Although tectonic plates are usually considered as rigid blocks, intraplate deformation such as lithospheric buckling or diffuse brittle deformation has been recognized for a long time. However, the origin of these deformations remains puzzling. Indeed, whereas the chronology of deformation at plate boundaries can be constrained by numerous methods (syntectonic sedimentary record, thermochronology, etc.), dating of brittle structures (faults, veins, and joints) in the far-field domains remains challenging, preventing a global interpretation of the system as a whole. In this contribution, we have combined a tectonic study with a synkinematical geochronological study of fault-related calcites of the Grands Causses intraplate domain, north of the Pyrenean orogeny. We show that these faults record a much longer history of deformation than previously thought. The Mesozoic extension, usually attributed to an early Jurassic Tethysian rifting event, probably lasted until the Barremian-Aptian epoch, in response to the Pyrenean basin’s opening. The so-called “Pyrenean deformation” of the Grands Causses domain, usually associated with the paroxysm of deformation in the belt during the late Eocene, began much earlier, around 100 Ma, and lasted for more than 60-70 Ma. This study demonstrates the high sensitivity of an intraplate domain (Grands Causses area) to record extensional or compressional deformations occurring at the edge of neighbouring plates.

1. Introduction

Tectonic inversion at plate boundaries has been known and studied since the beginning of the 20th century, for example, through the studies of Lamplugh [1] and Stille [2], although the term “inversion” only appeared in the early 80s [3]. Since then, many other works have highlighted the widespread character of this phenomenon at plate boundaries, such as Williams et al. [4] or more recently Graveleau et al. [5].

However, the consequences of tectonic inversion at plate boundaries in the intraplate domain (far-field deformation) have not been studied as much. This deformation, several hundred kilometres away from the orogeny, is generally expressed by (1) kilometre scale denudation generally associated with lithospheric folding, (2) reactivations of crustal scale inherited faults, and (3) small-scale fracturing including low-displacement faults and tectonic joints sets—all of them being able to trigger, for instance, significant fluid flow and potential mineralization. These far-field phenomena have been characterized with numerical or analogue modelling or field studies (see examples in Cloetingh et al. [6], Gerbault et al. [7], or Sokoutis et al. [8] for lithospheric folding, Ziegler et al. [9] for intraplate fault reactivations, and Navabpour et al. [10] for small-scale fracturing), but the question of their age and the time relationships with the building of the orogen itself has remained puzzling. Lithospheric buckling relative dating can be asserted using analogue modelling. Fernández-Lozano et al. [11] showed that such long wavelength folding occurs in the early stages of basin inversion. Nonetheless, dating of crustal scale brittle deformations in the intraplate domain is still challenging, due to the lack of (1) syntectonic sedimentary records, in relation to very low displacements, and (2) adequate materials to be dated in the fault zones or joints. In the absence of such information, their age has been deduced by comparison with the generally well-constrained tectonic timeline of the adjacent domains based on stress orientation from the structure directions, on syntectonic sedimentary record, on thermochronology, etc. [9, 10]—a comparison that could in some cases lead to erroneous age attributions (see Parrish et al. [12]).

Our study questions the nature, origin, and age of intraplate deformations using the example of the Pyrenean belt and its repercussion in the far-field domain, north of the northern foreland basin (France). The studied area is located at the crossroads between several geodynamic areas with contrasting evolutions (Liguro-Provençal basin, Pyrenees, Massif Central, and Alps; Figures 1 and 2). We have focused on the Grands Causses area, 150 km north-east of the north Pyrenean frontal trust (NPFT; Figure 1). The contractional deformation of this area is usually associated with the Eocene Pyrenean tectonic pulse, although a syntectonic sedimentary record associated with these structures is lacking [1316]. In the first part of this paper, the structural context of the Grands Causses domain has been highlighted. Then, we have detailed the petrographic and geochronological (U-Pb) results of synkinematical fault-related calcite analyses performed using High Resolution Inductively Coupled Plasma Mass Spectrometry coupled with a Laser Ablation System (LA-HR-ICPMS). Finally, we have proposed a global scenario linking intraplate fault activity/fluid circulations and plate geodynamics.

