{'fraction_non_alphanumeric': 0.0392006011593788, 'fraction_numerical': 0.02329492592857654, 'mean_word_length': 4.953030141655129, 'pattern_counts': {'":': 0, '<': 4, '<?xml version=': 0, '>': 3, 'https://': 0, 'lorem ipsum': 0, 'www.': 0, 'xml': 0}, 'pii_count': 0, 'substrings_counts': 0, 'word_list_counts': {'cursed_substrings.json': 1, 'profanity_word_list.json': 0, 'sexual_word_list.json': 0, 'zh_pornsignals.json': 0}}
Tectono‐sedimentary evolution of transverse extensional faults in a foreland basin: Response to changes in tectonic plate processes

Late Paleocene to Middle Eocene strata in the easternmost part of the Southern Pyrenees, up to 4 km thick, provide information on tectono‐sedimentary evolution of faults transversal to the Pyrenean chain. To know how changes in tectonic plate processes control the structural evolution of transverse faults and the synchronous thickness and lithological distribution of sedimentary strata in a foreland basin, field observations, interpretation of 2D seismic lines tied to lithostratigraphic data of exploration wells and gravity modelling constrains were carried out. This resulted in the following two tectono‐sedimentary phases in a foreland basin: first phase, dominated by transverse extensional faulting, synchronous with deposition of marine carbonates (ca. 57 to 51 Ma); and second phase, characterized by transverse contractional faulting, coeval to accumulation of marine and transitional siliciclastics (51 to 44 Ma). During the first phase, Iberia and Adria were moving to the east and west respectively. Therefore, lithospheric flexure in the easternmost part of the Iberian plate was developed due to that Sardinia was over‐thrusting Iberia. Consequently, activation of E‐dipping normal faults was generated giving rise to thick‐deep and thin‐shallow carbonate platform deposits across the hanging walls and footwalls of the transverse structures. During the second phase, a shearing interaction between Iberia and Sardinia prevailed re‐activating the transverse faults as contractional structures generating thin‐shelf and thick‐submarine fan deposits across the hanging walls and footwalls of the transverse structures. In the transition between the first and second phases, evaporitic conditions dominated in the basin suggesting a tectonic control on basin marine restriction. The results of our study demonstrate how thickness and lithology distribution, controlled by transverse faulting in a compressional regimen, are influenced by phases related to processes affecting motions and interactions between tectonic plates and continental blocks.

Several of the previous studies on transverse extensional faults discuss the change in stress directions during a single compressional event (periods of less than 20 Ma). However, the way in which changes in tectonic plate processes control the structural evolution of transverse faults and the synchronous thickness and lithological distribution of sedimentary strata in a foreland basin has yet to be considered. Addressing this issue should be relevant to consider the tectono-sedimentary evolution of transverse faults as key to establish kinematic histories of complex compressional zones. These histories are useful to better explain the present-day activity of seismic and volcanic regions with more than one orogenic salient; such as the Mediterranean region. The tectono-sedimentary evolution of transverse faults can also be used for understanding hydrocarbon and geothermal systems, as it can influence the distribution of thickness of source rocks by the generation of transverse depozones; as well as the quality of reservoirs and fluid migration by the formation of fractures and favouring pathways through these faults.
The Paleocene to Eocene succession of the easternmost part of the southern Pyrenees ( Figure 1a) provides an opportunity to investigate how multiphase transverse faulting in a foreland basin controls thickness and lithology distribution during tectonic plate processes. This area was located along the north-eastern margin of the Iberian Plate, which experienced continental collision with Eurasia, Corsica and Sardinia during the Eocene (Figure 1b) (Advokaat et al., 2014;Andreani, Loget, Rangin, & Pichon, 2010;Bestani, Espurt, Lamarche, Bellier, & Hollender, 2016;Lacombe & Jolivet, 2005). In the South-eastern Pyrenees, NNW-SSE striking faults are present (Figure 2a), transverse to the main W-E trend of this chain. These faults were active during the Early Eocene; interpreted as normal faulting, coeval to the Pyrenean compression (Estévez, 1970;Martínez et al., 1994;Santisteban & Taberner, 1979). However, there is still considerable uncertainty about the structural evolution, role on the stratigraphy and regional importance of these transverse faults.
In this study, we use a robust unpublished seismic reflection and well data set for the South-eastern Pyrenees in conjunction with field observations to: (a) characterize a foreland basin structure; (b) establish relations between transverse faults and thickness and lithology distributions through time and (c) record tectonic events. The Southeastern Pyrenees has been widely studied as a result of oil exploration. However, the number of boreholes is low and the quality of seismic data is poor due to the structural complexity and the existence of evaporite units (halite and anhydrite) with high-density contrasts (Figures 2b and 3). Therefore, a gravity analysis is integrated with the structural and lithostratigraphic data becoming a significant tool for validating the geological results. The aim of this work is to show how changes in tectonic plate motions and interactions control the structural evolution of transverse faults and the synchronous thickness and lithological distribution of sedimentary strata in a foreland basin. To achieve this, we study the relationship among transverse faults, thickness and lithology distributions through time with the geodynamic evolution of Iberia, Sardinia and Adria. By recording part of the geodynamic history of the Western Mediterranean region, we contribute to the knowledge of transverse extensional faults in foreland basins; specifically on 4-D structural and sedimentary evolution.

Highlights
• Two main tectonic phases are identified in the stratigraphic record of the South-eastern Pyrenees, related to the interactions between Iberia and Sardinia during Late Paleocene to Middle Eocene (57 to 44 Ma). • Tectono-sedimentary evolution of transverse faults in foreland basins can record changes of motion and interactions for tectonic plates and continental blocks. • The South-eastern Pyrenees is an exceptional area for interpreting the structural evolution of transverse extensional faults, which have been activated by the bending of foreland lithosphere surrounded by two orogenic salient.

