Geology of the Cabuérniga Fault System: evolution of a large Alpine structure with Variscan inheritance

ABSTRACT The Cabuérniga Fault System (CFS) extends for 150 km in E–W trend along the northern divide of the Cantabrian Mountains, near the Cantabrian Sea. This fault system was reactivated several times since the Paleozoic, and displays hydrothermal activity nowadays. We present here a detailed geological map of the whole CFS and six cross-sections distributed along strike, which allowed us to unravel the kinematic history of the individual fault segments comprising the CFS. The current CFS geometry mainly results from the reactivation of E–W Variscan thrusts and NW–SE Late-Variscan strike-slip faults during the Alpine Cycle. Reactivation took place in relation to both Permian and Mesozoic pre-orogenic extension and subsequent Cenozoic shortening. The eastern end of the CFS corresponds to a N–S structure arriving at the Cantabrian coast, while the western edge connects with the E–W trending Llanera fault, leading to an Alpine fault system extending 220 km.


Introduction
The Cabuérniga Fault System (CFS) is located in the Cantabrian Mountains, an E-W mountain range of Northern Iberia with elevations over 2600 m (see sketches 1-4 at the bottom of the geological map, Figure 1) that are the western extension of the Pyrenees. The CFS also displays a rough E-W direction, extending along 150 km between the localities of Arriondas (Asturias) and Ramales de la Victoria (Cantabria). East from Ramales de la Victoria, this fault system ends in a N-S fault that reaches the sea near the town of Santoña (Cantabria) (Sketch 2). To a great extent, the CFS conforms to the Cabuérniga fault as defined by Carreras Suárez et al. (1978). However, rather than a single structure, there is a complex system of anastomosing faults that we have named 'Cabuérniga Fault System'. From a geological point of view, the CFS belongs to the Alpine orogenic belt that gave rise on its westernmost termination to the Pyrenean-Cantabrian Mountains (Alonso et al., 1996;Martín-González & Heredia, 2011b;Muñoz, 2002Muñoz, , 2019and references therein). The northern and southern borders of the Cantabrian Mountains correspond to the Cantabrian Sea (Bay of Biscay) and the Cenozoic Duero foreland basin (Sketches 1 and 2). Phanerozoic rocks deformed by compressive and extensional structures developed during the Variscan and Alpine orogenic cycles crop out along this mountain belt.
The interest of the CFS lies on its large extension and long-lasting activity. The system contains inherited Variscan structures, and was active throughout the Alpine Cycle until recent times. Its complexity is evidenced by a wide range of interpretations on its kinematics, including: (i) south-dipping Mesozoic normal fault, which underwent Cenozoic inversion (García-Mondéjar et al., 1986), (ii) south-dipping Mesozoic normal fault cut by a north-dipping Cenozoic reverse fault (García-Espina, 1997) or (iii) vertical Mesozoic normal fault reactivated as strike-slip fault in Cenozoic times (Tavani, 2012).
In order to study the complex evolutionary history of this fault system and improve the construction of a coherent kinematic model, a detailed geological map has been carried out, distinguishing: (i) Variscan compressional structures (to improve the legibility of the map, these have been simplified), (ii) Late-Variscan faults and (iii) extensional-compressional structures from the Alpine cycle, indicating the age of the identified reactivations. The map aims the characterisation of the CFS along its complete extension along the Cantabrian Mountains, with special emphasis on the Paleozoic basement, where this fault system has been scarcely studied in previous works. Such characterisation includes also the control exerted by Alpine faults on the deposition of some lithostratigraphic units, the distribution of mineral deposits and the hydrothermal water springs location.