2. Geological and Tectonic Setting

2.1. Overall Geodynamic Context
2.1.1. Tethys Opening

From a geodynamic point of view, the dislocation of the Pangea that began in the Triassic period was at the origin of the opening from east to west of the Tethys Ocean in the south of France during the Jurassic period [17, 18]. The Liassic extension recorded in the Grands Causses domain is likely to be related to this general event at the scale of Western Europe. The dislocation of the Pangea induced the separation of two distinct continents: Gondwana in the south and Laurussia in the north, and an extensional synsedimentary phase on the southern margin of Western Europe, as evidenced by the presence of tilted blocks [19]. The subduction of this oceanic crust would be, during the Cenozoic period, at the origin of the collision of the Eurasian and Apulian plates, initiating the formation of the Alps chain.

2.1.2. Pyrenean History

The Pyrenean history follows a classical two-step scenario including the opening and later closure and inversion of basins [2024].

First, the divergence of the Iberian/Eurasian plates induced an extensional phase and a crustal thinning from the late Jurassic to the early Cretaceous. This phase would be linked to a movement of Iberia towards the south-west of several hundred kilometres and is at the origin of a magma-poor hyperextended rift, the Pyrenean rift system [21, 2533] (Figure 3). Then, according to paleomagnetic data, a counterclockwise rotation of the Iberian plate generated convergence between Iberian and Eurasian blocks, starting from the late Cretaceous [34]. This convergence induced basin inversion and continental collision in the Pyrenean domain.

The formation of the Pyrenean orogen is at the origin of many deformations of varying intensity depending on the distance from the chain [35]. In the inner Pyrenean belt, two deformation pulses have been identified [36]: the early stages of deformation occurred during the end of the late Cretaceous (Campanian-Maastrichtian; Figure 3) [37, 38]. The Paleocene was then marked by a slowdown in convergence related to the beginning of the continental collision [38, 39]. The orogenic paroxysm has been dated between the Ypresian and the Bartonian linked to an acceleration in the convergence rate [4042] (Figure 3). Thermochronology studies evidenced two cooling events confirming these two stages of contractional deformation (Campanian-Maastrichtian and Bartonian) and a denudation phase during the Oligocene period [43, 44].

In the foreland basins, the tectonic structures associated with the formation of the Pyrenean orogen are generally oriented east-west and derive from a stress field whose σ1 is oriented north-south [45, 46]. The age of the first deformations is younger than those recorded in the inner part of the belt according to some authors who have dated them as being early late Cretaceous (Cenomanian-Turonian) [47, 48] (Figure 3). In the Provence domain (Figure 1), the first deformations are dated from the Santonian to the Danian [49], as is the Aquitaine Basin, which underwent flexural deformation in the Maastrichtian [50]. Plaziat [40] highlighted a compressive phase in the Aquitaine Basin from the Ypresian to the Bartonian. Gély and Sztràkos [51] evidenced a 1st phase of deformation associated with the flexuration of the Western European lithospheric plate during the Paleocene and then three successive tightening episodes during the Eocene (Figure 3). The recent study from Ortiz et al. [52] showed a paroxysm of shortening in the Aquitaine Basin during the Priabonian (35.8 Ma). The end of the Pyrenees shortening occurred during the Chattian (27.1 to 25.2 Ma) [52].

In the intraplate field, Briais et al. [53] highlighted two deformation phases in the Paris Basin associated with the convergence of the African and European plates during the Paleocene and the uppermost Ypresian. During the Paleocene, an exhumation phase has also been evidenced from a thermochronology study in the Morvan basin [54]. Parrish et al. [12] demonstrated intraplate deformation in Southern England during the Eocene based on U-Pb dating of calcite veins (Figure 3).