| Sedimentary record
The South-eastern Pyrenees sedimentary succession ranges in age from Triassic to Pliocene and overlies a Paleozoic basement of granite and metamorphic rocks (Figures 2 and 3) Martínez et al., 1994;Muñoz et al., 1994;Vergés et al., 1994). During the early episodes of the continental collision, the sedimentary environments were characterized by non-and shallow-marine deposition. From the Early to Late Eocene, marine facies prevailed through an Atlantic Ocean connection. However, a period of isolation from the sea (not disconnect) from 50 to 46 Ma created conditions for the deposition of marine evaporites (Puigdefàbregas, Muñoz, & Marzo, 1986;Puigdefàbregas, Muñoz, & Vergés, 1992;Vergés, 1993;Vergés & Burbank, 1996).
For the purpose of this study, the most significant lithostratigraphic units are those of the Late Paleocene as well as the Early and Middle Eocene ages, which form a sedimentary cover up to at least 2,600 m in thickness ( Figure 3). The Late Paleocene to Early Eocene sequence is comprised of carbonates with a middle siliciclastic unit referred to as the Corones Formation. Overlaying this carbonate sequence is a substantial evaporite formation, known as the Serrat Evaporites, which closely marks F I G U R E 1 Regional setting. (a) Overview of the Western Mediterranean Region at present-day showing the location of the South-eastern Pyrenees (rectangle, Figure 2a), the main geographic features, structures (yellow lines) and basins (annotated Google Earth image). Structures are based on Vergés (1993) and Doglioni, Gueguen, Harabaglia, and Mongelli (1999). CCR, Catalan Coastal Range; NBFZ, Northern Balearic Fracture Zone. (b) Paleogeographic map (after Advokaat et al., 2014;Lacombe & Jolivet, 2005;Vacherat et al., 2017) showing the interactions between Iberia, Eurasia, Adria, Africa, Sardinia and Corsica during the Early Eocene. Paleogeographic information is based on Martínez et al. (1994), Christophoul et al. (2003), Carmignani et al. (2004), Andreani et al. (2010) and Costamagna and Schäfer (2017). Plate motions are extracted from Malusà et al. (2016) and Macchiavelli et al. (2017)  the boundary between the Early and Middle Eocene. The units deposited before the Serrat Evaporites are referred to herein as the Presalt group. Above these evaporites is the Vallfogona Formation, a unit dominated by siliciclastics and carbonates, and an upper discontinuous second evaporite unit known as the Beuda Gypsum Formation. Collectively, the Vallfogona and Beuda Formations will be referred to as the Suprasalt group. The stratigraphic sequence overlying this group is the Bellmunt sequence, comprised of siliciclastic units (Busquets, 1981;Calvet, Playà, Giménez-Montsant, & Permanyer, 2007;Carrillo, 2009;Carrillo, Rosell, & Ortí, 2014;Gimènez-Montsant & Salas, 1997;Martínez, Rivero, & Casas, 1997;Martínez et al., 2000;Puig et al., 2003;Puigdefàbregas et al., 1986Puigdefàbregas et al., , 1992).

| Structural features
The lithostratigraphic units are divided into two main structural zones: the autochthonous zone to the south and the allochthonous zone to the north. The boundary between these two zones is marked by the Middle-to-Late Eocene Vallfogona thrust (Figure 2a,b), a north-dipping frontal thrust (Cruset et al., 2018;Muñoz, Martínez, & Vergés, 1986;Ramos, Busquets, & Vergés, 2002).
The allochthonous zone is divided into two tectonic units known as upper and lower thrust sheets. The upper thrust sheets are located to the most western and eastern parts (Western and Eastern Upper Thrust sheets) of the South-eastern Pyrenees (Figure 2a), and are formed of sedimentary rocks ranging from Triassic to Late Eocene ages that were displaced between the Paleocene and Late Eocene (e.g. Martínez, Vergés, & Muñoz, 1988; Pujadas et al., 1989;Vergés, 1993;. The lower thrust sheets are characterized by the following subunits: the Cadí thrust sheet, whose major structure is a W-E trending fold known as Ripoll syncline; and the Serrat unit, underlying the Cadí thrust sheet, bounded by both a floor-androof thrust (e.g. Bello et al., 2008;Martínez et al., 1997;Muñoz et al., 1986;Souquet et al., 1975). Several studies have described three stratigraphic locations for regional décollement levels: the lower décollement, situated in the Late Cretaceous-Paleocene rocks Pujadas et al., 1989); the middle décollement, located in the Early Eocene units Muñoz et al., 1986) and the upper décollement, situated in the Serrat Evaporites (Carrillo, Koyi, & Nilfouroushan, 2017).
The Cadí thrust sheet deforms stratigraphic successions ranging from the Paleozoic (basement) to the Oligocene. This structural unit, displaced from the Middle to Late Eocene, is limited to the north by a regional backstop which was active before the Late Eocene ( Figure 2b) , 1989;Martínez et al., 1994;Muñoz et al., 1994;Ramos et al., 2002). The stratigraphic succession of the structural Serrat unit is incomplete (Bello et al., 2008;Vergés, 1993).
Both the autochthonous and allochthonous zones within the study area are affected by three major east-dipping faults ( Figure 2a) Muñoz et al., 1994;Pallí, Roqué, & Costa, 2011). In this study, from west to east, these faults are labelled as Western transverse fault (WTF), Central transverse fault (CTF) and Eastern transverse fault (ETF). It has been interpreted that some of these faults were generated as joints during the Late Carboniferous to Permian shearing stage and have been active since the Early Eocene (Estévez, 1970;Goula et al., 1999;Martínez et al., 1994;Santisteban & Taberner, 1979;Saula et al., 1994).

| Data set
This study uses 24 prestack time-migrated 2-D seismic reflection profiles (Figure 4a Geological map where the lithostratigraphic units, the main structures, exploration wells and stratigraphic and key sections are displayed (modify from Carrillo et al., 2017;Martínez et al., 2000;Muñoz et al., 1994;Pi et al., 2000). Exploration wells: A2, Ampurdan-2; B1, Bestrecà-1; B2, Banyoles-2; B4, Besalú-4; R1, Riudaura-1; R2, Riudaura-2; S1, Serrat- main structural trend) and W-E (parallel to the main structural trend). These profiles were acquired in 1985 by Unión Texas España Inc. and reprocessed by CEPSA. The seismic lines are presented in two-way time (TWT), having a sample rate of 4 ms and a record length of 6 s. An estimated average interval velocity of the sedimentary cover of ca. 4,850 m/s is derived from checkshot data from wells. Although the frequency of the peak amplitude in the time-frequency gathers is nowadays not available, considering the source (dynamite) and the frequency of the geophones (10 Hz) we can assume that the dominant frequency is about 30 Hz. This means that the vertical resolution of the seismic lines is of the order of 40 m. To reflect the relations between transverse faults and thickness and lithology distributions through time, structural and stratigraphic features observed are herein described using seconds in TWT, and meters considering this velocity and applying the velocities shown in Table 1 respectively.
Field descriptions and measurements (stratigraphic and structural features) were collected during field campaigns. This information enabled the updating of previous geological maps Mató et al., 1996;Muñoz et al., 1994;Pi et al., 2000) and the generation of eight key lithological sections (sections 1 to 8) in the Cadí thrust sheet (Figure 4b). Lithological information from the Serrat-1 and Bestrecà-1 wells was used to complete two of these sections (sections 2 and 4 in Figure 4b respectively).
A detailed gravity survey was conducted with a total of 844 data points measured ( Figure 4a) with a Lacoste-Romberg gravity meter, model G831, and referred to the IGSN-71 through the Spanish Gravity Net. Spatial positioning and height of the stations were obtained from global positioning systems (GPS) and benchmarks with an elevation precision of ±0.1 m. The new gravity measurements were integrated with data from previous surveys (Martínez et al., 1997;Rivero, 1993;Rivero, Pinto, & Casas, 2002).

