Geochemical characterization and isotopic age of Caradocian magmatism in the northeastern Iberian Peninsula: Insights into the Late Ordovician evolution of the northern Gondwana margin Marina Navidad a,*, Pedro Castiiieiras a, Josep M. Casas b, Montserrat Liesa d, Javier Fernandez Suarez a, Antoni Barnolas c, Jordi Carreras e, lnmaculada Gil-Peiia C • Departamento de Petro[ogfa y Geoqufmica-Instituto de Geologfa Econ6mica (CSIC) Universidad Complutense de Madrid,jose Antonio Novais 2, Madrid, 28040, Spain b Departament de Geodinamica i Geofisica-Institut de recerca GEOMODELS, Universitat de Barcelona, Martf i Franques s!n, Barcelona, 08028, Spain c Instituto Ge%gico y Minero de Espaiia, Rlos Rosas 23, Madrid, 28003, Spain d Departament de Geoqufmica, Petr%gia i Prospecci6 Ge%gica, Universitat de Barcelona, Mart! i Franques s!n. Barcelona, 08028, Spain • Departament de Geologia, Universitat AutOnoma de Barcelona, Bellaterra (Cerdanyo/a del Valles), 08193, Spain Keywords: Metaluminous and aluminous orthogneisses Pyroc1astic cale-alkaline metavoleanic rocks GeochemistlY Sr-Nd isotopes U-Pb dating Llte Ordovician magmatism Northeastern Iberia 1. Introduction ABSTRACT This paper presents new geochronological, geochemical and isotopic data for Late Ordovician magmatism (455 ± 2 Ma) from the Variscan massifs in the northeastern Iberian Peninsula. This magmatism is mainly composed of aluminous and metaluminous plutonic rocks and a sequence of calc-alkaline ignimbrites and epiclastic volcanic rocks, which were deformed and metamorphosed during the Variscan orogeny. The metaplutonic rocks are located in the lower part of a NeoproterozoicjEarly Cambrian metasedimentary sequence in the Canigo massif (Eastern Pyrenees). Aluminous orthogneisses have isotopic signatures that indicate a crustal origin, whereas metaluminous orthogneisses originate from a mixture of juvenile and crustal melts. The metavolcanic rocks define a cale-alkaline ignimbritic sequence of crustal origin located in the lower part of the Late Ordovician sequence in the Les Cavarres massif (Catalan Coastal Ranges). TDM ages for all the studied Late Ordovician rocks range from Mesoproterozoic to Paleoproterozoic (13 Ca, 1.5 Ca, 1.7 Ca and 2.2 Ca), suggesting a Neoproterozoic crust with Paleoproterozoic components and an arc signature, as the most probable source for the crustal melts. The isotopic similarities between the studied samples and other Early Ordovician magmatic rocks of the Central and Western Iberian Massif suggest repeated melt extractions from a common basement as the source for the Early and Late Ordovician magmatism in the Iberian realm. However, the studied massifs are more closely linked to other massifs from the Mediterranean realm. The most probable tectonic setting for the emplacement of the Late Ordovician volcanic and plutonic bodies is that of an extensional regime, which postdates a previous Ordovician contractional event and the opening of the Rheic Ocean. Variscan rocks that crop out in northeastern Iberia are concen­ trated in two areas: the Pyrenees and the Catalan Coastal Ranges (Fig. 1). The pre-Silurian succession of these areas consists of an up to 5000 m thick sequence that contains records of three different magmatic episodes: a Late Neoproterozoic-Early Cambrian event, an Early Ordovician event and a Late Ordovician event. The Late Neoproterozoic-Early Cambrian magmatism gave rise to small bodies of metaluminous and aluminous granitoids and syn-sedimentary bimodal volcanic rocks interbedded in the lower part of an upper Neoproterozoic-lower Cambrian sequence (Cocherie et al., 2005; Castifieiras et al., 2008). This magmatism resulted from the continued oblique subduction of the Iapetus or Prototethys oceans under the northern Gondwana margin from the Late Neoproterozoic to Early Cambrian, (Chen et al., 2000; Crowley et al., 2000; Fernandez-Suarez et al., 2000; Finger et al., 2000; Murphy et al., 2008a,b). The Early-mid-Ordovician magmatic episode emplaced large bodies of aluminous granites into the middle part of the pre-Silurian sequence (Deloule et al., 2002; Cocherie et al., 2005; Castifieiras et al., 2008), although minor sub-volcanic bodies have also been identified (Uesa et al., in press). This event, which is related to the break-up of the northern Gondwana margin, is widely exposed in other areas of the European Variscides, where it has been the subject of numerous studies (e.g. Pin and Marini, 1993; von Raumer, 1998; Helbing and Tiepolo, 2005; Casini and Oggiano, 2008; Melleton et al., 2010). * Corresponding author. Fax: + 34 91 522 2535. E-mail address: navidad@geo.ucm.es (M. Navidad). · . 0,-0 �� • • •..........•• � "" ....... - • .L'� ...... -� , �-�/ 500 km \ .z::::::.�'Fig.2 �� Sardinia Fig. 1. Variscan massifs in southwestern and central Europe with the location of Figs. 2 and 3 (stars). CZ, cantabrian Zone; WAI2, West Asturian Leonese Zone; GTOMZ, GaJicia-Tras­ os-Montes Zone; ell, Central Iberian Zone; Ol\t1Z, Ossa Morena Zone; SPZ, South Portuguese Zone; PAl, Pyrenean Axial Zone; eCR, catalan Coastal Ranges; NPF, Northern Pyrenean Fault; NPM, North Pyrenean Massifs; BAGLG, Beja-Acebuches-GaJicia-Lizard-Giesen oceanic suture Zone. After Robardet (2002). The Late Ordovician magmatic event is represented mainly by calc­ alkaline ignimbrites, andesites, volcaniclastic rocks, diorites and various types of granitic bodies. This magmatism is well known in the Catalan Coastal Ranges (Navidad and Bamolas, 1991; Bamolas and Garda-Sansegundo, 1992) and the Pyrenees (Robert and Thiebaut, 1974; Mart( et al., 1986; Calvet et al., 1988), where it is coeval with normal fault development. The event postdates a mid-(?) Ordovician phase of folding that accounts for an Upper Ordovician unconformity (Casas and Femandez, 2007; Casas, in press). Recent data (SHRIMP U­ Pb in zircon) provide Late Ordovician ages for the protolith of several orthogneisses in the Pyrenees (Casas et al., 201O-this issue). However, the voleanic rocks have not yet been dated, and detailed geochemical data for the volcanic and plutonic rocks are scarce. In this paper, we present new petrographic, geochemical, geochronological and isotopic data for Late Ordovician magmatic rocks of the Pyrenees (Canigo massif) and the Catalan Coastal Ranges (Les Gavarres massif). We discuss the significance of this magmatic event in relation to the Upper Ordovician magmatism of surrounding areas in order to elucidate the evolution of this portion of the peri­ Gondwana margin and the Rheic Ocean during the Late Ordovician. 2. Geological framework 2.1. Canigo massif The Canigo massif exposes one of the most complete pre-Silurian sequences of the Pyrenees (Fig. 2a). A thick (up to 3000 m) aluminous metagranitic body, the Canigo gneiss (Guitard, 1970), derived from an Early Ordovician intrusion (Deloule et al., 2002; Cocherie et al., 2005) divides the sequence into two parts (Fig. 2b and c). The lower series, underlying the gneiss, is up to 1500 m thick and comprises a succession of mica schists with marble, quartzite and metabasite intercalations (Balaig series; Guitard, 1953). The basal part of the sequence, which records high-grade metamorphic conditions (sillimanite-K-feldspar zone), is exposed in the central part of the massif because of a set of Neogene normal faults that crosscut the massif (Fig. 2b). The Casemi gneiss (Guitard, 1970) constitutes the most significant metaigneous intercalation in the Balaig series (Fig. 2c). It is a tabular body up to 1000 m thick mainly made up of fine-grained biotitic and amphibolic granitic gneisses. Previous geochronological data provide an age of 425 ± 18 Ma (Pb evaporation method on single zircon, Delaperriere and Soliva, 1992), but to date no other Silurian ages for magmatism have been found in the pre-Variscan rocks of the Pyrenees. Recent geochronological data indicates a Late Ordovician age for the protolith of this orthogneiss (445.9 ± 4.8 and 451.6 ± 4.8 Ma Ma, SHRIMP U-Pb in zircon, Casas et al., 2010-this issue). The Marialles metadiorite, which comprises metre-scale thick bodies interlayered in the Balaig mica schists, also yields a Late Ordovician age for the formation of its protolith (453 ± 4.4 Ma, SHRIMP U-Pb in zircon, Casas et al, 201O-this issue). The structurally lowest unit recognized in the Canigo massif is the Cadi gneiss (Guitard, 1970), which crops out beneath the Balaig series (Fig. 2c). It is an aluminous metagranitic body with lithological characteristics similar to those of the Canigo gneiss, although recent geochronological data provide a Late Ordovician age for its protolith (456.1 ± 4.8 Ma, SHRIMP U-Pb in zircon, Casas et al., 2010- this issue). The upper sequence, overlying the Canigo gneiss, is up to 3000 m thick and is mainly composed of metapelites and metagreywackes interbedded with numerous layers of marbles, quartzites and calc­ silicates at their base (Canavelles series of Cavet, 1957). It includes syn-sedimentary bimodal metavolcanic rocks (Ayora and Casas, 1986; Navidad and Carreras, 2002), which have yielded Late Neoproter­ ozoic-Early Cambrian ages (SHRIMP U-Pb in zircon; 581 ± 10 Ma, Cocherie et al., 2005; �540 Ma, Castifieiras et al., 2008). In its upper part, this sequence consists of a monotonous succession of shales and sandstones that is of unknown age due to its unfossiliferous character, but which has been classically assigned to the Cambro-Ordovician (Cavet, 1957). An erosive angular unconformity (Fig. 2c) separates the upper sequence from a well-dated Upper Ordovician succession (Cavet, 1957; Hartevelt, 1970; Santanach 1972; Garda-Sansegundo et al., 2004; Casas and Femandez, 2007) that includes syn-sedimentary cale-alkaline metavoleanic rocks of crustal origin (Cavet, 1957; Hartevelt, 1970; Robert and Thiebaut, 1976, Mufioz et al., 1983; Marti et aL, 1986). � • • • " 1 • '" u c • li D - § @ - D . . . rn � - Quartzl!8$ Sch;'ts with decelcified holM Btack phyllites _ ... Malbles and doIomites Metavok:llni<;s Silts1on85 11mb O-Ms mytmile!l Metapemes and metagrawaekes Metabas�es Ortllogoelsses (b) (a) VARISCAN BASEMENT • Gneisses and Paleozoic melasediments Granitoids D Post-Variscan cover sequences N E J � .:: ,---, late Neoproterozoic and L-..J Cambro-Ordovician metasediments _ Canig6 orthogneiss � Anatectic granites G::J Granodiorites o o I � '" w u z � '" " � 453±4.4 Ma Fig. 2. (a) Geological sketch map of the Variscan basement and cover sequences in the Eastern Pyrenees. (b) Geological map of the Canig6 massif. Modified after Guitard (1970). Santanach (1972). casas (1984). and Ayora and Casas (1986). (c) Synoptic stratigraphic column of the pre-Upper Ordovician rocks of the Canig6 massif. Data modified after Guitard (1970). Santanach (1972). casas (1984). Ayora and Casas (1986) and Casas et al. (2010-this issue). 2.2. The Catalan Coastal Ranges In the Catalan Coastal Ranges, the lower part of the pre-Upper Ordovician sequence is made up of a 2000 m thick heterogeneous succession of metapelites with intercalations of marble, quartzite, felsic gneiss and amphibolite (Fig. 3a). A Cambro-Ordovician age has been assigned to this succession (Bamolas and Garda-Sansegundo, 1992; Julivert and Duran, 1992) and U-Pb dating of the metaigneous rocks has yielded an Early Ordovician age (464 and 475 Ma, respectively; Gebauer, 1993). The upper part of the pre-Upper Ordovician sequence consists of a rhythmic shale series in which a Dobrotivian (early Late Ordovician) fossil assemblage has been described (GuW�rrez Marco et al., 2003). The Upper Ordovician sequence unconformably overlies the pre-Upper Ordovician sequence with a basal conglomeratic unit (Barnolas and Garda-Sansegundo, 1992) that exhibits similar lithological characteristics to those of the Central and Eastern Pyrenees. The conglomerates are overlain by the Els Metges unit, a succession of ignimbrites and epiclastic metavol- canic rocks up to 300 m thick (Navidad and Barnolas, 1991; Barnolas and Garda-Sansegundo, 1992). In the upper part of this unit, greywackes with an important volcanic component are located on top of the volcanic section. In other massifs of the Catalan Coastal Range, the greywackes contain a brachiopod faunal assemblage of Caradoc-Ashgill age (Villas et al., 1987). A sedimentary sequence (Mas Cals unit, Fig. 3b) composed of calcareous siltstones with limestones in its lower part overlies the Els Metges unit (Barnolas and Garda-Sansegundo, 1992). The abundant fossil content of the Mas Cals unit allows it to be assigned an Ashgillian age (Sarmiento et al., 1995). Late Visean to Serpukhovian polyphase Variscan deformation accompanied by high temperature-low pressure metamorphism affects the entire succession (Guitard, 1970; Zwart, 1979). A pervasive fold-related crenulation cleavage is the main Variscan deformational structure in the pre-Silurian rocks (Guitard, 1962; Hartevelt, 1970; Santanach, 1972), whereas south-directed thrust sheets are well developed in the overlying Silurian, Devonian and Carboniferous MONTNE (BAR FM) FITOI"tSLO.TES (ANSOYELl. FM) w.s CAL$SI..ATES(I:'STA/NO RIBBON �OSTONES INCLUDING SPOR.lDIC FINE lJO,VERED U"'ESTONES BEDS (JUJOLSFM) (a) CJ Gran�es !Hill Silurian-()eVQnian E2!ill F�or slates o Mas Cals slates , • Q Eis Metges volcan�es _ i I o Cambro-Ordovician _ Bagur limestones PALS rc-,-:-cc�--" • 8AGl.ll UMESfONES �O ........ �5 �===�'Okm Fig. 3. (a) Geological map of the Les Gavarres massif, and (b) Synoptic stratigraphic column of the pre-Upper Ordovician rocks of the Les Gavarres massif. After Barnolas and Garcia­ Sansegundo (1992). successions (Hartevelt, 1970; Domingo et al., 1988; Casas et al., 1989). Infra-Silurian rocks are affected by thrusts, but these are scarce. The thrusts cut the fold-related crenulation cleavage, but are themselves folded by south-verging cleavage-related folds. Thus, thrust and fold development are broadly synchronous. Syn- to late-orogenic (Mos­ covian-Kasimovian, Romer and Sol er, 1995) granitoids are intruded mainly into the upper levels of the succession, producing local contact metamorphism (Autran et al., 1970). 3. Field relations and petrographic data 3.1. The Casemf and Cadf orthogneisses The Casemi orthogneiss is a deformed tabular body intercalated in the Balaig mica schists. Together with the Cadi orthogneiss and several metadiorite bodies, it represents the structurally lowest metaigneous rock recognized in the Canigo massif (Fig. 2c). Both the Canigo massif and the Casemi gneiss exhibit a dome-like structure because of the interference of two Late Variscan folding episodes and/ or because of a fold linked to an underlying Alpine thrust. The Casemi orthogneiss forms the highest relief of the massif, and its best outcrops are located around Canigo peak (Fig. 2b). To the north, the base of the body is made up of felsic and biotitic orthogneisses that gradually change to amphibolic orthogneiss to the south. The Balaig mica schists are the host rocks and exhibit different metamorphic conditions at their lower and upper contacts (sillimanite-K-feldspar and andalusite zones, respectively). The Casemi biotitic orthogneiss shows a metamorphic paragenesis made up of quartz, microcline, biotite, albite-rich plagioclase and garnet, with muscovite, epidote and titanite as retrograde minerals. The amphibolic orthogneiss is mainly composed of Fe-hornblende and plagioclase (An35_40), whereas biotite and K-feldspar are present in the transitional facies to biotitic orthogneiss. Zircon, apatite and allanite are the main accessory minerals. Both orthogneiss types show an elongated granoblastic texture. The Marialles metadiorite crops out as a thin discontinuous layer interbedded in the Balaig metasediments along the southern contact of the Casemt gneiss. It shows a metamorphic paragenesis composed of plagioclase (oligoclase), Fe-hornblende. The Cadi orthogneiss (Fig. 2c) is located in the structurally lowest and more metamorphic layers of the Balaig series (sillimanite-K­ feldspar). It exhibits a porphyroclastic texture with a pervasive foliation. The mineral paragenesis is made up of quartz, biotite, plagioclase and K-feldspar. 3.2. EIs Metges volcanic sequence In Les Gavarres massif, the felsic volcanic and volcaniclastic rocks overlie the basal conglomerates of the Upper Ordovician sequence (Fig. 3b). Volcanic rocks are rhyodacitic tuffs and massive rhyolites with well preserved relic magmatic structures such as flow lamination in the ash beds, pumice shards, U-, V- and T-shaped glass shards, and perlitic and spherulitic textures. Pseudomorphs after the mafic minerals, embayed textures in quartz, and feldspar, apatite and ilmenite phenocrysts are also abundant in the tuff. The epiclastic ash beds differ significantly from the rhyodacitic tuff with the presence of numerous fragments of prismatic quartz, metamorphic rocks, oxide­ welded quartz and albite spherulites in the matrix. The volcanic rocks consequently provide ample evidence of their phreato-magmatic origin. 4. U-Pb zircon dating Zircon grains from the Els Metges pyroclastic tuff (sample NF-29) were dated by U-Pb 10-TIMS (isotope dilution-thermal ionization mass spectrometry) (Table 1). As noted above, zircons from the Casemi and Cadi orthogneisses and a metadiorite body were recently dated by SHRIMP U-Pb by Casas et al. (2010-this issue). U-Pb analytical work for sample NF-29 was carried out at the Institute of Geology, University of Oslo, Norway. Zircon grains selected for analysis were intensely air-abraded, producing a ca. 40 to 60% reduction in grain size. The abraded zircon grains were washed in 4 N HN03 on a hot plate and rinsed repeatedly with H20 and acetone (with ultrasonication after each rinsing step). After washing and Table 1 Results of ID-TIMS U-Pb dating of sample NF-29. (Les Gavarres massif). Sample, Weight· U Thl Pbcomc 206Pbd 206Pb" analysis [pg[ [ppm] U' [pg[ 204Pb B'U NF-29 Z1 327.1 0.53 1.3 1194.1 0.07357 NF-2922 4 435.0 0.26 5.0 1613.4 0.07286 NF-29 Z3 2 451.6 0.50 2.7 1554.2 0.07321 NF-29 Z4 5 223.9 0.25 1.5 3545.7 0.07316 • Weights are known with ca. 20 to 40% uncertainty. b Model ThjU ratio estimated from 2°OPbfo6Pb ratio and age of the sample. C Total common Ph in sample, including initial and blank Pb. d Measured ratio, corrected for fractionation and spike contribution . 20' 207Pb 20 207Pb 20 Apparent age [Ma] labs] 235U labs] 206Pb labs] 206Pb 207Pb 207Pb B'U mU 206Pb 0.00033 0.5699 0.0043 0.05618 0.00034 457.6 457.9 459.4 0.00036 0.5640 0.0035 0.05614 0.00022 453.4 454.1 457.8 0.00043 0.5629 0.0040 0.05577 0.00029 455.5 453.4 443.2 0.00029 0.5647 0.0027 0.05598 0.00014 455.2 454.6 451.5 • Corrected for spike, fractionation, blank and initial common Pb (Stacey and Kramers, 1975). r 20 uncertainty calculated by error propagation procedure that takes into account internal measurement statistics and external reproducibility as well as uncertainties in blank and common Pb correction. drying, the zircon grains were individually weighed in a precision balance. After the size reduction induced by air-abrasion, the analyzed grains weighed between 1 and 5 j.I.g (Table 1). A mixed 205PbP35U spike was added to the sample after weighing and transferring to the dissolution vessel. Zircon was dissolved in HF (+ HN03) in Teflon mini-bombs at approximately 185 cC for 5 days. Since all zircons weighed less than 5 j.I.g, no chemical separation of U­ Pb was performed and the whole sample taken in 6 N HCl was dried and loaded on out-gassed Rhenium filaments with H3P04 and silica geL Isotopic ratios were measured on a Finnigan-MAT 262 mass spectrometer by peak jumping on a secondary electron multiplier (ion counting mode). Total procedural blanks were less than 2-5 pg Pb and 0.1-0.3 pg U. The Stacey and Kramers (1975) model was used to subtract initial common Pb in excess of the laboratory blank. Decay constants are those of jaffey et al. (1971). Concordia plots and concordia age calculations were created using Isoplot 3.00 (Ludwig, 2003). Four single grain analyses were performed on zircons from sample NF-29. All zircons were idiomorphic stubby equant prisms with a pink hue. The results of the U-Pb analyses are given in Table 1 and are shown graphically in a concordia plot on Fig. 4. All four analyses are concordant and reproducible within error of individual measure- ments. The concordia age (Ludwig, 1998) is the most statistically robust age for this data set and yields a date of 455.2 ± 1.8 Ma (decay constant uncertainties included) as the best estimate for the crystallization of zircon. This age coincides within error to the base of the sixth stage of the Ordovician (455.8 ± 1.6 Ma), which corresponds to the middle Upper Ordovician (Katian), which is roughly equivalent to the middle Caradoc in the regional British series division (see Finney, 2005). 5. Whole-rock geochemistry and Sr-Nd isotopes Representative samples of Lite Ordovician orthogneisses, meta­ diorite, tuffs and epiclastic ashes from the Canigo and Les Gavarres massifs were selected for whole-rock and Sr-Nd isotopic analyses (Tables 2 and 3). 5.1. Analytical techniques Major and trace elements were measured in whole-rock samples using ICP-OES (Inductively Coupled Plasma-Optical Emission Spec­ trometry) for major and minor elements and ICP-MS (Inductively Coupled Plasma-Mass Spectroscopy) for trace elements at the Spectro­ chemical Laboratory of the Centre de Recherches en Petrographie et 0.0742 .,---------------------, ::> m M N .c 0.. l!i N NF-29 0.0738 Eis Metges 0.0734 0.0730 0.0726 0.0722 0.554 0.558 0.562 Concordia Age: 455.2 ±1.B Ma 95% confidence, decay-const. errs included MSWD (of concordance) = 0.84 Probability (of concordance) = 0.36 0.566 0.570 0.574 0.578 Fig. 4. U-Pb Concordia diagram for sample NF-29 (Les Gavarres massif). Table 2 Chemical analyses of samples from the canig6 massif. Chemical analyses of Eis Metges volcanic tuffs (VI) and epic1astic ashes (EA) from Catalonian coastal ranges (Les Gavarres massif). 100762 100766 100786 100768 C96940F CG-07-02' Felsic orthogneiss Bt orthogneiss CG-07-03' Bt orthogneiss C96767Q2b Et orthogneiss Bt orthogneiss CG-07-05' Amp orthogneiss Amp orthogneiss Amp orthogneiss 770288 Amp orthogneiss Marialles diorite Si02 Ti02 AhO} Fe20} MnO MgO caD Na20 K,O P20S H,O Total B, Be Co Q Cu G, Hr Nb Ni Rb Se St T, Th U V Y Zn Zt '" re Pt Nd Sm Eu Gd Th Ily Hu Et Tm Yb In 74.55 0.2 13.03 2.39
    <3 DAmp orthogneiss 0 a: .•••• Epiclastic ashes a: ti<:I Marialles diorite 0,1 0.1 I! Cadf orthogneiss Rb Th K Nb C. Nd HI Srn Y Rb Th K Nb C. Nd HI Srn Y Ba U Ta La Sr P Zr T; Ba U Ta La Sr P Zr T; 1000 1000 (C) les Gavarres massif (d) Canigo massif 2 2 � {5 100 {5 100 c "":"";;;;;;;;;;;;;""";;::�::;;;;;::::::::::�::; c 0 0 .c .c <.> � 32 <.> <.> 10 0 10 0 a: a: m ID ;;: " 'E E ::I '0 • �. ." E .0 0 m ID ;;: " 'E E ::I'OD>'·O ·"'" E.o 0 -' () Z (l)w�� o 0 UJ f- >- -' -' () Z U)UJ(9�OIW f->- -' "- J: CL Fig. 6. (a) and (b) Mantle-normalized trace element plots (normalization values after Wood et ai., 1979) from the Les Gavarres and Canig6 massifs. (c) and (d) Chondrite-normalized REE plots (normalization values after Evensen et aI., 1978) for the Les Gavarres and Canig6 and massifs. tuff (LaN/LuN = 5.9-7.6) than in the epiclastic ash (LaN/LuN = 8.5- 13.7) (Fig. 6c), and a negative Eu anomaly (Eu/Eu' = 0.44-0.52, in the volcanic tuff; Eu/Eu* = 0.20-0.30, in the epiclastic ash). The negative Eu anomaly indicates crystallization from evolved felsic magmas with plagioclase fractionation. In the Canigo samples, the chondrite­ normalized REE plot (Fig. 6d) shows a moderate slope (LaN/ LUN=4-7) with a negative Eu anomaly (Eu/Eu*=0.2s-0.76). The pattern for the Marialles metadiorite displays a similar slope (LaN/ LUN=4.23), although its REE content is lower and the Eu negative anomaly is less pronounced (Eu/Eu* = 0.93). All these REE patterns are comparable to those found in calc-alkaline magmas from active margin environments or from evolved crustal melts. To gain further insight into the possible tectonic setting for these rocks, the relation between characteristically immobile elements Nb and Y (Pearce et al., 1984) has been plotted (Fig. 7). In this diagram, all E Q. S .0 Z 1000 -.--�----------, Canigo massif • Bt orthogneiss I 100 10 VAG+ syn·COLG • • • • WPG / / ORG a / Casemf DAmp orthogneiss � Cad! orthogneiss Les Gavarres massif • Volcanic tuff o Epiclastic ashes 10 100 1000 Y (ppm) Fig. 7. Nb versus Y tectonic setting discrimination diagram for felsic rocks (Pearce et al.. 1984). samples fall into the area of anomalous rift, except for the volcaniclastic rocks, which are displaced into the arc setting. However, as Twist and Harmer (1987) point out, the trace element geochem­ istry in juvenile rocks contaminated by crustal rocks may preserve the tectonic setting signature of the crustal component. 5.3. Sr-Nd isotopes Thirteen samples were selected for Sr and Nd isotopic analysis, ten from the Canigo massif and three from the Les Gavarres massif (Table 3). Initial 87Srl6Sr values in some of the samples are very low «0.702), suggesting that this isotopic system was disturbed, ruling out its use. Some samples of volcanic tuff and felsic orthogneiss have plausible initial isotopic ratios (between 0.7023 and 0.7095), but these results were rejected owing to the sensitivity of the Rb-Sr system to disturbance. In contrast, the initial 143Nd;t44Nd ratios are more uniform. In the Casemi orthogneiss, this ratio varies between 0.5118 and 0.5119, and is 0.5120 in the Marialles metadiorite. The Cadi orthogneiss shows a Nd isotope ratio (0.5118) similar to the most differentiated Casemi samples (Table 3). The Els Metges volcanic tuff has 143Nd;t44Nd ratios between 0.5117 and 0.5118. The fractionation of Sm/Nd (fSm/Nd, DePaolo, 1988) from the chondritic uniform reservoir (CHUR) for crustal rocks shows a narrow range of values when averaged over volumes of typical igneous source areas (Haskin and Paster, 1979; Ben Othman et al., 1984). This fractionation index can be plotted versus c;Nd(T) to discriminate the tectonic setting (McLennan and Hemming, 1992). On such a plot (Fig. 8), all the studied samples cluster in a small area, although the Els Metges volcanic tuff and the Cadi orthogneisses are slightly more fractionated from CHUR than the Casemi orthogneiss. The Pyrenean samples have also been plotted on Fig. 8, together with an average for Early Ordovician igneous protoliths from the Central Iberian zone of the Iberian massif, ten orthogneisses from Guadarrama and five Ollo de Sapo metavolcanic rocks (data from Castro et al., 2003; Villaseca 0.2 1-------1---;::=::::� MORS 0.1 0.0 +------------+---''-'''-1 -0.1 " '" E -0.2 o fQ -0.3 +---------y/ -0.4 Precambrian Upper crust older than 1.6 Ga -0.5 Passive margin­ craton Active margin �--' -0.6 +---,--,----.--t--,--+--r----i -30 -20 -10 ENd (455) o Canigo massif • 8\ orthogneiss I � D Amp orthogneiss � 181 Marialles diorile � Cadf orthogneiss Les Gavarres massif • Eis Melges volcanic tuff 01her samples o 0110 de Sapo (eIZ) • Guadarrama (eIZ) 1 0 Fig. 8.fSm/Nd verslls£Nd(455) diagram for the analyzed samples showing fields used to discriminate tectonic setting (McLennan and Hemming, 1992). et al., 1998; and unpublished data). All data have been corrected assuming an age of 455 Ma. The uniformity in the jSm/Nd and c;Nd values suggests a common or similar composition of the crustal source for all of the samples. Furthermore, samples of crustal provenance (such as the Els Metges tuff, and the Cadi and Guadarrama orthogneisses) plot within the passive margin area, whereas samples with a juvenile signature (such as the Casemi orthogneiss and the MariaUes metadiorite) plot along the boundary between the passive and active margin areas, which is also where the OUo de Sapo volcanic tuff plots. These results suggest a transitional tectonic setting between an active margin (Casemi orthogneiss) and a passive margin (Els Metges tuff). Epsilon neodymium (c;Nd) values are negative in all of the samples. In the Canigo massif, the most depleted value (- 0.6) is displayed by the MariaUes metadiorite, whereas values range between -1.2 and -4.0 in the orthogneiss samples. The Els Metges volcanic tuff has the most enriched c;Nd values, between -4.8 and -5.2. In an age versus c;Nd diagram (Fig. 9), most of the Nd evolution lines intersect the depleted mantle curve at ages (TOM) ranging from 1.3 to 1.7 Ga, although two samples from the Canigo massif yield TOM ages of 2.0 and 2.2 Ga. The significance of these Mesoproterozoic TOM ages is difficult to evaluate in igneous felsic rocks without additional supporting data (e.g., zircon inheritance) and will be discussed in the following section. 6. Discussion Whole-rock geochemistry shows that the orthogneiss samples from the Canigo massif are compositionaUy uniform. As fractionation processes cannot be recognized, crustal contamination of mantle melts is the most probable process accounting for the formation of the various analyzed lithologies (biotite and amphibole orthogneisses, and metadiorite). By contrast, in the Les Gavarres massif, the negative trends of the samples suggest that the epiclastic ash resulted from the fractionation of the rhyolitic tuff. However, given the marked depletion of mobile elements (e.g., NazO and Rb), post-depositional 8 6 4 2 " Z 0 w CHUR -2 -4 -6 -8 -10 500 1 000 1500 2000 Age Fig. 9. Age versus eNd diagram for the canig6 and Les Gavarres massifs. Nd evolution lines are shown in grey. processes could also have modified the composition of the epiclastic ash. Nd isotopes suggest that the rocks originated in an active margin. In this setting, the composition of mantle-derived melts is strongly influenced by the crust, giving rise to negative c;Nd values (from - 1 to -4.8) - high negative c;Nd values ( - 3 to - 4) suggesting considerable crustal contamination. The isotopic signature of the Casemi gneiss (c;Nd between -1.4 to - 4) and the Marialles metadiorite (c;Nd = - 0.6) indicates that their protoliths were derived from mantle melts with no uniform crustal contamination. In contrast, the higher c;Nd negative values (- 4.2 to - 5.2) of the Cadi orthogneiss and Els Metges volcanic rocks indicate a crustal origin. Nd model ages (TOM) for most of the rocks are similar to those in samples from the Iberian massif (Femandez-Suarez et al., 2000; Castro et al., 2003; Rodrtguez Alonso et al., 2004; GuW�rrez Alonso et al., 2005; Murphy et al., 2008a). However, the strong crustal influence on the c;Nd values for felsic rocks rules out an unequivocal interpretation of the TOM ages. Recent geochronological studies in the Pyrenean massifs have reported Pan-African (800 Ma), Mesopro­ terozoic (1.0 Ga), Ebumean (1.8 to 2.0 Ga) and Archean (2.5 to 3.1 Ga) ages in inherited zircons from Neoproterozoic-Early Cambrian and Early Ordovician igneous rocks (Cocherie et al., 2005; Castifieiras et al., 2008). The absence of 1.4-1.7 Ga inherited zircons suggests that the TOM ages obtained in this study are the result of a mixture of an old component (probably Ebumean) and juvenile mantle melts (Ordovi­ cian). This favours a Neoproterozoic source mixed with Paleoprother­ ozoic components as the probable source for the melts that resulted in the Els Metges volcanic rocks, the Cadi aluminous orthogneiss and the source of contamination in the Casemi metaluminous orthogneiss. Moreover, the similarity of the isotopic signatures of some of the Early Ordovician aluminous magmatic protoliths of the Iberian massif (such as the OUo de Sapo volcanic Formation and the Guadarrama orthogneisses) suggests the repeated extraction of crustal melts from a common source during the Ordovician as previously proposed (e.g., Femandez-Suarez et al., 2000). This crustal recycling would account for the volcanic arc signature of the samples. This signature was probably inherited by melting of pre-existing Neoproterozoic­ Early Paleozoic calk-alkaline crust. Crustal recycling has been invoked to explain the voleanic arc affinity of the Early Ordovician Ollo de Sapo magmatic rocks (Dlez Montes et al., 201O-this issue). The age obtained for the Les Gavarres volcanic rocks (455.2 ± 1.8 Ma) is equivalent to the age of crystallization of the magmatic protoliths of the Canigo gneisses (Cadi and Casemi gneiss, and Marialles metadiorite 456- 446 Ma, Casas et al., 201 O-this issue). Thus, the Les Gavarres volcanic rod{s can be regarded as the volcanic equivalent of coeval Late Ordovician plutonic rocks. In the Pyrenees, this Late Ordovician magmatic event can be distinguished from the Early Ordovician event (477 and 467 Ma) that gave rise to the protoliths of large gneissic bodies such as Aston­ Hospitalet, Canigo, Roc de Frausa and Albera (Castifieiras et al., 2008; Oenele et al., 2009; Liesa et al., in press). Hence, two different Ordovician magmatic episodes can be identified: an Early Ordovician magmatism that does not seem to be related to any deformational or metamorphic episode, and a Late Ordovician episode, which is coeval with normal fault development. These two magmatic events are separated by a compressional episode, probably of mid-Ordovician age (Casas, in press), that is at odds with a continuous extensional regime related to the opening of the Rheic or the Rheic and Paleotethys oceans during the Ordovician and Silurian (e.g. Nance et al., 2010-this issue). This evolution differs from that observed in the rest of the Iberian Massif, where thick (up to 4500 m) Early Ordovician detrital sediments were deposited coevally with volca­ nism and plutonism (Perez-Estaun et al., 1990; Valverde-Vaquero et al., 2005; Dlez Montes et al., 2010-this issue), Late Ordovician magmatic activity is scarce (Valverde-Vaquero et al., 2007), and evidence of Ordovician deformation is limited (Martfnez Catalan et al., 1992). By contrast, Late Ordovician ages have also been obtained for magmatic rocks in several massifs of the Mediterranean realm: French Massif Central (Pont de Lam orthogneiss: 456 ± 3 Ma and Gorges d'Heric orthogneiss: 450 ± 6 Ma, Roger et al., 2004), Sardinia (Tanau­ nella orthogneiss: 458 ± 7 Ma and Lode orthogneiss: 456± 14 Ma, Helbing and Tiepolo, 2005) and Sicily (felsic porphyroids, 456-452 Ma, Trombetta et al., 2004). Late Ordovician (450-430 Ma) granitoids with an arc signature are also distributed in a longitudinal zone across all the Alpine realms (Zurbriggen et al., 1997; von Raumer, 1998) and in the Southern Alps. Bellieni and Sassi (1981) proposed an "Upper Ordovician rhyolitic plateau" formed by a granite-rhyolite association related to an anorogenic event favoured by extension. Late Ordovician (450-420 Ma, Rb-Sr isochron) crustal anatectic magmatism comprising tonalite to granite bodies of calc-alkaline composition and cogenetic acid tuffs have been described in the Eastern Alps, (Sassi and Zirpoli, 1979; Heinisch, 1981), and a 450-460 Ma magmatic event has also been reported in Penninic (Guillot et al., 2002) and Central Swiss Alps (Schaltegger et al., 2003). In addition, the occurrence of a mid-Ordovician tectono-metamor­ phic event also appears to be characteristic of the eastern part of the northern Gondwana margin. In the Alps (Sesia Zone), a mid-Ordovician compressional event associated with amphibolite fades metamorphism has been reported by Zurbriggen et al. (1997) and Handy et al. (1999). High-pressure metamorphism in the Aar Massif (468 Ma, Schaltegger et al.,2003) and in the Southern Alps (457 ± 5 Ma, Franzand Romer,2007) also provides evidence of an Ordovician orogenic cycle subsequently overprinted by Variscan deformation. This suggests that the Pyrenees and the Catalan Coastal Ranges were probably located in a different position on the northern Gondwana margin from that occupied by the rest of the Iberian Massif, and that both areas evolved differently following the Early Ordovician birth of the Rheic Ocean. In the Pyrenees and the Catalan Coastal Ranges, Late Ordovician magmatic activity was likely triggered by extensional tectonics affecting a previous thickened crust This evolution supports the proposal of Stampfli et aL (2002) and von Raumer et al. (2002), who suggested that the Late Ordovician extensional regime resulted from the collapse of a short-lived cordillera that formed during the mid-Ordovician. This cordillera, formed during a transient mid-Ordovician orogenic pulse started to collapse in a setting dominated by the Gondwana-directed subduction of a former peri­ Gondwanan ocean. 7. Conclusions Late Ordovician magmatism (ca. 455 Ma) in the northeastern Iberian Variscan massifs is represented by metaluminous and aluminous metagranites and metadiorites, and by sequences of cale­ alkaline metavolcanic rocks. The metaluminous metagranites were derived from juvenile mantle melts contaminated by crustal melts, whereas the aluminous metagranites and cale-alkaline metavolcanic rocks originated from crustal melting. This Late Ordovician magmatic event differs from the Early Ordovician event that gave rise to the large gneissic bodies (Aston-Hospitalet, Canigo, Roc de Frausa and Albera gneisses) that form the core of the Eastern Pyrenean metamorphic massifs. The emplacement of the volcanic and plutonic bodies was coeval with the development of extensional tectonics. A recycled Neoproterozoic-Early Paleozoic crust with a Paleopro­ terozoic signature and calk-alkaline affinity constitutes the most probable source of the crustal melts. Successive extraction of crustal melts from equivalent sources during the Ordovician would account for the geochemical similarities between the Early Ordovician and Late Ordovician magmatic rocks of the Iberian massif, the Pyrenees and the Catalan Coastal Range. The evolution suggests that the Pyrenees and the Catalan Coastal Ranges were located in different positions on the northern Gondwana margin than those occupied by the Iberian Massif, and that both areas evolved differently following the Early Ordovician birth of the Rheic Ocean. Acknowledgements This work was supported by the Spanish Commission for Science and Technology, project CGL-2007-66857C02-02.JMC acknowledges the support of the project Consolider-Ingenio 2010 programme, under CS02006-00041 "Topoiberia". Some of the analyses were funded by the Petrology Department of the Paris VI and Complutense (Madrid) universities. MN acknowledges Professor G. Guitard for his introduc­ tion to the study of the Canigo massif and his help in the transfer of samples, and the Paris VI University for analytical support. The detailed comments of Oamian Nance, Gabriel Gutierrez-Alonso and two anonymous referees greatly enhanced the original version of the manuscript. George von Knorring improved the English language. Appendix A. Supplementary data Figure 51. Figure 52. Note: The supplementary material accompanying this article is available at doi:1 O.1016jj.gr.2009.11.013. References Autran. A.. Fonteilles. M. Guitard. G .. 1970. Relations entre les intrusions de granitoides, l'anatexie, et le metamorphisme regional considerees principalment du point de vue de l'eau: cas de la Chaine hercynienne des Pyrenees Orientales. Bulletin Societe Geologique de France 7, 673-731. Ayora, C, casas, j.M., 1986. Stratabound As-Au mineralization in pre-Caradocian rocks form the Vall de Ribes, Eastern Pyrenees, Spain. Mineralium Deposita 21, 278-287. Barnolas, A., Garcia-Sansegundo, j., 1992. Caracterizaci6n estratigrafica y estructural del Paleozoico de Les Gavarres (Cadenas Costero Catalanas, NE de Espana). Boletin Geol6gico y Minero 103, 94-108. Bellieni, G., Sassi, F.P., 1981. In: Karamata, S., Sassi, F.P. (Eds.), New chemical data and a review on the South-Alpine "Pre-Hercynian Rhyolitic Plateau" in the Alps: IGCP no 5, Newsletter, vo!. 3, pp. 22-27. Ben Othman, D., Polve, M, Allegre, Cj.. 1984. Nd-Srisotopic composition of granulites and constraints on the evolution of the lower continental crust. Nature 307, 510-515. Brown, CG" Thorpe, R.S" Webb, p.c, 1984, The geochemical characteristics of granitoids in contrasting arcs and comments on magma sources, journal of the Geological Society of London 141, 413-426, Calvet, p" Lapierre, H., Charvet, j" 1988. Diversite du volcanisme Ordovicien dans la region de Pierrefitte (Hautes-Pyrenees): rhyolites calco-alcalines et basaltes alcalins, Comptes Rendus de I'Academie des Sciences de Paris 307, 805-812, Carignan, J" Hild, p" Mevelle, G" Morel, j" Yeghicheyan, D" 2001, Routine analyses of trace elements in geological samples using flow injection and low pressure on-line liquid chromatography coupled to ICP-MS: a study of geochemical reference materials BR. DR-N, UB-N, AN-G and GH. Geostandards Newsletter 25, 187-198, Casas, J,M" 1984, Estudi de la deformacio ens els gneiss del massis del Canigo, Ph,D, thesis, Barcelona University, Casas, J,M" in press, Ordovician deformations in the Pyrenees: new insights into the significance ofPre-Variscan ("Sardic") tectonics, Geological Magazine, doi: 10, 1017/ S0016756809990756, Casas, J,M, Fern.lndez, 0" 2007, On the Upper Ordovician unconformity in the Pyrenees: new evidence from the La Cerdanya area, Geologica Acta 5, 193-198. Casas, J,M" Domingo, E, Poblet, J" Soler, A, 1989, On the role of Hercynian and Alpine Thrusts in the Upper Paleozoic rocks of the Central and Eastern Pyrenees, Geodinamica Acta 3, 135-147, Casas,J,M, castifieiras, p" Navidad, M" Liesa, M" carreras,j" 2010, New insights into the Late Ordovician magmatism in the Eastern Pyrenees: U-Pb SHRIMP zircon data from the canigo massif. Gondwana Research 17 (2-3),317-324 (this issue), Casini, L, Oggiano, G" 2008, Late orogenic collapse and thermal doming of the northern Gondwana margin incorporated in the Variscan Chain: a case study from the Ozieri Metamorphic Complex, northern Sardinia, Italy, Gondwana Research 13, 396-406, Castifieiras, p" Navidad, M, liesa, M" Carreras, j" casas, J,M, 2008. U-Pb zircon ages (SHRIMP) for cadomian and Lower Ordovician magmatism in the Eastern Pyrenees: new insights in the pre-Variscan evolution of the northern Gondwana margin, Tectonophysics 461, 228-239, Castro, A, Corretge, LG" de la Rosa, j" Fernandez, C, Lopez, S,' Garcia Moreno, 0" Chacon, H., 2003, The appinite-migmatite complex of Sanabria, NW Iberian Massif, Spain, journal of Petrology 44, 1309-1344, Cavet, p" 1957, Le Paleozoique de la zone axiale des Pyrenees orientales fram;aises entre le Rousillon et I'Andorra (etude stratigraphique et paleontologique), Bulletin Service carte Geologique de France 55 (254), 303-518. Chen, E, Hegner, E, Todt, W" 2000, Zircon ages and Nd isotopic and chemical compositions of orthogneisses from the Black Forest, Germany, Evidence for a Cambrian magmatic arc. International journal of Earth Sciences 88, 791-802, Cocherie, A, Baudin, Th" Autran, A, Guerra, C, Fanning, CM, Laumonier, B" 2005, U-Pb zircon (ID-TIMS and SHRIMP) evidence for the Early Ordovician intrusion of metagranites in the late Proterozoic canaveilles Group of the Pyrenees and the Montagne Noire (France), Bulletin Societe Geologique de France 176, 269-282, Crowley, Q,G" Aoyd, PA, Winchester, jA, Franke, W" Holland, J,G" 2000, Early Palaeozoic rift-related magmatism in Variscan Europe: fragmentation of the Armorican Terrane Assemblage, Terra Nova 12, 171-180, Delaperriere, E" Soliva,J" 1992, Determination d'un age Ordovicien superieur - Silurien pour les gneiss de casemi (Massif du canigou, Pyrenees orientales) par la methode d'evaporation de plomo sur monozircon, Comptes Rendus de I'Academie des Sciences de Paris 314, 345-350, Deloule, E, Alexandrov, p" Cheilletz, A, Laumonier, B., Barbey, p" 2002, In-situ U-Pb zircon ages for Early-Ordovician magmatism in the eastern Pyrenees, France: the Canigou orthogneisses, International journal of Earth Sciences 91, 398-405, Denele, Y., Barbey, p" Deloule, E, Pelleter, E, Olivier, Ph" Gieizes, G" 2009, Middle Ordovidan U-Pb age of the Aston and Hospitalet orthogneissic laccoliths: their role in the Variscan evolution of the Pyrenees, Bulletin Societe Geologique de France 180, 209-221. DePaolo, DJ" 1988, Neodymium Isotope Geochemistry, Springer-Verlag, Berlin, Diez Montes, A, Martinez catalan, J,R., Bellido Mulas, E, 2010. Role of the 0110 de Sapo massive felsic volcanism of NW Iberia in the Early Ordovician dynamics of northern Gondwana, Gondwana Research 17 (2-3), 363-376 (this issue), Domingo, E, Mufioz, J,A, Santanach, p" 1988. Estructures d'encavalcament en les materials del socol hercinia del massis de la Tossa d'Alp, Pirineu oriental. Acta Geologica Hispanica 23, 141-153, Evensen, N,M" Hamilton, P,J" O'Nions, RJ" 1978. Rare earth abundances in chondritic meteorites, Geochimica Cosmochimica Acta 42, 1199-1212, Fernandez-Suarez,j" Gutierrez Alonso, G" jenner, G" Tubret, M, 2000, New ideas on the Protherozoic-Early Paleozoic evolution of NW Iberia: insights from U-Pb detrital zircon ages, Precambrian Research 102, 185-206, Finger, E, Hanzl, p" Pin, C, von Quadt, A, Steyrer, H.P" 2000, The Brunovistulian: Avalonian Precambrian sequence at the eastern end of the Central European Variscides? In: Franke, W" Haak, v., Oncken, 0" Tanner, D, (Eds,), Orogenic Processes, Quantification and Modelling in the Variscan Belt: Geological Society of London, Special Publication, voL 179, pp, 103-112, Finney, S,' 2005, Global series and stages for the Ordovician System. A progress report, Geologica Acta 3, 309-316, Franz, L, Romer, R.L, 2007, caledonian high-pressure metamorphism in the Strona­ Ceneri-Zone (Southern Alps of southern Switzerland and northern Italy), Swiss journal of Geosciences 100, 457-467, Garcia-Sansegundo, j" Gavalda, j" Alonso, J,L, 2004, Preuves de la discordante de 1'0rdovicien superieur dans la zone axiale des Pyrenees: exemple de dome de la Garona (Espagne, France), Comptes Rendus Geosciences 336, 1035-1040, Gebauer, D" 1993, Intra-grain zircon dating within the Iberian massif: 0110 de sapo Augengneiss, bimodal gneisses from the Massif de Guilleries (Girona), Graywacke of the Tentudia Group (Serie Negra, Sw Spain) and the HP/HT rock association at Cabo Ortegal (Galicia), XII Reuniao de Geologia do Oeste Peninsular Comunica\oes, Evora 2, 41-46, Guillot, E, Schaltegger, U" Bertrand, J,M" Deloule, E, Baudin, T., 2002, Zircon U-Pb geochronology of Ordovician magmatism in the polycyclic Ruitor Massif (Internal W-Alps), International journal of Earth Sciences 91, 964-978. Guitard, G" 1953, La structure du massif du canigou, Bulletin Societe Geologique de France 3, 907-924, Guitard, G" 1962, Sur l'existence de plis "anteschisteux"dans le Paleozolque inferieur et les gneiss du massif du canigou-Caran\a (Pyrenees Orientales), Bulletin de la Societe Geologique de France, Compte Rendu Sommaire des Seances de la Societe Geologique de France 4, 82, Guitard, G" 1970, Le metamorphisme hercynien mesozonal et les gneiss oeilles du massif du Canigou (Pyrenees-Orientales), Memoire Bureau Recherches Geologi­ ques et Minieres 63, Gutierrez Alonso, G" Fernandez Suarez, J" Collins, AS, Abad, t, Nieto, E, 2005, Amazonian Mesoprotherozoic basement in the core of the Ibero-Armorican Arc: 4°Arf9Ar detrital mica ages complement the zircon's tale, Geology 33, 637-640, Gutierrez Marco, J,C, Robardet, M" Rabano, t, Sarmiento, G,N" San Jose Lancha, MA, Herranz Araujo, p" Pieren Pidal, AP" 2003, Ordovician, In: Gibbons, W" Moreno, T. (Eds,), The Geology of Spain, The Geological Society, London, pp, 31-49, Handy, MR, Franz, L, Helier, E,janott, B" Zurbriggen, R, 1999, Multistage accretion and exhumation of the continental crust (Ivrea crustal section, Italy and Switzerland), Tectonics 18, 1154-1177, Hartevelt, jA, 1970, Geology of the Upper Segre and Valira valleys, Central Pyrenees, Andorra/Spain, Leidse Geologische Mededelingen 45, 167-236, Haskin, LA, Paster, T.P" 1979, Geochemistry and mineralogy of the rare earths, In: Gschneidner, KA (Ed,), Handbook on the Physics and Chemistry of the Rare Earths, North-Holland, Amsterdam. Hastie, AR., Kerr, AC, Pearce, jA, Mitchell, S,E, 2007, Oassification of altered volcanic island arc rocks using immobile trace elements: development of the Th-Co discrimination diagram. journal of Petrology 48, 2341-2357, Heinisch, H" 1981, Preliminary report on Early Paleozoic acidic volcanism in the Eastern and Southern Alps - a review. In: Karamata, S,' Sassi, EP, (Eds,), IGCP no 5: Newsletter, voL 3, pp, 80-88. Helbing, H., Tiepolo, M, 2005, Age determination of Ordovician magmatism in NE Sardinia and its bearing on Variscan basement evolution, journal Geological Society Inndon 162, 689-700, jacobsen, S,B" Wasseburg, G,J" 1980, Sm-Nd isotopic evolution of chondrites, Earth and Planetary Science Letters 50, 139-155, jaffey, AH., Flynn, KE, Glendenin, LE" Bentley, w.c, Essling, AM, 1971, Precision measurement of half-lives and specific activities of235U and 23BU, Physical Review C 4 (5), 1889-1906, julivert, M, Duran, H., 1992, El Paleozoico inferior de las Cadenas Costero Catalanas, In: Gutierrez Marco, J,C, Saavedra, j" Rabano, t (Eds,), Paleozoico Inferior de Ibero­ America, Universidad de Extremadura, pp, 607-613, Liesa, M, carreras, j" castifieiras, p" Casas, J,M, Navidad, M, Vila, M, in press, U-Pb zircon age of Ordovician magmatism in the Albera Massif (Eastern Pyrenees), Geologica Acta, Ludwig, KR., 1998. On the treatment of concordant uranium-lead ages, Geochimica et Cosmochimica Acta 62, 665-676, Ludwig, KR., 2003, User's manual for Isoplot 3,00, a geochronological toolkit for Microsoft Excel. Berkeley Geochronology Center Special Publication 4, 74 pp, Marti,j" Mufioz, jA, Vaquer, R., 1986, Les roches volcaniques de I'Ordovicien superieur de la region de Ribes de Freser-Rocabruna (Pyrenees catalanes): caracteres et signification, Comptes Rendus de I'Academie des Sciences de Paris 302, 1237-1242, Martinez Catalan, J,R., Rodriguez, M,P.H., Alonso, pv" Perez Estatin, A, Gonzalez Indeiro, E, 1992, Lower Paleozoic extensional tectonics in the limit between the West Asturian-Leonese and Central Iberian Zones of the Variscan Fold-Belt in NW Spain, Geologische Rundschau 81, 545-560, McLennan, S,M" Hemming, S,' 1992, Samarium-neodymium elemental and isotopic systematic in sedimentary rocks, Geochimica Cosmochimica Acta 561, 2015-2050, Melleton, J" Cocherie, A, Faure, M" Rossi, p" 2010, Precambrian protoliths and Early Paleozoic magmatism in the French Massif Central: U-Pb data and the North Gondwana connection in thewest European Variscan belt, Gondwana Research 17, 13-25, Mufioz, jA, Sabat, E, Santanach, p" 1983, Stratigraphic correlation form in Eastern Pyrenees, In: Sassi, EP" Szederkenyi, V. (Eds,), IGCP No 5: Newsletter, voL 5, pp, 102-107, Murphy, BJ" Gutierrez Alonso, G" Fernandez Suarez, j" Braid, j" 2008a, Probing crustal and mantle lithosphere origin through Ordovician volcanic rocks along the Iberian passive margin of Gondwana, Tectonophysics 461, 166-180, Murphy, BJ" Dostal,j" Keppie, D" 2008b, Neoprotherozoic-EarlyDevonian magmatism in the Antigonish Highlands, Avalon terrane, Nova Scotia: tracking the evolution of the mantle and crustal sources during the evolution of the Rheic Ocean, Tectonophysics 461, 181-261, Nance, R.D" Gutierrez-Alonso, G" Keppie,J,D" Linnemann, U" Murphy, J,B" Quesada, C, Strahan, RA, Woodcock, N.H., 2010, Evolution of the Rheic Ocean, Gondwana Research 17 (2-3), 194-222 (this issue), Navidad, M" Barnolas, A, 1991, El magmatismo (Ortogneises y volcanismo del Ordovicico Superior) del Paieozoico de los CataIanides, Boletin Geologico y Minero 102, 187-202, Navidad, M, Carreras, J" 2002, El volcanismo de la base del Paleozoico Inferior del macizo del Canigo (Pirineos Orientales), Evidencias geoquimicas de la apertura de una cuenca continental. Geogaceta 32, 91-94, Pearce,jA, Nigel, B., Harris, W" Tindle,AG" 1984, Trace element discrimination diagrams for the tectonic interpretation of granitic rocks, journal of Petrology 25, 956-983, Perez-Estatin, A, Bastida, E, Martinez-catalan, J,R., Gutierrez-Marco, J,c, Marcos, A, Pulgar, jA, 1990, West Asturias-Leonese Zone: stratigraphy, In: Dallmeyer, R.D" Martinez-Garcia, E (Eds,), Pre-Mesozoic Geology of Iberia, Springer Verlag, Berlin, pp, 92-102, Pin, e, Marini, E, 1993, Early Ordovician continental break-up in Variscan Europe: Nd­ Sr isotope and trace element evidence from bimodal igneous associations of the Southern Massif Central, France, Lithos 29, 177-196, Robardet, M, 2002 Alternative approach to the Variscan Belt in southwestern Europe: preorogenic peleobiogeographical constraints, In: Martinez C'atalan, JR, Hatcher, RD., Arenas, R, Diaz Garda, E (Eds,), Variscan-Appalachian Dynamics: The Building of the Llte Paleozoic Basement: Geological Society of America, Special Paper, voL 364, pp, 1-15, Robert, J,E, Thiebaut, j" 1976, Decouverte d'un volcanisme acide dans le C'aradoc de la region de Ribes de Fresser (Prov, de Gerone), Comptes Rendus de I'Academie des Sciences de Paris 282, 2050-2079, Rodriguez Alonso, MD" Peinado, M, L6pez Plaza, M, Franco, p" Carnicero, A, Gonzalo, J,e, 2004, Neoproterozoic-Cambrian synsedimentary magmatism in the Central Iberian Zone (Spain): geology, petrology and geodynamic significance, International journal Earth Sciences 93, 897-920, Roger, E, Respaut, J,P" BruneI, M, Matte, Ph" Paquette, J,L, 2004, Premiere datation U­ Pb des orthogneiss oeilles de la zone axiale de la Montgane Noire (Sud du Massif central): nouveaux temoins du magmatisme ordovicien dans la chalne varisque, Comptes Rendus Geoscience 336, 19-28. Romer, RL, Soler, A, 1995, U-Pb age and lead isotopic characterization of Au-bearing skarn related to the Andorra granite, Mineralium Deposita 30, 374-383, Santanach, p" 1972, Estudio tect6nico del Paleozoico inferior del Pirineo entre la Cerdana y el riD Tee Acta Geol6gica Hispanica 7, 44-49, Sarmiento, G,N" Sanz-L6pez,j" Barnolas, A, 1995, Conodontos delAshgill en las C'alizas de Madremanya, Les Gavarres (Girona), XIjornadas de Paleontologia, Tremp 161-163, Sassi, EP" Zirpoli, G" 1979, Pre-Variscan acidic magmatism in the Eastern Alps, IGCP no 5, Newsletter 1, 111-121, Schaltegger, U,' Abrecht,j" Corfu, E, 2003, The Ordovician orogeny in the Alpine basement: constraints from geochronology and geochemistry in the Ailr Massif (Central Alps), Schweizerische Mineralogische und Petrologische Mitteilunge 83, 183-195, Stacey, j.5., Kramers,J,D" 1975, Approximation of terrestrial lead isotope evolution by a two-stage model. Earth and Planetary Science Letters 26, 207, Stampfli, G,M" von Raumer, J,E, Borel, G,D" 2002, Paleozoic evolution of pre-Variscan terranes: from Gondwana to the Variscan collision, In: Martinez C'atalan, J,R, Hatcher, RD" Arenas, R, Diaz Garcia, E (Eds,J, Variscan-Appalachian Dynamics: The Building of the Llte Paleozoic Basement: Geological Society of America, Special Paper, voL 364, pp, 263-280, Steiger, RH., jager, E., 1977, Subcomission on geochronology: convention on the use of decay constants in geo- and cosmochronology, Earth and Planetary Science Letters 36, 359-362, Trombetta, A, Orrincione, R, Corfu, E, Mazzoleni, p" Pezzino, A, 2004, Mid-Ordovician U-Pb ages porphyroids in the Peloritan Mountains (NE Sicily): palaeogeographical implications for the evolution of the Alboran microplate, journal of the Geological Society London 161, 265-276, Twist, D" Harmer, REJ" 1987, Geochemistry of contrasting siliceous magmatic suites in the Bushveld Complex: genetic aspects and implications for tectonic discrimination diagrams, journal of Volcanology and Geothermal Research 32, 83-98. Valverde-Vaquero, p" Marcos, A, Farias, p" Gallastegui, G" 2005, U-Pb dating of Ordovician felsic volcanism in the Schistose Domain of the Galicia-Tras-os-Montes Zone near C'abo Ortegal (NW Spain), Geologica Acta 3, 27-37, Valverde-Vaquero, p" Farias, p" Marcos, A, Gallastegui, G" 2007, U-Pb dating of Siluro­ Ordovician volcanism in the Verin synform (Orense: Schistose Domain, Galicia­ Tras-os-Montes zone), Geogaceta 41, 247-250, Villas, E., Duran, H., julivert, M" 1987, The Upper Ordovician clastic sequence of the Catalan Coastal Ranges and its brachiopod fauna, Neues jahrbuch fUr Geologie und Palaontologie Abhandlungen 174, 55-74, Villaseca, e, Barbero, L, Rogers, G" 1998. Crustal origin of Hercynian peraluminous granitic batholiths of Central Spain: petrological, geochemical and isotopic (Sr, Nd) constraints, Lithos 43, 55-79, von Raumer,J,E, 1998. The Paleozoic evolution in the Alps: from Gondwana to Pangea, Geologische Rundschau 87, 407-435, von Raumer, J,E, Stampfli, G,M" Borel, G" Busy, E, 2002, Organization of pre-Variscan basement areas at the north-Gondwanan margin, International journal of Earth Sciences 91, 35-52, Winchester, jA, Floyd, PA, 1977, Geochemical discrimination of different magma series and their differentiation products using immobile elements, Chemical Geology 20, 325-343, Wood, DA, joron, J,L, Treuil, M, 1979, A re-appraisal of the use of trace elements to classify and discriminate between magma series erupted in different tectonic settings, Earth and Planetary Science Letters 45, 326-336, Zurbriggen, R, Franz, L, Handy, MR, 1997, Pre-Variscan deformation, metamorphism and magmatism in the Strona-Ceneri Zone (southern Alps of northern Italy and southern Switzerland), Schweizerische Mineralogische und Petrologische Mittei­ lungen 77, 361-380, Zwart, HJ" 1979, The geology of the Central Pyrenees, Leidse Geologische Mededelin­ gen 50, 1-74, Navidad et al 2010_Página_01 Navidad et al 2010_Página_02 Navidad et al 2010_Página_03 Navidad et al 2010_Página_04 Navidad et al 2010_Página_05 Navidad et al 2010_Página_06 Navidad et al 2010_Página_07 Navidad et al 2010_Página_08 Navidad et al 2010_Página_09 Navidad et al 2010_Página_10 Navidad et al 2010_Página_11 Navidad et al 2010_Página_12 Navidad et al 2010_Página_13