Historic Storms Detected in a Changing Environment over Recent Centuries in the Belle Henriette Lagoon

: This paper explores extreme storms with marine ﬂooding events in historical times. It focuses on the challenges associated with detecting these events with a degree of certainty in a changing environment. The paleoenvironmental changes which turned the Gulf of Pictons into the Belle Henriette lagoon have required the reconstruction of events to be based on cored sediments. This study aims to establish how sedimentological washovers can be detected in a depositional environment that has changed profoundly over recent centuries. The Belle Henriette lagoon site is unique. Despite the profound environmental changes that have taken place, the impact of storms was well preserved in the ﬁve cores extracted. Eighteen historical extreme storms with marine ﬂoods were detected in the sediment. By cross-referencing with historical archives, thirteen storms were estimated in 1990, 1940, 1896, 1876, 1859, 1838, 1820, 1811, 1751, 1711, 1645, 1469, and 1351. Five older storms dating back to 1090, 1036, 941, 809, and 581 will also be discussed. By conducting a thorough historical review, we can conclude that these extreme storms caused signiﬁcant damage and had a profound impact on the socioeconomic coastal communities.


Introduction
The risk of extreme storms has always represented a challenge for society and is a matter of growing concern today. Coastal storms and the marine flooding they entail are one of the principal natural disasters which have occurred all over the world throughout the history of mankind, and have been increasing over the last two centuries with the rise in sea level [1]. However, the explanation for this trend on extratropical coastlines is unclear [2]. Quantification is required as to whether (i) the hazard itself has increased due to climate change, (ii) people's exposure to the hazard has increased, or (iii) the definition of what constitutes an extreme storm has changed.
Environmental changes, whether natural or caused by humans, are a significant concern as they are occurring at a new stochastic rhythm and could increase natural hazards [3]. By studying the history of coastal risks and the sedimentology, we can construct a chronology of past storms. Coupled with ancient archives, this allows precision dating of the storms that occurred in previous centuries [4]. However, studying sedimentary records to pinpoint historical storms in a changing environment is a considerable challenge [5]. Detecting washover deposits is a complex issue requiring a pluridisciplinary approach [6]. It is all the more challenging since the environment is changing rapidly and the methods for washover detection must adapt to the different paleoenvironmental configurations of a particular site. This paper explores the ability to study extreme historical storms entailing marine flooding events. Marine flooding is caused by significant meteorological and oceanological parameters that produce intense morphogenic activity in coastal environments [6]. A storm surge occurs when the air pressure drops and there is significant wind, but also when the wave run-up is powerful. More precisely, the study looks at the uncertainties associated with changing environments and how this affects the identification of historical washovers. A significant storm surge can involve washovers, or allochthonous deposits that originate in the marine domain and are delivered to coastal depositional environments [7]. The sedimentological detection of past storm deposits from washovers has increased in the 2000s. The United States was the first to develop this characteristic method in the early 21st century [7,8] and it has since been used worldwide, with some notable examples in Oceania [9], Asia [10], Africa [11], South America [12], and Europe [13][14][15][16]. However, these studies explore environments that have remained constant over time.
This study underlines how new approaches can improve storm reconstructions and briefly highlights some issues associated with environmental changes and decision-support analysis in relation to such extreme events. It reveals a 1000-year storm chronology in the Belle Henriette lagoon (western Marais Poitevin), a highly dynamic environment that has changed considerably over time, based on stratigraphical records coupled with historical archives.

