# Structural style of the Languedoc Pyrenean thrust belt in relation with the inherited Mesozoic structures and with the rifting of the Gulf of Lion margin, southern France

Romain Hemelsdaël, Michel Seranne, Michel Séranne, Grégory Ballas

► To cite this version:

Romain Hemelsdaël, Michel Seranne, Michel Séranne, Grégory Ballas. Structural style of the Languedoc Pyrenean thrust belt in relation with the inherited Mesozoic structures and with the rifting of the Gulf of Lion margin, southern France. Bulletin de la Société Géologique de France, 2021, Orogen lifecycle: learnings and perspectives from Pyrenees, Western Mediterranean and analogues, 192, pp.46. 10.1051/bsgf/2021037. hal-03435967

HAL Id: hal-03435967

https://hal.science/hal-03435967

Submitted on 19 Nov 2021

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
# Structural style of the Languedoc Pyrenean thrust belt in relation with the inherited Mesozoic structures and with the rifting of the Gulf of Lion margin, southern France

Romain Hemelsdaël<fnref n="1" />,<fnref n="2" />, Michel Séranne<fnref n="1" />,*, Eglantine Husson<fnref n="3" /> and Grégory Ballas<fnref n="1" />

Received: 16 April 2021 / Accepted: 21 September 2021 / Publishing online: 18 October 2021

**Abstract** – The E-trending Pyrenean orogen results from the inversion of the mid-Cretaceous rift structures responding to an overall N-S extension, as suggested by the balanced and restored cross-sections oriented normal to the orogen. However, oblique convergence/divergence that involve strain partitioning and arcuate segments of the orogen prevent simple tectonic restorations. The Languedoc region (southern France) provides a case study of a complex polyphase deformation involving a range of reactivated structures and cross-cutting relationships, acquired in response to varying tectonic stresses with different orientations. We analyze and correlate the onshore-offshore structures of the Languedoc region, based on reassessment of existing and newly acquired subsurface data. New results in the poorly documented coastal area point to the existence of unrecognized major structures that improves onshore-offshore correlations. Our results show: (i) the part played by the Mesozoic (Early Jurassic, then mid-Cretaceous) extensional phases in the development and the localization of Pyrenean-related contractional structures; (ii) the control of inherited crustal structure on the later Oligocene rifting of the Gulf of Lion. This restoration of the Pyrenean shortening and Oligocene rifting, constructed along sections (approximately perpendicular to each other) indicates minimum shortening of 28 km and extension of 14 km, respectively, in the Languedoc foreland. Integration of the Pyrenean structural framework of Languedoc reveals a wide, NE-trending transfer zone linking the eastern Pyrénées to Provence.

**Keywords:** Languedoc / Pyrénées / oblique convergence / reactivation / structural restoration / transfer zone

**Résumé** – Style structural de la chaîne pyrénéenne en Languedoc, en relation avec les structures héritées du Mésozoïque et le rifting du Golfe du Lion (sud de la France). L'orogène pyrénéen, d'orientation générale E-W, résulte de l'inversion du rift crétacé « moyen », lui-même contrôlé par une extension de direction méridienne. Ceci est décrit dans les coupes structurales équilibrées et restaurées, perpendiculaires aux structures de l'orogène. Cependant, dans les cas de convergence ou divergence oblique et dans les cas de virgation, les restaurations tectoniques sont problématiques. Le Languedoc fournit un cas d'étude de déformations polyphasées, impliquant des réactivations ou des superpositions de structures, en réponse à des cinématiques distinctes et d'orientations différentes. Nous analysons et corrélons les structures à terre et en mer du Languedoc, sur la base de données anciennes réinterprétées et de données nouvellement acquises. Les résultats dans la zone littorale, jusqu'à présent très peu documentée, mettent en évidence de nouvelles structures majeures qui aident à la corrélation terre-mer. Nos résultats montrent: (i) le rôle des phases extensives mésozoïques (Jurassique inférieur et milieu du Crétacé) dans le développement et la localisation des structures compressives pyrénéennes et (ii) le contrôle de ces dernières sur le rifting Oligocène du Golfe du Lion. La restauration de coupes construites selon des directions choisies quasiment perpendiculaires entre elles donne $\ge$ 28 km de raccourcissement pyrénéen et $\ge$ 14 km d'extension Oligocène dans l'avant-pays languedocien. L'intégration au schéma structural pyrénéen dessine une large zone de transfert orientée NE-SW qui relie les Pyrénées franco-espagnoles à la Provence.

**Mots clés** : Languedoc / Pyrénées / convergence oblique / réactivation / restauration structurale / zone de transfert

<footnote n="1">¹ Géosciences Montpellier, Université de Montpellier-CNRS, Montpellier, France</footnote>

<footnote n="2">² GeoRessources, CNRS—Université de Lorraine, Vandœuvre-lès-Nancy, France</footnote>

<footnote n="3">³ BRGM, Orléans, France</footnote>

<footnote>*Corresponding author: michel.seranne@umontpellier.fr</footnote>
# 1 Introduction

The Pyrénées provide a much-studied example of inversion of a rift into an orogen. Geological cross-sections normal to the strike of the orogen display the compressional structures and also allow to decipher the pre-organic extensional structures that have been inverted (Gomez-Romeu et al., 2019; Jammes et al., 2009; Lescoutre and Manatschal, 2020; Saspiturry et al., 2020; Teixell et al., 2016, 2018). Inversion of rifted margins during later convergence is controlled by the inherited lithospheric structure (Espurt et al., 2019a; Jourdon et al., 2019; Lacombe and Bellahsen, 2016; Lacombe and Mouthereau, 2002) and lithosphere thermal state (Dielforder et al., 2019). The nature of the upper crust (crystalline basement versus layered sedimentary basins) and the distribution of sedimentary sequences are also critical (Duretz et al., 2020). Syntectonic sedimentary sequences preserved in upper crustal structures and in forelands provide insights into the kinematics of rift basin formation and their subsequent inversion during later contraction (Christophoul et al., 2016; Ford and Vergés, 2020; Ford et al., 2016; Lagabrielle et al., 2010; Lopez-Mir et al., 2016). These 2D analyses based on structural sections parallel to the direction of shortening and extension assume coaxial rifting and shortening (Espurt et al., 2019a; Teixell et al., 2016, 2018). However, oblique convergence/divergence with related strain partitioning, evidenced by strike-slip faults and lateral ramps, are difficult to constrain (Angrand et al., 2020; Olivet, 1996; Tavani et al., 2018; Vissers and Meijer, 2012).

After rift inversion and mountain building, the late orogenic stages of the Pyrénées are characterized by uplift and denudation (Fillon et al., 2020; Gunnell et al., 2009; Tavani et al., 2018). The present western and central Pyrénées experience N-S-directed extension characterized either by earthquake focal mechanisms and GPS measurements (Asensio et al., 2012; Nguyen et al., 2016; Rigo et al., 2015). Such extension is normal to the strike of the orogen, and the previous structures inherited from the rifting are reactivated (Fillon et al., 2020). However, the post-organic structures in the eastern Pyrénées display a different pattern, with complex interplay of reverse, strike-slip and normal faulting (Goula et al., 1999; Rigo et al., 2015), probably in relation with a distinct evolution of the deep lithosphere. Moreover, the deep root of the orogen is missing (Chevrot et al., 2018), through lithosphere thinning (e.g. Gunnell et al., 2009) or lateral delamination (e.g. Jolivet et al., 2020). Most authors agree on the link between removal of the roots beneath eastern Pyrénées and the later (Oligocene-Aquitanian) rifting of the Gulf of Lion (Chevrot et al., 2018; Diaz et al., 2018; Wehr et al., 2018). The NW-SE direction of extension is oblique to the N-S shortening in the Pyrénées, which raises questions on the possible reactivation of the E-W oriented orogenic structures during Oligocene-Aquitanian rifting. Such obliquity also requires complex step-by-step restoration in order to evaluate the direction and amount of orogenic shortening.

In addition, arcuate fold and thrust belts results from complex, non-cylindrical kinematics, involving strain partitioning, vertical axis rotation, transfer and strike-slip zones

(Macedo and Marshak, 1999; Marshak, 1988; Thomas, 1990). The rather linear Cantabric-Pyrenean orogen terminates to the east, on the shores of the NW Mediterranean (Jolivet et al., 2020). However, the fold and thrust belt is transferred northward and connects with Provence through a left-stepping accommodation zone, observed in the Corbières and Languedoc foreland (Arthaud and Séguret, 1981; Bestani et al., 2016; Lacombe and Jolivet, 2005; Séranne et al., 1995). This setting suggests a transfer zone (Frizon de Lamotte et al., 1995; Mascle et al., 1994; Tavani et al., 2018). In such areas of oblique convergence and strain partitioning, the restoration of structural sections may be an issue. The directions of successive extension, contraction and extensional collapse may be non-coaxial. Reactivation of inherited structures will depend on their more or less favorable orientations with respect to the stress regime (e.g. Martín-Gonzalez et al., 2021; Sibson, 1985; Soliva et al., 2019), especially in foreland, where inherited fault zones weaken an overall strong crust (e.g. Butler et al., 2006; Lescoutre and Manatschal, 2020).

In southern France, the Languedoc area recorded: (i) Early Jurassic Tethyan rifting and margin subsidence, (ii) mid-Cretaceous uplift of the Durancian ithmus, (iii) Pyrenean shortening and (iv) Oligocene Gulf of Lion rifting. These successive phases are characterized by specific and significantly different strain axes orientations. This results in a range of basement and cover structures with different orientations and cross-cutting settings. The area thus provides a case study of complex network of basement and cover faults, some of which have been reactivated once or more during a later stage. Deciphering the succession of deformation and associated kinematics are important for determining carbonate reservoir geometry and connectivity. Forelands display sedimentary sequences affected by folding and faulting, that provide reservoirs, covers and structural traps for oil and gas plays (e.g. Roure et al., 1994). In Languedoc, such methodology is applied to water resources exploration and management (Hemelsdaël et al., 2021; Ladouche et al., 2019; Pétré, 2020).

The main objectives of this paper are:

- to reassess the structural evolution of the northeastern foreland of the Pyrenean orogen, in a zone of oblique convergence, linking the Pyrénées and Provence, through construction of structural sections parallel to the different directions of shortening and extension, which affected the area from Early Jurassic to Miocene;

- to document multiple reactivations of a complex structural network, in response to successive non-colinear strain axes, through (i) analyses of the different stratigraphic records in footwalls and hangingwalls, and (ii) sequential restorations of the sections.

# 2 Structural and stratigraphic settings of Languedoc

The Languedoc area of southern France lies between the southeastern end of the Massif Central Variscan basement and the Mediterranean (Fig. 1). A several tens of kilometres wide strip of deformed Mesozoic and Paleogene sedimentary sequences runs parallel to the shoreline, from the eastern end of the Pyrénées to Provence.
<figure><img src="image_2.png" /></figure>

**Fig. 1.** Structural map of southern France showing the relation between the Pyrénées and Gulf of Lion margin; modified from Séranne et al. (1995). Na: Narbonne; Pe: Perpignan; Mo: Montpellier; Ma: Marseille; NPFZ: North Pyrenean Faut Zone. Offshore boreholes: AGM: Agde Maritime; AUT: Autan; CAL: Calmar; CIC: Cicindelle; MIS: Mistral; SIR: Sirocco; TRA: Tramontane. Directions of Pyrenean shortening from Gaviglio and Gonzales (1987), Lacombe et al. (1992), Arthaud and Laurent (1995) and Frizon de Lamotte et al. (2002) and Oligocene extension from Arthaud et al. (1977), Hippolyte et al. (1993) and Benedicto (1996).

