<<

Fluid systems in foreland -and-thrust belts : on overview from the Southern Pyrennees A. Travé, Pierre Labaume, J. Vergès

To cite this version:

A. Travé, Pierre Labaume, J. Vergès. Fluid systems in foreland fold-and-thrust belts : on overview from the Southern Pyrennees. O. Lacombe; F. Roure; J. Lavé; J. Vergés. Thrust belts and foreland basins, from fold kinematics to hydrocarbon systems, Springer, pp.93-115, 2007, ￿10.1007/978-3-540- 69426-7_5￿. ￿hal-00408001￿

HAL Id: hal-00408001 https://hal.archives-ouvertes.fr/hal-00408001 Submitted on 6 May 2020

HAL is a multi-disciplinary open access L’archive ouverte pluridisciplinaire HAL, est archive for the deposit and dissemination of sci- destinée au dépôt et à la diffusion de documents entific research documents, whether they are pub- scientifiques de niveau recherche, publiés ou non, lished or not. The documents may come from émanant des établissements d’enseignement et de teaching and research institutions in France or recherche français ou étrangers, des laboratoires abroad, or from public or private research centers. publics ou privés.

Distributed under a Creative Commons Attribution| 4.0 International License Fluid Systems in Foreland Fold-and-Thrust Belts: An Overview from the Southern

Travé, A. · Labaume, P. · Vergés, J.

Abstract. The analysis of three different regions of the South- The progressive increase of the 87Sr/86Sr ratio through time is Pyrenean fold-and-thrust belt reveals that during the due to the progressive uplift, exposure and erosion of the in- compression the hydrological system was compartmental- ternal Pyrenean Axial Zone. ised in time and space. During the early-middle , when the thrust front affected marine soft- in the Ainsa basin, the thrust zones were dominated by formation 1 Introduction fluids derived from Eocene marine waters trapped in the un- derlying Eocene marls, although influences of meteoric wa- Geofluids are related to the formation and deforma- ters were also present. During the middle-late Eocene, when tion of the fold-and-thrust belts and sedimentary fore- the thrust front emplaced marine rocks over continental red- land basins. The geodynamic evolution of the belt - beds in the eastern Catalan basin (L’Escala thrust), the thrust system moves geofluids and, in their fault zones were dominated by meteoric fluids. These fluids turn, geofluids condition the geodynamic evolution flowed preferentially along these faults, draining laterally the of the system by controlling the propagation of frac- meteoric fluids and acting as barriers hindering their flowing tures and faults and the dynamics of fault activity. The towards more external parts of the belt. During the Oligo- fluid sources are likely to change during this coupled cene, the most external part of the fold-and-thrust belt in the evolution, in particular because most foreland basins eastern Catalan basin developed on top of a salt detachment evolve from marine and largely underfilled to conti- horizon. The thrust front affected continental materials of late nenta overfilled and largely bypassing (e.g. Allen et al., Eocene- age. At this moment, the thrusts were con- 1986) (Fig. 1). Fractures and faults may have a major duits for meteoric fluids arriving from the surface and also for control on fluid distribution and migration pathways, evolved meteoric fluids migrating over short distance up- acting as barriers or as conduits. The etailed petrolo- wards after being in contact with the underlying evaporitic gy of the filling fractures, coupled with their dis- beds. tribution, geometry and microstructural study, is a Most of the fractures show a similar sequence of micro- key to understanding the coupling between fractures, fractures. Microfractures of stage 1 formed when the sedi- faults and fluid, in particular 1) the mechanisms and ment was poorly lithified. Microfractures of stage 2 represent kinematics of fluid-deformation relationships, 2) the the main episode of vein formation which developed when a fluid flow events and their relative timing, and 3) the progressive induration of the host occurred. During synchronism, or not, between fracturation and fluid microfracture stage 3, formed in an extensional regime, the flow. Thegeochemical composition of the -fill- host sediment was more indurated. The repetition of this se- ing minerals and their comparison with the geochem- quence of microfractures in different fracture generations of istry and mineralogy of their host rocks allow to assess the same outcrop indicates that the sediment induration was 1) the physico-chemical characteristics of the fluid, restricted to the vicinity of the vein. Away from the veins, the and therefore the fluid type and origin, 2) the degree sediment remained poorly lithified during the entire defor- of fluid/rock interaction (high versus low), 3) the type mation cycle. Calcite cement within the host rock precipitat- of fluid flow regime (closed versus open), and 4) the ed later than the syn-compressive veins, when the sediment fluid pathways (pervasive versus channelized). Finally, was more indurated. the integration of all these data within a specific geo- Elemental geochemistry and stable isotopes of the calcite logical setting constrains the relative timing of fluid veins indicates that early generation of microfractures is in- flow and the fluid driving force. filled by local fluids (either marine or meteoric), whereas ex- Themain fluid flow systems in foreland-orogen con- ternal fluids (meteoric or evolved meteoric) infilled the main texts are driven by (Oliver, 1986; Burkhard & compressive stage of microfractures. The hot temperature of Kerrich, 1988; Machel & Cavell, 1999) and topography these fluids (157ºC to 183ºC in the Atiart-Arro example) indi- (Garven 1985, 1989, 1995). Structural and microstruc- cates their circulation through deep parts of the thrust belt. tural study of faults related to compression indicates

1 Ge & Garven (1994) calculated that the total volume of fluid expelled during thrusting is in the order of 1% of the total volume of pore fluid that moved through the section after thrusting (i.e., by topographically-driv- en flow), and these fluids did not travel more than a few km into the foreland. Low fluxes of tectonically- expelled fluids are also indicated in the Rocky Moun- tain Foreland Basin (Machel & Cavell, 1999) and in the Variscan foreland because recrystallization and lithifi- cation of the host rocks occurred earlier than thrust- ing (Muchez et al., 1995). TheSouth-Pyrenean fold-and-thrust belt represents an ideal natural fi ld laboratory for the study of fluid migration during orogenic compression because the evolution of the foreland basin from marine under- filled to continental overfilled and the coupled thrust front propagation are particularly well recorded. The aims of this paper are to compare the type and ori- gin of the synkinematic fluids and the type of palaeo- hydrogeological regime at diff rent moments of the thrust front propagation. For this, we use data from three areas in three diff rent parts of the belt repre- sentative of the main situations of the belt front dur- ing the main tectonic events: the Atiart-Arro thrust- fold system in the Ainsa basin, and the L’Escala thrust Fig. 1. Evolution of a foreland basin from marine and largely and El Guix , both located in the eastern Cat- underfilled to continental. overfilled and largely bypassing. alan basin (Fig. 2). Data from Atiart-Arro and El Guix Post-orogenic vertical motions are mainly controlled by the in- are synthesised from detailed works already published terplay between tectonic and surficial denudation with con- separately elsewhere (Atiart-Arro: Travé et al., 1997, comitant isostatic rebound. 1998; El Guix: Travé et al., 2000), whereas data from L‘Escala are new data presented here for the first time. Details of analytical techniques and methods, as well that a high-fluid pressure highly favours thrusting as the detailed petrology in the cases of Atiart-Arro (Roure et al., 1994, 2005) although it is not a prereq- and El Guix, are not presented in this paper and the uisite (Dietrich et al., 1983; Cosgrove, 1993). Orogenic readers are referred to the previouly published papers. compression may drive fluid migration directly (Oli- We will also discuss the place of oil migration in ver, 1986), and significant fluid flow may migrate to- the fluid history from an example in the eastern Cat- wards the foreland coevally with thrusting (Ge and alan basin (Caja et al., 2006), and compare fluid flow Garven, 1992; Qing & Mountjoy, 1992). Bradbury & in the thrust front with that associated with the em- Woodwell (1987) inferred that each thrust sheet act- placement of basement thrust sheets in inner parts of ed as a separate hydrodynamic unit for expulsion of the belt (Grant et al., 1990; Banks et al., 1991; McCaig fluids into the foreland. Geochemical studies of cal- et al., 1995). cite cements in compressive fault zones indicate in some cases high fluid-rock interaction and relatively closed palaeohydrogeological regime (Dietrich et al., 2 South-Pyrenean Fold-and-Thrust Belt 1983; Budai, 1985; Kyser & Kerrich, 1990; Muchez et al., 1995; Calvet et al., 1996; Travé et al., 1997, 1998; The of Iberia and Europe pro- Ferket et al., 2004) but, in other cases calcite cements duced the formation of the Pyrenean orogen with a indicate low fluid-rock interaction and relatively open partial subduction of the Iberian lithosphere to the palaeohydrogeological regime (Burkhard & Kerrich, north (Choukroune et al., 1989; Roure et al., 1989; Mu- 1988; Marquer & Burkhard, 1992). Quantitative mod- ñoz, 1992; Beaumont et al., 2000). Convergence began elling of fluid migration in active foreland basins re- during the Campanian and most of compressional de- veals that fluid flow derived from basin sediment com- formation ceased during the early . Two large paction is subordinate to that resulting from the rising flexural foreland basins formed on both sides of the hinterland topographic relief of the orogen (Ge & Gar- growing mountain chain: the Aquitaine retroforeland ven, 1992; Bethke & Marshak, 1990; Bitzer et al., 1998). and the Ebro foreland basins (Fig. 2). After an early

2 Fig. 2. Structural sketch of the Pyrenean belt and location of the studied outcrops. A: Atiart- Arro (early-middle Eocene); B: L’Escala (middle-late Eocene); C: El Guix (Oligocene); D: outcrop zone of Armanciès Formation marls (lower Eocene); E: outcrop area of Gavarnie basement thrust discussed in the text.

period of low subsidence rate and mainly continental acter (e.g., Coney et al., 1996; García-Castellanos et al., or marine platform sedimentation, the South-Pyrene- 2003). an foreland basin developed as underfilled and ma- During Pyrenean shortening there was a gener- rine during the early-middle Eocene, with turbiditic al migration of tectonic activity towards the foreland troughs both in the eastermost and westernmost parts as well as of the sedimentary depocentres. Three main (the eastern Catalan and the Ainsa-Jaca basins, respec- periods of thrust sheet movements can be diff renti- tively). Subsequently, the basin became overfilled and ated with diff rent rates of motion (e.g. Vergés et al., continental from the middle-late Eocene to the end of 2002), namely: shortening during the late Oligocene (Puigdefàbregas et al., 1986, 1992; Vergés et al., 1995, 1998). At the mid- 1. An early period during the Late -Pa- dle-upper Eocene boundary (~37 Ma) the uplift of the leocene, related to the emplacement of the upper- Western Pyrenees triggered the end of the marine fore- most thrust sheets, with very low rates of about land basin stage and originated an intermountain ba- 0.5 mm/a; sin limited by the Pyrenees, the Catalan Coastal Rang- 2. A second period, during the early and middle Eo- es and the Iberian Range (Vergés & Burbank, 1996). cene, which is characterized by fast rates of thrust An internal fluvial network delivered sediments to the sheet emplacement in the order of 4.0 to 4.4 mm/a characterized by a large central lake (e.g., in submarine conditions; Anadón et al., 1979; Arenas & Pardo 1999). In the east- 3. During the third period (late Eocene – Oligocene) ern Catalan basin, the deposition of the Cardona salts there was a decrease of the rates of thrust emplace- (~37 Ma) coincides with this transition from marine to ment down to 1.5 to 2.6 mm/a. continental conditions. The end of deformation along the South-Pyrenean thrust front occurred during the Thesecond and third periods largely coincide with the late Oligocene as determined by magnetostratigraphy underfilled and overfilled stages of basin evolution, re- on growth strata (~24.7 Ma) (Meigs et al., 1996). Fis- spectively. sion track cooling ages in the Axial Zone of the central The three selected structural domains (Fig. 2) to Pyrenees, however, reveal that basement uplift in the conduct this geofluid study are representative of the back of the fold-and-thrust belt ended its major up- two last stages of the thrust front – foreland basin sys- lift at about 30 Ma (Fitzgerald et al., 1999), and that tem migration. The Atiart-Arro fold-and-thrust sys- younger movements occurred around 20 Ma along the tem represents the submarine stage, during the ear- southern edge of the Axial Zone (Sinclair et al., 2005; ly-middle Eocene. The L’Escala thrust represents an Campani et al., 2005). The Ebro basin infilled up to intermediate stage, during the middle-late Eocene in the late Miocene times when it was captured by Med- transitional marine-continental conditions. The El iterranean rivers and thus ending its endorheic char- Guix anticline represents the latest stage, during the Oligocene, with continental conditions (Fig. 3).

