Hostname: page-component-586b7cd67f-rdxmf Total loading time: 0 Render date: 2024-11-26T08:44:15.601Z Has data issue: false hasContentIssue false

Thermo-mechanical numerical model of the transition from continental rifting to oceanic spreading: the case study of the Alpine Tethys

Published online by Cambridge University Press:  03 October 2016

ANNA MARIA MAROTTA*
Affiliation:
Department of Earth Sciences ‘Ardito Desio’, Università degli Studi di Milano, via L. Mangiagalli 34 I-20134, Milan, Italy
MANUEL RODA
Affiliation:
Department of Earth Sciences ‘Ardito Desio’, Università degli Studi di Milano, via L. Mangiagalli 34 I-20134, Milan, Italy
KATYA CONTE
Affiliation:
Department of Earth Sciences ‘Ardito Desio’, Università degli Studi di Milano, via L. Mangiagalli 34 I-20134, Milan, Italy
MARIA IOLE SPALLA
Affiliation:
Department of Earth Sciences ‘Ardito Desio’, Università degli Studi di Milano, via L. Mangiagalli 34 I-20134, Milan, Italy
*
*Author for correspondence: [email protected]
Rights & Permissions [Opens in a new window]

Abstract

We develop a two-dimensional thermo-mechanical numerical model in which the formation of oceanic crust and serpentinite due to the hydration of the uprising mantle peridotite has been implemented, with the aim of discussing the behaviour of the lithosphere of the Alps and Northern Apennines during the transition from continental rifting to ocean spreading of the Alpine Tethys. The predictions of the model are compared with natural data related to the Permian–Triassic high-temperature – low-pressure (HT-LP) metamorphism affecting the continental lithosphere and data from the Jurassic P–T evolution of the oceanic lithosphere from the Alps and the Northern Apennines. Our analysis indicates that a thinned continental crust, an ocean–continent transition zone and an oceanic lithosphere characterize the final structure of the system in a poor magma rift pre-Alpine configuration. We also find that mantle serpentinization starts before crustal break-up and that denudation occurs before ocean spreading. The mantle denudation starts several million years before the gabbros/basalt formation, generating an ocean–continent transition zone from the passive continental margin to the oceanic lithosphere of size 160–280 km. The comparative analysis shows that the extension of a hot and weak lithosphere, which promotes the development of hyperextended Alpine margins, better agrees with the natural data. Finally, our comparative analysis supports the hypothesis that the lithospheric extension preceding the opening of the Alpine Tethys did not start in a stable continental lithosphere, but developed by recycling part of the old Variscan collisional suture.

Type
Original Articles
Copyright
Copyright © Cambridge University Press 2016 

1. Introduction

The aim of the present work is to discuss the thermo-mechanical behaviour of the lithosphere of the Alps and Northern Apennines during the transition from continental rifting to oceanization of the Alpine Tethys.

After the Variscan Orogeny, the future Alpine area (Fig. 1) underwent an extensional stage leading to the break-up of the Pangaea continental lithosphere and the opening of the Alpine Tethys Ocean (Lardeaux & Spalla, Reference Lardeaux and Spalla1991; Diella, Spalla & Tunesi, Reference Diella, Spalla and Tunesi1992; Dal Piaz, Reference Dal Piaz, Raumer and Neubauer1993; Bertotti et al. Reference Bertotti, Picotti, Bernoulli and Castellarin1993; Handy et al. Reference Handy, Franz, Heller, Janott and Zurbriggen1999; Schuster et al. Reference Schuster, Scharbert, Abart and Frank2001; Schuster & Stüwe, Reference Schuster and Stüwe2008; Marotta, Spalla & Gosso, Reference Marotta, Spalla, Gosso, Ring and Wernicke2009). The influence of the thermal, structural and compositional inheritance of the Variscan collision and collapse on the following extensional stage is still under debate (e.g. Marotta & Spalla Reference Marotta and Spalla2007; Von Raumer et al. Reference von Raumer, Bussy, Schaltegger, Schulz and Stampfli2013; Spalla et al. Reference Spalla, Zanoni, Marotta, Rebay, Roda, Zucali, Gosso, Schulmann, Martínez Catalán, Lardeaux, Janoušek and Oggiano2014).

Figure 1. Tectonic map of the Alps and Apennines (after Marotta & Spalla, Reference Marotta and Spalla2007; Handy et al. Reference Handy, Schmid, Bousquet, Kissling and Bernoulli2010) with the locations of: metamorphic rocks and main gabbro bodies of Permian–Triassic age occurring in the pre-Alpine continental crust of the Alps; and mantle rocks and oceanic gabbros from the Alps and Northern Apennines. The codes are defined in Tables 2–5.

The thermal structure of the Alpine lithosphere before the rifting event is poorly understood. Based on petrological analyses of the Ivrea crustal section, Handy et al. (Reference Handy, Franz, Heller, Janott and Zurbriggen1999) and Smye & Stockli (Reference Smye and Stockli2014) proposed that a series of thermal pulses after Carboniferous time affected the pre-rifting Alpine crust, potentially associated with and following the asthenosphere upwelling driven by hyperextension of the Adriatic margin during Late Triassic – Early Jurassic time. According to recent interpretations, the elevated temperatures in the deep crust of the Ivrea Zone may persist for millions of years, remaining close to the solidus for approximately 30 Ma (Klötzli et al. Reference Klötzli, Sinigoi, Quick, Demarchi, Tassinari, Sato and Günes2014) and thermally perturbing the Alpine lithosphere before the beginning of the rifting. Passive extension in the Europa-Adria system is thought to have been active during Triassic time (Handy & Zingg, Reference Handy and Zingg1991; Muntener, Hermann & Trommsdorf, Reference Muntener, Hermann and Trommsdorf2000; Muntener & Hermann, Reference Muntener, Hermann, Wilson, Whitmarsh, Taylor and Froitzheim2001; Montanini, Tribuzio & Anczkiewicz, Reference Montanini, Tribuzio and Anczkiewicz2006) or to have started during late Palaeozoic time (Lardeaux & Spalla, Reference Lardeaux and Spalla1991; Diella, Spalla & Tunesi, Reference Diella, Spalla and Tunesi1992; Dal Piaz, Reference Dal Piaz, Raumer and Neubauer1993; Marotta & Spalla, Reference Marotta and Spalla2007; Marotta, Spalla & Gosso, Reference Marotta, Spalla, Gosso, Ring and Wernicke2009). In contrast, based on data from Val Malenco, Corsica, Erro-Tobbio and Voltri Massif, Piccardo, Padovano & Guarnieri (Reference Piccardo, Padovano and Guarnieri2014) proposed that the pre-rift mantle lithosphere was equilibrated along an intermediate subcontinental geothermal gradient and the Permian high-temperature event(s) was followed by isobaric cooling, lasting more than 50 Ma (Manatschal, Reference Manatschal2004). This interpretation suggests the occurrence of a thermally equilibrated lithosphere before the rifting event (Lavier & Manatschal, Reference Lavier and Manatschal2006). This idea is supported by numerous rifting models proposed for the Alpine area (e.g. Beltrando, Rubatto & Manatschal, Reference Beltrando, Rubatto and Manatschal2010; Mohn et al. Reference Mohn, Manatschal, Beltrando, Masini and Kusznir2012) that were conceived based exclusively on data from the exploration of the Iberia passive margin (e.g. Boillot, Beslier & Girardeau, Reference Boillot, Beslier, Girardeau, Banda, Torné and Talwani1995; Hébert et al. Reference Hébert, Beaudoin, Rochon and Gardien2008).

Numerous numerical and analogue models of continental extension highlight the roles of different parameters and mechanisms in dictating the final geometry and style of the rifting and continental break-up. Among others, Reston & Morgan (Reference Reston and Morgan2004), Van Avendonk et al. (Reference Van Avendonk, Lavier, Shillington and Manatschal2009), Brune & Autin (Reference Brune and Autin2013) and Brune et al. (Reference Brune, Heine, Pérez-Gussinyé and Sobolev2014) focus on the role of the thermal state of the lithosphere; Buck (Reference Buck1991), Corti et al. (Reference Corti, Bonini, Sokoutis, Innocenti, Manetti, Cloetingh and Mulugeta2004), Nagel & Buck (Reference Nagel and Buck2004), Huismans, Buiter & Beaumont (Reference Huismans, Buiter and Beaumont2005), Huismans & Beaumont (Reference Huismans and Beaumont2011, Reference Huismans and Beaumont2014), Cloetingh et al. (Reference Cloetingh, Burov, Matenco, Beekman, Roure and Ziegler2013), Brune et al. (Reference Brune, Heine, Pérez-Gussinyé and Sobolev2014) and Liao & Gerya (Reference Liao and Gerya2015) analyse the effects of the composition, rheology and strength of the lower crust; Brune (Reference Brune2014) and Brune et al. (Reference Brune, Heine, Pérez-Gussinyé and Sobolev2014) investigate the role of the extensional velocity; Manatschal, Lavier & Chenin (Reference Manatschal, Lavier and Chenin2015) and Naliboff & Buiter (Reference Naliboff and Buiter2015) examine the role of structural and compositional inheritance of the system; and Escartín, Hirth & Evans (Reference Escartín, Hirth and Evans1997) and Pérez-Gussinyé et al. (Reference Pérez-Gussinyé, Reston, Phipps Morgan, Wilson, Whitmarsh, Taylor and Froitzheim2001, Reference Pérez-Gussinyé, Morgan, Reston and Ranero2006) focus on serpentinization. However, these works are not strictly focused on the opening of the Alpine Tethys and do not perform systematic comparisons between the model predictions and the natural data from the continental and oceanic crust of the Alps and Northern Apennines, including PT estimates, radiometric ages and lithological affinity.

This work represents an advancement of the work of Marotta, Spalla & Gosso (Reference Marotta, Spalla, Gosso, Ring and Wernicke2009) in which a rifting process following a continental collision and ending before the break-up of the continental crust is modelled. They implemented a bi-dimensional numerical geodynamic model to analyse the effects of an active extension during the Permian–Triassic period (300–220 Ma), assuming that the extension developed in a lithosphere already thermally and mechanically perturbed by a previous subduction-collision phase which occurred during the Variscan age, up to 300 Ma. Marotta, Spalla & Gosso (Reference Marotta, Spalla, Gosso, Ring and Wernicke2009)’s results supported the idea of an asymmetric rifting in which the Adriatic continental crust represented the hanging wall, although satisfactory and complete agreement with the natural geological data in terms of the coincidence of age, lithology and PT values was only obtained at the higher rate of forced extension (2 cm a−1).

Indeed, Marotta, Spalla & Gosso (Reference Marotta, Spalla, Gosso, Ring and Wernicke2009)’s model has two main limits: (a) the model does not evolve until the continental lithosphere break-up and subsequent ocean spreading; and (b) the conditions favourable to mantle partial melting have not been considered; the timing of the beginning of the new oceanic lithosphere was therefore not predicted. For these reasons, in the present work the transition from continental rifting to ocean spreading has been investigated using a two-dimensional (2D) thermo-mechanical numerical model in which the serpentinite formation due to the hydration of the upraising peridotite has been implemented and extension occurs over a mechanically unperturbed lithosphere. The predictions of the model have been compared with natural data related to the Permian–Triassic high-temperature – low-pressure (HT-LP) metamorphism affecting the continental lithosphere and data from the pre-Alpine Jurassic PT evolution of the oceanic crust from the Alps and the Northern Apennines.

2. Numerical model

2.a. Model set-up

To simulate the transition from rifting to oceanic spreading, a time-dependent 2D thermo-mechanical numerical model has been used in which the dynamics of the crust–mantle system have been investigated by numerical integration of the three fundamental equations of conservation of mass, momentum and energy:

(1) $$\begin{equation} \nabla \cdot \vec{u} & =& 0\\ \end{equation}$$
(2) $$\begin{equation}- \nabla P + \nabla \cdot \vec{\tau } + \rho \vec{g} & =& 0\\ \end{equation}$$
(3) $$\begin{equation} \rho {c_p}\left( {\frac{{\partial T}}{{\partial t}} + \vec{u} \cdot \nabla T} \right) & =& \nabla \cdot \left( {K\nabla T} \right) + \rho {H_d} \end{equation}$$

respectively, where $\vec{u}$ is the velocity, P is the pressure, $\vec{\tau }$ is the deviatoric stress, ρ is the density, $\vec{g}$ is the gravity acceleration, cp is the specific heat at constant pressure, T is the temperature, K is the thermal conductivity and Hd is the radiogenic heat production rate per unit mass.

Equations (1), (2) and (3) are solved using the 2D finite-elements code SubMar, which has been exhaustively described by Marotta, Spelta & Rizzetto (Reference Marotta, Spelta and Rizzetto2006). This numerical code uses the penalty function formulation to integrate the equation for the conservation of momentum and the streamline upwind/Petrov-Galerkin method to integrate the equation for the conservation of energy.

The marker in-cell technique has been used to compositionally differentiate crust and mantle rocks.

A viscous-plastic behaviour has been assumed for both materials. The effective viscosity is calculated as follows:

(4) $$\begin{equation} {\mu ^{{\rm{eff}}}} = {\mu _{{\rm{viscous}}}} = {\mu _{0,i}}\exp\left[ {\frac{{{E_i}}}{R}\left( {\frac{1}{T} - \frac{1}{{{T_0}}}} \right)} \right] \end{equation}$$

where μ 0, i and Ei are the reference viscosity at the reference temperature T 0 and the activation energy for the crust (i = c) and the mantle (i = m), respectively, with a maximum value defined by the plastic viscosity assumed equal to 1025 Pa s (Table 1).

Table 1. Material properties used in the 2-D numerical modelling

We account for a brittle behaviour of the crust to define the rheological conditions for mantle serpentinization only, as better specified in Section 2.b.

The material parameters are listed in Table 1. Initially, the lithosphere is mechanically unperturbed and laterally homogeneous. The initial thickness of the continental crust is assumed to be 30 km (Fig. 2), which is in agreement with models envisaging the beginning of extension-transtension onto a lithosphere characterized by a crust with a thickness of approximately 30 km (Muntener, Hermann & Trommsdorf, Reference Muntener, Hermann and Trommsdorf2000; Manatschal, Reference Manatschal2004).

Figure 2. (a) 2D geometry and numerical set-up of the model. (b) Thermal and rheological profiles at the beginning of the evolution for the hot (red lines) and cold (blue lines) models. The solid lines indicate the effective viscosity profiles, corresponding to the geotherms indicated by the dashed lines.

Two thermal settings are proposed here to satisfy two contrasting pre-rifting settings of the Alpine lithosphere characterized by different depths of the 1600 K isotherm: 80 and 220 km (Fig. 2b). We refer to these simulations as the hot and cold models, respectively. The initial thermal conditions correspond to an almost conductive thermal profile from 300 K at the surface to 1600 K at the base of the lithosphere; an initial homogeneous temperature of 1600 K is assumed below the lithosphere (Fig. 2b).

Boundary conditions are defined in terms of the temperature and velocity. A temperature of 300 K is fixed at the top of the crust and throughout the air–water layer, a temperature of 1600 K is fixed at the base of the model and zero flux is assumed through the lateral sides of the model. We apply an extension rate of 1.25 cm a−1 at both lateral sides of the model throughout the crustal thickness, resulting in a symmetric passive rifting with a total extension rate of 2.5 cm a−1, compatible with the magma-poor nature of the rift (e.g. Manatschal & Müntener, Reference Manatschal and Müntener2009). The 2D domain is closed vertically with shear-free conditions prescribed along the top and the bottom of the model domain, while both crust and lithospheric mantle are allowed to exit the model boundaries allowing the thinning of either crust and lithospheric mantle (Fig. 2a).

Considering that in the Alpine literature the onset of rifting is proposed to occur during 220–200 Ma (Müntener, Hermann & Trommsdorf, Reference Muntener, Hermann and Trommsdorf2000; Manatschal, Reference Manatschal2004; Montanini, Tribuzio & Anczkiewicz, Reference Montanini, Tribuzio and Anczkiewicz2006; Piccardo, Padovano & Guarnieri, Reference Piccardo, Padovano and Guarnieri2014) and that the oldest rocks of the ophiolitic associations belonging to the Liguria–Piemonte Ocean have been dated at 175–160 Ma (Tribuzio, Thirwall & Vannucci, Reference Tribuzio, Thirwall and Vannucci2004; Rossi et al. Reference Rossi, Lahondère, Cocherie, Caballero and Féraud2012; Kaczmarek, Müntener & Rubatto, Reference Kaczmarek, Müntener and Rubatto2008; Li et al. Reference Li, Faure, Lin and Manatschal2013, Reference Li, Faure, Rossi, Lin and Lahondère2015), we run the simulation over a time span of 60 Ma to match the time interval needed to generate the oldest gabbro intrusions in the natural system.

2.b. Conditions for mantle serpentinization

