Brigitte Van Vliet-Lanoë
email: brigitte.vanv-lietlanoe@univ-brest.fr

Christine Authemayou

Stéphane Molliex

Michael Hugh Field

Manfred Frechen

Pascal Le Roy

Julie Perrot

Valérie Andrieu-Ponel

Gwendoline Grégoire

Bernard Hallégouët

Middle Pleistocene seismically induced clay diapirism in an intraplate zone, western Brittany, France

Keywords: Clay diapirism, Liquefaction, Earthquake, Saprolite, Brittany, Variscan structures, Paleoseismology, Intraplate, Middle Pleistocene

The Brittany region of France is located in a low seismicity intraplate zone. Most of the instrumented earthquakes are limited to a shallow crustal depth without surface rupture. A paleoseismological analysis was performed on deposits on the Crozon Peninsula and in the Elorn estuary. We highlight hydroplastic deformations induced by liquefaction leading to clay diapirism, which were likely triggered by past earthquakes. This diapirism seems to be frequent in continental nonconsolidated sediments and to develop on the inherited tectonic structures, when a shallow water table and confining layers exist. Timing of deformation is dated using paleoenvironmental data, and electron spin resonance and infrared-stimulated luminescence dating methods. Two seismic periods were identified in western Europe during early Marine Oxygen Isotope Stage (MIS) 10 (~380 ka) and early MIS 8 (~280-265 ka). The lack of similar deformations affecting the Holocene tidal deposits in the Bay of Brest suggests that the magnitude of the triggering paleoearthquakes is probably higher (M w ~6) than the recent events (M w 5.4). These unusual intraplate major paleoearthquakes need specific factors affecting the far-field crustal stress loading to be triggered, such as a brief acceleration of the Africa-Eurasia lithospheric plate convergence, glacio-isostatic stress perturbations associated with the onset of major glaciations in northern Europe, or other processes induced by orbital forcing.

INTRODUCTION

Analyzing the paleoseismicity in a low-activity intraplate region is still challenging because most earthquakes are limited to a shallow crustal depth without evidence of rupture at the surface (Landgraf et al., 2016). In general, there is no long (>1 ka) paleoseismic record in stable continental regions. Western Europe is an intraplate region, influenced by both the collision between Eurasian and African plates (Ziegler and Dezès, 2006) and glacio-isostatic deformations [START_REF] Lambeck | Reply to comment by W. Fjeldskaar 'What about asthenosphere viscosity? Sea-level change, glacial rebound and mantle viscosity for northern Europe[END_REF][START_REF] Marotta | Combined effects of tectonics and glacial isostatic adjustment on intraplate deformation in central and northern Europe: applications to geodetic baseline analyses[END_REF]. Very few outcrops allow for the analysis of paleoseismicity in Quaternary deposits on continental western Europe. In eastern and central Brittany, Neogene and Quaternary deposits are preserved within small valleys or basins [START_REF] Vliet-Lanoë | Paleogeographic evolution of northwestern Europe during the Upper Cenozoic[END_REF]. In southern and western Brittany, such as the Crozon Peninsula and Bay of Brest, the low relief and "blocks and grabens" morphology inherited from the Neogene offer suitable conditions to preserve Quaternary deposits in depressions (Fig. 1). In this region, load casts and plastic deformations are commonly observed in Quaternary deposits and are interpreted as resulting from periglacial activity. Indeed, during the Quaternary glacial stages the region was occupied by discontinuous and/or continuous permafrost, allowing for the development of cryoturbation or the diagenesis of deformed sediments by differential frost heaving and solifluction processes [START_REF] Vliet-Lanoë | Distinguishing between tectonic and periglacial deformations of quaternary continental deposits in Europe[END_REF]; Table 1). However, frequent regional load cast figures do not correspond to periglacial criteria of formation, notably because of the absence of traces of ice segregation and of spatial regularity related to the classic control by desiccation polygons <1 m (Van [START_REF] Vliet-Lanoë | Distinguishing between tectonic and periglacial deformations of quaternary continental deposits in Europe[END_REF]. In Brittany, meter-scale load casts are present in the Aber Ildut valley to the west of Brest, in the Elorn River valley, and on the Crozon Peninsula (Fig. 1B). Seismic surveys and coring within the Bay of Brest did not permit us to observe deposits older than the Weichselian [START_REF] Grégoire | Control factors of Holocene sedimentary infilling in a semi-closed tidal estuarine-like system: the Bay of Brest (France)[END_REF]. However, the Trez Rouz and Pen Hat sections located on the Crozon Peninsula are exceptional sites exhibiting the most complete sections of middle and upper Quaternary units in Brittany. Studying these sites gives us an opportunity to investigate farther back in time (from ~400 ka to present). The Trez Rouz site is locally famous because of its wavelike peat deformations, traditionally attributed to periglacial activity [START_REF] Hallégouët | La presqu'île de Crozon: évolution géomorphologique[END_REF]. In addition, the entrance of the Elorn estuary (Fig. 1B) locally provided sediment traps, such as at the Kergleuz sections, which recorded and preserved an important part of the Plio-Quaternary geologic history. Sections presented in this study (Fig. 1B) were analyzed from stratigraphic and paleoseismic points of view before being compared to other dated sections with similar records: at Pénestin, close to the Vilaine River estuary (Southern Brittany), and at Le Gurp (Gironde River estuary, Aquitaine coast) (Figs. 1A and2).

In this article, we aim to define the tectonic, paleoenvironmental, and stratigraphic settings of several deformed Quaternary deposits to argue a mechanistic interpretation of the Quaternary deformation processes at the scale of western Brittany. We aim to determine whether nonperiglacial deformation might result from local or distal coseismic liquefaction of clayey saprolite, mud, and sand, because low magnitude earthquakes with no evidence of surface ruptures presently occur within the study area. Our data will provide a discussion of the processes responsible for the deformation structures, especially clay diapirism. We will also constrain the age of deformations with luminescence and paleoenvironmental methods. We rely on the age of these deformations with similar paleoseismic events recorded on (or close to) other Variscan structures in western Europe (Fig. 2). We also discuss these seismic events in the context of plate tectonics, glacio-isostatic mechanisms, and other less constrained processes, such as the impact of Earth's rotational velocity and mantle activity on lithosphere unloading.

MATERIAL AND METHODS

