The Betic Ophiolites and the Mesozoic Evolution of the Western Tethys

The Betic Ophiolites consist of numerous tectonic slices, metric to kilometric in size, of eclogitized mafic and ultramafic rocks associated to oceanic metasediments, deriving from the Betic oceanic domain. The outcrop of these ophiolites is aligned along 250 km in the Mulhacén Complex of the Nevado-Filábride Domain, located at the center-eastern zone of the Betic Cordillera (SE Spain). According to petrological/geochemical inferences and SHRIMP (Sensitive High Resolution Ion Micro-Probe) dating of igneous zircons, the Betic oceanic lithosphere originated along an ultra-slow mid-ocean ridge, after rifting, thinning and breakup of the preexisting continental crust. The Betic oceanic sector, located at the westernmost end of the Tethys Ocean, developed from the Lower to Middle Jurassic (185–170 Ma), just at the beginning of the Pangaea break-up between the Iberia-European and the Africa-Adrian plates. Subsequently, the oceanic spreading migrated northeastward to form the Ligurian and Alpine Tethys oceans, from 165 to 140 Ma. Breakup and oceanization isolated continental remnants, known as the Mesomediterranean Terrane, which were deformed and affected by the Upper Cretaceous-Paleocene Eo-Alpine high-pressure metamorphic event, due to the intra-oceanic subduction of the Jurassic oceanic lithosphere and the related continental margins. This process was followed by the partial exhumation of the subducted oceanic rocks onto their continental margins, forming the Betic and Alpine Ophiolites. Subsequently, along the Upper Oligocene and Miocene, the deformed and metamorphosed Mesomediterranean Terrane was dismembered into different continental blocks collectively known as AlKaPeCa microplate (Alboran, Kabylian, Peloritan and Calabrian). In particular, the Alboran block was displaced toward the SW to occupy its current setting between the Iberian and African plates, due to the Neogene opening of the Algero-Provençal Basin. During this translation, the different domains of the Alboran microplate, forming the Internal Zones of the Betic and Rifean Cordilleras, collided with the External Zones representing the Iberian and African margins and, together with them, underwent the later alpine deformation and metamorphism, characterized by local differences of P-T (Pressure-Temperature) conditions. These Neogene metamorphic processes, known as Meso-Alpine and Neo-Alpine events, developed in the Nevado-Filábride Domain under Ab-Ep amphibolite and greenschists facies conditions, respectively, causing retrogradation and intensive deformation of the Eo-Alpine eclogites. Geosciences 2017, 7, 31; doi:10.3390/geosciences7020031 www.mdpi.com/journal/geosciences Geosciences 2017, 7, 31 2 of 29


Introduction
The ultramafic and mafic rocks included in the Mulhacén Complex of the Betic Cordillera (SE Spain; Figure 1) have been interpreted as an eclogitized ophiolite association, i.e., metamorphosed remnants of a Mesozoic oceanic lithosphere since the 1970s (Puga [1], Puga and Díaz de Federico [2]).This hypothesis relating the Betic Ophiolite Association (hereafter named as BOA following Puga [3]) to a Mesozoic oceanic basin was criticized at the beginning, due to the paucity of outcrops conforming to the three layers model defined for Ophiolites in the Penrose Conference (1972), i.e., sections including peridotites/serpentinites, gabbros/basalts and deep sea sediments, such cherts, hypothetically corresponding to the stratigraphy of the oceanic lithosphere.In spite of the criticism, the Betic Ophiolites were identified as Mid-Oceanic-Ridge MOR-type by Puga [1,3], Bodinier et al. [4], Morten et al. [5], Puga et al. [6,7], who investigating BOA basic and ultramafic lithologies emphasized geochemical analogies with the rocks formed along the Atlantic Ridge and those included within the Jurassic Alpine-Apenninic Ophiolites.Moreover, later palaeogeographic, petrologic and geochemical studies [8][9][10][11] corroborated the ophiolitic origin of the BOA rocks.Nevertheless, Gómez-Pugnaire et al. [12] highlighted apparent petrological incongruence, emphasizing that the geochemical composition of Betic basalts/gabbros do not conform to Normal Mid Oceanic Ridge Basalt (N-MORB) features, i.e., tholeiitic character and extreme depletion in incompatible trace elements, and/or that the Betic ultramafic rocks tend to be less residual (i.e., lherzolite, scarcely depleted by partial melting) with respect to the abyssal peridotites having prevalent harzburgite composition.
Therefore, the aim of this contribution is to summarize the current understanding of BOA sequences, emphasizing analogies and differences with the adjoining rock associations and synthetizing possible geodynamic relationships.In this view, the BOA represents a key jigsaw of the Mesozoic Tethys puzzle and the definition of its paleogeographic and geodynamic evolution contributes to the understanding of the current geological configuration of the Central-Western Mediterranean area.

Geographical and Geological Setting of the Betic Ophiolites
The Betic Ophiolites (BOA), crop out discontinuously along some 250 km in the central and eastern sectors of the Betic Cordillera (SE Spain) (Figure 1).These ophiolites are represented by tectonic slabs ranging in size from metres to kilometres and made up of metamorphosed ultramafic and mafic lithologies sometimes associated to sedimentary rocks.Figure 1 shows the location of the most significant betic ophiolitic outcrops from which most of the petrological, geochemical and geochronological data mentioned in this study have been obtained.These ophiolites are exclusively located in the Mulhacén Complex (MC), in the central and eastern area of the Betic Cordillera, and are never found in the underlying Veleta Complex (VC) [13,16,23,37,38].This suggests that these two complexes pertained to different palaeogeographic domains before and during the Jurassic-Cretaceous period, when the ophiolites originated from exhumation of an oceanic floor [8,10,15,[17][18][19]39,40].
At present the rocks of the VC crop out forming several tectonic windows below the (MC), and both metamorphic complexes, forming the Nevado-Filábride Domain, are the deepest in the current tectonic pile forming the Betic Cordilleras (Figure 1).
Geosciences 2017, 7, 31 3 of 29 [8,10,15,[17][18][19]39,40].At present the rocks of the VC crop out forming several tectonic windows below the (MC), and both metamorphic complexes, forming the Nevado-Filábride Domain, are the deepest in the current tectonic pile forming the Betic Cordilleras (Figure 1).The BOA outcrops, deriving from the betic oceanic floor [8,13,19], allow the reconstruction of the pristine oceanic lithostratigraphy represented in Figure 2, which has been tectonically superimposed on the crustal rocks of its adjoining western continental margin, currently forming the Caldera crustal unit of the MC, and were overthrusted by the rocks deriving from its eastern continental margin, now forming the Sabinas crustal unit (Figure 2).Each of these crustal units consists of rocks deriving from a Paleozoic basement with abundant granitic and rhyolitic orthogneisses, covered by Triassic sediments [13,16,37,38,41,42] and, together with the Jurassic-Cretaceous Ophiolite Unit, constitute the MC. Figure 2 also shows, in its legend, the main lithologic differences existing between the Caldera and the Sabinas units, as well as the more notable differences characterizing the currently underlying VC by respect to the overlying MC, and their different tectonic units [13,23].
The reconstruction of the primitive relationship of the different rock types forming the Ophiolite Unit on the oceanic floor is mainly based on the present mutual location among the rocks forming the different outcrops of the metamorphosed and dismembered ophiolitic relics, and also The BOA outcrops, deriving from the betic oceanic floor [8,13,19], allow the reconstruction of the pristine oceanic lithostratigraphy represented in Figure 2, which has been tectonically superimposed on the crustal rocks of its adjoining western continental margin, currently forming the Caldera crustal unit of the MC, and were overthrusted by the rocks deriving from its eastern continental margin, now forming the Sabinas crustal unit (Figure 2).Each of these crustal units consists of rocks deriving from a Paleozoic basement with abundant granitic and rhyolitic orthogneisses, covered by Triassic sediments [13,16,37,38,41,42] and, together with the Jurassic-Cretaceous Ophiolite Unit, constitute the MC. Figure 2 also shows, in its legend, the main lithologic differences existing between the Caldera and the Sabinas units, as well as the more notable differences characterizing the currently underlying VC by respect to the overlying MC, and their different tectonic units [13,23].
by the absolute ages of many dated rocks presented in Table 1.The main rock types forming the Ophiolite Unit in the lithostratigraphic reconstruction shown in Figure 2 are, from bottom to top: Serpentinites, intruded by rodingitized dolerite dykes, eclogitized and amphibolitized gabbros, dolerite dykes and basalts.These rock types may be found covered by a metasedimentary sequence in which the intercalation of amphibolitized dolerite sills is common, probably partly deriving from thin basaltic levels.

Petrology and Metamorphic Evolution
The rocks included in the BOA lithostratigraphy are reported in Figure 2 as located in the original oceanic lithosphere, before being affected by Alpine tectonic activity sensu lato, and the more representative lithologies are illustrated in Figure 3a-c.These figures show the rocks on the The reconstruction of the primitive relationship of the different rock types forming the Ophiolite Unit on the oceanic floor is mainly based on the present mutual location among the rocks forming the different outcrops of the metamorphosed and dismembered ophiolitic relics, and also by the absolute ages of many dated rocks presented in Table 1.The main rock types forming the Ophiolite Unit in the lithostratigraphic reconstruction shown in Figure 2 are, from bottom to top: Serpentinites, intruded by rodingitized dolerite dykes, eclogitized and amphibolitized gabbros, dolerite dykes and basalts.These rock types may be found covered by a metasedimentary sequence in which the intercalation of amphibolitized dolerite sills is common, probably partly deriving from thin basaltic levels.

