Master Robin Honlet
email: robin.honlet@kuleuven.be

PhD Marta Gasparrini

Professor Hartmut Jäger

Professor Philippe Muchez

Professor Rudy Swennen

& Swennen

Precursor and ambient rock paleothermometry to assess the thermicity of burial dolomitization in the southern Cantabrian Zone

Keywords: paleothermometry, Cantabrian Zone, dolomitization, Rock-Eval pyrolysis, vitrinite reflectance Paleothermometry, Cantabrian Zone, dolomitization, Rock-Eval pyrolysis, vitrinite reflectance

Precursor and ambient rock paleothermometry to assess the thermicity of burial dolomitization in the southern

siliciclastic rocks experienced more diffuse fluid circulation and heat dissipation. The different paleothermometry datasets compiled for the study area indicate that the fluids circulating during Late-to Post-Variscan times, with associated fracture-related dolomitization and ore mineralization in carbonate precursors, are hydrothermal. The highest paleotemperatures were recorded in ambient and precursor rocks in the highly tectonized northern part of the study area, where several thrusts and faults allowed intense fluid circulation. Positive temperature anomalies within the precursor carbonates correlate well with the occurrence of dolomite geobodies and ore mineral deposits. Such anomalies could thus be used as an exploration tool for hydrothermal dolomite bodies in analogue sub-surface settings.

Powered by Editorial Manager® and ProduXion Manager® from Aries Systems Corporation

Introduction