A high-resolution swath bathymetric map of the Bay of Brest obtained by a multibeam echo-sounder survey coupled to a detailed topographic map of the Crozon Peninsula and surrounding areas (Litto3D program) are used to locate the most prominent faults controlling the topography. Altitudes are given in meters NGF (French Ordnance Survey, 5 m NGF corresponding to above sea level [asl]). A seismic profile (Fig. 3) was also performed to distinguish the offshore continuity of onshore faults at depth. This profile was acquired with a Sparker Source, a six channel SIG streamer, and a numerical DELPH2 acquisition system. The sparker emission energy ranges from 250 to 500 J, the shooting rate is 1.5 s, and the vessel bottom speed is 5 knots. Seismic processing (simple stacking and filtering, Butterworth, 200-800 Hz) was done to enhance data quality with SPW (Parallel Geosciences Corporation) and Seismic Unix softwares. Paleoseismological analyses were performed from field observations, descriptions, and regular photographic surveys of coastal sections made from the 1980s onward, but mostly prior to 1997, the recent access of the sections being restricted because of their location in the Crozon Natural Park (dating on the pictures). With the present-day rising sea level, coastal environments are either refreshed or degraded, and several successive surveys allow us to constrain the observed deformed structures in three dimensions. For the Trez Rouz site, first palynological analyses were performed in 1974 by Figure 2. (color online) Map of the main superficial and subsurface faults and shearing zones compiled by the French Bureau de recherches géologiques et minières (1/1,000,000; https://www.geoportail.gouv.fr/donnees/cartes-geologiques) and the British Geological Survey. Details for Aquitaine: [START_REF] Rolet | The concealed basement of Aquitaine[END_REF] and [START_REF] Gély | L'évolution paléogéographique et géodynamique du Bassin aquitain au Paléogène: enregistrement et datation de la tectonique pyrénéenne[END_REF]; details for northwestern France and Belgium: [START_REF] Mansy | Role of basement 22 Brigitte van Vliet-Lanoë et al. reactivation in the formation and inversion of the Weald-Boulonnais basin[END_REF]. Official seismic risks map compiled by [START_REF] Giardini | Seismic Hazard Harmonization in Europe (SHARE): Online Data Resource[END_REF]European Seismic Hazard Map [ESHM]

4

Brigitte van Vliet-Lanoë et al.

M.T. Morzadec-Kerfourn (Mullenders extraction methodology).

A second analysis was made by V. Andrieu-Ponel at the Marseille laboratory using a treatment with a heavy liquid (d = 2) in complement to the sampling for seeds and plant remains [START_REF] Field | Significance of two new Pleistocene plant records from western Europe[END_REF]. All these data are summarized in Figure 4. Dating of the Quaternary units was previously performed by electron spin resonance (ESR) methods on the Trez Rouz, Pen Hat, and Penestin sections (M. Laurent, in [START_REF] Vliet-Lanoë | Neotectonic and seismic activity in the Armorican and Cornubian Massifs: regional stress field with glacio-isostatic influence[END_REF][START_REF] Beck | Late Quaternary lacustrine paleo-seismic archives in north-western Alps: examples of earthquake-origin assessment of sedimentary disturbances[END_REF]. New infrared-stimulated luminescence (IRSL) ages on Marine Oxygen Isotope Stage (MIS) 6 deposits were performed on the Pen Hat section (Table 2). Sand-sized K-feldspar minerals were prepared using the standard procedure including dry sieving, treatments with HCl and H 2 O 2 , and density separation (d < 2.58). The prepared K-feldspar samples with 150-200 μm grain size were mounted on stainless steel discs using silicone oil. Each aliquot contained 10-20 grains. All luminescence measurements were conducted with a Riso TL/OSL (Thermo-Luminescence/Optically-Stimulated Luminescence) DA-20 reader, and the luminescence signal was detected through a combination of blue filters (7-59 and BG-39) with a neutral density filter (ND10). The dose rates were calculated by radioactivity measurements of the samples using high-resolution gamma spectrometry. The equivalent dose (D e ) of the samples were measured using elevated temperature post-IR IRSL (pIRIR), stimulated at 290°C after preheating at 320°C for 60 s and with the first IR stimulation at 200°C. By increasing the first stimulation temperature from the standard 50°C [START_REF] Thiel | Luminescence dating of the Stratzing loess profile (Austria)testing the potential of an elevated temperature post-IR IRSL protocol[END_REF] to 200°C, the signal can be considered to be free from anomalous fading (e.g., [START_REF] Ito | Post-IR IRSL dating of K-feldspar from last interglacial marine terrace deposits on the Kamikita coastal plain, northeastern Japan[END_REF]. Ten aliquots each were measured for samples PHD1-6, and six aliquots were measured for GWD1-2 and TRG1-2.

Mean D e values with 1σ standard uncertainty were used to determine the ages of the samples, except for GWD1-2, whose D e values are in saturation. For these two samples, the minimum ages were calculated using 2 times characteristic saturation dose (D 0 ) obtained from the dose response curve. To test the validity of the applied pIRIR protocol, a dose recovery test was performed. Six aliquots each for PHD1-6 were bleached for 4 hours in a solar simulator (Hönle SOL2). Three aliquots per sample were given a laboratory beta dose of 250 Gy, and we tested if this laboratory dose could be accurately determined. The residual dose originated from the remaining pIRIR signal after the 4 hour bleaching was subtracted from the recovered dose. The dose recovery ratio has a tendency to overestimation of the given dose; however, they are all Middle Pleistocene seismically induced clay diapirism in an intraplate zone, western Brittany, France consistent with the acceptable range (0.9-1.1 Gy) within 2σ uncertainty.

GEOLOGIC SETTING

The Bay of Brest forms the boundary between two Variscan domains of the Armorican massif, the Leon domain to the north, and the central Armorican domain to the south (Fig. 1A; [START_REF] Ballèvre | Palaeozoic history of the Armorican Massif: models for the tectonic evolution of the suture zones[END_REF]. Crozon Peninsula is located to the south of the Bay of Brest and is affected by a low-grade brittle deformation with folding and faulting, along several main tectonic structures, mostly inherited from the Devonian Variscan orogeny and reactivated during the Eocene Pyrenean orogeny, the Oligocene rifting phase, and lately the Miocene-Alpine orogeny (Ziegler and Dézes, 2006;[START_REF] Anell | Cenozoic uplift and subsidence in the North Atlantic region: geological evidence revisited[END_REF]. The later tectonic phase mostly occurred from the Langhian-Serravalian to the late Messinian under a 150°t rending shortening axis [START_REF] Bergerat | Stress fields in the European platform at the time of Africa-Eurasia collision[END_REF] and produced vertical motions of west-east elongated blocks (Le [START_REF] Roy | Cenozoic tectonics of the Western Approaches Channel basins and its control of local drainage systems[END_REF]. Further vertical motion of the Crozon Peninsula did not occur during the Plio-Quaternary, as proved by an Upper Miocene low-stand platform preserved in several places close to the present-day sea level, at Porz Kubu (Fig. 1B) Le Rhys, close to Douarnenez [START_REF] Vliet-Lanoë | Paleogeographic evolution of northwestern Europe during the Upper Cenozoic[END_REF], and at Pénestin (Vilaine River estuary; [START_REF] Vliet-Lanoë | Seismically induced shale diapirism: the Mine d'Or section, Vilaine estuary, southern Brittany[END_REF] (Fig. 2).

The morphology of the Crozon Peninsula mainly results from Late Hercynian pedimentation (Triassic, ~250 Ma), further reshaped by successive marine transgressions, from the Early Jurassic (203 Ma) to Late Cretaceous (~65 Ma; e.g., [START_REF] Bessin | The Armorican Massif (western France): a two times exhumed relief shaped by planation surfaces in response to Iberia-Eurasia relative movements[END_REF]. The intense weathering of emerged Brittany, which occurred during most of the Cretaceous and Paleogene, provided up to 30 m of argillaceous saprolite [START_REF] Wyns | Evolution tectonique du bâti armoricain oriental au Cénozoïque d'après l'analyse des paléosurfaces continentales et des formations géologiques associées[END_REF] preserved in basement compartments downwarped by Cenozoic inversion tectonics. et al (2015). Those of Pénestin are extracted from [START_REF] Vliet-Lanoë | Neotectonic and seismic activity in the Armorican and Cornubian Massifs: regional stress field with glacio-isostatic influence[END_REF] and [START_REF] Vliet-Lanoë | Seismically induced shale diapirism: the Mine d'Or section, Vilaine estuary, southern Brittany[END_REF]. ESR, electron spin resonance; IRSL, infrared-stimulated luminescence; MIS, Marine Oxygen Isotope Stage; NGF, French Ordnance Survey (5 m NGF corresponding to above sea level); SASZ, South Amorican Shear Zone. The study zone is located north of the main active structure of the region, the South Amorican Shear Zone (SASZ; Fig. 1B). The two main brittle faults around the Crozon Peninsula are the right-lateral 145°trending Kerforne and the left-lateral 045-070°trending Elorn Fault Systems that branch eastward to the North Armorican Shear Zone [START_REF] Plusquellec | Notice explicative de la feuille Châteaulin, 1/50 000[END_REF][START_REF] Darboux | Notice explicative, Carte géol[END_REF][START_REF] Gall | LiDAR offshore structural mapping and U/Pb zircon/monazite dating of Variscan strain in the Leon metamorphic domain[END_REF]Figs. 1 and 2). The morphology of the Bay of Brest outlet (also called the Goulet) is largely controlled by these faults.

The opening of the Bay of Brest occurred after the lower Oligocene (lower Rupelian; [START_REF] Hallégouët | Découverte d'un dépôt oligocène inférieur dans la haute vallée de l'Aber Ildut, au nord-ouest de Brest (Finistère)[END_REF]. The collapse of regional deep weathered blocks of basement, well marked in the topography and visible in the digital elevation model (Fig. 1B andC) allowed for the further excavation of the Goulet during the Neogene and Quaternary transgressiveregressive repetitive marine stages, especially by tidal scouring at the entrance of the bay. The Goulet embayment flooded at least after the Middle Miocene as coastal Langhian faluns have been discovered on the peri-Armorican shelf off Ouessant Island and the Penmarch Peninsula [START_REF] Andreieff | Contribution à l'étude géologiques des approches de la Manche[END_REF], to the west of the Bay of Brest (Fig. 2). The oldest deposits in the Elorn estuary are Lower Pliocene (lower Zanclean, Za1, ~4.8 Ma) and suggest a limited opening of the bay (tidal clays) during this period. A landward transgression followed and allowed the deposition of a stratigraphic sequence constituted by poorly preserved Piacenzian units, Gelasian units dated at 2.4 ± 0.3 Ma by ESR methods, and finally the Tiglian muds [START_REF] Vliet-Lanoë | Paleogeographic evolution of northwestern Europe during the Upper Cenozoic[END_REF][START_REF] Darboux | Notice explicative, Carte géol[END_REF]. From the Gelasian onward, the sea level globally lowered. South of the Bay of Brest, the graben of Trez Rouz was flooded during most of the Pliocene [START_REF] Chauris | Notice de la carte géologique 1/50 000[END_REF]Fig. 1B).

SEISMIC SETTING

Currently, a number of Variscan structures are seismically active faults that could be assigned to the far stress field produced by the Africa-Eurasia convergence and the mid-Atlantic ridge push, with a northwest-southeast compressional stress and a northeast-southwest extensional stress deduced from numerical modeling analyses [START_REF] Golke | Origin of the European regional stress field[END_REF], and GPS and seismological studies [START_REF] Mazabraud | Short scale stress determination in central and western France, an intraplate slow deforming region[END_REF][START_REF] Perrot | Analysis of the Mw 4.3 Lorient earthquake sequence: a multidisciplinary approach to the geodynamics of the Armorican Massif, westernmost France[END_REF][START_REF] Nocquet | Present-day kinematics of the Mediterranean: a comprehensive overview of GPS results[END_REF]. The stress field is mainly strike slip with an extensional component attributed either to the buoyancy forces arising since the last deglaciation and explaining the general apparent uplift in Brittany [START_REF] Bonnet | Large-scale relief development related to Quaternary tectonic uplift of a Proterozoic-Paleozoic basement: the Armorican Massif, NW France[END_REF] or to local variations in rheology favoring stress perturbation [START_REF] Mazabraud | Is earthquake activity along the French Atlantic margin favoured by local rheological contrasts?[END_REF].

The seismicity around the Bay of Brest is significant, with several earthquakes [START_REF] Arroucau | Sismicité du Massif armoricain: relocalisation et interprétation tectonique[END_REF][START_REF] Haugmard | Détermination non-linéaire des paramètres hypocentraux et structuraux: application à la sismicité intracontinentale du Massif armoricain[END_REF] with low magnitudes between 1.6 ≤ M w ≤ 3.2 and within a range of depth from 2 to 25 km, following the SI-Hex catalog from 1965 to 2009 [START_REF] Cara | SI-Hex: a new catalogue of instrumental seismicity for metropolitan France[END_REF]. However, the seismicity in the western part of France shows an occurrence of moderate earthquakes with M w ≥ 4, such as the Hennebont earthquake in 2002, M w 4.3, along the SASZ [START_REF] Perrot | Analysis of the Mw 4.3 Lorient earthquake sequence: a multidisciplinary approach to the geodynamics of the Armorican Massif, westernmost France[END_REF]. During the instrumental period (since 1965), only one stronger earthquake occurred, with M w 5.0, close to La Rochelle [START_REF] Cara | SI-Hex: a new catalogue of instrumental seismicity for metropolitan France[END_REF] in 1972.

The historical seismicity presents a similar pattern along the Atlantic margin of France, with 61 events of 4 > M w > 5, and 7 with M w ≥ 5 between 1286 and 1959 (FCAT-17 catalog; [START_REF] Manchuel | The French seismic CATalogue (FCAT-17)[END_REF]SISFRANCE www.sisfrance.net). Four of them are located in the Vendée and La Rochelle area including one of the biggest historical events, the M w 6.3 Bouin earthquake that occurred in the Marais Breton in 1799 [START_REF] Manchuel | The French seismic CATalogue (FCAT-17)[END_REF]. Two other events occurred close to Meucon (in 1286 and1930), and one M w 5.4 event occurred near Quimper in 1959 probably along the SASZ [START_REF] Manchuel | The French seismic CATalogue (FCAT-17)[END_REF]. Our study sites are located 55 km from Quimper, 105-115 km from Hennebont, and at a minimum distance of 25-40 km from the SASZ.

DATA Trez Rouz

The key section, Trez Rouz, is located on the Roscanvel peninsula of Crozon (Figs. 1B and5) at the entrance of a southwest-northeast trending shallow depression linking the Bay of Brest to the Atlantic Ocean. The formations that outcrop along the coast (48°18.31′N; 4°34.28′W) are located at the western termination of a 060°trending graben, which is bound by satellite faults of the 060°trending Elorn Fault system. The site is also under the influence of several branches of the 150-160°trending Kerforne Fault System (Fig. 1B). Indeed, northwest-southeast trending faulting has been observed into Paleozoic shale that mostly crop out at the southern end of the beach. This graben preserves 20-25 m of saprolite and supposed Pliocene sandy gravel deposits [START_REF] Chauris | Notice de la carte géologique 1/50 000[END_REF]. The overall geometry of the graben infill seems rather similar to the observation extracted from a Seismic Sparker profile located east of the site. This profile shows the structures extending the graben onshore to the east in the Roscanvel Bay, partly affected by a flooded river incision, with possible Pliocene relic sedimentary bodies (Fig. 3).

The northwest side of the Trez Rouz graben (Figs. 3 and5) is shaped by at least two erosional surfaces. The highest one is perched at about 15 m NGF, to the north, and affects the faulted and weathered sandstones constituting the Pointe du Pouldu, located in the graben (Fig. 3). Following the clay mineralogy of the sediments [START_REF] Melou | Contribution à l'étude sédimentologique du Quaternaire Sud-Finistérien[END_REF], the Pointe du Pouldu corresponds to the feeding zone of most of the periglacial slope deposits in the main section. These deposits are rich in kaolinite and sandstone clasts (TRZ9). The southern part of the graben is occupied by ~2 m of Quaternary sediments lying on a much more shallow basement (5 m close to the rivulet, black arrow in Fig. 5A, emerging at ~50 m to the south of the river). The second erosional surface corresponds to a flat surface that truncates the shale and can be interpreted as a raised platform at 7 m NGF, related to the last raised beach (possibly Eemian, MIS 5; Fig. 5E). Periglacial deposits are commonly labeled "heads" in western France and the United Kingdom. Their description has been completed from [START_REF] Vliet-Lanoë | Neotectonic and seismic activity in the Armorican and Cornubian Massifs: regional stress field with glacio-isostatic influence[END_REF] for the TRZ1 to TRZ3 deposits (Fig. 5). New observations have shown a succession of four main upward dissymmetric injections or folds (A-B-C-D; Fig. 5), from north to south. The units are described subsequently.

TRZ1: Tidal mud, stratified and faulted. Only visible in the southernmost major injections (D) (Fig. 5D).

TRZ2: Lower peat complex. It consists of two levels of black peat separated by a gyttja layer and includes large fragmented wood remains.

TRZ3: Beige clayey loam, including angular clasts of quartz and unweathered shale. This unit is derived from the local saprolite [START_REF] Melou | Contribution à l'étude sédimentologique du Quaternaire Sud-Finistérien[END_REF] and is faulted by 060-070°t rending fractures spaced about 70 cm apart, with vertically injected clasts of quartz and shales concentrated in structure B (Fig. 5B).

TRZ4-6: Main peat complex (Fig. 5-2). Above this unit, most of the deposits are derived from the weathered Landevennec sandstone to the north [START_REF] Melou | Contribution à l'étude sédimentologique du Quaternaire Sud-Finistérien[END_REF]. It is the sequence studied by [START_REF] Morzadec-Kerfourn | Variations de la ligne de rivage armoricaine au Quaternaire[END_REF] and [START_REF] Field | Significance of two new Pleistocene plant records from western Europe[END_REF].

TRZ4a: Bleached clayey loam, faintly stratified, laterally linked to an organic gyttja type accumulation, truncated by a stone line, and faulted by 060-070°trending fractures spaced about 70 cm apart (Fig. 5E).

TRZ4b: Dark gray organic loamy clay. TRZ5a: Peat, homogenous, dark brown, evolving in a water-laid accumulation of wood and bark fragments (Pinophyta, Alnus, Betula, shrubs, ferns), often in a vertical position (injection).

TRZ5b: Humic sand-sandy schorre soil evolving laterally into a tidal channel.

TRZ6a: Peaty loam, rich in wood fragments and seeds (Potamogeton), richer in intraclasts to the top; partly fragmented to the top.

TRZ6b: Stratified alternation of silty sand and gravels, supported by a silty matrix, corresponding to a tidal channel to the south of the main structure B (Fig. 5B). It partly overlaps structure B and is slumped.

TRZ6c: Humic bleached loam reworking sand and quartz gravel and onlapping TRZ6a and TRZ6b. To the south, a thick peat complex is preserved at its base.

TRZ7-8: Lower head complex sloping to the north. TRZ7 consists of a series of two heads, including several weakly developed hydromorphic paleosoils, and ends with very coarse clasts of shale and quartz. At the top of structure B, the upper layer presents gravels in stacking (slump) with fragmented organic remains, resulting in the bulging of structure B. On the southern side of the main deformation, between the two most important injections, the TRZ7 layer is overthickened and almost horizontal at its base; it consists of a 1-m-thick deformed aggradation of sandy to gravely series including more loamy layers, somewhat organic, with rare wood fragments. It seems to correspond to a dell infilling. TRZ8 consists of a humic loam, interstratified sandy loamy heads deformed by periglacial ploughing blocks (Fig. 6B), with some coarse clean sand at the base, and it corresponds to a gyttja. A thick gley complex buried this layer.

TRZ9-10: This is a second sequence of heads, partly reworking dune sand and humic material, buried under bleached heads, including large clasts. TRZ10, a clayey loamy head, is rubified by pseudo-gley paleosol on reworked saprolite.

TRZ11-13: Upper sequence of yellowish heads, mostly loamy. It includes a brown podzolic paleosol at its base that evolves to the south of the rivulet into a raised cobble beach (TRZ12) including stony patterned ground to the north. These heads truncate the whole sequence with a slope toward the south.

Paleoenvironmental reconstruction of the main peat complex at Trez Rouz

The presence of Azolla filiculoides (Morzadec-Kerfourn, 1974) indicates a pre-Eemian age for the main peat complex (TRZ4-6). Peat layers are formed under mild temperate oceanic conditions with thermophilous aquatic plants, indicating deposition during interglacial periods. The absence of Tertiary pollen discounts an early Pleistocene age (last occurrence or last appearance datum/LADs) and is thus middle Pleistocene. This peat was first attributed by palynological studies to the Holsteinian (MIS 11; Morzadec-Kerfourn, 1974, 1999; [START_REF] Vliet-Lanoë | Neotectonic and seismic activity in the Armorican and Cornubian Massifs: regional stress field with glacio-isostatic influence[END_REF][START_REF] Field | Significance of two new Pleistocene plant records from western Europe[END_REF]. The sandy layer, TRZ5b, intercalated within the peat complex has been dated to 470 ± 53 ka by ESR methods (MIS 13;[START_REF] Vliet-Lanoë | Middle Pleistocene raised beach anomalies in the English Channel: regional and global stratigraphic implications[END_REF]. Several interglacial stages that had slightly higher sea levels than today (MIS 9, MIS 11, MIS 13, and MIS 15;[START_REF] Rohling | Antarctic temperature and global sea level closely coupled over the past five glacial cycles[END_REF] have similar palynological characteristics [START_REF] Scourse | Growth, dynamics and deglaciation of the last British-Irish ice sheet: the deep-sea ice-rafted detritus record[END_REF]. The plant macrofossil assemblage does not contain any biostratigraphic markers that allow for a specific interglacial stage determination.

TRZ4-5a corresponds to a marshy area partly submerged under a slightly brackish, shallow water body with both woodland (Quercus, Corylus, Fraxinus, Ilex aquifolium), waterside, and damp ground habitats (Alnus glutinosa, Osmunda regalis in TRZ5a) present on the surrounding slopes. After this humiferous clayey sand deposition, poor in plant macrofossils, the dated TRZ5b clayey sands were deposited within the marsh, probably associated with a slight cooling and storminess. This may be related to the opening of a coastal ridge or most likely to the colluvial infill by dune sand mixed with clayey slopewash, isolating the depression.

TRZ6a peat corresponds to a renewed shallow brackish lagoon with Caulinia goretskyi, Potamogeton occidentalis, and Ruppia maritima, possibly tracing another high stand in a forested landscape. The absence of both reed swamp and damp woodland in the site suggests a narrow depression with some spring seepage, with most of the local saprolite being eroded. The presence of Vitis and Buxus in TRZ6a suggests rather mild conditions, initially warmer than TRZ4-5, followed by the onset of limited climate cooling. These conditions can be correlated with the Cromerian Interglacial IV recorded as MIS 13 in the Netherlands [START_REF] Zagwijn | The Cromerian Complex Stage of the Netherlands and correlation with other areas in Europe[END_REF]. As the ESR dating gives a maximum age of 523 ka, the whole main complex might fit either with MIS 13a or MIS13c. It is thus covering the thermal spike of MIS 13a (~490 ka; [START_REF] Railsback | An optimized scheme of lettered marine isotope substages for the last 1.0 million years, and the climatostratigraphic nature of isotope stages and substages[END_REF]. The entire main peat complex (TRZ4-5a) should have formed during MIS 13a to the eustatic maximum MIS 13c. More precisely, TRZ4-5a might fit the MIS 13c (~525 ka) because palynological content indicates in this unit an ongoing transgression (Fig. 4).

TRZ6c could correspond to the classic Holsteinian (MIS 11; [START_REF] Koutsodendris | Vegetation dynamics and climate variability during the Holsteinian interglacial based on a pollen record from Dethlingen (northern Germany)[END_REF], which is longer and warmer than MIS 13. The unit is associated with TRZ6b, which represents a shallow tidal channel, probably initially 25 m wide, with a lateral peaty salt marsh (Fig. 7, phase I). Units TRZ1 and TRZ2 should correspond to interglacial high stands older than MIS 13, probably the complex MIS 15, with THREE high stands (MIS 15e, 15c, and 15a; Cromerian Interglacial III: [START_REF] Railsback | An optimized scheme of lettered marine isotope substages for the last 1.0 million years, and the climatostratigraphic nature of isotope stages and substages[END_REF]. They represent the closing of a tidal channel and two episodes of brackish peat formation (TRZ2) separated by a gyttja.

The last humic layer, TRZ10, consists of lacustrine organic sediment deposited near the shore. TRZ10 was considered as an upper interglacial stage by [START_REF] Morzadec-Kerfourn | Variations de la ligne de rivage armoricaine au Quaternaire[END_REF] presenting "cold-temperate" vegetation (Corylus, despite presence of Ilex) in the lower part and a late temperate period in the upper part. The unit could correspond to part of the MIS 7, as TRZ12 is clearly connected to the south with the last raised beach (MIS 5). The presence of a thin loess cover at the top of TRZ13 heads (Fig. 5E), with periglacial patterned ground, suggests that the whole sequence is older than Eemian (MIS 5). Below, the pseudo-gley soil (TRZ10) is rubified, partly inherited from the Landevennec sandstone, which is related to some seasonal hydromorphism, consistent with the morphology. This unit may have form during MIS 7 or earlier.

The remaining organic soil, the base of TRZ8, should, by deduction, correspond to a MIS 9 paleosol stretched by solifluction. Unit TRZ9 is deposited by periglacial solifluction and colluviation, which implies available moisture and seasonal frost. The conditions favorable to these processes normally occurred during early glacial conditions, when the sea level was not too low and the climate was still oceanic (with the presence of snow). This unit could be correlated with MIS 8c. The injection of structure A occurred early during this substage, because jarosite, a sulfate complex inherited from the bacterial degradation of pyrite issued from brackish or marine waters, is present in the structure.

Structures and deformation

Four main structures are considered here, from north to south A-B-C-D, forming upward dissymmetric hydroplastic injections, oriented 060°(Figs. 5, 6A, and8).

The northernmost structure A is the largest (5 m wide). This structure vertically deforms the unit TRZ8 and is present exclusively in saprolites. The B, C, and D structures are thin (<2 m wide). They are mostly tilted to the southwest, in conformity with the present-day slope. They rework several generations of peats (MIS 11: TRZ6; MIS 13: TRZ5; MIS 15: TRZ2). Deformation B is the most spectacular and affects the dated sandy layer (TRZ5b). Structure C reworks the two oldest peats and is cored by a third peat. Structure D only reworks the two older peat units (TRZ2) and is injected by stratified clays (Fig. 5D,TRZ1). Some small sandy liquefaction features occurred within fold B in association with the brecciation of the lower part (TRZ5a; Figs. 7B [red star] and 8-1). Similar brecciation deformed the upper part of the injected peat structure C. In TRZ4, quartz vein clasts are vertically oriented but in a steplike stacked succession as in an injection (Fig. 5A). Fragments of roots are mostly still in connected positions in the TRZ6a peat. None of the ligneous remains are intact in the TRZ5a and TRZ2 peats, and most of them are fractured wood and bark fragments and dispersed throughout the entire sediment. These suggest flotation and liquefaction. Many wood fragments were reworked with oriented clasts in unit TRZ6c, attesting to mudflow and erosion of the main peat injection. One of the deformation features occurred in the upper part of TRZ6c when the water table level was still high (floated fragments). Some larger deformations clearly suggest fold shortening between injection folds B and C (Fig. 5, unit TRZ7) and to the south of structure A.

A set of fissures, mostly parallel and vertical, affects most of the upward injected features (structures A, B, C, and D; Figs. 5 and8-2). The fissures trend 060°parallel to the length of the graben and to the direction of major faults in the Paleozoic basement that crops out on the northern side of the graben (Figs. 5-2 and 8-2 [diagram]). Horizontal microstructures that look like slickensides are present on some fissure planes. They are similar to slickensides on basement faults that trend 060°( Fig. 3). Such slickensides are common in clayey soils after desiccation [START_REF] Irfan | Mineralogy, fabric properties and classification of weathered granites in Hong Kong[END_REF], but here they are likely oriented following the local stress field. This deformation is contemporaneous with a general consolidation of the fissured units.

At a smaller scale, the injected masses in fold B are affected at the same time by antithetic microfracturation, dipping at 90°in the TRZ4 organic layer, with scalelike normal successive collapses and in a form of 20°dipping microshear planes in the stretched loamy sand (TRZ5b; Fig. 8-2, C andD). These deformation structures do not affect unit TRZ6a. They seem related to substratum retraction, fissuring, and consolidation by drainage and desiccation. These microfaults in the peat are underlined by jarosite. Structure D is located at the crossing of two injections trending 150°and 060°with lateral injection blisters (Fig. 6A). These two directions correspond to those of the Kerforne Fault System and the Elorn Fault System, respectively. These two fault systems join at the Trez Rouz site (Fig. 1B). Periglacial deformations end the record of each set of units including frost creep bending and ploughing blocks (Fig. 6B).

Sequence of deformation events

Retro-deformation analysis allowed us to restore stratigraphic units to their original geometries and to recognize several seismic events (Fig. 7, phase I). From MIS 15 to MIS12, a succession of fluvial and coastal units (TRZ1 to TRZ7) incised the thick saprolite preserved in the graben (phase I, Fig. 7). The first main deformation occurred from the end of TRZ6c deposition, probably at the time of the transition between MIS 11 and the onset of the MIS 10 cooling event, ca. 380 ka (phase II, Fig. 7). This deformation is associated with the formation of the B, C, and D structures, further bent by solifluction. After that, periglacial slope deposits partly filled the depression between the uplifted zone and the northern border fault (phase III, Fig. 7). The next seismic event is attributed to structure A, which is younger than most of the TRZ9 unit (phase IV, Fig. 7). This event produces only local retraction fissures, the other structures being (TRZ1 to TRZ6) now consolidated. An early MIS 8c age could be attributed to the event that deformed the TRZ9 unit. The deformation is possible because of the presence of marine water that impregnated the saprolite (jarosite), indicating that the sea level was still high despite the periglacial colluvial nature of TRZ9. This may also be the result of a persisting high sea level or saturation by seismic pumping from a shallow marine water wedge (Figs. 3 and7). In addition, periglacial solifluction was still active, regulating the northern side of the graben in two steps separated by paleosols (phases V and VI).

Pen Hat section

The Pen Hat section is located (48°16.6′N; 4°37.17′W) between the Toulinguet and Pen Hir Points (Figs. 4 and9A). The two ridges are armored by Armorican sandstones, raised by a northeast-southwest trending faulted anticline leading to the outcrop of Brioverian quartz-rich shales. Pen Hat anticline is separated from the Trez Rouz graben by the Kerforne Fault System (Fig. 1B). The Pen Hat section is located along a fossil cliff incised in the faulted quartz-rich shale at the contact with the sandstones (Fig. 9A). The bedrock was abraded by marine erosion and is partly buried by two series of Pleistocene deposits. A Holocene dune and a beach cobble ridge occupy the northern part of the bay. This section was completely independent of the Bay of Brest fluvial system during low stands [START_REF] Hallégouët | Formation de la rade de Brest[END_REF].

Stratigraphy (Fig. 4)

The most complete section is visible in the northeastern part of the Pen Hat bay, at the foot of Toulinguet Point, with an approximately 20-m-thick layer of Pleistocene sediment. The upper levels crop out along the lower to central portion of the bay. The lower contact is highly weathered and is thus the source of dark-gray clayey saprolite from the quartz-rich shales. The base of the Pen Hat section (unit PH4) yielded an ESR age of 403 ± 70 ka or MIS 11 [START_REF] Vliet-Lanoë | Middle Pleistocene raised beach anomalies in the English Channel: regional and global stratigraphic implications[END_REF].

The following units may be distinguished (Fig. 4): PH1: Marine cobbles and gravels in a silty (kaolinitic crystalline silt) sandy matrix, resting on a platform cut into quartz-rich shale, ~5 m NGF.

PH2: Mixed marine pebbles and clasts in a loamy matrix. This unit corresponds to a shore deposit with shore ice rafting characteristics.

PH3: Marine sands with some gravels. A lower paleolithic industry is included, close to the foot of the paleocliff. PH4: Hydromorphic "podzolic" sandy clayish paleosol, buried by gray clayey sands with some sandstone and shales flakes. This sand yields an ESR age of 403 ± 70 ka BP [START_REF] Vliet-Lanoë | Middle Pleistocene raised beach anomalies in the English Channel: regional and global stratigraphic implications[END_REF]. This unit is interpreted as an oxidized organic gley in a salt marsh depression located between two beach ridges.

PH5: Level of brownish sand deformed by load casts, consolidated.

Units PH4 and PH5 are locally pierced by hollows or pipes filled by marine sands with local concentrations of heavy minerals, with tubular stratified "marine" sand bodies infilling the hollow in PH4 (Fig. 9C). Pressure fracturing is marked at the base of the hollows. Some partly reworked archaeological artifacts are preserved (Fig. 4).

PH6: Dark-gray stratified clayish mud including lateral heads, thinning to the south, strongly consolidated.

PH7a: Stratified sand body interlayered with gray silty clay; it forms a storm built tidal bar, invading a low marsh (washover fan) and evolving into beach on the top, at the contact of the 7 m platform (southern part of the section). Behind this ridge, a lagoon facies with developed tidal clayey layers and with some admixture of clasts toward the top is present. The ESR dating (Van Vliet-Lanoë et al., 2000) and the IRSL dating failed because there was too little initial bleaching of the clayey sandstone saprolite, with aberrant age between 126 and 143 ka. However, the basal stratified marine sands from the similar and nearby Trégana section (washover fan; 6 km to the north, on granite; see Figs. 1B and4) provided IRSL ages of 254 ± 18 and 271 ± 20 ka (early MIS 7; see Table 2).

PH7b: Eolian sand including gravels and clasts in upper horizons, dropped from the overhanging cliff. The IRSL dating yields an age of 161 ± 12 ka.

PH8a: Eolian sand, slightly humic and some frost wedges in the southern part of the section. The IRSL dating yields an age of 138 ± 10 ka.

PH 8b: Humic podzolic soil developed in eolian sand. PH9: Eolian sand including some dropped clasts from the cliff.

PH10 to PH19 corresponds to the upper Pleistocene heads following the Eemian (MIS 5e) paleosol. PH11 yields an IRSL age of 105 ± 7 ka. They are connected with a cryoturbated beach in the central part of the bay. The Holocene dune sand, PH19, includes Pomatia elegans dated at 5220 ± 320 cal yr BP at its base.

Interpretation

The stratigraphic interpretation and correlations with the Trez Rouz units are given in Figures 4 and10. We argue for seismic deformation produced the load casts in unit PH5 (Fig. 8C) and the erratic hollows filled by tubular bodies, themselves filled by stratified marine sands, with local concentrations of heavy minerals in PH4. These structures represent a water escape, under pressure and confined upward by the overconsolidated stratified clay of unit PH6 (Fig. 9A, B, andD). The existence of overconsolidation of the lower units at Pen Hat, of a level of load casts in unit PH5, and of younger oriented water escape features in unit PH4-5, confined in space by unit PH6, argue for a series of earthquakes. PH4 yielded an ESR age of 403 ± 70 ka, suggesting the largest liquefaction events are younger, but older than PH7, a MIS 7a tidal complex (Fig. 4). The load cast event (PH5) is thus slightly younger than 403 ± 70 ka.

Kergleuz section

The Kergleuz section is located at the outlet of the Elorn valley, on an abrasion terrace on the right bank of the Elorn estuary (48°23.62′N; 4°24.29′W; Fig. 1B), at an elevation of 30 m NGF. The Elorn Fault System morphologically controls this part of the valley. Most of the sediments are clays and are preserved in a 12-m-deep, 150-m-wide incision cut within quartz-rich shale. The fluvial infilling (Fig. 3) mostly consists of Zanclean and early Quaternary deposits [START_REF] Darboux | Notice explicative, Carte géol[END_REF].

Description

A complex deformation system is evident from the stratigraphy proposed by [START_REF] Hallégouët | Formation de la rade de Brest[END_REF]. Deformations mostly affect the Piacenzian and Gelasian sands (ESR dating to 2.4 ± 0.35 Ma) and terrace gravel that rest on Zanclean clayey infratidal sediments [START_REF] Darboux | Notice explicative, Carte géol[END_REF]. These deformation structures are truncated by the periglacial slope deposits resting on the lower unit of the Middle Terrace group (unit Fxa of the French geologic map; [START_REF] Darboux | Notice explicative, Carte géol[END_REF], Middle Pleistocene seismically induced clay diapirism in an intraplate zone, western Brittany, France deformed like the Gelasian unit in the western wall of the Kergleuz section (Fig. 10).

Clay layers are folded mostly synclinally, before being truncated by the 30 m asl thin alluvial terrace Fxa. Fold limbs are limited by upward vertically injected sand masses and by three to four vertical large mud escapes of about 5 m wide (in yellow in Fig. 10B). These sand-escape diapirs are similar in shape and size to those of the clayish structure A at Trez Rouz. Shortening structures, dehydration fissures, and clay layer brecciation provide evidence for liquefaction and possible load casts (Fig. 10). The periglacial head cover is not deformed. Frost wedges are often superimposed on diapirs. The basal Zanclean tidal clays are only subject to wavy deformation. Frost wedges (Fig. 10B), probably sand wedges, postdate most of the deformations.

Interpretation

The setting up of the Fxa Terrace basal surface likely occurred in the late Miocene (western and northern France; 16

Brigitte van Vliet-Lanoë et al.

Van [START_REF] Vliet-Lanoë | Paleogeographic evolution of northwestern Europe during the Upper Cenozoic[END_REF], during the opening of the Bay of Brest. The age of this terrace is not accurately known but could be considered to be ~400 ka by comparison with similar ESR-dated Fxa deposits in northern and northwestern France [START_REF] Cordier | Fluvial system response to Middle and Upper Pleistocene climate change in the Meurthe and Moselle valleys (Eastern Paris Basin and Rhenish Massif)[END_REF][START_REF] Bahain | ESR chronology of the Somme River Terrace system and first human settlements in Northern France[END_REF][START_REF] Voinchet | ESR chronology of alluvial deposits and first human settlements of the Middle Loire Basin (Region Centre, France)[END_REF] (Fig. 4). The first major development of frost wedges in Brittany occurred in MIS 8 (280-240 ka; Pénestin section: [START_REF] Vliet-Lanoë | Neotectonic and seismic activity in the Armorican and Cornubian Massifs: regional stress field with glacio-isostatic influence[END_REF][START_REF] Beck | Late Quaternary lacustrine paleo-seismic archives in north-western Alps: examples of earthquake-origin assessment of sedimentary disturbances[END_REF]. The upper boundary of the Fxa terrace unit is deformed; its deformation might have occurred around the same time (same event?) as in Trez Rouz and Pen Hat (Fig. 3).

The Gurp section (Gironde estuary, 45.4345°N, 1.1530°W)

Background

Observations performed on this site (Figs. 1A and5) occur in the vicinity (<30 km) of the SASZ, in the entrance to the Gironde estuary [START_REF] Mazabraud | Short scale stress determination in central and western France, an intraplate slow deforming region[END_REF][START_REF] Mazabraud | Is earthquake activity along the French Atlantic margin favoured by local rheological contrasts?[END_REF]) (Fig. 2). Present-day seismicity in this zone is significant, with historical earthquakes up to M w 6 (SI-Hex: [START_REF] Cara | SI-Hex: a new catalogue of instrumental seismicity for metropolitan France[END_REF][START_REF] Manchuel | The French seismic CATalogue (FCAT-17)[END_REF]. The stratigraphy of the section was described by [START_REF] Tastet | Le Pléistocène de la façade atlantique du Nord-Médoc (France): état des connaissances sur la lithologie et la chronostratigraphie des «Argiles du Gurp» s[END_REF], which describes several units of clay, sands, and peats. The stratigraphic position of the Gurp clay topped with peat (called "lignite"; Fig. 4) was attributed to the ends of MIS 11 or MIS 9 [START_REF] Diot | Le Pléistocène de la façade atlantique du Nord-Médoc (France): synthèse sur la palynologie des «Argiles du Gurp» s.l[END_REF][START_REF] O'brien | Early and Middle Pleistocene vegetation history of the Medoc region, southwest France[END_REF]. Recently, [START_REF] Sitzia | The Quaternary coversands of southwest France[END_REF], using ESR methods, dated the Amélie sands located above the peat to 217 ± 31 and 222 ± 31 ka (Fig. 3), which correlated this unit to MIS 8 (243-280 ka).

Description

The Holstenian (MIS 9) littoral deposits are mechanically deformed by successive events of load casts. Here, deformations in the clay (~2 m thick) show sliding planes stained with iron hydroxides and load casts, associated with upward clay injection into the peat, leading to macroscopic circular patterns called pseudo-pingo features [START_REF] Dubreuilh | Notice explicative, carte géologique de la France (1/50000), feuille St-Vivien de Médoc -Soulac-sur-Mer (730)[END_REF][START_REF] Dubreuilh | Contribution à l'étude sédimentologique du système fluviatile Dordogne-Garonne dans la région bordelaise. Les ressources en matériaux alluvionnaires du département de la Gironde[END_REF]. These features are affected by a retraction fissure set and brecciation, which suggests a postdeformational dehydration. Beneath the peat, the lower part of the Amélie sands display ground patterned by small cryogenic features, disturbed by clustered liquefaction features (Fig. 11) comparable to load casts in seismic settings [START_REF] Vliet-Lanoë | Distinguishing between tectonic and periglacial deformations of quaternary continental deposits in Europe[END_REF].

Interpretation

The ages of the Gurp unit should be similar to that of the Trez Rouz site with stacked and deformed coastal interglacial deposits close to today's sea level, although deformations are commonly attributed to periglacial features [START_REF] Tastet | Le Pléistocène de la façade atlantique du Nord-Médoc (France): état des connaissances sur la lithologie et la chronostratigraphie des «Argiles du Gurp» s[END_REF][START_REF] Sitzia | The Quaternary coversands of southwest France[END_REF]Andrieux et al., 2017). Despite the interpretations of [START_REF] Sitzia | The Quaternary coversands of southwest France[END_REF] that correlated the Amélie sands with MIS 8, we chose to use the stratigraphic position provided by the ESR dating (217 ± 31 and 222 ± 31 ka) consistent with the MIS 7. The estuarine sand, beneath the Amélie sands, separated by a major discontinuity with a paleosol, and above the MIS 11 peat, was correlated with MIS 9. According to our observations, the Gurp section shows at least two main events: (1) at the end of MIS 11, after the peat formation and before the estuarine sands deposition; and

(2) at the end of MIS 9, older than the Amelie sands but younger than the estuarine sands, late MIS 7-based ESR dating. The lower part of the Amélie sands also shows evidence of smaller and younger events during the middle MIS 7.

DISCUSSION

Arguments for seismically induced deformations

Before this study, deformation structures at Trez Rouz, Pen Hat, and Kergleuz were considered to be a result of periglacial processes [START_REF] Melou | Contribution à l'étude sédimentologique du Quaternaire Sud-Finistérien[END_REF][START_REF] Morzadec-Kerfourn | Variations de la ligne de rivage armoricaine au Quaternaire[END_REF][START_REF] Hallégouët | La presqu'île de Crozon: évolution géomorphologique[END_REF], despite the lack of syngenetic evidence of segregated ice. However, the numerous features shown allow us to argue that deformations took place during or just after interglacial stages (MIS 11 and 9; Fig. 4). With this timing, the deformations should have been produced by other processesseismic events being the most probable (see Table 1). Neither fault reactivation nor diagnostic evidence of frost mounds has been found at Trez Rouz or Kergleuz despite evidence of late permafrost during MIS 8a and 2. The only periglacial features present are scarce frost fissures, ploughing block involutions, clast or cobble frost jacking, and bending of features by frost creep.

Most of the younger deformed units highlighted in this study indicate that they formed at a sea level close to the present day, thus in temperate to boreal conditions (MIS 11 and MIS 9a deposits; Fig. 4). Indeed, the MIS 9e, MIS 11a, and MIS 13a high stands have regionally reached +20 to +30 m (corresponding to the Middle Terrace Fx in estuaries), whereas the MIS 9a and MIS 5e are closer to today's sea level along the coast of western Europe [START_REF] Vliet-Lanoë | Middle Pleistocene raised beach anomalies in the English Channel: regional and global stratigraphic implications[END_REF]. The age of the youngest units affected, deduced by ESR dating and the stratigraphic evidence, puts them at MIS 11, 9a/8c. Stratigraphic evidence suggests that paleoearthquakes occurred during the onset of a cooling period with a limited glacio-eustatic regression (Fig. 4). This scenario postulates that the sea level dropped only by 10 to 20 m below its interglacial maximum of +20 and -30 m during MIS 11a and 9a, respectively, for most of the events considered in our study. At Trez Rouz, a sea level close to today's is likely for the youngest seismic event and is supported by the occurrence of jarosite along structure A. This implies the likelihood of a persisting brackish water prism persisting at these sites when the deformation occurred. Furthermore, deformations observed at Trez Rouz above the thickest body of saprolite (25 m) needed large water storage at depth to have occurred (Figs. 3 and4). Deformations are probably enhanced by recurrent liquefaction and a downslope mass movement, occurring when the water table was high. This situation occurs rarely during glacial stages. Deformations affecting the thickest alluvial layer (15 m) in the Kergleuz incision also suggest a high elevation of the water table. Erosional features visible at Trez Rouz, Pen Hat, and Kergleuz indicate overconsolidation processes and superficial fissuring of the deformed zones. Large water escape structures also exist at Kergleuz (Fig. 11), but their age is hard to determine more accurately than lower or middle Quaternary. The occurrence of subsurface overconsolidation of the saprolite or in the clayey units in the three sites, as well as intraclast formation (Kergleuz) or brecciation, or mud injection or water escape, suggests a rapid liquefaction and deformation in agreement with an assumption involving earthquakes [START_REF] Vliet-Lanoë | Seismically induced shale diapirism: the Mine d'Or section, Vilaine estuary, southern Brittany[END_REF]. The development of large deformations and shortening with injections may be related to a slower motion of the affected units over at least several days after the seismic initiation because it progressively affects the base of the overlying units.

Coseismic liquefaction and sediment deformation by clay diapirism

Detailed observations made on the three studied sites allow us to determine the likely mechanisms of deformation. Similar superficial liquefaction features associated with clay diapirism were already attributed to middle Quaternary seismic activity south of the study zone in the Vilaine estuary (Fig. 2; [START_REF] Vliet-Lanoë | Seismically induced shale diapirism: the Mine d'Or section, Vilaine estuary, southern Brittany[END_REF].

Liquefaction susceptibility reflects the relative resistance of soils (sandy or silty) to loss of fabric and strength when subjected to ground shaking. As water is a noncompressible fluid, the alternating squeezing and suction produced by vibrations from the local water table in the faulted Paleozoic basement induces an excess of pore water pressure, exceeding the liquid limit [START_REF] Woith | Heterogenous response of hydrogeological systems to the Izmit and Düzce (Turkey) earthquakes of 1999[END_REF]. Oversaturation is reached by fabric destruction (unconsolidated sand or porous saprolite) and/or by coseismic water pumping [START_REF] Muir-Wood | Hydrologic signatures of earthquake strain[END_REF][START_REF] Sibson | Geofluids: Origin, Migration and Evolution of Fluids in Sedimentary Basins[END_REF], which leads to loss of soil strength and liquefaction, or even fluidization [START_REF] Sibson | Geofluids: Origin, Migration and Evolution of Fluids in Sedimentary Basins[END_REF][START_REF] Montgomery | Streamflow and Water Well Responses to Earthquakes[END_REF][START_REF] Derakhshandi | The effect of plastic fines on the pore pressure generation characteristics of saturated sands[END_REF]. With fluidization and a very low dynamic viscosity, buoyancy forces and hydraulic fracturing develop in the confining layer, with differential loading [START_REF] Bolton | Fluid-flow pathways in actively deforming sediments: the role of pore fluid pressure and volume change[END_REF]. During the hydrofracturing process, an inner hydraulic pressure gradient was probably involved in tension gashes opening and forming "clastic" dykes in peat [START_REF] Sibson | Geofluids: Origin, Migration and Evolution of Fluids in Sedimentary Basins[END_REF][START_REF] Levi | Earthquake-induced clastic dikes detected by anisotropy of magnetic susceptibility[END_REF]. These processes also allowed (1) an incipient swell or doming that stretched the overburden cover [START_REF] Powley | Shale Domes. Search and Discovery Article #60001[END_REF], (2) mud injected close to its apex within the confining layer, and (3) an upward stretching and thinning of the flanking units (Figs. 5 and 8A;[START_REF] Vliet-Lanoë | Seismically induced shale diapirism: the Mine d'Or section, Vilaine estuary, southern Brittany[END_REF]. [START_REF] Youd | Lessons learned and need for instrumented liquefaction sites[END_REF] showed that ground oscillations during an earthquake led to a continual rise of pore water pressure after ground shaking ceases, triggering the prolonged postshock deformation. When the hydraulic pressure is lowered by dewatering, the clay diapirism normally stops, and the processes of consolidation begin, followed by a later collapse [START_REF] Powley | Shale Domes. Search and Discovery Article #60001[END_REF]. Coseismic pumping also leads to overconsolidation of sediments at depth by water expulsion to the surface, although the superficial layers are first liquefied and further overconsolidated after desiccation [START_REF] Youd | Lessons learned and need for instrumented liquefaction sites[END_REF]. This could explain the late retraction fissures formed by desiccation in the zones suspected to have been liquefied, especially at the Trez Rouz and Kergleuz sites.

Coseismic liquefaction is commonly studied in sands (e.g., [START_REF] Obermeier | Liquefaction Evidence for One or More Strong Holocene Earthquakes in Wabash Valley of Southern Indiana and Illinois, with a Preliminary Estimate of Magnitude[END_REF]. Because kaolinite-rich saprolites behave as quick clays along shore cliffs exposed to the swash impact of waves, it is better to use the term "clay diapirism" instead of "shale diapirism" in our study. This mechanism, already well observed at Pénestin [START_REF] Vliet-Lanoë | Seismically induced shale diapirism: the Mine d'Or section, Vilaine estuary, southern Brittany[END_REF], could explain most of the mesoscopic deformation structures observed in saprolites at the Kergleuz and Trez Rouz sites, which includes the extension in sands (unit TRZ5) and retraction fissures developed by dewatering and desiccation in the injected saprolite or liquefied sediments. The collapse of the interstitial porosity by vibration also explains why any remobilization by a latter event is almost impossible in sands and in saprolites. The mechanisms are not fundamentally different in both materials. Furthermore, clay diapirism also explains shortening fabric in partly consolidated material in the same places.

The stretching by solifluction and/or thinning of the overburden alluvial body, upslope at Trez Rouz, limited the confining pressure for the injection in structure A. Moreover, postdiapirism consolidation could have prevented a renewed liquefaction in the same location. This is also true for structures B to C, leading to the transfer of the liquefaction zone to the northern part of the graben (structure A). The confining of the saprolite under the accumulation of TRZ9 slope deposits at the northern edge of the graben, and the absence of preexisting overconsolidation, facilitate the upward migration of the saprolite in structure A (Fig. 7, phase IV).

Impact of altered Variscan structures on deformation initiation

The main deformation process highlighted at Trez Rouz, Kergleuz, and Pénestin (Van Vliet-Lanoë et al., 2009) is a small-scale clay diapirism, initially triggered by coseismic water pressure. Further loading and lateral spreading on an already deformed and deeply weathered shale bedrock allowed the clay diapirism to develop.

The initial bulge shape, induced by liquefaction, is modified by lateral or oblique upward migration of fluidized clays and may generate a network of clay dykes or sand pipes, as happened at Pen Hat. In classical clay diapirism, the diapirs usually develop in the form of ring or conic faults, which enhance the rise of the diapir. Here, and in all other studied sites, the bulge shape is parallel to the water escape from the basement, leading to injection swarms roughly parallel to the basement faults. This localized overpressure source results in tension gash development in the confining stratified alluvial system, accumulated in a trough to the south of structure B at Trez Rouz (Fig. 3). The folds and injections are parallel to the main basement structures, but also to the orientation of the preexisting confined channel sedimentary pattern (TRZ6b; Fig. 7, phase I). The parallelism between sand fissures and channels is similar to the sand dikes along the Mississippi River caused by the New Madrid 1811 and 1812 earthquakes [START_REF] Obermeier | Using liquefaction induced features for palaeoseismic analysis[END_REF]. The absence of surface fault rupture and the subparallel characteristic of the injections related to the basement faults result from a mechanism of water injection from basement faults into the capping saprolite (Fig. 3A). Mechanically, the water table level must not be deeper than 5 m from ground level to create these structures [START_REF] Olson | Revised magnitudebound relation for the Wabash Valley seismic zone of the central United States[END_REF].

Clay diapirism is rather common in Europe, for instance, in the Oligocene Bovey-Tracey basin located in Devon in the United Kingdom along the Sticklepath Fault (Fig. 2; [START_REF] Vliet-Lanoë | Seismically induced shale diapirism: the Mine d'Or section, Vilaine estuary, southern Brittany[END_REF], and in the lignite grabens in Germany and Poland [START_REF] Eissmann | Lagerungsstörungen im Lockergebirge: Endogene und exogene Tektonik im Lockergebirge des nördlichen Mitteleuropas[END_REF]. Site effects [START_REF] Semblat | Site effects: basin geometry vs soil layering[END_REF] on unconsolidated sedimentary bodies certainly amplified the likelihood of diapirism in all the studied sites near Variscan structures. The minimal thickness of saprolite needed to cause diapirism is about 10 m, as evident at Pénestin and Kergleuz.

Seismo-tectonic significance

Currently, the focal depth for most of major instrumental earthquakes indicates that the accommodation of the regional plate deformation is recorded mostly along major shear zones (mostly Variscan shear zones; [START_REF] Mazabraud | Is earthquake activity along the French Atlantic margin favoured by local rheological contrasts?[END_REF]. The present-day stress field is northwest-southeast with a clear strike-slip signature and a slight west-east transtensive component in Brittany [START_REF] Olaiz | European continuous active tectonic strainstress map[END_REF]. Seismic analysis produced after the Lorient/Hennebont earthquake in 2002 (M w 4.3 or Lg 5.4 magnitude scales) along the SASZ indicates a main 170°trending regional stress field [START_REF] Perrot | Analysis of the Mw 4.3 Lorient earthquake sequence: a multidisciplinary approach to the geodynamics of the Armorican Massif, westernmost France[END_REF]. No fault displacement affecting the surface has been detected from this event or during the historical period. Generally, the regional seismicity is rather low. No mud volcanoes have ever been observed on the Holocene mudflats of the Bay of Brest after the recent earthquakes (e.g., October 11, 2013;December 11, 2013;October 30, 2015;December 8, 2016;December 11, 2016;January 24, 2018: M l 3.4 to 3.9, local magnitude according to Comissariat à l'Energie Atomique network) or after the recent, major Hennebont (in 2002; M w 4.3) and Quimper (in 1959; M w 5.4) earthquakes. Nevertheless, Holocene submarine slumps in estuarial clays are present at the entrance of the Elorn estuary [START_REF] Grégoire | Control factors of Holocene sedimentary infilling in a semi-closed tidal estuarine-like system: the Bay of Brest (France)[END_REF]Fig. 1B: S). Clay diapirism associated with liquefaction has only been observed in middle Pleistocene deposits in this region. This process needs specific conditions that do not currently occur in Brittany and surrounding areas, despite the fact that environmental conditions look similar between the middle Pleistocene and today, notably from sedimentary and sea level points of view (relative sea level close to the 0 NGF; see the "Trez Rouz" section). Therefore, as for middle Pleistocene peat deposits, Holocene mudflats in the Bay of Brest could easily be marked by mud ejection above basement faults in case of a major earthquake.

The only major regional faults able to generate strong earthquakes are the Elorn and the Kerforne Fault Systems and the SASZ (Figs. 1 and2). The major system, the SASZ, is the farthest away, >25 km, from our studied sites, while the other faults exist at depth beneath our studied sites (Fig. 1). To estimate the range of magnitudes possible to explain the paleoearthquakes, we considered regional magnitude-bound relations adapted from stable continental regions, based on relationships between earthquake magnitude and the most distal occurrences of level-ground liquefaction from the earthquake source [START_REF] Obermeier | Liquefaction Evidence for One or More Strong Holocene Earthquakes in Wabash Valley of Southern Indiana and Illinois, with a Preliminary Estimate of Magnitude[END_REF][START_REF] Olson | Revised magnitudebound relation for the Wabash Valley seismic zone of the central United States[END_REF], because no magnitude-bound relationships are available for Brittany. Based on the minimum distance between the studied sites and the SASZ (25-40 km), a minimum M s of 6.4 to 6.9 is required to produce liquefaction at the studied sites (M s values are similar to M w when they are greater than 5.5; [START_REF] Lolli | Empirical conversion between teleseismic magnitudes (mb and Ms) and moment magnitude (Mw) at the Global, Euro-Mediterranean and Italian scale[END_REF]. For the other two faults, located closer, a M s of 5.5 is enough. If we consider that a saprolite has high liquefaction susceptibility, because of its fragile porous fabric, the required minimum M s for an earthquake along the SASZ can be reduced to an M s of 6 to 6.4 and reduced to M s 5 for the closest faults (Kerforne and Elorn Faults).

Moreover, because earthquakes occurring during periods of relatively high sea levels and consequently high water tables, with site amplification the magnitude need not be so high for significant deformation. Indeed, the discharge created by a lowering sea level could open the fissures, promoting deeper water infiltration, reducing the friction and thus limiting the accumulation of strain, especially on shale or clay. Fluid pressure would propagate as slow stress waves along fissure paths and could enhance the effects of earthquakes [START_REF] Sibson | Geofluids: Origin, Migration and Evolution of Fluids in Sedimentary Basins[END_REF][START_REF] Segall | Dilatancy, compaction, and slip instability of a fluid-infiltrated fault[END_REF][START_REF] Mulargia | Anthropogenic triggering of large earthquakes[END_REF].

Geodynamic significance

Any M w 5 earthquake along the Kerforne and Elorn Faults, or M w 6 earthquake along the SASZ, is documented by historical and instrumental seismicity, but the possibility that such events occurred during the Quaternary cannot be excluded. Indeed, according to the length-magnitude scaling relationship [START_REF] Wells | New empirical relationships among magnitude, rupture length, rupture width, rupture area, and surface displacement[END_REF], the ~40-km-long Elorn and Kerforne faults and the ~100-km-long segments of the SASZ possibly could generate earthquakes of a maximum magnitude ~6.8 and ~7.4, respectively. Recently, Calais et al. (2016) suggested that large earthquakes in intraplate contexts such as in Brittany do not obey the seismic recurrence cycle, such as those located at the frontiers of plates. Moreover, the authors suggest that large earthquakes in stable continental regions are better explained by transient perturbations of local stress on fault strength that release elastic energy from a prestressed lithosphere.

From the stratigraphic analysis and the available dating, at least two periods of clusters of earthquakes are evident in western France: (1) at the onset of MIS 10, ~380 ka, and (2) at the onset of MIS8, ~280-265 ka (Fig. 3). No older or younger important deformations is evident, except some younger events in the Aquitaine during the Younger Dryas stade [START_REF] Schneider | Déformations co-sismiques à Cestas-Pot-aux-Pins et Larrusey[END_REF]. The second period (280-265 ka) is frequently associated with other potential paleoseismic events in western Europe. It could be the case for major falling down of the Sangatte paleocliff, the Cagny Epinette living-floor faulting in the Somme valley, and recurrent load casts in saprolite in central Brittany (Fig. 2; [START_REF] Vliet-Lanoë | Paleogeographic evolution of northwestern Europe during the Upper Cenozoic[END_REF]. Nevertheless, some deformations such as reverse faulting affecting the 160 ka loesses of Sangatte are now interpreted as thermokarstic features [START_REF] Vliet-Lanoë | Quaternary thermokarst and thermal erosion features in northern France: origin and palaeoenvironments[END_REF]. Sediment records in our study zone are incomplete because of lack of exposures [START_REF] Mccalpin | Paleoseismology[END_REF][START_REF] Caputo | Ground effects of large morphogenic earthquakes[END_REF][START_REF] Beck | Late Quaternary lacustrine paleo-seismic archives in north-western Alps: examples of earthquake-origin assessment of sedimentary disturbances[END_REF], so other earthquakes may have been missed. We highlight many discontinuities in the records, even when overconsolidation normally limits the erosion efficiency (see Fig. 3). Dating and paleoclimate records clearly indicate that the seismic events occurred during the onset of MIS 10 and 8, so we postulate that these two periods of seismic events could be triggered by of the rapid growth of major ice sheets (particularly MIS 10; [START_REF] Vliet-Lanoë | Neotectonic and seismic activity in the Armorican and Cornubian Massifs: regional stress field with glacio-isostatic influence[END_REF][START_REF] Vliet-Lanoë | Paleogeographic evolution of northwestern Europe during the Upper Cenozoic[END_REF]. We did not highlight any paleoseismic events at the onset of the glaciations that are more recent than MIS 8 in western France. The Gurp section, located closer to the Alpine sutures, is still active. However, the ice caps of the less extensive subsequent glaciations (MIS 6, 4, 2) could be too small and too far away from the study zone to produce a forebulge that adequately deforms the Armorican massif lithosphere and adequately increases the stress loading. Also, some cryogenic consolidation of the saprolite (permafrost) and old tidal muds during the cold Saalian I (MIS 8; 280-240 ka) could prevent new liquefaction during the onset of more recent glaciations.

Glacio-isostatic stress perturbation is not the only factor that could trigger a major earthquake in Brittany. Internal Earth dynamics such as, for example, changes in pressuredriven flow within the asthenosphere or slab sinking could also bring answers to the origin of these middle Pleistocene intraplate seismic events. Indeed, these paleoearthquakes could have been induced by the African-Eurasian plate convergence, and, more precisely, a brief acceleration of plate motion during the investigated period is possible [START_REF] Iaffaldano | Rapid plate motion variations through geological time: observations serving geodynamic interpretation[END_REF]. This convergence is also postulated as the 20

Brigitte van Vliet-Lanoë et al.

origin of other paleoseismic events during the same period, occurring to the south of the studied zone, closer to the plate boundary. For instance, on the northern foot slope of the Pyrenean mountains, superficial ruptures have been evidenced, deforming the upper terrace of the Adour River at Meilhan dated to <500 ka (Fig. 2; [START_REF] Carbon | Rupture de surface post-Pléistocène moyen dans le Bassin aquitain[END_REF]. Moreover, to the south of the Western Alps, coseismic ruptures dated to between 100 and 300 ka occurred along the Trévaresse Fault (Fig. 2; Aix-en-Provence: [START_REF] Chardon | First paleoseismological constraints on the strongest earthquake in France (Provence) in the twentieth century[END_REF]. Additionally, orbital forcing enhanced by solar activity and affecting Earth's rotation velocity could be envisaged as another source of lithospheric loading favoring seismic periods [START_REF] Bendick | Do weak global stresses synchronize earthquakes[END_REF]. This forcing also influences the climate and dynamic processes in the liquid metallic core of Earth controlling the magnetic field [START_REF] Hide | Earth's variable rotation[END_REF][START_REF] Worm | A link between geomagnetic reversals and events and glaciations[END_REF][START_REF] Encrenaz | The Solar System[END_REF][START_REF] Laut | Solar activity and terrestrial climate: an analysis of some purported correlations[END_REF]. Perturbations in the deepest part of Earth could also affect the stress field near the surface by possibly accelerating mantle convection and plume activity in the asthenosphere [START_REF] Larson | Mantle plumes control magnetic reversal frequency[END_REF][START_REF] Worm | A link between geomagnetic reversals and events and glaciations[END_REF]. Curiously, the highlighted seismic periods occurred during the onset of glaciations and at the same time as two magnetic low (geomagnetic low dipole) periods at 380 and 280 ka [START_REF] Worm | A link between geomagnetic reversals and events and glaciations[END_REF][START_REF] Ménabréaz | The geomagnetic dipole moment variation between 250 and 800 ka BP reconstructed from the authigenic 10 Be/ 9 Be signature in West Equatorial Pacific sediments[END_REF]

CONCLUSION

The combined examination of dating, geomorphology, and sedimentologic and paleoecological data from sites with Quaternary deformations in western Brittany allow for an analysis of the conditions that drive paleoseismicity in an intraplate region of Europe. From a mechanical point of view, clay diapirism development and liquefaction in regions with limited tectonic activity is linked to the occurrence of an unconsolidated sedimentary or saprolitic clay unit, with the presence of a rather shallow water table and a confining surface layer. The triggering seismicity is not necessarily important (not much higher than present day) and is probably associated with site effects. From a geodynamic point of view, two periods of coseismic events are evident: during early MIS 10 (~380 ka) and early MIS 8 (~280-265 ka).

They may be associated with the dynamics of Africa-Eurasia's plate convergence, in part accommodated in the far field by reactivation of major shear zones (Fig. 2). Still open to discussion is the possible glacio-isostatic impact on the stress field during the onset and decay of ice sheets together with their forebulge dynamics and the influences of orbital forcing on lithospheric loading.

Figure 1 .

 1 Figure 1. Location of the sections on the Crozon Peninsula and on the Elorn estuary. (A) Location of studied zones in the Armorican Massif framework. (B, C) LIDAR (light detection and ranging) image of the Bay of Brest entrance from the Litto3D LIDAR pictures. Completed for tectonic with data from the geologic map of the Department of Finistère, Bureau de recherches géologiques et minières 2008. P, Pénestin; R, sparker seismic profile Roscanvel (see Fig. 8); S, submarine slump. In panel C, major faults are in red; secondary faults are in black. (For interpretation of the references to color in this figure legend, the reader is referred to the web version of this article.)

Figure 3 .

 3 Figure 3. (color online) Trez Rouz graben geologic section with saprolite depth at scale (A) and a comparable Sparker seismic profile no. P4450, on its prolongation in the Bay of Roscanvel (B). Stippled lines below structures A to D in the upper diagram correspond to probable water-supplying faults preserved into the saprolite.

Figure 4 .

 4 Figure 4. (color online) Relative chronology for the study sites, with the stratigraphic position of the coseismic and periglacial deformations. Ages between brackets are estimated. The stratigraphy for the Gurp section follows O'Brien and Jones (2003) and Sitzia et al (2015). Those of Pénestin are extracted from Van Vliet-Lanoë et al. (1997) and Van Vliet-Lanoë et al. (2009). ESR, electron spin resonance; IRSL, infrared-stimulated luminescence; MIS, Marine Oxygen Isotope Stage; NGF, French Ordnance Survey (5 m NGF corresponding to above sea level); SASZ, South Amorican Shear Zone.

Figure 5 .

 5 Figure 5. Trez Rouz sites. (1) Location of the Trez Rouz section (in red) on the high-resolution digital elevation model. (2) Log of the section and detail of the section structures A, B, C, D, and E. Black arrow: rivulet. Note the shortening deformation figures of the units TRZ9 and TRZ7. The red square corresponds to the section investigated for plant macrofossil assemblage and palynology. (For interpretation of the references to color in this figure legend, the reader is referred to the web version of this article.)

Figure 6 .

 6 Figure 6. (color online) (A) Lateral injection in structure D, deforming unit TRZ2 (observation 2017). (B) Ploughing (soliflucted) block involutions between structures B and C, unit TRZ10; deforming TRZ8 organic soil (2004).

Figure 7 .

 7 Figure 7. Reconstitution of different phases of deformation of the Trez Rouz section. Red stars stand for injection. MIS-Marine Oxygen Isotope Stage; NGF-French Ordnance Survey (5 m NGF corresponding to above sea level). (For interpretation of the references to color in this figure legend, the reader is referred to the web version of this article.)

Figure 8 .

 8 Figure 8. Details of sites 8-1 and 8-2. Structure B details (Fig. 5). (A and B) Site 8-1 showing main injection, liquefaction (red stars in Fig. 7) retimed by stretching fissuring (slickensides in white), C-D to the north and B to the south; notice dehydration fissures below the injection in A. Site 8-2 showing stretching (stretching fissure) (C) and fissuration (retraction) (D) related to the dehydration of the liquefied saprolite. (E) Simple fissuring by dehydration, north of structure C. Notice the antithetical character of the stretching fissures on each side of the injection: C-D to the north and B to the south. The diagram plots the orientation of the basement faults and the two types of fissuration in the sediment. (F) Fault planes and associated slickenside (black arrows). (For interpretation of the references to color in this figure legend, the reader is referred to the web version of this article.)

Figure 10 .

 10 Figure 10. Deformation at Kergleuz during the 1995 excavation of the highway across a paleochannel. (A) Location of the section on a digital elevation model: sections in red. (B) Sections. yellow: main water-escape features. Blue: ca. 430 ka Fxa gravel terrace. Orange: Gelasian (~2.4 Ma). Please note the different vertical and horizontal scales. (For interpretation of the references to color in this figure legend, the reader is referred to the web version of this article.)

Figure 11 .

 11 Figure 11. (color online) Le Gurp study site. (A) Regional digital elevation model. (B) Coseismic involutions above the Le Gurp peat (courtesy J.L. Schneider). (Cp) Small, periglacial ground pattern in the Amélie Sands. (Cs) The same unit deformed by liquefaction.

Table 1 .

 1 Characteristics of cryoturbations and coseismic involutions (completed from[START_REF] Vliet-Lanoë | Distinguishing between tectonic and periglacial deformations of quaternary continental deposits in Europe[END_REF].

	Deformation	Process	Microfabric	Local conditions	Organization	Orientation
	Cryoturbation	Differential frost (ice	Platy to blocky	Depression, downslope,	Random adjacent	None or slope
	(periglacial)	segregation)	Silty caps,	moist to wet	or regular	direction
		heave, positive or negative	particle	Superficial liquefaction	(thermal or
		gradient of frost	translocation	(<20 cm)	desiccation net)
		susceptibility	Surface crusting		Decimeters,
		Small diapirs: injection	(pseudo-		regular in size
		(negative gradient of frost	liquefaction)
		susceptibility)
		Stone jacking
	Solifluction	Differential frost	Platy to granular,	Slope form 2°to steeper	Global from a	Globally
	(periglacial)	heave, positive or negative	Silty caps,	Massive form on drained	single surface or	slopeward
		gradient	aggregate	slope	in superimposed	Common except
		lateral and downward	rotation, particle	Stretched form on poorly	tongues	on south facing
		sliding during thaw	translocation	drained slope (snow		slopes
		consolidation and drainage	Surface crusting	patch)
		Stone jacking downstream	(pseudo-	Possible liquefaction at the
			liquefaction)	lower end
	Load cast	Stratified fast sedimentation	None specific,	Water saturated or	Random or local	None or surface
	(sedimentary)	with some contrast in	layer rupture	underwater	millimeters to	stretching in
		viscosity and density			decimeters	slope direction
	Load cast	Impact, stratified	None specific,	Water saturated or shallow	Local or along a	None
	(single	sedimentation with some	rather massive	underwater	shore
	shock)	contrast in viscosity and	layer rupture,		millimeters to
		density	locally kink		meters
		(e.g., swach, trampling)
	Load cast	Vibration at the surface of	None specific,	Water saturated or under	Patchy, in nodal	None to very well,
	(earthquake)	stratified sedimentation	rather massive	shallow water (<1 m) or	position,	unidirectional,
		with some contrast in	layer rupture	above a shallow water	In association	sinking in slope
		viscosity and density	wrinkling,	table (<5 m deep)	with water	direction
			rock clast		escape
			rotation,		millimeters to
			microinjections,		decimeters,
			overconsolidation		irregular in size
	Large	One or all the following	Stretching,	Depression, basin	Control by seismic
	involutions	processes:	overconsolidation	Wet to shallow water;	wave length or
	Folds trains,	earthquake	Secondary loads,	in series,	basement faults
	Diapirs	compression/tension	water injection,	above a shallow water	(or Riedel)
		slump (e.g., valley bulging)	shear, microfaults	table (<5 m deep)	meters to
					decimeters

Oriented

Middle Pleistocene seismically induced clay diapirism in an intraplate zone, western Brittany, France

Table 2 .

 2 Infrared-stimulated luminescence data and results for the Pen Hat section. "NGF" indicates French Ordnance Survey (5 m NGF corresponding to above sea level). "Lum Nr" indicates Luminescence Number.

													Error							Number
												Dose	in				Age	Dose	Error in	of
	Geologic	Sedimentology	Sample	Lum	U	Error	Th	Error	K	Error	Water	rate	dose	D e	Error	Age	uncertainty	recovery	dose	aliquots
	unit	remarks	code	Nr	(ppm)	U	(ppm)	Th	(%)	K	(%)	(Gy/ka)	rate	(Gy)	in D e	(ka)	(ka)	ratio	recovery	measured
	PH11	Weathered dune	PHD1	3728 1.54 0.08	4.36	0.22 1.54 0.08	6	2.67	0.15 280	8 105	7	1.18	0.04	10
		sand,
		with some gruss
	PH8a	Coarse grain dune	PHD2	3729 1.43 0.07	3.73	0.19 1.43 0.07	4	2.53	0.15 350 16 138	10	1.19	0.07	10
		sand,
		very slightly
		organic,
		frost wedge
	PH7b base Stratified sands,	PHD3	3730 1.27 0.06	3.43	0.17 1.27 0.06	5	2.29	0.14 369 15 161	12	1.11	0.05	10
		probable storm
		dune
	PH7a top Stratified sands,	PHD4	3731 1.36 0.07	3.17	0.16 1.36 0.07	3	2.43	0.15 306	6 126	8	1.12	0.08	10
		with gray clayish
		colluvial beds
		(washover in a
		marsh,
		summit +10 m
		NGF)
		insufficient
		bleaching
	PH7a		PHD5	3732 1.41 0.07	4.49	0.22 1.41 0.07	8	2.47	0.14 352 15 143	10	1.05	0.04	10
	middle
	PH7a base		PHD6	3733 1.46 0.07	6.01	0.30 1.46 0.07	9	2.59	0.14 370	7 143	8	1.19	0.09	10
	TR2	Washover tidal	TRG1	3737 1.51 0.08	4.29	0.21 1.51 0.08	4	2.66	0.16 678 24 254	18			6
	summit	marsh infilling,
		culminating
		at 9 m NGF
	TR2 base		TRG2	3738 1.51 0.08	2.83	0.14 1.51 0.08	4	2.59	0.16 702 30 271	20			6

https://www.cambridge.org/core/terms. https://doi.org/10.1017/qua.2018.63 Downloaded from https://www.cambridge.org/core. Worldwide Information Pty Ltd, on 31 Aug 2018 at 23:24:31, subject to the Cambridge Core terms of use, available at

https://www.cambridge.org/core/terms. https://doi.org/10.1017/qua.2018.63 Downloaded from https://www.cambridge.org/core. Worldwide Information Pty Ltd, on 31 Aug 2018 at 23:24:31, subject to the Cambridge Core terms of use, available at

Brigitte van Vliet-Lanoë et al. https://www.cambridge.org/core/terms. https://doi.org/10.1017/qua.2018.63 Downloaded from https://www.cambridge.org/core. Worldwide Information Pty Ltd, on 31 Aug 2018 at 23:24:31, subject to the Cambridge Core terms of use, available at

Brigitte van Vliet-Lanoë et al. https://www.cambridge.org/core/terms. https://doi.org/10.1017/qua.2018.63 Downloaded from https://www.cambridge.org/core. Worldwide Information Pty Ltd, on 31 Aug 2018 at 23:24:31, subject to the Cambridge Core terms of use, available at

ACKNOWLEDGMENTS

This research did not receive any specific grant from funding agencies in the public, commercial, or not-for-profit sectors. We first thank the "Reserve Naturelle Régionale de Crozon" in Brittany to have allowed the complementary sampling for new IRSL dating at Pen Hat. We thank especially Jamie Shulmeister and Roy Van Arsdale for their constructive recommendations to improve our paper, including enhancements to scientific clearness and language. We thank also an anonymous French reviewer for his judicious comments on the French seismic activity.


Alain Lescoat
email: alain.lescoat@chu-rennes.fr

MD Cristian Caimmi

MD Éric Hachulla

Patrick Jego

PhD Guillaume Coiffier

PhD Patrick Jégo

Imaging in Systemic Sclerosis make "US" great again?

Keywords: Systemic sclerosis, Scleroderma, Digital Ulcers, Power Doppler Ultrasonography, Ultrasound, Skin. Dear Editor

teaching and research institutions in France or abroad, or from public or private research centers.

Accepted Article

In this issue of Arthritis Care & Research, T. Santiago and colleagues [START_REF] Santiago | Ultrasonography for the assessment of skin in systemic sclerosis: a systematic review[END_REF] have conducted an informative systematic literature review on the relevance of ultrasound (US) for the assessment of skin involvement in systemic sclerosis (SSc). In this Original Article, the authors highlight the heterogeneous results concerning the reliability of US in studies published to date, but also underscore very promising data from recently published articles [START_REF] Sulli | Subclinical dermal involvement is detectable by high frequency ultrasound even in patients with limited cutaneous systemic sclerosis[END_REF]. Moreover, this work raises the issue considering US evaluation as an interesting alternative to mRSS for skin assessment in SSc. Therefore, beyond imaging considerations, this article questions the place of US evaluation in the global management of SSc. This topic is still a controversial matter as new assessment technics are emerging for microvascular, macrovascular and musculoskeletal evaluations in SSc [START_REF] Cutolo | Imaging of connective tissue diseases: Beyond visceral organ imaging?[END_REF][START_REF] Cutolo | Macrocirculation versus microcirculation and digital ulcers in systemic sclerosis patients: Macro-microcirculation and scleroderma[END_REF].

Nonetheless, a central question remains concerning these new evaluation tools for this disease: new imaging technologies may come, indeed, but for what purpose in the end?

In addressing this issue from a purely descriptive perspective, as the question of reframing and reshaping phenotypes in SSc is arising [START_REF] Ligon | Editorial: Scleroderma: Bringing a Disease From Black-and-White Into Technicolor[END_REF], US evaluation may offer a more precise description of vascular and musculoskeletal manifestations of the disease. This is all the more relevant as always more performant probes (>18Mhz) lead to always more precise descriptions and high-quality images. US characterization of these manifestations, combined with a precise clinical phenotyping [START_REF] Allanore | Review: Defining a Unified Vascular Phenotype in Systemic Sclerosis[END_REF] may help to enlarge the current view of SSc bringing the disease into the "technicolor" enriched vision of SSc suggested by Ligon and Wigley [START_REF] Ligon | Editorial: Scleroderma: Bringing a Disease From Black-and-White Into Technicolor[END_REF]. Beyond echocardiography and the detection of Pulmonary Arterial Hypertension, US may also have a part to play in the early detection of visceral manifestations of the disease such as lung fibrosis, as suggested by recent studies evaluating the diagnostic performances of US in this field [START_REF] Tardella | Ultrasound B-lines in the evaluation of interstitial lung disease in patients with systemic sclerosis: Cut-off point definition for the presence of significant pulmonary fibrosis[END_REF].

In a concrete and practical approach, US may be one step ahead of many assessment tools as it is both widely available and non-invasive/non-irradiating. The issue of reproducibility and observerdependent reliability may be easily addressed when US evaluation is performed by a trained rheumatologist [START_REF] Ruaro | Well-trained sonographers are worth their weight in gold: ultrasound in systemic sclerosis[END_REF]. As musculoskeletal US assessment is now included in many training programs for young Rheumatologists, switching from joints to other organs might be all the easier. Further studies are nonetheless needed to sketch simple and reproducible consensual evaluation criteria, constituting the first step for a coordinate initiative to homogenize practices concerning vascular, skin, lung and musculoskeletal US evaluation in SSc. This step is essential to foster future multi-centric studies.

From a rheumatologic perspective, US assessment may deserve a specific attention as it could have a part to play in the management of SSc in the future. Considering hand evaluation, one of the main challenges for clinicians dealing with SSc, is the multiplicity of causes that lead to hand functional disability [START_REF] Hachulla | SEDUCE study group). Efficacy of sildenafil on ischaemic digital ulcer healing in systemic sclerosis: the placebocontrolled SEDUCE study[END_REF][START_REF] Mouthon | Ischemic digital ulcers affect hand disability and pain in systemic sclerosis[END_REF]: skin roughness limiting mobility; Raynaud phenomenon, digital ulcers, calcinosis or acroosteolysis hindering precise gripping; synovial or tenosynovial involvement responsible for friction rubs and pain… Deciphering the cause of pain itself can also be challenging in SSc patients. Arthralgia or joint pain can result from numerous and co-existing causes, ranging from specific inflammatory synovitis to sclerotic tenosynovitis, mechanical digital ulcers, osteoarthritis or intra-articular calcinosis. Understanding the precise cause of pain and disability is also a central question as it may lead to a guided and more adapted therapeutic management. US offers the unique opportunity of combined and simultaneous dynamic evaluations of vascular, musculoskeletal and skin involvement: a single tool to assess a multicausal impairment [START_REF] Lescoat | Combination of Capillaroscopic and Ultrasonographic Evaluations in Systemic Sclerosis: Results of a Cross-Sectional Study[END_REF].

A Lescoat, C Caimmi and G Coiffier report no conflict of interest. E Hachulla received speaking fees from Actelion Pharmaceuticals Ltd, Glaxosmithkline and Bayer (<€10,000 each) outside of the current study. P Jégo received speaking fees from Actelion Pharmaceuticals Ltd and Bayer No specific funding was received to carry out the work described in this manuscript.

Accepted Article

In addressing the issue of new imaging technics in SSc, one must not forget that clinical evaluation must remain the keystone for the management of SSc patients. Nonetheless, clinical examination may rise unsolved questions that could be answered through US assessment.