Magmatic-poor margins formed by the rifting of continental crust are characterized by the occurrence of serpentinized peridotites within a broad continent-ocean transition (e.g. Pérez-Gussinyé et al. Reference Pérez-Gussinyé, Reston, Phipps Morgan, Wilson, Whitmarsh, Taylor and Froitzheim2001; Manatschal, Reference Manatschal2004) and by the lithostratigraphy of the ophiolitic sequences (Mevel, Caby & Kienast, Reference Mevel, Caby and Kienast1978; Lagabrielle & Cannat, Reference Lagabrielle and Cannat1990; Chalot-Prat, Reference Chalot-Prat, Foulger, Natland, Presnall and Anderson2005; Manatschal et al. Reference Manatschal, Sauter, Karpoff, Masini, Mohn and Lagabrielle2011; Li et al. Reference Li, Faure, Lin and Manatschal2013). The role of serpentinization of the lithospheric mantle during continental rifting and the transition to oceanic spreading has been extensively discussed by Pérez-Gussinyé et al. (Reference Pérez-Gussinyé, Reston, Phipps Morgan, Wilson, Whitmarsh, Taylor and Froitzheim2001, Reference Pérez-Gussinyé, Morgan, Reston and Ranero2006). They assume that mantle serpentinization occurs when the overlaying crustal layer is under brittle conditions, so that faults can cut across the crust allowing hydrous fluids to penetrate the mantle, and mantle matches the appropriate pressure and temperature conditions for the stability field of serpentine.

In order to implement mantle serpentinization and the consequent rheological and compositional changes, we check whether the pressure and temperature of each mantle-type marker match the stability field of serpentine and the overlying crustal layer is under brittle conditions. To define whether the overlying layer is under brittle conditions, we use a simplified formulation of Byerlee's law criterion:

(5) $$\begin{equation} {\sigma _{{\rm{brittle}}}} = \beta \cdot y\quad \quad {\rm{with}}\quad \quad \beta = 16\;{\rm{MPa}}\;{\rm{k}}{{\rm{m}}^{ - 1}} \end{equation}$$

where y is the depth, and compare the brittle strength σbrittle with the temperature- and pressure-based plastic strength:

(6) $$\begin{equation} {\rm{Regime}} \Leftarrow \min \left\{ {{\sigma _{{\rm{brittle}}}},{\sigma _{{\rm{plastic}}}}} \right\}. \end{equation}$$

2.c. Conditions for partial melting

During an active extension of a continental lithosphere, the temperature in the lithospheric mantle increases as a consequence of the upwelling asthenospheric flow. If the pressure and temperature conditions of peridotite solidus are matched, partial melting of the lithospheric mantle occurs and gabbroic-basaltic melts form. We assume here that once the mantle partial melting occurs in the system, the oceanic lithosphere starts to form. This implies an instantaneous transfer of the mantle melt to the surface in agreement with the estimates of the rate of magma ascent across the continental and oceanic crust (Clague, Reference Clague1987; Turner et al. Reference Turner, George, Evans, Hawkesworth and Zellmer2000).

In order to individuate the beginning of oceanic spreading and to identify the extension of the partially molten mantle region, we check the pressure and temperature conditions of each mantle-type marker during the evolution. When the PT conditions reach the dry solidus field of peridotite (Rogers et al. Reference Rogers, Blake, Widdowson, Parkinson and Harris2008):

(7) $$\begin{equation} {P_{\rm{m}}} \le \alpha \cdot {T_{\rm{m}}} + \beta , \end{equation}$$

where α = 0.00789792857 GPa K−1 and β = – 11.1071202 GPa, its typology is changed from mantle type into potential partially molten mantle type.

The extension of the partially molten mantle region is shown in Figure 3 distinguishing the oceanic lithosphere (dark green), formed by serpentinized mantle hosting gabbros and basalts that can be produced once partial melting conditions are attained, from the mantle that serpentinized before the occurrence of partial melting (light green).

Figure 3. Successive stages of the tectonic evolution predicted by the hot (panels ai) and cold (panels bi) models at different times after the beginning of the forced extension. Black and red dashed lines correspond to 800 K and 1500 K isotherms, respectively. Ages refer to the time span from the beginning of the simulations.

In the present form, the compositional and rheological changes consequent to partial melting are not implemented.

3. Model predictions

3.a. Structural configurations

Figure 3 shows the structural configurations of the crust and lithosphere at different stages of the evolution of the hot model (panels ai, left side) and the cold model (panels bi, right side). The crustal boundaries coincide with the envelope of the crustal type markers. The base of the lithosphere is thermally defined by the 1600 K isotherm (red dashed line).

During the initial phase of forced extension in both models a progressive thinning of the crust is predicted, mostly concentrated around the position of the future ridge. The crustal thinning occurs very early in the cold model, approximately 1 Ma after the beginning of the extension (Fig. 3b1), while more than 10 Ma passes before thinning becomes significant in the hot model (Fig. 3a1). After 15.4 Ma (hot model, Fig. 3a2) and 4.4 Ma (cold model, Fig. 3b2), the thermal thinning is localized around the position of the future ridge. These times mark the beginning of mantle serpentinization in both models. The long time interval necessary in the hot model for the beginning of deformation localization therefore seems related to the higher thermal state at the base of the continental crust, which results in a different rheological behaviour of the shallow lithospheric mantle as discussed at the end of Section 2.c.

The progression of forced extension leads to the occurrence of the crustal break-up which occurs at 31.4 Ma in the hot model, approximately 16 Ma after the beginning of the mantle serpentinization (Fig. 3a3), and at 7.4 Ma in the cold model, only 3 Ma after the beginning of the mantle serpentinization (Fig. 3b3). The occurrence of crustal break-up is followed by the exhumation of the serpentinized mantle replacing the thinned continental crust at the floor of the basin. In the hot model, the mantle serpentinization progresses even at the base of the continental crust for more than 250 km away from the ridge. In contrast, in the cold model deep-seated serpentinization does not occur. This is a consequence of the thinner continental crust characterizing the hot model, allowing the occurrence of mantle at shallower depths with respect to the cold model in which continental lithospheric mantle resides at greater depths and PT conditions are inappropriate for serpentinization.

The pressure and temperature conditions at the base of the lithosphere in the hot model are favourable for partial melting for a relatively short time after the crustal break-up (approximately 5 Ma, Fig. 3a4). In contrast, in the cold model favourable pressure and temperature conditions are reached after a relatively long time from the break-up (approximately 15 Ma, Fig. 3b4). Although there are different time intervals between the different stages (thinning, serpentinization and partial melting) in the hot and cold models, once crustal thinning starts comparable time spans pass before the beginning of the mantle partial melting in both models (approximately 21 Ma and 18 Ma for the hot and the cold models, respectively; Fig. 4).

Figure 4. Thermal and velocity fields of the system at different times after the beginning of the forced extension for the hot (panels ai) and cold (panels bi) models. Time scale in the centre of the figure represents the simulation duration for both models, on which the onset of serpentinization, crustal break-up and partial melting are located by arrows (red for the hot and blue for the cold models). Time intervals between serpentinization and partial melting onset are indicated by dashed lines.

At the beginning of the extension, the low strength of the lithosphere in the hot model reduces the efficiency of stress transmission up to shallow depths, making the localization of crustal thinning slower than that for the stronger lithosphere of the cold model. Once the crustal thinning is localized, the successive evolution is dominated by local thermal gradients because this stage is comparable in the two models. After their formation, the mantle, which is depleted by partial melting (dark yellow markers for the partial melting area and light yellow markers for the depleted mantle in Fig. 3a5, b5), moves laterally below the expanding oceanic crust, forming the oceanic lithosphere (dark green markers in Fig. 3a5, b5).

Figure 5. Surface horizontal crustal velocity predicted by the hot (a) and cold (b) models at different times after the beginning of the forced extension.

Here we assume that the oceanic lithosphere forms after the beginning of the mantle partial melting, when gabbros, basalts and part of serpentinized mantle contribute to the formation of the oceanic crust. Based on the structural configuration of the system at the time of partial melting (Fig. 3a4, b4) and after a given time span of 10 Ma (Fig. 3a5 and b5 for the hot model and the cold model, respectively), it is possible to estimate the variation in time of the width of the oceanic lithosphere, characterized by a gabbro-basalt-bearing crust, and of the serpentinized mantle denuded before melt generation. In particular, 10 Ma after the onset of the mantle partial melting our results indicate a total basin width ranging over 360–480 km as a function of the initial thermal state of the lithosphere (hot or cold models, respectively). In both models the oceanic lithosphere extends for approximately 200 km, while the denuded serpentinized mantle covers a width of approximately 160 km in the hot model and 280 km in the cold one.

3.b. Thermo-mechanical evolution

The structural configuration discussed above can be better understood if the thermo-mechanical evolution of the system is analysed. Figure 4 shows the thermal and velocity fields of the system at different stages of the evolution for the hot and the cold models (Fig. 4ai and bi, respectively). During the early stages of the evolution of both models, the velocity field is controlled by the far-field traction driving a predominantly horizontal velocity pattern though the lithosphere, being the intensity of the mantle upwelling below the future ridge lower than half the intensity of the far-field for the cold model (panel b1) and even negligible for the hot model (panel a1). Mantle upwelling increases during the evolution (Fig. 4a2) and reaches a magnitude comparable to or higher than that of the far-field traction only after the onset of the mantle partial melting (Fig. 4a3). For both models, the thermal thinning of the lithosphere localizes at the future ridge position when the mantle serpentinization conditions are matched, at approximately 15 Ma (Fig. 4a1) and 4.4 Ma (Fig. 4b1) after the beginning of the forced extension for the hot and cold models, respectively. Thermal thinning achieves its maximum when the pressure and temperature conditions are favourable for the mantle partial melting at 36.4 Ma (Fig. 4a3) and 22.4 Ma (Fig. 4b3) for the hot and cold models, respectively. The thermal thinning in the two models therefore differs only in the initial stages and, once the thermal destabilization starts, the time span needed to reach the maximum thinning is approximately independent of the initial thermal state (thick red and blue lines along the time scale in Fig. 4).

The different behaviours of the cold and hot models during the initial deformation phase and their similar behaviours at longer times can be ascribed to the differences in the lithosphere strength (Fig. 2b). During the initial phase, the lower strength (in terms of the lower effective viscosities of the lithospheric mantle levels) that characterizes the hot model (solid red line in Fig. 2b) may attenuate the transmission of the far-field traction up to the future ridge position, making the lithosphere thinning a very slow process compared with the cold model.

3.c. Velocity configuration

A more detailed analysis of the surface crustal horizontal velocity (Fig. 5a, b) indicates that both models are characterized by an initial phase during which the mean surface crustal velocity is stationary and increases linearly from the future ridge outwards until it reaches, at the boundary of the model, the velocity value constrained by the far-field traction. This stationary pattern lasts a short time (approximately 5 Ma) in the cold model (Fig. 5a) and a long time (approximately 30 Ma) in the hot model (Fig. 5b), ending when the crustal break-up occurs. This initial phase is followed by a transition phase lasting approximately 9 Ma and 3 Ma for the cold and hot models, respectively, during which the far-field and upwelling flow concur to the extension rate and the mean surface crustal horizontal velocities progressively increase until they reach the value of the far-field traction velocity. At the advanced stages of evolution the surface crustal horizontal velocity remains almost constant from the ridge to the margins for both models, with the exception of a 150–200 km wide region around the ridge where the mean surface horizontal velocity of the crust is mainly controlled by the very intense mantle upwelling flow and may overcome the far-field value.

The transmission of the spreading rate towards the periphery of the domain is likely limited by the fixed far-field prescribed velocity at the boundaries of the model producing, at first sight, an artificial local shortening.

Figure 6 shows the maps of the velocity modules through the lithospheric thickness at different time steps of the evolution, for the hot (Fig. 6ai) and cold (Fig. 6bi) models. Until crustal break-up occurs, the linear increase from the ridge outwards observed in the mean surface crustal velocity (Fig. 5) characterizes the entire lithosphere thickness (Fig. 6a1, a2 for the hot model and Fig. 6b1, b2 for the cold model). A small increase in the velocity is evident within the mantle where a decrease in the viscosity occurs because of the serpentinization. After the crustal break-up, the strongest velocity gradients localize within 200 km around the ridge and are associated with the increase in the magnitude of the upwelling mantle flow (Fig. 4). In the proximity of the ridge, in both models the extension rate at the base of the crust increases by a factor of 2–4 since the continental break-up to oceanic spreading (compare Fig. 6a2, b2 to Fig. 6a3, b3). Further away, the velocity gradients decrease significantly until they disappear, and the entire continental lithosphere moves as a rigid plate at a rate equal to the traction velocity after the onset of the mantle partial melting. Significant velocity gradients persist within the part of the basin floored by the serpentinized mantle and at the deep lithosphere levels below the ridge, where the intensity of the asthenosphere upwelling overcomes the far-field traction velocity (compare Fig. 4a3, b3 to Fig. 6a4, b4). At the advanced stages of the evolution after the onset of the mantle partial melting, as already observed in Figure 5 for the mean surface horizontal velocity, even the lithosphere velocity overcomes the far-field traction velocity up to a maximum value of 1.4 cm a−1 compared with the prescribed 1.25 cm a−1.

Figure 6. Maps of the velocity modules throughout the lithosphere at different times after the beginning of the forced extension for the hot (panels ai) and cold (panels bi) models.

4. Natural data

The natural data, compared with the model predictions in the following discussion, are derived from continental and oceanic rocks from the Alps and Apennines and include:

  1. 1. lithotypes: useful for inferring the continental, oceanic, crustal or mantellic provenance of each considered rock type as well as the raw characterization of the associable lithostratigraphic setting;

  2. 2. PT conditions: useful for individuating potential provenance regions along the model-predicted lithospheric cross-sections based on the occurrence of compatible thermal states; and

  3. 3. geochronological data: useful for checking the time correspondences of the mineral assemblage developments with the predicted sequence of thermal states from rifting to oceanic spreading.

The kinematics of the natural structures supported by mineral assemblages, which are used to infer the considered PT estimates, are also reported as a supplementary check for compatibility with the investigated lithospheric extensional regime.

The rock types, mineral assemblages and climax metamorphic conditions are described in Tables 2–5. In the pre-Mesozoic continental lithosphere of the Alps and Apennines, the metamorphic T max imprints and igneous activity compatible with the high thermal regime induced by mantle upwelling during lithospheric thinning have Permian–Triassic ages. These estimated conditions have therefore been selected for comparison with model predictions. On the contrary, the age data from the oceanic lithosphere selected for comparison with model predictions are concentrated around Middle Jurassic ages and are referred to gabbro emplacement and metamorphic re-equilibration at the ocean floor.

Table 2. Permian–Triassic metamorphic rocks from the continental crust of the Alps. Labels listed in the code column are reported in Figures 1 and 9–11 to indicate rock positions and duration of thermal fitting

A synthetic geological outline is given to help readers who are not familiar with the Alps and Apennines to place the samples that have provided the natural data in an adequate tectonic frame.

4.a. Geological outline

The Alps and Apennines developed during the closure of the Mesozoic Tethys along two opposite subduction zones during different time intervals, during Cretaceous–Oligocene time for the Alpine system and from Eocene time to the present day for the Apennines system (e.g. Dal Piaz, Reference Dal Piaz, Beltrando, Peccerillo, Mattei, Conticelli and Doglioni2010; Handy et al. Reference Handy, Schmid, Bousquet, Kissling and Bernoulli2010; Carminati & Doglioni, Reference Carminati and Doglioni2012). Different erosion percentages and exhumation efficiencies (Carminati & Doglioni, Reference Carminati and Doglioni2012) allowed the exposure of deep structural levels in the axial zone of the Alpine nappe; however, these levels are still buried in the Apennines. In contrast to the Apennines, denudation of the axial part of the belt makes the Alps an important site where it is possible to explore continuous sections of deep and intermediate Alpine and pre-Alpine continental crust and to investigate the metamorphic and igneous effects of the lithospheric-thinning-induced high thermal regime preceding Alpine convergence.

The Alps spread out from the Gulf of Genova to the Vienna Basin; the chain was truncated during the Neogenic opening of the Ligurian–Provençal–Algerian and Tyrrhenian basins (e.g. Bozzo et al. Reference Bozzo, Campi, Capponi and Giglia1992; Séranne, Reference Seranne, Durand, Jolivet, Horvath and Seranne1999; Federico et al. Reference Federico, Spagnolo, Crispini and Capponi2009; Turco et al. Reference Turco, Macchiavelli, Mazzoli, Schettino and Pierantoni2012). Alpine subduction-collision is responsible for the distribution of the continental pre-Alpine and Mesozoic oceanic rocks in four tectonic domains, individuated and explored at the lithospheric scale after the seismic integrated geophysical prospection projects of the whole belt (e.g. Polino, Dal Piaz & Gosso, Reference Polino, Dal Piaz and Gosso1990; Schmid et al. Reference Schmid, Fügenschuh, Kissling and Schuster2004; Cassinis, Reference Cassinis2006). From the internal to the external part of the belt, they are the Southalpine, Austroalpine, Penninic and Helvetic domains (Fig. 1).