The Cantabrian Zone (CZ) in Northern Spain represents the Variscan foreland fold-and-thrust belt (FFTB) on the Iberian Peninsula. This most external zone of the Iberian Variscan Massif has been affected by several diagenetic to epizonal thermal events during its complex geodynamic evolution. Pre-Variscan burial metamorphism is responsible for the recorded thermal isograds in the northern and western CZ [START_REF] Bastida | Tectono-thermal evolution in a region with thinskinned tectonics : the western nappes in the Cantabrian Zone (Variscan belt of NW Spain)[END_REF][START_REF] Brime | Transition from Diagenesis to Metamorphism Near the Front of the Variscan Regional Metamorphism (Cantabrian Zone, Northwestern Spain)[END_REF]). In the southern CZ, the burial metamorphism has been overprinted by two thermal events related to Lateto Post-Variscan orocline formation and subsequent lithospheric delamination [START_REF] Raven | Metamorphic fluids and transtension in the Cantabrian Mountains of northern Spain: an application of the conodont colour alteration index[END_REF][START_REF] Gutiérrez-Alonso | Orocline triggered lithospheric delamination[END_REF][START_REF] Aller | Superposition of tectono-thermal episodes in the southern Cantabrian Zone (foreland thrust and fold belt of the Iberian Variscides, NW Spain)[END_REF]. The latter induced the circulation of hot fluids in an extensional setting which caused massive fracture-related dolomitization of Paleozoic limestones, creating seismic-scale porous dolomite geobodies (Gasparrini et al. 2006b;[START_REF] Lapponi | Hydrothermal dolomitization in a complex geodynamic setting (Lower Palaeozoic, northern Spain)[END_REF][START_REF] Muñoz-Quijano | Hydrothermal dolomitization of Paleozoic successions in Northern Spain: petrophysical properties and structural control[END_REF]. These geobodies are now exposed and represent world class outcrop analogues for dolomitic subsurface oil and gas reservoirs elsewhere [START_REF] Phipps | Exploring for dolomitized Slave Point carbonates in northeastern British Columbia[END_REF][START_REF] Hurley | Albion-Scipio and Stoney Point fields[END_REF][START_REF] Cantrell | Genesis and characterization of dolomite, Arab-D Reservoir, Ghawar Field, Saudi Arabia[END_REF][START_REF] Feng | Evidences for hydrothermal dolomite of Sinian Dengying Formation in Gaoshiti-Moxi area, Sichuan Basin[END_REF].

A variety of paleothermometers has been applied to deduce the thermal history of the CZ, with the southern thrust units in particular. They include the conodont color alteration index (CAI; [START_REF] Raven | Metamorphic fluids and transtension in the Cantabrian Mountains of northern Spain: an application of the conodont colour alteration index[END_REF][START_REF] Bastida | Tectono-thermal evolution in a region with thinskinned tectonics : the western nappes in the Cantabrian Zone (Variscan belt of NW Spain)[END_REF][START_REF] Brime | Transition from Diagenesis to Metamorphism Near the Front of the Variscan Regional Metamorphism (Cantabrian Zone, Northwestern Spain)[END_REF][START_REF] Aller | Superposition of tectono-thermal episodes in the southern Cantabrian Zone (foreland thrust and fold belt of the Iberian Variscides, NW Spain)[END_REF], the Kübler Index (KI) of illite crystallinity [START_REF] Brime | Postdepositional transformations of clays in Palaeozoic rocks of Northwest Spain[END_REF][START_REF] Aller | Cleavage and its relation with metamorphic grade in the Cantabrian Zone (Hercynian of North-West Spain)[END_REF][START_REF] Marschik | Der Übergang von Diagenese zur sehr niedergradigen Metamorphose im externen Variszikum (Kantabrische Zone)[END_REF][START_REF] Bastida | Tectono-thermal evolution in a region with thinskinned tectonics : the western nappes in the Cantabrian Zone (Variscan belt of NW Spain)[END_REF][START_REF] Brime | Transition from Diagenesis to Metamorphism Near the Front of the Variscan Regional Metamorphism (Cantabrian Zone, Northwestern Spain)[END_REF], vitrinite reflectance analysis (VR; [START_REF] Colmenero | Coal basins in the Cantabrian Mountains, Northwestern Spain[END_REF][START_REF] Frings | Coalification history of the Stephanian Ciñera-Matallana pull-apart basin, NW Spain: Combining anisotropy of vitrinite reflectance and thermal modeling[END_REF][START_REF] Colmenero | Genesis and rank distribution of Upper Carboniferous coal basins in the Cantabrian Mountains, Northern Spain[END_REF], fluid inclusion (FI) microthermometry [START_REF] Paniagua | Tectonic setting, mineralogical characteristics, geochemical signatures and age dating of a new type of epithermal carbonatehosted, precious metal-five element deposits: the Villamanín area (Cantabrian Zone, northern Spain)[END_REF][START_REF] Ayllón | Re-equilibration of fluid inclusions in diagenetic-anchizonal rocks of the Ciñera-Matallana coal basin (NW Spain)[END_REF]Gasparrini et al. 2006a;[START_REF] Lapponi | Hydrothermal dolomitization in a complex geodynamic setting (Lower Palaeozoic, northern Spain)[END_REF][START_REF] Muñoz-Quijano | Hydrothermal dolomitization of Paleozoic successions in Northern Spain: petrophysical properties and structural control[END_REF] and apatite fission-track (AFT) thermochronology (Carrière 2006, Botor and[START_REF] Botor | Thermal history of the Sabero Coalfield (Southern Cantabrian Zone, NW Spain) as revealed by apatite fission track analyses from tonstein horizons: implications for timing of coalification[END_REF]. In the southern CZ fracture-related dolomites, occurring mostly in Carboniferous carbonates, formed at temperatures between 130 and 150 °C (based on FI microthermometry; Gasparrini et al. 2006a). The dolomitization has been interpreted as hydrothermal (i.e. formed at temperatures at least 10 °C higher than those of the precursor carbonate rocks; sensu [START_REF] Machel | Hydrothermal dolomitea product of poor definition and imagination[END_REF] based on CAI data from the Carboniferous precursors (Gasparrini et al. 2006b), indicating that they experienced maximum temperatures between 70 and 95 °C [START_REF] Raven | Metamorphic fluids and transtension in the Cantabrian Mountains of northern Spain: an application of the conodont colour alteration index[END_REF]. However, the temperatures derived from CAI data do not correspond to the common occurrence of ordered mixed-layer illite/smectite in ambient siliciclastic rocks of the southern CZ, indicating temperatures above 100 °C [START_REF] Eberl | Three zones for illite formation during burial diagenesis and metamorphism[END_REF][START_REF] Brime | Transition from Diagenesis to Metamorphism Near the Front of the Variscan Regional Metamorphism (Cantabrian Zone, Northwestern Spain)[END_REF][START_REF] Aller | Superposition of tectono-thermal episodes in the southern Cantabrian Zone (foreland thrust and fold belt of the Iberian Variscides, NW Spain)[END_REF]. These observations highlight the need for a coherent interpretation of the paleothermal data available for the southern CZ, in order to better constrain the thermal conditions of dolomitization, the diagenetic process of prime importance in this area. Up to date, most of the paleothermometers listed above have mainly been used to map isograds in order to understand the different metamorphic events which affected the CZ. A more exhaustive understanding of the thermal history experienced by the Carboniferous precursor and ambient rocks in the southern CZ is thus necessary in order to address the thermal constraints of the dolomitization process in such a complex geodynamic setting.

The present study extends the published paleothermometry database by performing Rock-Eval (RE) pyrolysis and VR analysis on organic-rich clay-and mudstones from the Bodón Unit, a tectonostratigraphic unit in the southern CZ. The samples were taken from ambient siliciclastic rocks stratigraphically under-and overlying the massive dolomitized carbonate units (Visean to Bashkirian) and the newly acquired data are integrated into the paleothermal framework of the southern CZ. The aim of this study is to deduce reliable maximum temperatures experienced by precursor and ambient rocks of the dolomite geobodies occurring in the study area, to improve the understanding of the paleothermal history of the southern CZ, as well as the thermal conditions (geothermal versus hydrothermal) which governed the occurrence of the massive dolomitization. Finally, the regional distribution of the thermal maturity of precursor and ambient rocks was compared with the regional distribution of hydrothermal dolomite geobodies, which allowed to correlate the spatial occurrence of dolomites with thermal anomalies in precursor and ambient rocks. This represents a potential new approach to detect hydrothermal subsurface geobodies in hydrocarbon and geothermal exploration plays.

Geological setting

Geodynamic history of the Cantabrian Zone (CZ)

The Cantabrian Mountains of Northern Spain consist of Paleozoic sedimentary rocks, locally covered by minor Mesozoic and Cenozoic successions. These Paleozoic rocks represent the FFTB of the Iberian Variscan Massif, referred to as the Cantabrian Zone [START_REF] Julivert | Décollement tectonics in the Variscan Cordillera of northwest Spain[END_REF]; figure 1A).

The Variscan Orogeny developed a N-S oriented east-verging thin-skinned mountain range, which evolved into a secondary thick-skinned orocline during the Upper Carboniferous (ca. 310-295 Ma;Guriérrez-Alonso et al. 2004[START_REF] Botor | Hydrothermal fluids influence on the thermal evolution of the Stephanian sequence, the Sabero Coalfield (NW Spain)[END_REF][START_REF] Muñoz-Quijano | Hydrothermal dolomitization of Paleozoic successions in Northern Spain: petrophysical properties and structural control[END_REF][START_REF] Weil | Kinematic constraints on buckling a lithospheric-scale orocline along the northern margin of Gondwana: A geologic synthesis[END_REF]. The lithospheric root underneath the orocline centre, roughly corresponding to the CZ, thickened, became unstable and lithospheric delamination occurred (ca. 295-285 Ma; Gutiérrez- [START_REF] Gutiérrez-Alonso | Orocline triggered lithospheric delamination[END_REF][START_REF] Gutiérrez-Alonso | Buckling an orogen: The Cantabrian Orocline[END_REF]. The delaminating lithosphere was replaced by upwelling asthenospheric material resulting in extension and increased crustal heat flow (Gutiérrez- [START_REF] Gutiérrez-Alonso | Orocline triggered lithospheric delamination[END_REF][START_REF] Weil | Kinematic constraints on buckling a lithospheric-scale orocline along the northern margin of Gondwana: A geologic synthesis[END_REF].

Precambrian rocks exposed in the core of the Narcea Antiform (figure 1A) separate the CZ from the metamorphic internal zones of the Variscan Orogen, located to the southwest. The CZ is subdivided into several tectonostratigraphic units [START_REF] Julivert | Décollement tectonics in the Variscan Cordillera of northwest Spain[END_REF]Pérez-Estaun et al. 1988; figure 1A): the Somiedo-Correcilla, La Sobia-Bodón, Aramo, Esla, Valsurvio, Ponga, Central Coal Basin, Picos de Europa and Pisuerga-Carrión units. The Somiedo-Correcilla and La Sobia-Bodón units are often subdivided into the Somiedo and Correcilla units and the La Sobia and Bodón units, respectively. This study is focused on the Bodón Unit, located south of the León Fault (a major tectonic lineament that crosscuts a large part of the CZ: figure 1A andB). Recently, the classical subdivision of the CZ has been modified by [START_REF] Alonso | Paleogeographic inversion resulting from large out of sequence breaching thrusts: The León Fault (Cantabrian Zone, NW Iberia). A new picture of the external Variscan Thrust Belt in the Ibero-Armorican Arc[END_REF], based on a new interpretation of the León Fault, being an out-of-sequence breaching thrust, duplicating tectonic units. However, in this contribution the former subdivision is used, to enable the correlation of the new thermal data obtained from this study with thermal data from previous studies. The Bodón Unit is subdivided into three structural nappes: the Gayo, Bodón and Forcada nappes [START_REF] Evers | Geology of the Leonides between the Bernesga and Porma rivers, Cantabrian Moutains, NW Spain[END_REF][START_REF]La tectónica de la Unidad de la Sobia-Bodón[END_REF]; figure 1B).

During the Cretaceous, an extensional episode occurred related to the opening of the northern Atlantic Ocean and the Bay of Biscay located northeast of the CZ [START_REF] Gong | The rotation of Iberia during the Aptian and the opening of the Bay of Biscay[END_REF]. Overall, the extensional tectonics affecting the CZ at the end of the Paleozoic and during the Mesozoic are not well constrained due to a lack of preserved sediments [START_REF] Martín-González | Interaction of two successive Alpine deformation fronts: constraints from low-temperature thermochronology and structural mapping (NW Iberian Peninsula)[END_REF].

The last tectonic event that affected the CZ is related to the Alpine Orogeny, leading to uplift and exhumation of the Variscan basement [START_REF] Alonso | Tertiary basins and Alpine tectonics in the Cantabrian Mountains (NW Spain)[END_REF][START_REF] Fillon | Alpine exhumation of the central Cantabrian Mountains, Northwest Spain[END_REF]. This created much of the present topography (i.e. the Cantabrian Mountains). The Alpine Orogeny caused minor internal deformation in the CZ, except for slight reactivation of faults and thrusts, block tilting and some further shortening [START_REF] Pulgar | La deformación alpina en al basamento varisco de la Zona Cantábrica[END_REF].

Stratigraphy of the Bodón Unit

In the Bodón Unit, as in most tectonostratigraphic units of the western and southern CZ, Paleozoic rocks were thrust along a décollement level located in carbonate rocks of the Lower Cambrian Láncara Fm. [START_REF] Julivert | Décollement tectonics in the Variscan Cordillera of northwest Spain[END_REF]). The CZ is subdivided into a pre-orogenic (pre-Carboniferous) and a synorogenic rock succession (Carboniferous), developed during the Variscan Orogeny [START_REF] Marcos | An approach to the tectonostratigraphic evolution of the Cantabrian Foreland thrust and fold belt, Hercynian Cordillera of NW Spain[END_REF]; figure 2). Due to the Variscan and Alpine compression, the rock successions in the Bodón Unit were tilted and presently dip subvertically, leading to excellent outcropping conditions, particularly along road cuts (figure 3).

The Láncara Fm. consists of limestones with a lower member composed of synsedimentary to early diagenetic dolomite [START_REF] Zamarreño | Las lithofacies carbonatadas del Cámbrico de la Zona Cantábrica (NW España) y su distribución paleogeográfica[END_REF]. Upper Cambrian to Silurian deposits are composed of siliciclastic rocks, with the Silurian Formigoso Fm. composed of dark shales. The Devonian is characterized by an alternation of carbonates and siliciclastic rocks, with the La Vid Group representing the thickest stratigraphic unit (figure 3A). The Ermita Fm. is a transgressive sandstone unit of a few meters thick, deposited during the Upper Devonian (figure 3A). The Vegamián Fm. at the base of the Lower Carboniferous is a thin succession (5 to 10 m) of black shales and siltstones with phosphate, manganese, markasite and chert nodules (figure 3A). Thin silt-to sandstone intercalations can occur, as well as carbonate lenses of a few cm in thickness. The Vegamián Fm. was deposited below wave base at low sedimentation rates under anoxic conditions [START_REF] Raven | Conodont biostratigraphy and depositional history of the Middle Devonian to Lower Carboniferous in the Cantabrian Zone (Cantabrian Mountains, Spain)[END_REF]. [START_REF] Seibert | Fazies und Paläogeographie des Unter-Karbon (Alba Formation) im Kantabrischen Gebirge[END_REF] correlates the deposition of the Vegamián Fm. with eustatic sea-level rise and subsidence. The Vegamián Fm. represents the transition to thick carbonate successions filling the Variscan foreland basin during the Lower to Middle Carboniferous [START_REF] Bahamonde | Pennsylvanian carbonate platforms adjacant to deltaic systems in an active marine foreland basin (Escalada Fm., Cantabrian Zone, NW Spain)[END_REF]. These successions are divided into three formations: the Alba, Barcaliente and Valdeteja fms. The Alba Fm. (Late Tournaisian to Late Serpukhovian) is usually less than 30 m in thickness and is composed of two nodular limestone packages, interpreted as condensed units deposited on a well oxygenated carbonate platform (figure 3A). A member composed of bedded cherts can be found between these limestone packages and represents an episode of pelagic sedimentation [START_REF] Wagner | Lithostratigraphic units of the lower part of the Carboniferous in Northern León, Spain[END_REF][START_REF] Sanchez De Posada | Carboniferous Pre-Stephanian Rocks of the Asturian-Leonese Domain (Cantabrian Zone)[END_REF]). The overlying Barcaliente Fm. (Late Serpukhovian to earliest Bashkirian) is composed of thinly bedded dark mudstones showing a remarkable consistency in both thickness (ca. 300 m) and depositional facies over most of the CZ. It reflects a period of uniform carbonate sedimentation in a shallow water, low energy, anoxic environment [START_REF] Evers | Geology of the Leonides between the Bernesga and Porma rivers, Cantabrian Moutains, NW Spain[END_REF][START_REF] Wagner | Lithostratigraphic units of the lower part of the Carboniferous in Northern León, Spain[END_REF][START_REF] Dietrich | Numerical modeling as a Means to Enhance Genetic Sedimentary Basin Interpretation[END_REF]. The Bashkirian Valdeteja Fm. (up to 1500 m thick) is made of massive carbonates deposited in a shallow marine platform environment [START_REF] Prins | Carboniferous Productidina and Chonetidina of the Cantabrian Mountains (NW Spain): Systematics, Stratigraphy and Palaeoecology[END_REF][START_REF] Wagner | Lithostratigraphic units of the lower part of the Carboniferous in Northern León, Spain[END_REF][START_REF] Dietrich | Numerical modeling as a Means to Enhance Genetic Sedimentary Basin Interpretation[END_REF][START_REF] Chesnel | Facies, geometry and growth phases of the Valdorria carbonate platform (Pennsylvanian, Northern Spain)[END_REF]; figure 3B), whose sedimentation is strongly influenced by the approaching Variscan orogenic front. The Valdeteja Fm. is diachronically overlain by siliciclastic rocks of the Late Bashkirian to Early Moscovian San Emiliano Fm. [START_REF] Lobato | Hoja del Mapa Geológico Nacional de España[END_REF][START_REF] Alonso | Hoja del Mapa Geológico Nacional de España (1:50000)[END_REF]; figure 3B), which filled up the foredeep of the Variscan Orogen [START_REF] Brouwer | La succession Carbonifère dans la partie méridionale des montagnes Cantabriques (Espagne Nord-Ouest)[END_REF]. It is composed of claystones (figure 3C), locally with abundant sandstone and carbonate lenses, and coal seams (van Ginkel 1965;[START_REF] Wagner | The position of the Bashkirian/Moscovian boundary in West European chronostratigraphy[END_REF][START_REF] Eichmüller | Die Valdeteja Formation: Aufbau und Geschichte einer oberkarbonischen Karbonatplattform (Kantabrisches Gebirge, Nordspanien)[END_REF][START_REF] Samankassou | Internal structure and depositional environment of Late Carboniferous mounds from the San Emiliano Formation, Cármenes Syncline, Cantabrian Mountains, Northern Spain[END_REF]. Unconformable Stephanian siliciclastic successions occur in local intramontane pull-apart basins (figure 1A), of which the Ciñera-Matallana (CM) and Canseco-Rucayo (CR) basins are located closest to the Bodón Unit [START_REF] Colmenero | Genesis and rank distribution of Upper Carboniferous coal basins in the Cantabrian Mountains, Northern Spain[END_REF]; figure 1B). Permian to Cretaceous sediments are absent in the southern CZ and crop out only locally in the northern CZ (Dalmeyer and Martínez-García 1990; figure 1A). They might have been initially deposited and eroded later on [START_REF] Alonso | Tertiary basins and Alpine tectonics in the Cantabrian Mountains (NW Spain)[END_REF][START_REF] Gómez-Fernández | Metallogenesis of Zn-Pb carbonate-hosted mineralization in the southeastern region of the Picos de Europa (Central Northern Spain) Province: geologic, fluid inclusion, and stable isotopes studies[END_REF][START_REF] Frings | Coalification history of the Stephanian Ciñera-Matallana pull-apart basin, NW Spain: Combining anisotropy of vitrinite reflectance and thermal modeling[END_REF][START_REF] Botor | Hydrothermal fluids influence on the thermal evolution of the Stephanian sequence, the Sabero Coalfield (NW Spain)[END_REF]).