Petrology and Metamorphic Evolution
The rocks included in the BOA lithostratigraphy are reported in Figure 2 as located in the original oceanic lithosphere, before being affected by Alpine tectonic activity sensu lato, and the more representative lithologies are illustrated in Figure 3a-c.These figures show the rocks on the outcrops, as well as their microscopic views, also specifying the outcrops of provenance of the different photos.Figure 3a represents the main lithologies of the ultramafic sequence, which are serpentinized spinel lherzolites, containing antigorite (Atg), diopside relics (Di) and neoformed tremolite (Tr) (Figure 3a-2).This rock-type presents gradual transition to metamorphic rocks consisting of acicular aggregates of olivine (Ol) and enstatite (En), with pseudo-spinifex texture, known as secondary harzburgites (Figure 3a-3).Both types of ultramafites contain numerous rodingitized and/or eclogitized boudinaged dykes (Figure 3a-4).Thin section microscopic view of a dolerite dyke transformed into eclogite, which is formed by almandine (Alm) omphacite (Omp) and rutile (Rt) paragenesis, partially amphibolitized with development of albite (Ab), epidote (Ep) and kataphoritic amphibole (Ktp).These rocks, of the mantle sequence, can be found directly overlayered by rocks of the sedimentary sequence (Figure 2), formed by near pure quartzites, probably deriving from radiolarites (Figure 3a-1), and/or tourmaline-rich micaschists and calc-schists.The reported photos correspond to rocks from the Almirez outcrop, but similar litohologies are found in the Santillana and Montenegro outcrops [10] (Figure 1).
Figure 3b corresponds to cumulitic troctolites and gabbros, intruded by dolerite dykes (Figure 3b-1), mainly preserving the plutonic structures and locally also the igneous paragenesis (Figure 3b-2,b-3), but more commonly transformed into coronitic eclogites, partly retrograded to amphibolites (Figure 3b-4).Figure 3b-2,b-3 show the local preservation at thin-section scale of igneous olivine (Ol), augitic clinopyroxene (Cpx) and calcic-plagioclase (Pl).Figure 3b-4 shows a coronitic eclogite originated from a cumulitic gabbro.In this metamorphic rock the igneous olivines have been pseudomorphosed by a garnet corona (Gr) surrounding an aggregate of omphacite (Omp) and amphibole (Amp), whereas the igneous plagioclases have been pseudomorphosed by an aggregate of albite (Ab) and clinozoisite (Czo).The preservation of some minerals of the igneous paragenesis, despite the later orogenic metamorphism is recognized only at thin section scale, and mainly in some outcrops like that of Cóbdar where eclogitization was less pervasive than in other sites, like that of Algarrobo where coronitic eclogites are very common (Figure 3b-4).
Figure 3c represents well preserved pillow structures in basaltic lithologies characterized by different degrees of alpine deformation (Figures 2 and 3c-1).The igneous parageneses are locally replaced by brown amphibole formed during oceanic-floor metamorphic processes (Figure 3c-3) or, more commonly, by the Eo-Alpine eclogitic parageneses partly retrograded to Meso-Alpine amphibolites (Figure 3c-4).In some basaltic levels the effects of the ocean-floor metamorphism has been locally well recorded by high temperature Ti-pargasite (Ti-Prg) brown amphibole (Figure 3c-3) in the groundmass and in millimetric veins, as observed in metabasalts from Cóbdar.Noteworthy, brown amphibole from veins of the Cóbdar metabasalts were dated with Ar-Ar laser probe, indicating a Jurassic age for the oceanic-floor metamorphic process (Table 1).
The Lugros outcrop (Figure 1) shows the best examples of preservation of igneous structures, such as pillows and lava flow-structures despite the complete transformation of the igneous paragenesis into eclogite facies paragenesis.In particular, Figure 3c-1 represents a good example of a minipillows aggregate in which the volcanic structure has been exceptionally well preserved, despite the complete eclogitization of the igneous paragenesis, as shown in Figure 3c-4.This microscopic photo shows a well preserved igneous variolitic texture, formed by a fibrous aggregate of volcanic plagioclase microcrystals later included in an almandinic garnet poeciloblast (Alm), which together omphacite (Omp) and barroisite (Bar) forms the eclogite paragenesis and its later local amphibolitization.
Both intrusive and volcanic basic rocks are locally covered by an oceanic meta-sedimentary sequence, similar to that previously described onto the ultramafic rocks, as shown in Figure 2. In this metasedimentary sequence calc-schists with centimetric ankerite nodules and some decimetric alternances of marbres and amphibolite sills are very common.Some relics of foraminifer probably cretaceous in age were found by Tendero et al. [39] into the calc-schists of the Cóbdar outcrop.Worth of note, in the Penninic ophiolite units of the Alps, ankeritic nodules similar to those observed in the BOA sedimentary sequence have also been interpreted as Cretaceous planktonic foraminifera, and have been used to establish the chronostratigraphy of metamorphic successions such as the "Schistes lustres" by Lemoine et al. [43] and Lagabrielle et al. [44].
On the whole, the metamorphic evolution inferred from the study of the igneous and sedimentary lithotypes constituting the Betic Ophiolites, together with the radiometric dating of some of their minerals, allow us to deduce their complex pressure-temperature-time evolution as schematized in Figure 4.
The range of ages corresponding to the different metamorphic events shown in Figure 4 can be seen in [13] and in Table 1 of this paper, and summarized as follows: (a) Ocean floor metamorphism and metasomatism induced serpentinization stage of the ultramafic rocks and rodingitization of the intruded doleritic dykes (Puga et al. [10,17,19]; Alt et al. [22]).
It also produced paragenesis of very high gradient in amphibolite facies, characteristic of ocean floor metamorphism, preserved in some gabbros and basalts [11,14,21,40].The P-T conditions of this oceanic metamorphism and metasomatism, represented in Figure 4, mainly correspond to the range of P-T calculated for the brown amphiboles of the Cóbdar metabasites, shown in Photo 3c-3, by Puga et al. [14].(b) Orogenic metamorphism, begining with subduction of the ocean floor during an initial metamorphic event known as Eo-Alpine, which attained eclogite facies conditions.This was followed by exhumation of part of this subducted ocean floor onto the continental margin, during and after which the Meso-Alpine and Neo-Alpine metamorphic events of Oligocene and Miocene ages took place [10,11,13,20].The successive parageneses of eclogite, Ab-Ep amphibolite and greenschist facies developed during these three metamorphic events, variously overprinting the pristine igneous parageneses and textures [11,[13][14][15]19,40].
The P-T trajectories shown in the upper inset of Figure 4 represent the temporal evolution of the orogenic metamorphism in the BOA, and other MC rocks, from the eclogite facies, with the peak P-T conditions approximately represented by M1, which were followed by Ab-Ep amphibolite and greenschist facies conditions, with a peak P-T conditions represented by M2.The more precise P-T conditions for the BOA rocks in these successive stages of alpine metamorphism are represented in Figure 4, into the eclogitic fields M1Ob and M1Ou, for basic and ultramafic rocks, respectively, and into the M2 field showing the conditions for the later partial transformation of eclogites into Ab-Ep amphibolites.A more precise explanation of the reactions selected to delimitate the successive metamorphic field conditions attained in the BOA rocks and in the other units forming part of the MC and VC can be seen in Puga et al. [10,11,13].
The P-T trajectories shown in the upper inset of Figure 4 represent the temporal evolution of the orogenic metamorphism in the BOA, and other MC rocks, from the eclogite facies, with the peak P-T conditions approximately represented by M1, which were followed by Ab-Ep amphibolite and greenschist facies conditions, with a peak P-T conditions represented by M2.The more precise P-T conditions for the BOA rocks in these successive stages of alpine metamorphism are represented in Figure 4, into the eclogitic fields M1Ob and M1Ou, for basic and ultramafic rocks, respectively, and into the M2 field showing the conditions for the later partial transformation of eclogites into Ab-Ep amphibolites.A more precise explanation of the reactions selected to delimitate the successive metamorphic field conditions attained in the BOA rocks and in the other units forming part of the MC and VC can be seen in Puga et al. [10,11,13].

Geochemical Characteristics and Tectonic Setting of the BOA Magmas
The chemical composition of representative rocks from selected outcrops of the Betic Ophiolites are shown in tectono-magmatic diagrams that discriminate the geodynamic context of distinct basaltic magmas, together with the mean compositions of these current magma types from different geodynamic settings for comparison (Figure 5).The symbols for the different types of ophiolitic rocks are explained in the legend of the figure, together with indication of the outcrops location.