3 a

b

c

Fig. 3. Cross-sections of the studied areas. a) Atiart-Arro fold-and-thrust system (schematic section before erosion of the inner part of the Cotiella ). The box corresponds to Fig. 4. b) L’Escala thrust. c) El Guix anticline. The box corresponds to Fig. 8

4 3 Early Development cene (Ypresian), and that deformation continued dur- of Submarine Thrust Front: ing the middle Eocene (Mutti et al., 1988). Atiart-Arro Fold-and-Thrust System The study of the fluid behaviour was carried out on four hectometric-scale outcrops, correspond- ing to three thrust-fault zones and one : the 3.1 Geological Framework Sampert, Atiart, and Los Molinos thrusts, from the NE to the SW (and structurally from top to bottom), The Atiart-Arro fold-and-thrust system is located in and the Arro syncline in the footwall of the Los Moli- the eastern part of the Ainsa sub-basin, a marine dep- nos thrust (Fig. 4) (Travé et al., 1998). The Atiart and ocentre in the inner part of the central South-Pyrenean Los Molinos thrusts probably root within the Cotiella basin. The Ainsa sub-basin is mainly filled up by low- thrust, whereas the Samper thrust is a rootless “out-of- er and middle Eocene marls with sandstone and con- the-syncline” thrust. The Arro syncline is asymmet- glomerate intercalations forming turbiditic channel- ric, with a steeply-dipping eastern limb and a gently- levee complexes, and comprises westward prograding dipping western limb. outer-shelf/slope facies intercalations in the eastern The studied mesostructures aff ct mainly marly li- part (e.g. Mutti et al., 1988). Th se deposits rest above thologies. The zones consist in a-few-me- Upper Cretaceous- shelf carbonates, Permo- tres thick intervals where deformation is marked by the marls and sandstones, and a base- association of bands in marls and calcite shear ment. The northern part of the basin-fill is thrust by veins. Calcite veins are also observed over a few tens of the Cotiella Nappe, detached in the Triassic sediments metres in the thrust-fault footwalls, but they are rare and mainly consisting of Upper Cretaceous-Paleocene or absent in the hangingwalls. Theshear bands feature carbonates (Fig. 4) (Séguret, 1972). The nappe moved S-C structures related to the development of pressure- southwards for at least 20 km and was probably root- solution . The most intensely strained shear ed in the Paleozoic basement of the Axial Zone of the bands have mineral assemblages characterised by the belt, although this root cannot be precisely identifi d occurrence of dickite, which is absent in the protolith due to subsequent deformation and erosion (e.g. Mar- as well as in the less strained domains of fault zones tinez-Peña & Casas-Sainz, 2003). TheAtiart-Arro fold- (Buatier et al., 1997). Many calcite veins show a three- and-thrust system is a SW-verging imbricate structure stage evolution of their microstructure with 1) diffuse developed in the basin-fill at the front of the nappe. veinlet networks associated with sand disaggregation, Erosional surfaces and unconformities related to the 2) crack-seal shear microstructures forming most of fold-and-thrust system show that the main movement the vein volume, and 3) late extensional openings. The of the Cotiella Nappe occurred during the early Eo- geometry of stage 1 and 2 microstructures is coher- ent with the syn-thrust shearing kinematics and their

Fig. 4. Detailed cross-section of the Atiart-Arro fold-and-thrust system.

5 chronology relates to the progressive induration of Hoefs, 1976; Hudson & Anderson, 1989), this deple- the sediment during vein formation, whereas stage 3 tion most probably resulting from burial diagenesis. corresponds to a syn- or post-shear boudinage of the Calcite cements from the three microscopic stages veins forming tabular bodies more competent than the of vein development, as well as the pre-, syn- and post- surrounding marls. Thecrack-seal microstructures at- syncline veins at Arro do not show significant diff r- test the episodic nature of fault slip and calcite precipi- ences of isotopic values. On each outcrop, the oxygen tation. Similar shear veins are also present in the Arro isotopic values of the calcite cements in veins are sys- syncline, with three generations that formed before, tematically depleted between 1.3 and 2.8‰ PDB in re- during and after folding, respectively. The most-abun- lation to the host marl calcite, probably indicating a dant ones are the syn-folding veins, which formed higher temperature of the fluid from which the vein along the bedding planes in relation to bedding-slip in calcite precipitated with respect to its adjacent host the steep limb of the syncline. In the Arro syncline, as rock (Marshall, 1992), or a deeper setting of precipi- well as in the Samper and Los Molinos thrusts, calcite tation. veins contain minor amounts of celestite. Other pe- The Atiart and Los Molinos thrusts, which are trographic features of the veins are described in Travé deep-rooted thrusts, also show depletion in 13C of et al., (1997). the calcite in veins with respect to the host marl cal- Macroscopic deformation is weak outside the thrust cite. This is interpreted as an input of an external 13C- fault zones and regional cleavage is absent. However, a depleted mineralising fluid channelized through the locality in the western limb of the Arro syncline fea- thrust-fault zones. On the contrary, in the Arro syn- tures oval burrows that attest about 20% of NE-SW cline and the Samper thrust, which have no roots to- horizontal shortening prior to sediment lithification wards deeper domains, the similar d13C in the calcite (Travé et al., 1998). Sand disaggregation observed in from veins and host marls shows that the d13C of the the earliest stages of shear vein formation also shows mineralising fluid was mainly controlled by that of the that deformation aff cted initially poorly-lithifi d sed- host sediment. iments. The isotopic values of calcite from the dickite-bear- ing marls in thrust zones range from -9.0‰ to -5.8‰ PDB in d18O, and from -2.6‰ to -1.5 ‰ PDB in d13C, 3.2 δ18O and δ13C of Host Rock and Calcite Veins respectively, i.e. without significant diff rences with respect to calcite from the less deformed marls. How- The d18O and d13C values from the host marls range ever, the calcite veinlets inside these dickite-bear- from -8.2‰ to -5.8‰, and from -3.3 ‰ to -0.7‰ PDB, ing marls show higher d18O values (from -7.3‰ to - respectively (Fig. 5, Table 1). Th y are thus lower than 4.4‰ PDB in d18O, and from -2.5‰ to -1.6‰ PDB in those of the Eocene marine carbonates which range d13C), probably resulting from in situ clay-water reac- from -4‰ to +2‰ and from -0.3‰ to +2.8‰ PDB, tions, i.e. the dickite formation through partial disso- respectively (Shackleton & Kennett, 1975; Veizer & lution of illite and drainage of potassium by fluid flow

Fig. 5. d18O versus d13C of calcite cements in the veins and host- rocks.

6 Table 1. δ18O, δ13C and 87Sr/86Sr ratio of calcite cements in the veins and host-rocks. The column dif δ18O indicates the difference between de δ18O of the calcite cement in the veins and its adjacent host-rock. For Atiart-Arro and El Guix, G1, G2 and G3 refer to the three stages of microstructure development, respectively. Structure/ δ18O δ13C dif 87Sr/86Sr Structure/ δ18O δ13C dif 87Sr/86Sr microstructure δ18O microstructure δ18O pre-thrust diaclase -8.5 -2.8 1.1 0.708179 pre-

G 2 -7.8 -0.5 1.7thrust pre-thrust diaclase -8.3 -1.9 0.9 G 2 -8.1 -1.8 2.0 bed to bed sliding -8.8 -2.1 1.4 G 2 -6.9 -0.6 0.8 thrust fault -8.3 -3.1 0.9 0.708165 G 2 -7.5 -0.9 1.4 thrust fault -8.6 -2.3 1.2

Host-rock -6.1 -0.7 0.707914thrust thrust fault -8.5 -2.2 1.1 0.708171 G 2 -8.2 -0.8 1.5 bed to bed sliding -8.5 -2.5 1.1 0.708189 G 2 -8.2 -0.8 1.5 bed to bed sliding -8.7 -2.4 1.3 G 2 -7.3 -0.9 0.6 sinistral fault -9.5 -0.6 2.1 L‘Escala Arro G 2 -8.4 -0.9 1.7 Host-rock -6.5 0 Syncline Host-rock -6.7 -0.3 NW-SE diaclase -12.6 -4 5.2 Host-rock -5.8 -1 0.707892 Host-rock -7.4 -0.8 Host-rock -5.8 -1.2 NW-SE diaclase -14.4 -3.1 7.3 0.70847 Host-rock -5.8 -1.2 Host-rock -7.1 -0.7 0.708952 post-thrust G 2 -8.1 -0.9 2.3 0.707862 sinistral fault -9.2 -0.5 2.5 0.708572 G 3 -7.6 -1.1 2.6 0.707912 Host-rock -6.7 -0.3 G 1 -7.3 -0.4 1.5 0.707886 sinistral fault -9.3 -0.6 2.6 G 2 -7.4 -0.9 1.4 0.707916 Host-rock -6.34 -5.46 - 0.70941 Host-rock -6.42 -5.35 - G 2 -9.2 -3.1 1.8 0.709173 G 2 -7.88 -5.01 1.5 0.70911 Host-rock -7.4 -2.4 G 2 -7.72 -4.97 1.3 0.70911 G 2 -9.6 -3.3 2.2 G 2 -8.41 -5.01 2.0 Host-rock -6.4 -2.4 G 2 -8.56 -5.16 1.5 Host-rock -6.9 -2.1 G 2 -8.49 -5.3 1.4 less deformed area

Los Molinos Host-rock -6.7 -1.9 G 2 -8.44 -5.24 1.4 Atiart-Arro G 2 -9.5 -3.3 2.3 Host-rock -7.07-5.66 - Host-rock -7.2 -2.4 G 2 -8.49 -4.57 3.8 0.70887 Host-rock -6.5 -1.9 G 2 -8.61 -4.61 3.9 fore- Host-rock -6.4 -1.5 thrust Host-rock -4.7 -3.52 - 0.70865 G 2 -5.9 -1.8 0.5 0.708153 G 1 -8.5 -4.36 2.5 G 2 -6.9 -2 0.5 G 2 -8.2 -4.3 2.2 0.70909 G 2 -7 -1.9 0.6 G 2 -8.42 -4.36 2.4 0.70909 G 2 -6.9 -1.9 0.5 G 2 -8.37 -4.36 2.4 0.70909 El Guix

Thrusts G 2 -7.9 -2.1 1.5 0.708247 G 2 -8.38 -4.32 2.4 0.70909

Atiart thrust Host-rock -7.1 -2 G 2 -8.34 -4.46 2.3 lower backthrust Host-rock -6.3 -1.7 G 2 -8.28 -4.47 2.3 Host-rock -6 -0.8 Host-rock -6.02 -4.72 G 2 -6.3 -1.6 0.3 G 2 -8.01 -4.58 1.9 0.70906 Host-rock -8.2 -3.3 G 3 -8.28 -4.23 2.1 Host-rock -7.8 -2.9 G 2 -8.32 -4.11 2.3 0.70906 G 2 -9.3 -2.9 1.1 Host-rock -6.16 -3.71 0.70878 G 2 -9.5 -2.5 1.3 G 2 -8.37 -4.63 2.0

G 2 -8.3 -3.1 0.5upper backthrust G 2 -8.22 -4.67 1.8

Samper G 2 -9.2 -3.1 1 Host-rock -6.42-3.82 G 2 -8.6 -2.4 0.4 G 2 -8.31 -4.17 1.9 G 2 -8 -3.1 0.2 G 2 -7.97 -4.2 1.6

up.b.foo. G 2 -8.25 -4.21 1.8

7 along the shear zones (Buatier et al., 1997). Altera- precipitation from an Eocene marine water (Travé et tion of rock-fragments and feldspars has been proved al., 1997). Th refore, the vein-forming fluid was prob- to produce an 18O increase in the fluid (Yeh & Savin, ably the interstitial water trapped in the Eocene ma- 1977; Longstaff , 1993). In the context of the Atiart- rine marls derived from the Eocene seawater, or a fluid Arro thrust system, the kinematics of dickite forma- equilibrated with the Eocene calcitic sediments. tion are inferred to have been promoted by strain and The siliciclastic detrital fraction in the marls has fluid flow (Buatier et al., 1997), rather than by temper- substantially higher 87Sr/86Sr ratios (0.71723–0.71894) ature as classically described in diagenetic or hydro- with respect to the Eocene seawater, giving rise to thermal environments. Comparing the values of the the more radiogenic 87Sr/86Sr ratios in the bulk marl four outcrops, both the host marl calcite and the cal- (0.70976–0.70987) than in the calcite fraction. This in- cite cements in veins show a progressive depletion in dicates that this detrital fraction originates from the both 18O and 13C from the Atiart thrust and Arro syn- erosion of Paleozoic rocks outcropping in the paleo- cline to Los Molinos thrust and to the Samper thrust Axial Zone. (Travé et al., 1997). The structural location and buri- In the Atiart and Los Molinos thrust-fault zones, al history of the four outcrops make it difficult to ex- the 87Sr/86Sr ratios in vein calcite range from 0.70815 plain the trend of isotopic depletion by diff rences in to 0.70927, i.e. more radiogenic than the Eocene sea- burial of the structures. A possible explanation could water values. Th se relatively high values indicate that be a higher meteoric influence in the structures locat- the vein calcite in fault zones precipitated from a flu- ed closer to the emerged part of the belt (i.e. the Sam- id with a diff rent composition than the fluid present per thrust). outside the fault zones and, consequently, that these zones acted as fluid conduits during shear deforma- tion. This more radiogenic fluid composition could 3.3 87Sr/86Sr Values of Host Rock and Veins originate from 1) interaction with the Paleozoic rocks or the Eocene marls existing in deeper settings, 2) In the Arro syncline, the calcite and celestite veins and in-situ mineralogical transformations of the Eocene the calcite fraction of the host marls, show 87Sr/86Sr ra- marls in the intensely strained sediments (e.g. the il- tio ranging between 0.70774 and 0.70795, without a lite to dickite transformation), or 3) the weathering by significant diff rence between the pre-, syn- and post- meteoric water of the Paleozoic rocks in the emerged syncline veins (Fig. 6, Table 1). This range of values is part of the belt. In the case that the source was the Eo- consistent with the 87Sr/86Sr ratio of the Eocene seawa- cene marls, either from deeper settings or from local ter (Katz et al., 1972; Burke et al., 1982; DePaolo & In- more deformed part of the thrust-fault zones, the de- gram, 1985), indicating that the source of strontium trital phyllosilicates submitted to higher temperature to precipitate the calcite and celestite in veins and the and/or pressure could have expelled their more radio- calcite cement in the host sediment was the same. The genic Sr during dissolution and mineralogical trans- later was controlled by the 87Sr/86Sr ratio of the Eo- formations. The three sources are possible and they cene seawater during the whole evolution of the Arro could have acted alone or simultaneously. syncline. The 34S values of the celestite also indicates