Geological setting
The stratigraphic succession of the studied area includes sedimentary rocks from the Cambrian to the Cenozoic, (see the chronostratigraphic scheme on the map). Martín-González and Heredia (2011a, 2011b) established a subdivision of the Cantabrian Mountains into two regions: (a) the Vasco-Cantábrica Region to the east, where Mesozoic rocks dominate, and (b) the Astur-Galaica Region to the west, where Paleozoic basement crop out extensively (Sketch 1). The CFS is located in the northern sector of both regions (Sketch 1). In the Vasco-Cantábrica Region, Alpine structures deform the Mesozoic-Cenozoic cover and are easier to recognise than across the Paleozoic outcrops of the Astur-Galaica Region. In the study area, the Variscan structures are Carboniferous in age (Marcos & Pulgar, 1982;Merino-Tomé et al., 2009;Pulgar et al., 1999) and belong to the Cantabrian Zone of the Variscan Orogen (Julivert, 1971). This zone represents the fold and thrust belt of the Variscan Orogen in the northwest of Iberian Peninsula and exhibits an arcuate shape convex to the west: the Ibero-Armorican Arc (see Weil, 2006 and references therein) (Sketch 2).
Between the Variscan and Alpine orogenic events, the whole Pyrenean-Cantabrian realm experienced strong lithospheric stretching in the Mesozoic, giving rise to the Bay of Biscay opening and the formation of the North-Iberian Margin (e.g. Pedreira et al., 2015 and references therein). During this time, several hyperextended rift basins were developed between the Iberian subplate and the European plate (Roca et al., 2011;Teixell et al., 2018;Tugend et al., 2014;Vergés & García-Senz, 2001). One of the most subsiding extensional basins of this domain was the Basque-Cantabrian Basin (e.g. Ábalos, 2016;Miró et al., 2020;Quintana et al., 2015;Rat, 1988), later incorporated to the orogenic building (Vasco-Cantabrica Region, Martín-González & Heredia, 2011a, 2011b. The CFS had a relevant role in this extensional period at the beginning of the Alpine cycle and was reactivated later, during the compressional deformation. The CFS affects both Paleozoic rocks, already deformed during the Variscan Cycle (lower Cambrian-lowest Permian) and those of the Permian-Mesozoic and Cenozoic (lower Permian-upper Eocene). The Paleozoic rocks in the study area belong to three tectonostratigraphic units of the Cantabrian Zone of the Variscan Orogen (Sketch 2). These units differ on their Carboniferous synorogenic successions (Colmenero et al., 2002) and are the following: (I) The Bodón-Ponga Unit (Alonso et al., 2009), which crops out in the NW sector of the study area and exhibits a Carboniferous succession (up to 3500 m thick) composed of interbedded siliciclastic and carbonate sedimentary rocks, the latter prevailing northwards. (II) The Picos de Europa Unit (Julivert, 1971), situated south and east from the former unit, covering most of the map and exhibiting a Carboniferous succession (up to 2100 m thick) eminently formed by carbonate rocks. Even though Alonso et al. (2009) included this unit within the Bodón-Ponga Unit, it is considered separately in this work, since its abundance of massive limestones also results in a very distinctive morphology in the Cantabrian Mountains.
(III) The Pisuerga-Carrión Unit (Julivert, 1971;Pérez-Estaún et al., 1988), situated south from the Picos de Europa Unit, it is present only on a small sector of the map, where its Carboniferous succession (up to 2000 m thick) is eminently siliciclastic. The Variscan pre-orogenic succession (mainly middle Cambrian-Early Ordovician) displays a reduced thickness (∼1000 m max.) in the Bodón-Ponga Unit, and it is absent in the southern part of the Picos de Europa Unit and in the Pisuerga-Carrión Unit (Aramburu & García-Ramos, 1993). This succession contains several formations that are grouped in the geological map and chronostratigraphic scheme (CD). The Variscan pre-and synorogenic successions were deposited in marine basins, related to the NW passive margin of Gondwana and a peripheral foreland basin, respectively.
The rocks related to the Alpine Cycle can be also grouped in a pre-orogenic and a synorogenic succession. The continental/marine pre-orogenic succession was deposited in the North-Iberian passive margin, is very thick (up to 8000 m) and has a lower Permian-Paleocene age. The Alpine synorogenic succession in the study area is very thin (less than 200 m) and was deposited in a continental foreland basin, although it can reach a thickness of 4000 m more to the south, in the Duero Basin (Sketches 1 and 2). The oldest synorogenic rocks date back to the upper Eocene, in contrast to the Pyrenees where the synorogenic sedimentation start in the later Upper Cretaceous (Muñoz, 2002). The Alpine pre-orogenic sedimentation is largely siliciclastic and mainly continental in its basal part (the Permian-Triassic succession), marine in the Lower and Middle Jurassic (extensive carbonate platform deposits), marine to continental siliciclastic in the Upper Jurassic-Lower Cretaceous and mainly calcareous and marine again in Upper Cretaceous times, although at this age the carbonate platform is more discontinuous (Ábalos, 2016;García-Senz et al., 2019;Heredia et al., 1990;Miró et al., 2020 and references therein).

Methods
A geological map at scale 1:25.000 was elaborated, based in previous data and new field data projected on ortophotographies and topographic maps from the Instituto Geográfico Nacional (IGN). This map was subsequently simplified to the broader scale 1:100.000 for publication, preserving the most outstanding cartographic information of the CFS. The map integrates previous data from: (i) published geological maps edited by the Spanish Geological Survey -Instituto Geológico y Minero de España (IGME), including 1:50.000 scale maps from the MAGNA 'Mapa Geológico Nacional/National Geological Map' (Carreras Suárez et al., 1974;Carreras Suárez et al., 1978;Julivert et al., 1979;Martínez-García et al., 1984;Portero García et al., 1975a, 1975b (Heredia et al., 1990) and several sheets of the more recent 'Mapa Geológico de Cantabria' at 1:25.000 scale, elaborated by the IGME in collaboration with the Regional Government of Cantabria (García-Senz & Cañas Fernández, 2013a, 2013bGarcía-Senz & Merino Tomé, 2011;Huerta Carmona, 2009;Larrondo Echevarria et al., 2008;Solé Pont et al., 2008a, 2008b, 2008c; (iii) geological maps published in scientific journals (Álvarez-Marrón et al., 1989;Marquínez, 1989); including also (iv) the location of the main mineral deposits and hydrothermal water springs for the study area, obtained from the BDMIN and BDAGUAS databases of the IGME, correcting in some cases these locations with new GPS coordinates. The new geological data were obtained during several fieldwork campaigns focused on the resolution and refinement of the cartographic expression of the CFS. All this information has been digitalised using ArcGis. A specific symbology has been composed, using different colours to describe different fault types and activity periods, and highlighting the Alpine structures, which are represented in red. Finally, for a detailed geometric and structural characterisation of the CFS, six N-S geological cross-sections were produced along the fault system and reproduced on top of the geological map.

Description of the structures
The structures of the CFS have been differentiated according to their age, kinematics and orientation. Reactivations of the structures are frequent and occurred in different times, along the entirety of their trace or in segments of favourable orientation. The geometric characteristics and the evolution of the SFC structures over time that are deduced from the geological map are described next.

Variscan structures
The main Variscan structures are thrusts and related folds (Julivert, 1971;Pérez-Estaún et al., 1988) that resulted in a shortening of ∼75% (Merino-Tomé et al., 2009). This tight imbricate system show variable orientations due to its location near the core of the Ibero-Armorican Arc (Sketch 2), and display differing tectonic transport directions and emplacement age. Accordingly, in the Bodón-Ponga Unit, the thrusts were emplaced in the upper Moscovian-Kasimovian (Heredia, 1998;Rodríguez-Fernández & Heredia, 1987), exhibit a general strike direction N-S that changes to E-W northwards, and record a tectonic transport direction toward the east (Weil, 2006). In the Picos de Europa Unit, the thrusts are Kasimovian-Gzhelian in age (Merino-Tomé et al., 2009), display strikes ranging ENE-WSW in their western portion to E-W in their eastern part, and present a generalised tectonic transport direction toward the S-SSE (Farias & Heredia, 1994). The latter thrusts accommodated to a great extent the closure of the Ibero-Armorican Arc by the end of the Variscan Orogeny. They also caused the rotation of the thrusts in the Bodón-Ponga Unit from a N-S direction to an E-W orientation in the northern sector of this unit (Julivert & Arboleya, 1986;Merino-Tomé et al., 2009;Pérez-Estaún et al., 1988;Rodríguez-Fernández & Heredia, 1988; among others) (Sketch 2). Finally, NW-SE strike-slip faults dominate in the fraction of the Pisuerga-Carrión Unit represented in this map.
The Picos de Europa Unit consists in an imbricate thrust system developed in a forward breaking sequence (Farias & Heredia, 1994;Merino-Tomé et al., 2009). In this system, the dominant E-W-trending thrusts represent frontal structures, while those with orientations near N-S (NNE-SSW to NNW-SSE) represent lateral and oblique structures. As a consequence of the stacking of thrusts with listric geometry, these thrusts are nearly vertical to overturned in the northern part of the Unit and dip moderately in the south, as revealed by their sinuous map traces (see geological map and cross-sections I-I' and II-II') (Merino-Tomé et al., 2009).
In the Bodón-Ponga Unit, the thrusts become subvertical to slightly overturned and present a N-S strike in the proximity of the western contact with the Picos de Europa Unit. This arrangement allowed their reactivation as lateral ramps during the southward emplacement of the Picos de Europa Unit. Evidencing this reactivation, unconformable Kasimovian deposits (CSs), syn-tectonic to the emplacement of the Picos de Europa Unit, are present in the hanging wall of the overturned thrusts of the Bodón-Ponga Unit.

Late-Variscan faults
Anastomosing strike-slip faults with a NW-SE to NNW-SSE direction cut through the Variscan thrusts accommodating displacements with a main dextral component. Some of these faults reactivated Variscan thrusts, resulting in nearly vertical reverse faults that still present a strike-slip component. These faults, which are sub-vertical, either strongly dipping toward the NE or, more rarely, toward the SW, were formed prior to the Permian deposits and hence they are considered as Late-Variscan faults. Moreover, the post-Variscan Permian deposits contain volcanic intervals in its lower part, whose plutonic counterparts are as old as ∼297 Ma (Gallastegui et al., 2004). Taking the aforementioned data into account, it is possible to infer that the development of the Late-Variscan faults comprised a short time span (4-5 Myr) between the Gzhelian (age of the youngest Variscan thrusts), and the lower Asselian (age of the oldest Permian deposits). This age is coincident with that estimated for the end of the Ibero-Armorican Arc closure, suggesting that these faults accommodated the last Variscan compressional stresses (López-Gómez et al., 2019;Rodríguez-Fernández & Heredia, 1987). The attitude of the Late-Variscan faults in the study area resembles that of the Ventaniella fault, a large Late-Variscan structure (Late-Hercynian according to Arthaud & Matte, 1977) that crosscuts the Cantabrian Zone west from the geological map (Sketch 2).

Alpine structures
We include under the definition of 'Alpine structures' those related to the whole Alpine cycle, including the normal and transfer faults formed during the Permian-Mesozoic extensional stages and the high-angle reverse faults, thrusts, folds and lateral structures linked to the development of the Cantabrian Mountains during the Cenozoic. Among the Alpine structures, we therefore, include those related to salt tectonics during the Mesozoic and their reactivation during the Alpine compression (e.g. Cámara, 2020; García-Senz et al., 2019 and references therein), represented in the study area by the El Mazo, Cabezón de la Sal, Sarón and Selaya diapirs.

Pre-orogenic extensional structures
They can be recognised because shortening achieved during the Alpine compressional event was moderate at upper crustal levels in this area, with maximum values close to 15% (García-Espina, 1997), although this shortening increases to the north and in depth (Gallastegui, 2000;Pedreira et al., 2015). The first extensional structures immediately postdated the Variscan Orogeny and formed in response to the final extensional collapse of the belt. These structures resulted from reactivation of Variscan and Late-Variscan faults and controlled the development of the small and isolated lower Permian (Cisuralian) basins extending along E-W and NW-SE preferred directions (López-Gómez et al., 2019). Sometimes, the faults controlling the formation of these basins are covered by Triassic deposits, as can be seen in the Deva river gorge, west of the La Hermida village (location in the central part of the map, close to 370.000 UTM line).
After a gap of 30 Myr, a second extensional event is identified in the Middle-Upper Triassic sedimentary record (López-Gómez et al., 2019). This event is related to the incipient rifting that led to the opening of the Bay of Biscay and the separation of the Iberian microplate from the Eurasian plate (Boillot et al., 1979;Rat, 1988). Usually, the faults related to this period were sealed by Jurassic marine deposits (JM1 on the map), although the base of the Jurassic was occasionally affected. These faults are sometimes spatially associated with subvolcanic basic rocks and controlled the deposition and thickness of the Triassic facies. They show different orientations: (i) E-W normal faults, such as the one controlling the thickness of the 'Buntsandstein Facies' (Tb) at the Peña Sagra range (377500-4782500 UTMs); (ii) NW-SE and E-W structures, such as the Aloños and Selaya faults, respectively (ALF and SEF in Sketch 3), which could have controlled the sedimentation of the 'Keuper Facies' (Tk) between the localities of Selaya (Pisueña river valley) and Puente Viesgo (Pas river valley), being the outcrops of this unit absent west and south from these faults. Eastwards, this last extensional period originated various hyperextended rift basins, which were situated between the Iberian microplate and the Eurasian plate (Roca et al., 2011;Teixell et al., 2018;Tugend et al., 2014;Vergés & García-Senz, 2001 and references therein) and accumulated over 10,000 m of sediments (García-Mondejar, 1989).
Overall, the strikes of the Permian-Mesozoic normal faults are commonly the same as those of the main Variscan structures, evidencing their frequent reactivation. In some cases, overturned Variscan thrusts were also reactivated as south-dipping normal faults this might be the case of the El Saigu thrust between the towns of Tresviso and Las Caldas de Besaya (see cross-sections III-III' and IV-IV'). On the other hand, Variscan thrusts tilted to dip angles larger than 40-50°commonly were reactivated as north-dipping Permian-Mesozoic normal faults, which are the most abundant.

Contractional structures
The structures related to the Alpine compressional deformation are thrusts and reverse faults, generally dipping to the north, and open folds that are clearly recognised in the Mesozoic cover but are scarce or absent in the more rigid Paleozoic basement (Alonso et al., 1996;Pulgar et al., 1999).
In the Paleozoic basement of the western area of the map, it is possible to see how the CFS thrusts developed mainly through the reactivation of previous faults, specifically, Variscan thrusts and/or Alpine normal faults with strike directions close to E-W and dipping toward the north. On the contrary, where the pre-existing structures are vertical or south dipping, they generally appear crosscut by the Alpine thrusts, which developed their own trajectory. Nonetheless, overturned Variscan thrusts occasionally reactivated as backthrusts, a process favoured where they previously reactivated as south-dipping normal faults. This appears to be the case for the Variscan thrusts located north from the El Saigu thrust, around Puentenansa, west from cross-section III-III'. Along the eastern and westernmost sectors of the map, the CFS develops imbricate thrust systems and high-angle reverse faults south from the main trace of the Cabuérniga fault (CAF in Sketch 3).
Frequently, the Late-Variscan faults reactivated as oblique and/or lateral ramps during the Alpine contraction, with an associated uplift of their northeastern block. Accordingly, the number of Alpine thrusts occasionally differs from one side of the Late-Variscan faults to the other (e.g. Cabañes and Mier faults; see geological map and Sketch 3).
Several E-W major (kilometre-scale) folds deform the rocks of the Mesozoic Alpine cover, thus being readily recognisable (see cross-sections III-III' to VI-VI'). From north to south, the following folds are present, exhibiting reduced lateral continuity and complex geometries: (I) The Panes synclinorium, north of the El Soplao fault (SOF in Sketch 3), has Cretaceous rocks in its core and is limited to the east by the N-S El Mazo fault (MAF in Sketch 3). East from this fault, the Panes synclinorium is replaced by the Cabezón de la Sal anticlinorium, truncated by a diapir of the same name. Moreover, the El Mazo fault, a normal Triassic fault reactivated in a short segment by the diapir of same name, separates two areas with a different number of thrusts, evidencing its role as a lateral structure during the Alpine shortening. (II) The Cabuérniga anticline, situated south of the El Soplao fault with Triassic rocks in its core, originated in the frontal ramp of the reverse Cabuérniga fault (CAF in Sketch 3). Eastwards, this anticline leads to a series of folds of smaller amplitude corresponding to the eastern termination of the main thrust-fault imbricate of the CFS.
(III) The Tudanca synclinorium is a kilometric-scale structure located south from the Cabuérniga fault. This fold appears displaced southwards, east from the Buelna paleohigh. This shift is linked to the development of the eastern imbricate thrust/fault system, which almost reaches the southern edge of the map (Sketch 3). Furthermore, this imbricate thrust system also resulted in both the southward displacement of the Cabuérniga anticline and the development of other structures such as the Selaya synclinorium, the Arredondo anticlinorium and the La Sía synclinorium (see cross-section VI-VI').
Between the Besaya and Pisueña rivers, the main N-S pre-orogenic transfer faults became reactivated as lateral structures also in the course of the Alpine compressional stage. For instance, the Cotillo fault (COF) constitutes the western termination of the eastern imbricate of the CFS. Similarly, the Besaya fault (BEF on the map) represents the eastern end of the Puentenansa backthrust (Sketch 3). At the terminations of the CFS, multiple NW-SE to E-W, northdipping Alpine pre-orogenic normal faults became reactivated as frontal and oblique ramps during the Alpine compression. Many of those extensional structures, in turn, may represent reactivations of Variscan or Late-Variscan structures (García-Espina, 1997;López-Gómez et al., 2019). Such is the case of the Tánago fault (TAF in Sketch 3), which delimits the Puentenansa backthrust in the west and may constitute the southeastern continuation of the Late-Variscan Mier fault (MIF in Sketch 3). East from Ramales de la Victoria, the CFS ends in another N-S transfer fault known as Ramales fault (RAF in Sketch 3).
In the western end, the CFS connects with the Llanera Thrust System (LLF in Sketch 2), which affected Cenozoic deposits formed in the Oviedo-Infiesto syn-orogenic basin in upper Eocene-Oligocene times, during the main pulse of orogenic uplift (Alonso et al., 1996;Fillon et al., 2016). A small outcrop of these Cenozoic synorogenic deposits appears in the northern part of the geological map. As a whole, the LLF and the CFS belong to a unique E-W Alpine imbricate thrust system extending along 220 km.
Associated with alpine faults of the CFS, six hydrothermal springs exist within the study area (Baeza et al., 2001). From west to east, they are Fuente Obar, La Hermida, La Brezosa, Las Caldas de Besaya, Puente Viesgo and Alceda (see Main map). The hydrothermal water temperature varies between 65°i n La Hermida and 27°in Alceda, then including both hypothermal and hyperthermal springs. The Alceda spring is located within the frontal portion of the eastern imbricate thrust/fault system linked to the CFS. Fuente Obar is coincident with the Late-Variscan San Carlos fault (SCF), while the La Hermida spring lies on the El Saigu thrust, both reactivated in the Alpine cycle. The Caldas de Besaya and Puente Viesgo upwellings are located on Variscan structures reactivated and connected to the CAF.

Fault-related mineralisations
Various mineral associations related with faults are present in the study area, including metallic and non-metallic elements (see the legend of the geological map). Mn and Fe ores usually show a sedimentary control, being located at brecciated zones and/or lateritic levels that result from subaerial exposure of carbonate rocks from either the Paleozoic (CUc, CLc in the map), Jurassic (JM1) or Cretaceous (Cru). The location of the Hg mineralisations tends to be associated with Variscan thrusts or lie in their proximity, suggesting a genetic link with these structures. On the contrary, Pb-Zn and F (associated to minor Cu) ore prospects seem mostly related to Alpine structures and/or Late-Variscan faults, either reactivated during the Alpine Cycle or not. However, these prospects are also present in areas relatively far from faults (at least on the surface), either hosted by Paleozoic limestones or Mesozoic units (specially the Urgonian Complex: units UC1 to UC4 in the map). The age of the F ores generally is Upper Triassic (ca. 206 Ma) (Symons et al., 2017), while Pb-Zn deposits mostly date back to the Lower Cretaceous (ca. 112 Ma) (Symons et al., 2015), thus suggesting that their origin was linked to the end of the main Alpine extensional stages recorded in the study area. Besides, the intense dolomitisation found in the limestones of the Picos de Europa Unit could have taken place mostly in Permian times, as in other areas of the Cantabrian Mountains (Gasparrini et al., 2006a(Gasparrini et al., , 2006b, and is linked to Variscan thrusts and Late-Variscan faults, either non reactivated or reactivated during the extensional collapse of the Variscan Orogen.

Conclusions
The detailed geological map presented here shows a new and clarifying image of the hitherto poorly known Alpine Cabuérniga Fault System (CFS). The faults associated with the CFS have been mapped over 150 km, through the Variscan Basement and the Alpine cover of northern Cantabrian Mountains, indicating the age and main type of movement for each fault or fault segment (normal, reverse or strike-slip). This map has made it possible to relate the genesis of the structures associated with the CFS with the sedimentation of certain lithostratigraphic units in the studied area, as well as to understand the complex kinematic evolution of these faults, identifying their successive reactivations. The role that each fault or fault segment has played over time provides the necessary framework to understand the relationship between tectonic processes and the different episodes of hydrothermalism and mineralisations that occurred in this part of the Cantabrian Mountains since the end of the Variscan cycle and throughout the entire Alpine cycle

Software
The geological map was originally constructed using Esri ArcGIS 10.5. The final draft was composed using Adobe Illustrator CS6.

Disclosure statement
No potential conflict of interest was reported by the author (s).

Funding
This work has been funded by the projects CGL2015-70970-P (FALPINO) from the I+D+i Spanish National Plan of the Ministerio de Ciencia, Innovación y Universidades and from the Principality of Asturias [GRUPIN-IDI/2018-000216], including both European Regional Development Funds. Additional funding was provided by the CANOA-321006500 project, included in the PLAN CARTOGRÁ-FICO 2017-2021 from the Spanish Geological Survey-Instituto Geológico y Minero de España (IGME).