Geochemical Characteristics and Tectonic Setting of the BOA Magmas
The chemical composition of representative rocks from selected outcrops of the Betic Ophiolites are shown in tectono-magmatic diagrams that discriminate the geodynamic context of distinct basaltic magmas, together with the mean compositions of these current magma types from different geodynamic settings for comparison (Figure 5).The symbols for the different types of ophiolitic rocks are explained in the legend of the figure, together with indication of the outcrops location.In spite of metamorphism, geochemical tracers such as rare earth and high field strength elements (REE, HFSE) still provide useful information on the original tectono-magmatic setting.BOA rocks show geochemical analogies with ophiolites from the Alps and the Apennines, whereas they differ from ophiolites widespread in the eastern Mediterranean, where supra-subduction affinities have been recorded [25,51,52].
Comparison in Figure 5A-D of the BOA magma compositions with those corresponding to the current magmas from different tectonic setting, indicate that the geochemical affinity of the different types of basic rocks in the Betic Ophiolites is mainly E-MORB or T-MORB, and locally N-MORB only for the case of magmas from Algarrobo outcrop (Figure 5C,D).
The basic rocks forming all the studied Betic Ophiolites outcrops, represented in Figure 1, are derived from Ti-rich tholeiitic magmas (Figure 5A).This magmatic affinity is characteristic of Ti-rich ophiolites, derived from oceanic floor created on oceanic ridges far from subduction zones (Beccaluva et al. [25]).
Most of the basic rocks from the Betic Ophiolites plotted in Figure 5C,D are consistent with the E-MORB average value of Sun and McDonough [53], and with the present-day volcanism at 63 • N, and less commonly at 45 • N, on the Mid-Atlantic Ridge.The enrichment in Th/Yb presented by some samples in Figure 5D, especially in rocks from Cóbdar and Algarrobo outcrops, can be related to assimilation of the thinned and rifted crustal rocks preexisting to the opening of the Betic Ocean, or to interaction with crustal rocks located at the continental margin.Tholeiitic magmas similar to those that formed the Betic Ophiolites are typical of slowly expanding oceanic ridges, especially in the initial stages of oceanization, when embryonic ridges evolve after continental breakup, i.e., during a transition between rifting and drifting before the attainment of a steady-state spreading regime, and therefore may be classified as "Continental margin ophiolites" according to the classification of [30,33] and Pearce [54] classification.
The geochemical and petrological characteristics of the Betic Ophiolites suggest an ultraslow-spreading ridge environment, similar to that hypothesized by Piccardo [55] for the Jurassic Ligurian Tethys.In fact, the current ultraslow ridges (such as the Gakkel Ridge in the Arctic region and the Southwest Indian Ridge) are typically characterized by relatively scarce and episodic magmatic activity, mainly originating in T-to E-MORB melts forming the "magmatic segments" distinguished by the abundance of basic rocks, linked by "amagmatic segments" in which serpentinized peridotite is the predominant lithology (Dick et al. [32]; Michael et al. [56]).This alternance of predominant basic and ultramafic rocks along the current ultra-slow ridges is similar to that characterizing the lithologies found in the Betic ophiolitic outcrops, along their extension of about 250 km in the Mulhacén Complex (Figure 1), despite the fact that the Betic oceanic-floor was subducted and subsequently exhumed onto the continental margin and, consequently, tectonically dislacerated during orogenic processes.
The enrichment of the BOA basic rocks in incompatible trace elements generating a predominant E-MORB and T-MORB character, by respect to N-MORB, could in theory be related to mantle sources activated by upraising mantle plumes comparatively enriched in trace elements (Gale et al. [34]; Herbrich et al. [35]), but in our view, for the BOA study case, is more likely attributed to magma genesis by (1) relatively low melting degree; and (2) associated with metasomatized lithospheric mantle sources.These inferences are perfectly compatible with both an OCT setting and an ultraslow-spreading ridge environment (see discussion below), where mantle sources were not depleted by prolonged magmatism.
The broad isotopic variation of the 87 Sr/ 86 Sr ratio in the Betic Ophiolites (Figure 5B) can be explained by the different degrees of metasomatism undergone by the basic and ultramafic rocks at the ocean floor.This process affected more intensively the basalts than the gabbros of each one of the studied outcrops (Lugros, Cóbdar and Algarrobo) and, in higher degree, the ultramafic rocks and their hosted dolerite dykes, which attained the highest Sr isotopic ratios, typical of Jurassic water (0.7072-0.7075, Jones et al. [57]), during their oceanic metasomatism, originating the rodingitisation and the coeval serpentinisation process [10].Similarly, abyssal peridotites from current oceanic fracture zones (plotted in Figure 5B) also show an increase in 87 Sr/ 86 Sr isotope values due to contamination by present-day seawater ( 87 Sr/ 86 Sr = 0.7092; Snow et al. [48]), which also affected this isotopic ratio for some of the Almirez serpentinites plotted in the same diagram.
Concluding, in our view the geochemical characteristics of the Betic Ophiolites are similar to those forming some of the outcrops of the Alpine-Apennine ophiolitic chain, which also derive from the Jurassic oceanic floor of the Western Alpine Tethys [19,25,58,59].

Radiometric Dating of the BOA
A relevant contribution of our researches has been the dating of the BOA rocks through the U-Pb isotopic analyses of zircons.These crystals although rare in the oceanic-floor basic rocks, have been separated from selected samples by processing about 20-30 kg of material from each sample.Note that in this contribution we preliminarily present some new igneous dating on zircons from rocks collected in the localities of Cóbdar and Algarrobo.These radiometric data implement those presented by Puga et al. [15,19] on rocks from the localities of Lugros and Almirez (Figure 1).
The cathodoluminescence images (CL) of four selected zircons from Lugros, Cóbdar, Algarrobo and Almirez outcops, are shown in Figure 6.The CL photos show three main microstructures in the studied zircons: (i) some crystals with idiomorphic rims and oscillatory zoning, which indicate the igneous origin of the zircon domains, together with (ii) locally recrystallized domains, which do not preserve the oscillatory zoning, probably due to oceanic-floor recrystallization processes, and (iii) crystal sectors presenting whitish color and clouded irregular forms that surround small rounded empty spaces (black in the CL photos).The latter bright and irregular crystal domains have been interpreted (by Puga et al. [15,19]) as developed during a fluid-assisted recrystallisation process of the zircon areas originally surrounding micropores that were previously filled by fluid inclusions.In the zircon crystal from Lugros (Figure 6A) these recrystallized whitish areas yield an age of Upper Cretaceous (83.7 ± 1 Ma), while the oscillatory zoned rim of the same crystal yield an igneous age of Lower Jurassic (189.4 ± 2.2 Ma); in the Almirez zircon crystal (Figure 6D), with igneous age of 186.7 ± 2 Ma, the recrystallysed area, surrounding a micropore, yield a Paleocene age (59.2 ± 0.8 Ma).This range of recrystallization ages from Upper Cretaceous to Paleocene would correspond to the Eo-Alpine metamorphic stage, which affected these rocks in eclogite facies conditions, during the subduction of the oceanic floor from which the Betic Ophiolites derive (Table 1).
contamination by present-day seawater ( 87 Sr/ 86 Sr = 0.7092; Snow et al. [48]), which also affected this isotopic ratio for some of the Almirez serpentinites plotted in the same diagram.
Concluding, in our view the geochemical characteristics of the Betic Ophiolites are similar to those forming some of the outcrops of the Alpine-Apennine ophiolitic chain, which also derive from the Jurassic oceanic floor of the Western Alpine Tethys [19,25,58,59].