| Structural and stratigraphic framework
Four main and relevant horizons were interpreted from all seismic profiles based on relationships obtained from surface and subsurface stratigraphy, VSP and synthetic seismograms ( Figure 5). These horizons are as following: (a) boundary between the Basement and overlying sedimentary cover; (b) top of the Presalt group; (c) top of the Serrat Evaporites and (d) top of the Suprasalt group. Another useful surface that provides a stratigraphic control of tectonic events through the time is the top of the Corones Formation, which was used to differentiate the Presalt group into two sub-groups: the Lower Presalt and the Upper Presalt. To match the structural interpretation between the autochthonous and allochthonous zones, 2-D and 3-D surface horizons and fault restorations were performed. According to Carrillo et al. (2017), a shortening of 30% and 15% of the Cadí thrust sheet and the structural Serrat unit, respectively, was taking into account for these restorations.
To highlight fault growth in 3-D and evaluate the impact of the transverse structures on sediment thickness in time and space, isochron and isopach maps (Figures 6 and 7) and a lithostratigraphic well correlation (Figure 8a) of the main horizons were generated based on seismic interpretation and surface mapping. In the case of the Suprasalt group and its overburden (the Bellmunt sequence), interval velocity is well constrained in the Cadí thrust sheet and the Empordà Basin, and it was converted from time to depth. An average velocity of 3,900 m/s, which corresponds to the average velocity of this group and overburden (Table 1), was used for this conversion. To identify the impact of transverse structures from field evidence, the eight lithological sections were correlated with a W-E orientation ( Figure 8b). Observations from detailed seismic profiles and outcrops ( Figure 9) were incorporated into our study, supporting the structural and stratigraphic interpretations.

| Gravity modelling
All the combined gravity data were reduced using the classical formulae of the Bouguer anomaly, where a series of corrections were applied to eliminate the non-geological causes related to gravity variations, including topographic correction. The Bouguer anomaly values were then interpolated by kriging to a 0.5 km × 0.5 km square grid and contoured.
As the Bouguer anomaly map integrates the effect of both long-and short-wavelength components, the regional factor was removed from the Bouguer anomaly to obtain a residual anomaly map. This last map was assimilated to a second-order surface whose orientation is consistent with the gravity map of the Pyrenees (Casas, Kearey, Rivero, & Adam, 1997).
To understand residual gravity anomalies and constrain the structural and stratigraphic interpretations, three representative cross-sections were selected and converted from time (in TWT) to depth ( Figure 10). The lithologies and interval velocities used for the conversion are listed in Table 1. The residual gravities of these sections were calculated, obtaining inversion models and compared with the measured ones (from the residual gravity map). The densities, assumed for the calculated gravities (Table  1), were extracted from previous works (Carrillo et al., 2014;Martínez et al., 1997;Rivero, Vilas, Pinto, & Casas, 2001).

| Basin and regional framework
The isopach maps were integrated with the stratigraphic correlations and results of the detailed seismic and outcrop interpretations as well as the results documented in previous works to produce three novel paleogeographic maps ( Figure 11). These maps correspond to tectono-sedimentary stage models. Each map shows structural features and depositional lithofacies, integrating both the autochthonous and allochthonous zones and accounting for the fore-mentioned shortenings established by Carrillo et al. (2017).
The paleogeographic maps were compared with regional tectonic studies to identify relationships between thickness and lithological distributions of sedimentary strata, influenced by transverse faults, and tectonic plate processes, such as Iberia and Sardinia. To illustrate these relationships, we constructed a transorogenic (350-km-long) cross-section at the scale of the lithosphere (Figure 12). This section shows the geodynamic evolution of Iberia and Sardinia during part of the continental collision stage (Late Paleocene to Middle Eocene), based on previous works, and the tectono-sedimentary features obtained in this study. Finally, to highlight the contribution of our work, a comparison between the tectono-sedimentary evolution of the studied transverse faults and similar structures in other regions was achieved.

| STRUCTURE
In this study, it is observed that the major transverse structures affect the Basement and sedimentary succession ranging from the Presalt group to the Bellmunt sequence in both the autochthonous and allochthonous zones with a high of at least F I G U R E 6 Structural maps in the autochthonous zone for the top Basement (a), in depth-time referred to a seismic datum, and the autochthonous and allochthonous zones for the top Serrat Evaporites (b), in height referred to the present-day sea level, where the main faults and folds (black and dashed white lines), the present-day location of the Vallfogona thrust (red lines) and exploration wells (circles) are displayed. Exploration wells: A2, Ampurdan-2; B1, Bestrecà-1; B4, Besalú-4; R1, Riudaura-1; R2, Riudaura-2; S1, Serrat-1; V1, Vallfogona-1 2.00 s TWT (ca. 4,850 m) (Figures 4b and 5). In the case of the WTF, the upper fault tip is located up to the lower part of the Bellmunt sequence (Figure 5a,b). These faults divide the study area into structural blocks ( Figure 6), with footwalls and hanging walls as well as secondary structures which are irrelevant in this study. To understand fault growth and thickness