Study Site
The Marais Poitevin is a fragile environment of great ecological interest. It is the second-largest wetland in France and like many wetlands around the world [17] it has a rich biodiversity. It has undergone a profound change in recent centuries. The Marais is an ancient gulf that has been progressively infilled [18] by rivers such as the Sèvre Niortaise and present-day Lay, formerly the Saint Benoist. On the western border of this marsh along the Atlantic Ocean, the Belle Henriette lagoon is a recently-formed four-kilometre-long lagoon covering nearly 200 hectares within a highly dynamic environment (Figure 1).
Cores were extracted from the Anse des Roullières precisely because this is the oldest part of the lagoon, according to historical cartography, and is, therefore, the most appropriate site for studying the paleoenvironmental evolution. This cove is composed of a mudflat colonised by salt marsh vegetation of between 0 m and 2 m NGF ('Nivellement Général de la France', a network of benchmarks in mainland France overseen by the Institut Géographique National www.ign.fr, accessed on 10 December 2021) according to topographic data from the OR2C (Observatoire Régional des Risques Côtiers) team [19]. The lagoon is bordered on each side by a dune ridge culminating at nearly 10 m NGF. The ridge to the southeast isolates the lagoon from the ocean and the one to the northeast isolates it from the anthropized marshes and Lay estuary (Figure 1). On the Atlantic coast, the rise in sea level has slowed since the last marine transgression ended at nearly 7500 BP. In this region, a notable point of inflection is observed in the sea-level curves between 7000 BP and 6000 BP (2.6 mm/year sea-level rise before and 0.8 mm/year after), indicating a relative stabilisation of the coastline [20]. Although the Belle Henriette site has undergone profound geomorphological changes over the centuries (Figure 2), this study shows that the method proposed by the OR2C to detect historical storms remains effective despite these changes.
By crossing sedimentological data with ancient cartography, we identified three main evolutions [18]: (i) the ancient Gulf of Pictons was filled up by the ocean during the transgression linked to the interglacial Holocene, (ii) the filling of the ancient bay that became a marsh as a result of alluvial deposits from the major rivers and the succinct progression of sandy spits along the shoreline; (iii) the current status of the marsh as a lagoon isolated from the ocean. In the Belle Henriette area, the sedimentary cores illustrated variations on a finer scale; (i) the former maritime bay and establishment of the sandy spit behind the lagoon are identified in the deepest facies; (ii) the temporary environment of the Cores were extracted from the Anse des Roullières precisely because this is the oldest part of the lagoon, according to historical cartography, and is, therefore, the most appropriate site for studying the paleoenvironmental evolution. This cove is composed of a mudflat colonised by salt marsh vegetation of between 0 m and 2 m NGF ('Nivellement Général de la France', a network of benchmarks in mainland France overseen by the Institut Géographique National www.ign.fr, accessed on 10 December 2021) according to topographic data from the OR2C (Observatoire Régional des Risques Côtiers) team [19]. The lagoon is bordered on each side by a dune ridge culminating at nearly 10 m NGF. The ridge to the southeast isolates the lagoon from the ocean and the one to the northeast isolates it from the anthropized marshes and Lay estuary (Figure 1). On the Atlantic coast, the rise in sea level has slowed since the last marine transgression ended at nearly 7500 BP. In this region, a notable point of inflection is observed in the sea-level curves between 7000 BP and 6000 BP (2.6 mm/year sea-level rise before and 0.8 mm/year after), indicating a relative stabilisation of the coastline [20]. Although the Belle Henriette site has undergone profound geomorphological changes over the centuries (Figure 2), this study shows that the method proposed by the OR2C to detect historical storms remains effective despite these changes.  (C) filling of the eastern part of the site due to the progradation of the sandy spit B and the continuous filling of alluvial deposits (estimated between the end of Antiquity and the Early Middle Ages); (D) start of the sandy spit C progradation (estimated since AD 750) and presentation of the polderization of the western Marais Poitevin since the 18th century; (E) the "Anse des Terres Noires" temporary silty environment established during the progradation of sandy spit C that has been breached near the Lay estuary (estimated near AD 1820); and (F) implementation of the BelleHenriette current lagoonal environment with the progradation of sandy spit D (estimated at the early 20th century). The map used as a background is named "Carte archéologique du Bas-Poitou" from the Regional Archives of Vendée.

Sedimentological Coring
Five sediment cores were collected over two sedimentary coring campaigns conducted in October 2018 and October 2019, following a transect from south to north (Figure 1). We used one of two methods, depending on the type of environment analysed in the Anse des Rouillères [6]. The first campaign in October 2018 resulted in the collection of two cores (BH-A and BH-B) from the protective dune separating the lagoon from the ocean. An Eijkelkamp motorised vibratory percussive core sampler was used to extract the core from the compact soil [6]. The BH-A core was 170 cm long and the BH-B core 100 cm long. During the second campaign in October 2019, we collected three cores (BH-C, BH-D, and BH-E) from the heart of the lagoon. We used an Eijkelkamp Beeker gravity corer, which is effective for sampling in wet clay or silty soils [21].

Laboratory Analysis
Each core was sliced lengthwise and each half-section photographed and described as it was opened, to preserve the colours of the sediment. To complement the sedimentological description, the Scopix© system (Gwent, UK) was used to take X radiographs [22]. A stratigraphic log was then constructed to describe the various visible facies. The sediments were characterised by lithostratigraphic visual analysis to identify major changes in facies in the sedimentary structure.
A sample was taken from every centimetre of each of the five cores to facilitate fine study of the grain size of all the facies, using a Malvern Mastersizer Hydro 2000G (Palaiseau, France). The calcined sediments were first placed in a solution of distilled water and sodium hexametaphosphate (NaPO 3 ) 6 at 5% for 48 h [23,24]. This allowed all the grains to disperse so the clays were no longer aggregated. For each sample, the software provided the grain size distribution according to 100 percentage grain-size classes, ranging from 0.02 µm to 2000 µm. The data were then transferred to the GRADISTAT V8 software, which transformed the gross particle size distribution into an ordered distribution according to grain size class (clay, silt, sand, coarse elements) [25]. GRADISTAT also calculated several statistical indicators (MGS (Mean Grain Size), sorting, skewness, and kurtosis) according to the Folk and Ward method, which is the most common, hence our choice [26]. The coarse sandy facies revealed an allochthonous marine sediment brought in by intense marine conditions and coarser than the clayey and silty deposits common in lagoon environments [27]; it may also be coarser than the finer sands from a sandy spit or bay environment and a useful proxy for stormy deposits [28].
High-resolution elemental analyses of the split sediment cores were carried out using an Avaatech© XRF core scanner (Dodewaard, The Netherlands). Element intensities were normalised by the total intensity (count per second of each spectrum: cps) [29,30], and element ratios [31,32] were built according to a principal component analysis (PCA) statistical analysis [13]. We used multivariate statistical analysis to investigate the geochemical structure of the sediments quantitatively (continental or marine context data), and to extract the geochemical markers characterising an extreme event of marine origin [13]. A PCA was carried out using R© software to sort these elements into two groups [33], revealing two types of deposit origin: continental and allochthonous (marine) elements. The software sorted the elements automatically following a dendrogram summary of the PCA. The geochemical ratios strontium/iron (Sr/Fe) and calcium/titanium (Ca/Ti) resulting from statistical analysis describe the dominance of exceptional marine conditions.