2.1.3. Oligocene Rifting of the Liguro-Provençal Basin

Geodynamic history was thereafter marked by the opening of the Liguro-Provençal oceanic basin during the Oligocene period [55, 56]. This extensional phase, associated with the northward movement of the African plate and the retreat of the subduction zone in the Mediterranean Sea, was at the origin of continental crust tearing in the back arc setting between Europe and the Corsica-Sardinia block [56, 57]. The structural style associated with this geodynamic episode was strongly influenced by the Pyrenean orogen that caused lithospheric thickening in the area [58]. From a structural point of view, the Gulf of Lion margin is currently controlled by a subparallel NE-SW normal fault network and many small synrift basins such as the Alès basin [59] showing an extension towards the SE [60].

2.2. Geology of the Studied Area

The geological history of the Grands Causses area is strongly influenced by the larger scale Western Europe history presented above. During the Mesozoic, this area was indeed under the influence of the opening of the Tethys Ocean associated with the dislocation of Pangea. It then became part of the distant foreland basin (intraplate domain) of the Pyrenean orogen since its establishment. Finally, the Grands Causses area, due to this proximity to the Liguro-Provençal basin, became an area likely to be marked by the opening of the Gulf of Lion during the Cenozoic.

Paleozoic basement rocks outcrop in the Massif Central and delimit the Grands Causses area to the north, east, and west, while the south of the basin is separated from the Languedoc by the Cévennes fault network (Figures 1 and 2). Mesozoic sedimentation began in the Triassic period with a transgressive episode inducing a continental/marine transition and the formation of evaporites and marls [61]. The Jurassic limestones, clays, and dolostones demonstrate a marine depositional environment. Nowadays, these marine formations are outcropping at elevations between 800 m and 1000 m amsl and form vast plateaus deeply incised by the river network. Above the Jurassic limestones and dolostones, very scattered Cretaceous remnants testify to a later but poorly constrained sedimentation history [62].

From a tectonic point of view, several episodes have been evidenced during the history of the Grands Causses area from the Jurassic to the present [1315], including extensional and compressional events. A first extensional event was at the origin of the formation of normal faults throughout the region during early Jurassic times. The Liassic age of this episode is well constrained due to the presence of kilometre-scale normal faults associated with progressive unconformity outcropping on the southern border of the Grands Causses domain [6365]. The compressive tectonic structures of the Grands Causses domain (either neoformed or normal fault reactivation) have been related by previous authors to the nearby orogen, i.e., the Pyrenean orogen building during the Eocene period [13, 45, 66]. Finally, a late third geodynamic episode, again extensional, has been highlighted by some authors such as Macquar [15]. However, the Liassic progressive unconformities are the only reliable argument to attribute an age to deformation. So far, the ages of later events—either extensional or compressional—have remained speculative due to the lack of post-Jurassic, potentially syntectonic, sedimentary records.

On a larger scale, vertical movements of the lithosphere deduced from the paleokarst studies of Camus [67], Bruxelles [68], and Husson [69] indicated (1) a regional uplift of the Grands Causses domain during the terminal early Cretaceous (>1600 m) and (2) an uplift during the Miocene whose amplitude varied from a few metres to 400 m. Thermochronological data from the Massif Central and the Grands Causses basin showed a significant exhumation episode (2000-2500 m of denudation) at the end of the early Cretaceous (110 Ma) suggesting erosion of a late Jurassic-early Cretaceous age sedimentary cover above the basement [7072].

3. Methodology

3.1. Tectonics

Tectonic structures (fault planes, stylolithic peak axes, and tension joints) were measured over the entire Grands Causses area in the Jurassic cover over a surface area of 2500 km2, from Ste-Enimie to the north to Navacelles to the south, and from Sévérac-le-Château to the west to Florac to the east (Figure 2). They have been inversed using WinTensor software [73] in order to determine the paleostress field orientation using the right dihedral inversion method [74]. Field observations also allowed us to determine relative chronology criteria in order to define the tectonic timeline in the Grands Causses area.