Radiometric Dating of the BOA
A relevant contribution of our researches has been the dating of the BOA rocks through the U-Pb isotopic analyses of zircons.These crystals although rare in the oceanic-floor basic rocks, have been separated from selected samples by processing about 20-30 kg of material from each sample.Note that in this contribution we preliminarily present some new igneous dating on zircons from rocks collected in the localities of Cóbdar and Algarrobo.These radiometric data implement those presented by Puga et al. [15,19] on rocks from the localities of Lugros and Almirez (Figure 1).
The cathodoluminescence images (CL) of four selected zircons from Lugros, Cóbdar, Algarrobo and Almirez outcops, are shown in Figure 6.The CL photos show three main microstructures in the studied zircons: (i) some crystals with idiomorphic rims and oscillatory zoning, which indicate the igneous origin of the zircon domains, together with (ii) locally recrystallized domains, which do not preserve the oscillatory zoning, probably due to oceanic-floor recrystallization processes, and (iii) crystal sectors presenting whitish color and clouded irregular forms that surround small rounded empty spaces (black in the CL photos).The latter bright and irregular crystal domains have been interpreted (by Puga et al. [15,19]) as developed during a fluid-assisted recrystallisation process of the zircon areas originally surrounding micropores that were previously filled by fluid inclusions.In the zircon crystal from Lugros (Figure 6A) these recrystallized whitish areas yield an age of Upper Cretaceous (83.7 ± 1 Ma), while the oscillatory zoned rim of the same crystal yield an igneous age of Lower Jurassic (189.4 ± 2.2 Ma); in the Almirez zircon crystal (Figure 6D), with igneous age of 186.7 ± 2 Ma, the recrystallysed area, surrounding a micropore, yield a Paleocene age (59.2 ± 0.8 Ma).This range of recrystallization ages from Upper Cretaceous to Paleocene would correspond to the Eo-Alpine metamorphic stage, which affected these rocks in eclogite facies conditions, during the subduction of the oceanic floor from which the Betic Ophiolites derive (Table 1).The Lugros eclogite histogram, Figure 7A, shows the average age corresponding to the Pliensbachian (187.4 ± 2.5 Ma), as well as a series of ages interpreted as corresponding to oceanic-floor recrystallisation between about 165 and 130 Ma, and Eo-Alpine metamorphic ages, comprised between 90 and 55 Ma.These latter ages are better represented in Lugros than in the partly eclogitized gabbros from Códar (Figure 7B) and Algarrobo (Figure 7C) outcrops, that were less pervasively recrystallized after their igneous crystallization (Figure 3b-2,b-3).The histogram of the dolerite dyke from Almirez outcrop (Figure 7D) shows, together with the Pliensbachian crystallisation age (with average value of 183.9 ± 1.6 Ma), another subordinate peak (with average age of 177.5 ± 2.0 Ma), corresponding to the more luminescent discontinuous rim of their zircons (Figure 6D), which is interpreted as the overgrowth of the igneous crystals during the oceanic rodingitisation process, as well as other recrystallization ages, about 170-150 Ma, which must correspond to the Upper-Middle Jurassic oceanic-floor metasomatism which strongly affected the ultramafic rocks and their hosted dolerite dykes.In this rock the Paleocene Eo-Alpine metamorphism, which is suggested by Figure 6D, is not registered in the histogram of Figure 7D.However, the complex evolution history, recorded by zircons of dolerite dykes intruded in the Almirez ultramafic rocks, which were rodingitized and later locally eclogitized, has been also delineated by the stable isotopes investigation of the hosting serpentinites (Alt et al. [22]).
It is important to note that all the dated rock types yielded, for the oscillatory zoned areas of their zircons, a common average magmatic age around 185 Ma ( [15,19]; Figure 7 and Table 1).This age is interpreted as the beginning of the oceanic crust development for the gabbros forming the basal layer of Lugros, Cóbdar and Algarrobo outcrops, and as the age of intrusion of the dolerite dykes into the ultramafic rocks forming the Almirez outcrop.
The Lugros eclogite histogram, Figure 7A, shows the average age corresponding to the Pliensbachian (187.4 ± 2.5 Ma), as well as a series of ages interpreted as corresponding to oceanic-floor recrystallisation between about 165 and 130 Ma, and Eo-Alpine metamorphic ages, comprised between 90 and 55 Ma.These latter ages are better represented in Lugros than in the partly eclogitized gabbros from Códar (Figure 7B) and Algarrobo (Figure 7C) outcrops, that were less pervasively recrystallized after their igneous crystallization (Figure 3b-2,b-3).The histogram of the dolerite dyke from Almirez outcrop (Figure 7D) shows, together with the Pliensbachian crystallisation age (with average value of 183.9 ± 1.6 Ma), another subordinate peak (with average age of 177.5 ± 2.0 Ma), corresponding to the more luminescent discontinuous rim of their zircons (Figure 6D), which is interpreted as the overgrowth of the igneous crystals during the oceanic rodingitisation process, as well as other recrystallization ages, about 170-150 Ma, which must correspond to the Upper-Middle Jurassic oceanic-floor metasomatism which strongly affected the ultramafic rocks and their hosted dolerite dykes.In this rock the Paleocene Eo-Alpine metamorphism, which is suggested by Figure 6D, is not registered in the histogram of Figure 7D.However, the complex evolution history, recorded by zircons of dolerite dykes intruded in the Almirez ultramafic rocks, which were rodingitized and later locally eclogitized, has been also delineated by the stable isotopes investigation of the hosting serpentinites (Alt et al. [22]).
It is important to note that all the dated rock types yielded, for the oscillatory zoned areas of their zircons, a common average magmatic age around 185 Ma ( [15,19]; Figure 7 and Table 1).This age is interpreted as the beginning of the oceanic crust development for the gabbros forming the basal layer of Lugros, Cóbdar and Algarrobo outcrops, and as the age of intrusion of the dolerite dykes into the ultramafic rocks forming the Almirez outcrop.The U/Pb radiometric dating by SHRIMP of distinct zircon domains distinguished in the studied zircons allowed us to determine with precision an igneous Lower Jurassic age, in the range of 190 to 180 Ma, with average age of 185 Ma, likely representing the beginning of the oceanic magmatism (Figures 6 and 7).Other groups of ages, although not as well defined, correspond to Lower to Middle Jurassic age, probably corresponding to ocean-floor metamorphism and metasomatism, and between the Upper Cretaceous to Paleocene, corresponding to the first orogenic Eo-Alpine recrystallization process that took place during the subduction and eclogitization of the oceanic floor (see Table 1).These two stages of recrystallization of the igneous zircons have not been ubiquitously developed, and/or preserved, in all the samples after to have undergone the Meso-Alpine and Neo-Alpine metamorphic events.They have been better preserved in the zircons from Lugros and Almirez outcrops as is shown in Figures 6 and 7 and mainly in Puga et al. [15,19].
The oceanic-floor metamorphism has been also well preserved in the brown amphibole (Ti-pargasite and Kaersutite) filling little veins in ophiolitic metabasites, which has been dated by 40 Ar/ 39 Ar laser probe as Middle Jurassic (Puga et al. [60]).Moreover, the glaucophane crystals formation, present in the last stage of some eclogites development [11], could correspond to the exhumation process of some oceanic floor slabs onto the continental margin.The crystallization of this Na-amphibole has been dated as near to 70 Ma by the K/Ar method in a Cóbdar amphibolite (Portugal Ferreira et al., [61]).The following metamorphic events, corresponding to the Meso-Alpine and Neo-Alpine stages, identified in BOA and in other Nevado-Filábride units, have produced scarce recrystallization in the zircons of the ophiolitic metabasites, although they have been radiometrically dated, mainly by Ar/Ar and K/Ar methods, corresponding to Eocene-Oligocene and Miocene ages, respectively (Table 1).
The most characteristic minerals developed during the Meso-Alpine metamorphism, in Ab-Ep amphibolite facies conditions, are the Na-Ca amphiboles, which replace the omphacite, garnet and glaucophane of the eclogitic paragenesis and, more locally, the high gradient brown amphibole relics formed during the Middle-Upper Jurassic ocean-floor metamorphic stage [11,13,14,60].These Na-Ca amphiboles have been dated by the Ar/Ar method between 50 and 40 Ma in partly amphibolitized eclogites (Monié et al. [62]; Puga et al. [13]).Other minerals of this same metamorphic event are the phengitic micas and the Fe-Mg amphiboles, which have been dated by Ar/Ar and K/Ar methods between 30 and 25 Ma [13, 61,62].This range of Meso-Alpine absolute ages include some of the U/Pb ages obtained on recrystallized zircon domains, especially in the Lugros eclogites (Puga et al. [15]).
Finally, the Neo-Alpine metamorphic event, developed mainly in greenschists facies conditions in all the units forming the Nevado-Filábride Complexes, has been dated beetwen 15 and 10 Ma. on different minerals and by different methods (Table 1).White micas and Fe-Mg amphiboles of amphibolitized eclogites and rodingites from several ophiolitic outcrops were dated by Ar-Ar and K-Ar as crystallized from 11 to 13 Ma.(Portugal Ferreira et al. [61]; Puga et al. [13], and zircon from a pyroxenite boudin, hosted by serpentinite from Almirez, was dated as 15 Ma.(López Sanchez-Vizcaíno et al. [63]).These later authors considered the dated zircons as formed during a Miocene high-pressure metamorphism meanwhile a discussion against this interpretation may be found in Puga et al. [15], pp.197-198.

Comparison of the Betic Ophiolites Radiometric Dating with Those of Alpine-Apenninic Ophiolites
The MOR-type ophiolites found at different locations in the Alps and Apennines (Figure 8) were dated using different radiometric techniques, of which the most reliable for inferring the age of magmatism are Sm/Nd and especially U/Pb applied to zircons from gabbros, complemented with dating of radiolarians preserved in the sedimentary sequences of some poorly metamorphosed ophiolites.Figure 8, by Bortolotti and Principi [66], shows the results of this dating and the techniques used on the ophiolites from each location.Some of the radiometric dating of Betic Ophiolites (shown in Table 1), coming from: Lugros, Almirez, Cóbdar and Algarrobo outcrops (Figure 1), which form part of the Ophiolite Unit of the MC, in the Nevado-Filábride Domain, have been also included in Figure 8 as Betic Cordillera for comparison.
As explained in the previous chapters, if we take into account only the dating carried out by U/Pb with SHRIMP on zircons, the mean absolute age for the magmatism originating from the Betic oceanic lithosphere ranges from 190 to 180 Ma, with an average age close to 185 Ma (Pliensbachian), for previously dated samples from Lugros and Almirez outcrops [15,19].These ages have been confirmed by the more recent dating of samples from Cóbdar and Algarrobo outcrops, shown in Figures 6 and 7 and Table 1 of this paper, whereas the SHRIMP U-Pb ages of gabbroic rocks from the Alps, Corsica and Apennines (with red color in Figure 8) yield comparatively younger ages, beginning at around 165 Ma.These Middle-Upper Jurassic ages for igneous zircons in gabbros and plagiogranites from the Alpine-Apennine Ophiolites (Figure 8) are also commonly accepted by other authors (Borsi et al. [67]; Costa and Caby [68]; Rubatto and Scambelluri [69]; Kaczmarek et al. [70]; Manatschal and Müntener [28]) and conform within the chronological interval 165-140 Ma.Thus, according to the SHRIMP zircon dating shown in Figures 6-8, and to the data compilation in Manatschal and Müntener [28], we can conclude that the break-up of Pangea and the development of the Western Tethys oceanic-floor began at the westernmost sector, from which the Betic Ophiolites derive (Betic Tethys) about 20 Ma earlier than in the northeastern extension (Ligurian and Alpine Tethys), from which the Alpine-Apennine Ophiolites originated.[66]), complemented by some representative zircons dating from the Betic Ophiolites, shown in Figures 6 and 7, and in Table 1.The comparable SHRIMP U/Pb dating on zircons from Betic and Alpine-Apennine ophiolitic gabbros are all represented with red color for more clarity.In green the Middle Jurassic age is represented, obtained for the oceanic floor metamorphism, applying Ar/Ar laser probe for dating millimetric veins of brown amphibole, which form part of the BOA basalts in the Cóbdar outcrop (Table 1).
As explained in the previous chapters, if we take into account only the dating carried out by U/Pb with SHRIMP on zircons, the mean absolute age for the magmatism originating from the Betic oceanic lithosphere ranges from 190 to 180 Ma, with an average age close to 185 Ma (Pliensbachian), for previously dated samples from Lugros and Almirez outcrops [15,19].These ages have been confirmed by the more recent dating of samples from Cóbdar and Algarrobo outcrops, shown in Figures 6 and 7 and Table 1 of this paper, whereas the SHRIMP U-Pb ages of gabbroic rocks from the Alps, Corsica and Apennines (with red color in Figure 8) yield comparatively younger ages, beginning at around 165 Ma.These Middle-Upper Jurassic ages for igneous zircons in gabbros and plagiogranites from the Alpine-Apennine Ophiolites (Figure 8) are also commonly accepted by other authors (Borsi et al. [67]; Costa and Caby [68]; Rubatto and Scambelluri [69]; Kaczmarek et al. [70]; Manatschal and Müntener [28]) and conform within the chronological interval 165-140 Ma.Thus, according to the SHRIMP zircon dating shown in Figures 6-8, and to the data compilation in Manatschal and Müntener [28], we can conclude that the break-up of Pangea and the development of the Western Tethys oceanic-floor began at the westernmost sector, from which the Betic Ophiolites derive (Betic Tethys) about 20 Ma earlier than in the northeastern extension (Ligurian and Alpine Tethys), from which the Alpine-Apennine Ophiolites originated.[66]), complemented by some representative zircons dating from the Betic Ophiolites, shown in Figures 6 and 7, and in Table 1.The comparable SHRIMP U/Pb dating on zircons from Betic and Alpine-Apennine ophiolitic gabbros are all represented with red color for more clarity.In green the Middle Jurassic age is represented, obtained for the oceanic floor metamorphism, applying Ar/Ar laser probe for dating millimetric veins of brown amphibole, which form part of the BOA basalts in the Cóbdar outcrop (Table 1).

Mesozoic Paleogeographic Reconstruction of the Betic Tethys
The new dating of zircons by SHRIMP for the Cóbdar and Algarrobo outcrops reinforce the hypothesis developed by Puga et al. [15,19], based on zircon dating from Lugros and Almirez outcrops, indicating that an important phase of formation of oceanic crust occurred between 190 and 180 Ma, all along the Betic Ocean domain (Figure 1).This suggest that the opening of the Betic sector of the neo-Tethys preceded that of the adjacent Ligurian and Alpine Tethys oceanic domains (Handy et al. [59], Bortolotti and Principi [66], Bill et al. [71]; Tribuzio et al. [72]; Rampone et al. [73]) of at least 20 Ma, as is shown in Figure 8.
This observation is important as it implies : (1) that the Betic Ophiolites would be the only preserved relics along the Alpine Chain deriving from the westernmost end of the Tethys, which is Pliensbachian in age as is shown in Figure 9a from Schettino and Turco [74]; and (2) that these ophiolites represent the birth of the westernmost Tethys Ocean and consequently the initial stage of the Pangaea break-up along the Betic Tethys, which was simultaneous to the opening of the Atlantic Ocean between 10 • N and 20 • N (Figure 9a).