| Western transverse fault (WTF)
In the Ebro Basin, the WTF is located eastward of the Riudaura-1 well ( Figure 4b). Here, this structure has a maximum throw of 0.45 s (ca. 1,100 m) at the top Basement with a normal geometry (Figure 5a). Based on a 3-D geometry of the seismic horizons, it is observed that this structure continues north below the Vallfogona thrust with a maximum throw at the top Basement of 0.83 s (ca. 2,010 m) and a normal geometry (Figure 6a). In the autochthonous zone, the length of the WTF is at least of 15 km.
In the Cadí thrust sheet, a significant NNW-SSE fault, exposed in the northwest region of the study area (Figure 4b), has been described in previous works (Martínez, Carrillo, Tallada, & Copons, 2015;Muñoz et al., 1994). Based on restorations, this fault links to the WTF in the autochthonous zone. Here, this fault has a reverse and dextral movement with a maximum throw at the top Serrat Evaporites of 0.37 s (ca. 900 m) (Figure 5e). To the southeast of Serrat-1, near the subsurface axis of the Ripoll syncline, the WTF is theorized to link to blind faults displaying NW-SE strikes ( Figure  6b) with normal displacements and throws at the top Serrat Evaporites of up to 0.13 s (ca. 315 m) (Figure 5b,d). The length of the WTF through the Cadí thrust sheet is at least of 12 km.

| Central transverse fault (CTF)
In a similar manner to the WTF, the CTF in the allochthonous zone fits with an ENE-dipping fault below the Vallfogona thrust (Figures 5d and 6a). This second fault corresponds to the CTF in the autochthonous zone having a normal geometry and throw at the top Basement of up to 0.54 s (ca. 1,310), which decreases northward. In the autochthonous zone, the length of the CTF is at least of 5 km.
In the Cadí thrust sheet, the CTF is superimposed by a NNW-SSE trending anticline with a length of at least of 12 km (Figures 5b and 6b). This anticline affects the stratigraphic successions ranging between the Presalt group and the Bellmunt sequence, including the sequences of the structural Serrat unit. Geometrically, the fold has a low curvature and its west limb is steeper than the east limb. In the western limb, the CTF is observed with a reverse movement and a throw at the top Serrat Evaporites of 0.15 s (ca. 360 m). In the northern part of the Cadí thrust sheet, the CTF steps westward and displays a normal kinematic sense (Figure 5d,e).

| Eastern transverse fault (ETF)
In the northern part of the ETF, the structure displays a reverse displacement of horizons in the autochthonous zone with a throw at the top Basement of 0.58 s (ca. 1,410 m) (Figure 6a). To the south, this fault has a normal geometry with a maximum throw of 0.88 s (ca. 2,135 m) (Figure 5c).

| Structural blocks
All the structural blocks in the study area contain low relief (<0.13 s depth; ca. 315 m) basement folds and faults. In the autochthonous zone, the footwall block of the WTF is mainly characterized by a NW-verging monocline (Figure 6a). By contrast, both the footwall and hanging wall blocks of the CTF present W-verging half-grabens. The hanging wall block of the ETF contains a S-verging graben monocline and an E-W trending blind fold-thrust, extending almost until the ETF (Figure 5c).
In the Cadí thrust sheet, the footwall of the WTF shows the geometry of the Ripoll syncline with a horizontal fold hinge (Figure 6b). Between the WTF and the CTF, the Serrat Evaporites are affected by a NNW-SEE trending depression (Figure 5d). Moreover, in the structural Serrat unit, an evaporite dome is noted in the north-eastern sector of the related block.
From the CTF to the ETF, two structural areas are recognized: the west side with a depression and the east side with a high.

LITHOLOGICAL DISTRIBUTION
Thickness and lithology variations in the lithostratigraphic units are identified in both the autochthonous and allochthonous zones ( Figure 5). These distributions vary markedly adjacent to the transverse faults as well as the low relief basement folds and faults (Figures 7 and 8). The variations are described separately for the Presalt group, the Serrat Evaporties and the Suprasalt group.

| Presalt group
In the autochthonous zone, the thickness of the Presalt group increases across the footwalls to the hanging walls of the transverse faults (Figure 7a). In the footwall of the WTF, a thickening of this group to the north-western sector is observed, varying from 0.14 to 0.32 s (ca. 370-840 m) thick. In the footwall of the CTF, a depocenter is recognized adjacent the WTF. This depocenter has a maximum thickness of 0.51 s (ca. 1,340), which dramatically decreases southward to 0.18 s where the throw of the WTF is less (Figure 6a). Toward the CTF, the thickness is reduced to 0.10 s (ca. 260 m). In the hanging wall of the CTF, the Presalt strata thickness is 0.39 s (ca. 1,020 m), which decreases slightly to the centre of the block. In the north-eastern sector of this block, a depocenter of up to 0.38 s thick is noted on an axial trace of a SW-NE trending basement syncline. In the southern sector, the Ampurdan-2 well crossed the entire Presalt Group (441 m thick), on a structural high (Figure 8a). In the hanging wall of the ETF, from outcrop the Presalt thickness is 2,200 m (stratigraphic section 8, Figure 8b). The greatest subsurface thickness in the study area (up to 1.10 s; ca. 2,670 m) is observed within the northern and central sectors of this block (Figures 5c and 7a). However, this thickness estimate is affected by structural thickening (up to ca. 30%), related to a blind fold-thrust (Figure 5c). The thickness decreases towards the southern sector of the block, where the throw of the ETF is less (Figure 6a).
In the Cadí thrust sheet, the thickness of the Lower Presalt increases greatly across the CTF and ETF from the footwall to hanging wall (Figure 8b). The same occurs for the Upper Presalt across the ETF; conversely, the thickness of the Upper Presalt decreases across the CTF. The structural block between the WTF and the CTF displays westward thickening of the Presalt group from 1,340 to 1,980 m and 1,110 to 1,340 m for the structural block between the CTF and ETF.
The thinnest (<500 m thick) successions of the Presalt group are dominated by about 80% of limestones and 20% of marls and sandstone in both the Ebro and Empordà basins and the allochthonous zone (Figure 8). On the contrary, the thickest (>500 m thick) successions are formed of 50% marls and shales and 50% of limestones and sandstones.