Southeast of the Cévennes Fault, the area displays E-W oriented folds and thrusts, associated with latest Cretaceous to mid-Eocene synorogenic continental sedimentation, related to Pyrenean mountain-building (Arthaud and Laurent, 1995; Arthaud and Séguret, 1981). Thin-skinned, north-verging thrusting dominates deformation in the northern foreland. In the seismically poorly documented land-marine transition area, the position of the basement ramps is uncertain and still discussed (Arthaud and Séguret, 1981; Séranne et al., 1995). Pyrenean deformation is also identified offshore, across the present-day Gulf of Lion (Lacombe and Jolivet, 2005; Mauffret and Gorini, 1996; Mauffret et al., 2001; Séranne et al., 1995). In this area, the Variscan basement has been thrusted, uplifted and eroded, prior to Neogene subsidence and sedimentation (Guennoc et al., 2000). This basement is considered to have been uplifted during the early stages of mountain building (Campanian – Maastrichtian) thereby forming a major source of clastic sediments for the north Pyrenean retrowedge (Ford et al., 2016; Odlum et al., 2019; Vinciguerra, 2020) and the Provence foreland (Guieu and Roussel, 1990). The basement of the present-day Gulf of Lion was therefore thickened during Pyrenean orogeny and maintained at significant elevation during Campanian to

Bartonian, in order to provide easterly-derived detrital sources to the retro-foreland (Ternois et al., 2019). Further east, this continental block was connected with Sardinia and Corsica, the collision of which with S. Europe resulted in the Provence Pyrenean fold and thrust belt (Bestani et al., 2016; Espurt et al., 2019b; Lacombe and Jolivet, 2005).

In map view, Pyrenean-related structures are cross-cut by NE-trending normal faults (Fig. 1) (Andrieux et al., 1971). However, in the onshore part of the margin, these extensional faults root at the basement-cover interface, thus reactivating the previous Pyrenean decollement (Séranne et al., 1995). NE-trending normal faults control the formation of half-grabens in response to NW-SE extension (Fig. 1) (Benedicto et al., 1999; Husson et al., 2018; Maerten and Séranne, 1995; Sanchis and Séranne, 2000; Serrano and Hanot, 2005). The extensional decollement ramps down into the basement through the major Nîmes Fault (Fig. 2). Further southeast, offshore, all extension is accommodated by basement faults (Benedicto et al., 1996; Guennoc et al., 2000; Séranne, 1999). Offshore, extensive seismic reflection survey allowed mapping of synrift basins controlled by thick-skinned basement faults. However, synrift deposits are noticeably thin on the margin (Gorini et al., 1993, 1994; Guennoc et al., 2000), as a result of the post-orogenic
<figure><img src="image_3.png" /><figcaption>**Fig. 2.** Geological map of the study area, simplified from BRGM 1/250000 map of Montpellier (Alabouvette et al., 2001) (location on Fig. 1). The trace of the offshore normal faults is extracted from and from Guennoc et al. (2000). Location of the seismic transects, of the cross-sections and the position of the deep boreholes.</figcaption></figure>
<figure><img src="image_4.png" /><figcaption>**Fig. 3.** Synthetic lithostratigraphic succession and main tectonic stages in the Montpellier-Sète area.</figcaption></figure>

elevated topography, which did not allow preservation of synrift sedimentary formations (Séranne et al., 1995). In more distal position, the rifted margin displays highly stretched and thinned continental crust and lower crust or upper mantle denudation (Jolivet et al., 2015; Séranne, 1999).

The sedimentary sequence in the study area (Fig. 3) consists of carbonate-dominated Mesozoic formations, thickening eastwards, toward the South-East Basin (Beaudrimont and Dubois, 1977; Debrand-Passard and Courbouleix, 1984). It was deposited during the Triassic-Liassic rifting and thermal subsidence of the Tethys margin (Bonjoly et al., 1996a; Le Pichon et al., 2010). The Variscan basement is unconformably overlain by 100 to 1000 m Triassic series that mainly consist of marls, sandy clays and minor amounts of gypsum. There is no evidence of salt-related structures such as those observed in the Corbières (Ford and Vergés, 2020; Menzer et al., 2019) and Provence areas (Espurt et al., 2019b). However, in Languedoc, the Triassic shales and gypsum acted as a major decollement level at the basement-cover interface, during the successive deformation phases (Gorini et al., 1991; Maerten and Séranne, 1995; Roure et al., 1992).

During mid-Cretaceous, thermal subsidence of the Tethyan margin stopped and gave way to uplift. This inversion produced uplift, limestone karstification and erosion of the sedimentary sequences in the onshore eastern Languedoc (Bardossy and Dercourt, 1990; Marchand et al., 2020), known as the “Durancian uplift” (Gignoux, 1926; Guyonnet-Benaize et al., 2010; Masse and Philip, 1976). Erosion and weathering occurred at distinct places of the area and at successive stages, analysed in detail elsewhere (Chanvry et al., 2020; Marchand, 2019; Marchand et al., 2020). Weathering of the thick sequences of the Valanginian marly limestone under wet tropical climate was followed by reworking, transport and deposition of bauxites above the karstified Late Jurassic limestones. Typically, in the study area, bauxite deposits mark an erosional hiatus, while continental sedimentation resumed in Late Cretaceous. Duration of the hiatus increases westwards and may cover the Lias-Santonian interval (Alabouvette et al., 2001). The bauxite event (Combes, 1990; Lajoinie and Laville, 1979) was recently dated of early to mid-Albian age in the Villeveyrac bauxite mine (Fig. 2) (Marchand et al., 2020). The driving force of the mid-Cretaceous event in Languedoc is still debated (Barbarand et al., 2001; Guyonnet-Benaize et al., 2010; Tavani et al., 2018; Wyns et al., 2003). Of interest to our purpose are the following facts:

- the regional bauxite surface constitutes a mid-Cretaceous regional reference surface (Lajoinie and Laville, 1979);

- in the study area, the bauxite deposits seals extensional faults, oriented NE-SW and E-W (Marchand et al., 2020);

- this mid-Cretaceous event records both uplift onshore and subsidence offshore of the study area, respectively;

- this event is temporarily correlated with rifting in the future Pyrénées (Debroas, 1990; Ford and Vergés, 2020) and offshore Provence (Fournier et al., 2016).

Continental to paralic sedimentation resumed in localized, Late Cretaceous basins, then extended across most of the study area during Maastrichtian.
<table><caption>**Table 1.** Reflection seismic surveys and profiles used in this study.</caption><thead><tr><th>Survey</th><th>Onshore/ Offshore</th><th>Date of acquisition</th><th>Reprocessing (BRGM)</th><th>Seismic quality</th><th>Studied profiles</th></tr></thead><tbody><tr><td>DEM</td><td>Onshore</td><td>2017</td><td></td><td>+</td><td>DEM01, DEM02</td></tr><tr><td>F</td><td>Onshore</td><td>1963</td><td>2017–2018</td><td>+</td><td>F01, F02, F05, F06, F10, F11, F12</td></tr><tr><td>84SW</td><td>Offshore</td><td>1984</td><td>2017</td><td>+++</td><td>84SW-10, 84SW-11, 84SW-22, 84SW23, 84SW24, 84SW25</td></tr><tr><td>84GL</td><td>Offshore</td><td>1984</td><td>2017</td><td>++</td><td>84GL-05, 84GL-13, 84GL-22</td></tr><tr><td>81GL</td><td>Offshore</td><td>1981</td><td>2017</td><td>++</td><td>81GL-02</td></tr><tr><td>H83</td><td>Onshore</td><td>1983</td><td>2007–2008</td><td>++</td><td>H83C, H83D, H83E, H83F, H83G, H83H, H83I, H83N, H83M, H83K</td></tr><tr><td>H84</td><td>Onshore</td><td>1984</td><td>2007–2008</td><td>++</td><td>H84B, H84C, H84D, H84E, H84F, H84G, H84I, H84T, H84W</td></tr><tr><td>LRM</td><td>Offshore</td><td>1996</td><td></td><td>+++</td><td>LRM18, LRM20, LRM22, LRM24, LRM26, LRM28, LRM30, LRM32, LRM34, LRM36, LRM05, LRM07, LRM09, LRM11, LRM13, LRM15, LRM17</td></tr></tbody></table>

The Pyrenean synorogenic series comprise Late Maastrichtian to lower Eocene continental breccia and silty marls, followed by Lutetian lacustrine deposits, and finally, Bartonian syntectonic alluvial-fan breccia that distally grade to silty marls. This sequence characterizes a two-stage compressional deformation, separated by a period of relative quiescence interval (Andrieux et al., 1971; Arthaud and Séguret, 1981).

Synorogenic Pyrenean sedimentation passes upward to an intermediate fluvial succession of Priabonian age, deposited in releasing bends and oversteps of left-lateral strike-slip faults, that affect a 50 km wide wrenching zone between the Cévennes and Nîmes faults, respectively (Séranne et al., 2021).

The onshore synrift late Rupelian to Aquitanian sequences consist of continental alluvial-fans, deposited along the active normal faults and grading distally to alluvial plain and/or lacustrine deposits (Benedicto et al., 1999; Crochet, 1984). Postrift marine sediments dated to early Burdigalian to middle Miocene onlaps onto the synrift and prerift sequences (Maerten and Séranne, 1995; Magné, 1978; Séranne, 1999). The postrift sequence is truncated by the Messinian erosional surface and covered by prograding Plio-Pleistocene sequence (Clauzon, 1979; Lofi et al., 2011 and references within).

## 3 Data

The study is based on the BRGM 1/50,000 geological maps, which have been locally reassessed (e.g. Séranne et al., 2021) and detailed field mapping in the key points of the Thau-Montpellier area (e.g. Hemelsdaël et al., 2021).

Subsurface data are provided by boreholes, deeper than 100 m, extracted from the “Banque du Sous-Sol” (BRGM), geothermal and water exploration boreholes (depth ranging from several 100’s to 1000 m) from BRGM internal reports (Dörfliger and Le Strat, 2001) and oil and gas exploration boreholes, both onshore and offshore (“Minergies” database) (Figs. 1 and 2).

The seismic reflection dataset originates from different geophysical surveys of which some critical seismic profiles have been reprocessed (Tab. 1). The onshore 2D seismic data acquired in 1963, NW of the Thau lagoon and offshore 2D seismic lines acquired in 1981–1984 between Agde and Montpellier have been reprocessed (Capar and Marc, 2018a, 2018b).

Two new seismic lines (DEM1, DEM2), 11 km-long, oriented E-W and NNW-SSE, have been acquired in 2017 by CGG-Veritas. These profiles link the exploration surveys (H83–H84) in the Hérault basin to the offshore surveys in the Gulf of Lion. The details of this new seismic acquisition and processing are provided in Coppo et al. (2018).

In addition, we have integrated previous structural and stratigraphic interpretations of the onshore Gulf of Lion margin, based on seismic surveys acquired in the mid-80’s (Benedicto et al., 1996; Husson, 2013; Husson et al., 2018; Maerten and Séranne, 1995; Serrano and Hanot, 2005). The upper margin has been documented by industrial and academic seismic lines (LRM, ECORS and ESP) previously interpreted (Bache et al., 2010; Duvail et al., 2005; Mauffret et al., 2001). Despite the diversity of the seismic surveys due to different acquisition parameters and processing techniques, it is possible to interpret and correlate different seismic units (Fig. 4).

The postrift sequence is well imaged in offshore profiles. The Pliocene-Pleistocene interval progrades basinward over the Messinian erosional surface, which truncates the parallel, aggrading Miocene postrift series. The underlying synrift (Oligocene-Aquitanian) succession is imaged on both onshore and offshore profiles, characterized by divergent reflectors and evidence of growth structures. Locally, they unconformably overlie Pyrenean synorogenic sediments. The Mesozoic strata are interpreted with more uncertainties depending on the quality of the seismic surveys. The Cretaceous strata are marked by parallel to subparallel and rather continuous reflectors of moderate to high amplitude. The reflectors of the Jurassic strata are characterised by low to moderate amplitude, and they can be chaotic and semi-transparent in some places. The different units within the Jurassic succession (Malm, Dogger, marly Lias, calcareous Lias) are well identified in onshore seismic profiles (Husson, 2013; Serrano and Hanot, 2005) below the Late Cretaceous (Fig. 4A, B). In profiles DEM1 and DEM2, it was only possible to trace the high amplitude reflectors at the interface between Triassic and Jurassic units. Interpretation of the top basement surface is rather uncertain in most of the profiles.

We have thus consistently identified the following units in all available surveys, from top to bottom: (1) the Plio-Pleistocene interval, (2) the undifferentiated Miocene postrift sequence truncated by the Messinian erosional unconformity,
<figure><img src="image_5.png" /><figcaption>**Fig. 4.** Examples of seismic reflection imaging of each survey used in this study (top panels) and the interpreted stratigraphy (lower panels). The different profiles and processing are listed in Table 1.</figcaption></figure>