In the Bodón Unit, fracture-related dolomitization mostly affected the Barcaliente and Valdeteja Fms. (Gasparrini et al. 2006b;[START_REF] Muñoz-Quijano | Hydrothermal dolomitization of Paleozoic successions in Northern Spain: petrophysical properties and structural control[END_REF]figure 3D). Smaller dolomite geobodies occur in the Láncara Fm. [START_REF] Lapponi | Hydrothermal dolomitization in a complex geodynamic setting (Lower Palaeozoic, northern Spain)[END_REF]) and in the top layers of the Alba Fm. (Gasparrini et al. 2006b). Detailed petrographic descriptions of dolomite exposed in the Bodón Unit have been made by Gasparrini et al. (2006b) and [START_REF] Lapponi | Hydrothermal dolomitization in a complex geodynamic setting (Lower Palaeozoic, northern Spain)[END_REF]. Both studies describe dolomite rocks with a typical succession of replacive and void-filling dolomite, followed by calcite cementation. Based on the work of Muñoz Quijano (2015) combined with results from aerial photographs and mapping, the occurrence of dolomite geobodies in Carboniferous carbonates has been indicated on figure 1B. Important ore mineralizations post-date dolomitization, as ore minerals occur in pores and cavities created by the dolomitization process [START_REF] Paniagua | Tectonic setting, mineralogical characteristics, geochemical signatures and age dating of a new type of epithermal carbonatehosted, precious metal-five element deposits: the Villamanín area (Cantabrian Zone, northern Spain)[END_REF]Gasparrini et al. 2006b).

Thermal history of the CZ