3.2. Analyses on Fault-Related Calcites

Our sampling strategy consisted of recovering synkinematic fault-related calcites that had undergone a relatively simple geological history, with a single deformation episode, in order to avoid multiple fluid circulation phases that could have induced a reopening of the isotopic system. For this reason, only calcites with a single family of slickensides were sampled for geochronological analyses.

3.2.1. Petrographic Observations

Calcite samples were analyzed in polished thin sections using optical microscopy in reflected and transmitted light and cathodoluminescence (CL). CL observations were carried out on an Olympus BX41 microscope (Olympus corporation, Tokyo, Japan) coupled to a cathodyne cold-cathode cathodoluminescence from NewTec (Nîmes, France) operating at 10 kV and 200-300 μA and a QICAM Fast 1394 digital camera (Teledyne QImaging, Surrey, Canada). The petrographic aim in this study was to highlight criteria for syndeformation fluid flow and to differentiate, as far as possible, the growth phases of calcitic minerals.

3.2.2. Calcite U-Pb Dating by LA-HR-ICPMS

U-Pb dating of some fault-related calcites by LA-HR-ICPMS was conducted using a High Resolution ICP-MS (Element XR from Thermo Scientific) coupled with a 193 nm ArF laser ablation system (Analyte Excite Excimer from Teledyne) at the GEOPS laboratory (Orsay, France). Since calcites commonly have low uranium concentrations and high common lead content, an elemental analysis of U and Pb was performed on a set of twenty-two samples in order to target samples with the highest U contents for U-Pb dating.

The calcite dating was carried out on polished samples. The laser beam diameter was set to 150 μm, using a frequency of 8 Hz and a fluence of 4 J/cm2. The analysis first consisted of the measurement of the background during 30 seconds before firing the laser for 30 seconds. Prior to U-Pb dating, a preablation was conducted at 155 μm, 8 Hz, and with a fluence of 4 J/cm2 in order to reduce potential contamination. For each sample, 40 to 90 ablation spots were analyzed. The NIST614 glass standard was used for the 207Pb/206Pb calibration, and the WC-1 calcite was used for 238U/206Pb calibration (; Roberts et al. [75]). Two secondary carbonate standards were also analyzed to validate the quality of the dating procedure (Duff Brown, —dated using the isotope dilution method, Hill et al. [76]; B6, 42.99 ± 0.98 Ma—dated from LA-ICPMS, Pagel et al. [77]). The uncertainties in brackets correspond to the propagated uncertainties (see Supplementary Materials for analytical details (available here)). Signals were processed with Iolite© software [78], and ages were obtained from the conventional Tera-Wasserburg diagram [79] using Isoplot© software [80]. Details on the method can be found in the Supplementary Materials.

4. Results

4.1. Tectonics
4.1.1. Stress Tensors

Almost 800 tectonic structures have been measured in the Grands Causses area (Figures 46, Table 1). Generally, extensional structures in this area corresponded to normal faults (Figure 4(a)) identified from tectoglyphs (slickensides on fault-related calcite) on fault mirrors or stratigraphic shifts. Compressive deformation had a much more varied range of tectonic structures: field observations showed the presence of reverse and strike-slip metric faults but also of folds with an east-west axis (Figure 4(b)). The tectoglyphs were also diverse: fault-related calcite, grooves, mechanical slickensides, stylolithic peaks, and secondary fault planes (Figures 4(c)4(f)).