Fig. 6. 87Sr/86Sr ratio versus d13C values of calcite cements in the veins and host-rocks.

8 3.4 Elemental Geochemistry mogenization temperatures ranging between 157 and 183ºC and salinities around 7 wt% eq. NaCl (Travé et al., 1998), reveal the presence of a hot, saline ascending The elemental composition of the fluid from which fluid restricted to the veins, where it was mixed with calcite precipitated was determined from the elemen- the local formation water (Travé et al., 1998). Th se tal geochemistry of the calcite cements (Travé et al., two types of fluids were drained towards the thrust- 1997). The results show that the elemental composi- fault zones where they acquired a higher 87Sr/86Sr ra- tion of the fluid from which all types of calcite precip- tio, probably related to local fluid-sediment reactions itated had a Mg/Ca ratio between 0.006 and 0.126, a (the dickite formation). The d13C depletion in vein cal- Sr/Ca ratio between 0.012 and 0.39, a Mg/Sr ratio be- cite from the Atiart and Los Molinos thrusts may re- tween 0.089 and 4.938, a Ca/Fe ratio between 79 and sult from the input along the thrust zones of 13C-de- 8286, and a Mn/Ca ratio between 0.1 × 10-4 and 4 × pleted mineralising fluids, whereas the d18O depletion 10-4. Most of these values are consistent with forma- in the calcite from the structurally highest/innermost tion water and not with unmodifi d marine waters or thrust-fault zones suggests also the influence of mete- oxidizing/reducing meteoric waters. Thesamples from oric water derived from the emerged part of the belt in the Atiart thrust and Arro syncline have higher val- these structures. ues of Sr/Ca and lower values of Mg/Sr than the sam- The earliest fluid regime in the Ainsa basin was an ples from the other two outcrops, reflecting original intergranular (porous) flow of the local fluid (compac- seawater composition to a certain extent (Travé et al., tional flow) allowing for a pervasive isotopic and ele- 1997). mental exchange with the marls prior to vein forma- tion. With the onset of compressional deformation, channelized flow along tectonic slip surfaces became 3.5 Fluid Circulation in the Early Development dominant, contributing to local sediment dewatering, of the Submarine Thrust Front but also to the input of external fluids along the thrust fault zones. The geochemical data do not record sig- Syn-kinematic fluid flow in the submarine thrust front nificant evolution of fluid composition during the suc- during the early-middle Eocene is evidenced macro- cessive episodes of vein formation along the slip sur- scopically by the abundance of calcite shear veins faces at a same locality. Such a fluid regime dominated within the marls involved in the Atiart-Arro fold-and- by the dewatering of local sediments is consistent with thrust system. the submarine context of deformation that aff cted re- Microstructures indicate that deformation aff cted cently deposited, poorly lithifi d sediments. This tec- initially poorly lithifi d sediment, and the crack-seal tonic context suggests that fluids may have been over- mechanism of formation of the shear veins attests the pressured. episodic nature of fault-slip and associated fluid flow in fractures. Distribution of the veins suggests that the main source of fluid was dewatering of the poorly per- 4 Intermediate Development meable marls from the thrust footwalls, probably relat- of the Fold-and-Thrust System ed to both (i) vertical due to burial under within a Marine-Continental thrust sheets, and (ii) tectonic horizontal shortening. Transitional Environment: L’Escala Thrust Th se fluids were drained upwards towards the thrust- fault zones, in which they migrated laterally towards 4.1 Geological Framework the thrust front due to the anisotropy of the fracture permeability in these zones. The L’Escala thrust is part of the Alpens-L‘Escala an- The 87Sr/86Sr ratios of the host marl calcite and of tiformal structure within the south-eastern Pyrene- the calcite and celestite in the veins away from the an foreland basin (eastern Catalan basin) (Muñoz et thrust-fault zones indicate that the original water al., 1986). The structure is constituted by thrusts and trapped interstitially in the marls was Eocene seawa- folds aff cting the middle and upper Eocene alluvial ter. This is corroborated by the d34S in the vein celestite fan and deltaic syn-tectonic deposits of the Bellmunt from the Arro syncline, which yielded values ranging and Puigsacalm Formations, respectively (Fig. 7). The from 18.3 to 21.7‰ (CDT) (Travé et al., 1997), consis- folded foreland is limited at the north by the Vallfogo- tent with those of the Eocene seawater sulphate (Clay- na thrust and, at the south, by the Bellmunt anticline pool et al., 1980). Thed18O and d13C values of the same from the subhorizontal sediments of the Ebro basin samples reveal a change of the porewater composi- (Fig. 3b). tion from marine to formation-water during the ear- The L’Escala anticline developed in the hanging- ly burial stage. Fluid inclusion microthermometry of wall of a thrust which places middle-upper Lutetian the celestite in the Arro syncline veins, indicating ho- continental sediments of the Bellmunt Formation over

9 Fig. 7. Stratigraphic framework of the south-eastern Pyrenean foreland basin with location of the El Guix outcrop (upper square), L’Escala outcrop (two intermediate squares), and Armancies Formation (lower square). lower Bartonian marine sediments of the Puigsacalm munt Formation. Calcite cements in the thrust faults Formation (Muñoz et al., 1986). and in the bedding planes precipitated synchronous- The Bellmunt Formation is constituted by a clastic ly to thrust activity due to successive shear openings, succession evolving from fluvio-deltaic to alluvial red as indicated by the occurrence of shear bands formed deposits prograding from the Pyrenean Axial Zone to- by the crack-seal mechanism. The post-thrust sinis- wards the south margin of the Lutetian basin (Ramos tral strike-slip faults and NW-SE fractures have been et al., 2002). The limit between the Bellmunt Forma- only observed in the footwall of the lower thrust in the tion and the Puigsacalm Formation is sharp and cor- Puigsacalm Formation. The former are filled by cal- responds to a marine transgression at a basin scale. cite cement precipitated synchronously to shear-fault The Puigsacalm Formation is constituted by grey mi- opening, whereas the NW-SE fractures correspond to croconglomerates and sandstones with glauconite extensional diaclases, parallel to the compression-re- prograding from the east and forming metric-thick lated cleavage planes, and filled by calcite cement not sequences interpreted as deltaic systems. aff cted by pressure-solution related to cleavage for- In the studied outcrop, the Bellmunt Formation is mation. formed by decimetric to metric intervals of clays, silt- The pre-thrust fractures are diaclases with sharp, ites, sandstones, and conglomerates interpreted as straight or undulated walls and are filled by calcite ce- channel and flood plain deposits. ThePuigsacalm For- ments constituting veins from 1 cm to 5 cm thick and mation in the studied outcrop basically consists of cal- up to several metres long. Th se fractures are filled by cite cemented sandstones. two generations of cement: a first generation constitut- The sampled section comprises three diff rent ed by a non isopachous rim of euhedral calcite crystals thrust fault zones. The main fault places materials and a second generation of anhedral sparry calcite ce- of the Bellmunt Formation on top of the Puigsacalm ment from few millimetres up to 2 centimetres thick, Formation, and the other two faults aff ct only the with deformed mechanical twin planes and undulat- Bellmunt Formation. Three sets of cross-cut calcite- ing extinction. filled fractures have been recognised: 1) extensional Thethrust fractures are the thrust faults and the as- pre-thrust fractures (cut and deformed by the thrust sociated stratification planes which moved synchro- plane), 2) thrust-related shear-fractures (thrust faults nously to the thrust movement due to bed to bed slid- and associated bed to bed sliding), and 3) post-thrust ing and are filled by calcite cements constituting veins fractures cutting the thrust planes. The latter consists from 1 to 10 cm thick. Theveins within the thrust frac- of: 3a) sinistral strike-slip faults, and 3b) NW-SE ex- tures are constituted by calcite cement, forming diff r- tensional fractures. The pre-thrust fractures have only ent shear bands (up to 6 successive shear bands have been observed in the Bellmunt Formation. Calcite ce- been observed in a single fracture) parallels to the frac- ments within the thrust faults have been observed ture wall. Locally, the limits between the shear bands in the three thrust faults and calcite cements in the are stylolitic surfaces. Each shear band is formed by a bedding planes have been only observed in the Bell- crack-seal mechanism. The calcite crystals are brown-