The Paleozoic rocks of the CZ have been affected by several thermal episodes related to the complex geodynamic history of the area [START_REF] Bastida | Tectono-thermal evolution of the Cantabrian Zone (NW Spain)[END_REF]. Peak burial prior to the onset of Variscan thrusting was the most important process causing the metamorphic grades recorded by the Paleozoic rocks in most of the CZ [START_REF] Bastida | Tectono-thermal evolution of the Cantabrian Zone (NW Spain)[END_REF][START_REF] Aller | Superposition of tectono-thermal episodes in the southern Cantabrian Zone (foreland thrust and fold belt of the Iberian Variscides, NW Spain)[END_REF]. This is indicated by the following: (1) metamorphic grades increase with rock age, (2) the distribution of metamorphic grades is not influenced by major thrusts and (3) metamorphic patterns are discontinuously inverted (i.e. rocks recording higher metamorphic grades are thrust over rocks recording lower metamorphic grades; [START_REF] Brime | Postdepositional transformations of clays in Palaeozoic rocks of Northwest Spain[END_REF][START_REF] Bastida | Tectono-thermal evolution in a region with thinskinned tectonics : the western nappes in the Cantabrian Zone (Variscan belt of NW Spain)[END_REF][START_REF] Aller | Superposition of tectono-thermal episodes in the southern Cantabrian Zone (foreland thrust and fold belt of the Iberian Variscides, NW Spain)[END_REF]. Burial resulted in diagenetic (up to 200 °C) to locally anchizonal (200 to 300 °C) conditions. The thermal peak occurred when the Paleozoic successions reached their maximum thickness (Bashkirian-Moscovian), preceding the emplacement of the tectonostratigraphic units [START_REF] Bastida | Tectono-thermal evolution of the Cantabrian Zone (NW Spain)[END_REF]. Based on CAI data from Devonian rocks in the Somiedo Unit, a paleogeothermal gradient of about 35 °C/km is estimated for this episode of peak burial [START_REF] Brime | Transition from Diagenesis to Metamorphism Near the Front of the Variscan Regional Metamorphism (Cantabrian Zone, Northwestern Spain)[END_REF].

Regional orogenic metamorphism, due to Variscan folding and thrusting, has only locally been recognized in the CZ [START_REF] Bastida | Tectono-thermal evolution of the Cantabrian Zone (NW Spain)[END_REF]. It is generally recorded closer to the internal areas of the Variscan Orogeny (e.g. in the northwestern part of the Somiedo Unit; [START_REF] Brime | Transition from Diagenesis to Metamorphism Near the Front of the Variscan Regional Metamorphism (Cantabrian Zone, Northwestern Spain)[END_REF]. It is of minor importance in the southern CZ and did not affect the Paleozoic rocks presently exposed in the Bodón Unit [START_REF] Bastida | Tectono-thermal evolution of the Cantabrian Zone (NW Spain)[END_REF][START_REF] Aller | Superposition of tectono-thermal episodes in the southern Cantabrian Zone (foreland thrust and fold belt of the Iberian Variscides, NW Spain)[END_REF].

Lithospheric delamination, initiated after the orocline formation, led to the upwelling of the asthenosphere and to a concomitant increase in crustal heat flow (Gutiérrez- [START_REF] Gutiérrez-Alonso | Orocline triggered lithospheric delamination[END_REF][START_REF] Weil | Kinematic constraints on buckling a lithospheric-scale orocline along the northern margin of Gondwana: A geologic synthesis[END_REF]. The upwelling also caused uplift of the CZ, accompanied by near-surface extension (Muñoz-Quijano and Gutiérrez-Alonso 2007). Two thermal events have been recognized in the CZ, corresponding respectively to the orocline formation and the subsequent lithospheric delamination, both accompanied by dominant extensional tectonics [START_REF] Bastida | Tectono-thermal evolution of the Cantabrian Zone (NW Spain)[END_REF][START_REF] Aller | Superposition of tectono-thermal episodes in the southern Cantabrian Zone (foreland thrust and fold belt of the Iberian Variscides, NW Spain)[END_REF].

The first event, often referred to as the Late-Variscan extensional event [START_REF] Bastida | Tectono-thermal evolution of the Cantabrian Zone (NW Spain)[END_REF]) and likely associated with the orocline formation, has been recorded in the southern part of the Central Coal Basin and in the Pisuerga-Carrión and Valsurvio units [START_REF] Aller | Deformación y metamorfismo en la parte Sur de la Cuenca Carbonífera Central (NO de España)[END_REF][START_REF] Aller | Cleavage and its relation with metamorphic grade in the Cantabrian Zone (Hercynian of North-West Spain)[END_REF][START_REF] Bastida | Tectono-thermal evolution of the Cantabrian Zone (NW Spain)[END_REF]. Here, the Paleozoic rocks record anchizonal to locally epizonal conditions (above 300 °C) typically associated with subhorizontal cleavage, as a result of extensional tectonics [START_REF] Aller | Cleavage and its relation with metamorphic grade in the Cantabrian Zone (Hercynian of North-West Spain)[END_REF]Marín 1997;[START_REF] García-López | Los episodios metamórficos de la Zona Cantábrica y su contexto estructural[END_REF]. The distribution of both anchizonal and epizonal metamorphism (affecting pre-Permian rocks only), as well as the cleavage development, do not correspond to the location of the main Variscan structures. Intrusions of small sills and stocks, which occur mainly in the western and eastern units of the CZ [START_REF] Corretgé | Igneous rocks of the Cantabrian/Palentine Zone[END_REF], accompanied this event. Close to these intrusions, contact metamorphism resulted in epizonal grades and the development of skarn deposits (Bastida et al. 2002 and references therein). Since subhorizontal cleavage associated with high metamorphic grades or intrusions has not been reported in the Bodón Unit [START_REF] Aller | Cleavage and its relation with metamorphic grade in the Cantabrian Zone (Hercynian of North-West Spain)[END_REF], the Late-Variscan extensional event was likely of minor importance.

The second thermal event, referred to as the Post-Variscan thermal event [START_REF] Bastida | Tectono-thermal evolution of the Cantabrian Zone (NW Spain)[END_REF], was responsible for the widespread occurrence of fracture-related dolomitization and ore mineralization in Paleozoic successions [START_REF] Martínez-García | Permian Mineralizations in the Cantabrian Mountains (North-West Spain)[END_REF][START_REF] Gómez-Fernández | Metallogenesis of Zn-Pb carbonate-hosted mineralization in the southeastern region of the Picos de Europa (Central Northern Spain) Province: geologic, fluid inclusion, and stable isotopes studies[END_REF][START_REF] Crespo | The Salamón gold deposit (León, Spain)[END_REF]Gasparrini et al. 2006b;[START_REF] Muñoz-Quijano | Modelo de evolución topográfica en el NO de la Península Ibérica durante la delaminación litosférica al final de la Orogenia Varisca[END_REF]. Several anomalies in CAI values occur within the Bodón and Picos de Europa units, in areas with dense fault patterns and abundant ore mineralization, and are interpreted to be caused by hot fluid movements during this thermal episode [START_REF] Raven | Metamorphic fluids and transtension in the Cantabrian Mountains of northern Spain: an application of the conodont colour alteration index[END_REF][START_REF] Bastida | Tectono-thermal evolution of the Cantabrian Zone (NW Spain)[END_REF]. In the southern CZ, anchi-to epizonal conditions recorded in the Stephanian pull-apart basins are probably related to intrusions of diorite during the same thermal episode [START_REF] Bastida | Tectono-thermal evolution of the Cantabrian Zone (NW Spain)[END_REF]. These intrusions resulted in reequilibration of fluid inclusions in quartz veins in the CM Basin [START_REF] Ayllón | Re-equilibration of fluid inclusions in diagenetic-anchizonal rocks of the Ciñera-Matallana coal basin (NW Spain)[END_REF]. The Post-Variscan thermal event is dated at 270 Ma based on K-Ar ages from illite [START_REF] Weh | The late orogenic history of the southeastern Cantabrian Mountains: illite-crystallinity and K-Ar data[END_REF]). These ages correspond with U-Pb and Pb-Pb absolute dating of dolomite-hosted ore minerals from the Bodón Unit [START_REF] Paniagua | Tectonic setting, mineralogical characteristics, geochemical signatures and age dating of a new type of epithermal carbonatehosted, precious metal-five element deposits: the Villamanín area (Cantabrian Zone, northern Spain)[END_REF]. Therefore, the Post-Variscan thermal event correlates temporally with lithospheric delamination following the orocline formation. Lithospheric delamination induced the circulation of hot and hypersaline marine-derived brines, modified through water-rock interactions, representing the source of Mg for massive dolomitization (Gasparrini et al. 2006a, b;[START_REF] Lapponi | Hydrothermal dolomitization in a complex geodynamic setting (Lower Palaeozoic, northern Spain)[END_REF]. Based on AFT thermochronology, [START_REF] Carrière | Neoproterozoic to Holocene tectonothermal evolution of the southern Cantabrian Mountains NW Iberia, revealed by apatite fission-track thermochronology[END_REF] and [START_REF] Botor | Thermal history of the Sabero Coalfield (Southern Cantabrian Zone, NW Spain) as revealed by apatite fission track analyses from tonstein horizons: implications for timing of coalification[END_REF] concluded that although this thermal event was likely short-lived, temperatures remained high within Stephanian siliciclastic successions (above 60 °C), at least until the Middle Triassic.

Materials and methods

Published paleothermometry database

In the following sections, the different paleothermometry datasets published for the Bodón Unit in previous studies are introduced. The stratigraphic intervals sampled for each paleothermometer are indicated on the stratigraphic column of figure 2. A brief explanation of the different paleothermometers is followed by an explanation on how the raw data have been converted to paleotemperatures. Most of the introduced paleothermometers are strongly affected by the duration of a thermal event, which has to be taken into account when comparing different paleothermal datasets. What also needs attention, is the fact that some datasets are available only from the precursor carbonate successions while others are derived from the siliciclastic ambient rocks. The calculated paleotemperatures have to be handled with care, because the empirical equations used for the different datasets have a specific accuracy and uncertainty regarding the parameters and preconditions applied, which might be a source of error and inaccuracy.

 Conodont Color Alteration Index (CAI)

The color alteration of conodonts depends on the carbonization of organic matter (OM) traces and is controlled by maximum temperature and heating time [START_REF] Epstein | Conodont color alterationan index to organic metamorphism[END_REF][START_REF] Rejebian | Conodont color and textural alteration: An index to regional metamorphism, contact metamorphism, and hydrothermal alteration[END_REF]. Recalculated paleotemperatures are strongly influenced by the duration of the heating event, as a relatively long heating time is required for the carbonization of organic matter in conodonts. Confining pressure does not affect the alteration process, excluding tectonics from influencing the color of conodont fragments [START_REF] Epstein | Conodont color alterationan index to organic metamorphism[END_REF]. Nevertheless, the presence of fluids under pressure, which is important in hydrothermal systems, can retard alteration. In such settings, OM carbonization is replaced by oxidation and volatilization of oxides, resulting in a loss of OM which affects the color alteration [START_REF] Epstein | Conodont color alterationan index to organic metamorphism[END_REF][START_REF] Rejebian | Conodont color and textural alteration: An index to regional metamorphism, contact metamorphism, and hydrothermal alteration[END_REF].

CAI data were obtained by [START_REF] Raven | Metamorphic fluids and transtension in the Cantabrian Mountains of northern Spain: an application of the conodont colour alteration index[END_REF] and [START_REF] Aller | Superposition of tectono-thermal episodes in the southern Cantabrian Zone (foreland thrust and fold belt of the Iberian Variscides, NW Spain)[END_REF] and are derived from conodonts found in Devonian and Carboniferous rock successions. The data have been converted to temperature using the Arrhenius plot constructed by [START_REF] Epstein | Conodont color alterationan index to organic metamorphism[END_REF]. A heating time of 10 Ma was used to convert the CAI data, based on the AFT studies of [START_REF] Carrière | Neoproterozoic to Holocene tectonothermal evolution of the southern Cantabrian Mountains NW Iberia, revealed by apatite fission-track thermochronology[END_REF] and [START_REF] Botor | Thermal history of the Sabero Coalfield (Southern Cantabrian Zone, NW Spain) as revealed by apatite fission track analyses from tonstein horizons: implications for timing of coalification[END_REF].

 Kübler index (KI) of illite crystallinity and clay mineralogy

The KI of illite crystallinity allows the determination of the kinetic clay mineral reaction progress in pelitic sequences. The crystallinity of illite is determined by measuring the full width at half maximum height of the illite (001) X-ray diffraction peak [START_REF] Kübler | La cristallinité de l'illite et les zones tout à fait supérieures du métamorphisme[END_REF][START_REF] Kübler | Evaluation quantitative du métamorphisme par cristallinité de l'illite[END_REF]). Measurements are performed on oriented mineral aggregate samples of the < 2 µm grain size fractions. Increasing crystallinity results in decreasing KI values. Despite the important role of temperature in the reaction process of phyllosilicates, numerous other factors (such as permeability, potassium availability, weathering and mechanical and hydrothermal alteration) affect the crystallinity of illite (e.g. [START_REF] Junfeng | Relationschip between illite crystallinity and temperature in active geothermal systems of New Zealand[END_REF][START_REF] Frings | Hydrothermally altered mudrock of the Ciñera-Matallana coal basin, Cantabrian Zone, northern Spain[END_REF][START_REF] Fukuchi | Changes in illite crystallinity within an ancient tectonic boundary thrust caused by thermal, mechanical, and hydrothermal effects: an example from the Nobeoka Thrust, southwest Japan[END_REF]). As a result, phyllosilicates in diagenetic conditions often reflect clay mineral reaction progress instead of thermodynamic equilibrium, which makes the use of the KI as a paleothermometer problematic [START_REF] Guggenheim | Report of the Association Internationale pour l'Etude des Argiles (AIPEA) Nomenclature Committee for 2001: Order, disorder and crystallinity in phyllosilicates and the use of the 'Crystallinity Index[END_REF].

Despite the many factors possibly influencing the crystallinity of illite, general correlation charts (e.g. [START_REF] Merriman | Patterns of very low-grade metamorphism in metapelitic rocks[END_REF][START_REF] Voldman | Conodont palaeothermometry of contact metamorphism in Middle Ordovician rocks from the Precordillera of western Argentina[END_REF] estimate that KI values lower than 0.42 correspond to temperatures above 200 °C (anchizone), whereas KI values lower than 0.25 correspond to temperatures above 300 °C (epizone). [START_REF] Marschik | Der Übergang von Diagenese zur sehr niedergradigen Metamorphose im externen Variszikum (Kantabrische Zone)[END_REF] and [START_REF] Aller | Superposition of tectono-thermal episodes in the southern Cantabrian Zone (foreland thrust and fold belt of the Iberian Variscides, NW Spain)[END_REF] analyzed the clay mineralogy and KI of clay-rich siliciclastic rocks in the Bodón and Correcilla units, ranging from Cambrian to Carboniferous in age. Corresponding temperatures were calculated using the empirical formula of [START_REF] Mukoyoshi | Quantitative estimation of temperature conditions for illite crystallinity: comparison to vitrinite reflectance from the Chichibu and Shimanto accretionary complexes, eastern Kyushu, Southwest Japan[END_REF]. This equation has been established based on a correlation between VR and KI data in an accretionary complex in Southwest Japan, and has been applied with extreme caution. [START_REF] Eberl | Three zones for illite formation during burial diagenesis and metamorphism[END_REF] identified 3 reaction zones for illite formation with the third zone corresponding to temperatures of at least 100 °C. Clay mineralogy of samples from the Bodón Unit suggests that all mixedlayer illite/smectite in the samples used for KI measurements appears to be ordered (e.g. [START_REF] Brime | Transition from Diagenesis to Metamorphism Near the Front of the Variscan Regional Metamorphism (Cantabrian Zone, Northwestern Spain)[END_REF][START_REF] Aller | Superposition of tectono-thermal episodes in the southern Cantabrian Zone (foreland thrust and fold belt of the Iberian Variscides, NW Spain)[END_REF]. This indicates that zone 3 of [START_REF] Eberl | Three zones for illite formation during burial diagenesis and metamorphism[END_REF] was reached [START_REF] Brime | Transition from Diagenesis to Metamorphism Near the Front of the Variscan Regional Metamorphism (Cantabrian Zone, Northwestern Spain)[END_REF][START_REF] Aller | Superposition of tectono-thermal episodes in the southern Cantabrian Zone (foreland thrust and fold belt of the Iberian Variscides, NW Spain)[END_REF].

 Vitrinite reflectance (VR) analysis

Due to changes in optical properties, the reflectance of vitrinite increases with increasing temperature in a predictable manner [START_REF] Karweil | Die Metamorphose der Kohlen vom Standpunkt der physikalischen Chemie[END_REF]Suggate 1959;[START_REF] Teichmüller | Application of coal petrological methods in geology including oil and natural gas prospecting[END_REF]. Although suppression or retardation of vitrinite reflectance can occur under certain circumstances (e.g. in case of overpressure; [START_REF] Schito | Overcoming pitfalls of vitrinite reflectance measurements in the assessment of thermal maturity: the case history of the lower Congo basin[END_REF], the analysis of vitrinite reflectance is considered a reliable paleothermometry technique to assess the thermal maturity of sediments in diagenetic settings [START_REF] Teichmüller | Organic material and very low grade metamorphism[END_REF].

Published VR data are only available for the Stephanian CM [START_REF] Frings | Coalification history of the Stephanian Ciñera-Matallana pull-apart basin, NW Spain: Combining anisotropy of vitrinite reflectance and thermal modeling[END_REF][START_REF] Colmenero | Genesis and rank distribution of Upper Carboniferous coal basins in the Cantabrian Mountains, Northern Spain[END_REF]) and CR [START_REF] Colmenero | Genesis and rank distribution of Upper Carboniferous coal basins in the Cantabrian Mountains, Northern Spain[END_REF]) basins. They have been converted to temperature by applying the formula proposed by [START_REF] Barker | Calculation of vitrinite reflectance from thermal histories and peak temperatures. A comparison of methods[END_REF] which, according to [START_REF] Frings | Coalification history of the Stephanian Ciñera-Matallana pull-apart basin, NW Spain: Combining anisotropy of vitrinite reflectance and thermal modeling[END_REF], yields the most realistic temperatures for samples from the Stephanian CM Basin.

 Fluid inclusion (FI) microthermometry

FIs are tiny vacuoles of fluid trapped in crystals during mineral precipitation [START_REF] Goldstein | Systematics of fluid inclusions in diagenetic minerals[END_REF]. They provide information on the composition and temperature of the motherfluids through microthermometric measurements. Of most importance for this study is the homogenization temperature (Th), which gives an estimation of the minimum trapping temperature of an inclusion. By applying a pressure correction, the true formation temperature of crystal precipitation can be derived from Th.

Primary FIs in fracture-related dolomites occurring in both Cambrian and Carboniferous precursors have been investigated in previous studies (Gasparrini et al. 2006a;[START_REF] Lapponi | Hydrothermal dolomitization in a complex geodynamic setting (Lower Palaeozoic, northern Spain)[END_REF][START_REF] Muñoz-Quijano | Hydrothermal dolomitization of Paleozoic successions in Northern Spain: petrophysical properties and structural control[END_REF], as well as FIs in late diagenetic calcite supposed to have precipitated from the same fluid system (Gasparrini et al. 2006a). FIs in quartz gangue minerals from ore deposits of the Providencia and Profunda mines (figure 1B) have been investigated by [START_REF] Paniagua | Tectonic setting, mineralogical characteristics, geochemical signatures and age dating of a new type of epithermal carbonatehosted, precious metal-five element deposits: the Villamanín area (Cantabrian Zone, northern Spain)[END_REF]. [START_REF] Ayllón | Re-equilibration of fluid inclusions in diagenetic-anchizonal rocks of the Ciñera-Matallana coal basin (NW Spain)[END_REF] performed a microthermometric study on quartz veins in the Stephanian CM Basin.

 Apatite fission-track (AFT) thermochronology

Fission track thermochronology is a suitable technique to investigate the low-temperature thermal history of sedimentary rocks [START_REF] Price | Fossil tracks of charged particles in mica and the age of minerals[END_REF][START_REF] Wagner | Fission-Track Dating[END_REF]. The technique is based on the natural decay of 238 U isotopes, producing radiation damage trails (so-called fission tracks) in U-bearing minerals. The spontaneous fission of 238 U isotopes produces fission tracks at a continuous rate throughout geological time. Fission-track thermochronology is mostly performed on apatite crystals (AFT thermochronology). Apatite fission tracks can undergo partial annealing (i.e. resetting) in a temperature range of 60 to 110 ± 10 °C, due to a measurable shortening of the fission tracks at these temperatures [START_REF] Gleadow | A natural long-term track annealing experiment for apatite[END_REF]. Complete annealing of fission tracks in apatite crystals occurs at temperatures higher than 110 ± 10 °C for 10 Ma. [START_REF] Carrière | Neoproterozoic to Holocene tectonothermal evolution of the southern Cantabrian Mountains NW Iberia, revealed by apatite fission-track thermochronology[END_REF] successfully deduced the thermal history for 2 Cambrian and 1 Stephanian sample in or close to the eastern part of the Bodón Unit. Recently, [START_REF] Botor | Thermal history of the Sabero Coalfield (Southern Cantabrian Zone, NW Spain) as revealed by apatite fission track analyses from tonstein horizons: implications for timing of coalification[END_REF] investigated the thermal history of the Stephanian Sabero Coalfield (located southeast of the study area) through a combination of AFT thermochronology and VR analysis.

Rock-Eval (RE) pyrolysis and vitrinite reflectance (VR) analysis

In addition to the published paleothermal datasets introduced above, the thermal maturity of a new set of samples was investigated through RE pyrolysis and VR analysis. A total of 78 samples of OM-rich clay-and mudstones from different stratigraphic levels were collected from road cuts and quarries. Most of the samples were collected in the Bodón Nappe (figure 1B). Other samples were collected in the Forcada and Gayo nappes, north of the León Fault (in the Central Coal Basin as well as in the Stephanian CR Basin) and in the northern part of the Correcilla Unit. Fresh material was preferentially sampled in order to minimize the influence of weathering (cfr. [START_REF] Copard | Erroneous maturity assessment given by Tmax and HI Rock-Eval parameters on highly mature weathered coals[END_REF].

RE pyrolysis is an analytical method used to characterize source rocks and evaluate their petroleum potential. It is based on the relative contents of free hydrocarbons and CO2 generated during pyrolysis, which is artificial heating of a sample [START_REF] Espitalié | Méthode rapide de caractérisation des roches mères, de leur potentiel pétrolier et de leur degré d'évolution[END_REF]1985a, b, c). Pyrolysis of the collected samples was performed using the Rock-Eval 6 apparatus [START_REF] Lafargue | Rock-Eval 6 applications in hydrocarbon exploration, production, and soil contamination studies[END_REF]) which employs a temperature-programmed heating of a powdered sample (± 100 mg) in an inert atmosphere (helium or nitrogen). The sample is first placed in a pyrolysis oven to measure the quantity of free hydrocarbons (revealed by the so-called S1 peak). Heating up to 650 °C allows the determination of the amount of hydrocarbons (S2 peak) and CO2 (S3 peak) produced during thermal cracking of kerogen. The sample is then transferred to an oxidation oven where the residual organic carbon, remaining in the sample after pyrolysis, is oxidized (and revealed by the S4 peak) at temperatures up to 850 °C. In this way, the Total Organic Carbon (TOC) content of the sample can be quantified. Of particular interest for this study is the Tmax parameter which corresponds to the S2 peak temperature in the pyrolysis oven, which is used to estimate the thermal maturity of rocks [START_REF] Tissot | Petroleum Formation and Occurrence. A new approach to oil and gas exploration[END_REF]. In this study, Tmax values were converted to equivalent VR values (referred to as Tmax-derived VR values) using the equation of [START_REF] Jarvie | Oil and Shale Gas from the Barnett Shale, Ft. Worth Basin, Texas[END_REF]. RE pyrolysis produced reliable results for 55 samples, which are listed in table 1. The other samples (including all samples from the Alba, Barcaliente and Valdeteja fms.) were excluded due to low TOC values (< 0.3 %) which have a high risk to produce unreliable RE pyrolysis data. From the 55 productive samples, 40 are from OM-rich layers of the San Emiliano Fm. In the Central Coal Basin, north of the León Fault, 5 samples are from the Lena and Sama Groups (local equivalents of the San Emiliano Fm.) and 2 samples are from the CR Basin. The remaining 8 samples were collected from the Formigoso (2 samples) and Vegamián fms. (4 samples), from the La Vid Group (1 sample) and from the Stephanian CM Basin (Roguera Fm.; 1 sample).

VR analysis, including specific maturity analysis of dispersed kerogen, was performed on 18 samples out of the 55 productive RE samples of the San Emiliano (15 samples) and Vegamián fms. (3 samples; table 1) through a method specifically applied at the GeoResources STC lab. The samples cover a wide range of TOC and Tmax values to evaluate the consistency between results from RE pyrolysis and VR analysis. Most samples are located within the Bodón Nappe and only 1 is located in the Forcada Nappe. Kerogen was isolated by the maceration technique (double acid treatment with HCl and HF) and concentrated and embedded in epoxy resin blocks, which were polished on the surface for reflectance analysis. VR analysis is based on digital image analysis, using high-resolution black/white images of vitrinite, obtained with a reflected-light microscope. Grey levels of the digital images represent the levels of reflectance and are recalculated to real VR values by specifically modified image analysis software. This enables high-resolution reflectance analysis of single vitrinite grains down to pixel-size (< 10 μm) providing highly accurate and reliable VR data, without any correction factor needed. It also supports the identification of in situ versus reworked and degraded vitrinite in mixed vitrinite assemblages and the separation of real vitrinite from vitrinite-like particles. In each sample 50 vitrinite particles were measured, if available. Additionally, one kerogen slide was made from each sample for optical analysis of kerogen composition, particularly the composition of the vitrinite assemblage, which was included in the detailed interpretation of the VR data. For the calculation of exact paleotemperatures from VR data several equations can be used, leading to different paleotemperatures, with increasing variation towards higher maturation (table 2). In this study the equation of [START_REF] Barker | Calculation of vitrinite reflectance from thermal histories and peak temperatures. A comparison of methods[END_REF] was used to provide the maximum comparability between the new data of this study and the previous data by [START_REF] Frings | Coalification history of the Stephanian Ciñera-Matallana pull-apart basin, NW Spain: Combining anisotropy of vitrinite reflectance and thermal modeling[END_REF].

Results

In the following sections the paleotemperatures obtained from literature datasets and from RE and VR measurements are reported. The regional distribution of paleotemperatures, calculated from each technique, is indicated on the geological map of the Bodón Unit (figures 4 and 5).

Published paleothermometry database

CAI data for Devonian and Carboniferous successions in the Bodón Unit range from 1 to 5 [START_REF] Raven | Metamorphic fluids and transtension in the Cantabrian Mountains of northern Spain: an application of the conodont colour alteration index[END_REF][START_REF] Aller | Superposition of tectono-thermal episodes in the southern Cantabrian Zone (foreland thrust and fold belt of the Iberian Variscides, NW Spain)[END_REF], corresponding to a wide range of diagenetic to anchizonal temperatures. The temperatures average around 60 °C for the Carboniferous precursor carbonates in the eastern part of the Bodón Unit, whereas precursor carbonates close to the dolomitized central and western parts are characterized by higher temperatures (up to 203 °C; figure 4A). These higher temperatures correlate with areas intensely affected by fracture-related dolomitization and ore mineralization.

Temperatures obtained from KI data from the ambient siliciclastic rocks are remarkably high compared to CAI temperatures (figure 4B). In the eastern part of the Bodón Unit temperatures approach 200 °C for samples from the San Emiliano Fm., while in the central part all KI data indicate temperatures between 200 and 250 °C, both for the San Emiliano Fm. and the pre-Carboniferous ambient rocks. In the central northern part of the Bodón Nappe, KI-derived paleotemperatures for pre-orogenic rock successions clearly increase towards the León Fault (figure 4B). The western part of the Bodón Unit is characterized by slightly lower temperatures, generally around 150 °C, but locally up to 200 °C. The southern part of the Central Coal Basin, north of the León Fault, is characterized by high temperatures between 250 and 300 °C.

VR data for both Stephanian CR and CM basins [START_REF] Frings | Coalification history of the Stephanian Ciñera-Matallana pull-apart basin, NW Spain: Combining anisotropy of vitrinite reflectance and thermal modeling[END_REF][START_REF] Colmenero | Genesis and rank distribution of Upper Carboniferous coal basins in the Cantabrian Mountains, Northern Spain[END_REF]) correspond to temperatures ranging from 120 to 240 °C (figure 5A). These relatively high temperatures are believed to be restricted to the Stephanian pull-apart basins because of their fault-controlled settings and their association with intrusions of sills and dykes [START_REF] Méndez | Estudio de la evolución de los carbones de la Cuenca Ciñera-Matallana[END_REF][START_REF] Frings | Coalification history of the Stephanian Ciñera-Matallana pull-apart basin, NW Spain: Combining anisotropy of vitrinite reflectance and thermal modeling[END_REF]).

FI trapping temperatures, derived from primary Th values (indicated on figure 4A), range from 130 to 150 °C for dolomite cement in Carboniferous rocks (Gasparrini et al. 2006a). Th values for dolomite cement in Cambrian rocks range between 85 and 105 °C (figure 4A) while true trapping temperatures lie approximately 5 to 10 degrees higher [START_REF] Lapponi | Hydrothermal dolomitization in a complex geodynamic setting (Lower Palaeozoic, northern Spain)[END_REF]. Late diagenetic calcite cement is characterized by trapping temperatures between 115 and 130 °C (Gasparrini et al. 2006a). Partial re-equilibration to decrepitation of FIs is described by [START_REF] Ayllón | Re-equilibration of fluid inclusions in diagenetic-anchizonal rocks of the Ciñera-Matallana coal basin (NW Spain)[END_REF] for quartz veins in the Stephanian CM Basin. Here, inclusions were trapped at temperatures between 110 and 120 °C and re-equilibrated by a short-lived thermal event recording temperatures as high as 300 °C. Ore mineralizations from the Bodón Unit are hosted by fracture-related dolomite and record trapping temperatures up to 200 °C [START_REF] Paniagua | Tectonic setting, mineralogical characteristics, geochemical signatures and age dating of a new type of epithermal carbonatehosted, precious metal-five element deposits: the Villamanín area (Cantabrian Zone, northern Spain)[END_REF]), but did not result in reequilibration of the FIs in the dolomite cements.

In one of the Cambrian AFT samples, as well as in the Stephanian AFT sample, [START_REF] Carrière | Neoproterozoic to Holocene tectonothermal evolution of the southern Cantabrian Mountains NW Iberia, revealed by apatite fission-track thermochronology[END_REF] observed partial annealing of fission tracks at about 270 Ma, indicating insignificant duration of peak thermal conditions to completely anneal the fission tracks (figure 4B). The other Cambrian AFT sample indicates complete annealing meaning that temperatures did exceed 110 ± 10 °C for more than 10 Ma. Most samples of the nearby Sabero Coalfield are completely annealed [START_REF] Botor | Thermal history of the Sabero Coalfield (Southern Cantabrian Zone, NW Spain) as revealed by apatite fission track analyses from tonstein horizons: implications for timing of coalification[END_REF].

RE pyrolysis and VR analysis

Tmax values obtained from RE pyrolysis show a very wide range from 449 to 608 °C (table 1), which corresponds to paleotemperatures of 130 to 250 °C (table 1; figure 5B). These data fall in a slightly cooler range of paleotemperatures as those derived from KI data, but are significantly higher compared to paleotemperatures deduced from CAI data in the study area. The 4 samples from the Vegamián Fm. show relatively low Tmax values with less variation (corresponding paleotemperatures: 129 -172 °C), while the many samples from the San Emiliano Fm. and the Lena and Sama Groups show the maximum variability (paleotemperatures: 129 -228 °C). The sample from the Devonian La Vid Group has a very high Tmax of 591 °C (paleotemperature: 236 °C). The maximum Tmax values, with very little variation, are recorded for the oldest sample, the Early Paleozoic Formigoso Fm. (Tmax: 602 °C; paleotemperature: 240 °C) and the youngest samples from the Stephanian CR Basin (Tmax: 603 -608 °C; paleotemperature: 241 -243 °C). Particularly the San Emiliano Fm. (and its equivalents) shows a very wide range of thermal maturity. Locations most intensely affected by fracture-related dolomitization show RE temperatures around 200°C (e.g. the central, western and southeastern parts of the Bodón Unit).

Additionally to Tmax values, organic maturity was studied trough VR analysis on 18 samples. Tmax-derived VR values obtained for these 18 samples (cfr. [START_REF] Jarvie | Oil and Shale Gas from the Barnett Shale, Ft. Worth Basin, Texas[END_REF], show a very wide range (0.92 -2.29 Rm%; table 1). During VR analysis, two populations of kerogen particles were observed in the kerogen slides of all 15 samples from the San Emiliano Fm., as well as in the distribution of measured VR values (figure 6A). The first population consists of less mature organic matter, identified by light to medium brown colors (figure 7A). The second population contains higher mature organic matter, identified by dark brown to black colors (figure 7A). For the current maturity analysis, only the less mature vitrinite population was analyzed. This led to a narrow range of reflectance in all samples from the San Emiliano Fm. (0.61 -0.93 Rm%; table 1), showing a much more homogeneous dataset compared to Tmax and Tmax-derived VR values. In the central part of the study area reflectance values range from 0.79 to 0.93 Rm% (corresponding paleotemperatures: 98 -130 °C). Reflectance values are lower in the eastern part (0.61 -0.72 Rm%; paleotemperatures: 96 -109 °C). At the very western end of the study area, reflectance is high (0.87 Rm%; paleotemperature: 124 °C), as well as in the southeastern end, where reflectance increases to 0.79 Rm% (paleotemperature: 116 °C). The total vitrinite assemblage for five samples from the San Emiliano Fm. has been analyzed to compare with the Tmax-derived VR values (table 1). Reflectance data from the total vitrinite assemblages are close to the VR values calculated from Tmax. Except for one sample, all measured VR values are slightly lower. The difference between calculated and measured VR values seems to increase towards higher maturity.

Samples from the Vegamián Fm. show strongly altered, poorly preserved organic matter (figure 7B). The alteration leads to lower reflectance values compared to the well preserved vitrinite in other samples. In one sample all organic matter, including vitrinite, is strongly altered. Therefore the reflectance value is much lower (0.45 %) compared to the other samples and was excluded from the maturity analysis. In the two other samples (underlined in blue in figure 5A) vitrinite was poorly to moderately preserved. Careful selection of the few relatively well preserved vitrinite particles, and avoiding the altered particles, led to usable VR values, which still might be slightly lower than in well preserved samples.

The paleotemperatures calculated from VR values of the less mature vitrinite population range from 96 to 130 °C, which is significantly lower compared to the temperatures derived from RE pyrolysis and previous VR data of the CM and CR basins (figure 5A; [START_REF] Frings | Coalification history of the Stephanian Ciñera-Matallana pull-apart basin, NW Spain: Combining anisotropy of vitrinite reflectance and thermal modeling[END_REF][START_REF] Colmenero | Genesis and rank distribution of Upper Carboniferous coal basins in the Cantabrian Mountains, Northern Spain[END_REF]. It has to be taken into account, however, that these Stephanian pull-apart basins experienced a different depositional and diagenetic history due to their fault-controlled setting [START_REF] Méndez | Estudio de la evolución de los carbones de la Cuenca Ciñera-Matallana[END_REF][START_REF] Frings | Coalification history of the Stephanian Ciñera-Matallana pull-apart basin, NW Spain: Combining anisotropy of vitrinite reflectance and thermal modeling[END_REF][START_REF] Botor | Thermal history of the Sabero Coalfield (Southern Cantabrian Zone, NW Spain) as revealed by apatite fission track analyses from tonstein horizons: implications for timing of coalification[END_REF].

Discussion

RE pyrolysis versus VR analysis

The less mature population of kerogen particles identified through optical kerogen analysis represents in situ organic matter recording the in situ basin maturity, while the higher mature population represents reworked organic matter showing the thermal overprint of previous depositional systems. Reflectance values of the total vitrinite population in the 5 analyzed samples from the San Emiliano Fm. are always higher compared to values of the corresponding in situ vitrinite population, but the differences show a significant variation (table 1). These differences seem to be controlled mainly by different proportions of highly mature reworked kerogen and less mature in situ kerogen, as observed in the kerogen slides (figure 7A). This can also explain the difference between the less mature in situ VR data and the higher mature Tmax-derived VR data, as Tmax (and Tmaxderived VR) is based on whole rock samples, including both the less mature in situ OM and the higher mature reworked OM. The differences between Tmax-derived VR data and in situ VR data are thus controlled by changes in the proportions of higher mature reworked vitrinite and less mature in situ vitrinite.

The occurrence of highly mature reworked OM in the samples from the San Emiliano Fm. supports the hypothesis of reworked detrital illite influencing the KI data recorded for synorogenic successions in the CZ, which has been proposed in several studies [START_REF] Bastida | Tectono-thermal evolution in a region with thinskinned tectonics : the western nappes in the Cantabrian Zone (Variscan belt of NW Spain)[END_REF][START_REF] Brime | Transition from Diagenesis to Metamorphism Near the Front of the Variscan Regional Metamorphism (Cantabrian Zone, Northwestern Spain)[END_REF][START_REF] Aller | Superposition of tectono-thermal episodes in the southern Cantabrian Zone (foreland thrust and fold belt of the Iberian Variscides, NW Spain)[END_REF]. Indeed, KIderived temperatures for the San Emiliano Fm. (figure 4B) are distinctly higher compared to the temperatures obtained from both in situ (VR) and bulk rock paleothermometers (RE), and are even higher compared to fluid inclusion trapping temperatures in dolomite and quartz gangue minerals. The KI-derived paleotemperatures clearly represent overestimations of the true thermal maturity, most likely due to incorporation of detrital illite. The synorogenic San Emiliano Fm. is composed of siliciclastic sediments filling the foredeep of the Variscan Orogen [START_REF] Brouwer | La succession Carbonifère dans la partie méridionale des montagnes Cantabriques (Espagne Nord-Ouest)[END_REF]. The highly mature detrital illite and the reworked kerogen most likely represent eroded sediments, derived from the deformed and metamorphosed Variscan Orogen to the southwest composed of older pre-orogenic rocks [START_REF] Julivert | Hercynian orogeny and Carboniferous palaeogeography in northwestern Spain: a model of deformation-sedimentation relationships[END_REF][START_REF] Bastida | Tectono-thermal evolution in a region with thinskinned tectonics : the western nappes in the Cantabrian Zone (Variscan belt of NW Spain)[END_REF]. The importance of sediment reworking in pre-orogenic siliciclastic successions cannot be assessed for in this study. Nevertheless, also for these successions, the KI data indicate very high paleotemperatures not in agreement with the real thermal maturity of the sediments. Next to sediment reworking, additional causes for the anomalously high KI-derived temperatures are related to data measurement and processing, and to the empirical equation used for paleotemperature calculation. The KI data of [START_REF] Marschik | Der Übergang von Diagenese zur sehr niedergradigen Metamorphose im externen Variszikum (Kantabrische Zone)[END_REF] were obtained by X-ray diffraction analysis of airdried samples [START_REF] Aller | Superposition of tectono-thermal episodes in the southern Cantabrian Zone (foreland thrust and fold belt of the Iberian Variscides, NW Spain)[END_REF]. This may result in an overestimation of the KI, as the presence of mixed-layer illite/smectite will modify the X-ray diffractograms in air-dried samples [START_REF] Junfeng | Relationschip between illite crystallinity and temperature in active geothermal systems of New Zealand[END_REF][START_REF] Bastida | Tectono-thermal evolution in a region with thinskinned tectonics : the western nappes in the Cantabrian Zone (Variscan belt of NW Spain)[END_REF]. Empirical equations used to calculate KI-derived temperatures are not widespread in literature, and most of these equations are bound by specific preconditions. The equation of [START_REF] Mukoyoshi | Quantitative estimation of temperature conditions for illite crystallinity: comparison to vitrinite reflectance from the Chichibu and Shimanto accretionary complexes, eastern Kyushu, Southwest Japan[END_REF] yields the most realistic paleotemperatures for the range of KI data of [START_REF] Marschik | Der Übergang von Diagenese zur sehr niedergradigen Metamorphose im externen Variszikum (Kantabrische Zone)[END_REF] and [START_REF] Aller | Superposition of tectono-thermal episodes in the southern Cantabrian Zone (foreland thrust and fold belt of the Iberian Variscides, NW Spain)[END_REF].

The incorporation of reworked OM and detrital illite limits the interpretation of RE and KI data of the Carboniferous formations for paleothermometry purposes, since they represent synorogenic successions [START_REF] Marcos | An approach to the tectonostratigraphic evolution of the Cantabrian Foreland thrust and fold belt, Hercynian Cordillera of NW Spain[END_REF]; figure 2). Maturity analyses and paleothermal studies based on bulk rock techniques such as RE pyrolysis and the KI of illite crystallinity have a high risk of failure and erroneous interpretations in most sedimentary systems due to sediment reworking. Ideally, additional data are needed to cross-check the consistency and reliability of these bulk rock techniques. In this study, the presence of two populations of OM with different maturity levels, as observed through optical kerogen analysis, identified the problematic effect of reworking. Next to optical analysis, reworking can also be observed in the pattern of CO2 produced during oxidation, following RE pyrolysis (figure 6B). Although RE pyrolysis can be used for fast screening of thermal maturity levels, additional maturity data are needed to verify the results from RE pyrolysis. This study indicates that a comparison with optical kerogen analysis of vitrinite is needed for more accurate interpretations, especially for the analysis of synorogenic sediments. Through VR analysis, reworked or altered vitrinite particles can be excluded, leading to more reliable and precise information on the thermal maturity of in situ vitrinite.

Comparison with paleothermal data from literature

The highest temperatures recorded for dolomitization and ore mineralization in the Variscan basement rocks of the Bodón Unit are derived from fluid inclusions in dolomite (130 -150 °C: Gasparrini et al. 2006a) and in quartz gangue minerals in the central part of the study area (200 °C;[START_REF] Paniagua | Tectonic setting, mineralogical characteristics, geochemical signatures and age dating of a new type of epithermal carbonatehosted, precious metal-five element deposits: the Villamanín area (Cantabrian Zone, northern Spain)[END_REF]. Trapping temperatures of fluid inclusions indicate the maximum temperatures of hydrothermal fluid flow related to the Post-Variscan thermal event, as they provide a snapshot of the temperature during dolomitization and ore mineralization (e.g. [START_REF] Middleton | Hot fluid flow events in Atlantic margin basins: an example from the Rathlin Basin[END_REF]. Hydrothermal fluid flow must have been most important in the central and western parts of the Bodón Unit, as well as in the southeastern part, given the abundance of dolomite and ore mineral deposits in these areas (figure 1B). The occurrence of high CAI-derived paleotemperatures in carbonate precursors close to dolomitized and mineralized areas (figure 4A) attests of the importance of hydrothermal fluid flow in these parts of the Bodón Unit. The episodes of fluid flow must have been sufficiently long to allow the conodont fragments to equilibrate. Likewise, RE-derived temperatures of the Early-Carboniferous Vegamián Fm. are higher in the central part of the Bodón Unit (151 -172 °C) compared to the eastern part (129 °C).

Paleotemperatures derived from the in situ VR dataset of the Vegamián and San Emiliano fms. are significantly lower compared to published data for the Stephanian CM and CR basins [START_REF] Frings | Coalification history of the Stephanian Ciñera-Matallana pull-apart basin, NW Spain: Combining anisotropy of vitrinite reflectance and thermal modeling[END_REF][START_REF] Colmenero | Genesis and rank distribution of Upper Carboniferous coal basins in the Cantabrian Mountains, Northern Spain[END_REF]; figure 5A), highlighting their different thermal histories. They have been more intensely affected by hydrothermal fluid flow and volcanic intrusions due to their fault-controlled setting [START_REF] Méndez | Estudio de la evolución de los carbones de la Cuenca Ciñera-Matallana[END_REF][START_REF] Frings | Coalification history of the Stephanian Ciñera-Matallana pull-apart basin, NW Spain: Combining anisotropy of vitrinite reflectance and thermal modeling[END_REF], and due to the relatively high porosity and permeability of their sedimentary successions compared to the tight Paleozoic basement rocks making up the Bodón Unit [START_REF] Duddy | Recognition of the thermal effects of fluid flow in sedimentary basins. In: Parnel J (ed) Geofluids: Origin, Migration and Evolution of Fluids in Sedimentary Basins[END_REF]. Althought the RE and KI-derived temperatures have to be treated with caution given the possible influence of reworked sediments, they seem to confirm the higher thermal maturity of the Stephanian basins (figures 4B and 5B).

Temperatures obtained from CAI data for precursor carbonates show baseline values around 60 °C in the eastern part of the study area, which was less affected by fracture-related hydrothermal dolomitization and ore mineralization (figures 1B and 4A). VR data yield temperatures around 100 °C for siliciclastics in the same area. This difference of approximately 40°C can have several causes related to either the alteration of conodonts or the reflectance of vitrinite. The variety of empirical equations to calculate paleotemperatures from VR data yields divergent paleotemperatures (table 2), which results in an error range regarding VR-derived paleotemperatures. The possible suppression of vitrinite reflectance [START_REF] Schito | Overcoming pitfalls of vitrinite reflectance measurements in the assessment of thermal maturity: the case history of the lower Congo basin[END_REF]) cannot be assessed for in this study. Regarding the CAI data, it is known that they, compared to other paleothermometers, might be less sensitive. [START_REF] Zhang | Late Ordovician-Early Silurian conodont biostratigraphy and thermal maturity, Hudson Bay Basin[END_REF] state that CAI values changing from 1 to 1.5 (49 to 62 °C) are not as sensitive as corresponding VR or KI data. Moreover, Raven and van der Pluijm (1986) report a general range of uncertainty of ± 0.5, corresponding to a temperature range of ± 25 °C. Secondly, the applied heating time of 10 Ma might be an overestimation of the true time of heating. This is unlikely, given the AFT results of [START_REF] Carrière | Neoproterozoic to Holocene tectonothermal evolution of the southern Cantabrian Mountains NW Iberia, revealed by apatite fission-track thermochronology[END_REF] and Botor and Anckiewicz (2015) which indicate that following the Post-Variscan thermal event, temperatures remained high for a significant period of time. Even if the heating time is decreased to 1 Ma on the Arrhenius plot of Eppstein et al. (1977), temperatures around 75 °C are obtained, which are still significantly lower compared to VR-derived temperatures. A third cause might be retardation of the color alteration process due to hydrothermal alteration. Hydrothermal alteration is usually accompagnied by indicative textures and patinas [START_REF] Voldman | Conodont palaeothermometry of contact metamorphism in Middle Ordovician rocks from the Precordillera of western Argentina[END_REF]), which have indeed been observed in the CZ [START_REF] Raven | Metamorphic fluids and transtension in the Cantabrian Mountains of northern Spain: an application of the conodont colour alteration index[END_REF]Bastida et al. 1999;Blanco-Ferrera et al. 2016).