The Southalpine domain consists of a south-verging thrust system that has been active since Cretaceous time, involving Palaeozoic continental basement and Permian–Cenozoic cover units, both only locally affected by very-low-grade Alpine metamorphism. The Austroalpine and Penninic domains have been deeply involved in the Alpine subduction and collision system, as accounted for by the high-pressure metamorphic imprints associated with the dominant Alpine fabrics. The Penninic domain consists of mingled crustal slices deriving from both pre-Alpine continental and Mesozoic oceanic lithosphere, the latter tectonically sampled from the subducted Tethys Ocean (e.g. Platt Reference Platt1986; Polino, Dal Piaz & Gosso, Reference Polino, Dal Piaz and Gosso1990; Stöckhert & Gerya, Reference Stöckhert and Gerya2005; Malatesta et al. Reference Malatesta, Crispini, Federico, Capponi and Scambelluri2012; Roda, Spalla & Marotta, Reference Roda, Spalla and Marotta2012; Malatesta et al. Reference Malatesta, Gerya, Crispini, Federico and Capponi2013). In contrast, the Austroalpine domain does not contain Mesozoic ophiolites but is infolded within them and the related Mesozoic sediments all along its external boundary. Finally, the Helvetic domain consists of a Europe-verging thrust system that includes basement and cover slices structured during the late stages of the Alpine continental collision since Tertiary. Due to its shallow-level Alpine tectonic history, the Helvetic and Southalpine units (Fig. 1) broadly preserve pre-Alpine metamorphic, structural and stratigraphic imprints (Fig. 7), whereas Cretaceous–Paleocene Alpine high-pressure rocks are confined within the axial part of the chain in a rootless crustal prism consisting of Penninic and Austroalpine units. These latter are bounded by the Penninic frontal thrust (PF, Fig. 1) towards the European foreland and by the Periadriatic lineament (PL, Fig. 1) towards the Adriatic hinterland (Handy & Oberhänsli, Reference Handy and Oberhänsli2004; Thöni et al. Reference Thöni, Miller, Blichert-Toft, Whitehouse, Konzett and Zanetti2008; Roda, Spalla & Marotta, Reference Roda, Spalla and Marotta2012). According to many authors the wide range of radiometric ages of high-pressure metamorphic imprints (Bousquet et al. Reference Bousquet, Engi, Gosso, Oberhänsli, Berger, Spalla, Zucali and Goffè2004; Goffé et al. Reference Goffé, Schwartz, Lardeaux and Bousquet2004; Handy et al. Reference Handy, Schmid, Bousquet, Kissling and Bernoulli2010; Lardeaux, Reference Lardeaux2014) suggests that they were buried and widely exhumed during subduction of the oceanic lithosphere (European lower plate), accompanied by tectonic erosion of the upper continental plate (Adria) before the onset of continental collision (e.g. Platt Reference Platt1986; Polino, Dal Piaz & Gosso, Reference Polino, Dal Piaz and Gosso1990; Spalla et al. Reference Spalla, Lardeaux, Dal Piaz, Gosso and Messiga1996; Gerya & Stöckhert Reference Gerya and Stöckhert2005; Roda, Marotta & Spalla, Reference Roda, Marotta and Spalla2010; Roda, Spalla & Marotta, Reference Roda, Spalla and Marotta2012; Rubatto et al. Reference Rubatto, Regis, Hermann, Boston, Engi, Beltrando and McAlpine2011). Here, pre-Alpine structural, metamorphic and igneous relics are also preserved even if pervasive Alpine structural and metamorphic reworking shapes them into small-sized and scattered lenses, inhibiting the correlation of pre-Alpine structures at the regional scale.

Figure 7. Timing of the metamorphic and magmatic events from Variscan to Jurassic time along the Alpine belt. Radiometric ages are plotted including analytical uncertainty intervals (dashed horizontal bars). Variscan evolutions are synthesized by data in Table 2. Permian–Triassic metamorphic and igneous data (continental gabbros) are synthesized by data in Table 3 and represented with their error margin. Ophiolites and mantle age data are listed in Tables 4 and 5. The grey stripe indicates the age of the earliest radiolarian cherts related to Alpine ophiolites (Cordey & Bailly, Reference Cordey and Bailly2007).

4.b. Lithology, structures, PT conditions and ages

Asthenospheric upwelling associated with lithospheric thinning causes high thermal regimes in the thinned continental lithosphere (e.g. Thompson, Reference Thompson1981; England & Thompson, Reference England and Thompson1984; Thompson & England, Reference Thompson and England1984; Sandiford & Powell, Reference Sandiford and Powell1986; Spear & Peacock, Reference Spear and Peacock1989; Beardsmore & Cull, Reference Beardsmore and Cull2001), and in the Alps the only igneous and metamorphic effects indicating the occurrence of such a thermal state before the Alpine convergence are of Permian–Triassic age. They are detectable along the whole belt, from the Ligurian Sea to the Pannonian Basin, even in domains strongly reworked by the Alpine tectonics and metamorphism. These records consist of a widespread emplacement of Permian–Triassic basic to acidic igneous activity and huge gabbro bodies (Tables 2 and 3; Fig. 1) associated with regional high-temperature – low-pressure (HT-LP) metamorphism, which postdate structures and metamorphic imprints widely developed during the Variscan subduction and collision (Fig. 7). These features are frequently associated with subcontinental peridotites (Table 4) and are mainly confined to the Austroalpine and Southalpine domains (e.g. Lardeaux & Spalla, Reference Lardeaux and Spalla1991; Bonin et al. Reference Bonin, Brändlein, Bussy, Desmons, Eggenberger, Finger, Graf, Marro, Mercolli, Oberhänsli, Ploquin, Quadt, Raumer, Schaltegger, Steyrer, Visonà, von Raumer and Neubauer1993; Bussy et al. Reference Bussy, Venturini, Hunziker and Martinotti1998; Rottura et al. Reference Rottura, Bargossi, Caggianelli, Del Moro, Visonà and Tranne1998; Schuster et al. Reference Schuster, Scharbert, Abart and Frank2001; Stahle et al. Reference Stahle, Frenzel, Hess, Saupe, Schmidt and Schneider2001; Rebay & Spalla, Reference Rebay and Spalla2001; Rampone, Reference Rampone2002; Peressini et al. Reference Peressini, Quick, Sinigoi, Hofmann and Fanning2007; Marotta, Spalla & Gosso, Reference Marotta, Spalla, Gosso, Ring and Wernicke2009; Spalla et al. Reference Spalla, Zanoni, Marotta, Rebay, Roda, Zucali, Gosso, Schulmann, Martínez Catalán, Lardeaux, Janoušek and Oggiano2014).

Table 3. Permian–Triassic gabbros emplaced in the pre-Alpine continental crust of the Alps and Apennines. Labels listed in the code column are reported in Figures 1 and 9–11 to indicate rock positions and duration of thermal fitting

Table 4. Subcontinental and oceanic mantle peridotites from the Alps and Apennines. Labels listed in the code column are reported in Figures 1 and 9–11 to indicate rock positions and duration of thermal fitting

Metamorphic Permian–Triassic imprints are widely recorded in the lower, intermediate and upper continental crust of the Austroalpine and Southalpine domains; only a few records have been recognized in the upper and intermediate Penninic crust of Western Alps, and they are never detected in the Helvetic domain (Fig. 1). In the Penninic domain, HT assemblages mainly developed in sillimanite-bearing metapelites and metaintrusives. They are generally derived from re-equilibration under low-pressure conditions and interpreted as Permian in age (Bouffette, Lardeaux & Caron, Reference Bouffette, Lardeaux and Caron1993; Table 2). Sapphirine-bearing HT- to UHT-IP granulites of the Gruf Complex have been recently petrologically and chronologically investigated, revealing an age of 260–290 Ma for T max conditions (Galli et al. Reference Galli, Le Bayon, Schmidt, Burg, Caddick and Reusser2011). The exhumation of Penninic HT Permian–Triassic rocks was generally associated both with cooling and heating (Desmons, Reference Desmons1992; Bouffette, Lardeaux & Caron, Reference Bouffette, Lardeaux and Caron1993) and occurred before the Alpine convergence, whereas the exhumation of UHT Gruf granulites occurred at the end of the Alpine convergence from the base of the internal European passive margin where they lay since the Permian lithospheric thinning (Galli et al. Reference Galli, Le Bayon, Schmidt, Burg, Caddick and Reusser2011). In the Austroalpine domain, Permian–Triassic HT metamorphism mainly developed in sillimanite- and biotite-bearing gneisses; it is associated with minor mafic granulites, amphibolites and high-grade marbles (Table 2). HT minerals locally mark mylonitic fabrics within discrete shear zones (Lardeaux & Spalla Reference Lardeaux and Spalla1991; Spalla et al. Reference Spalla, Lardeaux, Dal Piaz and Gosso1991), and exhumation paths can be characterized by cooling, heating or isothermal decompression (e.g. Dal Piaz, Lombardo & Gosso, Reference Dal Piaz, Lombardo and Gosso1983; Stöckhert, Reference Stöckhert1987; Vuichard, Reference Vuichard1987; Lardeaux & Spalla, Reference Lardeaux and Spalla1991; Spalla, Messiga & Gosso, Reference Spalla, Messiga and Gosso1995; Schuster et al. Reference Schuster, Scharbert, Abart and Frank2001; Manzotti & Zucali, Reference Manzotti and Zucali2013). Locally, the exhumation was accomplished up to very shallow structural levels (e.g. Rebay & Spalla, Reference Rebay and Spalla2001), suggesting that some Austroalpine units belonged to a thinned continental crust before being subducted during Cretaceous convergence. In the Southalpine basement, Permian–Triassic HT metamorphism developed in metapelites, mafic granulites, amphibolites and high-grade marbles and mainly re-equilibrated under granulite-amphibolite-facies conditions (Table 2). HT paragenesis marks a pervasive foliation that is locally mylonitic within up to kilometre-thick discrete shear zones that are often steepened by Alpine thrusting (e.g. Bertotti et al. Reference Bertotti, Picotti, Bernoulli and Castellarin1993; Gosso, Siletto & Spalla, Reference Gosso, Siletto and Spalla1997). Mylonitic belts developed under upper amphibolite-granulite- to greenschist-facies conditions are widespread in the central Southalpine domain, accounting for regional-scale pervasive extensional tectonics (Brodie, Rex & Rutter, Reference Brodie, Rex, Rutter, Coward, Dietrich and Park1989; Diella, Spalla & Tunesi, Reference Diella, Spalla and Tunesi1992; Bertotti et al. Reference Bertotti, Picotti, Bernoulli and Castellarin1993) and widely interpreted as related to regional-scale normal faults responsible for the exhumation of HT-LP complexes (e.g. Brodie, Rex & Rutter, Reference Brodie, Rex, Rutter, Coward, Dietrich and Park1989; Handy et al. Reference Handy, Franz, Heller, Janott and Zurbriggen1999) during a continuous evolution from Permian to Triassic or Jurassic time (e.g. Bertotti et al. Reference Bertotti, Picotti, Bernoulli and Castellarin1993). Intrusive stocks emplaced at shallow levels during Permian – Early Triassic time are associated with metamorphic aureoles (Povoden, Horacek & Abart, Reference Povoden, Horacek and Abart2002; Benciolini et al. Reference Benciolini, Poli, Visona and Zanferrari2006; Gallien, Abart & Wyhlidal, Reference Gallien, Abart and Wyhlidal2007). Exhumation paths were characterized by cooling or by increasing temperature during decompression (e.g. Brodie, Rex & Rutter, Reference Brodie, Rex, Rutter, Coward, Dietrich and Park1989; di Paola & Spalla Reference di Paola and Spalla2000) and were generally accomplished under a high thermal regime. HT Permian–Triassic mineral associations have never been detected in the pebbles of Permian conglomerates from the Orobic Alps, suggesting that HT rocks were not yet exposed in their Variscan source areas before the late Permian – Triassic period (Spalla et al. Reference Spalla, Zanoni, Gosso and Zucali2009; Zanoni, Spalla & Gosso, Reference Zanoni, Spalla and Gosso2010).

Permian–Triassic continental gabbros have been detected in the Austroalpine and Southalpine domains of the Alps (Bonin et al. Reference Bonin, Brändlein, Bussy, Desmons, Eggenberger, Finger, Graf, Marro, Mercolli, Oberhänsli, Ploquin, Quadt, Raumer, Schaltegger, Steyrer, Visonà, von Raumer and Neubauer1993; Rottura et al. Reference Rottura, Bargossi, Caggianelli, Del Moro, Visonà and Tranne1998; Stahle et al. Reference Stahle, Frenzel, Hess, Saupe, Schmidt and Schneider2001; Spiess et al. Reference Spiess, Cesare, Mazzoli, Sassi and Sassi2010; Spalla et al. Reference Spalla, Zanoni, Marotta, Rebay, Roda, Zucali, Gosso, Schulmann, Martínez Catalán, Lardeaux, Janoušek and Oggiano2014). The mafic products of the widespread Permian–Triassic igneous activity mainly consist of gabbroic bodies with subcontinental peridotites (Brodie, Rex & Rutter, Reference Brodie, Rex, Rutter, Coward, Dietrich and Park1989; Bonin et al. Reference Bonin, Brändlein, Bussy, Desmons, Eggenberger, Finger, Graf, Marro, Mercolli, Oberhänsli, Ploquin, Quadt, Raumer, Schaltegger, Steyrer, Visonà, von Raumer and Neubauer1993; Schuster et al. Reference Schuster, Scharbert, Abart and Frank2001; Stahle et al. Reference Stahle, Frenzel, Hess, Saupe, Schmidt and Schneider2001; Rampone Reference Rampone2002; Spalla et al. Reference Spalla, Zanoni, Marotta, Rebay, Roda, Zucali, Gosso, Schulmann, Martínez Catalán, Lardeaux, Janoušek and Oggiano2014). They occurred at different structural levels, and the country rocks vary from high-temperature – intermediate-pressure metamorphics to Triassic carbonatic sediments (Sills Reference Sills1984; Handy & Zingg Reference Handy and Zingg1991; Lardeaux & Spalla Reference Lardeaux and Spalla1991; Gallien, Abart & Wyhlidal, Reference Gallien, Abart and Wyhlidal2007; Miller et al. Reference Miller, Thöni, Goessler and Tessadri2011; Spalla et al. Reference Spalla, Zanoni, Marotta, Rebay, Roda, Zucali, Gosso, Schulmann, Martínez Catalán, Lardeaux, Janoušek and Oggiano2014). From a geochemical point of view most of the gabbros have a tholeiitic to alkaline signature, although they are generally considered to be generated from variably contaminated mantle sources in an extensional tectonic regime under a high thermal state associated with lithospheric thinning and rifting (Spalla et al. Reference Spalla, Zanoni, Marotta, Rebay, Roda, Zucali, Gosso, Schulmann, Martínez Catalán, Lardeaux, Janoušek and Oggiano2014 and references in Table 3). The ages of these gabbroic intrusions in the Austroalpine domain cluster around Permian (Table 3). The main Triassic magmatic signal is recorded in the Southalpine domain by Predazzo and Monzoni intrusives and in the Ivrea Zone by alkaline rocks (Table 3). Stahle et al. (Reference Stahle, Frenzel, Hess, Saupe, Schmidt and Schneider2001) interpret this igneous activity as the result of a long-lasting process of active rifting.

Few Permian–Triassic continental mantle slices have been found in the Alps (Table 4). In particular, the North Lanzo body in the Penninic domain is characterized by subcontinental lithospheric mantle protoliths and underwent progressive exhumation during pre-oceanic lithosphere extension and rifting. The South Lanzo body shows impregnation by mid-ocean-ridge basalt (MORB) melts rising from the underlying molten asthenosphere during the rifting stage of the Liguria–Piemonte Ocean (Piccardo & Guarnieri Reference Piccardo and Guarnieri2010). The Lanzo peridotites would therefore represent a lithosphere changing from a thinned continental plate to an ocean–continent transition zone (OCTZ) (Piccardo & Guarnieri Reference Piccardo and Guarnieri2010). In the Austroalpine domain, a small peridotite body has been detected in the Dent–Blanche Nappe (Nicot, Reference Nicot1977). Although no radiometric age is available for this mantle rock, its structural relation with continental rocks of the Valpelline Series and the analogy with similar rocks of the Ivrea Zone suggest a Permian age and therefore a continental affinity for this peridotite.

In the Alps, the transition from rifting to oceanic spreading is indicated by the deposition of post-rift sediments (172–165 Ma, Baumgartner et al. Reference Baumgartner, Bartolini, Carter, Conti, Cortese, Danelian, De Wever, Dumitrica, Dumitrica-Jud, Gorican, Guex, Hull, Kito, Marcucci, Matsuoka, Murchey, O'Dogherty, Savary, Vishnevskaya, Widz and Yao1995; Stampfli et al. Reference Stampfli, Mosar, Marquer, Marchant, Baudin and Borel1998; Bill et al. Reference Bill, O'Dogherty, Guex, Baumgartner and Masson2001; Handy et al. Reference Handy, Schmid, Bousquet, Kissling and Bernoulli2010) postdating syn-rift Triassic deposits (e.g. Gillcrist, Coward & Mugnier, Reference Gillcrist, Coward and Mugnier1987). This transition involved the exhumation and serpentinization of the subcontinental mantle at the Liguria–Piemonte Ocean margins (e.g. Desmurs, Manatschal & Bernoulli, Reference Desmurs, Manatschal, Bernoulli, Wilson, Whitmarsh, Taylor and Froitzheim2001; Manatschal, Reference Manatschal2004; Manatschal & Müntener, Reference Manatschal and Müntener2009). The radiometric ages of the ophiolitic gabbros (Fig. 7, Table 5) clustering around approximately 160 Ma (Mevel, Caby & Kienast, Reference Mevel, Caby and Kienast1978; Li et al. Reference Li, Faure, Lin and Manatschal2013), with older values of 166–183 Ma from the Apennines, Corsica and Erro–Tobbio ophiolitic units (e.g. Tribuzio, Thirwall & Vannucci, Reference Tribuzio, Thirwall and Vannucci2004; Rampone et al. Reference Rampone, Borghini, Romairone, Abouchami, Class and Goldstein2014; Li et al. Reference Li, Faure, Rossi, Lin and Lahondère2015), confine the beginning of the spreading of the Liguria–Piemonte Ocean. Ages of 198 ± 22 (Sm–Nd on gabbro WR) have been obtained by Costa & Caby (Reference Costa and Caby2001) in ophiolites from the Western Alps (Chenaillet) and interpreted by the authors as the signature of lithospheric extension announcing the oceanic spreading. Oceanic gabbros exclusively occur in the ophiolitic sequences of the Penninic domain (Fig. 1). They generally have a MORB affinity and are usually associated with serpentinized mantle but sometimes to volcanic sequences, such as lava flows, pillow basalts and pillow breccias, and oceanic sediments (e.g. Mevel, Caby & Kienast, Reference Mevel, Caby and Kienast1978; Ohnenstetter et al. Reference Ohnenstetter, Ohnenstetter, Vidal, Cornichet, Hermitte and Mace1981; Riccardi, Tribuzio & Caucia, Reference Riccardi, Tribuzio and Caucia1994; Martin, Tartarotti & Dal Piaz Giorgio, Reference Martin, Tartarotti and Dal Piaz Giorgio1994).