10 ish-greenish, anhedral or subeuhedral, locally elon- became open for fluid circulation and for precipitation gated, variable in size from 30 µm up to 2 mm, small- of calcite cements. er close to the crack-seal border and increasing in size The d18O of the calcite cement in the post-thrust farward. The bigger crystals show abundant mechani- strike-slip faults range from -9.5 to -9.2‰ PDB and cal twinning planes and undulating extinction. the d13C range from -0.6 to -0.5‰ PDB. Th se d18O The post-thrust fractures include two sets of frac- values are on average 2-3‰ more depleted than those tures: the sinistral faults and the NW-SE fractures. The of the Puigsacalm Formation bulk host rock, whereas veins filling the sinistral faults are constituted by cal- the d13C values are within the same range as those of cite crystals, with planar or interpenetrated borders. the host rock. The similarity between the d13C values The crystals are transparent and up to 5 mm large and between the host rock and the calcite cement within show mechanical twinning planes and undulating ex- the fractures indicates host-rock buff ring. The diff r- tinction. The NW-SE fractures are diaclases. The frac- ence in d18O between the bulk host rock and the calcite ture walls are sharp and straight. The fractures are cement within the pre-thrust and thrust-related frac- filled by euhedral sparry transparent calcite crystals tures and post-thrust strike-slip faults, is probably be- up to 5 millimetre large, with planar or interpenetrat- cause d18O of the bulk host rock is a mixture between ed borders, abundant deformed mechanical twinning a lower d18O of the interparticle calcite cement and planes and undulating extinction. a d18O closer to 0‰ of the marine bioclastic compo- nents, whereas the d18O of the calcite cements within the fractures reflect only the composition of the calcite 4.2 δ18O and δ13C of Host Rock cement which has a non marine origin, either precipi- and Calcite Veins tated from meteoric or from formation waters. The d18O of the calcite cement in the post-thrust The d18O of the Puigsacalm Formation host rock NW-SE fractures range from -14.4 to -12.6‰ PDB and range from -7.4 to -6.5‰ PDB and the d13C range the d13C range from -4.0 to -3.1‰ PDB. The d18O of from -0.8 to 0.0‰ PDB (Fig. 5, Table 1). Th se values the calcite cement in the post-thrust NW-SE fractures are more depleted in d18O and within the same range is in average 6-7‰ more depleted than the Puigsacalm for d13C than the isotopic values of the Eocene marine Formation host rock, and the d13C is in average 3-4‰ carbonates (from -4 to +2‰ PDB and from -0.3 and more depleted than the host rock. The very low d18O +2.8‰ PDB, respectively). Th se depleted d18O values values suggest a relatively hot fluid. The depleted d13C of the bulk host rock are probably a mixture of deplet- values indicate that this external fluid was probably, in ed d18O values of the calcite cement precipitated in the origin, a meteoric fluid enriched in 12C probably due interparticle porosity during relative burial conditions to the ingress of isotopically light soil-CO2 (Cerling, and a more enriched d18O values of the marine com- 1984; Cerling et al., 1989). This fluid did not interact ponents. The d13C values around 0‰ probably indi- with the Puigsacalm Formation host rock. cate that the main source for C to produce the calcite cements may have been the marine bioclastic compo- nents (Discocyclines, Nummulites, echinoderms, bi- 4.3 87Sr/86Sr Values of Host Rock valvia, bryozoans,..) which form 50% of the Puigsa- and Calcite Veins calm Formation host rock. The carbonate fraction of the Bellmunt Formation The 87Sr/86Sr ratio of the Puigsacalm Formation bulk host rock has not been analysed because the texture of host rock is 0.70895. In the calcite cements, values are the rock did not allow to analyse separately the calcite 0.70818 in the pre-thrust fractures, 0.70816 to 0.70817 cement from the abundant Paleozoic, and in the thrust fault zones, 0.70819 in the bedding-par- older Tertiary carbonate clasts. allel veins, 0.70857 in the post-thrust strike-slip faults, The d18O and the d13C values of the calcite cement and 0.70847 in the post-thrust NW-SE fractures in the pre-thrust fractures range from -8.5 to -8.3‰ (Fig. 6, Table 1). PDB and from -2.8 and -1.9‰ PDB, respectively. The The 87Sr/86Sr ratio of the Bartonian sea (age of the d18O and the d13C of the calcite cement in fractures Puigsacalm Formation) was 0.7077–0.7078 (DePaolo related to bed to bed sliding and in the thrust faults & Ingram, 1985; Hess et al., 1986; Koepnick et al., are very similar, ranging from -8.8 to -8.3‰ PDB and 1985; Mead & Hodell, 1995; Palmer & Elderfi ld, 1985) from -3.1 to -2.1‰ PDB, respectively. The similarity of and, therefore, all the studied cements are non marine the d18O and d13C between the calcite cements within in origin. pre-thrust and thrust-related fractures indicates that The 87Sr/86Sr ratio of the Puigsacalm Formation both precipitated from the same fluid or from diff r- bulk host rock (0.70895) reflects the mixture of the ent fluids with similar composition, which is less prob- bioclastic Eocene marine components and the silici- able. Hence, during thrusting, the pre-thrust fractures

11 clastic fraction derived from an external source, prob- The cement in the pre-thrust and thrust-related ably the interaction of fluids with Paleozoic rocks. fractures precipitated from a fluid consistent with a The 87Sr/86Sr ratio homogeneity of calcite cements meteoric water composition. in the pre-thrust fractures, thrust fault zones and bed- The cement in the Puigsacalm Formation host rock ding-parallel veins, indicates a common origin for and in the post-thrust strike-slip faults precipitated the fluid, as already suggested by the oxygen and car- from a fluid consistent with a formation water compo- bon isotopic values. This radiogenic fluid composition sition, probably because a deeper burial location of the could originate from 1) interaction with the Paleozo- structure after the thrust emplacement. ic rocks or the Paleozoic siliciclastic components of The cement in the NW-SE fractures precipitated the Puigsacalm Formation host rock existing in deep- from a fluid consistent either with a meteoric modifi d er settings, 2) in-situ mineralogical transformations of to a formation water composition. the Paleozoic siliciclastic components of the host rock, or 3) the weathering by meteoric water of the Paleozo- ic rocks in the emerged part of the belt. 4.5 Fluid Circulation During The calcite cements in the post-thrust strike-slip the Intermediate Development of faults and in the post-thrust NW-SE fractures show the Fold-and-Thrust System within a 87Sr/86Sr values closer to the Puigsacalm Formation Marine-Continental Transition Environment host rock values than the cements in the former frac- tures indicating a higher contribution either of the Pa- The diff rent stages of vein formations, characterised leozoic rocks or of the Paleozoic siliciclastic compo- by diff rent geochemical signatures, attest the evolu- nents of the host rock. In the case of the post-thrust tion of the fluid system during the tectonic history. strike-slip faults, the similar d13C values are more in agreement with a high contribution of the Paleozo- ic siliciclastic components of the host rock, whereas 4.5.1 Prior to the Main Tectonic Shortening in case of the post-thrust NW-SE fractures the low- er d18O and lower d13C of the calcite cement better ac- No evidences of fluid circulation exist during the de- count for an origin of the 87Sr/86Sr values by miner- velopment of the pre-thrust fractures because no pre- alogical transformations of the underlying Paleozoic cipitation of cement took place. Precipitation of calcite rocks, at higher temperature. cements within these fractures did not occurred until the main tectonic shortening. 4.4 Elemental Geochemistry 4.5.2 The Main Thrust Event The calcite veins filling the pre-thrust and thrust-re- lated fractures show a similar range of values charac- Compressive deformation and thrust emplacement terised by low Fe (between 405 and 2770 ppm), low Sr occurred while sedimentation was still active in the values (from below the detection limit up to 1735 ppm basin. The veins related to the main tectonic shorten- of Sr), Mn content from 370 up to 1815 ppm, and Mg ing occur either in the thrust faults or related to bed- content between 645 and 7945 ppm. ding slip during flexural folding. The concentration of The calcite veins filling the post-thrust strike-slip veins inside and in the few tens of metres below and faults are characterised by a very high Fe content (be- above the thrust faults and their absence far from the tween 2915 and 9105 ppm), between 365 and 1260 ppm thrust faults indicates that these zones were zones of of Mn, between 885 and 5715 ppm of Mg and from be- preferential fluid flow during compressive deforma- low the detection limit up to 1515 ppm of Sr. tion. As it also occurred in the Ainsa area, the crack- Thecalcite veins filling the post-thrust NW-SE frac- seal mechanism of formation of the calcite shear veins tures are characterised by very high Sr content (from indicates synchronicity between vein formation and 655 and 3535 ppm), between 345 and 1185 ppm of Mn, fault movement within the thrust-related fractures between 400 and 3195 ppm of Fe, and from below the and attests the episodic nature of fault-slip and asso- detection limit up to 1795 ppm of Mg. ciated fluid flow. The calcite cement filling the porosity within the A meteoric fluid, with relatively low Fe, low Sr and Puigsacalm Formation host rock has a composition high Mg content, d18O ranging from -8.8 to -8.3‰ similar to the post-thrust strike-slip faults, with Mn PDB and d13C ranging from -3.1 to -2.1‰ PDB, and content between 265 and 705 ppm, Fe content between 87Sr/86Sr ranging from 0.70816 to 0.70819, circulated 2595 and 6900 ppm, Mg content between 1230 and through the thrust-related fractures as well as through 3605 ppm and from below the detection limit up to the pre-thrust fractures during thrust development. 445 ppm of Sr. This fluid had no or very low interaction with the host

12 rock. Theidentical chemical composition of the calcite 4.5.4 The Late Extensional cement filling the three main thrust faults, the bed- ding planes and the previously formed fractures indi- cates that they belonged to the same interconnected The post-thrust NW-SE extensional fractures are the fluid system and therefore, the presence of a wide- latest observed features, formed by the opening of the spread distributed homogeneous fluid. pre-existing cleavage planes. Th y thus post-date the The meteoric fluids, derived from the emerged part compressional deformation. of the belt and arriving to the thrust fault zones, flowed The calcite cements in these NW-SE fractures, with preferentially along these zones, in a drainage pattern relatively low Mg and Fe and relatively high Sr content, dominated by higher fracture permeability parallel low d18O values ranging from -14.4 to -12.6‰ PDB, to the thrust faults. Due to this lateral drainage, the d13C ranging from -4.0 to -3.1‰ PDB and 87Sr/86Sr ra- thrust fault zones probably acted as barriers hindering tio of 0.70847, are interpreted as precipitated from a their flowing towards more external parts of the belt. hot meteoric fluid. This fluid had no or very low inter- The diff rences in the 87Sr/86Sr ratio and d13C be- action with the carbonate fraction of the Puigsacalm tween the calcite cement in the pre-thrust and thrust- Formation host rock, whereas the 87Sr/86Sr ratio indi- related fractures and the Puigsacalm Formation host cates contribution of the Palaeozoic rocks. The high rock indicates that the geochemistry of cement was temperature during calcite precipitation indicates a poorly controlled by that of the host rock and, there- deep-circulating fluid. fore, that the fluid was not buff red by the host rock Thus, at this stage, the paleohydrogeological system and an open palaeohydrological system during the remained open to external meteoric fluids arriving main thrust event. from more internal and emerged parts of the belt, but with rather long circulation in relatively deep struc- tures. 4.5.3 The Post-Thrust Strike-Slip Event The extensional character of these fractures is in- terpreted as a stage of compressive relaxation as- Calcite shear veins precipitated synchronously to sociated with tectonic uplift nd relief formation. strike-slip fault activity following thrusting. The el- emental geochemistry of the calcite cements in these faults is consistent with precipitation from a fluid with 5 Late Development of Continental Thrust a formation water composition. The calcite cements and Fold system: El Guix Anticline precipitated from this fluid are characterised by very high Fe content, d18O ranging from -9.5 to -9.2‰ PDB, 5.1 Geological Framework d13C ranging from -0.6 to -0.5‰ PDB and 87Sr/86Sr of 0.70857. Th se 87Sr/86Sr ratio and d13C values are clos- The El Guix anticline corresponds to the southern- er to the Puigsacalm Formation bulk host rock values most structure of the south-eastern Pyrenean fore- than the calcite cements in all the other fractures in- land basin, detached above the Priabonian Cardona dicating the higher intensity of fluid-rock interaction. salt separating the non-deformed strata from the fold- A high interaction between the fluid and the host rock ed and thrusted strata. Compressional activity in the seems to be restricted, in the studied section, at this area represents the latest stages of the Alpine compres- episode of the geodynamic evolution, and cementa- sional tectonics (Sans & Vergés, 1995). tion of the Puigsacalm Formation host rock proba- The El Guix anticline is a double structure com- bly occurred at this time. The d18O values of the Puig- posed of two (Fig. 3c). The sampled section sacalm Formation bulk host rock probably indicate a displays two sets of linked fractures rooted at diff rent mixture between the marine contribution of the bio- depths and aff cting diff rent . Structurally, clastic components and a more depleted contribution it is located in the northern anticline in which a major of the calcite cement precipitated during burial condi- north-directed thrust (backthrust) merges with sever- tions. Thed13C values indicate that the main source for al south-directed thrusts (forethrusts) (Fig. 8). Short- carbon to produce the calcite cements may have been ening across the anticline is approximately 21%, con- the marine bioclastic components which form 50% of sisting of 16% of layer-parallel shortening prior to the host rock. The 87Sr/86Sr ratio indicates the non ma- folding and thrusting, and only 5% during folding and rine origin of the fluid trapped in the sediment, and thrusting (Sans et al., 2003). the great influence of the siliciclastic components of The deformed sedimentary pile is made up of fine- the host rock. The 87Sr/86Sr anomaly in the calcite ce- grained, fluvio-lacustrine deposits, consisting of clays, ment in the fault probably results from local clay min- sandstones and limestones of late Eocene-Oligocene eral reactions. age, overlying the evaporite sequence. Th se rocks are located approximately 300 metres above the detach-