Stress tensors resulting from the inversion of extensional structures revealed various orientations of the stress field (Figure 5, blue lines). The variation of the stress field existed at the regional level (for instance, on site “33”—in the west of the studied area, north of Millau—the σ3 orientation is WNW-ESE compared to site “10”—north of the Causse Méjean, in the north of the studied area—the σ3 orientation is NE-SW (Figure 5)). These variations also existed at the local level (there were, for example, two different orientations of the stress field on sites “9” (σ3 ~ NW-SE and NNE-SSW)—south-east of the Causse du Larzac, in the south of the studied area—and “31” (σ3 ~ NE-SW and NW-SE)—in the east of the Causse du Larzac, north of St-Jean-du-Bruel).

Conversely, reverse and strike-slip fault inversion gave relatively uniform results (Figure 5, red lines). Indeed, the inversion of nearly 400 microtectonic data (corresponding mainly to fault plane measurements) showed a generally homogeneous N-S shortening direction at the outcrop level as well as at the regional level.

4.1.2. Relative Chronology

The overlap of tectonic structures was observed at the outcrop scale (relative chronology between structures) as well as at the fault mirror scale (superposition of slickensides with different orientations). We evidenced an early extensional episode whose stress field orientation had not been determined precisely because of the small number of fault mirrors on which this relative chronology was observed (Figure 6). The Grands Causses area was subsequently affected by a compressive regime that was responsible for the formation of a large network of reverse and strike-slip (sinistral and dextral) faults with σ1 generally oriented north-south (Figures 5 and 6). In the area, some of the reverse faults overlap strike-slip faults (Figure 6). However, due to a lack of microstructural data, it was difficult to accurately determine the relative chronology of these tectonic structures. Based on previous studies [13, 15, 46] and the homogeneity of the stress field throughout the Grands Causses region, we believe here that both strike-slip and reverse structures are associated with a single deformation episode. Finally, a last extensional episode was observed since some strike-slip faults were crosscut by normal faults (Figure 6). This last relative chronology highlighted an extensional episode associated with a σ3 orientation roughly NW-SE (based on only three measurements).

4.2. Petrographic Observations

Optical microscopic observation of fault-related calcites showed the presence of various markers of syndeformation crystallization, such as the inclusion trails (parallel to the surface of the host rock) and inclusion bands (oblique to the surface of the host rock) highlighted by Ramsay and Huber [81], Koehn and Passchier [82], Fagereng et al. [83], and Bons et al. [84] (Figure 7). The texture of calcites varied from one sample to another as well as within a sample. Indeed, fault-related calcites may be composed of rounded grains and elongated grains (with the axis of the elongation inclined to the host rock) as well as fibers (parallel to the host rock) (Figure 7).

Cathodoluminescence microscopy showed (1) a varying degree of luminescence from one fault-related calcite to another, unrelated to the nature of the fault: cathodoluminescence observations showed a dull-orange color for CA18J06-1 sample and brown for CA18J05-1 sample (Figures 7(c) and 7(d)); (2) a luminescence of fault-related calcites which was relatively uniform across the samples except at the location of some rounded grains which showed growth zonations alternating from bright to nonluminescent (Figure 7(b)).

4.3. U-Pb Dating

U-Pb data with 2σ errors are compiled in Figure 8. The uncertainties in brackets correspond to the propagated uncertainties. Among the selected samples, only 1/3 could be dated because of the high common lead content and the low amount of uranium in the samples (Table 2). Regarding the dating of secondary standards in this study, Duff Brown was dated at () and B6 at () (Figures 8(a) and 8(b)). Normal fault-related calcite CA18J03-1 was dated at [15.4] Ma with a MSWD of 1.5 (Figure 8(c)). For strike-slip fault-related calcites, CA18J09-4 was dated at [13.3] Ma (), CA18J03-2 was dated at [6.8] Ma (), and NAV01 at [7.8] Ma (Figures 8(d), 8(h), and 8(i)). For reverse fault-related calcites, an age of [7.9] Ma () was obtained for CA18J05-1, of [15.2] Ma () for CA18J06-1, and of [4.9] Ma () for CA18J03-4 (Figures 8(e)8(g)).

5. Discussion

5.1. Characteristics of the Fault-Related Calcites