Next to a lower sensitivity of the CAI technique, a secondary cause which might explain the observed differences between CAI-and VR-derived paleotemperatures, is the influence of deformation partitioning in different mechanical units of the Variscan FFTB. CAI data are obtained from the thick Cambrian and Carboniferous carbonate successions (figure 3D) representing rigid carbonate units. Hydrothermal fluid circulation was likely restricted to spaced fractures in these tight units, causing fracture-related dolomitization and ore mineralization, and resulting in a strong thermal overprint of nearby carbonate precursors. VR data on the other hand, are obtained from weak siliciclastic clay-rich successions which were more intensely deformed and might have allowed a more wide-spread circulation of hydrothermal fluids and a more effective heat dissipation in respons to the Post-Variscan thermal event.

Hydrothermal versus geothermal dolomitization

In the field of ore geology, hydrothermal activity refers to the action of hot waters directly related to igneous activity, in contrast with geothermal waters which are not related with igneous intrusions [START_REF] Robb | Introduction to Ore-Forming Processes[END_REF]. In dolomitization studies, hydrothermal generally applies to dolomitization by fluids with temperatures at least 10 °C higher than those experienced by the surrounding rocks, including both carbonate precursors and ambient siliciclastics [START_REF] Stearns | Thermal springs in the United States. Contributions to the hydrology of the United States[END_REF][START_REF] White | Thermal waters of volcanic origin[END_REF][START_REF] Machel | Hydrothermal dolomitea product of poor definition and imagination[END_REF]. The current dolomitization model for the dolomite geobodies in the Bodón Unit assumes that thermal convection of hydrothermal marine-derived brines resulted in fracture-related porous dolomites (Gasparrini et al. 2006b). This was based on the difference between FI trapping temperatures of 130 to 150 °C (Gasparrini et al. 2006a) and CAI-derived paleotemperatures of 70 to 95 °C in Carboniferous carbonate formations [START_REF] Raven | Metamorphic fluids and transtension in the Cantabrian Mountains of northern Spain: an application of the conodont colour alteration index[END_REF].

