Influence of basement rocks on fluid evolution during multiphase deformation: the example of the Estamariu thrust in the Pyrenean Axial Zone

Calcite veins precipitated in the Estamariu thrust during two tectonic events are studied in order to (i) decipher the temporal and spatial relationships between deformation and fluid migration in a long-lived thrust and (ii) determine the influence of basement rocks on the fluid chemistry during deformation. Structural and petrological observations constrain the relative timing of fluid migration and vein formation, whilst geochemical analyses (δ13C, δ18O, 87Sr/86Sr, clumped isotope thermometry, and elemental composition) applied to the related calcite cements and host rocks indicate the fluid origin, pathways, and extent of fluid–rock interaction. The first tectonic event, recorded by calcite cements Cc1a and Cc2, is attributed to the Alpine reactivation of the Estamariu thrust. Analytical data indicate that these cements precipitated from heated meteoric fluids (temperatures in the range of 50 to 100 C) that had interacted with basement rocks (87Sr/86Sr> 0.71) before upflowing through the thrust zone. The second tectonic event, attributed to the Neogene extension, is characterized by the reactivation of the Estamariu thrust and the formation of normal faults and shear fractures sealed by calcite cements Cc3, Cc4, and Cc5. Analytical data indicate that cements Cc3 and Cc4 precipitated from hydrothermal fluids (temperatures between 130 and 210 C and between 100 and 170 C, respectively) that had interacted with basement rocks (87Sr/86Sr > 0.71) and been expelled through fault zones during deformation. In contrast, cement Cc5 probably precipitated from meteoric waters that likely percolated from the surface through small shear fractures. The comparison between our results and already published data in other structures from the southern Pyrenees suggests that regardless of the origin of the fluids and the tectonic context, basement rocks have a significant influence on the fluid chemistry, particularly on the 87Sr/86Sr ratio. Accordingly, the cements precipitated from fluids that have interacted with basement rocks have significantly higher 87Sr/86Sr ratios (> 0.710) with respect to those precipitated from fluids that have interacted with the sedimentary cover (< 0.710), which involves younger and less radiogenic rocks.

Abstract. Calcite veins precipitated in the Estamariu thrust during two tectonic events are studied in order to (i) decipher the temporal and spatial relationships between deformation and fluid migration in a long-lived thrust and (ii) determine the influence of basement rocks on the fluid chemistry during deformation. Structural and petrological observations constrain the relative timing of fluid migration and vein formation, whilst geochemical analyses (δ 13 C, δ 18 O, 87 Sr/ 86 Sr, clumped isotope thermometry, and elemental composition) applied to the related calcite cements and host rocks indicate the fluid origin, pathways, and extent of fluid-rock interaction. The first tectonic event, recorded by calcite cements Cc1a and Cc2, is attributed to the Alpine reactivation of the Estamariu thrust. Analytical data indicate that these cements precipitated from heated meteoric fluids (temperatures in the range of 50 to 100 • C) that had interacted with basement rocks ( 87 Sr/ 86 Sr > 0.71) before upflowing through the thrust zone. The second tectonic event, attributed to the Neogene extension, is characterized by the reactivation of the Estamariu thrust and the formation of normal faults and shear fractures sealed by calcite cements Cc3, Cc4, and Cc5. Analytical data indicate that cements Cc3 and Cc4 precipitated from hydrothermal fluids (temperatures between 130 and 210 • C and between 100 and 170 • C, respectively) that had interacted with basement rocks ( 87 Sr/ 86 Sr > 0.71) and been expelled through fault zones during deformation. In contrast, cement Cc5 probably precipitated from meteoric waters that likely percolated from the surface through small shear fractures.
The comparison between our results and already published data in other structures from the southern Pyrenees suggests that regardless of the origin of the fluids and the tectonic context, basement rocks have a significant influence on the fluid chemistry, particularly on the 87 Sr/ 86 Sr ratio. Accordingly, the cements precipitated from fluids that have interacted with basement rocks have significantly higher 87 Sr/ 86 Sr ratios (> 0.710) with respect to those precipitated from fluids that have interacted with the sedimentary cover (< 0.710), which involves younger and less radiogenic rocks.