13 Fig. 8. Detail of the sampled structures in the El Guix area. ment horizon. During thrusting, the maximum burial Less deformed areas of the anticline also present determined by vitrinite reflectance was 2 km, indicat- small fractures filled with calcite. We studied fractures ing that the folding and thrusting occurred at depth located about 20 m to the north of the backthrust ar- (Vergés et al., 1998), which is consistent with the pre- ray aff cting subhorizontal beds of grey marls and fine growth nature of the deformed rocks. sandstones, and 200 m above the other samples aff ct- The backthrust is an intensely deformed, 2 m thick ing horizontal lacustrine limestones in the hanging- zone, which has an offset of 30 m, formed by an ar- wall of the backthrust array. ray of four north-directed faults. Only two of them In the El Guix anticline, calcite veins are mostly present calcite veins, and they henceforth are refered found where the host rocks are clays and limestones, as lower and upper mineralised backthrust. The foot- and are almost absent in the fractures aff cting sand- wall of the upper mineralised backthrust presents mi- stones. The calcite veins range from millimetres to nor fractures filled with calcite veins. The sediments centimetres in thickness and show . located between the main forethrust and the back- The evolution pattern of the calcite-sealed frac- thrust array are deformed by a set of south-directed tures consists of three stages (Fig. 9) which reflect an thrusts (forethrusts) which branch in layer parallel de- evolution of the microstructures very similar to that tachments and, at a deeper level, merge with the back- observed at Atiart-Arro. In each fracture, deforma- thrust (Fig. 8). Forethrusts I to IV only cut sandstone tion started with a network of discontinuous micro- beds, whereas forethrusts V and VI cut sandstone beds fractures (stage 1), slightly oblique to the sedimentary and a lower gypsum horizon. None of these forethrusts lamination and associated with a ductile deformation is mineralised and their slip ranges from half a metre of the host sediment. Later, crack-seal shear veins, lo- to a few metres. The footwall of forethrust V presents cally aff cted by micro-stylolitic surfaces, record the a set of cross-cutting south and north-directed minor main thrust activity (stage 2). Finally, extensional faults with an offset of several centimetres. Th se mi- fractures, cross-cut at a large angle the microfractures nor faults are filled with calcite veins. of stage 2 (stage 3).

14 values and 0.4 to 0.9 ‰ for d13C. The calcite veins fill- ing small fractures in the upper mineralised back- thrust footwall have d18O values between 1.6 and 1.9‰ and d13C values 0.4‰ more negative than the values of their adjacent host rock. Thecalcite veins fill- ing the lower mineralised backthrust have d18O val- ues between 2.2 and 2.5‰ more negative and d13C values between 0.3 and 0.4 more positive than the val- ues of their adjacent host rock. The calcite veins filling minor faults in the forethrust footwall have d18O val- ues between 3.8 and 3.9‰ and d13C values 1.1‰ more negative than the values of their adjacent host rock. The calcite veins filling small fractures in the less de- formed areas have d18O values between 1.3 and 2.0‰ more negative than the values of their adjacent host rock and d13C values between 0.4 and 0.5‰ more pos- itive than their adjacent host rock (Fig. 5, Table 1). The oxygen isotopic compositions of all the calcite veins in the El Guix anticline are consistent with pre- cipitation from a meteoric or an evolved meteoric flu- id (Hudson, 1977; Marshall, 1992). Thecarbon isotopic compositions of the calcite veins with values within the same range of values as the calcite fraction of the host rocks indicate that the carbon for these veins pro- ceeded from homogenising C derived from the host rock. Fig. 9. Petrology of the veins showing the three microfracture The values of the calcite cements in microfracture stages stages 1, 2, and 3 are all within the same range, be- tween -8.6‰ and -7.7‰ PDB for d18O and between - 5.3‰ to -4.1‰ PDB for d13C, indicating a common The cement filling microfracture stage 3 also fills fluid responsible for calcite precipitation during the the vuggy porosity aff cting the calcite cement filling diff rent stages. microfracture stage 2 and the intergranular porosity of the host rock. Further petrographic features of the veins are de- 5.3 87Sr/86Sr Values of Host Rock scribed in Travé et al. (2000). and Calcite Veins

The calcite fractions of the host rock have a 87Sr/86Sr 5.2 δ18O and δ13C Values ratio ranging from 0.70865 to 0.70941 (Fig. 6, Table 1). of Host Rock and Calcite Veins The calcite veins in the diff rent fractures have a 87Sr/ 86Sr ratio varying from 0.70887 to 0.70911. The ha- Thed18O values of the host rock (lacustrine limestones lite and gypsum samples from the underlying Cardo- and carbonate fraction in the marls) range from -7.1‰ na Salt Formation have values ranging from 0.70793 to -4.7‰ PDB, and the d13C values vary between - to 0.70797, which are slightly higher than the late Eo- 5.7‰ and -3.5‰ PDB (Fig. 5, Table 1), which are con- cene-early Oligocene marine signal ranging from sistent with freshwater limestones and meteoric calcite 0.7077 to 0.7079 (Hess et al., 1986; Hess et al., 1989; cements (Veizer, 1992). Denison et al., 1993). The oxygen and carbon isotopic composition of the Although the host rocks adjacent to the diff r- calcite cements filling the diff rent fractures (from – ent fractures have a diff rent 87Sr/86Sr ratio, the cal- 8.6 to -7.7‰ PDB and from -5.3 to -4.1‰ PDB, respec- cite veins in the diverse fractures plot within a nar- tively) plot within a much narrower range of values row range of values, averaging those of the host rocks, than the isotopic composition of the host rocks adja- as has been already obseved for the stable isotope val- cent to the fractures. ues. An 87Sr/86Sr versus d13C plot gives a well defined Thecalcite veins filling the upper mineralised back- correlation of the 7 values (r=0.91), with the two ex- thrust have more negative isotopic values than the val- treme values determined by the calcite fraction of the ues of their adjacent host rock by 1.8 to 2.3‰ for d18O host rocks and the intermediate values defined by the

15 calcite veins in the fractures. This correlation indicates a pathway for a fluid which has been in contact with that the 87Sr/86Sr isotopic composition of the calcite the underlying evaporites. veins in the fractures was sourced within the adjacent Diff rences between the three microfracture stages host rocks. Influences from the underlying Cardona are not evident in the Mg/Ca, Ca/Fe and Mn/Ca mo- Salt Formation are not reflected by the 87Sr/86Sr ratio lar ratios of the precipitating fluids. However, the pres- of the calcite veins in the fractures. ence of Sr and Na only in microfracture stage 2 indi- cates that during this stage the topographycally-driven fluid flow, arrived deeper, at the underlying Cardona 5.4 Elemental Geochemistry Salt Formation. The intergranular porosity of the host rock was The elemental composition of the fluid from which filled during precipitation of cement in microfracture calcite precipitated was determined from the elemen- stage 3. tal geochemistry of the calcite cements (Travé et al., 2000). The results show that most of the mineralising fluids for the analysed calcite veins have Mg/Ca molar 5.5 Fluid Circulation During ratios lower than 0.2, which is consistent with a forma- the Late Development of tion water composition. Only three samples from the Continental Fold-and-Thrust System small fractures in the upper mineralised backthrust footwall have a Mg/Ca molar ratio higher than 0.2, In the El Guix area, meteoric fluids, with high Fe/Mn which is consistent with a meteoric water composition. and Fe/Mg ratios, without Na and Sr, enriched with The Sr/Ca molar ratio of the mineralising fluid for cal- 13C and with low 87Sr/86Sr with respect to their host cite veins in the small fractures in the upper miner- rock were widely distributed in the structure through- alised backthrust footwall, in the lower mineralised out all its evolution, within a relatively open palaeohy- backthrust, and in the minor faults in the forethrust drological system. Evolved meteoric fluids, with lower footwall is higher than 0.007, which is consistent with Fe/Mn and Fe/Mg ratios, with Na and Sr, depleted in a formation water composition. The Sr content in the 13C and with high 87Sr/86Sr with respect to their host upper mineralised backthrust and in the small frac- rock were only present during thrust faults develop- tures in the less deformed areas is always below the ment within a relatively closed palaeohydrological detection limit. Most of the mineralising fluids for cal- system. The underlying evaporites acted as the lower cite veins in the upper mineralised backthrust have a boundary of the aquifer. Ca/Fe molar ratio higher than 1000, consistent with a meteoric water composition. By contrast, most of the mineralising fluids for calcite veins in the small frac- 6 Discussion and Conclusions tures in the upper mineralised backthrust footwall, in the lower mineralised backthrust, in the minor faults In this section we summarize the existing links be- in the forethrust footwall and in the small fractures in tween fluid history and tectonic evolution in the the less deformed areas have a Ca/Fe molar ratio low- South-Pyrenean foreland basin during the Tertiary. er than 1000, consistent with a formation water com- Th se relationships show a mutual interaction, also position. Most of the mineralising fluids for calcite detected in other foreland basins, which are working veins in all the studied samples have a Mn/Ca molar together during the forwards propagation of tecton- ratio higher than 0.0002, consistent with a meteoric ic stresses through the . The under- water composition. The presence of Na is restricted to standing of the role of fluids before, during and after the upper mineralised backthrust, the small fractures thrusting as well as the role of the thrust faults as con- in the upper mineralised backthrust footwall and the duits or seals for fluid migration are key questions to lower mineralised backthrust, indicating that only the be addressed. Understanding this interaction may help backthrust array allowed precipitation from a fluid de- to clarify and better explore the migration of hydro- rived from the underlying Cardona Salt Formation. carbons in fold-and-thrust belts and foreland basins. To sum up, the fluid precipitating the calcite veins in the upper mineralised backthrust was basical- ly consistent with a meteoric origin, whereas the flu- Timing of Thrusting, Fluid Migration, Thermal His- id precipitating the calcite veins in the other fractures tory and Relationships With Hydrocarbons was basically consistent with a formation water ori- gin. However, the presence of Na in all the fractures of Fluid migration has been detected through the entire the backthrust array including the upper mineralised Tertiary evolution of the South-Pyrenean thrust front backthrust indicates that these fractures also served as as revealed by the three analysed examples (Atiart-

16 Arro, L‘Escala and El Guix) with diff rent ages of de- Diff rently, in the L’Escala site, the sequence of mi- formation (Table 2). crofractures coeval with thrusting has not been recog- In Atiart-Arro and El Guix sites, a similar sequence nised. In this outcrop, the analysed structures corre- of microstructures, formed during the main compres- spond to pre-thrusting fractures, probably developed sive event, has been recognised. Microfracture stage 1 coevally to the first shortening events characterised by is characterised by sediment disaggregation showing intense spread layer-parallel shortening (Casas et al., that vein formation began in poorly lithifi d sediment. 1996; Sans et al., 2003), syn-thrusting fractures, and Microfracture stage 2 corresponds to the main episode post-thrusting fractures. As occurred in the other two of shear vein formation indicating a change in the de- areas, the main episode of bulk host-rock cementation formation mechanism due to progressive induration by calcite occurred after the main compressive defor- of the host sediment. Microfracture stage 3 indicates a mation, i.e., synchronously to the development of the decrease of intensity of compressive deformation and post-thrust strike-slip faults. could correspond to local extensional rebound in the The sequence of calcite veins in the Atiart-Arro veins that form tabular bodies more competent than fold-and-thrust system developed under submarine the host-sediment host sediment. thrusting during early-middle Eocene times. The se- Because veins of diff rent generations in a same quence of calcite veins in the L‘Escala thrust precipi- outcrop show a similar chronology of microfractures, tated during middle-late Eocene times under transi- we conclude that sediment induration was restricted tional marine-continental conditions. And finally, the to the vicinity of the vein, probably due to cement pre- calcite veins filling the fractures at the El Guix anti- cipitation in the vein, and that the sediment away from cline precipitated during early Oligocene times un- the veins remained poorly lithifi d during the whole der continental conditions (Fig. 10). Th se continen- deformation sequence. tal conditions were established in the entire basin at In both areas, Atiart-Arro and El Guix, geochem- the Priabonian (late Eocene), when the uplift of the ical characteristics of the calcite cements in the three western Pyrenees closed up the foreland connection to microfracture stages are similar, and the main episode the Atlantic (e.g., Vergés et al., 1995; Serra-Kiel et al., of bulk host-rock cementation by calcite occurred af- 2003a, b). ter compressive deformation. This occurred because in A good correlation for the 87Sr/86Sr values of the this foreland basin the front of deformation advanced calcite veins in the Arro fold (Atiart-Arro fold-and- within the soft ediment early after its deposition. thrust system) with the Eocene marine waters is ob- served (Fig. 6). However, the samples from L‘Escala