Core Dating
The age models of the BH-B, BH-C, BH-D, and BH-E cores were derived from a combination of two radionuclides adapted to the multi-decadal time scales. One, 210 Pb (t 1/2 = 22.2 years), has been widely used to calculate short-term (decades) sediment accumulation rates in continental and oceanic environments in the last 40 years [34]. 210 Pb is a naturallyoccurring radioisotope, the atmospheric inputs of which are rapidly incorporated into the sediment. This 210 Pb excess ( 210 Pb xs ) decreases with the depth of the sediment according to its radioactive period and the sedimentation rate. The profiles of 210 P bxs in sediments allow us to estimate sedimentation rates and thus date the sediments. The accuracy of the 210 Pbbased age model must be tested using an independent time-stratigraphic marker, such as 137 Cs (t 1/2 = 30 years), an artificial radioisotope with a pulsed input function: atmospheric nuclear tests (early 1950, maximum fallout in 1963 for the northern hemisphere) and the Chernobyl accident (1986) [35].
Determining the activities of the radioelements of interest ( 210 Pb, 226 Ra, and 137 Cs) is based on non-destructive measurements by direct gamma spectrometry on about 2 g of dry sediment. The detector used was a high-efficiency, broad-energy gamma spectrometer (BeGE TM , Mirion Technologies, Meriden, CT, USA). Detector calibration was carried out using IAEA-certified sediments (RGU). Sediment accumulation rates (SARs) were calculated from 210 Pb xs profiles with depth in sediments. The CF:CS (constant flux, constant sedimentation) model was used in our study as this method is widely employed for lagoon and lake depositional environments [36,37]. The age of sediment layers was obtained by dividing the depth of sediment layers by the SAR. We then estimated the year of deposition for each sediment layer by considering the age of the top core as corresponding to the sampling date (2018 for BH-A and BH-B and 2019 for BH-C, BH-D, and BH-E).
Finally, to estimate the age of the deepest facies where 210 Pb activities are no longer detectable, a marine shell was taken from the BH-C core and dated using the radiocarbon ( 14 C) method. This shell enabled calibration of the sandy facies of the three lagoon cores BH-C, BH-D, and BH-E between 50 and 80 cm deep, since it was collected in an environment common to the three cores. The element extracted was a common Cerastoderma edule shell detected at 62 cm in the BH-C core, dated at 1810 ± 30 years by accelerator mass spectrometry (AMS) [38]. The analysis was carried out by the BETA Analytic laboratory (Dublin, Ireland). The date was corrected using the two atmospheric (IntCal13) and marine (Marine13) curves according to the Reimer reference [39] and the Marine Correction Database based on the work of Tisnérat-Laborde [40].