The intraplate deformation of the Grands Causses area has been characterized from the observation and analysis of tectonic structures that have undergone, one by one, a relatively simple geological history. In fact, the choice to consider small, only metric-sized, structures and to analyze single-phase fault-related calcites was aimed at reducing the risk of reopening geochemical systems, which may have led to results that would not have been representative of the considered tectonic phase.

Some previous studies have focused on the morphology, petrography, and growth of fault-related calcites (slickenfibre calcite), highlighting some criteria of syndeformation such as inclusion bands and inclusion trails and the growth direction of elongated grains associated with the σ1 main stress [8187]. Bons et al. [84] described inclusion bands as “arrays of wall-rock inclusions, dust, secondary minerals or fluid inclusions that were included or precipitated at the vein growth front.” Inclusion trails corresponded to arrays of inclusions which were “assumed to follow the opening trajectory of a vein” [82]. Although these two types of morphology associated with fault-related calcites may allow us to determine the history of the opening of the slickenfibre calcite, their interpretation is often not straightforward [82].

In this study, the petrography of the fault-related calcites confirms the syntectonic character of their crystallization due to the presence of these criteria (Figure 7). Moreover, cathodoluminescence observations showed that the studied objects crystallized during a single growth phase, except for some zones associated with rounded grains that may have growth zonations.

5.2. Unraveling the Complexity of the Extensional Footprint

Fieldwork showed that the formation of some of the normal faults occurred prior to the compressive episode which was associated with the Pyrenean phase and with the northward movement of Africa according to some authors such as Macquar [15], De Charpal et al. [14], and Constantin et al. [13] (Figure 6). According to this latter, the extension record of the whole area was mainly Jurassic in age as a consequence of the Tethys Ocean opening. This geodynamic event was highlighted in the literature by progressive unconformity outcropping [63, 64].

In this study, the U-Pb results show that some normal faults of the Grands Causses area resulted from a much later extension, most likely at the early Cretaceous in age (Barremien; Figure 8(c)), relatively simultaneous to the opening of the Pyrenean basins to the south [24]. In addition, the normal fault CA18J03-1 ( [15.4] Ma) is related to a NNE-SSW extension. By extrapolating this U-Pb age to the other normal faults of this family, the Cretaceous extension would therefore correspond to a stress field such as σ3 which is approximately oriented N-S (Figure 5). This result is consistent with the geometry of the Pyrenean basin opening [88, 89] (Figure 3). Consequently, the Mesozoic extension in the Grands Causses domain probably lasted much longer than previously proposed in response to far-field stress transmission either related to the Tethys or Pyrenean basin opening episodes.

In summary, we provide new evidence that the Grands Causses domain underwent two different extensional events during Mesozoic times: a first one during the Jurassic in response to the Tethys basin opening and a second one during the early Cretaceous as a consequence of Pyrenean basins opening along a N-S direction. The transition between both events may have been gradual, implying that the Grands Causses domain has been subjected to a long-lived extension episode since the Lias to the early Cretaceous in response to plate boundaries’ deformation. Moreover, based on the overlap of normal slickensides above strike-slip slickensides (Figure 6), a late extensional geodynamic episode was recorded in the Grands Causses area after the Pyrenean compressive episode (see discussion in the next section). Although no calcite related to these faults was dated, this episode could be interpreted as a consequence of the opening of the Liguro-Provençal basin during Oligocene times [46]. The NW-SE orientation of the stress field associated with these structures is indeed consistent with this geodynamic episode [56].

Concerning the stress field associated with the extension in the Grands Causses area, the microtectonic data showed a significant heterogeneity of its orientation (Figure 5). The early extension during the Cretaceous could be associated with a stress field with σ3 north-south while the late extension to the Cenozoic could come from a stress field σ3 NW-SE. However, this heterogeneity in the set of stress tensors can be explained by (1) a superposition of several deformation episodes since the Mesozoic, (2) the compaction of formations at the origin of extensional structures without preferential orientation, and (3) fault planes formed from preexisting planes (fractures, reverse, and strike-slip fault planes for late extension).