Cretaceous Paleogeographic Reconstruction of the Betic Tethys and Eo-Alpine Metamorphic Evolution
Figure 10 shows the palaeogeographic reconstruction for the Early Cretaceous of the central-western area of the Mediterranean by Guerrera et al. [8], with slight modifications.In this figure, N-F Ocean indicates the area of oceanic floor where the BOA rocks, described in this paper, could have originated SE of the Iberian-European margin for the original paleogeographic continental domains of the tectonic units within which these ophiolites were intercalated.From West to East we have coloured the bands indicating the location of the following Mesozoic Betic Domains: the Subbetic margin, the VC, and the different original domains corresponding to the continental tectonic units forming part of the Mulhacén Complex: Caldera and Sabinas.These two latter units In Figure 9a the red-coloured bands represent the zones of fracture and crustal thinning (rifting) marking the break-up of the Pangaea supercontinent from the Triassic, where the Eurasian and African plates separated in this area.The first ocean floors of the Central Atlantic and Western Tethys, marked here in dark blue, began to develop following these rifting patterns approximately 185 Ma ago in the Pliensbachian (Lower Jurassic).The area marked BT (Betic Tethys) between the tectonic microplates located to the SE of Iberia-European plate during the Mesozoic, would correspond to the ridge where the Betic Ophiolites originated.This hypothesis is firmly supported by U/Pb radiometric dating with SHRIMP of the zircons separated from eclogitized gabbros and dolerites of the Betic Ophiolites, whose absolute magmatic correspond to the Pliensbachian (Puga et al. [15,19]; Figures 6 and 7, and Table 1 of this paper).The oceanic floors continued to develop throughout the Jurassic in both the Central Atlantic and the Western Tethys, as shown by the paleogeographic reconstruction of Figure 9b for the Tithonian (Late Jurassic).This stage of the paleogeograhic reconstruction agrees with some radiometric and petrologic evidence suggesting continuation of the spreading period of the Betic Tethys up to the Upper Jurassic or even to the Lower Cretaceous, as suggested by sedimentary evidences such as the development of an overlying sedimentary sequence, with cretaceous fossils relics [39], containing some sills with similar mineralogy and chemical composition to the underlying volcanic ophiolites.

Cretaceous Paleogeographic Reconstruction of the Betic Tethys and Eo-Alpine Metamorphic Evolution
Figure 10 shows the palaeogeographic reconstruction for the Early Cretaceous of the central-western area of the Mediterranean by Guerrera et al. [8], with slight modifications.In this figure, N-F Ocean indicates the area of oceanic floor where the BOA rocks, described in this paper, could have originated SE of the Iberian-European margin for the original paleogeographic continental domains of the tectonic units within which these ophiolites were intercalated.From West to East we have coloured the bands indicating the location of the following Mesozoic Betic Domains: the Subbetic margin, the VC, and the different original domains corresponding to the continental tectonic units forming part of the Mulhacén Complex: Caldera and Sabinas.These two latter units would correspond to both continental margins of the Betic Ocean, or Nevado-Filábride Ocean.The relative provenance of some Betic ophiolitic outcrops is indicated by the letters A to E, according to their current geographical location and differential petrologic and geochemical peculiarities.These characteristics, which are mainly shown in Figure 5D, suggest a probable origin of Cóbdar (D) and Algarrobo (E) outcrops closer to the western continental margin of the Jurassic Betic Ocean than those of Lugros (A), Almirez (B) and Caniles (C) outcrops.Moreover, eastwards the corresponding domains of the Alpujárride (AC) and Maláguide (MaC) Complexes are represented, as well as the emerged areas of the Mesomediterranean Terrane (MT) in the early Cretaceous.According to this reconstruction, at the start of the Cretaceous, the Betic Ocean could have connected the Piemontese-Ligurian (Pi-Li) and the Maghrebian (Mag) Oceans with the Central Atlantic Ocean.
The Jurassic extensional tectonic phase, which originated the opening of the Betic Ocean from the Pliensbachian onward, must have also affected other Betic Domains such as the Alpujárride and the Subbetic, whose relative location during the Mesozoic period is shown in Figure 10.This hypothesis is based on the tectonic setting and ages of the basic magmatism in each of these domains.This extensional tectonic phase affecting the Alpujárride Domain would explain the intrusion of some gabbros and basaltic lenses into the crustal rocks of the Ojen Nappe, identified as eclogites by Tubia et al. [75], which currently host the Ronda peridotites containing pyroxenite dykes.Some zircons from these eclogites and pyroxenites, dated by Sanchez Rodriguez and Gebauer [76], yielded ages from Pliesbachian to early Cretaceous, similar to the basic rocks forming the betic oceanic floor, although their originating magmas, also from asthenospheric origin, were scarcely developed, suggesting aborted rift tectonic conditions, and less effective extension with respect to what occurred in the the Nevado-Filábride domain during the same period.Despite these local tectonic differences, the coincidence of Mesozoic zircon crystallisation ages for two of the main domains of the Betic Cordillera, such as the Nevado-Filábride and the Alpujárride, suggests that in the Betic realm Pangaea actually began to break up in the Pliensbachian ( [19] and present paper), and not in the Middle Jurassic as generally reported in the literature (Principi et al. [77]; Bortolotti and Principi [66]; Manatschal and Müntener [28] and references herein).
Algarrobo (E) outcrops closer the western continental margin of the Jurassic Betic Ocean than those of Lugros (A), Almirez (B) and Caniles (C) outcrops.Moreover, eastwards the corresponding domains of the Alpujárride (AC) and Maláguide (MaC) Complexes are represented, as well as the emerged areas of the Mesomediterranean Terrane (MT) in the early Cretaceous.According to this reconstruction, at the start of the Cretaceous, the Betic Ocean could have connected the Piemontese-Ligurian (Pi-Li) and the Maghrebian (Mag) Oceans with the Central Atlantic Ocean.The Jurassic extensional tectonic phase, which originated the opening of the Betic Ocean from the Pliensbachian onward, must have also affected other Betic Domains such as the Alpujárride and the Subbetic, whose relative location during the Mesozoic period is shown in Figure 10.This hypothesis is based on the tectonic setting and ages of the basic magmatism in each of these domains.This extensional tectonic phase affecting the Alpujárride Domain would explain the intrusion of some gabbros and basaltic lenses into the crustal rocks of the Ojen Nappe, identified as eclogites by Tubia et al. [75], which currently host the Ronda peridotites containing pyroxenite dykes.Some zircons from these eclogites and pyroxenites, dated by Sanchez Rodriguez and Gebauer [76], yielded ages from Pliesbachian to early Cretaceous, similar to the basic rocks forming the betic oceanic floor, although their originating magmas, also from asthenospheric origin, were scarcely developed, suggesting aborted rift tectonic conditions, and less effective extension with respect to what occurred in the the Nevado-Filábride domain during the same period.Despite these local Extensional tectonics during the Jurassic affected also the Subbetic Domain along which a fissural basic magmatism developed, deriving from a lithospheric mantle.This magmatism was dated by K/Ar between the Lower Jurassic up to Lower Cretaceous (Puga [78]; Puga and Portugal Ferrerira [79]; Puga et al. [80][81][82]; Portugal Ferreira et al. [83]).Some of these ages have been interpreted as partial resetting of the original igneous ages during the alpine metamorphism in prehnite-pumpellyite and pumpellyite-actinolite facies, which was developed mainly in dolerite bodies intruded in Triassic formations, especially for those closer to the Betic Internal Zones, which collided with the External Zones during the Miocene (Morten and Puga [84]; Puga et al. [85,86]; Aguirre et al. [87]).
According to Guerrera et al.
( [8] and references herein) the Nevado-Filábride and Ligurian-Piemontese Oceans as well as the northern margin of the Mesomediterranean microplate were structured from Late Cretaceous to the Eocene to create an Eo-Alpine Chain The U-Pb data of zircons from the Betic Ophiolites presented in Puga et al. [15,19] and in Figures 6 and 7, and Table 1 of this paper, showing Upper Cretaceous to Paleocene recrystallization domains on the Jurassic igneous zircon crystals, support the development of this first Alpine tectonic and metamorphic event, which originated the transformation of the basic and ultramafic rocks forming the oceanic-floor, into eclogites, during the Eo-Alpine subduction process as is shown in Figure 11.Eo-Alpine eclogites have been also identified in many other ophiolite outcrops from the Alps, where they have been interpreted as originated by Upper Cretaceous-Paleocene intra-oceanic subduction of the Ligurian sections of the Neothetys, located to the East of Iberian and the west of Adrian microplates (Handy et al. [59] and references herein).

Neogene Geodynamic Evolution the Mesomediterranean Terrane
Subsequently to the Eo-Alpine reestructuration of the Mesomediterranean Terrane, mainly affecting its northern margin, took place, during the late Oligocene and earliest Miocene, the dislaceration of this microplate into the different blocks collectively known, according to Bouillin et al. [88], as AlKaPeCa (Alboran, Kabylian, Peloritan and Calabrian).These blocks were displaced toward the south to attain their respective present locations due to subduction toward the NW and roll-back under them of different oceanic plates, genetically related with the development of the new Mediterranean oceanic basins (Handy et al. [59]; Chalouan and Michard [89]; Lonergan and White [90], Rosebaum et al. [91]; Lustrino et al. [92]).In particular, the Alboran block, or microplate, was displaced to the SW to occupy its present setting between the Iberian and African plates mainly due to the Neogene opening of the Algero-Provençal Basin.During this translation, the different domains of the Alboran microplate, forming the Internal Zones of the Betic and Rifean Cordilleras, collided with their External Zones belonging to the South Iberian and African margins and, together with them, underwent the later alpine stages of deformation and metamorphism.These are characterized by local differences of P-T conditions between the different domains of the Internal Betic Zones, probably due to differential ascent of the asthenosphere below them during the Miocene.These Neogene metamorphic processes, known as Meso-Alpine and Neo-Alpine events, developed in the Nevado-Filábride Domain in Ab-Ep amphibolite and greenschist facies conditions, respectively, overprinting the Eo-Alpine eclogites that were partly retrograded and tectonically deformed.