Introduction
Deformation associated with crustal shortening is mainly accommodated by thrust faulting and related fault zone structures (Mouthereau et al., 2014;Muñoz, 1992;Sibson, 1994). Successive faulting may occur, and favorably oriented structures may undergo reactivation during different tectonic events in a long-lived orogenic belt (Cochelin et al., 2018;Sibson, 1995). The reactivation of faults may produce changes in the hydraulic behavior of fault zones as well as in the origin and regime of fluids circulating through them Barker and Cox, 2011;Cantarero et al., 2018;Cruset et al., 2018;Lacroix et al., 2018;Travé et al., 2007). Consequently, constraining the timing of deformation and fluid migration is essential to better understand the main factors leading to the current configuration of a mountain belt, its evolution through time, and the mobilization of different fluids during successive deformation events (Baqués et al., 2012;Crespo-Blanc et al., 1995;Faÿ-Gomord et al., 2018;Fitz-Diaz et al., 2011;Lacroix et al., 2014). Understanding basin-scale fluid flow is of primary importance to reconstruct the diagenetic history of a sedimentary basin, as fluids take part in a wide range of geological processes including precipitation of new mineral phases, dolomitization, and petroleum migration, among others (Barker et al., 2009;Foden, 2001;Fontana et al., 2014;Gomez-Rivas et al., 2014;Martín-Martín et al., 2015;Mozafari et al., 2019;Piessens et al., 2002). Due to economic interest in these processes, in particular related to oil and ore deposit exploration, CO 2 sequestration, seismic activity, and water management, many researchers have addressed the relationship between deformation and fluid migration (Beaudoin et al., 2014;Breesch et al., 2009;Cox, 2007;Dewever et al., 2013;Gasparrini et al., 2013;Muñoz-López et al., 2020;Suchy et al., 2000;Travé et al., 2009;Voicu et al., 2000;Warren et al., 2014).
In the Pyrenees, the basement rocks from the Axial Zone are affected by numerous fault systems related to the Variscan orogeny (late Paleozoic) but reactivated during the Pyrenean compression (Late Cretaceous to Oligocene) (Cochelin et al., 2018;Poblet, 1991). However, no real consensus exists about the influence of the Alpine deformation on the basement rocks, and the age of basement-involved structures is still debated (Cochelin et al., 2018;García-Sansegundo et al., 2011). As a consequence, the relationships between deformation and fluid flow have widely focused on structures from the Mesozoic and Cenozoic cover (Beaudoin et al., 2015;Crognier et al., 2018;Cruset et al., 2016Cruset et al., , 2018Cruset et al., , 2020aLacroix et al., 2011Lacroix et al., , 2014Martinez Casas et al., 2019;Muñoz-López et al., 2020;Nardini et al., 2019;Travé et al., 1997Travé et al., , 1998Travé et al., , 2000, where the timing of deformation and thrust emplacement is well-constrained (Cruset et al., 2020b;Vergés, 1993;Vergés and Muñoz, 1990). In contrast, few studies, only concentrated along the Gavarnie thrust system, have examined the relationship between deformation and fluid migration in the Paleozoic basement (Banks et al., 1991;Grant et al., 1990;Henderson and McCaig, 1996;McCaig et al., 1995McCaig et al., , 2000aRye and Bradbury, 1988;Trincal et al., 2017). Another important aspect of studying the faultfluid system is related to the heat flow and the influence of faults on the development of geothermal systems (Faulds et al., 2010;Grasby and Hutcheon, 2001;Liotta et al., 2010;Rowland and Sibson, 2004). Particularly, in the NE part of the Iberian Peninsula (including the Pyrenees and the Catalan Coastal Range), high-permeability Neogene extensional faults, acting as conduits for upward migration, provide efficient pathways for hydrothermal fluids to flow from deeper to shallower crustal levels (Carmona et al., 2000;Fernàndez and Banda, 1990;Taillefer et al., 2017Taillefer et al., , 2018. In this sense, understanding the fault-fluid system evolution and the relative timing of hydrothermal fluid migration is of great importance to characterize the potential geothermal resources of this area. In this contribution, we report the temporal and spatial relationships between deformation and fluid migration in a long-lived Variscan thrust deforming basement rocks in the Pyrenean Axial Zone. For this purpose, we combine structural, petrological, and geochemical analyses of calcite veins precipitated in the Estamariu thrust during two reactivation episodes related to the Alpine compression (Late Cretaceous to Oligocene) and the Neogene extension. Structural and petrological observations allow us to unravel the relative timing of fluid migration and vein formation in relation to the involved tectonic events. The geochemistry of the vein cements and related host rocks provides information on the fluid origin, pathways, and extent of fluid-rock interaction during deformation. Therefore, the main objectives of this paper are the following: (i) to constrain the relative timing of vein formation and fluid migration; (ii) to determine the fluid origin and pathways during successive compressional and extensional deformation phases; (iii) to assess the influence of basement rocks on the chemistry of fluids circulating during deformation; and (iv) to provide insights into the fluid flow at regional scale in the NE part of the Iberian Peninsula, where the presence of hydrothermal fluids has been reported from Neogene times to present.

Geological setting
The Pyrenees constitute an asymmetric and doubly verging orogenic belt that resulted from the Alpine (Late Cretaceous to Oligocene) convergence between the Iberian and European plates (Choukroune, 1989;Muñoz, 1992;Roure et al., 1989;Sibuet et al., 2004;Srivastava et al., 1990;Vergés and Fernàndez, 2012). The Pyrenean structure consists of a central antiformal stack of basement-involved rocks from the Axial Zone, flanked by two oppositely verging fold-andthrust belts and their associated foreland basins (Muñoz, 1992;Muñoz et al., 1986) (Fig. 1A). The Pyrenean Axial Zone has been deformed by successive Variscan, Alpine, and Neogene phases (Saura and Teixell, 2006). In the eastern Axial Pyrenees, an E-W to ENE-WSW fault system developed during the Neogene extension (Roca, 1996;Roca and Guimerà, 1992;Vergés et al., 2002). The main fault, La Tet fault, has associated a set of E-W extensional basins such as La Cerdanya, Conflent, La Seu d 'Urgell, and Cerc (Cabrera et al., 1988;Roca, 1996). The Cerc basin consists of a Stephano-Permian accumulation of volcanic rocks discordantly overlying Cambro-Ordovician materials. This basin is thrusted in its eastern limit by the Estamariu thrust, whereas the northern and southern boundaries correspond to two Neogene extensional faults, La Seu d'Urgell fault and the Ortedó fault, respectively (Hartevelt, 1970;Roca, 1996;Saura, D. Muñoz-López et al.: Influence of basement rocks on fluid evolution during multiphase deformation 2259 2004) (Fig. 1B, C). In the NW part of the basin, the limit between the Stephano-Permian unit and the upper Ordovician sequence corresponds to a Stephano-Permian extensional fault formed coevally with the deposition of the volcanic sequence (Saura, 2004). This fault was reactivated during the latest stages of the Neogene extension (Saura, 2004) and is referred to here as the Sant Antoni fault (Fig. 1C).
The Estamariu thrust is a basement-involved reverse fault originated during the Variscan orogeny with a minimum displacement of 27 km (Poblet, 1991). However, in its southwestern termination, it juxtaposes the Devonian Rueda Formation against the Stephano-Permian Erill Castell Formation. The Erill Castell Formation developed during the late to post-orogenic collapse of the Variscan belt (Lago et al., 2004;Martí, 1991Martí, , 1996Ziegler, 1988), evidencing the reactivation of the Estamariu thrust during the Alpine orogeny (Poblet, 1991;Saura, 2004). Rocks cropping out around the Estamariu thrust and the Cerc basin range from upper Ordovician to Miocene (Fig. 1C). However, due to the complex structural setting, the stratigraphic record is discontinuous and only upper Ordovician, Devonian, Stephano-Permian, and Neogene rocks are present in the study area. The basement lithologies consist of upper Ordovician and Devonian metasedimentary rocks affected by multiscale folds and related pervasive axial plane regional foliation (Bons, 1988;Casas et al., 1989;Cochelin et al., 2018;Zwart, 1986). This deformation is linked to low-grade metamorphic conditions developed during the Variscan orogeny (Hartevelt, 1970;Poblet, 1991;Saura, 2004). The upper Ordovician succession includes an alternation of shales, sandstones, conglomerates, quartzites, and phyllites, and the Devonian sequence consists of an alternation of limestones and black slates (Rueda Formation) (Mey, 1967). The Stephano-Permian sequence developed during the late to post-orogenic extensional collapse of the Variscan belt and in the study area is represented by a volcanic and volcanoclastic unit (the Erill Castell Formation) (Martí, 1991;Mey et al., 1968) involving tuffs and ignimbrites at the base and andesites in the upper part (Martí, 1996;Saura and Teixell, 2006). Finally, the Neogene sequence is constituted by detrital and poorly lithified sediments, mainly shales, sandstones, and conglomerates deposited during the Neogene extension associated with the opening of the NW Mediterranean Sea (Roca, 1996).

Methods
This study integrates a field compilation of structural data as well as petrological and geochemical analyses of calcite cements and related host rocks. The structural data include the orientation of bedding, foliations, and fractures in addition to crosscutting relationships and kinematics. Such data were plotted in equal-area lower-hemisphere projections, and different fracture sets were established according to their type, strike, mineral infillings, and relative age deduced from crosscutting relationships. All these data were integrated in a schematic map and a cross section of the Estamariu thrust and the Cerc basin (Figs. 2a,b and 3). Samples considered representative of the involved host rocks and all calcite vein generations observed in the different fracture sets and fault-related structures were selected for petrological and geochemical analyses. Thin sections of these samples were prepared and studied under a Zeiss Axiophot optical microscope and a cold cathodoluminescence (CL) microscope model 8200 Mk5-1 operating between 16-19 kV and 350 µA gun current.
The geochemical analyses have been performed in calcite cements and related host rocks in order to determine the origin, composition, and temperature of the vein-forming fluids in addition to the extent of fluid-rock interaction. These analyses include (i) stable isotope analyses in 37 calcite cements and the carbonate portion of the Devonian rocks, (ii) the elemental composition of 12 samples analyzed using high-resolution inductively coupled plasma mass spectrometry (HR-ICP-MS), (iii) the 87 Sr/ 86 Sr and 143 Nd/ 144 Nd ratios of eight representative samples, and (iv) clumped isotope thermometry of calcite cements (three samples), the results of which were converted to temperatures by applying the calibration method of Kluge et al. (2015). Calculated δ 18 O fluid values are expressed in per mill (‰) with respect to the Vienna Standard Mean Ocean Water (VSMOW). Details on methods and procedures can be found in the Supplement.

Structure and associated calcite cements
The Estamariu thrust strikes N-S to NW-SE and dips between 40 and 70 • towards the NE. It has a displacement of a few hundred meters and juxtaposes a Devonian alternation of limestones and shales in the hanging wall against Stephano-Permian andesites in the footwall (Poblet, 1991) (Figs. 2-4). The main slip plane is undulose and generates a 2-3 m thick thrust zone affecting both the hanging wall and footwall, but it is thicker in the hanging wall, up to 2.5 m thick. In the footwall the thrust zone is less than 1 m thick and has associated minor restricted thrust zones developed as subsidiary accommodation structures related to the main thrust fault (Fig. 2a,  b). All kinematic indicators, including S-C structures and slickenlines, indicate reverse displacement towards the west.
The mesostructures and microstructures observed in the study area are described below according to their structural position in relation to the Estamariu thrust, that is, hanging wall, thrust zone, and footwall (Figs. 3, 4). The relative timing of the different mesostructures and microstructures has been determined by means of crosscutting relationships and microstructural analysis.  . Lower-hemisphere equal-area stereoplots of the Devonian bedding (S 0 ), regional foliation (S 1 ), and thrust zone foliation affecting the hanging wall and footwall (S 2 ); magmatic layering and the different faults and veins observed in the study area are also included. Location in Fig. 1B.

Hanging wall
In the studied outcrops, the Devonian Rueda Formation is characterized by a well-bedded alternation of dark to light grey limestones with subordinate dark grey shales (S0) (Fig. 5a). Limestones are made up of encrinites, which consist of a bioclastic packstone formed essentially of crinoid stems (Fig. 5b). Under cathodoluminescence, en-crinites show dark to bright orange colors (Fig. 5c). Devonian rocks form a decametric anticline oriented NW-SE with a well-developed axial plane foliation (S 1 ) concentrated in the pelitic intervals. S 1 , which is oriented NNW-SSE, is a pervasive regional foliation dipping 30 to 55 • towards the E and NE and is generally between 2 and 5 cm spaced. In the hinge of the anticline, bedding (S0) dips towards the SE and forms a high angle with S 1 (Figs. 2b and 5a), whereas in its eastern limb, the regional foliation (S 1 ) dips steeper than S0. These geometric relationships between bedding and foliation have been used to determine the fold type at large scale (i.e., as shown in Fig. 2b).

Thrust zone
The thrust zone consists of a deformation zone affecting both the hanging wall and the footwall. Within the hanging wall, the Devonian host rocks are still recognizable, but the intensity of deformation progressively increases towards the main thrust plane. This deformation consists of a penetrative thrust zone foliation (S 2 ), two generations of stylolites (e1, e2), and three generations of calcite veins (V0, V1a, and V1b) (Figs. 3, 6). These structures are described below in chronological order.
The foliation within the thrust zone affecting the Devonian hanging wall (S 2 ) strikes NW-SE and dips 40-50 • NE, similar to the regional foliation (S 1 ), but it is more closely spaced, generally between 0.2 and 1 cm (Fig. 6a, b). This observation points to a progressive transposition of the regional foliation within the thrust zone during thrusting. At mesoscale, S 2 has related shear surfaces (Ci) defining centimetric S-Ctype structures, again indicating reverse kinematics (Fig. 6a). Stylolites e1 have a wave-like shape and trend subparallel to the thrust zone foliation (S 2 ) (Fig. 6b, c). When present, these stylolites are very systematic, exhibiting a spacing of 1-2 mm (Fig. 6c).
The first calcite vein generation (V0), only observed at microscopic scale (Fig. 6b, c), corresponds to up to 1 cm long and less than 1 mm thick veins cemented by blocky to elongated blocky calcite crystals featuring a dark brown luminescence (cement Cc0). Veins V0 and stylolites e1 are perpendicular between them and show ambiguous crosscutting relationships. These microstructures are concentrated into discontinuous fragments of the Devonian host rocks within the thrust zone. Calcite veins V1a crosscut the previous vein generation (V0) as well as the stylolites e1 and are developed within S 2 surfaces (Figs. 3, 6d). These veins are the most abundant, exhibit a white to brownish color in hand samples, and are up to 10 cm long and 1 cm thick. The vein cement (Cc1a) is formed of anhedral crystals up to 3-4 mm in size displaying a blocky texture and a dark brown luminescence (Fig. 6e). Stylolites e2, more abundant than stylolites e1, are up to 10 cm long and show spacing between 0.5 and 2 cm (Fig. 6d, f). These stylolites mainly correspond to sutured areas developed between the host rock and the calcite veins V1a and between foliation surfaces S 2 .
Calcite veins V1b, up to 1 cm long and less than 1 mm thick, were also identified at microscopic scale ( Fig. 6d, f). The vein cement (Cc1b) consists of up to 0.1 mm calcite crystals with a blocky texture and a bright yellow luminescence. These veins postdate the previous V0 and V1a generations and trend perpendicular to stylolites e2.
Towards the fault plane, the thrust zone foliation S 2 is progressively more closely spaced, and stylolites e2 become more abundant (showing millimeter spacing) and exhibit ambiguous crosscutting relationships with veins V1b (Fig. 6f). The main slip surface corresponds to a discrete plane that contains calcite slickensides (veins V2). The vein cement (Cc2) is milky white in hand samples and consists of up to 3 mm blocky to elongated blocky crystals (Fig. 6g) with a dull to bright orange luminescence (Fig. 6h). Deformation in the footwall is concentrated within the main thrust zone and subsidiary thrust zones and corresponds to the thrust zone foliation (S 2 ) and calcite veins V3 (Fig. 3). This foliation (S 2 ) strikes NW-SE, dips towards the NE, and is millimeter to centimeter spaced (Fig. 7a). Calcite veins V3 are generally 1-2 cm thick and strike NW-SE. They are parallel or locally branch off, cutting the S 2 planes in the subsidiary thrust zone (Fig. 7a, b). Outside the thrust zone, veins V3 are locally present but have a NE-SW strike. These veins are mostly less than 1 m long and are spaced between a few centimeters and 50 cm. The vein cement (Cc3) is made up of a milky white calcite characterized by up to 3 mm long fibrous crystals oriented perpendicular to the vein walls (Fig. 7c). Locally, anhedral blocky crystals ranging in size from 0.1 to 1 mm are also present. This cement displays a bright yellow to bright orange luminescence (Fig. 7d).

Footwall
In the footwall, the Stephano-Permian Erill Castell Formation comprises massive, dark greenish andesitic levels showing a rhythmic magmatic layering (Fig. 7e), which corresponds to a fluidal structure of the host rock. The local presence of pyroclastic and brecciated volcanoclastic levels is also ubiquitous, mainly in the lower part of this sequence. Andesites are characterized by a porphyritic texture defined by a dark fine-grained spherulitic matrix partially devitrified, large zoned crystals of plagioclase (Fig. 7f) up to 2-3 cm long, and less abundant biotite and hornblende. These mafic phenocrystals are systematically pseudomorphosed by clay minerals and frequently show evidence of oxidation and chloritization. Andesites are affected by E-W-striking open joints (J1) dipping indistinctively towards the north and south (Fig. 7e). These joints locally trend parallel to the magmatic layering (Fig. 3).  Finally, as described above, the northern and southern limits of both the Cerc basin and the Estamariu thrust correspond to two Neogene extensional faults, La Seu d'Urgell and the Ortedó fault systems (Fig. 1C). These faults are subvertical or steeply dip towards the north. In the northern part, the slip plane of the La Seu d'Urgell fault has not been observed, and the limit between the Stephano-Permian rocks and the Neogene deposits is not well-constrained due to the poor quality of the Neogene outcrops and the presence of Quaternary deposits. In the southern part, the Ortedó fault generates a dark greyish to brown fault zone that is several meters thick and characterized by the presence of clay-rich incohesive fault rocks developed at the contact between Stephano-Permian and upper Ordovician rocks. Related to these main fault systems, mesoscale normal faults commonly affect the andesites within the Cerc basin. These faults are mainly E-W and locally NE-SW, are subvertical, and dip indistinctly towards the N and S. Fault planes are locally mineralized with calcite cement (veins V4) and exhibit two striae set generations indicating dip-slip and strike-slip movements (Fig. 8a). The calcite cement (Cc4) consists of up to 2 mm blocky to elongated blocky crystals (Fig. 8b) with a homogeneous dark orange luminescence (Fig. 8c). On the other hand, the main Estamariu thrust zone is locally displaced by shear fractures (Fig. 8d) and a later set of shear bands (Cn) (Fig. 8e), both having an overall NNW-SSE to NNE-SSW strike (Fig. 3) that indicates a minor normal displacement. Shear fractures are locally mineralized with calcite (veins V5). The vein cement consists of a greyish microsparite calcite cement (Cc5) (Fig. 8f, g).

Geochemistry of calcite cements and host rocks
The geochemistry (δ 18 O, δ 13 C, δ 18 O fluid , 87 Sr/ 86 Sr, 143 Nd/ 144 Nd, and elemental composition) and the calculated temperature of precipitation of the different calcite cements Cc1a to Cc5 are described below. Veins V0 and V1b were only observed at microscopic scale and their calcite cement Cc0 and Cc1b could not be sampled to perform these geochemical analyses.
The δ 18 O and δ 13 C isotopic composition of the carbonate fraction of the Devonian hanging wall and the different calcite cements (Cc1a to Cc5) are summarized in Table 1 and represented in Fig. 9. The micritic matrix of the Devonian packstone ranges in δ 18 O values between −10.5 ‰ VPDB and −8.4 ‰ VPDB and in δ 13 C values between +1.5 ‰ VPDB and +2.8 ‰ VPDB, whereas the calcite cements have a broader range of values depending on the cement generation (Fig. 9).
The obtained 47 values from clumped isotope thermometry were converted into temperatures and δ 18 O fluid (Table 1 and Fig. 10) using the equations of Kluge et al. (2015) and Friedman and O'Neil (1977), respectively. The calculated T and δ 18 O fluid for calcite cement Cc2 range between 50 and 100 • C and between −6.4 ‰ VSMOW and −0.3 ‰ VSMOW, respectively. Cement Cc3 shows T and δ 18 O fluid in the range of 130 to 210 • C and +4.3 ‰ VSMOW to +12.1 ‰ VSMOW, respectively. Cc4 exhibits T and δ 18 O fluid in the range of 100 to 170 • C and +0.9 ‰ VSMOW to +8.1 ‰ VSMOW, respectively.
The 87 Sr/ 86 Sr ratios of calcite cements Cc1a to Cc5 and host rocks are reported in Table 1 and Fig. 11. Devonian limestones from the hanging wall have a 87 Sr/ 86 Sr ratio of 0.710663, whilst the Stephano-Permian andesites in the footwall exhibit a more radiogenic 87 Sr/ 86 Sr ratio of 0.743983. The calcite cements have more radiogenic 87 Sr/ 86 Sr ratios than the Devonian limestones but lower radiogenic values than the Stephano-Permian andesites. This ratio is 0.713018 to 0.714092 in Cc1a, 0.718294 for Cc2, 0.714619 for Cc3, 0.717706 for Cc4, and 0.716923 for Cc5. These results are compared with already published data from synkinematic veins and deformed rocks in other Pyrenean structures developed in the basement and in the sedimentary cover during the Pyrenean compression (Fig. 11). This comparison shows that values obtained in this study are (1) significantly more radiogenic than the values of marine carbonates and synkinematic veins precipitated in the sedimentary cover (i.e., in the south Pyrenean fault and thrust belt) and (2) within the same range of values of synkinematic veins and deformed rocks in the Pyrenean basement (Axial Zone).
The analyzed samples for 143 Nd/ 144 Nd ratios in calcite cements and host rocks are reported in Table 1. However, due to the general low Nd concentrations in most of the analyzed calcite cements and the limited number of powdered samples that were available, only calcite cement Cc5 and the andesite host rock (footwall) could be measured. Cc5 has a 143 Nd/ 144 Nd ratio of 0.512178, which is similar that of its footwall host rocks: 0.512196.  The obtained elemental composition broadly varies among the different calcite cements and related host rocks (Table 2 and Fig. 12). In the thrust zone affecting the hanging wall, calcite cement Cc1a shows a similar trend to that of the Devonian limestones, both having high Sr, intermediate to high Mg and Fe, and low Mn contents (Fig. 12)   lasted from late Paleozoic (Variscan orogeny) to Neogene times.
The Paleozoic metasedimentary rocks from the Pyrenean basement are broadly affected by multiscale folds and axial plane regional foliation developed during the main Variscan deformation phase (Bons, 1988;Cochelin et al., 2018;Zwart, 1986). Similar structures, a decametric-scale anticline, and pervasive axial plane foliation (S 1 ) are found in the Devonian sequence located in the thrust hanging wall (Fig. 2b). We therefore consider them to have developed during the Variscan compression, contemporaneous with the main activity of the Estamariu thrust. Veins V0 are perpendicular to stylolites e1 and show ambiguous crosscutting relationships between them. Thus, they are interpreted as originat- Figure 11. Simplified geological map of the south-central Pyrenees showing the location of structures where 87 Sr/ 86 Sr analyses have been carried out. Below are 87 Sr/ 86 Sr ratios from this study compared to results from other structures involving either cover units (1-10) or basement rocks (11-16). The thick blue line refers to the 87 Sr/ 86 Sr range of Phanerozoic seawater, and the dashed brown line represents the 87 Sr/ 86 Sr limit value between basement and cover structures. (1) El Guix anticline , (2) Puig Reig anticline (Cruset et al., 2016), (3) L'Escala thrust , (4) Vallfogona thrust , (5) Ainsa basin (Travé et al., 1997), (6) Ainsa-Bielsa area (McCaig et al., 1995), (7) Minor Bóixols thrust (Muñoz-López et al., 2020), (8) Bóixols anticline (Nardini et al., 2019), (9) lower Pedraforca thrust (Cruset et al., 2020a), (10) upper Pedraforca thrust (Cruset, 2019), (11) Gavarnie thrust (McCaig et al., 1995), (12) Pic de Port Vieux thrust (Banks et al., 1991), (13) Pic de Port Vieux thrust (McCaig et al., 2000b), (14) Plan de Larri thrust (McCaig et al., 1995), (15) La Glere shear zone (Wayne and McCaig, 1998), (16) Trois Seigneurs Massif (not on the map) (Bickle et al., 1988). ing coevally. Both microstructures are concentrated into discontinuous fragments of the Devonian host rocks and are therefore considered inherited microstructures likely developed in Variscan times. However, as pointed out above, in the study area the Estamariu thrust affects late to post-Variscan Stephano-Permian andesites, thus confirming its reactivation during the Alpine orogeny. Accordingly, the structures that are strictly attributed to the Alpine reactivation of the thrust are those structures indicating reverse kinematics or associated with a compressional stress, which are found within the thrust zone deformation at the contact between Devonian and Stephano-Permian units. Contrarily, the magmatic lay-ering and joints J1 are broadly present in the andesitic footwall, outside the thrust zone, and in other Stephano-Permian basins; they are therefore considered inherited fluidal and cooling structures, respectively. For this reason, the calcite veins V1a and V2 (and related cements Cc1a and Cc2), exclusively occurring within the thrust zone, have been associated with the reactivation of the Estamariu thrust. During this period, and associated with ongoing deformation and progressive shortening, stylolites e2 developed as sutured areas between host rock and veins V1a and between foliation surfaces, coevally with the development of veins V1b, as denoted by their crosscutting relationships and orientations.  (Ca, Mg, Fe, Mn, Sr) of the calcite cements Cc1a to Cc5 and host rocks from the hanging wall (HW) and footwall (FW). The qualitative scale in green indicates different contents (for each element, the darkest green points to the highest concentration and vice versa). Other structures present in the study area, such as veins V3 to V5 and related cements Cc3 to Cc5, are attributed to the Neogene extension. Veins V3 precipitated in the subsidiary thrust zone developed in the footwall of the Estamariu thrust. These veins strike parallel to the thrust zone foliation (S 2 ) (Fig. 7a, b) but are characterized by calcite fibers growing perpendicular to the vein walls and to the foliation surfaces (Fig. 7c), thus evidencing their extensional character. The presence of extensional calcite veins opened along previously formed foliation surfaces in a thrust zone has been reported in other structures in the Pyrenees and has been considered to postdate the thrust activity (Lacroix et al., 2011(Lacroix et al., , 2014. Veins V4 precipitated in subvertical and E-W mesoscale faults affecting the Stephano-Permian andesites (Fig. 8) outside the thrust zone (Fig. 3). The fault orientation and dip as well as the two striae set generations observed on the fault planes are compatible with the Neogene extensional faults that bound the Cerc basin and postdate the Esta-mariu thrust (Cabrera et al., 1988;Roca, 1996;Saura, 2004). Calcite cements Cc3 and Cc4, occluding veins V3 and V4, have a similar geochemical composition (Figs. 9-12), supporting the idea that their precipitation occurred during the same tectonic event associated with a similar fluid regime (i.e., although these cements precipitated in different structures, they are likely contemporaneous). Finally, veins V5 (and related cement Cc5) precipitated locally in shear fractures crosscutting and postdating the thrust-related deformation (Figs. 3, 4, and 8d). These veins strike parallel to the shear bands (Cn) located in the main thrust zone (Fig. 8e), exhibiting normal slip kinematics, postdating the reverse structures, and therefore indicating reactivation of the Estamariu thrust during the Neogene extension.

Fluid system during the Alpine reactivation of the Estamariu thrust
As veins V1a and V2 are consistent with the Alpine reactivation of the Estamariu thrust, the geochemistry of their related calcite cements Cc1a and Cc2 is interpreted to record the fluid system during this tectonic event.
Cements Cc1a and Cc2 are characterized by high 87 Sr/ 86 Sr ratios (from 0.713 to 0.714 for Cc1a and 0.718 for Cc2), significantly higher radiogenic than ratios of Phanerozoic seawater (between 0.7070 and 0.7090) (McArthur et al., 2012). This may reflect the incorporation of radiogenic Sr from a fluid that derived from or interacted with Rb-rich and/or Sr-rich basement rocks such as those underlying the Estamariu thrust. In addition, Cc1a has a narrow range of δ 13 C, between +0.91 ‰ VPDB and +2 ‰ VPDB, consistent with values of the Devonian marine limestones from the hanging wall (between +1.54 ‰ VPDB and +2.75 ‰ VPDB) and within the range of Devonian marine carbonate values (Veizer et al., 1999). Likewise, the elemental composition of Cc1a follows a similar trend to that of its Devonian host, both having high Mg and Sr and low Mn contents with respect to the other calcite cements (Fig. 12). These geochemical similarities indicate significant fluid-rock interaction and buffering of the carbon and elemental composition of the precipitating fluid by the Devonian carbonates (Marshall, 1992). Calcite cement Cc2 has slightly lower δ 13 C, lower Mg and Sr, and higher Mn contents with respect to both Cc1a and the Devonian host. On the other hand, the temperature and the δ 18 O composition of the vein-forming fluids, calculated from clumped isotope thermometry of Cc2, range between 50 and 100 • C and between −6.4 ‰ SMOW and −0.3 ‰ SMOW, respectively. These values are interpreted as the involvement of meteoric fluids that were probably heated at depth and enriched in radiogenic Sr during their flow and interaction with basement rocks. These fluids flowed preferentially along the thrust zone (Fig. 13a), as evidenced by the exclusive presence of calcite in this area, likely due to the enhanced permeability associated with the thrust discontinuity (McCaig et al., 1995;Trincal et al., 2017). As Cc1a and Cc2 precipitated in the thrust zone during the same tectonic event, they likely precipitated from the same fluids, which progressively increased the extent of fluid-rock interaction from the thrust plane (Cc2) towards the hanging wall (Cc1a), as indicated by the host-rock-buffered composition of the latter. Previous studies already reported syntectonic migration of fluids that had interacted at depth with basement rocks before upflowing along thrust zones in other structures from the Pyrenean basement, such as the Gavarnie thrust and the related Pic de Port Vieux thrust (McCaig et al., 1995).

Fluid system during the Neogene extension
Calcite veins V3 to V5 are attributed to the Neogene extension, and the geochemistry of their related calcite cements Cc3 to Cc5 characterizes the fluid system during this period.
Cc3 and Cc4 have considerably high 87 Sr/ 86 Sr ratios (0.714619 and 0.717706, respectively), similar to the ones reported for Cc1a and Cc2 (Fig. 11), indicating interaction with basement rocks. The δ 18 O fluid calculated from clumped isotopes, between +4.3 ‰ SMOW and +12.1 ‰ SMOW for Cc3 and between +0.9 ‰ SMOW and +8.1 ‰ SMOW for Cc4, falls within the range of metamorphic and/or formation brines (Taylor, 1987). The δ 18 O-depleted values of these cements (around −14 ‰ VPDB) are due to the high temperatures of the fluids (between 130 and 210 • C for Cc3 and between 100 and 170 • C for Cc4). Assuming a normal geothermal gradient of 30 • C, these temperatures would have been reached at a minimum depth of 3-5 km. However, these veins have never reached such a burial depth, since during the Neogene extension the studied structure acquired its current configuration (Saura, 2004) and was only buried under the Devonian sequence (hanging wall), which has a maximum thickness of several hundred meters (Mey, 1967). This assumption evidences the hydrothermal character of the circulating fluids, which probably migrated rapidly enough through normal faults to maintain their high temperatures and to be in thermal disequilibrium with the surrounding rocks. Similarly, the high Mn content of Cc3 and Cc4 (around 7700-8300 and 4000 ppm, respectively), responsible of their bright luminescence (Figs. 7d, 8c), is consistent with hydrothermal waters (Pfeifer et al., 1988;Pomerol, 1983;Pratt et al., 1991). On the other hand, the δ 13 C-depleted values of these cements are indicative of the influence of organic-derived carbon (Cerling, 1984;Vilasi et al., 2006). The most probable source for these low δ 13 C values is the Silurian black shales that do not crop out in the study area but acted as the main detachment level during the Variscan compression and locally during the Alpine compression in the Pyrenean Axial Zone (Mey, 1967). These black shales have significant organic carbon contents (TOC around 2.3 %), and around the Gavarnie thrust, they exhibit syntectonic carbonate veins yielding δ 13 C values between −2 ‰ VPDB and −8 ‰ VPDB (McCaig et al., 1995). Thus, cements Cc3 and Cc4 precipitated from hy- Figure 13. Tectonic and geochemical evolution of the study area (not to scale) and relationships with the evolution of the fluid system. (a) During the Alpine reactivation of the Estamariu thrust, a meteoric fluid (red arrows) interacted at depth with basement rocks and then migrated along the fault plane towards the hanging wall, precipitating cements Cc1a and Cc2. (b) During the Neogene extension, basementderived hydrothermal fluids (orange arrows) flowed upwards through newly formed and reactivated fault zones. This fluid precipitated calcite cements Cc3 and Cc4. Finally, during ongoing deformation, meteoric fluids (green arrows) percolated in the system and precipitated Cc5, revealing a main change in the fluid regime. drothermal fluids derived from and/or equilibrated with basement rocks expelled through newly formed and reactivated fault zones during deformation (Fig. 13b). The hydrothermal character of the fluids involved in this deformation event, and their relatively high temperatures (up to 210 • C), could have altered the clumped isotope composition of the previous calcite cement generations (Cc1a and Cc2). Clumped isotopes may be reset by recrystallization and by solid-state isotopic exchange reactions and diffusion within the mineral lattice (Shenton et al., 2015;Stolper and Eiler, 2015). However, in the studied vein samples there is no evidence of calcite recrystallization (for instance, grain coarsening linked to grain boundary migration). Therefore, we conclude that although it is possible, there is little evidence for clumped isotopes of cements Cc1a and Cc2 to have been reset during precipitation of Cc3 and Cc4.
Finally, the isotopic signature of Cc5, ranging between −8.1 ‰ VPDB and −5.7 ‰ VPDB for δ 18 O and between −8.2 ‰ VPDB and −3.8 ‰ VPDB for δ 13 C, falls within the range of meteoric carbonates (Travé et al., 2007;Veizer, 1992). The similar tendency in the elemental composition of this cement and the Stephano-Permian volcanic rocks, both having the highest Mg and Fe and the lowest Mn and Sr contents with respect to the other cements and host rocks, reveals significant fluid-rock interaction with the footwall rocks. The significant water-rock interaction is also demonstrated by the Nd isotopic composition of Cc5 (0.512178), yielding values similar to the volcanic host (0.512196). This fact indicates that this cement precipitated from percolation of meteoric fluids, the geochemistry of which was controlled by the volcanic host rock. Studies focused on infiltration of meteoric fluids and subsequent upflowing along the La Tet fault during the Neogene extension have shown that meteoric waters in the area infiltrate at high altitudes of around 2000 m and low temperatures of around 5 • C (Krimissa et al., 1994;Taillefer et al., 2018). Therefore, we suggest that Cc5 probably precipitated during the latest stages of extension after the fluid regime changed from upward fluid migration to percolation of cold meteoric waters, as also occurred in the Barcelona Plain (Catalan Coastal Range) (Cantarero et al., 2014).
In conclusion, the high 87 Sr/ 86 Sr ratios of the different calcite cements (Cc1a to Cc5) evidence the interaction between the vein-forming fluids and Paleozoic basement rocks with a more radiogenic signature. The geochemical evolution of these cement generations also highlights the progressive change in the fluid regime and composition during successive compressional and extensional tectonic events (Fig. 13). The continuous increase in precipitation temperatures and enrichment in δ 18 O fluid from calcite cements Cc1a and Cc2 (Alpine) to cements Cc3 and Cc4 (Neogene) is probably linked to a higher extent of fluid-rock interaction with basement rocks. By contrast, during the latest stages of extension, the infiltration of meteoric fluids likely indicates a more significant change in the fluid regime from upward to downward fluid migration (Fig. 13).

Influence of Paleozoic basement rocks on the 87 Sr/ 86 Sr ratios of the fluids circulating during deformation
In this section, we assess the influence of basement rocks on the chemistry of fluids circulating during deformation in the Pyrenees. The comparison between previous studies and the new data provided in this contribution evidences that fluids migrating through basement or cover units have a different geochemical signature, which is recorded in the 87 Sr/ 86 Sr ratios of the vein cements. The high 87 Sr/ 86 Sr ratios (0.713 to 0.718) of the analyzed cements, originated during successive compressional and extensional tectonic events, indicate that regardless of the origin of the fluids and the tectonic context, basement rocks have a significant influence on the fluid chemistry. Accordingly, cements precipitated from fluids that have circulated through basement rocks have significantly high 87 Sr/ 86 Sr ratios (> 0.710) (Fig. 11), reflecting the interaction between the vein-forming fluids and rocks with a more radiogenic signature. Similar radiogenic 87 Sr/ 86 Sr ratios have also been attributed to basementderived fluids in the Glarus nappe (Swiss Alps) (Burkhard et al., 1992). By contrast, vein cements precipitated from fluids that have circulated through the Mesozoic-Cenozoic sedimentary cover in the Pyrenees (i.e., through younger rocks with a different radiogenic signature) have significantly lower 87 Sr/ 86 Sr ratios (< 0.710). Such lower values may be similar to Phanerozoic seawater values, evidencing interaction between the vein-forming fluids and marine carbonate units, or they may be higher, evidencing interaction with siliciclastic rocks Travé et al., 2007). A previous study, focused on fluid flow along the Gavarnie thrust in the central-western Axial Pyrenees, used this limit value ( 87 Sr/ 86 Sr = 0.710) to differentiate between the unaltered limestone protolith and the Cretaceous thrust-related carbonate mylonite affected by fluids carrying radiogenic Sr (Mc-Caig et al., 1995). In addition, Williams et al. (2015) also attributed 87 Sr/ 86 Sr ratios > 0.710 to the input of basementderived fluids in the Rio Grande rift, USA.

5.5
Fluid flow at regional scale: the NE part of the Iberian Peninsula during the Neogene extension As pointed out above, the structural and geochemical data indicate that calcite cements Cc3 and Cc4 (veins V3 and V4, respectively) precipitated from hydrothermal fluids (up to 210 • C) that interacted at depth with basement rocks before ascending through newly formed and reactivated structures during the Neogene extension. These interpretations are consistent with the presence of several hydrothermal springs (temperatures of 29 to 73 • C) currently upwelling aligned along the La Tet fault and related Neogene deformation in the Pyrenean Axial Zone (Krimissa et al., 1994;Taillefer et al., 2017Taillefer et al., , 2018. Several studies indicate the origin of these hotwater springs as meteoric fluids infiltrated at high-elevated reliefs above 2000 m, warmed at great depths by normal geothermal gradients, and migrated upwards along permeability anisotropies related to fault zones (Taillefer et al., 2017(Taillefer et al., , 2018. The 87 Sr/ 86 Sr ratios of these springs, ranging between 0.715 and 0.730 (Caballero et al., 2012), are within the range of values obtained in this study and also account for interaction between circulating fluids and basement lithologies. Studies based on numerical models suggest that the La Tet fault and the involved basement rocks are still permeable down to 3 km of depth (Taillefer et al., 2017(Taillefer et al., , 2018, although the fault has been dormant since the Mio-Pliocene (Goula et al., 1999). It has also been suggested that the footwall topography, which induces high hydraulic gradients and produces fluid advection, is the major factor controlling the infiltration of meteoric fluids, the circulation depth, and the maximum temperature reached by the migrating fluids (Taillefer et al., 2017). A similar geological context and fluid regime evolution to those explained above are found in the Barcelona Plain and the Vallès Basin, located in the northeast part of the Catalan Coastal Range (CCR) (Fig. 1A). Consequently, the comparison between the two geological contexts allows us to give insights into the fluid circulation in extensional basins at regional scale (in the NE part of Iberia). In these locations of the CCR, the main fault system associated with the Neogene extension acted as a conduit for hydrothermal fluid circulation at temperatures between 130 and 150 • C during synkinematic periods (Cantarero et al., 2014;Cardellach et al., 2003) and is also responsible for the present-day circulation of hotwater springs up to 70 • C (Carmona et al., 2000;Fernàndez and Banda, 1990). In both cases, fluids would have been topographically driven from elevated areas to great depths (Cantarero et al., 2014), where they circulated through basement rocks, acquiring high 87 Sr/ 86 Sr ratios (> 0.712) and high temperatures (Cardellach et al., 2003) before ascending through fault discontinuities. However, in the Penedès basin, which corresponds to the southwestern termination of the Neogene structure in the CCR, basement lithologies do not crop out and the extensional faults only involve Neogene deposits filling the basin and a Mesozoic sedimentary substrate. In this location, the main fault system acted as a conduit for several episodes of meteoric fluids percolation during the Neogene extension, and evidence of hydrothermal fluid circulation has not been reported in the area (Travé and Calvet, 2001;Travé et al., 1998;Baqués et al., 2010Baqués et al., , 2012. This fact agrees with previous studies that highlight the fact that hydrothermal activity, in particular the occurrence of hot-water springs in the Pyrenees and in the CCR, is preferably concentrated in basement rocks, which constitute the elevated footwall of the main extensional fault systems (Taillefer et al., 2017;Carmona et al., 2000).
All these observations indicate an open fluid system in the NE part of the Iberian Peninsula associated with the Neogene extensional deformation. This extensional fault system acted as a conduit for the circulation of hot fluids in Neogene times and does so in the present. This fault-controlled fluid flow could have been continuous through time or could be related to intermittent pulses. Fault control on upflowing of hot fluids along fault systems is a common process in different geological settings and has been reported in the Great Basin, USA (Faulds et al., 2010;Nelson et al., 2009), in western Turkey (Faulds et al., 2010), in the southern Canadian Cordillera (Grasby and Hutcheon, 2001), and in southern Tuscany, Italy (Liotta et al., 2010).

Conclusions
This study assesses the influence of basement rocks on the fluid chemistry during deformation in the Pyrenees and provides insights into the fluid regime in the NE part of the Iberian Peninsula. Our data indicate that regardless of the fluid origin and the tectonic context, the fluids that have interacted with basement rocks have a higher 87 Sr/ 86 Sr ratio (> 0.710) than those that have circulated through the sedimentary cover (< 0.710). On the other hand, extensional deformation structures in both the eastern Pyrenees and the northeastern part of the Catalan Coastal Range acted as conduits for hydrothermal fluid migration in Neogene times and do so in the present. These fluids likely interacted with basement rocks before ascending through fault zones and related structures.
The studied thrust resulted from a multistage late Paleozoic (Variscan) to Neogene tectonic evolution. In the study area, it places a Devonian pre-Variscan unit against a Stephano-Permian late to post-Variscan sequence, and therefore the structures present within the thrust zone, affecting both sequences, are attributed to the Alpine compression (Late Cretaceous to Oligocene) and subsequent Neogene extension. During the Alpine compression, the reactivation of the thrust resulted in the transposition of the Variscan regional foliation within the thrust zone and in the formation of a subsidiary thrust zone affecting the andesites in the footwall. Geochemical analyses indicate that meteoric fluids circulated through the thrust during this period at temperatures between 50 and 100 • C. These fluids progressively increased the extent of fluid-rock interaction from the thrust plane towards the hanging wall, as suggested by the hostrock-buffered composition of the calcite in the hanging wall deformation structures. During the Neogene extension, the Estamariu thrust was likely reactivated, and normal faults and shear fractures were formed. Geochemical data reveal that basement-derived fluids circulated through these structures at temperatures up to 210 • C. Finally, during the latest to post-stages of extension and uplift of the structure, the fluid regime changed to percolation of meteoric fluids that were buffered by the volcanic host rocks. Data availability. All data are available in Tables 1 and 2 and in the Supplement.
Author contributions. DML, GA, DC, IC, and AT were responsible for conceptualization. DML handled data curation. DML, GA, DC, IC, CMJ, and AT conducted the formal analysis, carried out the investigation, contributed to the methodology, and participated in review and editing. AT was responsible for funding acquisition. DML was responsible for writing and original draft preparation.
Competing interests. The authors declare that they have no conflict of interest.
Special issue statement. This article is part of the special issue "Faults, fractures, and fluid flow in the shallow crust". It is not associated with a conference.
Acknowledgements. We acknowledge the constructive comments from Owen Callahan and Brice Lacroix, who helped to improve the quality of the paper, as well as the editorial guidance of Randolph Williams. We also thank Josep Maria Casas from the Universitat de Barcelona for showing us the location of the studied outcrops. Carbon and oxygen isotopic analyses were carried out at the Centre Científics i Tecnològics of the Universitat de Barcelona. Strontium and neodymium analyses were performed at the CAI de Geocronología y Geoquímica Isotópica of the Universidad Complutense de Madrid. The elemental composition was analyzed at the geochemistry facility labGEOTOP of Geosciences Barcelona (GEO3BCN-CSIC). Clumped isotope thermometry was carried out at the Imperial College London. This research was carried out within the framework of the DGICYT Spanish project PGC2018-093903-B-C22 (Ministerio de Ciencia, Innovación y Universidades/Agencia Estatal de Investigación/Fondo Europeo de Desarrollo Regional, Unión Europea) and the Grup Consolidat de Recerca "Geologia Sedimentària" (2017-SGR-824). The PhD research of DML is supported by the FPI2016 (BES-2016-077214) Spanish program from MINECO.
Review statement. This paper was edited by Randolph Williams and reviewed by Brice Lacroix and Owen Callahan.