Table 2. Evolution pattern of the structures, microstructures, fluids, hydrologic regime, and fluid type. The “*” indicate the main thrusting event. The arrows indicate chronology.

Evolution pattern Structures Fractures Fluids Hydrologic system Fluid type Atiart-Arro 3 thrusts faults 3 stages of local and deep closed (in the syn- Seawater to formation water 1 footwall syncline microfractures derived hot fluids cline) and more (in the syncline) in syn-thrust (in the syncline) open* (in the and shear veins + thrusts) Deep derived + meteoric + larger distance* local transformations* (in the thrusts) (in the thrusts) L’Escala pre-thrust diaclases 3 stages of no fluid recorded macrofractures open* Meteoric* 3 thrust faults (pre, syn and larger distance* post-thrust) post-thrust sinistral closed Evolved meteoric faults local open Hot meteoric or deeply post-thrust larger distance derived NW-SE diaclases El Guix forethrusts 3 stages of local closed Meteoric backthrusts microfractures in syn-thrust larger distance* open* Evolved meteoric* shear veins

local closed Meteoric

17 Fig. 10. Model of fluid flow pattern during the geodynamic evolution of the Southern Pyrenean fold-and-thrust belt. The different studied areas are drawn in the different sections although they did not develope at the same time or within the same section. They should be regarded as “type” outcrops.

and El Guix show a progressive more radiogenic 87Sr/ unit is a reflection of the rate of flow of formation wa- 86Sr ratio and more depleted d13C values. ter through the lithologic unit or aquifer. Relatively The 87Sr/86Sr ratios in lake water and in non ma- rapid water flow tends to make 87Sr/86Sr ratios homo- rine carbonates reflect the ages and Rb/Sr ratios of the geneous, whereas stagnant conditions caused by low rocks exposed to weathering in the drainage basin in permeability or a low hydrostatic gradient permit local which carbonates were deposited (Neat et al., 1979). variations to develop in the isotopic composition of Sr Thus, these 87Sr/86Sr ratios may increase by exposure of calcite cements (Stanley & Faure, 1979). of crystalline basement rocks to weathering (Brass, The progressive increase in the 87Sr/86Sr ratio ob- 1976), and the stratigraphic variations of the 87Sr/86Sr served in the studied samples (Fig. 6), from the Arro- ratios of non marine carbonate rocks can be used to de- Atiart thrust system to L’Escala thrust and, finally, to tect changes in the geology or hydrology of the drain- the El Guix anticline, is interpreted to have been pro- age basin (Neat et al., 1979). Similarly, the 87Sr/86Sr ra- duced by a progressive increase of exposure of crys- tios of calcite cements have been attributed to reflect talline basement rocks to weathering and progressive the isotopic composition of Sr release into the pore flu- erosion of the Axial Zone of the belt. An increase of id by diff rent Rb-bearing rocks and minerals (Stan- the basement clasts, reaching up to about 70% at the ley & Faure, 1979). On the other hand, the homogenei- top of a middle Eocene succession, has been observed ty of 87Sr/86Sr ratios of calcite cements in a lithological from the lower Lutetian to middle Bartonian succes-

18 sion of the Ripoll syncline, in the northern part of the which developed coeval to Bellmunt Formation depo- eastern Catalan basin (Ramos et al., 2002). sition, and which were buried about 3 km (Vergés et All the host rock samples of the three studied re- al., 1998). So, with a geothermal gradient of 30ºC km- gions show a great overlap of the d18O values, around - 1, that is slightly high for foreland basins, the observed 6.5‰ PDB (from -8.2 to -4.7‰) (Fig. 5). Thecalcite ce- temperatures for the fluids are higher than the ones re- ments within the fractures in the three studied regions sulting from only burial depths. also show a great overlap of the d18O values, around So, we interpret that the thrust-fault zones acted as -8.5‰ PDB (from -9.6 to -5.9‰). Systematically, the drains for the fluids derived from both the emerged calcite veins of each locality show a shift around -2 ‰ part of the belt (meteoric fluids) and deeper settings (always lower than -3.9‰) of the d18O values with re- of the thrust systems (footwall dewatering and hot as- spect to their adjacent host rock. Conversely, the host cending fluids). rock samples of the three studied regions show dis- This fault control on the fluid system was relat- tinct d13C values (from -5.7 to -3.5‰ at El Guix, from ed to the architecture of the fold-and thrust belt, it- -3.3 to -0.7‰ at Atiart-Arro and from -0.8 to 0‰ at self controlled by the organisation of the stratigraphy. L’Escala). Thed13C values of the calcite cements within Indeed, the main sole thrust of the wedge is located the fractures in the three studied regions are all com- within evaporites (salt or gypsum) horizons, either in prised from -5.3 to -0.4‰. In Atiart-Arro and L’Escala the Triassic (Ainsa), the middle Eocene (l’Escala) or regions they show a major overlap of the values, where- the upper Eocene (El Guix), according to the distri- as at El Guix the calcite cements show more depleted bution of the various depocentres (Séguret, 1972; Ver- d13C values than in the former areas (Figs. 5 and 6). In gés et al.1992). Th se evaporite detachment levels were all the cases the maximum diff rence between the cal- barriers for deep fluids migrating from the underly- cite cement in the fracture and their adjacent host rock ing basement, whereas thrusts within the wedge may is lower than 2.3‰ for d13C. have drained laterally deep fluids originating from the In the three studied regions, calcite cement within hinterland thrust basement units forming the Axial the host rock precipitated later than in the syn-com- Zone. Thus, the hot and high radiogenic fluids that ar- pressive veins, in a more evolved stage of the defor- rived at the thrust front most probably correspond to mation history, when the sediment was more indurat- originally meteoric fluids involved in topographical- ed. That means that in the earlier stage of deformation ly-driven flow between the Axial Zone relief and the the fractures permeability was higher than the host thrust front, and that could interact at depth with the rock permeability even when the porosity within Axial Zone thrust basement before migrating upward the host rock was still important. In the Atiart-Arro through the frontal thrusts. and L’Escala outcrops the d18O and d13C values of the Other than the thrust fault zones, the only strati- host rock probably is a mixture of the marine com- graphic level potentially capable to drain the flu- ponents and the calcite cements precipitated after the ids at a regional scale is the Paleocene to middle Eo- main compressive event, with more negative d18O val- cene limestone unit below the detritic basin-fill (cf. ues. The lower d13C values of the El Guix (both host Fig. 3A for Ainsa, and Fig. 7 for L’Escala and El Guix). rocks and calcites cements in veins) could be indica- Th se limestones are sealed upward by marly or evap- tive of the increase of meteoric waters involved in the orite horizons and may have drained deep fluids to- system, as also deduced from the 87Sr/86Sr ratios. This ward the foreland, preventing them to go further up- progressive depletion of the d13C values of host rock ward in the tectonic wedge. Nevertheless, in the Ainsa and calcite cements have also been observed in the area, these limestones are also present in the Cotiella Atiart-Arro region when comparing the four diff rent nappe, where they could drain meteoric fluids toward outcrops, with possible higher meteoric influences in the frontal thrusts in the detritic basin-fill (Fig. 3A). the structures located closer to the emerged parts of By contrasts, the detritic units, which constitute the the belt (Travé et al., 1998). main part of the wedge and its foreland are constitut- Hot temperature of fluids moving through the frac- ed by lenticular bodies of sandstones and conglomer- tures during thrusting seems to be coherent with two ates embedded within marly units. Th se bodies may additional results. In the Atiart-Arro example, micro- have acted as local drains, but could not compete with thermometry of fluid inclusions give temperatures be- the thrust faults for long-distance circulation within tween 157ºC and 183ºC (Travé et al., 1998), which are the wedge. too high for the constrained burial depth between a The close interaction between thrusting and fluids few hundred of metres and 3 km (Travé et al., 1998). can be used to decipher the timing of fluid migration. Temperatures between 120ºC and 145ºC have also As in many fold-and-thrust belts, the thrust system been reported in calcite veins from the lower Eo- propagates towards the foreland involving younger cene Armancies carbonates (Caja et al., 2006) locat- undeformed strata. We suggest that the thrusts are not ed in the northern flank of the Ripoll syncline (Fig. 2), only the primary conduits for fluids moving from up-

19 Factors controlling the fluid dynamics Table 3. Factors controlling the fluid Type of Type of Existence of Type of dynamics in the three studied structure microstructure high relief host rocks areas. Atiart-Arro ******** L’Escala ***-**- El Guix ********** lifted internal areas of the belt to buried external re- of the Vallfogona thrust during deposition of the up- gions of the foreland, but also that the fluids also mi- per part of the Bellmunt redbeds. Consequently, oil grate through time along the propagating thrust generation probably spanned a short period of time system, carrying hot fluids from deeper areas of the before significant thrusting took place and oil migra- thrust system. tion probably occurred soon after its generation and Diff rent structures, mesostructures and micro- on limited distance, most of the oil shows being locat- structures were embedded by fluids of diff rent origin ed in the same Armancies formation. and composition (i.e., dewatering of connate waters, The frontal thrusts studied in this paper were con- diff rent degree of interaction with meteoric waters, nected to basement thrusts in the inner part of the belt. hot ascending fluids, fluids having interacted with In the central Pyrenees, the major Gavarnie basement salts, etc.), showing an important coupling of thrust thrust deformed and uplifted the inner part of the geometry and fluid composition. Specifically, the evo- Cotiella thrust sheet from the middle Eocene. In the lution of the nature of the fluids and of their interac- Pic de Port Vieux-Plan de Larri area where the Gavar- tions with the host rocks are consistent with the evo- nie thrust places Palaeozoic strata (Silurian to Devoni- lution of the South-Pyrenean thrust front from deep an sediments) above Permo-Triassic redbeds (Fig. 2), submarine (Atiart-Arro) to mixed shallow marine - an hypersaline Sr-rich brine was pumped by fault ac- continental (L‘Escala), and eventually to intra-moun- tivity from the underlying redbeds (Grant et al., 1990; tainous continental (El Guix), coevally to the increase Banks et al., 1991; McCaig et al., 1995). The flow pat- of relief forming and Palaeozoic rock and tern was highly organised and unidirectional, paral- erosion in the Axial Zone of the belt (Table 3). The lel to the limestone in the thrust zone. Flu- type of structure/microstructure and the existence or id extraction from the footwall was over short (metric) not of a high relief seems to be the main factors con- distances, but flow along the was rapid and at trolling the fluid dynamics (Table 3), whereas the host a kilometre-scale, veins also attesting episodic upward rock type do not seem to have played a specific role. escape of fluid from the fault zone into the hanging- Evidences of oil seep or hydrocarbon bearing flu- wall (McCaig et al., 1995). By contrast, no evidences id inclusions have not been observed in the three stud- were found for a significant input of either surface or ied regions, however an interesting point in the South- metamorphic fluids during thrusting. Synchronous- ern Pyrenean belt is in the northern flank of the Ripoll ly, in the Pineta thrust complex that aff cted the Up- syncline where the lower part of the Armancies For- per Cretaceous-Paleocene cover in the hangingwall of mation is a good hydrocarbon source rock and pres- the Gavarnie thrust, the fluid flow system was dom- ents many oil shows along more than 100 km of E-W inated by two fluids, an external fluid resulting from trending outcrops. A recent study of the calcite ce- metamorphic devolatilisation reactions in the un- ments filling fractures in the Armancies Formation re- derlying silicate rocks (Permo-Triassic slates or Her- veals the presence of two main generations of cements cynian granodiorites), and a descending connate flu- (Caja et al., 2006). The first is attributed to the syn- id from the overlying carbonates, both transported at thrusting event whereas the second, which presents oil least several kilometres (Bradbury & Woodwell, 1987; inclusions, evidences hydrocarbon migration through Rye & Bradbury, 1988). the Armancies Formation coeval with this late calcite Th se works and our results thus show that thrust cement precipitation. The petrological and geochemi- fault zones are likely to focus long distance transport cal characters of this second generation of calcite ce- of fluids pumped from the surrounding rocks and dis- ments allow us to correlate them to the post-thrusting placed from the basement to the foreland basin. calcite cements filling the late extensional fracture- sat L‘Escala, and ascribed to a late shortening or post- shortening event. TheArmancies Formation went into Acknowledgements the oil under burial of about 3 km during de- position of Bellmunt redbeds (Vergés et al., 1998), and This work was carried out within the framework of was subsequently incorporated into the hanging-wall DGICYT grant CGL2006-04860. We acknowledge the