Our study indicates that in the eastern part of the study area, CAI-derived paleotemperatures of undolomitized precursor carbonate rocks average around 60 °C, which is likely an underestimation of their true thermal maturity. Siliciclastic ambient rocks are characterized by VR-derived temperatures around 100 °C. Since dolomitization and ore mineralization were more important in the central and western parts of the study area, one can assume that these lower temperatures are the result of Pre-Variscan burial and high crustal heat flow during the Post-Variscan thermal event. It confirms the hydrothermal nature of dolomitization (130 to 150 °C) and ore mineralizations (up to 200 °C) which are clearly related to positive temperature anomalies in the central and western parts of the study area.

Conclusions

The thermal maturity of Paleozoic rocks in the southern Cantabrian Zone (the Iberian Variscan foreland foldand-thrust belt) has been investigated to better constrain the thermal conditions of Late-Variscan fluid flow which caused widespread dolomitization and ore mineralization. Maximum temperatures experienced by both carbonate precursor and ambient siliciclastic rocks have been assessed based on a compilation of published paleothermal data (CAI, KI, VR, FI) combined with newly acquired RE pyrolysis and VR analysis data.

Both RE and KI data from synorogenic successions yield paleotemperatures that are anomalously high compared to in situ VR data and CAI data. Optical kerogen analysis of VR samples revealed an important amount of reworked highly mature kerogen, giving evidence for significant input of reworked sediments in these synorogenic units. The presence of reworked sediments, most likely eroded material from higher metamorphic zones of the Variscan Orogen, critically affects maturity analysis through bulk rock techniques like RE pyrolysis, KI and VR analysis of total vitrinite populations. This is a general problem in many sedimentary systems, where maturity analyses based on bulk rock data are limited by uncertainties and the risk of erroneous misinterpretations. During optical kerogen analysis, in situ and reworked kerogen populations can be identifed. Using only the reflectance values of in situ populations significantly increases the accuracy and reliability of the maturity data and minimizes the risk of failure of maturity analyses. CAI-derived paleotemperatures for carbonate precursor rocks not affected by intense fluid circulation average around 60 °C and are attributed to burial and increased crustal heat flow following lithospheric delamination. Paleotemperatures within these carbonate rock units increase towards dolomitized and mineralized zones, which formed during Post-Variscan fluid flow at temperatures up to 200 °C [START_REF] Paniagua | Tectonic setting, mineralogical characteristics, geochemical signatures and age dating of a new type of epithermal carbonatehosted, precious metal-five element deposits: the Villamanín area (Cantabrian Zone, northern Spain)[END_REF]. Siliciclastic ambient rocks affected by burial and high crustal heat fluxes have a higher average thermal maturity, corresponding to temperatures around 100°C (based on in situ VR data). The different thermal maturity of carbonate and siliciclastic rocks is likely caused by a difference in analytical techniques used to obtain CAI and VR data and their corresponding paleotemperatures. Conodont fragments were likely not as sensitive to the thermal overprint compared to VR, and the color alteration process might have suffered from hydrothermal alteration. A secondary cause to explain the difference in thermal maturity might be a different response to Variscan compression and orocline formation of rigid precursor carbonates and softer ambient siliciclastics which caused the hot fluids to circulate in a different manner. Circulation was confined along spaced fractures through the carbonates, and was more widespread through the siliciclastics, resulting in more efficient heat transport.