Genetic and Evolutionary History of the Betic Ophiolites and the Nevado-Filábride Complexes
Figure 11 shows the origin and evolution of the Betic Ophiolites based on the study of the successive parageneses and textures preserved in them and their thermodynamic formation conditions, together with the radiometric dating of some of their minerals, particularly zircons, and the geochemistry of their different lithologies.Further information is provided by data relative to the tectonic units in which these ophiolites are currently intercalated.This model, based on the previously exposed data in this paper and on other data previously published ; Nieto [41]; Nieto et al. [42]; Ruiz Cruz et al. [20,21]; Aerden et al. [23]), shows the schematic geodynamic evolution of the VC and MC from the Carboniferous to the Miocene.It also includes the evolution of the various units comprising the latter complex (MC), one of which is the Ophiolite Unit (Figure 2).Development in time has been divided into several episodes in Figure 11, numbered 1 to 6, to which correspond successive geological stages affecting both the VC and MC, which are graphically shown and enumerated on the right of the diagram.
The oceanic lithosphere from which the Betic Ophiolites derived was formed in episode 3 by asthenospheric rise, under extensional conditions that persisted throughout the Jurassic and Cretaceous, along a rifting zone on the continental crust of the MC, located between the Caldera (W) and Sabinas (E) units.Among other Pre-Alpine metamorphic rocks, this continental crust contained orthogneisses derived from Hercynian syn-collisional granits (red color) and Permian-Triassic intra-plate rhyolites (blue color) formed earlier in episodes 1,2.In episode 4, the Jurassic-Cretaceous oceanic floor, and part of its continental margins, were subducted toward the SE as a result of a first convergent movement between of the Iberian-European and African-Adria plates beginning in the Late Cretaceous, which caused the Eo-Alpine metamorphism in eclogite facies conditions (Table 1).Genetically related to this subduction episode there are only scarce magmatic products, mainly represented by pyroclastic andesites which, melanged with continental carbonates and gypsum deposits, would gave origin to the Soportújar Formation (SF) ([13] and Figure 2).During episode 5 a small part of the subducted rocks, metamorphosed at depths from 50 up to more than 100 km ([10,11] and Figure 4), were exhumed as tectonic slabs, onto the Caldera domain (its western continental margin), forming the Ophiolite Unit, which was also overthrusted by the Sabinas unit (its eastern continental margin) and by the SF, forming an initial stacking involving the three units of the MC and probably also of the MC onto the VC.During this compressive episode, the Meso-Alpine metamorphic event developed, affecting both complexes in Ab-Ep amphibolite facies conditions during Paleogene (Table 1).Finally, episode 6 took place along the Neogene, during which the Alboran microplate, dismembered from the western part of the Mesomediterranean Terrane, was displaced to the SW, where the metamorphic Internal Zones emplaced on top of the External Zones of the Betic and Rift domains and underwent, together with them, the continental collision between the Iberian and African plates.This collision originatedthe Neo-Alpine metamorphic event, (mainly developed in greenschist facies conditions in the MC and VC), which was followed by the overstacking of the Alpujárride and Maláguide onto the previously empiled Nevado-Filábride Complexes.
Geosciences 2017, 7, 31 20 of 29 conditions, together with the radiometric dating of some of their minerals, particularly zircons, and the geochemistry of their different lithologies.Further information is provided by data relative to the tectonic units in which these ophiolites are currently intercalated.This model, based on the previously exposed data in this paper and on other data previously published ; Nieto [41]; Nieto et al. [42]; Ruiz Cruz et al. [20,21]; Aerden et al. [23]), shows the schematic geodynamic evolution of the VC and MC from the Carboniferous to the Miocene.It also includes the evolution of the various units comprising the latter complex (MC), one of which is the Ophiolite Unit (Figure 2).Development in time has been divided into several episodes in Figure 11, numbered 1 to 6, to which correspond successive geological stages affecting both the VC and MC, which are graphically shown and enumerated on the right of the diagram.The oceanic lithosphere from which the Betic Ophiolites derived was formed in episode 3 by asthenospheric rise, under extensional conditions that persisted throughout the Jurassic and Cretaceous, along a rifting zone on the continental crust of the MC, located between the Caldera (W) and Sabinas (E) units.Among other Pre-Alpine metamorphic rocks, this continental crust contained orthogneisses derived from Hercynian syn-collisional granits (red color) and Permian-Triassic intra-plate rhyolites (blue color) formed earlier in episodes 1,2.In episode 4, the Jurassic-Cretaceous oceanic floor, and part of its continental margins, were subducted toward the SE as a result of a first
Palinspastic restoration of BOA rocks is not simple and necessarily requires the tentative back-stripping and positioning of the adjoining Alboran domain and continental blocks such as those currently outcropping in crystalline massifs of Kabylie (Algeria) Calabria and Peloritani (Italy) that were to the BOA sections before (and during) the neo-Tethys opening, probably forming a unique microplate known as Mesomediterranean Terrane (Guerrera et al. [8]).This microplate, located to the SE margin of the Iberia-European plate, underwent the effects of the Eo-Alpine Upper Cretaceous Orogeny, by convergence of the Iberia-European and Adria-African plates, with subduction toward the south of the Betic oceanic-floor, followed by Paleocene exhumation of some eclogitized relics on the northern continental margin.This Mesomediterranean microplate was subsequently dismembered into its different continental blocks, collectively known as AlKaPeCa (Alboran, Kabylies, Peloritani, Calabria) [88], from the Oligocene onward due to subduction towards the NW under this terrane of a new oceanic microplate (Guerrera et al. [8]; Handy et al. [59]; Chalouan and Michard [89]; Lonergan and White [90]; Rosebaum et al. [91]; Lustrino et al. [92]; Perez-Lopez and Perez Valera [94]; Critelli et al. [95]).
Consistently with this interpretation, the BOA rocks and those of the surrounded crystalline basements, which formed the continental margins of the Betic Ocean, were involved in the oceanic subduction and the following collision of major continental plates such as Iberia-Europe and Africa, recording similar P-T-t metamorphic paths during the Eo-Alpine orogenic phase that led to the formation of eclogite facies in the MC.This first alpine metamorphic event was followed during the Neogene by the Meso-Alpine and Neo-Alpine orogenic events, which produced the re-crystallisation in the Ab-Ep amphibolite and greenschist facies conditions of the Eo-Alpine eclogites.
The remarkable petrological analogy between the BOA and Alpine-Apennine ophiolites has to be emphasized; it indicates that during the Jurassic they were part of the same oceanic system elongated in a NE-SW direction, that in our view, was westward connected with the embryonic Atlantic Ocean (see Figures 9 and 10).This ocean was unrelated to the oceanic systems represented by the Dinarides and Hellenides ophiolites, which developed eastward and concomitant with active subduction processes.This hypothesis is not based only on ophiolite investigations, but is supported by facies analyses and paleontological studies of coeval sedimentary sequences (Tendero et al. [39]; Lemoine et al. [43]; Lagabrielle et al. [44]; Principi et al. [77]; Perez-Lopez and Perez-Valera [94]; Critelli et al. [95]; Hanne et al. [96]; Martin-Rojas et al. [97]; Somma et al. [98]; Zaghloul et al. [99]).
The described scenario is well represented in the paleogeograhic reconstruction of the Western Tethys of Figure 9, which is perfectly in agreement with the model proposed by Schettino and Turco [74].According to our radiometric zircon dating on the BOA rocks (Table 1), the Betic oceanic floor initiated its development in the Pliensbachian (from about 190 to 180 Ma), and continued its metasomatic and metamorphic transformations, under oceanic floor conditions, affecting most of the constituent lithologies, from the ultramafic sequence up to the volcanic one, from the Toarcian to the Tithonian (about 177 to 152 Ma for different outcrops).This was followed by a late magmatic activity originating basaltic sills and dykes having chemical affinity with that of the previous MORB magmas forming the oceanic volcanic and plutonic underlying sequences.These sills are interlayered in the overlying sedimentary sequence, which could be Cretaceous in age (Tendero et al. [39] and references herein) and the scarce dykes cut across some pillowed volcanic layers and more rarely the overlying metasediments.The radiometric dating presented in Table 1, together with the paleogeographic reconstruction of Schettino and Turco [74] for the westernmost segment of the Neotethys (Figure 9), suggest the development of the Betic Ocean for a period of at least 40-50 Ma.On the other hand, the eastward extension of the Betic Tethys Ocean along the Ligurian and Alpine Tethys, that took place from the Middle Jurassic onward, beginning at about 165 Ma, developed up to the Upper Tithonian for about 20-30 Ma, according to SHRIMP dating of zircons in the Western and Ligurian Alps and the radiolarians dating of the sedimentary sequence from Tuscan Apennines (see Figures 8 and 9).
The Betic Ocean and its prolongation in the Piemonte-Ligurian domain evolved in slow spreading mode and they were dissected by many transform zones that displaced the ridges settings, leaving significant continental blocks isolated from the adjoining macro-plates.
Major plates re-organization caused cessation of oceanic spreading and incipient consumption of the Neo-Tethys since the Late Cretaceous.As observed in the Alps-Apennine systems, and also in the Betic realm, the subduction developed during several stages characterized by changes in polarity (Lonergan and White [90]; Martin-Martin et al. [93]; Carminati et al. [100]; Verges and Fernandez [101]).The Eo-Alpine subduction, active during Late Cretaceous-Paleocene, was dipping southeastward and plausibly continued up to the continental collision between the south Iberian-European margin and the Mesomediterranean microplate, ultimately leading to obduction of some BOA eclogitized sequences on the western continental margin of the Betic Ocean, formed by the present Caldera Unit of the MC.Then, plate convergence reactivated during Eocene with a reversed polarity characterized by active subduction with remarkable slab passive sinking and roll-back, which allowed extension on the Betic margin originating exhumation of deep mantle rocks and favouring, from Oligocene to Late Miocene, the development of tholeiitic dykes, was followed by extensive calcalkaline and ultrapotassic magmatism (Carminati et al. [100]; Verges and Fernandez [101]; Torres-Roldan et al. [102]; Doblas et al. [103]).The concomitant back arc extension occurring in this Neogene tectonomagmatic phase dismembered the Mesomediterranean microplate, whose remnants progressively migrated to their current settings (Alboran, Kabilies, Peloritani massif and Calabria), mainly due to the opening of the oceanic Algerian-Provençal basin and Tyrrhenian Sea.