20 contribution of the Grup Consolidat „Geologia Sedi- A Core Workshop. Soc. Econ. Paleont. Miner. Core Work- mentària“ 2005/SGR00890. We thank J.A. Muñoz for shop, 7: 377–407, Tulsa. showing us the l‘Escala outcrop, the „Serveis cientifi- Burke, W.H., Denison, R.E., Hetherington, E.A., Koepink, R.B., Nelson, H.F., Orro, J.B. (1982) Variation of sea-water 87Sr/ cotècnics“ of Barcelona University for the facilities in 86Sr throughout Phanerozoic time. Geology 10: 516–519. the electron microprobe and stable isotope analyses, Burkhard, M., Kerrich, R. 1988. Fluid regimes in the deforma- the „C.A.I. de Geocronología y Geoquímica Isotópi- tion of the , Switzerland, as inferred from ca of the Universidad Complutense de Madrid“ for the stable isotope data. Contrib. Miner. Petrol, 99: 416–429. 87Sr/86Sr analyses. We greatly appreciate the construc- Caja, M.A., Permanyer, A., Marfil, R., Al-Aasm, I.S., Martín- tive comments of Rudy Swennen and François Roure, Crespo, T. (2006) Fluid flow record from fracture-fill calcite which have improved the manuscript. in the Eocene limestones from the South-Pyrenean Basin (NE Spain) and its relationship to oil shows. Journal of geo- chemical exploration 89: 27–32. Campani M., Jolivet M., Labaume P., Brunel M., Monié P., Arn- References aud N., 2005. Denudation kinematics of an orogenic prism: integrated thermochronology and tectonic study in the W- Allen, P.A., Homewood, P., Williams, G.D. (1986) Foreland Ba- Central Pyrenees (France-Spain). Earth Science Meet- sins: an introduction. In: Allen, P.A., Homewood, P. (eds.) ing « Thrust Belts and Foreland Basins », Société Géologique Foreland Basins. Special Publications of the International de France – Sociedad Geologica de Espana, 14–16 décembre Association of Sedimentologists 8: 3–12. 2005, Rueil-Malmaison. Anadón, P., Colombo, F., Esteban, M., Marzo, M., Robles, S., Calvet, F., Travé, A., Roca, E., Soler, A., Labaume, P. (1996) Santanach, P., Solé Sugrañes, L. (1979) Evolución tectonoes- Fracturación y migración de fluidos durante la evolución tratigráfica de los Catalánides. Acta Geol. Hispánica 14: tectónica neógena en el Sector Central de las Cadenas 242–270. Costero Catalanas. Geogaceta 20: 1715–1718. Arenas, C., Pardo, G. (1999) Latest Oligocene-Late Miocene la- Casas, J. M., Durney, D., Ferret, J., Muñoz, J. A. (1996) Determi- custrine systems of the north-central part of the Ebro Basin nación de la deformación finita en la vertiente sur del Piri- (Spain): sedimentary facies model and palaeogeographic neo oriental a lo largo de la transversal del río Ter. Geogace- synthesis. Palaeogeography, Palaeoclimatology, Palaeoecol- ta 20: 803–805. ogy 151: 127–148. Cerling, T.E. (1984) The stable isotopic composition of modern Banks, D.A., Davies, G.R., Yardley, B.W.D., McCaig, A.M., soil carbonate and its relationship to climate. Earth planet. Grant, N.T. (1991) The chemistry of brines from an Alpine Sci. Lett. 71: 229–240. thrust system in the Central Pyrenees: an application of flu- Cerling, T.E., Quade, J., Wang, Y. Bowman, J.R. (1989) Carbon id inclusion analysis to the study of fluid behaviour in oro- isotopes in soils and palaeosols as ecology and palaeoecolo- genesis. Geochim. Cosmochim. Acta 55: 1021–1030. gy indicators. Nature 341: 138–139. Bartrina, M.T., Cabrera, L., Jurado, M.J., Guimerà, J., Roca, E. Choukroune, P., and ECORS team (1989) The ECORS Pyrene- (1992) Evolution of the central Catalan margin of the Valen- an deep seismic profile reflection data and the overall struc- cia through (Western Mediterranean). Tectonophysics 203: ture of an . Tectonics 8: 23–39. 219–247. Claypool, G.E., Holser, W.T., Kaplan, I.R., Sakai, H., Zak, I. Beaumont, C., Muñoz, J. A., Hamilton, J., Fullsack, P. (2000) (1980) The age curves of sulfur and oxygen isotopes in ma- Factors controlling the Alpine evolution of the central Pyre- rine sulfate and their mutual interpretation. Chem. Geol. nees inferred from a comparison of observations and geody- 28: 199–260. namical models. Journal of Geophysical Research 105: 8121– Coney, P. J., Muñoz, J. A., McClay, K. R., Evenchick, C. A. (1996) 8145. Syntectonic burial and post-tectonic exhumation of the Bethke, C.M., Marshak, S. (1990) Brine migrations across southern Pyrenees foreland fold-thrust belt. Journal Geo- North America — the of groundwater. An- logical Society, London 153: 9–16. nual Review of Earth and Planetary Sciences 18: 287–315. Cosgrove, J. W. (1993) The interplay between fluids, folds and Bitzer, K., Travé, A., Carmona, J.M., Calvet, F. (1998) Fluid flow thrusts during deformation of a sedimentary succession. in foreland basins during emplacement of thrust sheets: Journal of 15: 491–500. modelling the south-Pyrenean Ainsa basin. Bull. Soc. Géol. Denison, R.E., Koepnick, R.B., Fletcher, A., Dahl, D.A., Baker, France 5: 627–634. M.C. (1993) Re-evaluation of early Oligocene, Eocene, and Bradbury, H.J., Woodwell, G.R. (1987) Ancient fluid flow with- Paleocene seawater strontium isotope ratios using outcrop in foreland terrains. In: Goff, J.C., Williams, B.P.J. (eds.). samples from the U.S. Gulf Coast. Paleoceanography 8: 101– Fluid flow in sedimentary bsins and aquifers. Geological So- 126. ciety Special Publication 34: 87–102. DePaolo, D. J., Ingram, B.L. (1985) High-resolution stratigra- 87 86 Brass, G.W. (1976) The variation of the marine Sr/ Sr ratio phy with strontium isotopes. Science 227: 938–941. during Phanerozoic time: interpretation using flux model. Dietrich, D., McKenzie, J.A., Song, H. (1983) Origin of Calcite Geochimica et Cosmochimica Acta 40: 721–730. in syntectonic veins as determined from carbon-isotope ra- Buatier, M., Travé, A., Labaume, P. Potdevin, J.L. (1997) Dickite tios. Geology 11: 547–551. related to fluid-sediment interaction and deformation in Ferket H., Swennen R., Ortuño-Arzate S., Cacas M.C., Roure F., Pyrenean thrust-fault zones. Eur. J. Mineral. 9: 875–888. 2004. Hydrofracturing in the Laramide foreland fold-and- Budai, J.M. (1985) Evidence for rapid fluid migration during thrust belt of Eastern mexico. In Swennen R., Roure F., deformation, , and Over- Granath J., eds., Deformation, fluid flow and reservoir ap- thrust Belt. In: Longman, M.W., Shanley, K.W., Lindsay, R. praisal in foreland fold-and-thrust belts, AAPG Hedberg F., Eby, D.E. (eds.). Rocky Mountain Carbonate Reservoirs – Series, Memoir, 1, 133–156.