(3) the Oligocene-Aquitanian synrift sequence, which truncates (4) the Pyrenean syntectonic sequence of Paleocene-Eocene age and (5) the undifferentiated Jurassic interval. The underlying Triassic sediments and Variscan basement remain poorly constrained.

## 4 Structure of the top of basement

The present-day top of the Variscan basement records deformation resulting from all tectonic phases that successively affected the Languedoc area. We compiled the map of the top of basement (Fig. 5) using different published studies. Arthaud et al. (1981) compiled the top of Palaeozoic basement for the onshore Languedoc, from analyses of geological map and borehole data (oil and gas exploration). Onshore seismic data allowed to image the top of the basement: N and W of Montpellier (Benedicto, 1996), in the Hérault basin (Couëffé, 2011; Maerten and Séranne, 1995; Serrano and Hanot, 2005); in Camargue (Benedicto et al., 1996); and in the Gardon and Alès basin (Sanchis, 2000; Sanchis and Séranne, 2000). In the offshore Gulf of Lion, (Guennoc et al., 2000) compiled a structural map of the top of prerift sequence, interpreted from extensive seismic reflection data tied to exploratory offshore wells. It has been shown that, in wide areas of the Gulf of Lion, the Palaeozoic basement appears to be unconformably covered by the Oligocene-Miocene sequences related to the rifting (Arthaud et al., 1981; Cravatte et al., 1974; Guennoc et al., 1994). In these areas, the post-Palaeozoic interval thickness has been determined by depth conversion of the seismic data. The precise distribution and thickness of preserved Mesozoic units are still debated (Gorini, 1993), but there is evidence for preserved Mesozoic in the south of Camargue and in the

downfaulted grabens of the western Gulf of Lion (Guennoc et al., 2000; Mauffret and Gorini, 1996).

Gravity maps acquired by the GOCE project provide a regional, low resolution image at the scale of the western Mediterranean (Dufréchou et al., 2019). In this study, we used published Bouger and gravity residual maps (Canva, 2018; Canva et al., 2020; Guennoc et al., 2004), to localize the main density contrasts, interpreted as Variscan basement highs and lows.

The top of basement map (Fig. 5A) displays a complex structure with a range of elevations varying from more than 1.5 km (altitude of the Aigoual summit) down to depth exceeding 9 km (N. of Nîmes). Such basement depth corresponds to the western end of the South-East Basin (Beaudrimont and Dubois, 1977; Le Pichon et al., 2010). Highs and lows are generally oriented along a NE-SW trend, with higher gradients normal to this trend. Superposition of the simplified structural sketch (Fig. 5B) suggests that the present-day top of basement elevation is controlled by the Pyrenean orogeny and later rifting of the Gulf of Lion. However, several structures are unaccounted for, especially in coastal areas, where geological structures are masked by Plio-Pleistocene cover, and unexplored by seismic reflection (SE of Béziers, SW of Montpellier, W of the CIC coastal borehole). This suggests that: (i) some structures are not correctly documented, (ii) some other structures may result from earlier tectonic events, and (iii) interactions of the successive stages of deformation are not fully understood.

## 5 Pre-Pyrenean structures

The study area, located in the westernmost South-East Basin of France (Fig. 1) (Beaudrimont and Dubois, 1977;
<figure><img src="image_6.png" /><figcaption>**Fig. 5.** A. Isobath map of the top of the Paleozoic basement, compiled from various sources (see text for references). B. Same map with superimposed structural framework of Séranne et al. (1995). Main cities: Al: Alès; Be: Béziers; Mi: Millau; Mo: Montpellier; Na: Narbonne; Ni: Nîmes; Ag: Aigoual summit (1567 m). Offshore boreholes: AGM: Agde Maritime; CAL: Calmar; CIC: Cicindelle; SIR: Sirocco. Location of this map on Figure 1.</figcaption></figure>

Dubois and Delfaud, 1989) is characterised by large thickness variations within the Mesozoic successions, documented by seismic and borehole data (Alabouvette et al., 2001; Benedicto et al., 1996; Bonijoly et al., 1996a; Maerten and Séranne, 1995; Sanchis and Séranne, 2000). The thickening of the Jurassic series towards the SE, across NE-trending normal margin faults has also been described (Beaudrimont and Dubois, 1977; Bonijoly et al., 1996a, 1996b; Dubois and Delfaud, 1989; Merzeraud and Colombie, 1999). In the study area (Fig. 2), boreholes and field sections show large thickness variations over short distances, which are interpreted as evidence for differential block faulting during margin development (Fig. 6A). The marly Lias unit is 100 to 500 m thick, but reaches 1000 m in the northern area (VLF well, Fig. 2). The Dogger unit gradually thickens towards the south (about 1300 m of Dogger in GAR well) but then drastically thins in the Sète area where the Jurassic series are significantly reduced. The presence of reduced Tethyan series in the Sète area is here explained by a N-dipping normal fault (Sète F., Fig. 6A) to the north of the Sète peninsula, although we have poor control on its strike. This normal fault correlates with the horst of the Gardiole structural high (Fig. 2), also identified along-strike in the Montpellier area beneath the Neogene cover (Husson, 2013). This fault accommodated differential subsidence during the Dogger and possibly also during the Liassic times. Another feature of the Sète area is the presence of the Palaeozoic basement at shallow depth (<1500 mbsl). Deep oil and gas exploration boreholes show that the Palaeozoic basement top lies at 1583 mbsl in GAR well and abruptly plunges toward the west below the Hérault Basin (about 4500 mbsl).

The Gardiole Massif is bordered by normal faults, oriented NE-SW, dipping between 45° and 60°. Other minor faults are characterised by low dips toward the SE (15-30°), and they pre-date the mid-Cretaceous bauxitic event (Combes, 1965; Gottis et al., 1967). The lack of calcareous Liassic unit in the GAR well (Fig. 6A) is explained by the presence of a low-angle south-dipping normal fault (Arthaud and Durand, 1989; Combes, 1965). Reduced Mesozoic interval in Sète also implies a horst structure during Tethyan margin formation. The differential thickening of the Jurassic series and the presence of basement structural highs beneath Gardiole and Sète areas strongly suggest the presence of inherited horsts and sub-basins limited by NE-trending faults. These newly described features are used for the construction of the regional cross-sections.

The pre-Pyrenean structures also include mid-Cretaceous Durancian faults, which are sealed by the bauxite erosional surface (Combes, 1990; Lajoinie and Laville, 1979). Northeast of the study area, Valanginian to Hauterivian series are preserved in downfaulted blocks, whereas the uplifted blocks (bounded by EW-oriented and NE-trending faults) are eroded and karstified down to Late Jurassic limestone (Husson, 2013). Later deformation inverted or masked these faults so they are now only evidenced by the distinct stratigraphic record in the footwall and hangingwall, respectively (Fig. 6B). The area was therefore subjected to mid-Cretaceous faulting coeval with the bauxitic event, which developed in relation to the opening of the Pyrénées-South Provence mid-Cretaceous rift and the uplift of the southern Massif Central (Chanvry et al., 2020; Marchand, 2019; Marchand et al., 2020).
<figure><img src="image_7.png" /><figcaption>**Fig. 6.** A. Correlation diagram of selected boreholes along a N–S transect (location on Fig. 2) and providing evidence for Jurassic (“Tethyan”) fault-controlled subsidence. Note that the northernmost section (Sumène) is measured in the field. B. Correlation diagram of stratigraphic sections measured in the field (location on Fig. 2) and showing evidence of mid-Cretaceous (“Durancian”) faulting. BU: Buntsandstein; MU: Muschelkalk; KE: Keuper; RH: Rhetian; HE: Hettangian; SI: Sinemurian; PL: Pleinsbachien; DO: Domerian; TO: Toarcian; BAJ: Bajocian; BAT: Bathonian; CA: Callovian; OX: Oxfordian; KI: Kimmeridgian; TI: Tithonian; BE: Berriasian; VA: Valanginian; MA: Maastrichtian.</figcaption></figure>

## 6 Pyrenean thrusts in the Montpellier-Thau area

The Montpellier thrust constitutes the major onshore compressional structure related to the Pyrenean orogeny in the study area (Arthaud and Séguret, 1981). Thrust geometry is constrained by the Murviels borehole data (MUR, Fig. 2), which documents lower Liassic limestone above syntectonic Paleocene-Maastrichtian sediments, indicating a minimal offset of 4 km (Andrieux *et al.*, 1971). According to fold geometry in the hangingwall and regional paleostress reconstruction, the structure responds to a NNE-SSW shortening (Arthaud and Laurent, 1995). It is associated with

NNE-trending lateral ramps, involving a left-lateral strike-slip component of movement (Arthaud and Durand, 1989; Benedicto *et al.*, 1996). This indicates shortening oblique to the surface trace of the thrust, with a component of left-lateral strike-slip movement. Seismic reflection interpretation across the whole area allowed 3D structural modelling of the Montpellier Thrust (Husson, 2013; Husson *et al.*, 2018).

The sections in Figure 7 show the geometry at depth of the Montpellier Thrust. Section C is close to the NNE direction of transport, and it corresponds to the transect of a previous regional cross-section displaying Pyrenean shortening (Arthaud and Laurent, 1995). The NW-SE sections (A and B) are close to normal to transport direction and display the fault as lateral or oblique ramp. The Montpellier Thrust
<figure><img src="image_8.png" /><figcaption>**Fig. 7.** Structural sections in the Montpellier–Thau area. Position given in the upper right corner and in Figure 2. See text for method and comments. Note that the fault activities related to the main tectonic stages are indicated (color-coded).</figcaption></figure>

hangingwall consists of a Jurassic to Late Cretaceous
monocline and the main thrust ramps down into the basement
further south.

Tectonic-sedimentation relationships at the thrust front
allow the analysis of the chronology of Montpellier Thrust
(Arthaud and Séguret, 1981). Initial folding of the hangingwall
and deposition of syntectonic formations at the thrust front
occurred between the latest Cretaceous and Early Eocene.
Alluvial-fan breccia made of carbonate clasts in a distinctive
matrix of red silty marls commonly display microcodium:
the “Vitrolian” facies (Freytet and Plaziat, 1982). These
continental formations interfinger northwards with fan-deltas
with a marine calcareous matrix that yielded early Palaeocene
foraminifera (Combes et al., 2007). A second phase of
shortening occurred after deposition of the mid-Eocene
lacustrine limestones (classically of Lutetian age). The latter
is observed to unconformably seal the folded thrust front, and it
is also folded and thrusted within alluvial fans that rapidly pass
distally (i.e. northward) to silty marls with pedogenesis that
characterize an alluvial plain environment (Andrieux et al.,
1971). This second syntectonic formation is dated to the
Bartonian (Arthaud and Séguret, 1981). Across the Montpellier

Thrust, the amplitude of this shortening is moderate compared to
the previous shortening phase, whereas Bartonian shortening
corresponds to the major tectonic phase in the Saint-Loup thrust,
located some 20 km northward (Philip *et al.*, 1978; Séranne *et al.*,
2021).

The Thau Thrust is a newly described major structure of
the study area, which is entirely covered by Neogene
sediments and the Thau lagoon (Fig. 7B, C). Several lines
of evidence converge to confirm its existence:

- the upper Jurassic carbonates on both sides of the Thau
lagoon dip 25°–30° to the SSW, which allows to interpolate
a vertical offset in the range of 400 m. Such structural
setting could result from either a north-verging thrust or a
north-facing normal fault;