Table 5. Ophiolitic gabbros from the Alps, Apennines and Corsica. The lack of P–T estimates for Corsica gabbros inhibits their comparison with the model predictions. Labels listed in the code column are reported in Figures 1 and 9–11 to indicate rock positions and duration of thermal fitting

Oceanic mantle rocks coming from the Penninic domain (Table 4) are variably serpentinized peridotites. Based on the relict texture, mineralogy and structural relations with gabbros and rodingites, the serpentinized peridotite of Mt Avic (Table 4) is considered to have an oceanic affinity (Fontana, Panseri & Tartarotti, Reference Fontana, Panseri and Tartarotti2008). Although the gabbroic bodies of the Lanzo Massif are considered to have originated during the opening of the Jurassic–Piedmont Ligurian ocean (Lagabrielle, Fudral & Kienast, Reference Lagabrielle, Fudral and Kienast1989; Pognante, Rösli & Toscani, Reference Pognante, Rösli and Toscani1985) or during the earliest stages of the formation of the embryonic oceanic crust (Kaczmarek, Müntener & Rubatto, Reference Kaczmarek, Müntener and Rubatto2008), the associated peridotites show a more complex affinity as already discussed. The Mg-rich gabbroic rocks of Erro-Tobbio complex in Voltri Massif (Ligurian Alps) is considered as representative of syn-rift melt intrusions in thinned lithospheric mantle exhumed at ocean–continent transition domains (Rampone et al. Reference Rampone, Borghini, Romairone, Abouchami, Class and Goldstein2014). They represent the oldest gabbroic bodies of the Alpine Tethys (201–163 Ma, Rampone et al. Reference Rampone, Borghini, Romairone, Abouchami, Class and Goldstein2014). Because rocks from the Alps are widely deformed and metamorphosed during Alpine subduction and collision, we add ophiolites from the Northern Apennines that escaped pervasive subduction-related metamorphic re-equilibration and the unique case of mafic granulite from the continental crust to the collection of natural data (Tables 3–5). The Northern Apennines are characterized by oceanic and continental units (Fig. 1). Oceanic units are divided into two different groups of thrust nappes – the Internal and External Ligurian units (e.g. Marroni & Pandolfi, Reference Marroni and Pandolfi2007) – strongly deformed under low-grade metamorphic conditions (Marroni & Pandolfi, Reference Marroni and Pandolfi2007; Donatio, Marroni & Rocchi, Reference Donatio, Marroni and Rocchi2013). An ophiolitic sequence of Jurassic age and a sedimentary cover ranging in age from Late Jurassic to Paleocene characterize the Internal Ligurian units (Marroni & Pandolfi, Reference Marroni and Pandolfi2007). In the External Ligurian units, Late Cretaceous sedimentary melanges containing slide-blocks of ophiolites occur at the base of the Upper Cretaceous carbonate turbidites (Helminthoid Flysch; Marroni & Pandolfi, Reference Marroni and Pandolfi2007). The Ligurian units were thrust onto the Tuscan nappes during the Oligo-Miocene post-collisional convergence. The successions of the Tuscan units belong to the Adria passive continental margin, recording in sequence rifting-, cooling- and subsidence-related imprints during the opening of the Liguria–Piemonte Ocean (Marroni & Pandolfi, Reference Marroni and Pandolfi2007). Underlying Tertiary metamorphic units are exposed in rare tectonic windows (Fig. 1).

Granulitic gabbros of the External Liguride Unit associated with felsic granulites locally intrude mantle peridotites (Marroni & Tribuzio, Reference Marroni and Tribuzio1996; Marroni et al. Reference Marroni, Molli, Montanini and Tribuzio1998). The emplacement of the gabbroic protoliths occurs at deep crustal levels of late Carboniferous – early Permian age (approximately 290 Ma), and their country rocks are felsic granulites of the lower continental crust. Most likely in association with the subcontinental mantle, mafic and felsic granulites underwent a multistage exhumation beginning during Permian–Triassic time and ending during Late Triassic – Middle Jurassic time, when they were finally exhumed to shallow levels by extensive brittle faulting (Marroni et al. Reference Marroni, Molli, Montanini and Tribuzio1998). The External Liguride Unit is interpreted as an ocean–continent transition zone (Marroni et al. Reference Marroni, Molli, Montanini and Tribuzio1998).

Some Jurassic oceanic gabbros have been detected in the Apennines (Table 5; Riccardi, Tribuzio & Caucia, Reference Riccardi, Tribuzio and Caucia1994; Tribuzio, Riccardi & Ottolini, Reference Tribuzio, Riccardi and Ottolini1995; Tribuzio, Riccardi & Messiga, Reference Tribuzio, Riccardi and Messiga1997; Rebay, Riccardi & Spalla, Reference Rebay, Riccardi and Spalla2015). The oldest gabbro bodies (169–179 Ma) belong to the External Liguride Unit in the Northern Apennines. Despite their N-MORB affinity, Tribuzio, Thirwall & Vannucci (Reference Tribuzio, Thirwall and Vannucci2004) interpreted these gabbros as having developed during an intermediate stage of the rifting process that led to the opening of the Ligurian Tethys and not strictly related to the oceanic spreading.

The peridotite of the External Liguride Unit has been interpreted as the mantle of an ocean–continent transition zone (Marroni et al. Reference Marroni, Molli, Montanini and Tribuzio1998), while the serpentinized peridotite of the Internal Liguride Unit is suggested to be representative of the Jurassic oceanic lithosphere of the Liguria–Piemonte Ocean (Marroni & Pandolfi Reference Marroni and Pandolfi2007; Donatio, Marroni & Rocchi, Reference Donatio, Marroni and Rocchi2013).

5. Comparison between the model predictions and the natural data

Before proceeding with the comparison between the model predictions and the natural data, it is necessary to reference the relative time of the numerical simulation with respect to the natural ages. To achieve this goal, we decided to match the beginning of the mantle partial melting in the model with the ages of the gabbros available in the literature (Table 5). In particular, we chose the three absolute ages of 160 Ma, 170 Ma and 185 Ma interpreted by different authors (Tribuzio, Riccardi & Ottolini, Reference Tribuzio, Riccardi and Ottolini1995; Tribuzio, Thirwall & Vannucci, Reference Tribuzio, Thirwall and Vannucci2004; Li et al. Reference Li, Faure, Lin and Manatschal2013, Reference Li, Faure, Rossi, Lin and Lahondère2015) as the oldest gabbro ages of Liguria–Piemonte Ocean; they can therefore be considered the temporal markers of the early oceanic spreading. Figure 8 depicts the time referencing of both models with respect to the three chosen absolute ages.

Figure 8. Time references from the hot and cold models with respect to the age of the oldest gabbros (160 Ma, 170 Ma and 185 Ma), based on the literature for the Northern Apennines and Western Alps (Tribuzio, Riccardi & Ottolini, Reference Tribuzio, Riccardi and Ottolini1995; Li et al. Reference Li, Faure, Lin and Manatschal2013; Rampone et al. Reference Rampone, Borghini, Romairone, Abouchami, Class and Goldstein2014). Green circles identify the times when the mantle partial melting occurs in the hot model (36.4 Ma) and in the cold model (22.4 Ma).

The data from the continental crust (Tables 2, 3) are compared with the continental markers; the data from the mantle (Table 4) are compared with the dry or serpentinized mantle according to the estimated PT conditions and rock assemblages. In the following, the process starting from the exhumation of serpentinized mantle and proceeding through the oceanic spreading is referred to as ‘oceanization’. The oceanic spreading leading to the formation of the Liguria–Piemonte Ocean starts when the conditions favourable for mantle partial melting are attained in the system. We compare the pressure–temperature values predicted for markers belonging to the oceanic lithosphere with the PT estimates available for gabbros belonging to Alpine and Apennine ophiolitic complexes (Table 5), as well as those available for the oceanic mantle (Table 4).

To verify if the simulated geodynamic context can reproduce a thermal state that is compatible with that recorded by the lithosphere during the Permian–Triassic period, in Figure 9 we show the duration of the agreement between the predictions and the natural data for the three different absolute ophiolite ages in terms of lithology and coincident PT values compared to the radiometric (black thick segments) and geologic (grey thick segments) ages of the natural data. In the following discussion, we refer to a ‘complete fit’ when there is agreement between the model predictions and the natural data in the lithology, PT values and ages. The fit is considered ‘partial’ if age coincidence is lacking. The comparative analysis takes into account the natural data with their error margins (Tables 2–5).

Figure 9. Duration of the agreement between the predictions and the natural data as well as number of fitting markers (colours) for the three different ages of the oldest gabbros (160 Ma, 170 Ma and 185 Ma, dashed black lines) in terms of lithological affinity and coincident PT values compared to the radiometric (black thick segments) and geologic (grey thick segments) ages of the natural data. In the text, we refer to the result as a ‘complete fit’ when the model predictions and the natural data of the lithological affinity, PT values and ages agree. The fit is considered ‘partial’ if age agreement is lacking. Panel (a) refers to the hot model, whereas panel (b) refers to the cold model. The light grey area represents the duration of the numerical simulation. The codes are defined in Tables 2–5.

The predictions of both the hot and cold models show complete fits with the same data, although the number of fitting markers is in general lower in the cold model than in the hot model. Complete fits are realized with a maximum of 13 data points out of the available 44: 10 are of the oceanic lithosphere type (Pp5a, Pp5b, Pp6, Pg1b and Pg2a, Penninic domain; and APN1a, APN1b, APN1c, APN1d and APN3, Apennine domain), and 2 are of the pre-oceanic lithosphere type (Ap5a, Pg3b and Pg3c, Austroalpine domain). Good agreement is obtained for all of the oceanic lithosphere type markers with the natural data from the Alpine and Apennine ophiolites for all of the ages proposed in the literature.

The predictions do not show a complete fit with the continental crust data. Nevertheless, the hot model reproduces thermo-barometric conditions that are compatible with those of most of the continental crust data in the interval during 220–150 Ma, according to the time of the oldest ophiolitic ages. A similar situation holds for the cold model, although with a very low number of markers. However, the radiometric ages of these data are older than the oldest model predictions.

Some data, that is, Sp1 (Eisacktal), Sp3 (Dervio-Olgiasca) and Sp5a and Sp5b (Ivrea) which are from the Southalpine domain, Ap17 (Sopron) and Ap7 and Ap8 (Valpelline) which are from the Austroalpine domain and Pp4 (Gruf) which is from the Penninic domain, do not show any agreement, even thermo-barometric, because the estimated temperatures are higher than the predicted temperatures.

Figures 10 and 11 depict the spatial distributions of the markers that exhibit complete fits with the data for the hot and the cold models, respectively. With the exception of the Pg3c and Ap5a (continental mantle, Austroalpine domain) data, the PT conditions that are compatible with those of the other natural data are predicted only in a 200 km wide region centred at the ridge and only after the formation of the oceanic lithosphere (Figs 10d, e, 11d, e). The Pg3c datum fits the mantle-type markers at shallow structural levels only after a significant thermo-mechanical thinning of the lithosphere. This occurs after approximately 5.4 Ma in the hot model, affecting a large amount of the continental lithospheric mantle (Fig. 10b). In the cold model, the fitting starts very early and affects only a 60 km wide area around the future ridge, concurrently with the thermal thinning localization. For both models, the fitting persists until the late stages of the evolution and affects the shallow structural levels of the entire non-serpentinized lithosphere. The hot model predicts PT conditions that are compatible with the Ap5a datum (continental mantle, Austroalpine domain) starting at the beginning of the simulation (Fig. 10a), while in the cold model the favourable conditions occur only after significant heating of the system (Fig. 11b). This result supports the idea that Ap5a is a representative slice of the continental mantle under the perturbed PT conditions rather than under the thermal regime of a stable lithosphere. Pg1b and Pg2a, both of which belong to the gabbroic rocks of the Lanzo Massif for both the hot and cold models, exhibit a complete fit with markers of the oceanic lithosphere type but at different structural levels. Pg1b agrees at deeper structural levels (Figs 10d, 11d), consistent with a gabbroic melt impregnation of the lithosphere (Compagnoni, di Brozolo & Sandrone, Reference Compagnoni, di Brozolo and Sandrone1984). Instead, Pg2a agrees at shallower levels (Figs 10d, 11d), in agreement with a decompression stage of an older subcontinental lithospheric mantle under hydrated conditions (Pognante, Rösli & Toscani, Reference Pognante, Rösli and Toscani1985). Fits with Pp6 (Mont Avic serpentinized peridotite), APN1a (gabbros of North Apennines ophiolites) and APN1b (gabbros of Apennines Ophiolites) data start at the beginning of the oceanization (Figs 10d, 11d) and can be compatible with either the oceanic lithosphere or the ocean–continent transition zone. During oceanic spreading, the complete fit localizes at the marginal portions of the oceanic basin (Figs 10d, e, 11d, e). Finally, the complete fit with the APN1c, APN1d and APN3 (gabbros and mantle of Apennines ophiolites) and Pp5a and Pp5b (Chenaillet) data occurs during oceanic spreading and at the appropriate litho-structural levels of the oceanic lithosphere (Figs 10e, 11e).

Figure 10. Spatial distribution of the markers that guarantee a complete fit of the data for the hot model at different stages of the evolution. The dashed black lines indicate the 800 K and 1500 K isotherms. The colours identify the markers showing a complete fit with the natural data as specified in the legend; t r indicates the time span from the beginning of the simulation and t a indicates the absolute age constrained choosing the oldest gabbro age at 170 Ma (see discussion in the text). The codes are defined in Tables 2–5.

Figure 11. Spatial distribution of the markers that guarantee a complete fit of the data for the cold model at different stages of the evolution. The dashed black lines indicate the 800 K and 1500 K isotherms. The colours identify the markers showing a complete fit with the natural data specified in the legend; t r indicates the time span from the beginning of the simulation and t a indicates the absolute age constrained choosing the oldest gabbro age at 170 Ma (see discussion in the text). The codes are defined in Tables 2–5.

6. Discussion

The model of crustal extension presented here, characterized by a weak lower crust and mantle serpentinization, results in symmetric rifting of the continental lithosphere and exhumation of a serpentinized lithospheric mantle. Our results support the idea that the occurrence of serpentinization of the mantle can favour the exhumation of the lithospheric mantle before the oceanic spreading in agreement with Lagabrielle & Cannat (Reference Lagabrielle and Cannat1990) and Pérez-Gussinyé et al. (Reference Pérez-Gussinyé, Morgan, Reston and Ranero2006).

The onset of lithospheric thinning localized around the future ridge strongly depends on the initial lithosphere thermal state: for a cold and strong lithosphere, the thinning is very rapid (after approximately 4.4 Ma) with respect to a hot and weak lithosphere (after approximately 15.4 Ma). Similarly, the time span between the onset of thinning and the occurrence of crustal break-up is shorter for a cold lithosphere (approximately 3 Ma) than for a hot lithosphere (approximately 16 Ma). These dynamics are attributable to the concurrent roles of the prescribed far-field traction and mantle upwelling flow. In the hot model, the contribution of the upwelling mantle flow to the lithosphere extension becomes efficient only in the advanced stages of the evolution, after the onset of the mantle partial melting. In contrast, for the cold model, both forces concur to the extension dynamics from the early stages of the evolution. These results agree with the models by Brune & Autin (Reference Brune and Autin2013) and Manatschal, Lavier & Chenin (Reference Manatschal, Lavier and Chenin2015) in which the break-up of a hotter and weaker lithosphere occurs later than in a colder and stronger lithosphere. In the case of a higher thermal state of the pre-rifting lithosphere, the viscous crustal layer is thicker than the brittle portion; consequently, the brittle strain softening is less efficient at focusing the deformation into discrete shear zones (Brune & Autin, Reference Brune and Autin2013). Lavier & Manatschal (Reference Lavier and Manatschal2006) and van Avendonk et al. (Reference Van Avendonk, Lavier, Shillington and Manatschal2009) suggest the opposite behaviour when the strong gabbroic lower crust is taken into account. In their models, a cold and strong lithosphere results in a longer rifting duration. On the other hand, the occurrence of a strong lower crust for the Alpine pre-rifting lithosphere is in contrast to the lithostratigraphy of the pre-Alpine continental crust. Given the number and size of Permian–Triassic gabbroic intrusions in the Alpine crust, the amount of detectable gabbroic rocks is less than 5 % of the total pre-Alpine lower continental crust actually exposed along the whole Alpine belt, which can be estimated from the tectonic map of the Alps (Bigi et al. Reference Bigi, Castellarin, Coli, Dal Piaz, Sartori, Scandone and Vai1990; Schmid et al. Reference Schmid, Fügenschuh, Kissling and Schuster2004) taking into account both units deeply involved in or escaping the Alpine subduction. Even considering a lower crust of the pristine passive margin that is richer in gabbroic rocks, it is reasonable to predict that, during the Alpine convergence, a selective tectonic sampling of the gabbro-poor lower crust does not occur. A coherent gabbroic lower crust for the Alpine pre-rifting lithosphere therefore seems unlikely.