Historical Sources
The reconstruction of storms that affected western France was also based on ancient primary and secondary sources from historical archives [41][42][43][44], mainly descriptive texts recording data and observations. It involved the manual consultation of 19,888 manuscripts held in various archives: (i) 19,691 unpublished historical manuscripts such as narrative documents (chronicles, annals, journals, and marginal notes in parish registers) instrumental data, old maps and practice documents such as records of accounts, proceedings, repairs, and chancery; and (ii) 196 published articles. Historical archives provide qualitative data (sea state observations, descriptions of damage or societal reactions) and quantitative data (atmospheric pressure measurements, water height measurements) which must be evaluated and analysed using the historical-critical method [45,46]. The earliest marine weather hazard relatively well documented in written sources was from the winter of 1351-1352 (n.st) [47], so the reconstruction could only go back as far as the latter half of the 14th century. In addition, the reconstruction from historical sources ended at the end of the 18th century, when instrumental data was developed and came into general usage. The historical reconstruction of marine storms and submersions in western France from the mid-14th century to the late 18th century confirmed the occurrence of 128 impacting events [46].
These 128 hazards do not represent the reality of the storm activity on the spatiotemporal framework studied, or the totality of the censuses carried out, for several reasons: (1) Historical sources record the occurrence of storms only where their impact disrupts the seigneurial levy, compromises the harvest, or entails unexpected expenses (for repairs) which are subject to tax relief [48]. The research involved in reconstructing historical storms and marine submersions in a particular area is therefore based on retrieving isolated information from documents that were not written or structured for this purpose.
(2) The major documentary gaps due to the loss or destruction of archives (religious wars, the French Revolution, the Second World War) mean that the work of researchers cannot be exhaustive.
(3) Strict criteria such as the date and description of the damage (any events that are undated or include a too-brief mention of the damage must be excluded) were applied for rigorous reconstruction of the events [45]. From the 18th century, newspapers and the Meteo-France database (http://tempetes.meteofrance.fr, accessed on 10 December 2021) were used to complete the storm chronology.

Results
3.1. Stratigraphic Chronologies 3.1.1. Dating of BH-C, BH-D, and BH-E Cores Taken from the Heart of the Anse des Rouillères The three cores BH-C, BH-D, and BH-E had a silty surface about 30-40 cm deep, which allowed detailed profiling of 210 Pb xs and 137 Cs. Significant 210 P bxs activities were measured to depths of 30 cm depending on the site. This enabled calculation of sediment accumulation rates of the same order of magnitude: 0.17 cm/year for the BH-C core, 0.32 cm/year for BH-D, and 0.215 cm/year for BH-E ( Figure 3). This enabled dating of the first few decimetres of the cores, the age of the top of the cores being 2018-the year the cores were taken. To test the dating, the 137 Cs profile of each core was plotted as a function of the time scale based on 210 Pb ( Figure 3). For cores where significant levels of 137 Cs were detected (BH-C, BH-D, and BH-E), there were obvious peaks in 137 Cs activity corresponding to the maximum 137 Cs weapon-test atmospheric fallout in 1963. In addition, the higher sedimentation in core BH-C gave a higher resolution 137 Cs profile, enabling identification of a secondary peak in 1986, ascribed to the Chernobyl nuclear accident. This test based on 137 Cs therefore validates the 210 Pb xs -based dating of the BH-C, BH-D, and BH-E cores, and indicates that there are fairly constant mean sediment accumulation rates at these sites. BH-E cores, and indicates that there are fairly constant mean sediment accumulation rates at these sites. The marine shell extracted at 62 cm from the BH-C core and 14 C dated to 1810 +/− 30 years was used to date the deepest facies of the three cores. This gave an average sedimentation rate of 0.45 cm/year between the base layers dated by 210 Pbxs and the level where the shell was extracted from the BH-C core. This sedimentation rate was higher than that of the superficial facies due to the high proportion of sands present. Since the three cores were extracted from the same environment and had comparable sedimentation rates for the upper facies, we used this average sedimentation rate to extrapolate the date of the deepest facies of the BH-D and BH-E cores ( The marine shell extracted at 62 cm from the BH-C core and 14 C dated to 1810 +/− 30 years was used to date the deepest facies of the three cores. This gave an average sedimentation rate of 0.45 cm/year between the base layers dated by 210 Pb xs and the level where the shell was extracted from the BH-C core. This sedimentation rate was higher than that of the superficial facies due to the high proportion of sands present. Since the three cores were extracted from the same environment and had comparable sedimentation rates for the upper facies, we used this average sedimentation rate to extrapolate the date of the deepest facies of the BH-D and BH-E cores (Figure 3).

Dating of BH-A and BH-B Cores Taken from the Sandy Spit D
The two cores BH-A and BH-B collected from the sandy spit isolating the lagoon from the ocean mainly comprised sand topped by a 10 cm thick layer of silty clay (Figure 1). The silty surface facies of the BH-B core were selected to date these two cores because it was the thickest. The 210 Pb xs activities were quite low and decreased rapidly with depth, reaching negligible levels towards 8-10 cm. This indicated a lower sedimentation rate (0.11 cm/year) than those determined for the cores collected at the heart of the Anse des Rouillères. In the absence of shells that could be used for additional dating, we used this average sedimentation rate to extrapolate the BH-A and BH-B core ages (Figure 3), as their log was homogeneous and comparable (Figures 7 and 8). Note that the dating of the deep layers of these two cores is less accurate than for the BH-C, BH-D, and BH-E cores because the 14 C date added to the BH-C, BH-D, and BH-E depth/age model improved their accuracy. As the dating carried out on these two cores was established only by 210 Pb from the first ten centimetres of the BH-B core, there remains a high level of uncertainty for the deepest layers of BH-A and BH-B.

BH-C, BH-D, and BH-E Cores
The three cores BH-C, BH-D, and BH-E collected in 2019 at the heart of the Anse des Rouillères in a flat, tidal environment ( Figure 1) had a similar lithostratigraphic structure (Figures 4-6). A first facies, common to all three logs, was identified on the surface structure comprising the first 16 centimetres of the cores. This silty layer, rich in organic matter and corresponding to the current surface of the Anse des Rouillères substrate, had an average grain size of about 40 µm heterogeneous (sorting index 5 to 6). This first set was interrupted by several layers of coarse silts at 4 cm (BH-C, dated approximately 1999), 6-7 cm (BH-D dated approximately 1999 and BH-E approximately 1970), and 11 cm (BH-C, dated approximately 1940). Several extreme events occurred in these years in western France [13]. A second facies B was only identified in the log of the BH-C core, composed of medium homogeneous sands of about 400 µm (sand rate up to 90%) between 16 cm (dated 1916) and 22 cm (dated 1908).
The three stratigraphic logs therefore presented comparable parameters between 22 cm (dated around 1896) and 44 cm for BH-E, 46 cm for BH-D, and 48 cm for BH-C (dated between 1820 and 1840). Although the overall composition of facies C seemed similar to facies A, we noted a much higher level of clays, up to 10% in the C facies, while it remained fairly constant and less than 5% in the superficial layer A. An initial C1 facies testified to calm morphogenic conditions with a silty sediment poor in sands, moving towards a more dynamic C2 environment alternating between fine silts, coarse silts, and fine sands.
We noted here a double progression in the dynamics of the C2 layer, with a much more marked heterogeneity in the BH-E core log. We also noted an evolution between the three cores. The thickness was finer for the C2 layer in the BH-E core log. Moving away from the heart of the lagoon, we noted greater heterogeneity and a decrease in the thickness of this layer, which then changed more rapidly into the sandy facies D, estimated around 1820.
Finally, the last facies D showed a deep, sandy layer between 44-48 cm (dated 1920-1840) and the base of the cores at 76 cm, dated 1780-1800. Sands dominated broadly, at more than 90% of the generally well-sorted sediment composition, except for the occurrence of finer sediments (silts or fine sands) observed at two points: 62-68 cm in BH-D and 54-56 cm in BH-E. A distinction can be made between the sandy composition of these three logs. There was a higher proportion of coarse and very coarse sands in the BH-C core, up to 50% of the total sediment composition.

BH-A and BH-B Cores
The two cores BH-A and BH-B collected in 2018 on the sandy spit produced comparable logs (Figures 7 and 8). The logs revealing the structure of the three cores taken from the tidal flat environment. The superficial facies A is a layer dominated by coarse silts (with an Atmosphere 2022, 13, 151 9 of 20 average grain size of between 40 and 60 µm) to a depth of 10-12 cm (dated 1900-1920). This facies was similar to the A layers in BH-C, BH-D, and BH-E, but thinner. This first facies was interrupted by a layer of medium sand (approx. 350 µm) between 5 and 7 cm (dated around 1940). The presence of this sand can be explained by a one-time extreme event already observed in the region [13]. With the exception of this isolated marine event, the first facies A, dominated by silts, was composed of heterogeneous sediments (sorting index of about 5). It was also rich in organic matter since the wet silt sampling site supported submerged salt-marsh vegetation at high tide.

BH-C, BH-D, and BH-E Cores
The three cores BH-C, BH-D, and BH-E collected in 2019 at the heart of the Anse des Rouillères in a flat, tidal environment ( Figure 1) had a similar lithostratigraphic structure (Figures 4-6). A first facies, common to all three logs, was identified on the surface structure comprising the first 16 centimetres of the cores. This silty layer, rich in organic matter and corresponding to the current surface of the Anse des Rouillères substrate, had an average grain size of about 40 µm heterogeneous (sorting index 5 to 6). This first set was interrupted by several layers of coarse silts at 4 cm (BH-C, dated approximately 1999), 6-7 cm (BH-D dated approximately 1999 and BH-E approximately 1970), and 11 cm (BH-C, dated approximately 1940). Several extreme events occurred in these years in western France [13]. A second facies B was only identified in the log of the BH-C core, composed of medium homogeneous sands of about 400 µm (sand rate up to 90%) between 16 cm (dated 1916) and 22 cm (dated 1908).  From 10-12 cm down to the deepest levels, fine to coarse sands dominated broadly in a much more homogeneous sedimentological set (sorting index around 1). The sandy composition reached proportions higher than 90% for most of these facies. Within this set, a sub-distinction was identified between a first facies B, composed of sand gradually passing from fine to medium grain at 20-25 cm (dated 1800-1820), with some clear variations in the sorting index and average grain size (18 cm for BH-A, 14-22 cm for BH-B). The structure of this facies was comparable with the B-facies sediment detected in BH-C. We then identified a second sandy surface, D1, composed of medium sand up to 400 µm down to 135 centimetres (dated 750 AD with extreme uncertainty). A third level, D2, was composed of coarser elements (up to 600 µm), but was also richer in organic matter and marked by a short rupture with a fine facies of silts (147 cm) that was found down to the bottom of the BH-A core. The grain sizes of the 'very coarse' and 'coarse' sand in the D2 facies were higher, between 135 and 170 cm (10% and 50%, respectively) than in the D1 facies between 25 cm and 135 cm (less than 5% and 30%, respectively) and the B facies between 12 cm and 25 cm (very coarse grains almost zero and coarse grains around 30%). The sedimentological characteristics of the D1 facies were comparable with those of the D layers for the three cores BH-C, BH-D, and BH-E. They showed very close similarities to the D facies of the BH-C core, which comprised coarser sands than the D facies of BH-D and BH-E.  The three stratigraphic logs therefore presented comparable parameters between 22 cm (dated around 1896) and 44 cm for BH-E, 46 cm for BH-D, and 48 cm for BH-C (dated between 1820 and 1840). Although the overall composition of facies C seemed similar to facies A, we noted a much higher level of clays, up to 10% in the C facies, while it re- Figure 6. BH-E core log with photography, radiography, grain size results, and geochemical signals.
silts (with an average grain size of between 40 and 60 µm) to a depth of 10-12 cm (dated 1900-1920). This facies was similar to the A layers in BH-C, BH-D, and BH-E, but thinner. This first facies was interrupted by a layer of medium sand (approx. 350 µm) between 5 and 7 cm (dated around 1940). The presence of this sand can be explained by a one-time extreme event already observed in the region [13]. With the exception of this isolated marine event, the first facies A, dominated by silts, was composed of heterogeneous sediments (sorting index of about 5). It was also rich in organic matter since the wet silt sampling site supported submerged salt-marsh vegetation at high tide.  From 10-12 cm down to the deepest levels, fine to coarse sands dominated broadly in a much more homogeneous sedimentological set (sorting index around 1). The sandy composition reached proportions higher than 90% for most of these facies. Within this Figure 8. BH-B core log with photography, radiography, grain size results, and geochemical signals.

Detection of Old Storms by Crossing Sedimentological Indices
Based on statistical analyses, the grain size (mean grain size (MGS), sorting index, and sand rate) and geochemical ratios (Sr/Fe (strontium/iron) and Ca/Ti (calcium/titanium)) were used to determine historical washovers (Figure 9). The scientific community considers concomitant peaks in particle size and geochemical indices as good storm proxies because they indicate a notable marine impulse [4,6]. The most notable event detected in the lagoon sequence was the 1896 event, dated 1899, 1878, and 1880-1910 in the three lagoon cores BH-C, BH-D, and BH-E. In BH-C, a significant sandy marine incursion having sharp contact with the surrounding silts was detected. The MGS, tenth decile, and sand percentages increased at 22 cm ( Figure 9). These increases are much greater than the proxy progressions linked to the other storms. The sand content shows a progression from 40% to 95%, with increases in tenth decile and MGS from 20 µm to 200 µm and from 40 µm to 500 µm, respectively. This is concomitant with significant pulses of S (from 0.55 to 1), Si (from 5 to 40), and Sr (from 1 to 2.3). In BH-D, a less significant grain size marine pulse was recorded, but the geochemical signal remained significant since Si and Sr increased respectively from 7 to 10 and from 7 to 10.5. BH-E was also highly impacted with a succession of sandy pulses (from 40 to 150 µm of MGS) dated between 1880 and 1910. This is characteristic of a high morphogenic environment originating from stormy conditions.

Socioeconomic Impacts of Events Detected in Sedimentological and Historical Records
As developed by Athimon (2021) [45], historical archives enabled us to identify 128 ancient storms between the 14th and 18th centuries. No storm information was identified before 1340 in these sources. Between 1340 and 1400, 8 storms were identified, and 22 occurrences were mentioned during the 15th century. However, the data are not exhaustive and the level of representation of the sources is qualified as 'random' for the 14th and 15th centuries. In the 16th, 17th, and 18th centuries, 31, 33, and 34 different storms were mentioned, respectively, but the sources are only partially representative and not exhaustive for these three centuries. Since the 18th century, local newspapers and the Meteo-France database have recorded nearly a hundred storms that partially or totally affected western France. The sections below relate to the main storms that impacted the Belle Henriette environment, the isotopic dating of which correspond to historical dates found in ancient and recent written archives.

Post-1900 Events
The first storm recorded in the archives refers to the event of 26-28 February 1990. From the historical archives, we learned that three storms hit the study area in 1990, causing a lot of damage. One of them recorded a strong tide gauge of 104 at the port of Pouliguen (SHOM) during the event of 26-28 February 1990. The port is located two kilometres south of the Traicts du Croisic. These three successive high-profile storms caused more than 100 deaths in France, with winds up to 176 km/h in the west. Numerous reports of flooded houses and broken dikes were also recorded. The wind damage was considerable: fallen trees, pylons, walls, cranes, and billboards; collapsed church steeples, torn roofs, crushed cars, and numerous traffic accidents [13].
The event of 16-17 November 1940 was also evident in the cores; historical data reveals the strength of a storm that killed 30 people in 1940 recorded in the sediments. The winds resulted in considerable damage: trees uprooted, towers and dikes destroyed, and major damage to boats and roofs. Numerous marine submersions were also identified.

XIXth-Century Events
Six storms that occurred during the XIXth century were recorded in the five cores. One of the most significant is associated with the event of 4 December 1896. Numerous breaches were identified along the Atlantic coastline. The event resulted in 33 deaths in France due to numerous storm-generated floods. It was combined with a strong tide gauge of about 95 according to the SHOM. We hypothesize that this significant event of 1896 affected much of the French Atlantic coast, since three cores from the Belle Henriette lagoon were affected, and all nine cores from the Petite Mer de Gâvres and Traicts du Croisic lagoons [13]. In the Petite Mer de Gâvres, the event may have breached the coastal spit that isolates the lagoon from the ocean. Through a paleoenvironmental reconstruction of the study area, we can assume that the event contributed to the progradation of the coastal spit D, which resulted in the establishment of the Belle Henriette lagoon.
The five other impacts are less evident in the sediments but were nonetheless considerable according to historical sources. A storm occurred on 31 December 1876 during a syzygy tide impacting the Loire-Atlantique: breach of a 2 km stretch of dike with 6 km of the embankment needing to be rebuilt, submerged marshes, state intervention, and repairs to the salt marshes of Guérande, Mesquer, and Pouliguen totalling FRF 380,000. The cost of the damage caused by the storm at Croisic and Guérande was estimated at over FRF 3 million; we discovered that residents required help to rebuild and it was still being discussed in 1880. During the two-day storm that started on 25 October 1859 with an estimated coefficient of 100, winds reached over 118 km/h. The affected areas were the departments bordering the Atlantic Ocean and the English Channel. We estimated around 800 deaths (90 damaged vessels and 133 individual drownings). Waves invaded fields and salt marshes with the equivalent of 6 years' loss of work [45,49]. On 24 February 1838, a storm called 'vimer' was reported in Charente-Maritime with heavy rains in Loire-Atlanique in the preceding days, which resulted in the Loire flooding [45]. On 18 October 1820, a three-day storm with significant winds affected Brittany, Charente-Maritime, Vendée, and Loire-Atlantique. Several accidents at sea and deaths through drowning (four boat fatalities) were reported, and many boats were destroyed [43,45]. The last storm of the 19th century occurred late in the afternoon on 24 February 1811, with a flood lasting more than 12 h at high tide. Two deaths were reported, along with breaches of the dikes at Le Croisic and Bouin. Half the island was submerged on Bouin, with the loss of salt and salt pans which were under water for more than 12 h. The loss of salt was estimated at around 1/3 of the total production [45].

Older Significant Events
The five older impacts refer to the storms of 17 March 1751, 9 December 1711, 28 January 1645, 27-28 January 1469, and the winter of 1351-1352. All five storms pro-foundly impacted the socioeconomic coastal communities; details of these storms were reported by Pouzet and Maanan (2020) [4].
During the 1751 storm, marine flooding was reported on the former island of Bouin and in Les Sables-d'Olonne. Wind impacts were recorded in Nantes and Thouars and Legé was severely affected, according to the Francheteau family records. In La Ploueze, the damage was so severe that people thought it was accompanied by an earthquake. The storm also had a significant impact off shore; at the request of the consulates, the Admiralty of Nantes issued an order to all boat owners on 17 March 1751 to unload their cargo within 24 h and set sail from Paimboeuf as quickly as possible. Their mission was to assist ships stranded by the storm and recover the merchandise. According to this source, the Admiralty's aim was to reduce the damage to trade caused by the storm. The parish register for Nantes reports that only three of the 60 ships remained unharmed and that many people were injured or died.
During the 1711 event, a strong storm surge was observed in La-Faute-sur-Mer [50]. Historical records report significant damage, such as the destruction of dikes and flooding of salt marshes. The total cost is estimated at 1.19 million euros in modern currency, mainly because the repairs to the dikes took so long. Ré Island was severely flooded. After the storm, salt production tools were discovered near the forecourt of the Portes en Ré church. According to local parish records, the storm and its associated marine flooding profoundly affected coastal families for centuries. One written source called this event a 'hurricane' after the Thouars bell tower collapsed into the church vault.
In 1645, marine flooding was reported in the marshes of the Ré, Aix, and Oléron islands and in the cities of La Rochelle, Marennes, Arvert, and Saint-Sornin. A large vessel was destroyed in the lowlands at an altitude of 3 m NGF and nearly 30 vessels sank. Reported by Athimon (2021), these impacts are recorded in the regional archives of Charente, Charente-Maritime, Gironde, and Sarthe and in local historical books. The storm is referred to as 'cloudy', 'horrible', and 'terrible', and a large proportion of the residents affected referred to it as 'exceptional'. In terms of the economic cost, it took the people of Ré Island three months to repair the damage that affected agricultural activities. The overall spatial extent of the storm was considerable and included much of the French Atlantic coast.
Concerning the event of 1469, the tidal coefficient is estimated as being around 106 (spring tide) on 28 January. The storm in 1469 caused severe damage to the dikes and salt marshes of what was the island of Bouin, which was submerged during the night. Historical data indicates that 1,500 tonnes of salt were likely lost, with 30-m long breaches of the salt marsh pavements. Numerous ridges and roads were destroyed and a salt merchant was drowned. In the Retz region, maritime floods can be detected through the records of taxes that merchants were required to pay to their lord (seigneurial taxes). The land had become barren due to the brutal marine incursions and prices fell as a result. Its effects were mainly felt on the coast (from central Vendée to the south of Brittany).
The last event recorded relates to the winter of 1351-1352. Many lowlands on Noirmoutier Island were flooded for nearly half a century. The marshes were also destroyed in Olonne, and salt production was still affected 15 years later. This event was defined as a major submersion of Ré Island and 'probably apocalyptic'. Athimon (2021) described the storm as 'one of the most violent and dramatic weather events of the last millennium along the French Atlantic coast'. Sedimentological and historical evidence highlights its impact in Brittany to the Pertuis Charentais. They also suggest that this storm may have hit the British Isles, meaning that much of the European Atlantic Coast was affected [4].

Challenges Associated with Crossing Sedimentological Data and Historical Storm Archives
As presented in Figure 9, eighteen marine facies were extracted from the five cores. However, when crossing their isotopic dating with the precise date of storm occurrences recorded in historical archives, only thirteen of them (1990,1940,1896,1876,1859,1838,1820,1811,1751,1711,1645,1469, and 1351) were confirmed as past storm events. The key sedimentological parameters used to detect marine layers are high mean grain size, tenth decile, sand, or coarse sand rates and high geochemical values for Sr/Fe and Ca/Ti. The historical data presented in the results section proves that each of these storms entailed significant marine conditions causing coastal flooding. Significant socioeconomic impacts were also reported, attesting to the significance of these historical disasters. The five other dates (1090, 1036, 941, 809, and 581) go back further than the period covered by the study of historical sources; they cannot be definitively confirmed as historical storms, but the hypothesis can be discussed. According to Bryant (2005) [51], tsunamis affecting the eastern Atlantic coast account for about 1.5% of total tsunamis recorded in the seas and oceans worldwide. The hypotheses of past tsunami signatures cannot be disregarded as they have already been discussed by some European studies, but tsunamis are less likely to occur than storms [52,53].
In addition, there are several challenges associated with crossing these two methods. Reconstructions are based on information gathered, and gaps in the information, especially documentaries, therefore influence the reconstruction, the results, and the interpretation. For example, although sedimentological impacts of past storms can be detected between 1420 and 1470, the historical archives from this period are insufficient to confirm several hypotheses [45]. Sedimentological archives also have wider margins of error than historical ones. Added to this, extratropical events are lower in intensity than tropical cyclones and therefore harder to identify in sedimentary records. A combination of parameters for intense marine-weather events is therefore required to identify impacts on coastal macrotidal depositional environments [6].

Spatial Extent of Storms
Other studies carried out by the OR2C reflect the often vast spatial extent of events detected in the sedimentological archives. The storms of 1940, 1896, 1751, 1469, and 1351 were detected at the Belle Henriette, but also left their mark in the Traicts du Croisic and the Petite Mer de Gâvres [4,13]. According to these two studies, the event of 1990 detected at the Belle Henriette also impacted these two environments, but to a much lesser extent (few cores were impacted). The events of 1711 and 1645 detected in the five cores from the Anse des Rouillères were also detected elsewhere, but only at the Traicts du Croisic, leaving a significant mark [4]. In terms of the events of 1876, 1838, and 1811 at the Belle Henriette, hypotheses were proposed at the Traicts du Croisic during these periods [13]. However, the dating is not accurate enough to support these assumptions. Finally, sedimentological analyses suggest that both the 1859 and 1820 events seem to have affected only the Belle Henriette region, but historical archives insist on how violent it was. Since the Belle Henriette site is located further south than the two others, the two storms may have taken a more southerly route through the country than the other historical events recorded. In terms of the other older hypotheses for 1090, 941, 581, and 809, these can be linked respectively with the hypothetical events detected at the Traicts of Croisic in 1083 ± 30, 973 ± 44, and 522 ± 60, and at the Petite Mer de Gavres in 768 ± 230 [4]. However, the event dated 1036 does not match the work previously carried out by our research team. On a wider scale, the event of 581 can be linked to the 1590 cal yBP storm phase detected on Yeu Island, and all, or nearly all, the events of 581, 809, and 941 are included in an EASE (European Atlantic Stormy Event) [54]. Historical records generally provide limited information on storm damage linked to a precise location. Consequently, it is difficult to document storm tracks from historical sources.

Storms Detected in a Changing Environment
As Figure 2 shows, the Belle Henriette site has undergone profound changes in recent centuries: (i) what was once the Gulf of Pictons was filled by the ocean during the transgression linked to the interglacial Holocene; (ii) it was then filled by alluvial deposits from major rivers and the succinct progression of sandy spits along the shoreline, until the present time, where the environment is (iii) a lagoon cut off from the ocean. The sedimentological approach based on statistical analyses of the geochemical characteristics of the sediments enabled storm events to be tracked in different sedimentary environments (gulf, bay, and lagoon). For this study, we used two ratios associated with grain size: Sr/Fe and Ca/Ti (Sr and Ca are elements of maritime origin; Fe and Ti are from the bedrock) [55][56][57][58].
Several studies have shown that these extreme events of the past are detectable in bays [59][60][61], due to the impulse of significant marine levels revealed by coarser sediments than those in bays. This configuration is identifiable from the 20 cm of cores BH-A and BH-B and the 50 cm of cores BH-C, BH-D, and BH-E. The configurations of the Anse des Terres Noires and the Belle Henriette lagoon are therefore conducive to the detection of washovers since the environment is muddy or silty and extreme events are detected by the incidence of sandy and marine elements (such as Sr and Ca) [6,32,62]. The distinction is therefore much clearer between a lagoon environment and an allochthonous marine facies (rich in Fe and Ti). From these five cores, the lagoon environment was detected in the first 20 centimetres of the BH-A and BH-B cores and the first 50 centimetres of the BH-C, BH-D, and BH-E cores. Regardless of the sedimentary environment, the peaks in these two reports describe an extreme event. This result is confirmed by historical archives. Despite the significant changes to the landscape, the method of detecting historical storms proposed by the Regional Coastal Risk Observatory (OR2C) works with different types of environments.

Conclusions
Despite significant environmental changes in the study area, the proposed approach has produced satisfying results in different types of coastal environments. In addition to providing information on historical storms, the sedimentological archives illustrate the various successive paleoenvironmental phases. From the cores collected at the Belle Henriette lagoon, thirteen ancient storms were detected, mainly through precise analysis of grain size and geochemical proxies. By cross-referencing with historical archives, these storms can be estimated as having occurred in 1990, 1940, 1896, 1876, 1859, 1838, 1820, 1811, 1751, 1711, 1645, 1469, and 1351. Hypotheses can also be discussed for older events dated 1090, 1036, 941, 809, and 581. Historical records are unable to confirm these, but several links with other storm phases have been detected on a wider scale (EASE). By cross-referencing with previous studies carried out at our laboratory, the storms of 1990, 1940, 1896, 1751, 1469, and 1351 were also identified in the Petite Mer de Gâvres (southern Brittany) and Traicts du Croisic, (Loire-Atlantique), those of 1876, 1838, 1811, 1711, and 1645 were also identified at the Traicts du Croisic, and those of 1859 and 1820 were only recorded at the Belle Henriette lagoon. The spatial extent of these historical events can therefore be discussed. Researchers involved in developing risk and uncertainty methods should consider the wider environmental context carefully. Crossing sedimentological data with historical archives is a challenge and is harder to adopt in changing environments. Decision makers should understand why different approaches may give different estimates of risk assessment.