21 Fitzgerald, P. G., Muñoz, J. A., Coney, P. J., Baldwin, S. L. (1999) Longstaffe, F.J. (1993) Meteoric water and sandstone diagenesis Asymmetric exhumation across the Pyrenean orogen: im- in the Western Canada Sedimentary Basin. In: Horbury, plications for the tectonic evolution of a collisional orogen: A.D., Robinson, A.G. (eds.). Diagenesis and basin develop- Earth And Planetary Science Letters 173: 157–170. ment. American Association of Petroleum Geologists. Stud- García-Castellanos, D., Vergés, J., Gaspar-Escribano, J., Cloet- ies in Geology, 36: 49–68. ingh, S. (2003) Interplay between tectonics, climate, and flu- Machel, H.G., Cavell, P.A. (1999) Low-flux, tectonically-in- vial transport during the evolution of the Ebro Ba- duced squeegee fluid flow (“hot flash”) into the Rocky sin (NE Iberia): Journal of Geophysical Research 108: B7, Mountain Foreland Basin. Bulletin of Canadian Petroleum 2347, doi: 10.1029/2002JB002073. Geology 47: 510–533. Garven, G. (1985) The role of regional fluid flow in the genesis Marquer, D., Burkhard, M. (1992) Fluid circulation, progressive of the Pine Point deposit, Western Canada sedimentary ba- deformation and mass-transfer processes in the upper crust: sin. Economic Geology 80: 307–324. the example of basement-cover relationships in the External Garven, G. (1989) A hydrogeologic model for the formation of Crystalline Massifs, Switzerland. Journal of Structural Ge- the giant oil sands deposits of the Western Canada Sedi- ology 14: 1047–1057. mentary Basin. Ameriacan Journal of Science 289: 105–166. Marshall, J. D. (1992) Climatic and oceanographic isotopic sig- Garven, G. (1995) Continental scale groundwater flow and geo- nals from the carbonate rock record and their preservation. logic processes. Annual Review of Earth and Planetary Sci- Geological Magazine 129: 143–160. ences 23: 89–117. Martínez-Peña, B., Casas Sainz, A. M. (2003) Cretaceous-Ter- Ge, S., Garven, G. (1992) Hydromechanical modeling of tecton- tiary tectonic of the Cotiella Basin (southern Pyr- ically driven groundwater flow with application to the Ar- enees, Spain): International Journal of Earth Sciences (Geol koma foreland basin. J. Geophys. Res. 97, B6: 9119–9144. Rundsch), 92: 99–113. Ge, S., Garven, G. (1994) A theoretical model for thrust-in- McCaig, A.M., Wayne, J.M., Marshall, J.D., Banks, D., Hender- duced deep groundwater expulsion with application to the son, I. (1995) Isotopic and fluid inclusion studies of fluid Canadian . J. Geophys. Res. 99, B7: 13851– movement along the Gavarnie thrust, Central Pyrenees: re- 13868. action fronts in carbonate mylonites. Am. J. Sci. 295: 309– Grant, N.T., Banks, D.A., McCaig, A.M., Yardley, B.W.D. (1990) 343. Chemistry, source, and behavior of fluids involved in Alpine Mead, G.A., Hodell, D.A. (1995) Controls on the 87Sr/86Sr com- thrusting of the Central Pyrenees. J. Geophys. Res. 95: 9123– position of seawater from the middle Eocene to Oligocene: 9131. Hole 689B, Maud Rise, Antarctica. Paleoceanography 10: Hess, J., Stott, L.D., Bender, M.L., Kennett, J.P., Schilling, J.-G. 327–346. (1989). The Oligocene marine microfossil record: age assess- Meigs, A. J., Vergés, J. Burbank, D. W. (1996) Ten-million-year ments using strontium isotopes. Paleoceanography 4: 655– history of a thrust sheet: Geological Society of America Bul- 679. letin 108: 1608–1625. Hess, J.; Bender, M.L.; Schilling, J.G. (1986) Evolution of the ra- Muchez, P., Slobodnik, M., Viaene, W. A., Keppens, E. (1995) tio of strontium-87 to strontium-86 in seawater from Creta- Geochemical constraints on the origin and migration of pa- ceous to present. Science 231: 979–984. laeofluids at the northern margin of the Variscan foreland, Hubbert, M. K., Rubey, W. W. (1959) Role of fluid pressure in soutern Belgium. Sedimentary Geology 96: 191–200. the mechanics of overthrust faulting. Geological Society of Muñoz, J. A. (1992) Evolution of a continental collision belt: America Bulletin 70: 115–205. ECORS-Pyrenees crustal balanced section, in McClay, K. R. Hudson, J.D. (1977) Stable isotopes and limestone lithification. (ed.) : London, Chapman & Hall, p. 235– J. Geol. Soc. London 133: 637–660. 246. Hudson, J. D., Anderson, T. F. (1989) Ocean temperatures and Muñoz, J. A., Martínez, A., Vergés, J. (1986) Thrust sequence in isotopic compositions through time. Transactions of the the Spanish eastern Pyrenees. Journal of Structural Geology Royal Society of Edinburgh 80: 183–192. 8: 399–405. Juez-Larré, J., Andriessen, P.A.M. (2006) Tectonothermal evo- Mutti, E., Séguret, M. and Sgavetti, M. (1988) Sedimentation lution of the northeastern margin of Iberian since the break- and deformation in the Tertiary sequences of the southern up of Pangea to present, revealed by low-temperature fis- Pyrenees: AAPG Mediterranean Basins Conference, p. 169. sion-track and (U-Th)/He thermochronology. A case history Neat, P. L., Faure, G., Pegram, W.J. (1979) The isotopic compo- of the Catalan Coastal Ranges. Earth and Planetary Science sition of strontium in non-marine carbonate rocks : The Letters 243: 159–180. Flagstaff ormation of Utah. Sedimentology 26: 271–282. Katz, A., Sass, E., Starinsky, A., Holland, H.D. (1972) Strontium Oliver, J. (1986) Fluids expelled tectonically from orogenic behaviour in the aragonite–calcite transformation: an ex- belts: their role in hydrocarbon migration and other geolog- perimental study at 40–98°C. Geochim. Cosmochim. Acta ic phenomena. Geology 14: 99–102. 36: 481–496. Palmer, M.R.; Elderfield, H. (1985) Sr isotope composition of Koepnick, R.B.; Denison, R.E.; Burke, W.H.; Hetherington, E. sea water over the past 75 Myr. Nature, 314: 526–528. A.; Nelson, H.F.; Otto, J.B.; Waite, L.E. (1985) Construction Puigdefàbregas, C., Muñoz, J. A. Marzo, M. (1986) Thrust belt of the seawater 87Sr/86Sr curve for the Cenozoic and Creta- development in the Eastern Pyrenees and related deposi- ceous: Supporting data. Chem. Geol., Isot. Geosci. Sect., 58: tional sequences in the southern foreland basin, in Allen, P. 55–81. A. Homewood, P. (eds.), Foreland basins, Special Public. of Kyser, T.K., Kerrich, R. (1990) Geometry of fluids in tectonical- the International Association of Sedimentologists 8: 229– ly active crustal regions. In: Nesbitt, B.E (ed.). Short Course 246. on Fluids in Tectonically Active Regimes of the Continental Puigdefàbregas, C., Muñoz, J. A., Vergés, J. (1992) Thrusting Crust. Mineralogical Association of Canada, 18: 133–230, and Foreland Basin Evolution in the Southern Pyrenees, in Vancouver.

22 McClay, K. R. (ed.), Thrust Tectonics, Chapman & Hall, p. Stanley, K.O., Faure, G. (1979) Isotopic composition and source 247–254. of strontium in sandstone cements: The High Sequence of Qing, H., Mountjoy, E.(1992) Large-scale fluid flow in the Mid- Wyoming and Nebraska. J. Sed. Petr. 49: 45–54. dle Presqu’ile barrier, Western Canada Sedimen- Travé, A., Labaume, P., Calvet, F., Soler, A. (1997) Sediment de- tary Basin. Geology 20: 903–906. watering and pore fluid migration along thrust faults in a Ramos, E., Busquets, P., Vergés, J. (2002) Interplay between lon- foreland basin inferred from isotopic and elemental geo- gitudinal fluvial and transverse alluvial fan systems and chemical analyses (Eocene south-Pyrenees, Spain). Tecto- growing thrusts in a (SE Pyrenees). In Mar- nophysics 282: 375–398. zo, M., Muñoz, J. A., Vergés, J. (eds.) Sedimentary Geology Travé, A., Labaume, P., Calvet, F., Soler, A., Tritlla, J., Buatier, on Growth Strata, Sedimentary Geology 146: 105–131. M., Potdevin, J.L., Séguret, M., Raynaud, S., Briqueu, L. Roure, F., Choukroune, P., Berastegui, X., Muñoz, J. A., Villien, (1998) Fluid migration during Eocene thrust emplacement A. M., Thereon, P., Bareyt, M., Séguret, M., Cámara, P., in the south Pyrenean foreland basin (Spain): an integrated Déramond, J. (1989) ECORS deep seismic data and balanced structural, mineralogical and geochemical approach. In: cross sections: Geometric constraints on the evolution of Mascle, A., Puigdefàbregas, C., Luterbacher, H.P., Fernàn- the Pyrenees: Tectonics 8: 41–50. dez, M. (eds.), Cenozoic Foreland Basins of Western Europe. Roure, F., Carnevali, J.O., Gou, Y. and Subieta, T. (1994) Geom- Geol. Soc. London, Spec. Publ. 134: 163–188. etry and Kinematics of the North Monagas thrust Belt (Ven- Travé, A., Calvet, F., Sans, M., Vergés, J., Thirlwall, M. (2000) ezuela). Marine and Petroleum Geology 11: 347–362. Fluid history related to the Alpine compression at the mar- Roure, F., Swennen, R., Schneider, F., Faure, J.L., Ferket, H., gin of the south-Pyrenean Foreland basin: the El Guix anti- Guilhaumou, N., Osadetz, K., Robion, P., Vandeginste, V. cline. Tectonophysics 321: 73–102. (2005) Incidence and importance of tectonics and natural Veizer, J. (1992) Depositional and diagenetic history of lime- fluid migration on reservoir evolution in foreland fold-and- stones: stable and radiogenic isotopes. In: Clauer, N., Chaud- thrust belts. Oil & Gas Science and Technology – Rev. IFP huri, S. (eds.), Isotopic Signatures and Sedimentary Records. 60: 67–106. Lecture Notes in Earth Sciences 43: 13–48. Springer, Ber- Rye, D.M., Bradbury, H.J. (1988) Fluid flow in the crust: an ex- lin. ample from a Pyrenean thrust ramp. American Journal of Veizer, J., Hoefs, J. (1976) The nature of O18/O16 and C13/C12 sec- Science 288: 197–235. ular trends in sedimentary carbonate rocks. Geochimica et Sans, M., Vergés, J. (1995) Fold development related to contrac- Cosmochimica Acta 40: 1387–1395. tional salt tectonics: southeastern Pyrenean thrust front, Vergés, J., Muñoz, J.A., Martinez, A. (1992) South Pyrenean Spain, in Jackson, M. P. A. Roberts, D. G. Snelson, S. (eds.), fold-and-thrust belt : Role of foreland evaporitic levels in AAPG Memoir 65 on Salt Tectonics: a global perspective, thrust geometry, in McClay, K.R. (ed.), Thrust Tectonics, Chapter 18: 369–378. London, Chapman et al., p. 255–264. Sans, M., Vergés, J., Gomis, E., Parés, J.M., Schiatarella, M., Vergés, J., Millán, H., Roca, E., Muñoz, J. A., Marzo, M., Cirés, Travé, A., Calvet, F., Santanach, P., Doulcet, A. (2003) Layer J., den Bezemer, T., Zoetemeijer, R., Cloetingh, S. (1995) parallel shortenning in salt-detached folds: constraint on Eastern Pyrenees and related foreland basins: Pre-, syn- and cross-. Tectonophysics 372: 85–104. post-collisional crustal-scale cross-sections. In Cloetingh, Séguret, M. (1972) Etude tectonique des nappes et séries décol- S., Durand, B., Puigdefàbregas, C. (eds.) Marine and Petro- lées de la partie centrale du versant sud des Pyrénées, Série leum Geology 12: 903–916. géologie structurale, n°2, Montpellier, Publications de l‘ Vergés, J., Burbank, D. W. (1996) Eocene-Oligocene Thrusting Université des Sciences et Techniques du Languedoc (Uste- and Basin Configuration in the Eastern and Central Pyre- la). nees (Spain), in Friend, P. F. Dabrio, C. J. (eds.), Tertiary Ba- Serra-Kiel, J., Travé, A., Mató, E., Saula, E., Ferràndez, C., Bus- sins of Spain. The Stratigraphic Record of Crustal Kinemat- quets, P., Tosquella, J., Vergés, J. (2003a) Marine and transi- ics, Cambridge University Press. World and Regional tional Middle/Upper Eocene Units of the Southeastern Py- Geology 6: 120–133. renean Foreland Basin (NE Spain). Geologica Acta 1: Vergés, J., Marzo, M., Santaeulària, T., Serra-Kiel, J., Burbank, 177–200. D. W., Muñoz, J. A., Giménez-Montsant, J. (1998) Quanti- Serra-Kiel, J., Mató, E., Saula, E., Travé, A., Ferràndez, C., Àlva- fied vertical motions and tectonic evolution of the SE Pyre- rez-Pérez, G., Busquets, P., Samsó, J.M., Tosquella, J., Fr- nean foreland basin., In Mascle, A., Puigdefàbregas, C., naquès, J., Romero, J., Barnolas, A. (2003b) An inventory of Luterbacher, H. P., Fernàndez, M. (eds.), Cenozoic Foreland the Marine and Transitional Middle/Upper Eocene Depos- Basins of Western Europe, London, Geological Society Spe- its of the Southeastern Pyrenean Foreland Basin (NE Spain). cial Publications 134: 107–134. Geologica Acta 1: 201–232. Vergés, J., Fernàndez, M. Martínez, A. (2002) The Pyrenean Shackleton, N. J., Kennett, J. P. (1975) Paleo-temperature histo- orogen: pre-, syn-, and post-collisional evolution, in Rosen- ry of the Cenozoic and the initiation of Antarctic glaciation: baum, J. G. Lister, G. S. (eds.), Reconstruction of the evolu- oxygen and carbon isotope analyses in DSDP Sites 277, 279 tion of the Alpine-Himalayan Orogen, Journal of Virtual and 281. Initial Reports of the Deep Sea Drilling Project, 29: Explorer 8: 55–84. 743–755. Yeh, H.W., Savin, S.M. (1977) Mechanisms of burial metamor- Sinclair, H.D., Gibson, M., Naylor, M., Morris, R.G. (2005) Asy- phism of argillaceous sediments: 3. O-isotope evidence. metric growth of the Pyrenees revealed through measure- Geological Society of America Bulletin 88: 1321–1330. ment and modelling of orogenic fluxes. American Journal of Science 305: 369–406.

23