For both chosen initial thermal configurations of the lithosphere, the exhumation of the serpentinized mantle starts before the oceanic spreading and the mantle partial melting (considered in our study coincident with gabbros formation), making the model compatible with the magma-poor rifting suggested for the Alpine case (e.g. Lavier & Manatschal, Reference Lavier and Manatschal2006; Pérez-Gussinyé et al. Reference Pérez-Gussinyé, Morgan, Reston and Ranero2006; Manatschal & Müntener, Reference Manatschal and Müntener2009; Manatschal, Lavier & Chenin Reference Manatschal, Lavier and Chenin2015), developing an ocean–continent transition zone similar to Galizia margin (e.g. Boillot, Girardeau & Kornprobst, Reference Boillot, Girardeau and Kornprobst1989; Manatschal, Reference Manatschal2004; Hébert et al. Reference Hébert, Beaudoin, Rochon and Gardien2008). The exhumation of a serpentinized lithospheric mantle before the oceanic spreading of the Liguria–Piemonte Ocean is also suggested based on the geochemical analysis of the syn-rift Alpine crust and sediments (Pinto et al. Reference Pinto, Manatschal, Karpoff and Viana2015).

For both the hot and cold models, the mantle partial melting does not appear immediately after the crustal break-up but after approximately 5 Ma in the hot model and after approximately 15 Ma in the cold model. This result is dependent on the thermal field predicted when the crustal break-up occurs. The hot model predicts more suitable thermo-barometric conditions for mantle melting a few million years after the crustal break-up. In contrast, in the cold model a long time is required to increase the thermal state to satisfy the pressure–temperature conditions that are suitable for mantle melting. Despite the different partial timings of each stage in the two models, once the serpentinization starts and the deformation localizes around the future ridge, the time required for the mantle partial melting and the beginning of oceanic spreading is comparable in the two models: 18 Ma for the cold model and 21 Ma for the hot model.

The continental crust thickness sensibly decreases during the extension but with different rates. In the hot model, the crustal thickness decreases from 30 km to approximately 18 km at the margin and to approximately 5 km close to the OCTZ within 31.4 Ma after the beginning of the evolution. In the cold model, it decreases from 30 km to approximately 22 km at the margin and to approximately 20 km close to the OCTZ within 21.4 Ma after the beginning of the evolution. Afterwards, in both models no further significant thinning occurs within the continental crust. The hyperextended margin envisaged by the hot model satisfies the model of an Alpine Tethys hyperextended system (Manatschal, Lavier & Chenin, Reference Manatschal, Lavier and Chenin2015). Our results indicate that in the proximity of the ridge, for both models, the extension rate at the base of the crust increases by a factor of 2–4 from the continental break-up to oceanization, in agreement with Whitmarsh, Manatschal & Minshull (Reference Whitmarsh, Manatschal and Minshull2001). Furthermore, at the advanced stages of the extension after the beginning of oceanic spreading, the extension rate may overcome the far-field traction.

A thinned continental crust (passive margin), an ocean–continent transition zone and an oceanic lithosphere characterize the final structure of the system from the periphery to the centre of the model domain for both thermal states. Our results show that the formation of the OCTZ starts 5–15 Ma before the partial melting of the mantle and develops with a size ranging over 160–280 km (according to the initial thermal configuration of the lithosphere), which is compatible with the observations of Galicia Margin (e.g. Boillot, Beslier & Girardeau, Reference Boillot, Beslier, Girardeau, Banda, Torné and Talwani1995; Hébert et al. Reference Hébert, Beaudoin, Rochon and Gardien2008) and similar to the interpretation for the Alpine rift (Manatschal & Müntener, Reference Manatschal and Müntener2009). This suggests that if the estimate of the oceanic basin width is based simply on the coincidence between the continental crust break-up and the onset of gabbros emplacement, as stated in several models proposed for the Alpine domain (e.g. Li et al. Reference Li, Faure, Lin and Manatschal2013), the effective basin width will be underestimated. In particular, Li et al. (Reference Li, Faure, Lin and Manatschal2013) estimated a basin width of 300 km after 10 Ma of extension at a full spreading rate of 3 cm a−1. In contrast, our model indicates that the extension of the basin would range over 360–480 km after the same time of 10 Ma for a full extension rate of 2.5 cm a−1.

The comparison between the natural data and the model predictions shows good agreement with all of the oceanic data for both the hot and cold models. In contrast, the comparison with the data from the continental crust lacks a complete fit because ages are not coincident. The lithological and thermal fits predicted by the hot and cold models with significant delays (from the Permian–Triassic to the Late Triassic – Late Jurassic periods; Figs 7, 9) suggest that, according to both models, the effects of a positive thermal anomaly should be recorded in the continental crust of both passive margins at 220–150 Ma. These effects have not yet been detected in the pre-Alpine continental crust of the Alps and Apennines. Similarly, the mantle partial melting occurs in both models 36.4 Ma and 22.4 Ma after the beginning of the extension in the hot and cold models, respectively, that would correspond to the early gabbros at 185 Ma, which is significantly younger than the Permian–Triassic continental gabbro emplacements. This time misfit supports the interpretation predicting a thermo-mechanical perturbation of the continental lithosphere in this portion of the Tethys due to the Variscan collision and the late orogenic extension (Spalla et al. Reference Spalla, Zanoni, Marotta, Rebay, Roda, Zucali, Gosso, Schulmann, Martínez Catalán, Lardeaux, Janoušek and Oggiano2014). In addition, the more favourable thermal regime predicted by the hot model could be due to a previously perturbed system either by the Variscan orogenic collapse or by an already thermally eroded and softened lithosphere. However, the sole orogenic collapse mechanism is not sufficient to reproduce the thermo-barometric conditions of the HT-LP Permian–Triassic metamorphism and intense igneous activity recorded in the continental lithosphere (Marotta & Spalla, Reference Marotta and Spalla2007; Marotta, Spalla & Gosso, Reference Marotta, Spalla, Gosso, Ring and Wernicke2009).

Finally, although the symmetry of the predicted thermal anomaly around the future ridge, HT natural parageneses from the Permian–Triassic continental lithosphere of the Alps are concentrated in the Austroalpine and Southalpine domains, while they are totally lacking in the Helvetic domain. This distribution of HT-LP metamorphic assemblages supports the interpretation of an asymmetric rifting (e.g. Lardeaux & Spalla, Reference Lardeaux and Spalla1991; Marotta, Spalla & Gosso, Reference Marotta, Spalla, Gosso, Ring and Wernicke2009).

7. Conclusions

We developed a 2D thermo-mechanical numerical model of passive rifting to investigate the evolution of the lithosphere of the Alps and the Northern Apennines during the transition from continental rifting to oceanic spreading of the Alpine Tethys. The model accounts for the crustal extension of a weak lower crust and mantle serpentinization, and results in symmetric rifting and denudation of the serpentinized lithospheric mantle.

A thinned continental crust (passive margin), an ocean–continent transition zone and an oceanic lithosphere characterize the final structure of the system. The thickness of the passive margin decreases over time from 30 km to 18 km (hot model) or 22 km (cold model) at the model boundaries and to 5 km (hot model, hyperextended margin configuration) or 20 km (cold model) close to the ocean–continent transition zone.

The mantle serpentinization starts before the crustal break-up, and the denudation occurs before the oceanic spreading. In addition, a hot and weak lithosphere evolves to oceanization slower than a cold and strong lithosphere, with a comparable time interval after the onset of serpentinization.

Our results indicate that, if the estimated basin width is based simply on the coincidence between the continental crust break-up and the onset of the gabbros emplacement, the effective width of the basin domain will be underestimated. The mantle denudation starts several million years before partial melting, generating an ocean–continent transition zone from the passive continental margin to the oceanic lithosphere with a size ranging over 160–280 km in a magma-poor rift pre-Alpine configuration.

The thermo-barometric predictions with their modelled timing were compared with the natural data derived from continental and oceanic rocks from the Alps and Apennines of Permian–Jurassic age, and the predictions from the hot model, which also promotes the development of hyperextended Alpine margins, agree with natural data better.

Our results support the idea that the Tethyan rifting should begin in a perturbed continental lithosphere, likely ascribable to the previous Variscan subduction-collision, as widely supported by the occurrence of HP metamorphic relics in the different Alpine structural domains. In fact, if rifting developed in a stable lithosphere, Triassic–Jurassic HT-LP metamorphism is predicted together with gabbro-basalt production younger than 185 Ma instead of the observed Permian–Triassic metamorphic and igneous records. Indeed, this has never been detected in the Alpine continental crust.

In addition, the distributions of the Permian–Triassic continental gabbros and the high-temperature metamorphism in the Austroalpine and Southalpine domains support the idea that it was asymmetric rifting in which the lithospheric signature of the Variscan subduction-collision can be a constraining inheritance for the successive geometry.

These ideas could be further confirmed by a new model that accounts for the previous history and the thermo-rheological consequences of Variscan Orogeny as initial configuration (as in Marotta, Spalla & Gosso, Reference Marotta, Spalla, Gosso, Ring and Wernicke2009) and evolves though continental break-up and the successive oceanization as in the present study.

Acknowledgements

The research was supported by PRIN 2011 (2010AZR98L) (Birth and death of oceanic basins: geodynamic processes from rifting to continental collision in Mediterranean and circum-Mediterranean orogens). The authors thank T. Gerya and U. Ring for fruitful suggestions. We also thank the editor and the two anonymous reviewers for their thorough review. All figures were created using GMT plotting software (Wessel & Smith, Reference Wessel and Smith2001) and Adobe Illustrator®.

References