Massive fracture-related dolomitization of Cambrian and Carboniferous carbonates of the southern Cantabrian Zone occurred at temperatures between 90 and 150 °C (based on FIs). Undisturbed ambient siliciclastic rocks experienced temperatures around 100°C, whereas this temperature might have been lower for precursor carbonates. Dolomitization was thus caused by hydrothermal fluids, circulating in thermal disequilibrium with the precursor carbonate rocks. The positive temperature anomalies recorded in the precursor carbonate rocks by CAI data broadly correlate with the occurrence of dolomitization and ore mineralization. This implies that in similar settings, mapping the thermal maturity anomalies within carbonate precursors could be used as an exploration tool for porous hydrothermal dolomitic geobodies and ore mineralization.

FIGURE CAPTIONS

Fig. 1 (A) Geological map of the CZ with indication of the main fault systems and tectonostratigraphic units (modified after [START_REF] Dallmeyer | Pre-Mesozoic Geology of Iberia[END_REF]. Unconformable Stephanian and Mesozoic-Cenozoic successions have been indicated. The Bodón Unit is encircled in red. The León Fault is indicated with black arrows. (B) Geological map of the Bodón Unit modified from [START_REF] Lobato | Hoja del Mapa Geológico Nacional de España[END_REF], [START_REF] Alonso | Hoja del Mapa Geológico Nacional de España (1:50000)[END_REF][START_REF] Rodríguez | Hoja del Mapa Geológico Nacional de España (1:50000), Memoria 102 (Los Barrios de Luna). IGME, Madrid Suggate RP[END_REF] and [START_REF] Alonso | Paleogeographic inversion resulting from large out of sequence breaching thrusts: The León Fault (Cantabrian Zone, NW Iberia). A new picture of the external Variscan Thrust Belt in the Ibero-Armorican Arc[END_REF]. The Gayo, Bodón and Forcada nappes are indicated, together with Carboniferous carbonate outcrops and the location of dolomite geobodies (based on Muñoz Quijano 2015) and ore mineralizations. (Gasparrini et al. 2006a;[START_REF] Lapponi | Hydrothermal dolomitization in a complex geodynamic setting (Lower Palaeozoic, northern Spain)[END_REF]) and (B) paleotemperatures based on KI data [START_REF] Marschik | Der Übergang von Diagenese zur sehr niedergradigen Metamorphose im externen Variszikum (Kantabrische Zone)[END_REF][START_REF] Aller | Superposition of tectono-thermal episodes in the southern Cantabrian Zone (foreland thrust and fold belt of the Iberian Variscides, NW Spain)[END_REF]) and locations of AFT thermochronology [START_REF] Carrière | Neoproterozoic to Holocene tectonothermal evolution of the southern Cantabrian Mountains NW Iberia, revealed by apatite fission-track thermochronology[END_REF].