Conclusions
The BOA and the related lithologies, mainly acidic gneissic levels in different units of the MC, dated from 278 to 215 Ma ([13,37,41,42]) testify that the Pangea continent was affected by pervasive rifting accompanied by acidic volcanism since the Permian/Triassic, which preceded the continental breakup and development of the Jurassic oceanic basin.Rifting and oceanic drifting was plausibly influenced by pre-existing tectonic structures inherited by the Varisic orogeny (Manatschal et al. [104]; Gilard et al. [105]), giving rise to a very complex paleogeography with many ridge offset and short zigzags accommodating seafloor accretion as observed in modern oceanic sectors characterized by slow spreading (Grindlay et al. [106]; Lazar et al. [107]; Kamesh Raju et al. [108]; Murton and Rona, [109]; Mercier de Lèpinay et al. [110]).
In this framework, sharing similar features, the Betic ocean was connected eastward with other oceanic strands currently representing the Liguride ophiolites in the Alpine and Apennine domains, and westward with the embryonic Atlantic Ocean.
According to most recent interpretations, this oceanic corridor having limited width (few hundreds of kilometers) and dissected by many transform zones, was punctuated by Oceanic Core Complexes (OCC) that are dome shaped, with corrugated features, developed along detachment faults (Lagabrielle et al. [31]; Lagabrielle [33]; Manatschal et al. [111]; Li et al. [112]), alternating magmatic sectors characterized by the outcrop of pillow basalts, and amagmatic sectors characterized by the exposure on the seafloor of ultramafic rocks, sometimes intruded by gabbros or dolerite dykes.These associations rarely fit the "Penrose three layer" model constituted by the regular coexistence of ultramafic rocks, basic magmatic rocks and sediments, showing a particular paucity of sheeted dyke complexes.On the contrary, Ligurian Ophiolites are often associated with remnants of the continental crust, and analogously the Betic Ophiolites locally present some indices of crustal assimilation, probably deriving from basement rocks of their continental margin, indicating in both cases (1) vicinity of the continental margin at least during the first stages of the oceanic-floor magmatism, and (2) intermittent spreading generally attributed to ultra-slow spreading ridges [19] or to OCT setting [30].
These tectonomagmatic features, favouring the persistence of isolated continental blocks within the oceanic domain, can explain the close association of the BOA and Ligurian oceanic rocks with continental counterparts, most of which were agglomerated within the microplate known as Mesomediterranean Terrane from which, after the Eo-Alpine orogeny, the AlKaPeCa blocks derived.
The Triassic-Jurassic paleogeography delineated above, obviously influenced the following convergence processes that initiated since the Late Cretaceous following major plate reorganization driven by opening of the South Atlantic basin and anticlockwise rotation of the African plate.The subduction during the Eo-Alpine cycle was dipping to the SE and S direction in the Betic and Alpine sectors, These subduction processes became progressively blocked in ca.20 Ma, probably in relation to continental collision between the above mentioned major plates (Iberia, Europe, Africa) and the interposed Mesomediterranean and Adria microplates, leading to obducion of the Betic and Liguride ophiolite belts.As indicated by the regional distribution of subduction-related volcanism (Wilson and Bianchini [113]; Beccaluva et al. [114][115][116]; Bianchini et al. [117]) and by geophysical studies (Lucente et al. [118]; Vignaroli et al. [119]; Mancilla et al. [120]) that provide a snapshot of the current geodynamic configuration, further complexity of the Betic and Alpine-Apennine orogenic systems was subsequently induced by Cenozoic subduction processes that were characterized by a flip of polarity, with active subduction dipping toward the N-NW in the Apennine sector, Northward in Maghrebian sector and NE in the circum-Alboran realm.

Figure 1 .
Figure 1.Geological sketch map of the central-eastern sector of the Betic Cordilleras, showing the relationships between the Veleta, Mulhacén and Alpujárride metamorphic complexes and the post-orogenic volcanic and sedimentary formations, modified from Figure 1 in Puga et al. [19].The upper inset shows a simplified version of the present tectonic relations of the Nevado-Filábride complexes, with Veleta (VC) beneath and Mulhacén (MC) on top, on which the Alpujárride (AC) and Maláguide (MaC) complexes and the External Betic Zones (EZ) are superposed.Black four-pointed stars represent the largest outcrops of the BOA located along the MC of the Nevado-Filábride Domain.Red encircled stars indicate the location of the main ophiolitic outcrops from which most of the studied rocks, considered in this revision, derive.The studied outcrops are the following sites: Lugros, Camarate, Loma Pelada, Almirez, Santillana and Montenegro (in Sierra Nevada), Cóbdar and Lubrín (Sierra de Filabres), Algarrobo (Sierra del Algarrobo) and Caniles (Sierra de Baza).Letters V and M in blue, in the western area of the map, indicate the approximate geographic setting of the Veleta and Mulhacén Pics.

Figure 1 .
Figure 1.Geological sketch map of the central-eastern sector of the Betic Cordilleras, showing the relationships between the Veleta, Mulhacén and Alpujárride metamorphic complexes and the post-orogenic volcanic and sedimentary formations, modified from Figure 1 in Puga et al. [19].The upper inset shows a simplified version of the present tectonic relations of the Nevado-Filábride complexes, with Veleta (VC) beneath and Mulhacén (MC) on top, on which the Alpujárride (AC) and Maláguide (MaC) complexes and the External Betic Zones (EZ) are superposed.Black four-pointed stars represent the largest outcrops of the BOA located along the MC of the Nevado-Filábride Domain.Red encircled stars indicate the location of the main ophiolitic outcrops from which most of the studied rocks, considered in this revision, derive.The studied outcrops are the following sites: Lugros, Camarate, Loma Pelada, Almirez, Santillana and Montenegro (in Sierra Nevada), Cóbdar and Lubrín (Sierra de Filabres), Algarrobo (Sierra del Algarrobo) and Caniles (Sierra de Baza).Letters V and M in blue, in the western area of the map, indicate the approximate geographic setting of the Veleta and Mulhacén Pics.

Figure 2 .
Figure 2. Simplified lithostratigraphic column showing the main rock types of the different tectonic units and formations making up the Veleta and Mulhacén Complexes, together with the lithostratigraphy (not to scale) of the overlying Alpujárride Complex, modified from Puga et al. [19].The main lithologies forming the different tectonic units and formations are represented in Figure 2 by different screens and colors, whose meaning can be seen in the legend of this figure.Abbreviations on the right border of Figure 2 have the following meaning, Pz: Palaeozoic, PTr: Permian-Triassic, T: Triassic, J: Jurassic, Cr: Cretaceous, Pl: Paleocene, and indicate the radiometrically dated ages, or the geologically more probable ages.

Figure 2 .
Figure 2. Simplified lithostratigraphic column showing the main rock types of the different tectonic units and formations making up the Veleta and Mulhacén Complexes, together with the lithostratigraphy (not to scale) of the overlying Alpujárride Complex, modified from Puga et al. [19].The main lithologies forming the different tectonic units and formations are represented in Figure 2 by different screens and colors, whose meaning can be seen in the legend of this figure.Abbreviations on the right border of Figure 2 have the following meaning, Pz: Palaeozoic, PTr: Permian-Triassic, T: Triassic, J: Jurassic, Cr: Cretaceous, Pl: Paleocene, and indicate the radiometrically dated ages, or the geologically more probable ages.

Figure 4 .
Figure 4. P-T diagram showing the stability fields of the different metamorphic events registered in the BOA (basic and ultramafic rocks) compared with other units of Veleta and Mulhacén complexes, modified from Puga et al. [13].Meaning of abbreviations indicating these fields are V and brown color: Veleta; M and green or blue color: Mulhacén; P: Pre-Alpine; 1: Eo-Alpine; 2: Meso-Alpine and Neo-Alpine; C: Caldera Unit; S: Sabinas Unit; Ou: Ultramafic rocks of the Ophiolite Unit; Ob: Basic rocks of the Ophiolite Unit.The schematic P-T-t path in the inset of this figure shows the climax conditions corresponding to the Pre-Alpine, Eo-Alpine and Meso-Alpine events represented, respectively, by the position of Mp, M1 and M2 for the MC, and by those of Vp, V1 and V2 for the VC.The position of the arrow with the label "Oceanic metamorphism" indicates the range of low pressure and high temperature conditions deduced for this process in the Cóbdar metabasalts (Puga et al. [14]).Labeled metamorphic reactions as in Puga et al. [13].

Figure 4 .
Figure 4. P-T diagram showing the stability fields of the different metamorphic events registered in the BOA (basic and ultramafic rocks) compared with other units of Veleta and Mulhacén complexes, modified from Puga et al. [13].Meaning of abbreviations indicating these fields are V and brown color: Veleta; M and green or blue color: Mulhacén; P: Pre-Alpine; 1: Eo-Alpine; 2: Meso-Alpine and Neo-Alpine; C: Caldera Unit; S: Sabinas Unit; Ou: Ultramafic rocks of the Ophiolite Unit; Ob: Basic rocks of the Ophiolite Unit.The schematic P-T-t path in the inset of this figure shows the climax conditions corresponding to the Pre-Alpine, Eo-Alpine and Meso-Alpine events represented, respectively, by the position of Mp, M1 and M2 for the MC, and by those of Vp, V1 and V2 for the VC.The position of the arrow with the label "Oceanic metamorphism" indicates the range of low pressure and high temperature conditions deduced for this process in the Cóbdar metabasalts (Puga et al. [14]).Labeled metamorphic reactions as in Puga et al. [13].