| Serrat Evaporites
In the autochthonous zone, the thickness of the Serrat Evaporites increases across the footwalls to hanging walls of the CTF and in the southern part of the ETF (Figures 5b and  7b). In the footwall block of the WTF, the average thickness of the Serrat Evaporites is 0.20 s (ca. 520 m). However, the thickness increases to 0.27 s along the Serrat Evaporites-Vallfogona thrust contact, probably due to structural thickening (salt tectonics). The footwall of the CTF shows a peculiar E-W F I G U R E 1 2 Crustal-scale cross-sections across the South-eastern Pyrenees from the low-relief thrust system, in Iberia, to the basement high, in Sardinia, during the first (a) and second tectonic phases (b) showing the interactions between Iberia and Sardinia from Late Paleocene to Middle Eocene. The sections illustrate the study area's lithological, paleogeographical and paleoenviromental context as well as main structures. Apatite U-Th/He (AHe) and apatite fission track (AFT) ages are also show (from Rushlow et al., 2013 andMalusà et al., 2016). In the Sardinian domain, the geological information is extracted from Carmignani et al. (2004) and Costamagna and Schäfer (2017) trending "salt" wall with thicknesses up to 0.34 s (ca. 885 m) below the Vallfogona thrust (Figure 6a). To the north and south of this wall, the evaporitic succession is thinnest (Figure 5b,d) with, locally, low values of 0.05 s (ca. 130 m) ("salt" welds). In the footwall block of the ETF, the Serrat Evaporites unit thickens from 0.1 s (ca. 260 m) thick in the central region to 0.25 s (ca. 650 m) thick towards the CTF. In the south-eastern sector of this block, near the contact between the Serrat Evaporites and the Vallfogona thrust, the thickness is around 0.29 s (ca. 750 m). However, eastward of the Besalú-4 well, the thickness decreases to 0.07 s (ca. 180 m) on a basement high (Figure 5c). The unit thicknesses in the Ampurdan-2, S-43 and Banyoles-2 wells are 240, 160 and 220 m respectively (Figure 8a). In the hanging wall of the ETF, the evaporites are not recognized outcropping northward (section 8 in Figure 8b), where the ETF has a reverse geometry (Figure 6a). On the other hand, to the south, a thin (0.07 s thick) succession of the Serrat Evaporites is recognized with erosional truncations (Figure 9a). This succession increases in thickness southward, up to 0.37 s (ca. 960 m) thick, where a significant depocenter is present.
In the allochthonous zone, different areas of prominent thickness of the Serrat Evaporites are observed (Figure 5b,d).
In the Cadí thrust sheet, these areas are located in the north limb of the Ripoll syncline, westward between the WTF and CTF, with at least 0.17 s (300 m thick; Figure 8b), and in the eastern part between the CTF and ETF, at least 0.20 s (ca. 520 m) thick. By contrast, thin (<100 m thick) successions of the Serrat Evaporites are recognized eastward of the footwalls related to the WTF and CTF. In the structural Serrat unit, an area with a marked thickness of at least 0.78 s (ca. 2,030 m) is observed on the footwall of the WTF ( Figure  5e). This thickness decreases progressively in the footwall of the CTF. In the anticline superimposing this structure, another area with a significant thickness of at least 0.51 s (ca. 1,330 m) is noted (Figure 5b).
The lithology distribution of the Serrat Evaporites has the following features: (a) in the Ebro Basin, anhydrite and carbonate layers dominate the structural highs ( Figure 8a

| Suprasalt group
In the autochthonous zone, the thickness of the Suprasalt group increases from the footwall to hanging wall across both the WTF and ETF (Figures 5a and 7c). In the footwall of the WTF, 100 m of the Suprasalt group was measured in the Riudaura-1 well. Adjacent the WTF hanging wall, a thickness up to 0.36 s (ca. 860 m) is identified (Figure 5b). Just east of the CTF, within the CTF hanging wall, this group reaches at least 0.25 s (ca. 600 m). This thickness decreases eastward, adjacent the ETF, ranging from 80 to 150 m thick ( Figure  8a). In the hanging wall of the ETF, the Suprasalt group displays up to 220 m thick in the northern part (Figure 8b). There, in some seismic profiles, reflectors corresponding to the Suprasalt group truncate the Serrat Evaporites (Figure 9a). In the south-eastern part of the ETF hanging wall, this group has a depocenter of up to 0.20 s (ca. 480 m) thick (Figure 5c).
On the northern limb of the Ripoll syncline in the Cadí thrust sheet, the thickness of the Suprasalt group decreases from the footwalls to hanging walls across the WTF and CTF (Figure 8b). By contrast, the thickness increases across the faults on the synclinal southern limb (Figure 9b). Therefore, in the north-eastern and south-western parts of the related structural blocks, depocenters ranging from 500 to 1,000 m thick for the Suprasalt group are identified (Figure 7c). While, in the north-western and south-eastern parts, the thickness decreases abruptly with values lower than 500 m. Seismic reflectors of the Suprasalt group onlaping the Serrat Evaporites are recognized in the south-western depocenter towards the transverse faults (Figure 9b).
The thinnest (up to 250 m thick) successions of the Suprasalt group are mainly dominated by carbonate, siltstone and sulphate layers (Figure 8). These successions are identified in the Ebro and Empordà basins as well as the local eastern part of the Cadí thrust sheet. By contrast, the thickest (>250 m thick) successions are formed of siltstone, sandstone and sulphate layers.

| GRAVIT Y CONSTRAINTS
In the residual gravity map, values ranges from 13 to −10 mGal. A series of significant anomalies and variations on residual gravity with NNW-SSE, NW-SE and NE-SW directions are observed. In the inversion models (Figure 10), part of these variations are located around the transverse faults. The models are referred to herein as "Central model" (Figure 10a), related to Figure 5b, "Eastern model" (Figure  10b), for Figure 5c, and "Northern model" (Figure 10c), linked to Figure 5d. While small variations in the thicknesses and densities used in the models can render similar residual gravity trends, no alternative structural interpretation has been found that matches the information available. sector. We relate this anomaly to the existence of a 1.8-kmthick deposit of pure anhydrite (>80% in anhydrite content with an average density of 2.90 g/cm 3 ) of the Serrat Evaporites, forming part of a relative thin (up to 3.0 km thick) sedimentary cover. Residual gravity decreases to −1 mGal at the axis of the Ripoll syncline, indicating a thickening (up to 4.1 km) of this sedimentary cover and thinning of the anhydrite deposit. Residual gravity increases up to 1 mGal at the CTF, signifying a thin sedimentary cover with thickening of the anhydrite deposit. From the CTF towards the Vallfogona thrust, residual gravity initially decreases to −1 mGal, however, continuing east it increases to a positive anomaly of up to 5 mGal. The negative value is attributed to a salt body (2.10 g/cm 3 ) up to 0.8 km thick in the structural Serrat unit. The positive anomaly is associated with thickening (up to 2.1 km) of the anhydrite deposit within a thin (up to 2.5 km thick) sedimentary cover. Eastward of the Vallfogona thrust, a negative residual anomaly of up to −8 mGal is recognized. We relate this anomaly to the following three features: (a) presence of the thinnest (0.7 km thick) sedimentary cover; (b) low thickness (<0.4 km thick) of the anhydrite deposit and (c) thickening of siliciclastics (2.45-2.60 g/cm 3 ) of the Bellmunt sequence.

| Eastern model
In the Eastern model (Figure 10b), between the Ampurdan-2 well and the intersection point with Figure 10a, the residual gravity gently reduces from −7 to −9 mGal. These negative values are associated with the same three features listed for the easternmost section of the Central model. These facts, combined with the Ampurdan-2 lithologic descriptions (Figure 8a), suggest a change in the basement lithology, from a high to low concentration of schist and increment of granite (2.72 and 2.64 g/cm 3 respectively). From the Central model intersection point to the hanging wall of the ETF, the residual gravity rises abruptly up to −3 mGal. This change is mainly due to two features: (a) a variation in both lithology and thickness of the Serrat Evaporites, from anhydrite and carbonate (2.83 g/cm 3 ) at ~0.4 km thick to pure anhydrite at ~1.0 km thick; and (b) a variation in lithology of the Basement, from granite to schist. Towards the north-eastern end of the Eastern model, residual gravity increases up to −2 mGal due to reduced thickness of the siliciclastic Bellmunt sequence.

| Northern model
In the Northern model (Figure 10c), from southeast to northwest, the residual gravity increases from 0 to 4 mGal in a NW direction. This positive trend is attributed to the existence of a thick deposit (1.0 km) of pure anhydrite in the frontal part of the structural Serrat unit. Passing northeast from the Cadi thrust sheet to the CTF, residual gravity reduces to −8 mGal, indicating a decrease in anhydrite content in the passage towards a more shale and salt prone section in the structural Serrat unit. A salt dome was added to the model to explain the lowest residual gravity point. Continuing northeast, the residual gravity rises to −5 mGal, suggesting a thick deposit (1.2 km) of pure anhydrite below the Vallfogona thrust.
Residual gravity values then drop up to −8 mGal, which is associated with a reduction in anhydrite thickness.

EVOLUTION OF THE TRANSVERSE FAULTS
We propose three paleogeographic maps for the Late Paleocene to Middle Eocene tectono-sedimentary evolution of the main transverse faults addressed in this study (WTF, CTF and ETF) ( Figure 11). In addition, these maps include the evolution of N, NNE, NW and NE dipping faults also identified in the study area. The N and NNE dipping faults, located in the northern and north-eastern parts of the study area, have been interpreted in previous works (Martínez et al., 1989;Pujadas et al., 1989). The maps are synchronous to the sedimentation of the Presalt group, Serrat Evaporites and Suprasalt group, corresponding to the following sedimentary stages: Stage 1 (Figure 11a

Eocene (57 to 51 Ma)
Thickening of the Presalt strata across the footwalls to the hanging walls of the transverse faults (Figures 7a and 8) suggests that these structures worked as normal faults during Stage 1 (Figure 11a). This fact is in agreement with Estévez (1970) for the ETF. Southward thinning of the strata along these structures indicates changes in throw.
Based on the relationships between thickness distribution, as controlled by extensional transverse faulting, and lithology of the Presalt strata (Figures 7a and 8), the following tectono-sedimentary features are interpreted for the present stage: limestones, deposited in structural highs (footwalls); and shales and marls, deposited in structural depressions (hanging walls) (Figure 11a). In agreement with previous works (Gimènez-Montsant & Salas, 1997; Martínez et al., 1988), the shale and marl correspond to deep platform environments, and the limestone to shallow platform. In turn, this study interprets that NNW-SSE depositional-environmental belts were formed within each of the structural blocks, characterized by deep platform environments to the west and shallow platform environments to the east. According to Gimènez-Montsant and Salas (1997), the north-eastern part of the study area was a shallow detrital environment, attributed to a delta plain deposited on an uplifted eastern margin.

| Stage 2: Early to Middle Eocene (51 to 49 Ma)
During the second stage, a westward thickening of the Serrat Evaporites between the WTF and CTF in the Cadí thrust sheet (Figure 8b) suggests that normal faulting was active in the northern portion of the WTF (Figure 11b). Determining topographic and tectonic configuration along the southern portion of the WTF is problematic, due to the high degree of present-day deformation observed for the Serrat Evaporites along this structure (Figure 7b). However, normal faulting is also proposed for the southern sectors of the CTF and ETF as they demonstrate the expected stratigraphic thickening across the footwall to the hanging wall blocks in the autochthonous zone ( Figure 10). On the other hand, the northward thinning along the ETF hanging wall and the unconformities within the Serrat Evaporites ( Figure  9a) indicate contractional faulting along the northern portion of the ETF. In the case of the salt dome modelled in the structural Serrat unit between the WTF and CTF ( Figure  10c), we interpret that this dome was transported from the prekinematic hanging wall of the CTF to the present-day allochthonous zone.
Based on the lithological distribution of the Serrat Evaporites (Figures 8 and 10) and the paleo-topography, influenced by faulting, the following tectono-sedimentary features are described for the present stage: anhydrite, limestone and dolostone, accumulated in structural highs; and halite or shale, deposited in structural depressions (Figure 11b). In agreement with Carrillo et al. (2014), the anhydrite, carbonate and dolostone correspond to sulphate-carbonate shelves, pure anhydrite to selenitic wedges or basins and the halite to salt deep basin. Therefore, in this study, it is interpreted that selenitic and salt deep basins were mainly concentrated along the hanging walls adjacent the active transverse extensional faults.

| Stage 3: Middle Eocene (49 to 44 Ma)
Thickening of the Suprasalt strata across the footwalls to the hanging walls of the transverse faults in the autochthonous zone and the southern part of the Cadí thrust sheet ( Figure 7c) suggests that the southern portion of these structures worked as normal faults during Stage 3 (Figure 11c). However, northern and north-western thinning of the units within the same blocks (Figure 8b) indicates contractional faulting along the northern portions of the transverse structures. In the case of the ETF, the erosional truncations of the Suprasalt group on the Serrat Evaporites in the northern sector of the Empordà Basin (Figure 9a) suggest that the contractional zone migrated southward from that in Stage 2, and was located further south than the contractional zones of the WTF and CTF in Stage 3.
Based on the relationships between lithology of the Suprasalt strata ( Figure 8) and thickness distribution, as controlled by transverse faulting, the following tectono-sedimentary features are interpreted for the current stage: carbonate, siltstone and sulphate layers, deposited in structural highs (footwalls); and siltstone, sandstone and sulphate layers, deposited in structural depressions (hanging walls) ( Figure 11a). In agreement with previous works (Carrillo et al., 2014;Costa, 1989), the carbonate, siltstone and sulphate layers correspond to shelves, and the siltstone, sandstone and sulphate layers to slope/submarine fans. Thus, in this study, it is interpreted that slope/submarine fans concentrated in the northern and south-western parts of the structural blocks.

| Relationships between tectonosedimentary evolution of transverse faults and tectonic plate processes
The results described above, collaborated with descriptions of eastern Pyrenees tectonic processes, open a discussion on how changes in tectonic plate motions and interactions control the structural evolution of transverse faults and the synchronous thickness and lithological distribution of sedimentary strata in a foreland basin.
The structural evolution of the transverse faults and the synchronous thickness and lithology distribution across these structures, as observed in this work (Figure 11), indicates a progressive change from stretching to contractional mechanisms migrating from north to south along faults and east to west across structural blocks. This change occurred during the ending of the Early Eocene (ca. 51 Ma) as a response of tectonic processes affecting Iberia and Sardinia. To evaluate this response, we present a cross-section ( Figure  12) which assumes the following two points: (a) Sardinia and Corsica were an independent continental block, not a part of Iberia; and (b) during the Eocene, the southernmost part of Sardinia was tectonically interacting with the easternmost part of Iberia. These points are in agreement with previous works ( CARRILLO et AL. et al., 2016;Horner & Lowrie, 1981;Lacombe & Jolivet, 2005).
The present-day thickness of the Iberian continental crust is characterized by eastward thinning, from ca. 45 km in the central Pyrenees to ca. 20 km in the Empordà Basin (Chevrot et al., 2018). This thinning is due to the Miocene to Pliocene extensional stage, which displaced Sardinia with a counter-clockwise rotation to the present-day position (e.g. Roca, Sans, Cabrera, & Marzo, 1999). The two continental blocks of Sardinia and Corsica have a maximum crustal thickness of 34 km (Egger, Demartin, Ansorge, Banda, & Maistrello, 1988), and they thin up to 25 km along the margins (Gailler, Klingelhoefer, Olivet, Aslanian, & Sardinia-group, 2009;Prada et al., 2013). According to Bestani et al. (2016), the south-eastern part of Eurasia and Corsica had maximum crustal thickness of 60 km during the Eocene. Therefore, in this study, it is assumed that during the same epoch, the easternmost part of the Iberian as well as the Sardinian crust would have been thicker than the present day, between 45 and 60 km.
According to previous works (Macchiavelli et al., 2017;Malusà, Danišík, & Kuhlemann, 2016), two tectonic phases related to motions of Adria and Iberia, with respect to Eurasia, are distinguished from the Late Paleocene to Middle Eocene. These phases are illustrated in Figure 12, where the first phase takes place from 57 to 51 Ma and a second phase between 51 and 44 Ma. The relationships between these phases and the tectono-sedimentary evolution of the South-eastern Pyrenees and southern part of Sardinia are discussed below.

| First phase (57 to 51 Ma)
During the first phase, Iberia was moving east and Adria to the west, relative to Eurasia (Macchiavelli et al., 2017;Malusà et al., 2016). Consequently, it is known that during the Early Eocene, Sardinia and Corsica were over-thrusting Iberia and Eurasia forming an active N-S striking mountain range (Figure 12a) (Andreani et al., 2010;Bestani et al., 2016;Lacombe & Jolivet, 2005). It is assumed in this study that the boundary between the Iberian plate and the Sardinian block was marked by a main thrust known as the present-day Northern Balearic Fracture Zone (NBFZ; Figure 1a).
In the Sardinian block, we identify three domains which prevailed during the first phase, from west to east: (a) a high-relief thrust system, verging to the east and deforming basement units and a sedimentary Mesozoic cover; (b) a Tethyan-influenced marine piggy-back basin and (c) a basement high with a sedimentary Mesozoic cover ( Figure  12a). Apatite U-Th/He (AHe) with cooling ages between 80 and 57 Ma in the southernmost part of Sardinia, analysed by Malusà et al. (2016), and E-verging thrusts affecting Mesozoic series (Barca & Costamagna, 1997) support the existence of the first domain. To the north, previous works (Carmignani, Funedda, Oggiano, & Pasci, 2004;Costamagna & Schäfer, 2017) recognize Early Eocene shallow marine carbonate deposits indicating the influence from the Tethys sea within the second domain. In the central-east of Sardinia, thermochronological interpretations by Zattin, Massar, and Dieni (2008) suggest uplifting from 140 Ma to Oligocene, supporting the basement high of the last domain.
In the Iberian plate, we also recognize three domains during the first phase, from east to west: (a) a high-relief W-verging thrust system, deforming basement units and a sedimentary Mesozoic cover; (b) an Atlantic-influence marine foreland basin controlled by N-S striking normal faults and (c) a low-relief thrust system, verging to the east and involving sedimentary Mesozoic cover (Figure 12a). West verging structures involving metamorphic and Mesozoic units have been well-documented (e.g. Carreras, 2001;Druguet, 2001;Fleta et al., 1994), supporting the existence of the first domain. We propose that during the first phase a load and lithospheric flexure in the easternmost part of the Iberian plate developed due to the Sardinian over-thrusting. Consequently, activation of E-dipping normal faults for both pre-and early orogenic structures was generated, giving rise to thick/deep and thin/shallow carbonate platform deposits across the hanging walls and footwalls of the transverse structures. Also, a forebulge transversal to an active W-E Pyrenean range was formed in the central part of this basin. The paleogeographic features of our Stage 1 (Figure 11a) are consistent with this proposal. Apatite fissions tracks (AFT) with cooling ages between 59 and 48 Ma for Mesozoic units in the Western Upper Thrust Sheet (Figure 2a), documented in Rushlow, Barnes, Ehlers, and Vergés (2013), suggest the low relief in the last Iberian domain.
In the Sardinian block, molasse facies (alluvial fan and fluvial systems) in the south, documented by Costamagna and Schäfer (2017) with a Middle Eocene age, suggest that the uplifting of the Sardinian high-relief thrust system persisted, and the piggy-back basin disconnected from the Tethys sea (Figure 12b). In the basement high domain, Middle Eocene deformation affected the Mesozoic cover with minor fold and thrusts (Arragoni, Maggi, Cianfarra, & Salvini, 2016).
In the Iberian plate, syn-orogenic conglomerates in the frontal parts of the Eastern and Western upper thrust sheets Pi et al., 2000), indicate that the high-and low-relief thrust systems continued their uplifting to at least the Late Eocene (Figure 12b). A change from extensional to contractional kinematic of the transverse faults occurred in the northern part of the foreland basin, giving rise to thin/shelf deposits in the hanging walls as it is observed in our Stages 2 and 3 (Figure 11b,c). By contrast, thick/submarine fan deposits were generated in the footwalls. We interpret that the structural change from extensional to contractional was due to the new displacement of Iberia to the northwest, where the active Pyrenan chain was acting as backstop. The east-to-west migration of contraction is consistent with the existence of Eulerian counter-clockwise poles to the central part of Iberia (Tavani et al., 2018 and referred herein). The transition between the first and second phases is also marked by a lithological change from carbonates (Presalt Group) to evaporites (Serrat Evaporites and the Suprasalt Group) and siliciclastics (Suprasalt Group) in the foreland basin suggesting a tectonic control on basin marine restriction and sedimentary conditions.

| Comparisons with other regions
Geometries and stratigraphic variations through transverse faulting in foreland basins have been reported in the central and southern Apennines (Doglioni, 1995;Tavani et al., 2015), the Eastern Maghrebides (Bianca et al., 1999;Billi et al., 2006;Gutscher et al., 2015;Torelli et al., 1998) and the Carpathian Bend Zone (Tărăpoancă et al., 2003). All these examples acted with two stages of different stress direction: a first extensional stage; and a subsequently reverse and/or strike-slip motion (e.g. Gutscher et al., 2015;Tărăpoancă et al., 2003;Tavani et al., 2015). The main tectonic and stratigraphic features of these foreland basins are highlighted and compared with our case study.
In the Apennines, Tavani et al. (2015) have provided a kinematic evolution of a fault system transversal to a paleo-subducting front. The same authors have interpreted that the transverse extensional faulting was attributed as a response of a syncline arching and flexure of a subducted plate forming a non-cylindrical forebulge. The sedimentary succession deposited during the motion of the fault system is only formed of siliciclastics. Fault-throws at the prefaulting horizons and synchronous stratigraphic thickness have not been reported in the Apennines.
In the Carphatian Bend Zone, Tărăpoancă et al. (2003) have shown foreland deposits within hanging walls of faults transversal to a thrust front. This work displays thickness variations, although it is based only on seismic data.
In the Eastern Maghrebides, Torelli et al. (1998) have recognized a stratigraphic record deposited coevally with normal faulting transversal to the Maghrebian Thrust Belt. A forebulge, transversal to this front, has also been identified (Bianca et al., 1999;Billi et al., 2006;Torelli et al., 1998). According to Billi et al. (2006), the normal faults were developed due to a flexure of a subducted plate associated with a lateral loading effect by a growing accretionary wedge. Thickness variations in 2D across the faults and a diversity of lithologies have been recognized in the related foreland basin (Hyblean basin; Torelli et al., 1998). However, this information is based on off-shore subsurface data.
From a tectonic point of view, we consider that the generation of the transverse extensional faults in the Eastern Maghrebides (described above) is analogue to the studied faults in this work. Furthermore, again similarly to the Eastern Maghrebides, the South-eastern Pyrenees underwent lateral loading from an additional orogenic salient (the Sardinian block; Figure 12a) forming flexural extension. By contrast, the cases from the Apennines and the Carphatian Bend Zone were controlled by a single orogenic salient in an along-strike stretching setting.
Our case study in the South-eastern Pyrenees provides the complete structural, stratigraphic and lithologic features related to a tectono-sedimentary evolution of transverse extensional faults in a compressional regimen. Here, we have an on-shore case which is supported by field, seismic and well data. Moreover, thickness distribution patterns and lithology variations (carbonates, evaporites and siliciclastics) during the structural evolution are recognized. In Tavani et al. (2015), the orientation of the fault dips and the impact of the structural evolution on thickness and lithology distribution remains uncertain. In Tărăpoancă et al. (2003), the relationships between lithology variations and structural evolution of transverse faults, supported by field and/or well evidences, have not been provided. In the Eastern Maghrebides, field examples have not been reported yet. All of these remaining features make the South-eastern Pyrenees an exceptional area to understand the 4-D structural and sedimentary evolution of transverse extensional faults in foreland basins.

| 1409
1. Two main tectono-sedimentary phases can be distinguished: (a) first phase with deep and shallow marine carbonate accumulation controlled by extensional faulting, transverse to an active Pyrenean chain, synchronoulsy to east displacement of Iberia and frontal collision of this plate with Sardinia (ca. 57 to 51 Ma); and (b) second phase with marine evaporitic and siliciclastic deposition influenced by re-activation of the transverse structures as contractional faults coevally to northwest motion of Iberia and transpressive stress regimen between Iberia and Sardinia (ca. 51 to 44 Ma). 2. Our study reveals how the tectono-sedimentary evolution of transverse faults in foreland basins record changes on motions and interactions of tectonic plates and continental blocks. These changes have an influence on structural evolution of transverse faults and the synchronous thickness and lithology distributions. Moreover, this work highlights the importance to analyse relationships between stratigraphic sequences, affected by transverse faults in orogenic chains, and tectonic processes as key to understanding kinematic histories of complex compressional and/or subduction zones. 3. The tectonic origin of transverse extensional faults in the South-eastern Pyrenees is similar to the Eastern Maghrebides where a lateral loading occurred presenting two orogenic salients. The South-eastern Pyrenees is an exceptional area to understand the structural evolution of transverse extensional faults active by bending foreland lithosphere. It provides field and subsurface evidences where thickness and lithology distributions are observed through the time as a response to the evolution of these structures. Moreover, the Late Paleocene to Middle Eocene stratigraphic record of this easternmost part of Iberia contributes to better understanding the complex geodynamic history of the Western Mediterranean region.

DATA AVAILABILITY STATEMENT
The data that support the findings of this study are available on request from the corresponding author. The data are not publicly available due to privacy or ethical restrictions.