- beneath the Neogene cover, several shallow water wells
along the northern shore of the Thau lagoon ("Banque du
Sous-Sol" BRGM database) have encountered reddish
continental sandy marls, containing breccia of Jurassic
limestones, which can be correlated with either outcrops of
synrift Oligocene or syntectonic Paleocene-Eocene. Since
Oligocene is not documented in this area (Gottis et al.,
1967), we interpret this interval as synorogenic Pyrenean
sediments deposited north of an emerging Pyrenean thrust. Correlative continental detritals of early Eocene age are exposed 15 km west of Sète (Berger et al., 1981);

- the thrust documented in Valensac borehole (VAL, Fig. 2) (Arthaud and Laurent, 1995), east of the Thau lagoon, extends along-strike northeast-ward and correlates with the newly described structure. Further east, past the shoreline, this thrust connects with the Nîmes Fault (Husson et al., 2018). Symmetrically, the Thau Thrust can also be considered as a compressive splay of the major Nîmes Fault;

- beneath the Neogene cover, the Thau Thrust trace is constrained by the results of a passive seismic acquisition (Bourgeois, 2019) providing a detailed map of the top Jurassic surface in the Sète-Balaruc area. The Thau Thrust subcrop is marked by topographic anomalies of the top Jurassic at shallow depth (<200 m), consistent with the E-W orientation of the Pyrenean thrusts.

Therefore, the Sète peninsula is an isolated and uplifted block corresponding to the hangingwall of the E-W oriented Pyrenean Thau Thrust, in a complex structural zone interacting with NNE-trending left-lateral wrenching (Nîmes Fault). In addition, the Thau Thrust is superimposed onto a previous north-facing, extensional fault of Tethyan age. The lack of exposure prevents any direct dating of the Thau thrust activity. It is however likely that the Thau thrust records the same tectonic history as the Montpellier Thrust.

# 7 Rifting structures of the Gulf of Lion margin

The study area underwent NW-SE extension during the Oligocene to Aquitanian times. Detailed analyses in the onshore part of the rifted margin have documented an extensional system, consisting of NE-striking extensional faults that detached in the Variscan basement-Mesozoic cover interface, and that controlled half-graben formation (Benedicto, 1996; Husson, 2013; Maerten and Séranne, 1995; Serrano and Hanot, 2005; Séranne et al., 1995). In the offshore part, the extensional faults affect the basement and control the deposition of thin synrift sequence (Gorini et al., 1994; Guennoc et al., 2000).

The established seismic stratigraphy (Fig. 4) enables us to construct two 80 km long seismic regional cross-sections linking the onshore and offshore domains (Fig. 8). These sections are parallel to the extension direction, so that the imaged extensional structures are significant for the rifting history.

In the offshore part of the profiles, the Miocene to Plio-Quaternary succession forms a thick wedge as a result of postrift subsidence, compared with the onshore thin and discontinuous sequence (Hérault basin). The Hérault Fault corresponds to the Oligocene top-of-basement extensional decollement. It propagates landwards and emerges along the inherited, basement-involved, Cévennes Fault, of which, it reactivates the shallowest part (Maerten and Séranne, 1995).

Similar structural setting is also observed along-strike in the synrift Alès basin: the Alès Fault activates the top-of-basement decollement during Oligocene rifting. Accordingly, the Alès Fault propagates upward and reactivates the Cévennes Fault, close to the surface (Sanchis and Séranne, 2000). Therefore, the Hérault Fault bounding the synrift Hérault basin does not affect the Palaeozoic basement. On the opposite, the Nîmes Fault affects the basement and bounds the major synrift basins, which are characterised by diverging reflectors. Its vertical synrift offset decreases toward the SW, from more than 2sTWTT offset in the onshore Camargue graben (Benedicto et al., 1996) to less than 1 s TWTT and 0.2 s TWTT (Fig. 8A and B, respectively). The offset is transferred to the antithetic, NW-facing, Vistrenque Fault, thus forming a relay zone.

One of the most striking features is a structural high in the coastal zone, with basement top at about 1.5 s TWT (Fig. 8). This basement high is located between the Nîmes Fault and the Thau Thrust, thus forming a complex zone of interaction between Pyrenean thrusts and Oligocene extensional faults. The reactivated Nîmes Fault connects with the Thau Thrust system below 2.0 s TWT. As argued above, the basement offset is explained by the combined effects of the north-dipping Sète Fault, which forms a Tethyan horst, and later Pyrenean thrusting across the Thau Thrust system. Unfortunately, the poor quality of the reprocessed line F02 in the Hérault Basin (Figs. 2 and 8B) does not allow satisfactory seismic imaging of the geometry of the Thau Thrust system and associated deformed units.

A key point to address is the offshore continuity of the Mesozoic units. Only two offshore boreholes, close to the coastline, encountered the Mesozoic substratum (CAL, CIC; Fig. 2). The other deep boreholes drilled the Variscan basement directly underlying the synrift or postrift sediments, displaying Palaeozoic metasediments in AUT, MIS and TRA wells, or granitic pluton in SIR well. Based on our seismic interpretations and the offshore borehole data, we propose that Mesozoic sequences are present north and northwest of the Mistral-Sirocco high (MIS, SIR, Fig. 1). Seismic data display parallel reflectors, interpreted as Mesozoic stratified sequences, below the synrift Vistrenque graben (Fig. 8). These seismic transects also provide clues to interpret the geometry and kinematics of the Pyrenean structures affecting the prerift basement. Despite the strong obliquity of the profiles with respect to the direction of compression, some south-dipping reflectors in the offshore seismic profiles are attributed to basement thrusts, such as between the Vistrenque and Petite Camargue faults (Fig. 8A), as well as below the coastal structural high (Fig. 8A, B). It thus appears that this coastal basement high results from the Pyrenean Thau Thrust. Similarly, although not clearly imaged by seismics, we suggest that a major NNE-verging thrust system controls the uplift and erosion of the Sirocco basement high. Likewise, the Mesozoic strata would have been preserved in the thrust footwall, but there is neither borehole or seismic evidence for Mesozoic strata SE of the Petite Camargue Fault (eroded or never deposited?). Later extensional faults, well imaged on the seismic data, developed horsts atop the thrust hangingwall culmination. This structural setting is used to build the structural cross sections presented below.
<figure><img src="image_9.png" /><figcaption>**Fig. 8.** Two composite seismic reflection transects across the study area (upper panels) and interpreted sections (lower panels). The NW–SE orientation of the sections, i.e. parallel to the direction of extension, displays the rifting structures in the upper crust, and their onshore-offshore correlation. Note that vertical scales are in seconds two-way-travel-time. Position of the two transects given in Figure 2.</figcaption></figure>

# 8 Restorations of sections

## 8.1 NW-SE cross-section: The extensional structures

A first section, oriented NW-SE, runs from the Cévennes Fault to the offshore Sirocco basement horst (64 km long; Fig. 9) and is used to restore the extensional structures of the proximal Gulf of Lion rifted margin. This cross section is parallel to both onshore-offshore seismic cross-sections (Fig. 8A), parallel to the direction of extension (Séranne, 1999) and crosses the main extensional structures, thus allowing restoration of the Oligocene rifting structures. When possible, interpreted offshore seismic sections were converted to depth, using velocity stacks provided by the seismic profiles, and controlled by the few available boreholes on the line (GAR and SIR). In addition, interval velocities were determined on the neighbouring profiles that are tied to the other offshore wells. Onshore, the Jurassic stratigraphic intervals display fairly well constrained thickness variation across faults.

Offshore, between the Nîmes and Vistrenque faults, seismic reflection shows layered seismic units above the Variscan basement, interpreted as undifferentiated Jurassic succession. Dogger carbonates were also found in CAL well (Fig. 2). The prerift structure between the antithetic Nîmes and Vistrenque faults, that form a relay zone, remains poorly constrained.

On the section (Fig. 9), we distinguish the Oligocene rifting faults, the Pyrenean compressional structures and the Tethyan normal faults. The balanced cross sections were built and restored on the basis of surface conservation, and deformation of normal fault hangingwall by antithetic shearing along a 60° angle (Benedicto, 1996; Withjack and Peterson, 1993). Restoration of the Oligocene rifting structures was performed assuming a progressive increase of topography toward the south, up to 1000 m representing a Pyrenean orogenic topography (Benedicto, 1996; Séranne et al., 1995, 2021). Restoration to pre-Oligocene rifting shows about 14 km of total extension, of which, the Nîmes and the Vistrenque Faults
<figure><img src="image_10.png" /><figcaption>**Fig. 9.** Geological cross-section oriented parallel to the extension direction, in order to show the rifting structures (note identical vertical and horizontal scales). See location on Figure 2. This section is primarily derived from the seismic transect and then restored section in a post-Pyrenean state. See text for justification of topography.</figcaption></figure>

accommodates 5 km and 3 km, respectively. The prerift Nîmes Fault reveals a reverse offset of the basement-cover interface, with similar amplitude to that of the Montpellier Thrust. Considering the strong obliquity (circa 80°) of the section with respect to the direction of Pyrenean shortening, the 3 to 5 km offset measured across these structures reflects the reverse component of movement, which was associated with a significant off-section, left-lateral, strike-slip component. The Nîmes Fault therefore accommodated oblique (combined reverse and left-lateral) offset during the Pyrenean shortening.

Restoration to the prerift, post-Pyrenean state also highlights the presence of an inherited structural high between the Montpellier Thrust and the Nîmes Fault, in the Gardiole Massif. There is no evidence of Pyrenean shortening within the Gardiole Massif, which is affected by low angle normal faults of Tethyan age. We thus interpret the structural high as passively uplifted in the hangingwall of the Montpellier Thrust, during Pyrenean compression. Indeed, the Nîmes Fault, which bounds the Gardiole massif to the SE is known to be active during Tethyan margin formation, NE of the study area (Le Pichon et al., 2010; Séguret et al., 1997). The Gardiole high persisted during the Oligocene rifting, as flexural uplift in the footwall of the Nîmes Fault, which prevented preservation of Oligocene synrift deposits. Consequently, only the postrift Miocene transgressive deposits are recorded in the Gigean Basin (Figs. 2 and 9).

## 8.2 NNE-SSW cross-section — The inherited Pyrenean structures

Large scale Pyrenean structures are associated with horizontal stress σ1, oriented NNE-SSW, as derived from microtectonics markers (Arthaud and Laurent, 1995). Following

these authors, we assume a generalized NNE-SSW shortening in the study area, and consequently, chose to analyse the Pyrenean deformation along a NNE-SSE regional cross-section (77 km long, Fig. 10). This cross-section includes the Cévennes Fault to the north and the Saint-Loup, Montpellier, Thau and Agde thrusts, respectively from north to south. Since this cross-section does not intercept significant Oligocene extensional fault, the pre-Pyrenean state can be restored without interference with rifting deformation. The Pyrenean thrusts, mid-Cretaceous Durancian faults, and Tethyan faults are distinguished on the section. The geometry of the Montpellier Thrust is consistent with the deformation style observed at surface level. The Jurassic series of the hangingwall is affected by folding, including subvertical axial planes running parallel to the thrust trace, over a 5 km wide area. Such folding is associated with subordinate thrusts (not represented on this section, but in Fig. 7—section A) locally bringing the marly Liassic unit at surface, close to MUR borehole (Fig. 2). South of the Sainton Fault, the amount of deformation in the hangingwall decreases significantly, with monocline structures and long-wavelength folding involving late Cretaceous sequences. This abrupt decrease of deformation of the hangingwall is interpreted as the expression of the footwall ramp gently dipping to the south and then into the Palaeozoic basement, while the wide Late Cretaceous, syncline results from varying dip of the basement ramp.

The depth of the hangingwall flat is constrained by borehole data (MUR borehole, Fig. 2). The Sainton Fault connects with the main thrust at depth, at the transition from flat-ramp geometry (Fig. 7). While the flat-ramp geometry of the Montpellier Thrust is consistent with both surface and borehole data, there is no topographic culmination to explain such basement structure. Evidence for limited hangingwall erosion comes from bauxite remnants preserved on the present
<figure><img src="image_11.png" /><figcaption>**Fig. 10.** **A.** Geological cross-section oriented parallel to the Pyrenean shortening direction, in order to show the Pyrenean thrusting and folding in the upper crust. See location on **Figure 2. B**: Restored section in a late Cretaceous/pre-Pyrenean state, showing the localization of the EW oriented thrusts above structures inherited from the Durancian uplift event and/or Late Cretaceous differential subsidence.</figcaption></figure>

topography. This thus requires a significant footwall flexure of the autochthonous Palaeozoic basement, driven by the thrust load, in order to accommodate the thickness of the formations in the hangingwall.

The Montpellier foreland mostly comprises Palaeocene to Early Eocene syntectonic sediments while the Saint-Loup foreland consists of younger Middle Eocene syntectonic series. The thin-skinned Saint-Loup Thrust merges at depth with the Montpellier Thrust and propagates northward at the basement-cover interface. Its emergence is controlled by an inherited, north-facing normal fault of Durancian age, sealed by the bauxite surface.

To the south, the Jurassic series is buried below the Late Cretaceous series. The Thau Thrust, presently entirely covered by Neogene sediments, is identified by a deep borehole (VAL well), which documents several fault segments, folded strata and variable thickness within the Jurassic series. We interpret this structure as a detachment fold and imbricate thrusts rooted in the Triassic decollement level. Further to the south, the Adge Thrust represents the most important thrust in the region, although it is neither imaged on seismic or sampled by borehole. Its location beneath the Neogene cover is shown by a sharp gravimetric gradient (Guennoc et al., 2004) and it is evidenced by correlating the deep boreholes CAA and AGM, onshore and offshore, respectively (Fig. 2). Palaeocene syntectonic series overlies a thick (>1000 m) Late Cretaceous sequences in the onshore footwall of the Adge Thrust, while offshore, the Miocene postrift unconformably overlies the Variscan basement at 1100 m depth (Guennoc et al., 2000). Such hiatus suggests that the basement was exhumed by thrusting during the Pyrenean compression, and possibly formed a significant topography in the hinterland. The thin-skinned Thau Thrust branches on the basement ramp of the Adge Thrust, located close to the southern basin border.

Borehole data suggests thickening of the Mesozoic basin towards the centre of the basin, north of the VAL well. Both Tethyan and Durancian, north-facing, normal faults are likely to have accommodated such thickening (Fig. 10B).

The results of our restoration show that the Montpellier and Thau-Agde thrusts accommodated 13 and ≥15 km of shortening, respectively (Fig. 10). The displacement

accommodated by the Adge Thrust is only a minimal value, as the amount of created hangingwall topography cannot be estimated, since it was later eroded prior to Miocene postrift sedimentation. Therefore, the total shortening along this transect is at least 28 km.

The restoration to the pre-Pyrenean state exhibits a Jurassic high between two distinct Cretaceous basins. To the north, an Early Cretaceous (Neocomian) basin is bounded by the Saint-Loup Fault, which is sealed by the bauxite surface (Fig. 6B). To the south, a Late Cretaceous basin developed above the bauxite surface, prior to the main Pyrenean phase. The driving force behind the formation of this southern basin is questionable. It could reflect the flexural response to an early phase of the Pyrenean orogeny (e.g. Ternois et al., 2019). Alternatively, and more likely, it may represent the thermal recovery of the lithosphere, thinned during the Aptian-Albian rifting (Grool et al., 2018; Rougier et al., 2016). Better resolution of the chrono-stratigraphy of this Late Cretaceous basin is required in order to fuel this discussion. Future paleothermal analyses of Late Cretaceous cores from the boreholes would also help to understand the origin of the Late Cretaceous basin.

The restored cross-section depicts a repeated pattern of localisation of the Pyrenean thrusts above inherited horsts and normal faults corresponding to mid-Cretaceous Durancian and/or Jurassic Tethyan extensional faults. The Saint-Loup Thrust emerged above a N-facing Durancian fault, informally named the Saint-Loup Fault. Similarly, the Montpellier thrust emerges above the Durancian Montpellier Fault (Fig. 6B). The Thau Thrust is localized above the N-facing extensional Sète Fault of Tethyan age (Figs. 7B and 8B). Thus, the Montpellier-Sète structural high has been formed since at least the mid-Cretaceous Durancian uplift. This structural high controls the transition from thin-skinned to thick-skinned thrusting.

## 9 Discussion

Restoration of 2D structural sections involves several simplifications. The estimated shortening must be computed along a section, which is parallel to the shortening direction. In this study, we determine the Pyrenean shortening direction by using the compressional stress trajectories derived from
<figure><img src="image_12.png" /><figcaption>**Fig. 11.** Fault activity chronological diagram, highlighting multiple reactivations according to (i) their orientation and (ii) the successive tectonic stages recorded in the area. Faults are identified on Figure 2.</figcaption></figure>

microtectonic structures (Arthaud and Laurent, 1995). This does not account for possible decoupling of the sedimentary cover from the basement, including possible rotations such as those, documented south of the study area in the Corbières (Frizon de Lamotte et al., 1995; Rouvier et al., 2012). Decoupling of the Mesozoic cover from the Variscan basement, within the Keuper shales and evaporite (Late Triassic) has been evidenced in Corbières (Ford and Vergés, 2020) and in Provence (Espurt et al., 2019b). In the study area, possible salt tectonics may have affected the Mesozoic sequences especially in areas where the Keuper interval is thick (≥800 m in the lower Hérault valley or towards the South-East Basin, Fig. 2). However, in the central part of the study area, the Keuper sequence is much thinner (about 100 m) and there are no field or seismic evidence of intense salt

tectonics involving diapirs. It appears that the Keuper interval just behaved as a low friction interface, that allowed decoupling of the cover.

The shortening estimate derived from restoration of sections gives only a minimal value, as it does not take into account possible lateral movements along strike-slip faults parallel to the contraction direction (e.g. Krézsek et al., 2013). The post-Late Jurassic left-lateral displacement of the Cévennes Fault has been estimated up to 17 km (Bodeur, 1980). The post-Pyrenean and pre-rifting (Priabonian) left-lateral displacement between Cévennes and Nîmes faults has recently been estimated to 5 km (Séranne et al., 2021). Unfortunately, the strike-slip offset across the Nîmes fault during Pyrenean shortening is not constrained. The amount of shortening across thrusts provides a reasonable estimate of the
<figure><img src="image_13.png" /><figcaption>**Fig. 12.** A. Simplified Pyrenean structural sketch of the study area, superimposed onto the top-of-basement map (shades of grey). The overall structure corresponds to a left-lateral wrenching accommodating oblique convergence; the thrusts appear as restraining bends or splays of major, NE-striking basement strike-slip faults. B. Simplified rifting structural sketch of the study area. The synrift isopach map (shades of pink overlay) highlights the major extensional fault. The NW-trending transfer zones did not appear in the previous stage: they are thus newly formed during rifting. Note that synrift structures in the SW corner of the area are poorly documented in the present study.</figcaption></figure>

offset across their lateral or oblique ramps; accordingly, the 28 km shortening measured along the section sub-parallel to the major Cévennes and Nîmes faults (Fig. 10) suggests a Pyrenean left-lateral displacement of similar amplitude, distributed between the two faults.

Similarly, computing synrift amount of extension provides only minimal values. The exact amount of extensional displacement along inverted Pyrenean thrusts cannot be determined because the offset of the initial thrust is unknown. For example, the Sainton Thrust has been completely inverted during Oligocene extension and presently displays a normal fault setting. In contrast, the basement ramp of the Montpellier Thrust kept its thrust geometry, even though it may have been reactivated as an extensional fault during rifting (Fig. 9).

Within the northern foreland of the Pyrénées-Provence thrust belt, the Agde Thrust is located in a similar structural setting as the North Pyrenean Frontal Thrust (Grool et al., 2018) and the Cap Sicié thrust in Provence (Bestani et al., 2016). We find discrepancies in the amount of shortening in the northern foreland of the orogen. The late Cretaceous to Eocene 28 km of shortening across the Languedoc section (Fig. 10) is higher than the 13.3 km of shortening across the northern foreland of eastern Pyrénées (Grool et al., 2018), and lower than the 40 km of shortening across central Provence (Bestani et al., 2015). Such discrepancies could be accounted for by the different width of the restored orogenic wedge. More likely, the deformed Pyrénées-Provence foreland is segmented by oblique strike-slip faults, such as Cévennes, Nîmes, Durance Faults, which separates blocks accommodating varying

amounts of shortening (e.g. Lacombe and Jolivet, 2005). Unlike the Eastern Pyrénées and Provence characterized by east-west oriented thrusts, Pyrenean deformation style in Languedoc is characterized by a complex association of NNE-trending strike-slip faults associated with N-verging thrusts located between the major Cévennes and Nîmes Faults. The Pyrenean northern thrust front is therefore displaced some 60 km northward across the Corbières-Cévennes transfer zone (Mascle et al., 1994; Rouvier et al., 2012). In Languedoc, the decollement of the Mesozoic cover is not associated with significant salt tectonics, in contrast to what has been described in the neighboring areas of Corbières and Provence (Bestani et al., 2015; Espurt et al., 2019b; Ford and Vergés, 2020). This is probably due to a thinner Triassic sequence and dominance of siliciclastic sediments over the Keuper evaporites (Lopez, 1992).

The faults identified in this study reveal different histories of tectonic activity (Fig. 11). The NNE trending regional faults (Cévennes and Nîmes Faults) were initiated as strike-slip faults during the late Variscan plate reorganization (Arthaud and Matte, 1975), and they have later been reactivated during each of the successive tectonic phase, as extensional or strike-slip faults. Secondary faults of similar orientation have only been reactivated once or twice, as normal faults or left-lateral strike-slip faults. Other secondary faults (see names on Fig. 2) of either EW (i.e. Cambellies Fault) or N-S orientation (i.e. Mosson Fault) have been active only during the Durancian uplift. Finally, north-verging thrusts of E-W orientation are not reactivated structures, they were formed during Pyrenean
deformation. The tectonic chronology chart of Figure 11 shows that the Pyrenean shortening phase stands out, with the association of newly formed, E-W thrusts and reactivated NNE-trending structures as strike-slip faults. The general structural pattern of the Pyrenean deformation in onshore Languedoc is that of a sinistral wrench zone along inherited NNE-trending crustal faults, with newly formed E-W thrusts in restraining bends and right-stepping relays zones (Christie-Blick and Biddle, 1985) (Fig. 12A). The Languedoc area therefore corresponds to a transpressional transfer zone between Pyrénées (SW of the study area) and Provence (to the E) segments of the orogen (Lacombe and Jolivet, 2005). The study area is a northeastward extension of the Corbières transfer (Frizon de Lamotte et al., 1995; Mascle et al., 1994; Rouvier et al., 2012), in line with the Cévennes Fault (Arthaud and Matte, 1975), which corresponds to the SE margin of the Massif Central Variscan basement. During Pyrenean shortening, this oblique transpressional zone corresponded to a weak lithosphere, which had been thinned during Tethyan and mid-Cretaceous rifting phases. It separated the thicker and more rigid Massif Central lithosphere from a continental block, extending from the present-day Gulf of Lion to Sardinia (Advokaat et al., 2014; Bestani et al., 2016; Vinciguerra, 2020).

Oligocene extensional faults mostly reactivated Pyrenean left-lateral strike-slip faults, which detach above the Triassic decollement in the onshore domain (Séranne et al., 1995). Offshore, SE of the Nîmes Fault, extensional faults affect the basement and seem to have initiated during this event. Another specific feature of the structural framework developed during rifting of the Gulf of Lion, is the formation of transfer zones (Sétoise and Arlésienne transfer zones; Fig. 12B) parallel to the extension direction that affected the whole Gulf of Lion margin (Gorini et al., 1994; Maillard et al., 2020; Séranne et al., 1995). The Sétoise transfer zone separates two distinct synrift depocenters in the onshore Hérault Basin; it also bounds to the SW the basement high offshore Camargue, and acts as lateral ramp for the main synrift depocenters (Vistrenque and Petite Camargue grabens). Such transfer zones correspond to an area of distributed deformation in the upper crust and sedimentary cover that accommodates mid- or lower crustal strike-slip faulting. The major Catalan transfer zone, south of the study area (Jolivet et al., 2021; Maillard et al., 2020), is interpreted as the response of the overriding plate above slab tears, in a retreating subduction (Romagny et al., 2020). However, it is unlikely that the minor transfer zones documented here depict multiple slab tears, instead, they probably represents partitioning of the extensional deformation in the proximal part of the overriding plate, driven by shearing at the base of the lithosphere, and leading to lower continental crust exhumation in the distal part of the margin (Jolivet et al., 2020).

## 10 Conclusion

1. This structural analysis of the Languedoc area demonstrates complex 3D interaction of inherited upper-crustal structures, expressed by multiple fault reactivations through time, under different successive stress orientations. Pyrenean structures and polyphase tectonics are deciphered by linking both onshore and offshore domains across the Gulf of Lion proximal margin. Integration of

new acquisitions and reappraisal of published data give rise to newly described Pyrenean and pre-Pyrenean structures in the coastal Montpellier-Sète area.

2. In Languedoc, two Mesozoic extensional phases: Tethyan (Triassic-Liassic rifting and subsequent thermal subsidence) followed by Durancian (mid-Cretaceous Pyrenean rifting and northern rift-shoulder uplift) reactivated NE-SW and E-W oriented structures, respectively, as extensional faults. Both Tethyan and Durancian Mesozoic basins were bounded southward by a non-subsiding continental block extending in the present-day Gulf of Lion.

3. During oblique convergence of the Gulf of Lion-Sardinia continental block with the NE-trending margin of the Variscan Massif Central, the NE-SW faults were reactivated as left-lateral strike-slip faults, and E-W oriented thrusts were generated. The latter acted as contractional splays of the formers, and they emerged above the previous E-W Durancian faults. The Languedoc and Corbières areas therefore correspond to a NE-trending transfer zone that accommodates oblique convergence, and link the Iberian Pyrénées with Provence.

4. The Montpellier Thrust is interpreted as a basement-involved structure with a double flat-ramp geometry. Further to the south, in the coastal area, the Thau Thrust is a newly described major structure that links eastwards with the Nîmes Fault. Southwestwards, the thin-skinned Thau Thrust branches on the basement ramp of the Agde Thrust.

5. The ≥28 km amount of shortening, estimated in the Languedoc foreland, is a minimal value of the Pyrenean shortening, and must be used with caution as it does not include the complete left-lateral offset across the Corbières-Languedoc transfer zone. Our section parallel to the NNE directed oblique convergence does not show significant reactivation during later NW-SE rifting; consequently, this single step restoration enables us to reconstruct the pre-orogenic setting, which is now buried below the Gulf of Lion syn- and postrift sequences.

6. Oligocene rifting structures formed in response to a NW-SE extension, which reactivates the previous NE-trending strike-slip faults. Onshore thin-skinned extensional tectonics ramps down into the basement, through the reactivation of the Nîmes Fault and of the inherited Montpellier and Agde basement thrusts, respectively. Finally, the distinctive NW-trending transfer zones responsible for the segmentation of the Gulf of Lion margin represent newly-formed rifting structures, which accommodate differential extension rates.

*Acknowledgements.* This work is part of the Dem'Eaux Thau Project, funded by: European Community, French State, Région Occitanie, Agence de l'Eau RMC, SMBT, City of Balaruc, and Montpellier Méditerranée Métropole. This contribution benefited from discussions within the Dem'Eaux Thau research group and the Far-field Working Group of the OROGEN Project. The authors are grateful for the comments and suggestions provided by Johanna Lofi and Pierre Labaume.

We are indebted to Armin Dielforder and Olivier Averbuch for their insightful and thorough reviews, which improved the manuscript. We extend our thanks to the Guest Editors Sylvain Calassou and Olivier Lacombe who provided constructive comments and helpful suggestions.
References

Advokaat EL, Hinsbergen DJJV, Maffione M, Langereis CG, Vissers RLM, Cherchi A, et al. 2014. Eocene rotation of Sardinia, and the paleogeography of the western Mediterranean region. *Earth and Planetary Science Letters* 401: 183–195.

Alabouvette B, Berger G, Demange M, Guérangé-Lozes J. 2001. Montpellier - Carte géologique de la France à 1/250 000⁶. BRGM.

Andrieux J, Mattauer M, Tomasi P, Martinez C, Reille JL, Matte P, et al. 1971. Carte géologique de la France au 1/50 000. In: Feuille de Montpellier n° 990. BRGM.

Angrand P, Mouthereau F, Masini E, Asti R. 2020. A reconstruction of Iberia accounting for W-Tethys/N-Atlantic kinematics since the late Permian-Triassic. *Solid Earth* 11(4): 1313–1322.

Arthaud F, Matte P. 1975. Les décrochements tardi-hercyniens du SW de l'Europe. Géométrie et essai de reconstitution des conditions de la déformation. *Tectonophysics* 25: 139–171.

Arthaud F, Mégard F, Séguret M. 1977. Cadre tectonique de quelques bassins sédimentaires. *Bulletin des Centres de Recherches Exploration – Production Elf-Aquitaine* 1: 147–188.

Arthaud F, Ogier M, Séguret M. 1981. Géologie et géophysique du Golfe du Lion et de sa bordure nord. *Bulletin du BRGM* 1(2): 175–193.

Arthaud F, Séguret M. 1981. Les structures pyrénéennes du Languedoc et du Golfe du Lion (Sud de la France). *Bulletin de la Société Géologique de France* XXIII: 51–63.

Arthaud F, Durand V. 1989. Relations entre chevauchements et décrochements dans la tectoniques pyrénéenne du « pli de Montpellier » (Languedoc, France). *C R Acad Sci Paris* 303, Serie II: 1387–1393.

Arthaud F, Laurent P. 1995. Contraintes, déformation et déplacement dans l'avant-pays Nord-pyrénéen du Languedoc méditerranéen. *Geodynamica Acta* 8: 142–157.

Asensio E, Khazaradze G, Echeverria A, King RW, Vilajosana I. 2012. GPS studies of active deformation in the Pyrenees. *Geophysical Journal International* 190: 913–921.

Bache F, Olivet J-L, Gorini C, Aslanian D, Labails C, Rabineau M. 2010. Evolution of rifted continental margins: The case of the Gulf of Lions (Western Mediterranean Basin). *Earth and Planetary Science Letters* 292: 345–356.

Barbarand J, Lucazeau F, Pagel M, Séranne M. 2001. Burial and exhumation history of the south-eastern Massif Central (France) constrained by an apatite fission-track thermochronology. *Tectonophysics* 335: 275–290.

Bardossy G, Dercourt J. 1990. Les gisements de bauxites téthysiens (Méditerranée, Proche et Moyen Orient) ; cadre paléogeographique et contrôles génétiques. *Bulletin de la Société Géologique de France* 6: 869–888.

Beaudrimont AF, Dubois P. 1977. Un bassin mésogéen du domaine péri-Alpin: le sud-est de la France. *Bulletin des Centres de Recherches Exploration – Production Elf-Aquitaine* 1: 261–308.

Benedicto A. 1996. Modèles tectono-sédimentaires de bassins en extension et style structural de la marge passive du Golfe du Lion (Partie Nord), sud-est de la France. Université de Montpellier II, 235 p.

Benedicto A, Labaume P, Séguret M, Séranne M. 1996. Low-angle crustal ramp and basin geometry in the Gulf of Lion passive margin: the Oligocene-Aquitanian Vistrenque graben, SE France. *Tectonics* 15: 1192–1212.

Benedicto A, Séguret M, Labaume P. 1999. Interaction between faulting, drainage and sedimentation in extensional hanging-wall syncline basins: example of the Oligocene Matelles Basin, Gulf of Lion margin (S.E. France). In: Durand B, Jolivet L, Horváth F, Séranne M, eds. The Mediterranean Basins: Tertiary extension

within the Alpine Orogen. *The Geological Society, London, Special Publication* 156: 81–108.

Berger G, Feist R, Freytet P. 1981. Carte géologique de la France au 1/50 000. In: Feuille de Pézenas. BRGM.

Bestani L, Espurt N, Lamarche J, Floquet M, Philip J, Bellier O, et al. 2015. Structural style and evolution of the Pyrenean-Provence thrust belt, SE France. *Bulletin de la Société Géologique de France* 186: 223–241.

Bestani L, Espurt N, Lamarche J, Bellier O, Hollender F. 2016. Reconstruction of the Provence Chain evolution, southeastern France. *Tectonics* 35: 1506–1525.

Bodeur Y. 1980. Kimmeridgien supérieur et Portlandien du Languedoc. *Geobios Mémoire* (4): 77–83.

Bonjoly D, Perrin J, Roure F, Bergerat F, Courel L, Elmi S, et al. 1996a. The Ardèche paleomeargin of the South-East Basin of France: mesozoic evolution of a part of the tethyan continental margin. *Marine and Petroleum Geology* 13: 607–623.

Bonjoly D, Sureau JF, Steinberg M. 1996b. GPF Ardèche: a scientific drilling programme for fluid-rock interaction studies along an extensional palaeo-margin. *Marine and Petroleum Geology* 13: 605–606.

Bourgeois T. 2019. Caractérisation du paleokarst sous couverture du Démonstrateur DEM'EAUX Thau, Bassin de Montbazin-Gigean et de l'extrémité Est de l'étang de Thau (Hérault) : approche par sismique passive. Montpellier, 50 p.

Butler RWH, Tavarnelli E, Grasso M. 2006. Structural inheritance in mountain belts: An Alpine Apennine perspective. *Journal of Structural Geology* 28: 1893–1908.

Canva A. 2018. Modélisation structurale du Golfe du Lion à partir de données géophysiques magnétiques et gravimétriques. Master 2. ENSEGID-Bordeaux INP, 64 p.

Canva A, Thinon I, Peyrefitte A, Couëffé R, Maillard A, Jolivet L, et al. 2020. The Catalan magnetic anomaly: Its significance for the crustal structure of the Gulf of Lion passive margin and relationship to the Catalan transfer zone. *Marine and Petroleum Geology* 113: 104174.

Capar L, Marc S. 2018a. Retraitement de lignes sismiques 2D dans la région de l'étang de Thau. Rapport final. BRGM/RP-68318-FR, 23 p.

Capar L, Marc S. 2018b. Retraitement de 3 lignes sismiques 2D dans la région de Florensac. Rapport final. BRGM/RP-68603-FR, 19 p.

Chanvry E, Marchand E, Lopez M, Séranne M, Le-Sout G, Vinches M. 2020. Tectonic and climate control on allochthonous bauxite deposition. Example from the mid-Cretaceous Villeveyrac basin, southern France. *Sedimentary Geology* 407.

Chevrot S, Sylvander M, Diaz J, Martin R, Mouthereau F, Manatschal G, et al. 2018. The non-cylindrical crustal architecture of the Pyrenees. *Scientific Reports* 8: 9591.

Christie-Blick N, Biddle KT. 1985. Deformation and basin formation along strike-slip. In: Christie-Blick N, Biddle KT, eds. Strike-Slip Deformation, Basin Formation, and Sedimentation. *SEPM, Special Publication* 37: 1–34.

Christophoul F, Ford M, Grool A, Rougier G, Hemmer L. 2016. Tectono-stratigraphic evolution of the northeastern Pyrenean Foreland. In: International Geological Congress, paper 3188.

Clauzon G. 1979. The Messinian Durance Canyon (Provence, France): palaeogeographic evidence for a desiccated deep basin. *Palaeogeography Palaeoclimatology Palaeoecology* 29: 15–40.

Combes P-J. 1965. Remarques sur un type particulier de failles normales dans le massif de la Gardiole (Hérault). *C R Somm Soc Géol Fr* 2: 55–56.

Combes P-J. 1990. Typologie, cadre géodynamique et genèse des bauxites françaises. *Geodinamica Acta* 2: 91–109.
Combes PJ, Peybernès B, Fondecave-Wallez M-J, Séranne M, Lesage J-L, Camus H. 2007. Latest-Cretaceous/Paleocene karsts with marine infillings from Languedoc (South of France); paleogeographic, hydrogeologic and geodynamic implications. *Geodinamica Acta* 20: 301-326.

Couëffé R. 2011. Caractérisation géologique du secteur de l'Agglomération de Béziers à partir de la réinterprétation des profils sismiques et définition de zones favorables pour l'implantation d'ouvrages de reconnaissance. BRGM, 142 p.

Coppo N, Jacob T, Bitri A, Darnet M, Bretaudieu F, Wawrzyniak P. 2018. Investigations géophysiques (CSEM, gravimétrie et sismique) sur le secteur de Balaruc-les-Bains (Hérault, France) – Rapport d’acquisition. BRGM, RP-68318-F, 77 p.

Cravatte J, Dufaire P, Prim M, Rouaix S. 1974. Les forages du Golfe du Lion. Stratigraphie, sédimentologie. *Notes et Mémoires, Compagnie française des pétroles* 11: 209-274.

Crochet JY. 1984. Géologie et paléontologie de la partie septentrionale du fossé oligocène des Matelles (Hérault, sud de la France). *Géologie de la France* 1-2: 91-104.

Debrand-Passard S, Courbouleix S. 1984. Synthèse géologique du Sud-Est de la France – Stratigraphie et paléogéographie. Orléans, France: BRGM, 615 p.

Debroas E-J. 1990. Le flysch noir albo-cénomanien temoin de la structuration albienne à sénonienne de la Zone nord-pyreneenne en Bigorre (Hautes-Pyrenees, France). *Bulletin de la Société Géologique de France* VI: 273-285.

Diaz J, Vergés J, Chevrot S, Antonio-Vigil A, Ruiz M, Sylvander M, et al. 2018. Mapping the crustal structure beneath the eastern Pyrenees. *Tectonophysics* 744: 296-309.

Dielforder A, Frasca G, Brune S, Ford M. 2019. Formation of the Iberian-European convergent plate boundary fault and its effect on intraplate deformation in Central Europe. *Geochemistry, Geophysics, Geosystems* 20: 2395-2417.

Dörfliger N, Le Strat P. 2001. Étude du pourtour est de l'étang de Thau (Hérault) – phase II – Définition du modèle géologique et inventaire des phénomènes karstiques, volume 1. BRGM/RP-50786-FR, 136.

Dubois P, Delfaud J. 1989. Le Bassin du Sud-Est. In: Technip, ed. Dynamique et méthodes d'étude des bassins sédimentaires. Paris: Association des sédimentologistes français, pp. 277-297.

Dufréchou G, Martin R, Bonvalot S, Bruinsma S. 2019. Insight on the western Mediterranean crustal structure from GOCE satellite-gravity data. *Journal of Geodynamics* 124: 67-78.

Duretz T, Asti R, Lagabrielle Y, Brun JP, Jourdon A, Clerc C, et al. 2020. Numerical modelling of Cretaceous Pyrenean Rifting: the interaction between mantle exhumation and syn-rift salt tectonics. *Basin Research* 34: 652-667.

Duvail C, Gorini C, Lofi J, Le Strat P, Clauzon G, Dos Reis AT. 2005. Correlation between onshore and offshore Pliocene-Quaternary systems tracts below the Roussillon Basin (eastern Pyrenees, France). *Marine and Petroleum Geology* 22: 747-756.

Espurt N, Angrand P, Teixell A, Labaumed P, Ford M, Blanquat MDS, et al. 2019a. Crustal-scale balanced cross-section and restorations of the Central Pyrenean belt (Nestes-Cinca transect): Highlighting the structural control of Variscan belt and Permian-Mesozoic rift systems on mountain building. *Tectonophysics* 764: 25-45.

Espurt N, Wattellier F, Philip J, Hippolyte JC, Bellier O, Bestani L. 2019b. Mesozoic halokinesis and basement inheritance in the eastern Provence fold-thrust belt, SE France. *Tectonophysics* 766: 60-80.

Fillon C, Mouthereau F, Calassou S, Pik R, Bellahsen N, Gautheron C, et al. 2020. Post-orogenic exhumation in the western Pyrenees: evidence for extension driven by pre-orogenic inheritance. *Journal of the Geological Society* 178.

Ford M, Hemmer L, Vacherat A, Gallagher K, Christophoul F. 2016. Retro-wedge foreland basin evolution along the ECORS line, eastern Pyrenees, France. *Journal of the Geological Society* 173. https://doi.org/10.1144/jgs2015-1129.

Ford M, Vergés J. 2020. Evolution of a salt-rich transtensional rifted margin, eastern North Pyrenees, France. *Journal of the Geological Society* 178: 2019-2157.

Fournier F, Tassy A, Thinon I, Münch P, Cornée J-J., Borgomano J, et al. 2016. Pre-Pliocene tectonostratigraphic framework of the Provence continental shelf (eastern Gulf of Lion, SE France). *Bulletin de la Société Géologique de France* 187: 187-216.

Freytet P, Plaziat JC. 1982. Continental carbonate sedimentation and pedogenesis – Late Cretaceous and early Tertiary of southern France. *Contributions to Sedimentology* 12: 1-213.

Frizon De Lamotte D, Guezou J-C, Averbuch O. 1995. Distinguishing lateral folds in thrust-systems; examples from Corbières (SW France) and Betic Cordilleras (SE Spain). *Journal of Structural Geology* 17: 233-244.

Frizon de Lamotte D, Souque C, Grelaud S, Robion P. 2002. Early record of tectonic magnetic fabric during inversion of a sedimentary basin. Short review and examples from the Corbières transfer zone (France). *Bulletin de la Société Géologique de France* 172: 461-469.

Gaviglio P, Gonzales J-F. 1987. Fracturation et histoire tectonique du bassin de Gardanne (Bouches-du-Rhone). *Bulletin de la Société Géologique de France* 8: 675-682.

Gignoux M. 1926. Géologie stratigraphique. Paris: Masson & Cie, 588 p.

Gomez-Romeu J, Masini E, Tugend J, Ducoux M, Kuznir N. 2019. Role of rift structural inheritance in orogeny highlighted by the Western Pyrenees case-study. *Tectonophysics* 766: 131-150.

Gorini C, Viaillard P, Deramond J. 1991. Modèle d'inversion négative: la tectonique extensive post-nappe du fossé de Narbonne-Sigean (Corbières, sud de la France). *CR Acad Sci Paris* 312: 1013-1019.

Gorini C. 1993. Géodynamique d'une marge passive : le Golfe du Lion (Méditerranée Occidentale). Univ. Paul Sabatier-Toulouse 3, 256 p.

Gorini C, Le Marrec A, Mauffret A. 1993. Contribution to the structural and sedimentary history of the Gulf of Lions (western Mediterranean), from the ECORS profiles, industrial seismic profiles and well data. *Bulletin de la Société Géologique de France* 164: 353-363.

Gorini C, Mauffret A, Guennoc P, Le Marrec A. 1994. Structure of the Gulf of Lions (Northwest Mediterranean Sea): A review. In: Mascle A, ed. Hydrocarbon and petroleum geology of France. Special Publication of the European Association of Petroleum Geology No. 4. Berlin: Springer-Verlag, pp. 223-243.

Gottis M, Denizot G, Combes PJ, Bertrand JM, Seguret M, Gérard M, et al. 1967. Carte Géologique de la France au 1/50 000. In: Feuille de Sète n° 1016. BRGM.

Goula X, Olivera C, Fleta J, Grellet B, Lindo R, Rivera LA, et al. 1999. Present and recent stress regime in the eastern part of the Pyrenees. *Tectonophysics* 308: 487-502.

Grool AR, Ford M, Vergés J, Huismans RS, Christophoul F, Dielforder A. 2018. Insights into the crustal-scale dynamics of a doubly vergent orogen from a quantitative analysis of its forelands: A case study of the Eastern Pyrenees. *Tectonics* 37: 450-476.

Guennoc P, Debeglia N, Gorini C, Le Marrec A, Mauffret A. 1994. Anatomie d'une marge passive jeune (Golfe du Lion - Sud France): Apports des données géophysiques. *Bulletin des Centres de Recherches Exploration - Production Elf-Aquitaine* 18: 33-57.
Guennoc P, Gorini C, Mauffret A. 2000. Histoire géologique du Golfe du Lion et cartographie du rift oligo-aquitanien et de la surface messinienne. *Géologie de la France* 3: 67-97.

Guennoc P, Martelet G, Gorini C, Mauffret A, Bache F. 2004. Structure du socle anté-rift du Golfe du Lion : synthèse et analyse des données gravimétriques et sismiques. Pau, France: GDR Marges - Atelier Golfe du Lion.

Guieu G, Roussel J. 1990. Arguments for the pre-rift uplift and rift propagation in the Ligurian-Provençal basin (Northwestern Mediterranean) in the light of Pyrenean Provençal orogeny. *Tectonics* 9: 1113-1142.

Gunnell Y, Calvet M, Brichau S, Carter A, Aguilar JP, Zeyen H. 2009. Low long-term erosion rates in high-energy mountain belts: Insights from thermo- and biochronology in the Eastern Pyrenees. *Earth and Planetary Science Letters* 278: 208-218. https://doi.org/10.1016/j.epsl.2008.1012.1004.

Guyonnet-Benaize C, Lamarche J, Masse J-P., Villeneuve M, Viseur S. 2010. 3D structural modelling of small-deformations in polyphase faults pattern. Application to the Mid-Cretaceous Durance uplift, Provence (SE France). *Journal of Geodynamics* 50: 81-93. https://doi.org/10.1016/j.jog.2010.1003.1003.

Hemelsdaël R, Séranne M, Caritg S, Courrioux G, Husson E, Jacob T, et al. 2021. Dem'Eaux Thau - Construction d'un modèle géologique 3D du secteur Montpellier-Sète et de l'étang de Thau. BRGM. RP-70789-FR,94.

Hippolyte JC, Angelier J, Bergerat F, Nury D, Guieu G. 1993. Tectonic-stratigraphic record of paleostress time changes in the Oligocene basins of the Provence, southern France. *Tectonophysics* 226: 15-35.

Husson E. 2013. Interaction géodynamique/karstification et modélisation géologique 3D des massifs carbonatés : implication sur la distribution prévisionnelle de la karstification. Exemple des paléokarsts crétacés à néogènes du Languedoc montpelliéran. Doctorat. Université Montpellier 2, 314 p.

Husson E, Guillen A, Séranne M, Courrioux G, Couëffé R. 2018. 3D Geological modeling and gravity inversion of a structurally complex carbonate area: Application for karstified massif localization. *Basin Research* 30: 766-782.

Jammes S, Manatschal G, Lavier L, Masini E. 2009. Tectonosedimentary evolution related to extreme crustal thinning ahead of a propagating ocean: Example of the western Pyrenees. *Tectonics* 28. https://doi.org/10.1029/2008TC002406.

Jolivet L, Gorini C, Smit J, Leroy S. 2015. Continental breakup and the dynamics of rifting in back-arc basins: The Gulf of Lion margin. *Tectonics* 34. https://doi.org/10.1002/2014TC003570.

Jolivet L, Romagny A, Gorini C, Maillard A, Thinon I, Couëffé R, et al. 2020. Fast dismantling of a mountain belt by mantle flow: late-orogenic evolution of Pyrenees and Liguro-Provençal rifting. *Tectonophysics* 776: 228312.

Jolivet L, Menant A, Roche V, Pourhiet LL, Maillard A, Augier R, et al. 2021. Transfer zones in Mediterranean back-arc regions and tear faults. *Bulletin de la Société Géologique de France* 192. https://doi.org/10.1051/bsgf/2021006.

Jourdon A, Laeticia LP, Mouthereau F, Masini E. 2019. Role of rift maturity on the architecture and shortening distribution in mountain belts. *Earth and Planetary Science Letters* 512: 89-99.

Krézsek C, Lăpădat A, Mațenco L, Arnberger K, Barbu V, Olaru R. 2013. Strain partitioning at orogenic contacts during rotation, strike-slip and oblique convergence: Paleogene-Early Miocene evolution of the contact between the South Carpathians and Moesia. *Global and Planetary Change* 103: 63-81.

Lacombe O, Angelier J, Laurent P. 1992. Determining paleostress orientations from faults and calcite twins: a case study near the

Sainte-Victoire Range (southern France). *Tectonophysics* 201: 141-156.

Lacombe O, Mouthereau F. 2002. Basement-involved shortening and deep detachment tectonics in forelands of orogens: Insights from recent collision belts (Taiwan, Western Alps, Pyrenees). *Tectonics* 21: 1030.

Lacombe O, Jolivet L. 2005. Structural and kinematic relationships between Corsica and the Pyrenees-Provence domain at the time of the Pyrenean orogeny. *Tectonics* 24: TC1003.

Lacombe O, Bellahsen N. 2016. Thick-skinned tectonics and basement-involved fold-thrust belts: insights from selected Cenozoic orogens. *Geological Magazine* 153: 763-810.

Ladouche B, Lamotte C, Hemelsdaël R, Pétré M-A., Dewandel B, Léonardi V, et al. 2019. DEM'EAUX Thau - Synthèse et valorisation préliminaire des données sur l'hydrosystème de Thau (34) - Rapport final. BRGM. BRGM/RP-68483-FR,313.

Lagabrielle Y, Labaume P, De Saint Blanquat M. 2010. Mantle exhumation, crustal denudation, and gravity tectonics during Cretaceous rifting in the Pyrenean realm (SW Europe): Insights from the geological setting of the lherzolite bodies. *Tectonics* 29: TC4012, 4011-4026.

Lajoinie J-P., Laville P. 1979. Les formations bauxitiques de la Provence et du Languedoc. Dimensions et distribution des gisements. *Mémoire du BRGM* 100: 146.

Le Pichon X, Rangin C, Hamon Y, Loget N, Lin JY, Andreani L, et al. 2010. Geodynamics of the France Southeast Basin. *Bulletin de la Société Géologique de France* 181: 477-501

Lescoutre R, Manatschal G. 2020. Role of rift-inheritance and segmentation for orogenic evolution: example from the Pyrenean-Cantabrian system. *Bulletin de la Société Géologique de France* 191: 18.

Lofi J, Deverchère J, Gaullier V, Gillet H, Gorini C, Guennoc P, et al. 2011. Seismic Atlas of the “Messinian Salinity Crisis” Markers in the Mediterranean and Black Seas. Paris: Commission for the Geological Map of the World and Mémoires de la Société Géologique de France, 179, 79 p.

Lopez M. 1992. Dynamique du passage d'un appareil terrigène à une plate-forme carbonatée en domaine semi-aride: le Trias de Lodève, sud de la France. Univ. Montpellier 2, 403 p.

Lopez-Mir B, Munoz JA, Garcia-Senz J. 2016. 3D geometric reconstruction of Upper Cretaceous passive diapirs and salt withdrawal basins in the Cotiella Basin (southern Pyrenees). *Journal of the Geological Society of London* 173: 616-627.

Macedo J, Marshak S. 1999. Controls on the geometry of fold-thrust belt salients. *Geological Society of America Bulletin* 111: 1808-1822.

Maerten L, Séranne M. 1995. Extensional tectonics of the Oligo-Miocene Hérault Basin (S France), Gulf of Lion margin. *Bulletin de la Société Géologique de France*: 739-749.

Magné J. 1978. Etudes microstratigraphiques sur le Néogène de la Méditerranée nord-occidentale 2: Le Néogène du Languedoc méditerranéen. Toulouse: Laboratoire Géologie Méditerranéenne - Univ Paul Sabatier, 433 p.

Maillard A, Jolivet L, Lofi J, Thinon I, Couëffé R, Canva A, et al. 2020. Transfer zones and associated volcanic province in the eastern Valencia Basin: Evidence for a hot rifted margin? *Marine and Petroleum Geology* 119: 104419.

Marchand E. 2019. Rôle des interactions tectonique-sédimentation sur l'évolution et la variabilité spatiale d'un gisement de bauxite karstique : exemple du bassin de Villeveyrac (sud de la France). Doctorat. Université de Montpellier - IMT Mines Alès, 303 p.

Marchand E, Séranne M, Bruguier O, Vinches M. 2020. LA-ICP-MS dating of detrital zircon grains from the Cretaceous allochthonous
bauxites of Languedoc (south of France): Provenance and geodynamic consequences. *Basin Research* 33(1): 270-290.

Marshak S. 1988. Kinematics of Orocline and arc formation in thin skinned orogens. *Tectonics* 7: 73-86.

Martín-Gonzalez F, Fernandez-Lozano J, De Vicente G, Crespo-Martín C, Heredia N. 2021. Role of multiple inherited basement structures on orogen geometry and evolution: Insights from analogue modelling. *Journal of Structural Geology* 144: 104267.

Mascle A, Jacquart G, Deville E. 1994. The Corbières transverse zone of the Pyrenees-Provence thrust belt (South France) — Tectonic history and petroleum. In: 6th Conference, European Association of Petroleum Geoscientists & Engineers, Vienna, P556.

Masse J-P., Philip J. 1976. Paléogéographie et tectonique Crétacé moyen en Provence: Révision du concept d'Isthme Durancien. *Revue de Géographie Physique et Géologie Dynamique* XVIII: 49-66.

Mauffret A, Gorini C. 1996. Structural style and geodynamic evolution of Camargue and Western Provençal basin, southeastern France. *Tectonics* 15: 356-375.

Mauffret A, Durand De Grossouvre B, Dos Reis AT, Gorini C, Nercessian A. 2001. Structural geometry in the eastern Pyrenees and western Gulf of Lion (Western Mediterranean). *Journal of Structural Geology* 23: 1701-1726.

Menzer L, Séranne M, Labaume P, Couëffé R. 2019. Sedimentary record of pre-Pyrenean salt tectonic inheritance in the Narbonne-Sigean synrift basin, Gulf of Lion margin. In: OROGEN Meeting CNRS-Paris, Total-BRGM-INSU.

Merzeraud G, Colombie C. 1999. Évolution morphologique des profils de dépôts dans le Sinémurien de la marge cévenole (région de Lodève). *Comptes rendus de l'Académie des Sciences de Paris* 329: 779-786.

Nguyen HN, Vernant P, Mazzotti S, Khazaradze G, Asensio E. 2016. 3-D GPS velocity field and its implications on the present-day post-orogenic deformation of the Western Alps and Pyrenees. *Solid Earth* 7: 1349-1363.

Odlum ML, Stocklia DF, Capaldi TN, Thomson KD, Clark J, Puigdefàbregas C, et al. 2019. Tectonic and sediment provenance evolution of the South Eastern Pyrenean foreland basins during rift margin inversion and orogenic uplift. *Tectonophysics* 765: 226-248.

Olivet JL. 1996. Cinématique de la plaque Ibérique. *Bulletin des Centres de Recherche Exploration-Production Elf Aquitaine* 20: 131-195.

Pétré M-A. 2020. Hydraulic and geochemical impact of occasional saltwater intrusions through a submarine spring in a karst and thermal aquifer (Balaruc peninsula near Montpellier, France). *Hydrology and Earth System Sciences* 24: 5656-5672.

Philip H, Mattauer M, Bodeur Y, Séguret M, Puech JP, Mattei J. 1978. Carte géologique de la France au 1/50 000. In: Feuille de St Martin de Londres n° 963. BRGM.

Rigo A, Vernant P, Feigl KL, Goula X, Khazaradze G, Talaya J, et al. 2015. Present-day deformation of the Pyrenees revealed by GPS surveying and earthquake focal mechanisms until 2011. *Geophysical Journal International* 201: 947-964.

Romagny A, Jolivet L, Menant A, Bessière E, Maillard A, Canva A, et al. 2020. Detailed tectonic reconstructions of the Western Mediterranean region for the last 35 Ma, insights on driving mechanisms. *Bulletin de la Societe Geologique de France* 191: 37.

Rougier G, Ford M, Christophoul F, Bader A-G. 2016. Stratigraphic and tectonic studies in the central Aquitaine Basin, northern Pyrenees: Constraints on the subsidence and deformation history of a retro-foreland basin. *Comptes Rendus Géosciences* 348: 224-235.

Roure F, Brun JP, Colleta B, Van Den Driessche J. 1992. Geometry and kinematics of extensional structures in the Alpine foreland basin of southeastern France. *Journal of Structural Geology* 14: 503-519.

Roure F, Brun JP, Coletta B, Vially R. 1994. Multiple extensional structures, fault reactivation, and petroleum plays in the Alpine foreland basin of southeastern France. In: Mascle A, ed. Hydrocarbon and petroleum geology of France, Special Publication of the European Association of Petroleum Geology No. 4. Springer-Verlag, pp. 245-268.

Rouvier H, Henry B, Le Goff M. 2012. Regional rotations evidenced by the paleomagnetism in the Corbieres orocline (France). *Bulletin de la paleomagnetisme de la Corbieres* 183: 409-424.

Sanchis E. 2000. Méthodologie d'imagerie 3D des bassins, des exemples dans le bassin du Sud Est. BRGM/Univ-Montpellier 2. BRGM GEOFR3D 98/04, 70.

Sanchis E, Séranne M. 2000. Structural style and tectonic evolution of a polyphase extensional basin of the Gulf of Lion passive margin: the Tertiary Alès Basin, southern France. *Tectonophysics* 322: 243-264.

Saspiturry N, Allanic C, Razin P, Issautier B, Baudin T, Lasseur E, et al. 2020. Closure of a hyperextended system in an orogenic lithospheric pop-up, Western Pyrenees: The role of mantle buttressing and rift structural inheritance. *Terra Nova* 32: 253-260.

Séguret M, Benedicto A, Séranne M. 1997. Structure profonde de la région du Gard Rhodanien — Apport du retraitement et de la réinterprétation de données sismiques régionales. In: Journées Scientifiques Gard, Bagnols/Cèze (Gard), pp. 31-32.

Séranne M, Benedicto A, Truffert C, Pascal G, Labaume P. 1995. Structural style and evolution of the Gulf of Lion Oligo-Miocene rifting: Role of the Pyrenean orogeny. *Marine and Petroleum Geology* 12: 809-820.

Séranne M. 1999. The Gulf of Lion continental margin (NW Mediterranean) revisited by IBS: an overview. In: Durand B, Jolivet L, Horváth F, Séranne M, eds. The Mediterranean Basins: Tertiary extension within the Alpine Orogen. *The Geological Society, London, Special Publication* 156: 15-36.

Séranne M, Couëffé R, Husson E, Baral C, Villard J. 2021. The transition from Pyrenean shortening to Gulf of Lion rifting in Languedoc (South France) — A tectonic-sedimentation analysis. *BSGF – Earth Sciences Bulletin* 192. https://doi.org/10.1051/bsgf/2021017.

Serrano O, Hanot F. 2005. Le bassin Oligo-Miocène de l'Hérault : un exemple de rétrocharriage des structures pyrénéennes. Implications hydrogéologiques. BRGM. BRGM/RP-53733-FR.

Sibson RH. 1985. A note on fault reactivation. *Journal of Structural Geology* 7: 751-754.

Soliva R, Maerten F, Maerten L, Mattila J. 2019. Fault slip envelope: a new parametric investigation tool for fault slip based on geomechanics and 3-D fault geometry. *Solid Earth* 10: 1141-1154.

Tavani S, Bertok C, Granado P, Piana F, Salas R, Vigna B, et al. 2018. The Iberia-Eurasia plate boundary east of the Pyrenees. *Earth-Science Reviews* 187: 314-337.

Teixell A, Labaume P, Lagabrielle Y. 2016. The crustal evolution of the west-central Pyrenees revisited: Inferences from a new kinematic scenario. *Comptes Rendus Géosciences* 348: 257-267.

Teixell A, Labaume P, Ayarza P, Espurt N, De Saint Blanquat M, Lagabrielle Y. 2018. Crustal structure and evolution of the Pyrenean-Cantabrian belt: A review and new interpretations from recent concepts and data. *Tectonophysics* 724-725: 146-170.

Ternois S, Odlum M, Ford M, Pik R, Stockli D, Tibari B, et al. 2019. Thermochronological Evidence of Early Orogenesis, Eastern Pyrenees, France. *Tectonics* 38: 1308-1336.

Thomas WA. 1990. Controls on locations of transverse zones in thrust belts. *Eclogae Geologicae Helvetiae* 83: 727-744.

Vinciguerra C. 2020. Reconstruction des paléo-drainages des bassins précoces péri-orogéniques (Crétacé terminal-Paléocène) à partir des dépôts fluviatiles dans le système pyrénéen oriental. Doctorat. Bordeaux 3, 307 p.
Vissers RLM, Meijer PT. 2012. Mesozoic rotation of Iberia: Subduction in the Pyrenees? *Earth-Science Reviews* 110: 93–110.

Wehr H, Chevrot S, Courrioux G, Guillen A. 2018. A three-dimensional model of the Pyrenees and their foreland basins from geological and gravimetric data. *Tectonophysics* 734–735: 16–32.

Withjack MO, Peterson ET. 1993. Prediction of normal fault geometries – A sensitivity analysis. *American Association of Petroleum Geologists Bulletin* 77: 1860–1873.

Wyns R, Quesnel F, Simon-Coinçon R, Guillocheau F, Lacquement F. 2003. Major weathering in France related to lithospheric deformation. *Géologie de la France* 1: 79–87.

**Cite this article as:** Hemelsdaël R, Séranne M, Husson E, Ballas G. 2021. Structural style of the Languedoc Pyrenean thrust belt in relation with the inherited Mesozoic structures and with the rifting of the Gulf of Lion margin, southern France, *BSGF - Earth Sciences Bulletin* 192: 46.