Figure 5 .
Figure 5. Distribution of scarcely mobile trace elements (A), or their ratios (D), REE values normalized to chondrites (C) and isotopic ratios of Sr vs. Nd (B) of the different lithologies forming the Betic Ophiolites compared with MORB, OIB and IAT basic magmas, modified from [19] in order to include the following outcrops: Cóbdar, Lugros, Algarrobo and Almirez.(A) Ti/1000 vs. V diagram (Shervais [45]); (B) Nb/Yb vs. Th/Yb diagram (Pearce [46]); (C) Chondrite-normalised REE diagram for basic and ultramafic rocks, using normalising values of Boynton [47]; (D) 87Sr/86Sr vs. 143Nd/144Nd isotope correlation diagram for basic and ultramafic rocks of the Betic Ophiolites, compared with abyssal peridotites (data from Snow et al. [48]); dashed lines showing mantle-seawater mixing curves for the abyssal peridotites are also reported by Snow et al. [48]); Sr and Nd isotopic values for EPR (Eastern Pacific Rise), MAR (Mid-Atlantic Ridge), SWIR (South Western Indian Ridge) and OIB values (Oceanic Island basalts) from Hoffman [49] ; CFB values (Continental flood basalts) from Marzoli et al. [50].Note that in Figure 5A,B,D, are also plotted, with black color, the rocks radiometrically dated, deriving from each of the four selected outcrops, to show the chemical similarity between the dated rocks and the other rocks from the same outcrops.

Figure 5 .
Figure 5. Distribution of scarcely mobile trace elements (A), or their ratios (D), REE values normalized to chondrites (C) and isotopic ratios of Sr vs. Nd (B) of the different lithologies forming the Betic Ophiolites compared with MORB, OIB and IAT basic magmas, modified from [19] in order to include the following outcrops: Cóbdar, Lugros, Algarrobo and Almirez.(A) Ti/1000 vs. V diagram (Shervais [45]); (B) Nb/Yb vs. Th/Yb diagram (Pearce [46]); (C) Chondrite-normalised REE diagram for basic and ultramafic rocks, using normalising values of Boynton [47]; (D) 87Sr/86Sr vs. 143Nd/144Nd isotope correlation diagram for basic and ultramafic rocks of the Betic Ophiolites, compared with abyssal peridotites (data from Snow et al. [48]); dashed lines showing mantle-seawater mixing curves for the abyssal peridotites are also reported by Snow et al. [48]); Sr and Nd isotopic values for EPR (Eastern Pacific Rise), MAR (Mid-Atlantic Ridge), SWIR (South Western Indian Ridge) and OIB values (Oceanic Island basalts) from Hoffman [49] ; CFB values (Continental flood basalts) from Marzoli et al. [50].Note that in Figure 5A,B,D, are also plotted, with black color, the rocks radiometrically dated, deriving from each of the four selected outcrops, to show the chemical similarity between the dated rocks and the other rocks from the same outcrops.

Figure 6 .
Figure 6.Cathodoluminescence images of selected zircons of BOA eclogitized gabbros and dolerites from the representative outcrops of Lugros, Cóbdar, Algarrobo and Almirez dated by U/Pb with SHRIMP (Sensitive High Resolution Ion Micro-Probe).

Figure 6 .
Figure 6.Cathodoluminescence images of selected zircons of BOA eclogitized gabbros and dolerites from the representative outcrops of Lugros, Cóbdar, Algarrobo and Almirez dated by U/Pb with SHRIMP (Sensitive High Resolution Ion Micro-Probe).

Figure 7 .
Figure 7. Tera-Wasserburg diagrams and probability histograms with mean ages for the dated zircons of four samples from Lugros (A), and Almirez (D) outcrops, according to [19], complemented with samples from Cóbdar (B) and Algarrobo (C).

Figure 7 .
Figure 7. Tera-Wasserburg diagrams and probability histograms with mean ages for the dated zircons of four samples from Lugros (A), and Almirez (D) outcrops, according to [19], complemented with samples from Cóbdar (B) and Algarrobo (C).

Figure 8 .
Figure 8.Comparison of igneous ages between Alpine-Apennine Ophiolites and the BeticOphiolites from Lugros and Almirez[19], together with those from Cóbdar and Algarrobo outcrops added in this paper.This figure includes the chronological data obtained by different methods on the Alpine-Apennine Ophiolites, deriving from the Jurassic Western Tethys (according to Bortolotti and Principi[66]), complemented by some representative zircons dating from the Betic Ophiolites, shown in Figures6 and 7, and in Table1.The comparable SHRIMP U/Pb dating on zircons from Betic and Alpine-Apennine ophiolitic gabbros are all represented with red color for more clarity.In green the Middle Jurassic age is represented, obtained for the oceanic floor metamorphism, applying Ar/Ar laser probe for dating millimetric veins of brown amphibole, which form part of the BOA basalts in the Cóbdar outcrop (Table1).

Figure 8 .
Figure 8.Comparison of igneous ages between Alpine-Apennine Ophiolites and the Betic Ophiolites from Lugros and Almirez [19], together with those from Cóbdar and Algarrobo outcrops added in this paper.This figure includes the chronological data obtained by different methods on the Alpine-Apennine Ophiolites, deriving from the Jurassic Western Tethys (according to Bortolotti and Principi[66]), complemented by some representative zircons dating from the Betic Ophiolites, shown in Figures6 and 7, and in Table1.The comparable SHRIMP U/Pb dating on zircons from Betic and Alpine-Apennine ophiolitic gabbros are all represented with red color for more clarity.In green the Middle Jurassic age is represented, obtained for the oceanic floor metamorphism, applying Ar/Ar laser probe for dating millimetric veins of brown amphibole, which form part of the BOA basalts in the Cóbdar outcrop (Table1).

Figure 9 .
Figure 9. Paleogeographic reconstruction of the Western Tethys and Central Atlantic from Schettino and Turco [74] at 185 Ma (Pliensbachian) (a).and 147.7 Ma (chron 21, Tithonian) (b).The strips represented in red in this Figure, along some zones of the Atlantic, Ligurian and Alpine future oceans, corresponded during the Pliensbachian to rifting zones, that evolved to continental break-up that preceded oceanic floor formation, from the Tithonian onward.In (a) we have indicated the probable location of the Betic Tethys (BT), on the initial oceanic-floor branch of the Western Tethys (represented in dark blue color), according with the Pliensbachian U/Pb zircon ages of the Betic ophiolitic gabbros of Lugros and Almirez outcrops [19] complemented, in this paper, by similar dating from Cóbdar and Algarrobo outcrops.In (b) are represented the synthetic isochrons relative to the Tithonian opening of the Tethys oceans, modified by us to locate, in this paper, the Betic Ocean to the western of the Ligurian and Alpine Tethys.GiF: Gibraltar Fault; NPF: North Pyrennean Fault.

Figure 9 .
Figure 9. Paleogeographic reconstruction of the Western Tethys and Central Atlantic from Schettino and Turco [74] at 185 Ma (Pliensbachian) (a).and 147.7 Ma (chron 21, Tithonian) (b).The strips represented in red in this Figure, along some zones of the Atlantic, Ligurian and Alpine future oceans, corresponded during the Pliensbachian to rifting zones, that evolved to continental break-up that preceded oceanic floor formation, from the Tithonian onward.In (a) we have indicated the probable location of the Betic Tethys (BT), on the initial oceanic-floor branch of the Western Tethys (represented in dark blue color), according with the Pliensbachian U/Pb zircon ages of the Betic ophiolitic gabbros of Lugros and Almirez outcrops [19] complemented, in this paper, by similar dating from Cóbdar and Algarrobo outcrops.In (b) are represented the synthetic isochrons relative to the Tithonian opening of the Tethys oceans, modified by us to locate, in this paper, the Betic Ocean to the western of the Ligurian and Alpine Tethys.GiF: Gibraltar Fault; NPF: North Pyrennean Fault.

Figure 10 .
Figure 10.Early Cretaceous paleogeography of the Central-Western Mediterranean region and distribution of the oceanic basins around the Mesomediterranean Terrane (M.T.), mainly the Nevado-Filábride (NF) Ocean, connecting the Piemontese-Liguride (Pi-Li) and the Maghrebian (Mag) oceans with the central Atlantic Ocean (Guerrera et al. [8]).Colored strips represent the different Betic domains, from which the Veleta and Mulhacén complexes derive, which surround the Jurassic-Cretaceous NF Ocean, or Betic Ocean [19].Letters A, B, C, D and E, indicate the relative provenance of Lugros, Almirez, Caniles, Cóbdar and Algarrobo ophiolitic outcrops along the Betic Ocean, according to their current geographical locations and differential geochemical characteristics.Key for other abbreviations in this figure can be seen in Guerrera et al. [8].

Figure 10 .
Figure 10.Early Cretaceous paleogeography of the Central-Western Mediterranean region and distribution of the oceanic basins around the Mesomediterranean Terrane (M.T.), mainly the Nevado-Filábride (NF) Ocean, connecting the Piemontese-Liguride (Pi-Li) and the Maghrebian (Mag) oceans with the central Atlantic Ocean (Guerrera et al. [8]).Colored strips represent the different Betic domains, from which the Veleta and Mulhacén complexes derive, which surround the Jurassic-Cretaceous NF Ocean, or Betic Ocean [19].Letters A, B, C, D and E, indicate the relative provenance of Lugros, Almirez, Caniles, Cóbdar and Algarrobo ophiolitic outcrops along the Betic Ocean, according to their current geographical locations and differential geochemical characteristics.Key for other abbreviations in this figure can be seen in Guerrera et al. [8].

Figure 11 .
Figure 11.Genetic and evolutionary pattern of the Betic Ophiolites and the continental units of the Nevado-Filábride Complexes, among which they are tectonically intercalated, and Paleogene overthrusting of the Mulhacén Complex onto the Veleta Complex.

Figure 11 .
Figure 11.Genetic and evolutionary pattern of the Betic Ophiolites and the continental units of the Nevado-Filábride Complexes, among which they are tectonically intercalated, and Paleogene overthrusting of the Mulhacén Complex onto the Veleta Complex.

Table 1 .
Radiometric dating of the different genetic processes originating from the Betic Ophiolites.