Fig. 5

Structural sketch map of the Bodón Unit with indication of temperatures based on (A) VR analysis [START_REF] Frings | Coalification history of the Stephanian Ciñera-Matallana pull-apart basin, NW Spain: Combining anisotropy of vitrinite reflectance and thermal modeling[END_REF][START_REF] Colmenero | Genesis and rank distribution of Upper Carboniferous coal basins in the Cantabrian Mountains, Northern Spain[END_REF] this study) and (B) RE pyrolysis (this study). [START_REF] Jarvie | Oil and Shale Gas from the Barnett Shale, Ft. Worth Basin, Texas[END_REF]) and corresponding temperatures (cfr. [START_REF] Barker | Calculation of vitrinite reflectance from thermal histories and peak temperatures. A comparison of methods[END_REF]. HI and OI roughly correspond to the H/C and O/C atomic ratios for each sample. VR results are expressed in Rm%. The first column of VR data refers to the total vitrinite assemblage. The second and third columns of VR data represent in situ kerogen and corresponding temperatures (cfr. [START_REF] Barker | Calculation of vitrinite reflectance from thermal histories and peak temperatures. A comparison of methods[END_REF]. authors cannot conclude if their measurements have been performed on real and/or larger vitrinite grains, so we prefer not to specifically mention this in the manuscript.

 You should absolutely mention that the occurrence of dolomite geobodies as indicated on the figures (i.e. figures 1B, 4 & 5) is based on the work of Muñoz-Quijano (2015).

The occurrence of dolomite geobodies indicated on the figures is indeed based on the work of [START_REF] Muñoz-Quijano | Hydrothermal dolomitization of Paleozoic successions in Northern Spain: petrophysical properties and structural control[END_REF]. This is mentioned in the manuscript (lines 167-169) and we have now also added a citation in the caption of figure 1 (lines 846-848).

 Several minor grammatical errors were highlighted by Thilo Bechstaedt. These errors (on lines 138, 197, 330, 473 and 595) have been corrected in the revised manuscript.

Fig. 2

 2 Fig. 2 Stratigraphic column of the Bodón Unit (modified from Aller et al. (2005) and Pérez-Estaún et al. (1988)) showing Carboniferous precursor carbonate successions in light blue and dolomite geobodies in pink. Stephanian siliciclastics are indicated in purple. Pre-orogenic and synorogenic successions are distinguished, as well as what is here referred to as dolomite precursor carbonate rocks and ambient siliciclastic rocks. The stratigraphic intervals sampled for the different published paleothermometry datasets and sampled for this study are shown on the right.

Fig. 3

 3 Fig. 3 (A) Exposure of Upper Devonian and Lower Carboniferous successions in the Las Chábanos quarry in the northern central part of the Bodón Unit. Formation boundaries are indicated with white dotted lines. Siliciclastic rocks of the Devonian La Vid Group are exposed on the far left, followed by brownish sandstone beds of the Ermita Fm. The black succession in the centre of the picture represents the Lower Carboniferous Vegamián Fm. which is overlain by nodular limestones of the Alba Fm. On the far right side, the central chert member of the Alba Fm. can be seen. Note professor for scale. (B) Diachronic contact between Valdeteja limestones (left) and San Emiliano clay-and sandstones (right) in the central part of the Bodón Unit. Note utility pole on top for scale. (C) Typical exposure of San Emiliano Fm. claystones, characterized by the occurrence of diffuse tectonic cleavage. Note hammer for scale. (D) Example of dolomitized limestones of the Valdeteja Fm. in the southeastern part of the Bodón Unit. Light gray rocks are the precursor limestones, whereas dark gray rocks are dolomites. Note utility poles on lower left side for scale.

Fig. 4

 4 Fig. 4 Structural sketch map of the Bodón Unit with indication of (A) paleotemperatures based on conodont CAI data (Raven and van der Pluijm 1986; Aller at al. 2005) and FI homogenization temperatures(Gasparrini et al. 2006a;[START_REF] Lapponi | Hydrothermal dolomitization in a complex geodynamic setting (Lower Palaeozoic, northern Spain)[END_REF]) and (B) paleotemperatures based on KI data[START_REF] Marschik | Der Übergang von Diagenese zur sehr niedergradigen Metamorphose im externen Variszikum (Kantabrische Zone)[END_REF][START_REF] Aller | Superposition of tectono-thermal episodes in the southern Cantabrian Zone (foreland thrust and fold belt of the Iberian Variscides, NW Spain)[END_REF]) and locations of AFT thermochronology[START_REF] Carrière | Neoproterozoic to Holocene tectonothermal evolution of the southern Cantabrian Mountains NW Iberia, revealed by apatite fission-track thermochronology[END_REF]).

Fig. 6

 6 Fig. 6(A) Histogram illustrating VR data expressed as Rm(%) for sample 46, showing a large population of in situ material and two groups of reworked higher mature material. (B) CO2 oxidation pattern from RE analysis for sample 46, showing the occurrence of two types of OM (black arrows). The third peak to the right corresponds to the carbonate content of the sample.

Fig. 7

 7 Fig. 7 Kerogen photomicrographs. (A) Two populations of kerogen can be recognized in this samples from the San Emiliano Fm. Light brown colors (white arrow) characterize the least mature population, interpreted as in situ, whereas dark brown to black colors (black arrow) characterize the reworked most mature population. (B) Strongly altered OM in a sample from the Vegamián Fm.

TABLE CAPTIONS Table 1

 CAPTIONS1 Results from RE pyrolysis and VR analysis. RE Tmax values are indicated, as well as hydrogen and oxygen indices (HI and OI), TOC, Tmax-derived VR values (cfr.

Table 2

 2 Paleotemperatures in °C calculated with common equations for increasing vitrinite reflectance levels. The difference in paleotemperatures obtained from different equations increases with increasing reflectance level.

	FIGURE 1

Acknowledgements

This study was funded by Research Foundation -Flanders (FWO). Additional funding was also provided by the IAS 2014 Grant Scheme. Rock-Eval measurements were performed at IFPEN (France). Daniel Pillot and Géremie Letort are kindly acknowledged for assistence with the analyses. We thank Elisabeth Bemer (IFPEN), head of the "water-rock interaction" project, who provided additional funding for vitrinite reflectance analyses.

This document provides a detailed overview of changes applied to the revised manuscript in response to the remarks of reviewer #1 (Prof. Emer. Thilo Bechstädt).

Reviewer #1: Prof. Emer. Thilo Bechstädt:

 The vitrinite reflectance measurements performed in this study are not classical measurements, but a method specifically applied at the GeoResources SCT lab in Heidelberg. Can you provide citations to publications of this specific method or mention this more clearly?

To the knowledge of the authors, there are currently no publications on this specific method of VR analysis performed at the Georesources SCT lab, but we adapted the methodological section of the manuscript by specifying that these measurements are measurements of dispersed kerogen (line 339), different from classical VR analysis. We also added that the measurements have been performed at the GeoResources SCT lab (line 341). The method itself is extensively described in the methodology (lines 339-357).

 During VR analysis, two populations of kerogen particles were observed in the kerogen slides. The first population is less mature and interpreted as "in situ", while the second population is more mature and interpreted as "reworked". This is correct but it is an interpretation and should be discussed in the discussion.

We fully agree with Thilo Bechstädt. The "in situ vs. reworked" hypothesis should be introduced in the discussion section, since it is an interpretation and not a result. We have now modified the results section where we only refer to a "less mature" and a "more mature" kerogen population. In de discussion section, we discuss the interpretation of "in situ vs. reworked" material (lines 443-445).

 The study indicates that a comparison with data from vitrinite reflectance analysis is needed to verify the results from Rock-Eval pyrolysis. You should mention that it is not classical VR analysis but optical kerogen analysis of vitrinite reflectance, in order to avoid confusion. This has been adapted (lines 484-486).

 Older studies reporting data on vitrinite reflectance analysis have measured real, larger vitrinite grains in the coal basins, so an input of older material can be excluded in these studies.

The older studies referred to (i.e. [START_REF] Frings | Coalification history of the Stephanian Ciñera-Matallana pull-apart basin, NW Spain: Combining anisotropy of vitrinite reflectance and thermal modeling[END_REF][START_REF] Colmenero | Genesis and rank distribution of Upper Carboniferous coal basins in the Cantabrian Mountains, Northern Spain[END_REF]) performed VR analysis on samples from Stephanian coal basins. This is mentioned several times in the manuscript (e.g. lines 280-281 & 378-379). However, based on these publications, the