5.3. So-Called “Pyrenean Deformation” of the Grands Causses Domain

The Grands Causses area is characterized by many reverse and strike-slip fault planes resulting from a rather homogeneous stress field with a N-S shortening orientation (Figure 5). These compressive structures are usually considered as contemporaneous with the Pyrenean belt building during Eocene times [13, 14], but it remains a purely speculative affirmation in the absence of a syntectonic sedimentary record.

The U-Pb ages have shed a new light on this commonly accepted idea. Indeed, although the age uncertainties are relatively high (, Figure 8) due to low uranium content and the relative abundance of common lead in calcites, the U-Pb ages showed the record of a continuous brittle deformation over several tens of millions of years, from Albian to Priabonian times (older age: [13.3] Ma; younger age: [7.8] Ma). This indicates that the intraplate brittle deformation that occurred in response to basin closure at the Iberia/Europe plate boundary lasted much longer than previously thought and is not simply the record of the paroxysmal phases of mountain building of the terminal Eocene age [45].

Furthermore, reconstruction of Iberian and European plate kinematics indicates that the beginning of compression occurred around 85 Ma [24, 41, 90]. This value is within the uncertainties of three of the ages obtained on compression-related structures (sample CA18J05-1: [7.9] Ma, sample CA18J09-4: [13.3] Ma, and sample CA18J06-1: [15.2] Ma). As a consequence, it appears that part of the intraplate brittle deformation in the Grands Causses area also recorded the very early phases of tectonic plate convergence, probably associated with the flexuration and denudation of the lithosphere [44].

6. Conclusion

Until now, the lack of sedimentary records associated with deformation in the Grands Causses area has led to a poorly constrained geodynamic calendar that appeared to be controversial in various studies [1315, 46]. In this study, the combination of a structural study over the whole area with a petrographic and geochronological study on relatively simple objects such as fault-related calcites allows us to establish a more detailed calendar of the brittle intraplate deformation.

This study demonstrates the high sensitivity of an intraplate domain (Grands Causses area) to record deformations occurring at the edge of tectonic plates. Indeed, not only does the intraplate domain record the main relatively brief events occurring at the plate boundary such as the Eocene Pyrenean pulse, but the far-field deformation is also long and continuous over several tens of millions of years since the Mesozoic. The tectonic structures of the Grands Causses domain have formed in response to both the extensional deformations that occurred in Southern France (opening of the Tethys Ocean in the Jurassic, of the Pyrenean basins in the early Cretaceous, and, possibly, of the Liguro-Provençal basin in the Oligocene) and also the Pyrenean compression, from the late Cretaceous to the Eocene.

Strikingly, although this domain seems to record all the geodynamic episodes that occurred in Southern France during the Meso-Cenozoic period, there is no evidence in the studied area of the closure of the Tethys Ocean and subsequent formation of the alpine belt 300 km to the west. Many aspects of stress transmission from plate boundaries to the intraplate domains thus still remain to be investigated, which will require a well constrained calendar of the intraplate domain deformation—a challenging task to establish.

Data Availability

The analytical and field data used to support the findings of this study are included within the article and within the supplementary information file.

Conflicts of Interest

The authors declare that they have no conflicts of interest.


The authors are grateful to Thomas Blaise and Benjamin Brigaud (GEOPS laboratory, Paris-Saclay University) for fruitful exchanges on U-Pb calcite dating. U-Pb analyses were performed on the HR-ICPMS which is part of the PANOPLY platform (GEOPS-LSCE). O.P. benefits from a PhD Grant from the Paris-Sud University. Valérie Godard is thanked for thin-section production.

Supplementary Materials

Supplementary materials contains a sample description table as well as the detailed U-Pb method and data tables of the dating results. (Supplementary Materials)