Arcay, D., Tric, E. & Doin, M. P. 2005. Numerical simulation of subduction zones. Effect of slab dehydration on the mantle wedge dynamics. Physics of the Earth and Planetary Interiors 149, 133–53.Google Scholar
Baumgartner, P. O., Bartolini, A., Carter, E. S., Conti, M., Cortese, G., Danelian, T., De Wever, P., Dumitrica, P., Dumitrica-Jud, R., Gorican, S., Guex, J., Hull, D. M., Kito, N., Marcucci, M., Matsuoka, A., Murchey, B., O'Dogherty, L., Savary, J., Vishnevskaya, V., Widz, D. & Yao, A. 1995. Middle Jurassic to Early Cretaceous radiolarian biochronology of Tethys based on Unitary Associations. In Middle Jurassic to Lower Cretaceous Radiolaria of Tethys: Occurrences, Systematics, Biochronology (eds InterRad Jurassic-Cretaceous Working Group), pp. 23. Memoires de Geologie, Lausanne, Switzerland.Google Scholar
Beardsmore, G. R. & Cull, J. P. 2001. Crustal Heat Flow: A Guide to Measurement and Modelling. Cambridge: Cambridge University Press, 321 pp.Google Scholar
Beltrando, M., Rubatto, D. & Manatschal, G. 2010. From passive margins to orogens: The link between ocean-continent transition zones and (ultra) high-pressure metamorphism. Geology 38, 559–62.Google Scholar
Benciolini, L., Poli, M. E., Visona, D. & Zanferrari, A. 2006. Looking inside Late Variscan tectonics: structural and metamorphic heterogeneity of the Eastern Southalpine basement (NE Italy). Geodinamica Acta 19, 1732.CrossRefGoogle Scholar
Bertotti, G., Picotti, V., Bernoulli, D. & Castellarin, A. 1993. From rifting to drifting: tectonic evolution of the South-Alpine upper crust from the Triassic to the Early Cretaceous. Sedimentary Geology 86 (1–2), 5376.Google Scholar
Best, M. G. & Christiansen, E. H. 2001. Igneous Petrology. London: Blackwell Science, 458 pp.Google Scholar
Biagini, L., Bistacchi, A., Gosso, G., Magistroni, C., Rossetti, I., Spalla, M. I. & Tognoni, A. 1995. The II DK, HT mega-relic in the Sesia-Lanzo Zone: Late Variscan collision or Permo-Triassic rifting? In IOS International Ophiolite Symposium, Pavia, 1823 September 1995, pp. 22.Google Scholar
Bigi, G., Castellarin, A., Coli, M., Dal Piaz, G.V., Sartori, R., Scandone, P. & Vai, G.B. 1990. Structural Model of Italy, sheets 1-2: CNR, Progetto Finalizzato Geodinamica.Google Scholar
Bill, M., O'Dogherty, L., Guex, J., Baumgartner, P. O. & Masson, H. 2001. Radiolarite ages in Alpine-Mediterranean ophiolites: Constraints on the oceanic spreading and the Tethys-Atlantic connection. Geological Society of America Bulletin 113 (1), 129–43.Google Scholar
Bocquet, J., Delaloye, M., Hunziker, J. C. & Krummenacher, D. 1974. K-Ar and Rb-Sr dating of blue amphiboles, micas, and associated minerals from the Western Alps. Contributions to Mineralogy and Petrology 47 (1), 726.Google Scholar
Boillot, G., Beslier, M. O. & Girardeau, J. 1995. Nature, structure and evolution of the ocean-continent boundary: the lesson of the west Galicia margin (Spain). In Rifted Ocean-Continent Boundaries (eds. Banda, E., Torné, M. & Talwani, M.), pp. 219–29. Dordrecht, Netherlands: Springer.Google Scholar
Boillot, G., Girardeau, J. & Kornprobst, J. 1989. Rifting of the Galicia Margin: crustal thinning and emplacement of mantle rocks on the seafloor. In Proceedings of the Ocean Drilling Program, pp. 741–56. College Station, Texas, Scientific Results no. 103.Google Scholar
Bonin, B., Brändlein, P., Bussy, F., Desmons, J., Eggenberger, U., Finger, F., Graf, K., Marro, C., Mercolli, I., Oberhänsli, R., Ploquin, A., Quadt, A. von, Raumer, J. von, Schaltegger, U., Steyrer, H. P. & Visonà, D. 1993. In Pre-Mesozoic Geology in the Alps (eds. von Raumer, J. F. & Neubauer, F.), pp. 327–44. Berlin, Heidelberg: Springer.Google Scholar
Boriani, A. & Burlini, L. 1995. Carta Geologica della valle Cannobina Scala 1:25000. Milano: Dipartimento di Scienze della Terra “Ardito Desio” Università degli Studi di Milano.Google Scholar
Boriani, A., Colombo, A. & Macera, P. 1985. Radiometric geochronology of Central Alps. Rendiconti della Società Italiana di Mineralogia e Petrologia 40, 139–86.Google Scholar
Borsi, S., Ferrara, G., Paganelli, L. & Simboli, G. 1968. Isotopic age measurements of the M.Monzoni intrusive complex. Mineralogica et Petrographica Acta 14, 171–83.Google Scholar
Bouffette, J., Lardeaux, J. M. & Caron, J. M. 1993. Le passage des granulites aux éclogites dans les méta-pélites de l'unité de la Punta Muret (Massif Dora-Maira, Alpes occidentales). Comptes Rendus de l'Académie des Sciences 317, 1617–24.Google Scholar
Bousquet, R., Engi, M., Gosso, G., Oberhänsli, R., Berger, A., Spalla, M. I., Zucali, M. & Goffè, B. 2004. Explanatory notes to the map: metamorphic structure of the Alps transition from the Western to the Central Alps. Mitteilungen der Österreichischen Mineralogischen Gesellschaft 149, 145–56.Google Scholar
Bozzo, E., Campi, S., Capponi, G. & Giglia, G. 1992. The suture between the Alps and Apennines in the Ligurian sector based on geological and geomagnetic data. Tectonophysics 206 (1–2), 159–69.CrossRefGoogle Scholar
Brodie, K. H., Rex, D. & Rutter, E. H. 1989. On the age of deep crustal extensional faulting in the Ivrea zone, Northern Italy. In Alpine Tectonics (eds Coward, R. G., Dietrich, M. P. & Park, D.), pp. 203–10. Geological Society, London, Special Publication no. 45.Google Scholar
Brune, S. 2014. Evolution of stress and fault patterns in oblique rift systems: 3-D numerical lithospheric-scale experiments from rift to breakup. Geochemistry, Geophysics, Geosystems 15 (8), 3392–415.Google Scholar
Brune, S. & Autin, J. 2013. The rift to break-up evolution of the Gulf of Aden: Insights from 3D numerical lithospheric-scale modelling. Tectonophysics 607, 6579.Google Scholar
Brune, S., Heine, C., Pérez-Gussinyé, M. & Sobolev, S. V. 2014. Rift migration explains continental margin asymmetry and crustal hyper-extension. Nature Communications 5, article no. 4014.Google Scholar
Buck, W. R. 1991. Modes of continental lithospheric extension. Journal of Geophysical Research: Solid Earth 96 (B12), 20161–78.CrossRefGoogle Scholar
Buergi, A. & Kloetzli, U. 1990. New data on the evolutionary history of the Ivrea Zone (Northern Italy). Bulletin of the Swiss Association of Petroleum Geology and Engineering 56 (130), 4970.Google Scholar
Bussy, F., Venturini, G., Hunziker, J. & Martinotti, G. 1998. U-Pb ages of magmatic rocks of the western Austroalpine Dent-Blanche-Sesia Unit. Schweizerische Mineralogische Und Petrographische Mitteilungen 78, 163–8.Google Scholar
Carminati, E. & Doglioni, C. 2012. Alps vs. Apennines: The paradigm of a tectonically asymmetric Earth. Earth-Science Reviews 112 (1–2), 6796.Google Scholar
Cassinis, R. 2006. Reviewing pre-TRANSALP DFF models. Tectonophysics 414, 7986.Google Scholar
Chalot-Prat, F. 2005. An undeformed ophiolite in the Alps: Field and geochemical evidence for a link between volcanism and shallow plate tectonic processes. In Plates, Plumes and Paradigms (eds Foulger, G. R., Natland, J. H., Presnall, D. C. & Anderson, D. L.), pp. 750–80. Geological Society of America, Special Paper no. 388.Google Scholar
Chopra, P. N. & Peterson, M. S. 1981. The experimental deformation of dunite. Tectonophysics 78, 453–73.CrossRefGoogle Scholar
Clague, D. A. 1987. Hawaiian xenolith populations, magma supply rates, and development of magma chambers. Bulletin of Volcanology 49 (4), 577–87.Google Scholar
Cloetingh, S., Burov, E., Matenco, L., Beekman, F., Roure, F. & Ziegler, P. A. 2013. The Moho in extensional tectonic settings: Insights from thermo-mechanical models. Tectonophysics 609, 558604.Google Scholar
Colombo, A. & Tunesi, A. 1999. Pre-Alpine metamorphism of the southern Alps. Schweizerische Mineralogische und Petrographische Mitteilungen 79, 6377.Google Scholar
Compagnoni, R., di Brozolo, F. R. & Sandrone, R. 1984. Kaersutite-bearing mylonitic gabbro from the Lanzo-peridotite (western Italian Alps). Geologie en Mijnbouw 63 (2), 189–96.Google Scholar
Cordey, F. & Bailly, A. 2007. Alpine ocean seafloor spreading and onset of pelagic sedimentation: new radiolarian data from the Chenaillet–Montgenèvre ophiolite (French–Italian Alps). Geodinamica Acta, 20, 131–8.Google Scholar
Corti, G., Bonini, M., Sokoutis, D., Innocenti, F., Manetti, P., Cloetingh, S. & Mulugeta, G. 2004. Continental rift architecture and patterns of magma migration: A dynamic analysis based on centrifuge models. Tectonics 23 (2), TC2012.Google Scholar
Costa, S. & Caby, R. 2001. Evolution of the Ligurian Tethys in the Western Alps: Sm/Nd and U–Pb geochronology and rare-earth element geochemistry of the Montgenevre ophiolite (France). Chemical Geology 3–4 (175), 449–66.Google Scholar
Dal Piaz, G. V. 1993. Evolution of Austro-Alpine and Upper Penninic Basement in the Northwestern Alps from Variscan Convergence to Post-Variscan Extension. In Pre-Mesozoic Geology in the Alps (eds Raumer, J. F. & Neubauer, F.), pp. 327–44. Berlin, Heidelberg: Springer.Google Scholar
Dal Piaz, G. V. 2010. The Italian Alps: a journey across two centuries of Alpine geology. In The Geology of Italy: Tectonics and Life along Plate Margins (eds Beltrando, M., Peccerillo, A., Mattei, M., Conticelli, S. & Doglioni, C.), pp. 1108. Conder, Australia: Journal of the Virtual Explorer, Electronic Edition, 36.Google Scholar
Dal Piaz, G. V., De Vecchi, G. P. & Hunziker, J. C. 1977. The Austroalpine layered gabbros of the Matterhorn and Mt. Collon-Dents de Bertol. Schweizerische Mineralogische und Petrographische Mitteilungen 57, 5988.Google Scholar
Dal Piaz, G. V., Lombardo, B. & Gosso, G. 1983. Metamorphic evolution of the Mt. Emilius klippe, Dent Blanche nappe, Western Alps. American Journal of Science 283A, 438–58.Google Scholar
De Capitani, L., Carnevale, M. & Fumagalli, M. 2007. Gamma-ray spectroscopy determination of radioactive elements in late-Hercynian plutonic rocks of Val Biandino and Val Trompia (Lombardy, Italy). Periodico di Mineralogia 76 (1), 2539.Google Scholar
Del Moro, A. & Visonà, D. 1982. The epiplutonic Hercynian Complex of Bressanone (Brixen, Eastern Alps, Italy). Petrologic and radiometric data. Neues Jahrbuch fur Mineralogie - Abhandlungen 145, 6685.Google Scholar
Desmons, J. 1992. The Briançon basement (Pennine Western Alps): mineral composition and polymetamorphic evolution. Schweizerische Mineralogische und Petrographische Mitteilungen 72, 3755.Google Scholar
Desmurs, L., Manatschal, G. & Bernoulli, D. 2001. The Steinmann Trinity revisited: mantle exhumation and magmatism along an ocean-continent transition: the Platta nappe, eastern Switzerland. In Non-Volcanic Rifting of Continental Margins: A Comparison of Evidence from Land and Sea (eds. Wilson, R.C.L., Whitmarsh, R.B., Taylor, B. & Froitzheim, N.), pp. 235–66. Geological Society, London, Special Publication no. 187.Google Scholar
di Paola, S. & Spalla, M. I. 2000. Contrasting tectonic records in pre-Alpine metabasites of the Southern Alps (lake Como, Italy). Journal of Geodynamics 30 (1–2), 167–89.Google Scholar
Diella, V., Spalla, M. I. & Tunesi, A. 1992. Contrasted thermo-mechanical evolutions in the South Alpine metamorphic basement of the Orobic Alps (Central Alps, Italy). Journal of Metamorphic Geology 10, 203–19.CrossRefGoogle Scholar
Donatio, D., Marroni, M. & Rocchi, S. 2013. Serpentinization history in mantle section from a fossil slow-spreading ridge sequence/evidences from Pomaia quarry (Southern Tuscany, Italy). Ofioliti 38 (1), 1528.Google Scholar
Dubois, J. & Diament, M. 1997. Géophysique. Paris: Masson, 205 pp.Google Scholar
England, P. C. & Thompson, A. B. 1984. Pressure-temperature-time paths of regional metamorphism I. Heat transfer during the evolution of regions of thickened continental crust. Journal of Petrology 25 (4), 894928.Google Scholar
Escartín, J., Hirth, G. & Evans, B. 1997. Effects of serpentinization on the lithospheric strength and the style of normal faulting at slow-spreading ridges. Earth and Planetary Science Letters 151 (3–4), 181–9.Google Scholar
Ewing, T., Hermann, J. & Rubatto, D. 2013. The robustness of the Zr-in-rutile and Ti-in-zircon thermometers during high-temperature metamorphism (Ivrea-Verbano Zone, northern Italy). Contributions to Mineralogy and Petrology 165 (4), 757–79.CrossRefGoogle Scholar
Federico, L., Spagnolo, C., Crispini, L. & Capponi, G. 2009. Fault-slip analysis in the metaophiolites of the Voltri Massif: constraints for the tectonic evolution at the Alps/Apennine boundary. Geological Journal 44 (2), 225–40.Google Scholar
Ferry, J. M., Wing, B. A., Penniston-Dorland, S. C. & Rumble, D. 2002. The direction of fluid flow during contact metamorphism of siliceous carbonate rocks: new data for the Monzoni and Predazzo aureoles, northern Italy, and a global review. Contributions to Mineralogy and Petrology 142 (6), 679–99.Google Scholar
Fontana, E., Panseri, M. & Tartarotti, P. 2008. Oceanic relict textures in the Mount Avic Serpentinites, Western Alps. Ofioliti 33 (2), 105–18.Google Scholar
Gaidies, F., Abart, R., De Capitani, C., Schuster, R., Connolly, J. A. D. & Reusser, E. 2006. Characterization of polymetamorphism in the Austroalpine basement east of the Tauern Window using garnet isopleth thermobarometry. Journal of Metamorphic Geology 24 (6), 451–75.CrossRefGoogle Scholar
Galli, A., Le Bayon, B., Schmidt, M. W., Burg, J.-P., Caddick, M. J. & Reusser, E. 2011. Granulites and charnockites of the Gruf Complex: Evidence for Permian ultra-high temperature metamorphism in the Central Alps. Lithos 124 (1–2), 1745.Google Scholar
Galli, A., Le Bayon, B., Schmidt, M. W., Burg, J.-P., Reusser, E., Sergeev, S. A. & Larionov, A. 2012. U–Pb zircon dating of the Gruf Complex: disclosing the late Variscan granulitic lower crust of Europe stranded in the Central Alps. Contributions to Mineralogy and Petrology 163 (2), 353–78.Google Scholar
Gallien, F., Abart, R. & Wyhlidal, S. 2007. Contact metamorphism and selective metasomatism of the layered Bellerophon Formation in the eastern Monzoni contact aureole, northern Italy. Mineralogy and Petrology 91, 2553.Google Scholar
Gardien, V., Reusser, E. & Marquer, D. 1994. Pre-Alpine metamorphic evolution of the gneisses from the Valpelline series (Western Alps, Italy). Schweizerische Mineralogische und Petrographische Mitteilungen 74, 489502.Google Scholar
Garuti, G., Bea, F., Zaccarini, F. & Montero, P. 2001. Age, geochemistry and petrogenesis of the ultramafic pipes in the Ivrea zone, NW Italy. Journal of Petrology 42 (2), 433–57.Google Scholar
Gerya, T. V. & Stöckhert, B. 2005. Two-dimensional numerical modeling of tectonic and metamorphic histories at active continental margins. International Journal of Earth Sciences 95 (2), 250–74.Google Scholar
Giacomini, F., Messiga, B., Tribuzio, R. & Braga, R. 1999. The Sondalo gabbroic complex and its country rocks: new geological and petrological data. Tübinger Geowissenschaftliche Arbeiten. Reihe A, Geologie, Palaeontologie, Stratigraphie, vol. 52, pp. 156.Google Scholar
Gillcrist, R., Coward, M. & Mugnier, J. L. 1987. Structural inversion and its controls: examples from the Alpine foreland and the French-Alps. Geodinamica Acta 1, 534.Google Scholar
Goffé, B., Schwartz, S., Lardeaux, J. M. & Bousquet, R. 2004. Explanatory notes to the map: metamorphic structure of the Alps Western and Ligurian Alps. Mitteilungen der Österreichischen Mineralogischen Gesellschaft 149, 125–44.Google Scholar
Gosso, G., Siletto, G. & Spalla, M. I. 1997. H-T/L-P metamorphism and structures in the South-Alpine basement near Lake Como, Orobic Alps; intracontinental imprints of the Permo-Triassic rifting. Ofioliti 22, 133–45.Google Scholar
Gregnanin, A. 1980. Metamorphism and magmatism in the western Italian Tyrol. Rivista Italiana di Mineralogia e Petrologia 36, 4964.Google Scholar
Haas, R. 1985. Zur Metamorphose des Suedlichen Oetz-talkristallins unter Besonderer Beruecksichtigung der Matscher Einheit (Vintschgau/Suedtirol). Ph.D. thesis, University of Innsbruck. Published thesis.Google Scholar
Habler, G. & Thöni, M. 2001. Preservation of Permo-Triassic low-pressure assemblages in the Cretaceous high-pressure metamorphic Saualpe crystalline basement (Eastern Alps, Austria). Journal of Metamorphic Geology 19, 679–97.Google Scholar
Handy, M. R., Franz, L., Heller, F., Janott, B. & Zurbriggen, R. 1999. Multistage accretion and exhumation of the continental crust (Ivrea crustal section, Italy and Switzerland). Tectonics 18 (6), 1154–77.Google Scholar
Handy, M. R. & Oberhänsli, R. 2004. Explanatory notes to the map: metamorphic structure of the Alps age map of the metamorphic structure of the Alps: tectonic interpretation and outstanding problem. Mitteilungen der Österreichischen Mineralogischen Gesellschaft 149, 201–25.Google Scholar
Handy, M. R., Schmid, S. M., Bousquet, R., Kissling, E. & Bernoulli, D. 2010. Reconciling plate-tectonic reconstructions of Alpine Tethys with the geological-geophysical record of spreading and subduction in the Alps. Earth-Science Reviews 102, 121–58.Google Scholar
Handy, M. R. & Zingg, A. 1991. The tectonic and rheological evolution of an attenuated cross section of the continental crust: Ivrea crustal section, southern Alps, northwestern Italy and southern Switzerland. Geological Society of America Bulletin 103 (2), 236–53.Google Scholar
Hansmann, W., Muntener, O. & Hermann, J. 2001. U-Pb zircon geochronology of a tholeiitic intrusion and associated migmatites at a continental crust-mantle transition, Val Malenco, Italy. Schweizerische Mineralogische und Petrographische Mitteilungen 81 (1), 239–56.Google Scholar
Hébert, R., Beaudoin, G., Rochon, M. & Gardien, V. 2008. Metamorphic evolution and oxygen isotope geochemistry of rift-precursor amphibolites from Hole 1067A ODP Leg 173 off West Iberian Galicia Bank rifted margin. Lithos 101, 162–76.Google Scholar
Henk, A., Franz, L., Teufel, S. & Oncken, O. 1997. Magmatic underplating, extension, and crustal reequilibration: insights from a cross-section through the Ivrea Zone and Strona-Ceneri Zone, northern Italy. Journal of Geology 105 (3), 367–77.Google Scholar
Hermann, J. & Rubatto, D. 2003. Relating zircon and monazite domains to garnet growth zones: age and duration of granulite facies metamorphism in the Val Malenco lower crust. Journal of Metamorphic Geology 21 (9), 833–52.Google Scholar
Hoke, L. 1990. The Altkristallin of the Kreuzeck Mountains, SE Tauern Window, Eastern Alps: Basement crust in a convergent plate boundary zone. Jahrbuch des Geologischen Bundesantall 133 (1), 587.Google Scholar
Honda, S. & Saito, M. 2003. Small-scale convection under the back-arc occurring in the low viscosity wedge. Earth and Planetary Science Letters 216, 703–15.Google Scholar
Huismans, R. S. & Beaumont, C. 2011. Depth-dependent extension, two-stage breakup and cratonic underplating at rifted margins. Nature 473 (7345), 7478.Google Scholar
Huismans, R. S. & Beaumont, C. 2014. Rifted continental margins: The case for depth-dependent extension. Earth and Planetary Science Letters 407, 148–62.Google Scholar
Huismans, R. S., Buiter, S. J. H. & Beaumont, C. 2005. Effect of plastic-viscous layering and strain softening on mode selection during lithospheric extension. Journal of Geophysical Research: Solid Earth 110 (B2), 117.CrossRefGoogle Scholar
Hunziker, J. C., Desmon, J. & Hurford, A. J. 1992. Thirty-two years of geochronological work in the Central and Western Alps: a review of seven maps. Mémoires de Géologie, Lausanne 13, 159.Google Scholar
Hunziker, J. C. & Zingg, A. 1980. Lower Paleozoic amphibolite to granulite facies metamorphism in the Ivrea Zone (southern Alps, northern Italy). Schweizerische Mineralogische und Petrographische Mitteilungen 60, 181213.Google Scholar
Kaczmarek, M.-A., Müntener, O. & Rubatto, D. 2008. Trace element chemistry and U-Pb dating of zircons from oceanic gabbros and their relationship with whole rock composition (Lanzo, Italian Alps). Contributions to Mineralogy and Petrology 155 (3), 295312.Google Scholar
Klötzli, U. S., Sinigoi, S., Quick, J. E., Demarchi, G., Tassinari, C. C. G., Sato, K. & Günes, Z. 2014. Duration of igneous activity in the Sesia Magmatic System and implications for high-temperature metamorphism in the Ivrea-Verbano deep crust. Lithos 206, 1933.CrossRefGoogle Scholar
Lagabrielle, Y. & Cannat, M. 1990. Alpine Jurassic ophiolites resemble the modern central Atlantic basement. Geology 18 (4), 319–22.Google Scholar
Lagabrielle, Y., Fudral, S. & Kienast, J.-R. 1989. La couverture océanique des ultrabasites de Lanzo (Alpes occidentales): arguments lithostratigraphiques et pétrologiques. Geodinamica Acta 4 (1), 4355.Google Scholar
Langone, A. & Tiepolo, M. 2015. U-Th-Pb ‘multi-phase’ approach to the study of crystalline basement: application to the northernmost sector of the Ivrea-Verbano Zone (Alps). Periodico di Mineralogia 84 (3B), 633–55.Google Scholar
Lardeaux, J. M. 1981. Evolution Tectono-metamorphique de la Zone Nord du Massif de Sesia-Lanzo (Alpes Occidentales): Un Exemple d’Éclogitisation de Croute Continentale. Ph.D. thesis, University of Paris VI. Published thesis.Google Scholar
Lardeaux, J.-M. 2014. Deciphering orogeny: a metamorphic perspective. Examples from European Alpine and Variscan belts: Part I: Alpine metamorphism in the western Alps. A review. Bulletin de la Societe Geologique de France 185 (2), 93114.Google Scholar
Lardeaux, J. M. & Spalla, M. I. 1991. From granulites to eclogites in the Sesia zone (Italian Western Alps): a record of the opening and closure of the Piedmont Ocean. Journal of Metamorphic Geology 9 (1), 3559.Google Scholar
Lavier, L. L. & Manatschal, G. 2006. A mechanism to thin the continental lithosphere at magma-poor margins. Nature 440 (7082), 324–8.CrossRefGoogle ScholarPubMed
Li, X.-H., Faure, M., Lin, W. & Manatschal, G. 2013. New isotopic constraints on age and magma genesis of an embryonic oceanic crust: the Chenaillet Ophiolite in the Western Alps. Lithos 160–161, 283–91.Google Scholar
Li, X.-H., Faure, M., Rossi, P., Lin, W. & Lahondère, D. 2015. Age of Alpine Corsica ophiolites revisited: Insights from in situ zircon U-Pb age and O-Hf isotopes. Lithos 220–223, 179–90.Google Scholar
Liao, J. & Gerya, T. 2015. From continental rifting to seafloor spreading: Insight from 3D thermo-mechanical modeling. Gondwana Research 28 (4), 1329–43.Google Scholar
Lu, M., Hofmann, A. W., Mazzucchelli, M. & Rivalenti, G. 1997. The mafic-ultramafic complex near Finero (Ivrea-Verbano Zone), II. Geochronology and isotope geochemistry. Chemical Geology 140 (3–4), 223–35.Google Scholar
Malatesta, C., Crispini, L., Federico, L., Capponi, G. & Scambelluri, M. 2012. The exhumation of high pressure ophiolites (Voltri Massif, Western Alps): Insights from structural and petrologic data on metagabbro bodies. Tectonophysics 568–569, 102–23.Google Scholar
Malatesta, C., Gerya, T., Crispini, L., Federico, L. & Capponi, G. 2013. Oblique subduction modelling indicates along-trench tectonic transport of sediments. Nature Communications 4, 16.Google Scholar
Manatschal, G. 2004. New models for evolution of magma-poor rifted margins based on a review of data and concepts from West Iberia and the Alps. International Journal of Earth Sciences 93 (3), 432–66.Google Scholar
Manatschal, G., Lavier, L. & Chenin, P. 2015. The role of inheritance in structuring hyperextended rift systems: Some considerations based on observations and numerical modeling. Gondwana Research 27 (1), 140–64.Google Scholar
Manatschal, G. & Müntener, O. 2009. A type sequence across an ancient magma-poor ocean-continent transition: the example of the western Alpine Tethys ophiolites. Tectonophysics 473 (1–2), 419.Google Scholar
Manatschal, G., Sauter, D., Karpoff, A. M., Masini, E., Mohn, G. & Lagabrielle, Y. 2011. The Chenaillet Ophiolite in the French/Italian Alps: An ancient analogue for an Oceanic Core Complex? Lithos 124 (3–4), 169–84.Google Scholar
Manzotti, P. & Zucali, M. 2013. The pre-Alpine tectonic history of the Austroalpine continental basement in the Valpelline unit (Western Italian Alps). Geological Magazine 150 (1), 153–72.Google Scholar
Marotta, A. M. & Spalla, M. I. 2007. Permian-Triassic high thermal regime in the Alps: Result of late Variscan collapse or continental rifting? Validation by numerical modeling. Tectonics 26, 127.Google Scholar
Marotta, A. M., Spalla, M. I. & Gosso, G. 2009. Upper and lower crustal evolution during lithospheric extension: numerical modelling and natural footprints from the European Alps. In Extending a Continent: Architecture, Rheology and Heat Budget (eds Ring, U. & Wernicke, B.), pp. 3372. The Geological Society, London, Special Publication no. 321.Google Scholar
Marotta, A. M., Spelta, E. & Rizzetto, C. 2006. Gravity signature of crustal subduction inferred from numerical modelling. Geophysical Journal International 166, 923–38.Google Scholar
Marroni, M., Molli, G., Montanini, A. & Tribuzio, R. 1998. The association of continental crust rocks with ophiolites in the Northern Apennines (Italy): implications for the continent-ocean transition in the Western Tethys. Tectonophysics 292 (1–2), 4366.Google Scholar
Marroni, M. & Pandolfi, L. 2007. The architecture of an incipient oceanic basin: a tentative reconstruction of the Jurassic Liguria-Piemonte basin along the Northern Apennines-Alpine Corsica transect. International Journal of Earth Sciences 96 (6), 1059–78.Google Scholar
Marroni, M. & Tribuzio, R. 1996. Gabbro-derived granulites from external liguride units (northern Apennine, Italy): implications for the rifting processes in the western Tethys. Geologische Rundschau 85 (2), 239–49.CrossRefGoogle Scholar
Martin, S., Tartarotti, P. & Dal Piaz Giorgio, V. 1994. The Mesozoic ophiolites of the Alps: a review. OGS/Bollettino di Geofisica Teorica e Applicata 36 (141–144), 175220.Google Scholar
Mayer, A., Mezger, K. & Sinigoi, S. 2000. New Sm–Nd ages for the Ivrea–Verbano Zone, Sesia and Sessera valleys (Northern Italy). Journal of Geodynamics 30 (1–2), 147–66.Google Scholar
Mevel, C., Caby, R. & Kienast, J.-R. 1978. Amphibolite facies conditions in the oceanic crust: example of amphibolitized flaser-gabbro and amphibolites from the Chenaillet ophiolite massif (Hautes Alpes, France). Earth and Planetary Science Letters 39 (1), 98108.Google Scholar
Miller, C. & Thöni, M. 1997. Eo-alpine eclogitisation of Permian MORB-type gabbros in the Koralpe (Eastern Alps, Austria): new geochronological, geochemical and petrological data. Chemical Geology 137 (3–4), 283310.Google Scholar
Miller, C., Thöni, M., Goessler, W. & Tessadri, R. 2011. Origin and age of the Eisenkappel gabbro to granite suite (Carinthia, SE Austrian Alps). Lithos 125 (1–2), 434–48.Google Scholar
Mohn, G., Manatschal, G., Beltrando, M., Masini, E. & Kusznir, N. 2012. Necking of continental crust in magma-poor rifted margins: Evidence from the fossil Alpine Tethys margins. Tectonics 31 (1), TC1012.Google Scholar
Monjoie, P. 2004. The Mont Collon mafic complex (Austroalpine Dent Blanche nappe): Permian evolution of the Western European mantle. Ph.D. thesis. Université de Lausanne. Published thesis.Google Scholar
Monjoie, P., Bussy, F., Lapierre, H. & Pfeifer, H.-R. 2005. Modeling of in-situ crystallization processes in the Permian mafic layered intrusion of Mont Collon (Dent Blanche nappe, western Alps). Lithos 83 (3–4), 317–46.Google Scholar
Montanini, A., Tribuzio, R. & Anczkiewicz, R. 2006. Exhumation history of a garnet pyroxenite-bearing mantle section from a continent-ocean transition (northern Apennine Ophiolites, Italy). Journal of Petrology 47 (10), 1943–71.CrossRefGoogle Scholar
Mundil, R., Brack, P. & Laurenzi, M. A. 1996. High resolution U–Pb single-zircon age determinations: new constraints on the timing of Middle Triassic magmatism in the Southern Alps. In 78a Riunione Estiva, Geologia delle Dolomiti, pp. 1, Società Geologica Italiana. San Cassiano, 16–18 September 1996.Google Scholar
Muntener, O. & Hermann, J. 2001. The role of lower crust and continental upper mantle during formation of non-volcanic passive margins: evidence from the Alps. In Non-Volcanic Rifting of Continental Margins: A Comparison of Evidence from Land and Sea (eds. Wilson, R.C.L., Whitmarsh, R.B., Taylor, B. & Froitzheim, N.), pp. 267–88. Geological Society, London, Special Publication no. 187.Google Scholar
Muntener, O., Hermann, J. & Trommsdorf, V. 2000. Cooling history and exhumation of lower-crustal granulite and upper mantle (Malenco, eastern Central Alps). Journal of Petrology 41 (2), 175200.Google Scholar
Nagel, T. J. & Buck, W. R. 2004. Symmetric alternative to asymmetric rifting models. Geology 32 (11), 937–40.Google Scholar
Nagy, G., Draganits, E., Demény, A., Pantó, G. & Árkai, P. 2002. Genesis and transformations of monazite, florencite and rhabdophane during medium grade metamorphism: examples from the Sopron Hills, Eastern Alps. Chemical Geology 191 (1–3), 2546.Google Scholar
Naliboff, J. & Buiter, S. J. H. 2015. Rift reactivation and migration during multiphase extension. Earth and Planetary Science Letters 421, 5867.Google Scholar
Nicot, E. 1977. Les roches meso and catazonales de la Valpelline (nappe de la Dent Blanche; Alpes Italiennes). Ph.D. thesis. Université Pierre et Marie Curie, Paris VI. Published thesis.Google Scholar
Ohnenstetter, M., Ohnenstetter, D., Vidal, P., Cornichet, J., Hermitte, D. & Mace, J. 1981. Crystallization and age of zircon from Corsican ophiolitic albitites: consequences for oceanic expansion in Jurassic times. Earth and Planetary Science Letters 54 (3), 397408.Google Scholar
Pennacchioni, G. & Cesare, B. 1997. Ductile-brittle transition in pre-Alpine amphibolite facies mylonites during evolution from water-present to water-deficient conditions (Mont Mary nappe, Italian Western Alps). Journal of Metamorphic Geology 15 (6), 777–91.Google Scholar
Peressini, G., Quick, J. E., Sinigoi, S., Hofmann, A. W. & Fanning, M. 2007. Duration of a large mafic intrusion and heat transfer in the lower crust: a SHRIMP U-Pb zircon study in the Ivrea-Verbano Zone (western Alps, Italy). Journal of Petrology 48 (6), 1185–218.Google Scholar
Pérez-Gussinyé, M., Morgan, J. P., Reston, T. J. & Ranero, C. R. 2006. The rift to drift transition at non-volcanic margins: Insights from numerical modelling. Earth and Planetary Science Letters 244 (1–2), 458–73.Google Scholar
Pérez-Gussinyé, M., Reston, T. J. & Phipps Morgan, J. 2001. Serpentinization and magmatism during extension at non-volcanic margins: the effect of initial lithospheric structure. In Non-Volcanic Rifting of Continental Margins: A Comparison of Evidence from Land and Sea (eds. Wilson, R.C.L., Whitmarsh, R.B., Taylor, B. & Froitzheim, N.), pp. 551–76. Geological Society, London, Special Publication no. 187.Google Scholar
Piccardo, G. B. & Guarnieri, L. 2010. Alpine peridotites from the Ligurian Tethys: an updated critical review. International Geology Review 52 (10–12), 1138–59.Google Scholar
Piccardo, G. B., Padovano, M. & Guarnieri, L. 2014. The Ligurian Tethys: Mantle processes and geodynamics. Earth-Science Reviews 138, 409–34.Google Scholar
Pin, C. 1986. Datation U–Pb sur zircon à 285 Ma du complexe gabbro dioritique du Val Sesia - Val Mastallone et age tardi hercynien du métamorphisme granulitique de la zone Ivrea-Verbano (Italie). Compte Rendu Academie des Sciences de Paris 303, 827–30.Google Scholar
Pinarelli, L. & Boriani, A. 2007. Tracing metamorphism, magmatism and tectonics in the southern Alps (Italy): constraints from Rb-Sr and Pb-Pb geochronology, and isotope geochemistry. Periodico di Mineralogia 76, 524.Google Scholar
Pinto, V. H. G., Manatschal, G., Karpoff, A. M. & Viana, A. 2015. Tracing mantle-reacted fluids in magma-poor rifted margins: The example of Alpine Tethyan rifted margins. Geochemistry, Geophysics, Geosystems 16 (9), 3271–308.Google Scholar
Platt, J. P. 1986. Dynamic of orogenic wedges and the uplift of high-pressure metamorphic rocks. Geological Society of America Bulletin 97, 10371053.Google Scholar
Pognante, U., Rösli, U. & Toscani, L. 1985. Petrology of ultramafic and mafic rocks from the Lanzo peridotite body (Western Alps). Lithos 18, 201–14.Google Scholar
Polino, R., Dal Piaz, G. V. & Gosso, G. 1990. Tectonic erosion at the Adria margin and accretionary processes for the Cretaceous orogeny of the Alps. Mémoires de la Societé Géologique de France 156, 345–67.Google Scholar
Povoden, E., Horacek, M. & Abart, 2002. Contact metamorphism of siliceous dolomite and impure limestones from the Werfen formation in the eastern Monzoni contact aureole. Mineralogy and Petrology 76, 99120.Google Scholar
Quick, J. E., Sinigoi, S., Negrini, L., Demarchi, G. & Mayer, A. 1992. Synmagmatic deformation in the underplated igneous complex of the Ivrea-Verbano zone. Geology 20 (7), 613–6.Google Scholar
Quick, J. E., Sinigoi, S., Snoke, A. W., Kalakay, T. J., Mayer, A. & Peressini, G. 2002. Geologic map of the southern Ivrea-Verbano zone, northwestern Italy. In Geologic Investigation Series Map I-2776 and booklet, pp. 22. Reston, Virginia: US Geological Survey, US Government Printing Office.Google Scholar
Rampone, E. 2002. Mantle dynamics during Permo-Mesozoic extension of the Europe-Adria lithosphere: insights from the Ligurian ophiolites. Periodico di Mineralogia 73, 215–30.Google Scholar
Rampone, E., Borghini, G., Romairone, A., Abouchami, W., Class, C. & Goldstein, S. L. 2014. Sm–Nd geochronology of the Erro-Tobbio gabbros (Ligurian Alps, Italy): Insights into the evolution of the Alpine Tethys. Lithos 205, 236–46.Google Scholar
Rampone, E., Hofmann, A. W. & Raczek, I. 2009. Isotopic equilibrium between mantle peridotite and melt: Evidence from the Corsica ophiolite. Earth and Planetary Science Letters 288 (3–4), 601–10.Google Scholar
Ranalli, G. & Murphy, D. C. 1987. Rheological stratification of the lithosphere. Tectonophysics 132 (4), 281–95.Google Scholar
Rebay, G., Riccardi, M. P. & Spalla, M. I. 2015. Fluid rock interactions as recorded by Cl-rich amphiboles from continental and oceanic crust of italian orogenic belts. Periodico di Mineralogia 84 (3B), 751–77.Google Scholar
Rebay, G. & Spalla, M. I. 2001. Emplacement at granulite facies conditions of the Sesia-Lanzo metagabbros: an early record of Permian rifting? Lithos 58 (3–4), 85104.Google Scholar
Reston, T. J. & Morgan, J. P. 2004. Continental geotherm and the evolution of rifted margins. Geology 32 (2), 133–6.Google Scholar
Riccardi, M. P., Tribuzio, R. & Caucia, F. 1994. Amphibole evolution in the metagabbros from East Ligurian ophiolites (Northern Apennines, Italy): Constraints on the ocean-floor metamorphism. Memorie della Società Geologica Italiana 48, 203–8.Google Scholar
Roda, M., Marotta, A. M. & Spalla, M. I. 2010. Numerical simulations of an ocean-continent convergent system: Influence of subduction geometry and mantle wedge hydration on crustal recycling. Geochemistry, Geophysics, Geosystems 11 (5), 121.Google Scholar
Roda, M., Spalla, M. I. & Marotta, A. M. 2012. Integration of natural data within a numerical model of ablative subduction: a possible interpretation for the Alpine dynamics of the Austroalpine crust. Journal of Metamorphic Geology 30 (9), 973996.Google Scholar
Rogers, N., Blake, S. K. B., Widdowson, M., Parkinson, I. & Harris, N. 2008. An Introduction to Our Dynamic Planet. New York: Cambridge University Press, 398 pp.Google Scholar
Rossi, P., Cocherie, A., Lahondère, D. & Fanning, C. M. 2002. La marge Européenne de la Téthys Jurassique en Corse: datation de trondhjémites de Balagne et indices de croûte continentale sous le domaine Balano-Ligure. Comptes Rendus Geoscience 334 (5), 313–22.Google Scholar
Rossi, P., Lahondère, J.-C., Cocherie, A., Caballero, Y. & Féraud, J. 2012. Notice explicative, Carte geol. France (1/50 000), feuille Bastelica (1118). Edition du Bureau de Recherches Géologiques et Minières, Orléans, 134 pp.Google Scholar
Rottura, A., Bargossi, G. M., Caggianelli, A., Del Moro, A., Visonà, D. & Tranne, C. A. 1998. Origin and significance of the Permian high-K calc-alkaline magmatism in the central-eastern Southern Alps, Italy. Lithos 45 (1–4), 329–48.Google Scholar
Rubatto, D., Regis, D., Hermann, J., Boston, K., Engi, M., Beltrando, M. & McAlpine, S. R. B. 2011. Yo-yo subduction recorded by accessory minerals in the Italian Western Alps. Nature Geoscience 4 (5), 338–42.Google Scholar
Rybach, L. 1988. Determination of heat production rate. In Handbook of Terrestrial Heat-Flow Density Determination (eds Haenel, R., Stegena, L. & Rybach, L.), pp. 125142. Solid Earth Sciences Library, Kluwer Academic Publishers.Google Scholar
Sanders, C. A. E., Bertotti, G. & Tommasini, S. 1996. Triassic pegmatites in the Mesoizoic middle crust of the Southern Alps (Italy): fluid inclusions, radiometric dating and tectonic implications. Eclogae Geologicae Helvetiae 89 (1), 505–25.Google Scholar
Sandiford, M. & Powell, R. 1986. Deep crustal metamorphism during crustal extension: modern and ancient examples. Earth and Planetary Science Letters 79, 151–8.Google Scholar
Schmid, S. M., Fügenschuh, B., Kissling, E. & Schuster, R. 2004. Tectonic map and overall architecture of the Alpine orogen. Eclogae Geologicae Helvetiae 97 (1), 93117.Google Scholar
Schuster, R. & Frank, W. 1999. Metamorphic evolution of the Austroalpine units east of the Tauern Window: indications for Jurassic strike slip tectonics. Mitteilungen der Gesellschaft der Geologie und Bergbaustudenten in Österreich 42, 3758.Google Scholar
Schuster, R., Scharbert, S., Abart, R. & Frank, W. 2001. Permo-Triassic extension and related HT/LP metamorphism in the Austroalpine-Southalpine realm. Mitteilungen der Gesellschaft der Geologie und Bergbaustudenten in Österreich 45, 111–41.Google Scholar
Schuster, R. & Stüwe, K. 2008. Permian metamorphic event in the Alps. Geology 36, 603–6.Google Scholar
Seranne, M. 1999. The Gulf of Lion continental margin (NW Mediterranean) revisited by IBS: an overview. In The Mediterranean Basins: Tertiary Extension within the Alpine Orogen (eds. Durand, B., Jolivet, L., Horvath, F. & Seranne, M.), pp. 1536. Geological Society, London, Special Publication no. 156.Google Scholar
Sills, J. D. 1984. Granulite facies metamorphism in the Ivrea Zone, NW Italy. Schweizerische Mineralogische und Petrographische Mitteilungen 64, 169–91.Google Scholar
Smye, A. J. & Stockli, D. F. 2014. Rutile U-Pb age depth profiling: A continuous record of lithospheric thermal evolution. Earth and Planetary Science Letters 408, 171–82.Google Scholar
Spalla, M. I., Lardeaux, J. M., Dal Piaz, G. V. & Gosso, G. 1991. Metamorphisme et tectonique a la marge externe de la zone Sesia-Lanzo (Alpes occidentales). Memorie di Scienze Geologiche 43, 361–9.Google Scholar
Spalla, M. I., Lardeaux, J. M., Dal Piaz, G. V., Gosso, G. & Messiga, B. 1996. Tectonic significance of the Alpine eclogites. Journal of Geodynamics 21 (3), 257–85.Google Scholar
Spalla, M. I., Messiga, B. & Gosso, G. 1995. LT-Alpine overprint on the HT rifting-related metamorphism in the steep belt of the Languard-Campo Nappe. The Cima Rovaia and Scisti del Tonale Units represent two different extents of Alpine re-equilibration. In International Ophiolite Symposium, Pavia, 18–23 September 1995, pp. 148.Google Scholar
Spalla, M. I., Zanoni, D., Gosso, G. & Zucali, M. 2009. Deciphering the geologic memory of a Permian conglomerate of the Southern Alps by pebble P–T estimates. International Journal of Earth Sciences 98 (1), 203–26.Google Scholar
Spalla, M. I., Zanoni, D., Marotta, A. M., Rebay, G., Roda, M., Zucali, M. & Gosso, G. 2014. The transition from Variscan collision to continental break-up in the Alps: insights from the comparison between natural data and numerical model predictions. In The Variscan Orogeny: Extent, Timescale and the Formation of the European Crust (eds Schulmann, K., Martínez Catalán, J. R., Lardeaux, J. M., Janoušek, V. & Oggiano, G.), pp. 363400. Geological Society, London, Special Publication no. 405.Google Scholar
Spear, F. S. & Peacock, S. M. 1989. Metamorphic Pressure‐Temperature‐Time Paths. Washington, DC: American Geophysical Union, 102 pp.Google Scholar
Spicker, G. & Huckenholz, H. G. 1986. Petrography and geochemistry of the Monzoni Intrusion, northern Italy. Fortschritte der Mineralogie, Beiheft 64 (1), 172.Google Scholar
Spiess, R., Cesare, B., Mazzoli, C., Sassi, R. & Sassi, F. P. 2010. The crystalline basement of the Adria microplate in the eastern Alps: a review of the palaeostructural evolution from the Neoproterozoic to the Cenozoic. Rendiconti Lincei Scienze Fisiche Naturali 21, 3150.Google Scholar
Stahle, V., Frenzel, G., Hess, J. C., Saupe, F., Schmidt, S. T. & Schneider, W. 2001. Permian metabasalt and Triassic alkaline dykes in the northern Ivrea zone: clues to the post-Variscan geodynamic evolution of the Southern Alps. Schweizerische Mineralogische Und Petrographische Mitteilungen 81, 121.Google Scholar
Stampfli, G., Mosar, J., Marquer, D., Marchant, R., Baudin, T. & Borel, G. 1998. Subduction and obduction processes in the Swiss Alps. Tectonophysics 296 (1–2), 159204.Google Scholar
Stöckhert, B. & Gerya, T. V. 2005. Pre-collisional high pressure metamorphism and nappe tectonics at active continental margins: a numerical simulation. Terra Nova 17, 102–10.Google Scholar
Stöckhert, B. 1987. Das Uttenheimer Pegmatitfeld (Ostalpines Altkristallin, Suedtirol) Genese und alpine Ueberpraegung. Erlanger Geologische Abhandlungen 114, 83106.Google Scholar
Tenczer, V., Powell, R. & Stuwe, K. 2006. Evolution of H2O content in a polymetamorphic terrane: the Plattengneiss Shear Zone (Koralpe, Austria). Journal of Metamorphic Geology 24 (4), 281–95.Google Scholar
Thompson, A. B. 1981. The pressure–temperature (P,T) plane viewed by geophysicists and petrologists. Terra Cognita 1, 1120.Google Scholar
Thompson, A. B. & England, P. C. 1984. Pressure-temperature-time paths of regional metamorphism II. Their inference and interpretation using mineral assemblages in metamorphic rocks. Journal of Petrology 25 (4), 929–55.Google Scholar
Thöni, M. & Jagoutz, E. 1992. Some new aspects of dating eclogites in orogenic belts: Sm-Nd, Rb-Sr, and Pb-Pb isotopic results from the Austroalpine Saualpe and Koralpe type-locality (Carinthia/Styria, southeastern Austria). Geochimica et Cosmochimica Acta 56 (1), 347–68.Google Scholar
Thöni, M. & Miller, C. 2000. Permo-Triassic pegmatites in the eo-Alpine eclogite-facies Koralpe complex, Austria: age and magma source constraints from mineral chemical, Rb-Sr and Sm-Nd isotope data. Schweizerische Mineralogische und Petrographische Mitteilungen 80 (2), 169–86.Google Scholar
Thöni, M. & Miller, C. 2009. The ‘Permian event’ in the Eastern European Alps: Sm–Nd and P–T data recorded by multi-stage garnet from the Plankogel unit. Chemical Geology 260 (1–2), 2036.Google Scholar
Thöni, M., Miller, C., Blichert-Toft, J., Whitehouse, M. J., Konzett, J. & Zanetti, A. 2008. Timing of high-pressure metamorphism and exhumation of the eclogite type-locality (Kupplerbrunn–Prickler Halt, Saualpe, south-eastern Austria): constraints from correlations of the Sm–Nd, Lu–Hf, U–Pb and Rb–Sr isotopic systems. Journal of Metamorphic Geology 26 (5), 561–81.Google Scholar
Thöni, M., Mottana, A., Delitala, M. C., De Capitani, L. & Liborio, G. 1992. The Val Biandino composite pluton: a Late Hercynian intrusion into the South Alpine metamorphic basement of the Alps (Italy). Neues Jahrbuch fur Mineralogie - Abhandlungen 12, 545–54.Google Scholar
Tribuzio, R., Riccardi, M. P. & Messiga, B. 1997. Amphibolitization of Mg- and Fe-rich gabbroic dykes within mantle-derived serpentinites from Northern Apennine ophiolites: Evidence for high-temperature hydration of the oceanic lithosphere. Ofioliti 22 (1), 7180.Google Scholar
Tribuzio, R., Riccardi, M. P. & Ottolini, L. 1995. Trace element redistribution in high-temperature deformed gabbros from East Ligurian ophiolites (Northern Apennines, Italy): constraints on the origin of syndeformation fluids. Journal of Metamorphic Geology 13 (3), 367–77.Google Scholar
Tribuzio, R., Thirwall, M. F. & Messiga, B. 1999. Petrology, mineral and isotope geochemistry of the Sondalo gabbroic complex (Central Alps, Northern Italy): implications for the origin of post-Variscan magmatism. Contributions to Mineralogy and Petrology 136, 4862.Google Scholar
Tribuzio, R., Thirwall, M. F. & Vannucci, R. 2004. Origin of the gabbro-peridotite association from the Northern Apennine Ophiolites (Italy). Journal of Petrology 45 (6), 1109–24.Google Scholar
Turco, E., Macchiavelli, C., Mazzoli, S., Schettino, A. & Pierantoni, P. P. 2012. Kinematic evolution of Alpine Corsica in the framework of Mediterranean mountain belts. Tectonophysics 579, 193206.Google Scholar
Turner, S. P., George, R. M. M., Evans, P. J., Hawkesworth, C. J. & Zellmer, G. F. 2000. Time-scales of magma formation, ascent and storage beneath subduction-zone volcanoes. Philosophical Transactions of the Royal Society A: Mathematical, Physical and Engineering Sciences 358 (1770), 1443–64.Google Scholar
Van Avendonk, H. J. A., Lavier, L. L., Shillington, D. J. & Manatschal, G. 2009. Extension of continental crust at the margin of the eastern Grand Banks, Newfoundland. Tectonophysics 468 (1–4), 131–48.Google Scholar
Vavra, G., Gebauer, D., Schmid, R. & Compston, W. 1996. Multiple zircon growth and recrystallization during polyphase Late Carboniferous to Triassic metamorphism in granulites of the Ivrea Zone (Southern Alps): an ion microprobe (SHRIMP) study. Contributions to Mineralogy and Petrology 122 (4), 337–58.Google Scholar
Vavra, G., Schmid, R. & Gebauer, D. 1999. Internal morphology, habit and U-Th-Pb microanalysis of amphibolite-to-granulite facies zircons: geochronology of the Ivrea Zone (Southern Alps). Contributions to Mineralogy and Petrology 134 (4), 380404.Google Scholar
Visonà, D. 1995. Polybaric evolution of calc-alkaline magmas: the dioritic belt of the Bressanone-Chiusa igneous complex (NE Italy). Memorie di Scienze Geologiche 47, 111–24.Google Scholar
Visonà, D. 1997. The Predazzo multipulse intrusive body (Western Dolomites, Italy). Field and mineralogical studies. Memorie di Scienze Geologiche 49, 117–25.Google Scholar
Visonà, D., Fioretti, A. M., Poli, M. E., Zanferrari, A. & Fanning, M. 2007. U-Pb SHRIMP zircon dating of andesite from the Dolomite area (NE Italy): geochronological evidence for the early onset of Permian Volcanism in the eastern part of the southern Alps. Swiss Journal of Geosciences 100 (2), 313–24.Google Scholar
von Raumer, J. F., Bussy, F., Schaltegger, U., Schulz, B. & Stampfli, G. M. 2013. Pre-Mesozoic Alpine basements: Their place in the European Paleozoic framework. Geological Society of America Bulletin 125 (1–2), 89108.Google Scholar
Voshage, H., Hunziker, J. C., Hoffmann, A. W. & Zingg, A. 1987. A Nd and Sr isotopic study of Ivrea zone, Southern Alps, N-Italy. Contributions to Mineralogy and Petrology 97, 3142.Google Scholar
Vuichard, J. P. 1987. Conditions P-T du métamorphisme anté-alpin dans la seconde zone diorito-kinzigitique (Zone Sesia-Lanzo, Alpes occidentales). Schweizerische Mineralogische und Petrographische Mitteilungen 67, 257–71.Google Scholar
Wessel, P. & Smith, W. M. F. 2001. New improved version of generic mapping tools released. EOS Transactions of the American Geophysical Union 79, 579.Google Scholar
Whitmarsh, R. B., Manatschal, G. & Minshull, T. A. 2001. Evolution of magma-poor continental margins from rifting to seafloor spreading. Nature 413 (6852), 150–4.Google Scholar
Wogelius, R. A. & Finley, B. G. 1989. Subsolidus emplacement history of the Lanzo massif, northern Italy. Geology 17 (11), 995–9.Google Scholar
Zanoni, D., Spalla, M. I. & Gosso, G. 2010. Vestiges of lost tectonic units in conglomerate pebbles? A test in Permian sequences of the Southalpine Orobic Alps. Geological Magazine 147 (1), 98122.Google Scholar
Zucali, M. 2001. La correlazione nei terreni metamorfici: due esempi dall'Austroalpino occidentale (Zona Sesia-Lanzo) e centrale (Falda Languard-Campo/ Serie del Tonale). Ph.D. thesis, Università degli Studi di Milano. Published thesis.Google Scholar
Figure 0

Figure 1. Tectonic map of the Alps and Apennines (after Marotta & Spalla, 2007; Handy et al.2010) with the locations of: metamorphic rocks and main gabbro bodies of Permian–Triassic age occurring in the pre-Alpine continental crust of the Alps; and mantle rocks and oceanic gabbros from the Alps and Northern Apennines. The codes are defined in Tables 2–5.

Figure 1

Table 1. Material properties used in the 2-D numerical modelling

Figure 2

Figure 2. (a) 2D geometry and numerical set-up of the model. (b) Thermal and rheological profiles at the beginning of the evolution for the hot (red lines) and cold (blue lines) models. The solid lines indicate the effective viscosity profiles, corresponding to the geotherms indicated by the dashed lines.

Figure 3

Figure 3. Successive stages of the tectonic evolution predicted by the hot (panels ai) and cold (panels bi) models at different times after the beginning of the forced extension. Black and red dashed lines correspond to 800 K and 1500 K isotherms, respectively. Ages refer to the time span from the beginning of the simulations.

Figure 4

Figure 4. Thermal and velocity fields of the system at different times after the beginning of the forced extension for the hot (panels ai) and cold (panels bi) models. Time scale in the centre of the figure represents the simulation duration for both models, on which the onset of serpentinization, crustal break-up and partial melting are located by arrows (red for the hot and blue for the cold models). Time intervals between serpentinization and partial melting onset are indicated by dashed lines.

Figure 5

Figure 5. Surface horizontal crustal velocity predicted by the hot (a) and cold (b) models at different times after the beginning of the forced extension.

Figure 6

Figure 6. Maps of the velocity modules throughout the lithosphere at different times after the beginning of the forced extension for the hot (panels ai) and cold (panels bi) models.

Figure 7

Table 2. Permian–Triassic metamorphic rocks from the continental crust of the Alps. Labels listed in the code column are reported in Figures 1 and 9–11 to indicate rock positions and duration of thermal fitting

Figure 8

Figure 7. Timing of the metamorphic and magmatic events from Variscan to Jurassic time along the Alpine belt. Radiometric ages are plotted including analytical uncertainty intervals (dashed horizontal bars). Variscan evolutions are synthesized by data in Table 2. Permian–Triassic metamorphic and igneous data (continental gabbros) are synthesized by data in Table 3 and represented with their error margin. Ophiolites and mantle age data are listed in Tables 4 and 5. The grey stripe indicates the age of the earliest radiolarian cherts related to Alpine ophiolites (Cordey & Bailly, 2007).

Figure 9

Table 3. Permian–Triassic gabbros emplaced in the pre-Alpine continental crust of the Alps and Apennines. Labels listed in the code column are reported in Figures 1 and 9–11 to indicate rock positions and duration of thermal fitting

Figure 10

Table 4. Subcontinental and oceanic mantle peridotites from the Alps and Apennines. Labels listed in the code column are reported in Figures 1 and 9–11 to indicate rock positions and duration of thermal fitting

Figure 11

Table 5. Ophiolitic gabbros from the Alps, Apennines and Corsica. The lack of P–T estimates for Corsica gabbros inhibits their comparison with the model predictions. Labels listed in the code column are reported in Figures 1 and 9–11 to indicate rock positions and duration of thermal fitting

Figure 12

Figure 8. Time references from the hot and cold models with respect to the age of the oldest gabbros (160 Ma, 170 Ma and 185 Ma), based on the literature for the Northern Apennines and Western Alps (Tribuzio, Riccardi & Ottolini, 1995; Li et al.2013; Rampone et al.2014). Green circles identify the times when the mantle partial melting occurs in the hot model (36.4 Ma) and in the cold model (22.4 Ma).

Figure 13

Figure 9. Duration of the agreement between the predictions and the natural data as well as number of fitting markers (colours) for the three different ages of the oldest gabbros (160 Ma, 170 Ma and 185 Ma, dashed black lines) in terms of lithological affinity and coincident PT values compared to the radiometric (black thick segments) and geologic (grey thick segments) ages of the natural data. In the text, we refer to the result as a ‘complete fit’ when the model predictions and the natural data of the lithological affinity, PT values and ages agree. The fit is considered ‘partial’ if age agreement is lacking. Panel (a) refers to the hot model, whereas panel (b) refers to the cold model. The light grey area represents the duration of the numerical simulation. The codes are defined in Tables 2–5.

Figure 14

Figure 10. Spatial distribution of the markers that guarantee a complete fit of the data for the hot model at different stages of the evolution. The dashed black lines indicate the 800 K and 1500 K isotherms. The colours identify the markers showing a complete fit with the natural data as specified in the legend; tr indicates the time span from the beginning of the simulation and ta indicates the absolute age constrained choosing the oldest gabbro age at 170 Ma (see discussion in the text). The codes are defined in Tables 2–5.

Figure 15

Figure 11. Spatial distribution of the markers that guarantee a complete fit of the data for the cold model at different stages of the evolution. The dashed black lines indicate the 800 K and 1500 K isotherms. The colours identify the markers showing a complete fit with the natural data specified in the legend; tr indicates the time span from the beginning of the simulation and ta indicates the absolute age constrained choosing the oldest gabbro age at 170 Ma (see discussion in the text). The codes are defined in Tables 2–5.