Exploring Offshore Sediment Evidence of the 1755 CE Tsunami (Faro, Portugal): Implications for the Study of Outer Shelf Tsunami Deposits

Outer shelf sedimentary records are promising for determining the recurrence intervals of tsunamis. However, compared to onshore deposits, offshore deposits are more difficult to access, and so far, studies of outer shelf tsunami deposits are scarce. Here, an example of studying these deposits is presented to infer implications for tsunami-related signatures in similar environments and potentially contribute to pre-historic tsunami event detections. A multidisciplinary approach was performed to detect the sedimentary imprints left by the 1755 CE tsunami in two cores, located in the southern Portuguese continental shelf at water depths of 58 and 91 m. Age models based on 14C and Pbxs allowed a probable correspondence with the 1755 CE tsunami event. A multi-proxy approach, including sand composition, grain-size, inorganic geochemistry, magnetic susceptibility, and microtextural features on quartz grain surfaces, yielded evidence for a tsunami depositional signature, although only a subtle terrestrial signal is present. A low contribution of terrestrial material to outer shelf tsunami deposits calls for methodologies that reveal sedimentary structures linked to tsunami event hydrodynamics. Finally, a change in general sedimentation after the tsunami event might have influenced the signature of the 1755 CE tsunami in the outer shelf environment.


Introduction
Research related to tsunamis has increased significantly in the last two decades, focusing not only on retrieving return periods, but also on understanding tsunami event dynamics [1]. Defining the coastal hazard of tsunamis by studying their sedimentary records benefits from their extensive chronology compared to instrumental or historic documentation [2][3][4]. While onshore tsunami deposits are easy are used for the first time to study offshore tsunami deposits in the mid to outer continental shelf environment. The MT analysis is implemented speculatively to explore its potential for tsunami deposits identification in deeper realms. Expanding the possibilities to detect offshore tsunami sediments will support offshore tsunami sediment identification and will provide for additional insights into the dynamics of tsunami sedimentation. Also, tsunami recurrence intervals derived from sedimentary records will be more reliable with a larger set of proven proxies for tsunami signatures in sedimentary records.
Ultimately, exploring offshore tsunami sediments will aid assessment of tsunami hazard and risk. Studying offshore tsunami deposits related to the 1755 CE Lisbon tsunami in the southern Portuguese shelf has advantages as it has well-studied onshore deposits [34,[40][41][42][43][44][45][46] and the fact that, out of seven verified historical tsunamis hitting Portugal, only the 1755 CE tsunami left geological evidence onshore [47]. Hence, this event is not only recognized as the most damaging tsunami that hit Europe since Roman times, but also, its geological imprint is remarkable and the object of numerous studies. This provides a good framework for exploring offshore tsunami deposits that can have small variations to the background sedimentation and small deposit thicknesses.

Study Site
The two studied gravity cores MW14-GC-106 (referred to as M106) (36.89° N, −8.07° W) and MW14-GC-107 (referred to as M107) (36.93° N, −7.99° W) were retrieved in 2014 (during the MOWER project sea mission [48]) from 91 m and 58 m below sea level (bsl), respectively. Core locations lie in the southern Portuguese continental shelf and correspond to the Eastern Algarve proximal continental margin about 15 km and 7 km offshore Quarteira/Faro for core M106 and M107, respectively ( Figure 1). Figure 1. Bathymetry map of the study area in 3D-perspective and core locations. The contour lines are at 100 m intervals. Cores M106 (red, 91 m bsl) and M107 (green, 58 m bsl) were analyzed in the present study, while core VC2B (yellow, 96 m bsl) was analyzed for tsunami imprints by Quintela et al. [19]. CSV = Cape St. Vincent; CSM = Cape Santa Maria. The map was produced using GMT v6.0 [49].
In the study area, the outer shelf passes over to the continental slope at depths of approximately 100 m. The shelf has a slope gradient of 0.3° to 0.4° in the outer continental shelf (core M106) and up to 2° in the inner continental shelf (core M107). The narrow width (ca. 5 km off Faro) and relatively high slope gradients of the continental shelf in proximity to the Cape Santa Maria (CSM) are created by the erosive action of the Mediterranean outflow water (MOW) [50]. The CSM divides the western Figure 1. Bathymetry map of the study area in 3D-perspective and core locations. The contour lines are at 100 m intervals. Cores M106 (red, 91 m bsl) and M107 (green, 58 m bsl) were analyzed in the present study, while core VC2B (yellow, 96 m bsl) was analyzed for tsunami imprints by Quintela et al. [19]. CSV = Cape St. Vincent; CSM = Cape Santa Maria. The map was produced using GMT v6.0 [49].
In the study area, the outer shelf passes over to the continental slope at depths of approximately 100 m. The shelf has a slope gradient of 0.3 • to 0.4 • in the outer continental shelf (core M106) and up to 2 • in the inner continental shelf (core M107). The narrow width (ca. 5 km off Faro) and relatively high slope gradients of the continental shelf in proximity to the Cape Santa Maria (CSM) are created by the Minerals 2020, 10, 731 4 of 30 erosive action of the Mediterranean outflow water (MOW) [50]. The CSM divides the western and eastern continental shelves with distinct surface water currents generated by two gyres, the cyclonic eddy off Cape São Vicente (CSV) and a coastal counter current (CCC), respectively [51]. With strong easterlies, the CCC bifurcates off CSM and invades the western part of the shelf [52]. Further, in the study area 70% of the waves originate from SW and 30% from SE with maximal significant wave heights of 7 m, although the wave climate is characterized by a predominantly smooth and moderate sea state [53]. The main contributors of sediments to the shelf off Algarve are cliff erosion [54] and the input of rivers, mainly the Guadiana River [55]. The red sand cliffs (Pleistocene) at Vale de Lobo consist mainly of quartz and clay minerals with accessory goethite, feldspars, and Ti-oxides. Illite and kaolinite dominate the clay fraction [56]. Generally, only the fine particles are deposited on the shelf because sands are directly transported eastwards by alongshore currents. However, the fine-grained material is to 86% transported further beyond the outer shelf, while the remaining 14% is deposited on the shelf [57]. The present high stand stage hydrodynamic conditions (since ca. 4.0 ky) favor depositional processes and add sediments to the continental shelf sedimentary cover [50]. Cores M106 and M107 are located in the outer shelf and inner shelf sub-domains, respectively [50]. The inner shelf is mostly covered by quartz-bioclastic sands with varying grain-size, while the mid to outer shelf is covered by mud and sandy-mud [58]. The main source of these distal aggradational deposits is interpreted to be point sources far from this area and controlled by along-shelf currents and a reduced local fluvial supply [59].

Previous Study
Some results of a previous study within the ASTARTE project (www.astarte-project.eu) are described and have been used to analyze specific layers ( Figure 2) in order to study possible 1755 CE tsunami deposits in a more sophisticated analysis (present study). All results that are not shown here can be found in the ASTARTE deliverable report 2.43 [60].
Both cores under study show a general grain size fining upward trend indicative of a transgressive succession with decreasing energy depositional regimes from the bottom to the top ( Figure 2).
In core M106, several lithologic units can be identified, a slightly clayey sandy silt unit between 354 cm and 130 cm below sea-floor (bsf), slightly sandy slightly clayey silt between 130 and 9 cm bsf, and slightly sandy clayey silt until the core top ( Figure 2a, lithology classes as defined by [61]). Within these units, the sediment is very uniform, with no visible structures or laminations that, however, are observed in the X-rays imagery (not shown here, see [60]). Frequently, shell fragments, mainly bivalves, are visible throughout the core. The mean grain size ranges from ca. 35 µm at the base to ca. 11 µm at the upper core section. This general decrease in mean grain size entails marked and successive decreases at specific core levels (at ca. 316, 263, 215, 133, and 73 cm bsf) with minor oscillations between these levels. The sand percentage mimics the mean grain size variation and varies between 44% at the base and 5% at the top of the core. The silt fraction is always the textural predominant type with a mean of 62% (max: 75%; min: 44%), exhibiting several peaks in abundance, in general, coincident with mean grain size decreases. Clay is more constant along the core, with an average of 14% (max: 20%; min: 9%).
Core M107 is coarser than M106 but similarly does not show clear visible structures or laminations (Figure 2b), although they can be observed at X-ray images (not shown here, see [60]). In core M107, the laminations are not as well defined as in core M106 and it seems intensely bioturbated with abundant shells and shell fragments. The main feature is a moderate yellowish brown (10YR 5/4) top layer with finer sediments until 20-25 cm depth, whose base is not well defined, and there is a gradual transition to a coarser sand and to dark yellowish brown (10YR 4/2) sediments. Lithologically, core M107 is characterized by alternating layers of slightly gravelly sandy mud, very slightly gravelly sandy mud, and muddy sand. Above a depth 38 cm bsf, some sandy mud layers appear while the top is exclusively constituted by an alternation of sandy mud and slightly sandy mud. The mean grain size is, in general, much coarser than core M106, ranging from 212 µm at the base to ca. 14 µm at the top of the section. The upward fining trend is intercepted between 35 and 30 cm bsf by a peak that reaches 216 µm. Sand fraction is the predominant textural type from the base until 20 cm bsf, where it is possible to observe a gradual decrease (from 75% to 42%). Above this depth, the sand decrease is marked until the top, where it reaches 21%. Silt has a profile almost inverse of the sand's profile, varying from 15% at the base until 59% at the top (with an average of 32%). Clay presents a similar profile to silt, with a lower percentage (average 11%; min 6% and max 25%). Gravel is also present along the core with minor fluctuations and an average percentage of 3% and a maximum of 14%.
X-ray fluorescence (XRF) on cores M106 and M107 was performed with a 1-cm resolution and the identified elements were: Al, Si, S, Cl, K, Ca, Ti, V, Cr, Mn, Fe, Ni, Zn, Ga, Ge, As, Br, Rb, Sr, Y, Zr, Cs, Ba, La, Ce, Nd, Sm, and Pb. Globally, both cores XRF-logs show general increasing Si/Ca trends from the base to the top, indicating increased carbonate facies at the base, also shown on the carbonate contents ( Figure 2). This general trend reflects the increased shells content in the very coarse sand chaotic facies, also noticeable in the Ca and Sr profiles (not shown here, see [60]). In both cores, this trend occurs with two distinct gradients. At the top of the cores (up to 176 cm bsf in M106 and up to 27 cm bsf in M107), this decrease in carbonate content is more pronounced, being more subtle below these depths ( Figure 2). The elements considered to indicate detrital source (Si, Zr, Rb, Al, Fe, and K) show, in general, low variability throughout the cores, indicating a relatively constant source of the sediments into the shelf. Elements considered to be involved in diagenetic (oxidation) processes (Mn, V, Co, S) are all consistent and show increased reactivity in the top 30-40 cm (not shown here, see [60]). Pb, which is usually related with anthropogenic sources, reach higher values at the top and at very discrete layers within the very coarse sand with chaotic facies of abundant shells. Combined analysis of the geochemical (the detrital source elements and the Si/Ca trends) and lithological properties of the cores allowed us to define discrete layers that represent anomalous behavior in relation to the general trends along the cores. In particular, events associated with high energy or abnormal sedimentary sources are being investigated with more detail to evaluate their potential relationship with past tsunami events. Several of these events were identified, and while some correspond to 1-2 cm thick layers, others correspond to 10-20 cm thick layers (see Figure 2). . From left to right: core image, computed tomography scan (only for core M106), core log, textural type percentages, mean grain size (MGS) (μm), CaCO3-content (%), magnetic susceptibility (MS) (10 −4 SI), log transformed XRF elemental counts/ratios (cps). These data originate from a previous study (see Section 1.2,[60]) and show anomalous layers that are interpreted as possible high energy event layers (green horizontal bars). Also, the cores sections of the present study are depicted (core section within the red lines). core image, computed tomography scan (only for core M106), core log, textural type percentages, mean grain size (MGS) (µm), CaCO3-content (%), magnetic susceptibility (MS) (10 −4 SI), log transformed XRF elemental counts/ratios (cps). These data originate from a previous study (see Section 1.2,[60]) and show anomalous layers that are interpreted as possible high energy event layers (green horizontal bars). Also, the cores sections of the present study are depicted (core section within the red lines).
Magnetic susceptibility profiles are quite similar, being possible to divide each core in two main domains, with limits at 175 cm bsf for core M106 and at 30 cm bsf for core M107 (Figure 2). The lower most domain is characterized by lower and quite stable magnetic susceptibility values (around 1 × 10 −7 m 3 ·kg −1 ). Above, values increase significantly, reaching 4 × 10 −7 m 3 ·kg −1 at 100 cm bsf for M106 and at 15 cm bsf for core M107). At the top of the cores, it is still possible to observe a gentle increase with some oscillations. The uppermost part defines a peak, reaching the highest values 6 × 10 −7 m 3 ·kg −1 at 10 and 4 cm bsf for cores M106 and M107, respectively.
Based on the variations in grain-size, XRF, carbonate content and magnetic parameters, several anomalous layers were identified on the basis of significant coincident changes in several properties. Five of these layers were identified on core M107 and seven layers were identified in core M106, all defined by: i) marked grain-size variations (increase or decrease in grain-size, usually expressed in the mean grain-size); ii) increases in terrigenous proxies (Fe/Ca and Si/Ca) or by detrital source proxies (Si, Zr, Rb, Al, Fe, and K) in XRF data; and iii) magnetic mineral phases variations and environmental magnetic parameters (not shown here, see [60]). Most of these layers were 210 Pb and 14 C-dated and the most recent of such layers was dated in both cores, indicating an age that is compatible with the AD 1755 tsunami event. Therefore, this layer was selected for further detailed investigation, the focus of this work.

Materials and Methods
The focus was given on the specific sections of the two cores, M106 and M107, that incorporate the possible 1755 CE tsunami layer and levels below and above the event layer (pre-and post-tsunami) ( Figure 2). For the detection of differences between the possible 1755 CE tsunami layer and the background sedimentation, MS, SC, and MT were examined, including multivariate statistical analysis.

Supporting Data
Supporting data, i.e., data obtained from previous works include grain size, carbonate content and inorganic geochemistry for cores M107 and M106 from a previous study (Section 1.2). The data are available for the whole cores at a minimum resolution of 1 cm.
Grain size and carbonate content analyses were performed at IPMA-Tavira geological lab (Tavira, Portugal). Grain size distribution was determined by laser diffraction with a Malvern Mastersizer Hydro 2000G, following the necessary protocol and the adequate software (v.5.61). The resulting grain size data are a volume percentage for all the textural distributions ranging from 0.01 µm to 2 mm in 1 cm intervals. In the coarser core M107, a 500 µm sediment separation was made beforehand to yield more representative and reliable grain size analysis by analyzing the fraction <500 µm by laser diffraction and >500 µm by classical sieving. Textural type percentages and mean grain size (geometric methods of moments) were calculated using Gradistat (v8.0) [62]. Carbonate content was analyzed with the volumetric method of Scheibler as described by Tatzber et al. [63] using an Eijkelkamp calcimeter.
MS was continuously measured on whole-round cores with a GEOTEK Multisensor Core Logger (MSCL) at the University of Vigo with a frequency of 1/cm and in a discrete way (mass normalized magnetic susceptibility χ LF ) through cylindrical samples (one centimeter of height and diameter) measured in the magnetic susceptibility bridge MFK1 (AGICO), at the Laboratory of Paleomagnetism and Geomagnetism of Instituto Dom Luiz (University of Lisbon, Lisbon, Portugal).
XRF data were collected every 10 mm down-core over a 15 mm 2 area by using generator settings of 30 kV, 55 mA, and 20 s directly at the split core surface of the archive half with an ITRAX Core Scanner at the University of Vigo (Vigo, Spain). The split-core surface was covered with a 4 µm thin SPEXCerti Prep Ultralene1 foil to avoid contamination of the XRF measurement unit and desiccation of the sediment. The acquired data have been processed with the software CoreScanner 8.6.3 using the Q-spec 8.6.0 package. Identified elements with sufficient abundance cover Al, Si, S, Cl, K, Ca, Ti, V, Cr, Mn, Fe, Ni, Zn, Ga, Br, Rb, Sr, Y, Zr, Cs, and Nd, and were recorded in counts per second (cps). Counts were normalized over cps to correct surface geometry variations and grain size. The ratio of incoherent over coherent scattering (MoR) is used as an indicator of organic matter content [64]. For terrestrial input, the elemental count of Al + Si is used as a proxy of the Pleistocene red sand (with high quartz and Illite/Kaolinite content), Ti/Ca as a sediment source dependent proxy, and Fe of the input of fine material [17].
Computed tomography (CT) data for core M106 was obtained using a HITACHI ECLOS 16 Multislice CT from the Veterinary University Hospital Rof Codina of Lugo (University of Santiago de Compostela, Santiago de Compostela, Spain) following the acquisition protocol of Mena et al. [65]. The final resolution per voxel is of 0.2 mm × 0.2 mm × 0.65 mm (x, y, z). Images from the CT-scan were obtained using the free software "MRIcro" [66]. The radio-density profile was extracted, using the free software "anidoC" [65], from a centered cylinder with 5 cm of diameter, trying to avoid the interaction between the plastic material of the core and the sediment.

Sand Composition and Microtextural Analysis of Quartz Grains
Samples for SC and MT analysis were wet sieved using a 63 µm mesh, deflocculated (with sodium hexametaphosphate, Na 6 O 18 P 6 ), dried (max. 40 • C), and sieved again in a 1 Φ interval (logarithmic scale after Krumbein [67]). The SC analysis incorporates 100 grains of the 250-125 µm fraction of each sample, counted, and identified under a binocular microscope (SZX7, OLYMPUS, Tokyo, Japan). A preliminary study was performed to statistically determine the count number and best-suited fraction for the SC analysis [68]. The grains are classified as terrigenous, biogenic, or non-identified and subdivided into 11 components (Table 1). Terrigenous components are composed of quartz, mica, opaque, aggregate, and other terrigenous components. In this group, aggregates are purely terrigenous without biogenic elements. Biogenic components are subdivided in planktonic and benthic foraminifera, as well as mollusk, terrestrial biogenic (continental origin, e.g., charcoal, wood), and other biogenic components (such as echinoderms). Not identified are components where visual identification with the binocular microscope cannot yield a clear classification for the terrigenous or biogenic group. Aggregated components consisting of biogenic and terrigenous elements are classified as non-identified.
Quartz grains of the size between 250-125 µm were randomly handpicked from 9 samples of core M107 using a binocular microscope representing the 3 different chronological regimes pre-tsunami, tsunami, and post-tsunami. A population of 30 grains for each sample was chosen to gain a meaningful sedimentological interpretation based on an empirical test conducted by Costa et al. [34]. The grains were mounted on adhesive carbon stickers and high-resolution images were photographed with a scanning electron microscope (SEM) (TM4000Hitachi, Tokyo, Japan). The quartz grain shape was evaluated based on the grey level SEM-images using the free software ImageJv.1.49 [69,70]. An image processing routine was established to separate the foreground (grain) from the background based on Sobel edge detection and the Yen-thresholding method [71]. The binarized output image (Figure 3c) was used to determine grain perimeter, grain area, circularity, aspect ratio, roundness, and roughness factor. Surface microtextures including dissolution, adhering particles, percussion marks, and fresh surfaces were identified visually by marking areas of occurrence with polygons on the grain surface ( Figure 3b). Their total percentages relative to the total surface area occupation were calculated and the results of automated routine for grain shape measurements were processed using R [72]. The automated image processing and mathematical morphology calculations, as well as the coding of the grain-image names prior processing reduces the human operator subjectivity. Shape parameters, including perimeter, circularity, and aspect ratio, were obtained using "ImageJ", allowing the calculation of the grain roundness parameter R following the methodology and equation introduced by Takashimizu and Iiyoshi [73]. The roughness factor (FR) in this study denotes the irregularity of the grain outline resembling an angularity scale as that of Powers [74] and the methodology is based on the work of Lira and Pina [75], in which a more detailed description can be found. Following the semi-quantitative approach of Costa et al. [34], the focus was given to 5 microtextural families of grain surface features of chemical and physical nature: angularity, fresh surfaces, percussion marks, adhering particles and dissolution ( Figure 3). This allows possible discrimination with statistical reliability due to the greater sampling compared to the analysis of single microtextures. To classify the microtextural area of individual grains each textural surface feature was scaled as 0 (absent), 1 (0-10% of grain surface), 2 (10-25% of grain surface), 3 (25-50% of grain surface), 4 (50-75% of grain surface), and 5 (>75% of grain surface).

210. Pb and 14 C Dating
The age models of both cores are based on accelerator mass spectrometry (AMS) 14 C-dating at Beta Analytic Inc. (Miami, FL, USA) and Actlabs (Canada) and 210 Pb-dating checked by 137 Cs activity at University Bordeaux (Pessac, France) ( Table 2).
Ages obtained from the age estimation analysis were projected on core depths using Bayesian statistics with the help of the "rbacon" package (v.2.4.2) [76] implemented in the software "R" [72]. This methodology divides the core into small vertical sections (here 5 cm) and uses Markov chain Surface microtextures including dissolution, adhering particles, percussion marks, and fresh surfaces were identified visually by marking areas of occurrence with polygons on the grain surface ( Figure 3b). Their total percentages relative to the total surface area occupation were calculated and the results of automated routine for grain shape measurements were processed using R [72]. The automated image processing and mathematical morphology calculations, as well as the coding of the grain-image names prior processing reduces the human operator subjectivity. Shape parameters, including perimeter, circularity, and aspect ratio, were obtained using "ImageJ", allowing the calculation of the grain roundness parameter R following the methodology and equation introduced by Takashimizu and Iiyoshi [73]. The roughness factor (FR) in this study denotes the irregularity of the grain outline resembling an angularity scale as that of Powers [74] and the methodology is based on the work of Lira and Pina [75], in which a more detailed description can be found. Following the semi-quantitative approach of Costa et al. [34], the focus was given to 5 microtextural families of grain surface features of chemical and physical nature: angularity, fresh surfaces, percussion marks, adhering particles and dissolution ( Figure 3). This allows possible discrimination with statistical reliability due to the greater sampling compared to the analysis of single microtextures. To classify the microtextural area of individual grains each textural surface feature was scaled as 0 (absent), 1 (0-10% of grain surface), 2 (10-25% of grain surface), 3 (25-50% of grain surface), 4 (50-75% of grain surface), and 5 (>75% of grain surface).

210. Pb and 14 C Dating
The age models of both cores are based on accelerator mass spectrometry (AMS) 14 C-dating at Beta Analytic Inc. (Miami, FL, USA) and Actlabs (Canada) and 210 Pb-dating checked by 137 Cs activity at University Bordeaux (Pessac, France) ( Table 2).
Ages obtained from the age estimation analysis were projected on core depths using Bayesian statistics with the help of the "rbacon" package (v.2.4.2) [76] implemented in the software "R" [72]. This methodology divides the core into small vertical sections (here 5 cm) and uses Markov chain Monte Carlo (MCMC) iterations to estimate the sedimentation rates for each section. Thus, no constant sedimentation rate between dated levels must be assumed, resulting in better sedimentation rate estimations [77]. Radiocarbon ages were corrected using the calibration curve "Marine13" [78] coupled with a calculated marine reservoir effect value (∆R) for the core location ( Figure 1). ∆R was established by using the database of "CALIB" (http://calib.org/marine/, 6 September 2019) entering the approximate core location and choosing the 20 closest points. Out of the 20 listed results, only points from the south coast of Portugal were chosen, excluding also points originating from freshwater influenced areas as the Ria Formosa lagoon. This yield ∆R = 268 ± 122 based on the average of 9 points out of this database ( [79], http://calib.org/marine/references.php). Table 2. Radiocarbon ages obtained from the 2 cores M106 and M107. Calibrated radiocarbon ages are minimum and maximum values of the calibrated ages with "rbacon" [77] (see text in chapter 2.3. for more details). BP = before present (i.e., before 1950 CE); CE = common era; BCE = before common era; prob = probability. Excess 210 Pb ( 210 Pb xs ) was calculated by subtracting the supported fraction, 226 Ra from the total 210 Pb activity. Measurements were done downcore using a low-background, high efficiency gamma detector until reaching negligible excesses. Then, 210 Pb xs profiles of each cores were used to calculate mean sedimentation rates using the CFCS (constant flux constant sedimentation) model [80]. This age-depth model was computed using R package "serac" (v.0.1.0) [81]. In contrast to 14 C-ages, 210 Pb-ages were used as absolute dates within the final age-depth model, with error estimates yielded from the CFCS model, and then combined in the same age model using "rbacon". Previously detected layers, that are interpreted to correspond to possible high energy events (see Section 1.2, Figure 2) were speculatively considered as sections of the instantaneous deposition of sediments. Hence, these layers were excised before modeling to further enhance the age-depth relations, because ignoring any age gaps may cause even more error in the age model [82]. Also, parameters for the age model vary from the basic settings of the "rbacon"-function (see Figure 4 for the values). Age-depth relation using 210 Pbxs and 14 C method ("rbacon" output) including calibrated 14 C-ages (blue). The red curve shows a single "best" model based on the weighted mean age for each depth. Darker grey indicates CE-ages that are more likely and dashed grey lines represent the 95% confidence intervals (max. and min. CE-ages). Horizontal grey bars represent possible high energy event layers (abrupt sedimentation) detected by a previous study (see Section 1.2,[60]). Age model parameters are shown for better reproducibility.

Statistical Analysis
Besides visualization of the supporting data, a linear discriminant analysis (LDA) was conducted incorporating selected geochemical and sedimentological data as an independent set of quantitative variables and the chronological regime grouping of samples as a response variable that was obtained independently from the explanatory variables. LDA can be used to determine to what degree the set of quantitative variables can explain the sample grouping [83]. Data of the LDA were checked for multivariate homogeneity of within-group covariance matrices (permutational test, pvalue > 0.05) and for distinct means of the explanatory variables (Wilks' lambda test, p-values < 0.05), as well as centered and scaled. The LDA was performed using package "vegan" (v.2.5-6) [84] implemented in "R" [72]. A principal component analysis (PCA) using the function "rda" (basic settings) of the package "vegan" was conducted incorporating a hierarchical cluster analysis (HCL) (Ward's method) and including samples of both cores. The PCA data matrix was standardized and consists of the supporting proxies along with textural type percentages. Age-depth relation using 210 Pb xs and 14 C method ("rbacon" output) including calibrated 14 C-ages (blue). The red curve shows a single "best" model based on the weighted mean age for each depth. Darker grey indicates CE-ages that are more likely and dashed grey lines represent the 95% confidence intervals (max. and min. CE-ages). Horizontal grey bars represent possible high energy event layers (abrupt sedimentation) detected by a previous study (see Section 1.2,[60]). Age model parameters are shown for better reproducibility.

Statistical Analysis
Besides visualization of the supporting data, a linear discriminant analysis (LDA) was conducted incorporating selected geochemical and sedimentological data as an independent set of quantitative variables and the chronological regime grouping of samples as a response variable that was obtained independently from the explanatory variables. LDA can be used to determine to what degree the set of quantitative variables can explain the sample grouping [83]. Data of the LDA were checked for multivariate homogeneity of within-group covariance matrices (permutational test, p-value > 0.05) and for distinct means of the explanatory variables (Wilks' lambda test, p-values < 0.05), as well as centered and scaled. The LDA was performed using package "vegan" (v.2.5-6) [84] implemented in "R" [72]. A principal component analysis (PCA) using the function "rda" (basic settings) of the package "vegan" was conducted incorporating a hierarchical cluster analysis (HCL) (Ward's method) and including samples of both cores. The PCA data matrix was standardized and consists of the supporting proxies along with textural type percentages.
The Hellinger-transformed (abundance values divided by sample total abundance with following square-root transformation) count data [85] of the SC were further analyzed using PCA for each core (function "rda", basic settings).
P-values for the MT analysis were calculated to test if there are statistical differences in mean values of the 3 populations (pre-tsunami, tsunami, post-tsunami) aiding interpretation of results by using Kruskal-Wallis test [86] from the base "R"-package "stats" (v.4.0).

Chronology
Based on the depth decrease in 210 Pb xs activities and the CFCS model, a mean sedimentation rate of 0.92 mm was estimated in the uppermost 20 cm of core M107, and 2.75 mm·yr −1 in the uppermost 10 cm of core M106.
In general, sedimentation rates varied over time in both cores with lower rates in lower core depths and higher rates in the upper core sections (Figure 4b,c). Notwithstanding, at the base of core M107, the sedimentation rate is apparently high. Moreover, the removing of the possible high energy events (abrupt sedimentation) in the age model is speculative, and their limits must be investigated more precisely. In the top cores, the artificial 137 Cs had permitted to validate the chronology. The onset of 137 Cs is observed early 50s, as expected, and the maximum activities in subsurface are in agreement with the known inputs events (maximum nuclear weapon test atmospheric fallout in 1963 CE and Chernobyl in 1986 CE), the low sedimentation rate in core M107 do not permit to identify the two peaks (Figure 4a,d). Thereafter, the age model presents increasing uncertainties with greater distance from dated levels which is due to the relatively low number of radiocarbon ages. Yet, core M107 includes a longer sedimentary sequence with suggested ages until 3785 BCE compared to core M106 that covers mean ages until 1192 BCE. Consequently, the mean sedimentation rate of the entire core length is ca. 1.09 mm·yr −1 for core M106 and 0.21 mm·yr −1 for core M107. Estimated ages for the possible 1755 CE tsunami layers is 1769 CE single best model age with a 95% confidence interval between 1628 CE and 1871 CE for core M106 and single best model age of 1589 CE for core M107, while the 95% confidence interval comprises ages between 1326 CE and 1794 CE.

Multi Proxy Data
Both cores show a general grain size fining upward trend indicative of a transgressive succession with decreasing energy depositional regimes from the bottom to the top of the studied core sections ( Figure 5). The mean grain size of core M106 is ranging from 19 µm at the base to ca. 13 µm at the upper core section. The mean grain size of core M107, which is generally coarser than core M106, ranges from 200 µm at the base to ca. 40 µm at top of the section. In core M107, sand is the predominant textural type, while in core M106 silt is predominant with percentages of 40-70% ( Figure 2; Section 1.2). Gravel is absent in core M106. This is coherent with the results of Lobo et al. [59], where core M106 represents the distal muddy belt environment and core M107 the infralittoral propagating wedge. The base of the pre-tsunami layer includes the highest values in mean grain size in both cores. The tsunami layer contains mean grain size peaks in both cores. However, it is subtle in core M106 (at core depth 71 cm and 69 cm) and more pronounced in core M107 (at core depth 23 cm and 20 cm), where it is directly linked with the increase of gravel content. In the post-tsunami layer, mean grain size is decreasing with minor variations. Carbonate content is varying between 10% to 11.5% and 12% to 20% in cores M106 and M107, respectively ( Figure 5). Peaks of carbonate content in the tsunami layer are visible in both cores. In contrast to core M107, core M106 has higher levels of carbonate content in the upper tsunami layer which remains high until 58 cm core depth in the post-tsunami section. MS is generally higher in core M106 compared to core M107 with values ranging from 4.2 to 5 (10 −4 SI) and 1 to 4.5 (10 −4 SI), respectively ( Figure 5). In core M106 between 70 cm to 55 cm, MS is increased with a distinct peak at depth 63 cm, which is the uppermost sample of the tsunami layer. In core M107, shows a tenuous increase in the upper tsunami layer and a pronounced peak in the post-tsunami layer around depth 11 cm. Considering mass-normalized magnetic susceptibility χLF for both cores (including all samples) ( Figure 6) it is verified for core M106 that pre to post tsunami depths are within an interval that mostly defines a relative plateau of magnetic susceptibility values between 100 and 32 cm compared to adjoining lower and upper core sections. Despite some oscillations, such an interval seems to interrupt the increase of susceptibility observed below 100 cm and above 32 cm. Within this range, it is possible to see some abrupt variations that are distinguished MS is generally higher in core M106 compared to core M107 with values ranging from 4.2 to 5 (10 −4 SI) and 1 to 4.5 (10 −4 SI), respectively ( Figure 5). In core M106 between 70 cm to 55 cm, MS is increased with a distinct peak at depth 63 cm, which is the uppermost sample of the tsunami layer. In core M107, shows a tenuous increase in the upper tsunami layer and a pronounced peak in the post-tsunami layer around depth 11 cm. Considering mass-normalized magnetic susceptibility χ LF for both cores (including all samples) ( Figure 6) it is verified for core M106 that pre to post tsunami depths are within an interval that mostly defines a relative plateau of magnetic susceptibility values between 100 and 32 cm compared to adjoining lower and upper core sections. Despite some oscillations, such an interval seems to interrupt the increase of susceptibility observed below 100 cm and above 32 cm. Within this range, it is possible to see some abrupt variations that are distinguished from low amplitude oscillations, such as those observed between 93 cm and 87 cm, the positive peak at 63 cm, and the negative one at 54 cm. For core M107, the pre-tsunami range is characterized by very stable, gradually increasing magnetic susceptibility values, the tsunami one by a sharp increase, and the post-tsunami by oscillations starting at 15 cm that interrupt the general increase observed below. from low amplitude oscillations, such as those observed between 93 cm and 87 cm, the positive peak at 63 cm, and the negative one at 54 cm. For core M107, the pre-tsunami range is characterized by very stable, gradually increasing magnetic susceptibility values, the tsunami one by a sharp increase, and the post-tsunami by oscillations starting at 15 cm that interrupt the general increase observed below.
XRF data show a significantly larger range in core M107 compared to core M106, except for the Al+Si elemental counts, where variability is similar ( Figure 5). Especially, the Ti/Ca ratio is higher in core M106 compared to core M107 which might be related to the higher mean grain size values in core M107. Out of all XRF proxies used in the present study, the Al + Si elemental counts show increased mean values in the tsunami layer relative to the post-and pre-tsunami units. In core M106 mean values are 6.15, 6.21, and 6.17 for post-tsunami, tsunami, and pre-tsunami units, respectively. In core M107, post-tsunami, tsunami, and pre-tsunami units have mean values of 6.14, 6.28, and 6.24, respectively. CT-scan data of core M106 show a decreasing radio-density trend (from lighter to darker colors) ( Figure 5). For the studied section (80-52 cm), CT shows a decreasing radio-density in the pre-tsunami layer (blue colors, 900-840 HU) and a relatively abrupt change in the tsunami layer, from low radiodensity materials (ca. 780 HU) to increasing amount of higher radio-density materials (ca. 860 HU) (Figure 7). The post-tsunami layer is characterized by material of relatively constant high radiodensity material (ca. 865 HU) with a scarce presence of low radio-density material (down to 810 HU). XRF data show a significantly larger range in core M107 compared to core M106, except for the Al + Si elemental counts, where variability is similar ( Figure 5). Especially, the Ti/Ca ratio is higher in core M106 compared to core M107 which might be related to the higher mean grain size values in core M107. Out of all XRF proxies used in the present study, the Al + Si elemental counts show increased mean values in the tsunami layer relative to the post-and pre-tsunami units. In core M106 mean values are 6.15, 6.21, and 6.17 for post-tsunami, tsunami, and pre-tsunami units, respectively. In core M107, post-tsunami, tsunami, and pre-tsunami units have mean values of 6.14, 6.28, and 6.24, respectively.
CT-scan data of core M106 show a decreasing radio-density trend (from lighter to darker colors) ( Figure 5). For the studied section (80-52 cm), CT shows a decreasing radio-density in the pre-tsunami layer (blue colors, 900-840 HU) and a relatively abrupt change in the tsunami layer, from low radio-density materials (ca. 780 HU) to increasing amount of higher radio-density materials (ca. 860 HU) (Figure 7). The post-tsunami layer is characterized by material of relatively constant high radio-density material (ca. 865 HU) with a scarce presence of low radio-density material (down to 810 HU). For the multivariate statistical analyses a data matrix was created consisting of the supporting data set including magnetic susceptibility (log(MS)), mean grain size (log(MGS)), textural type percentages (sand, silt, clay) carbonate content (CaCO3), organic matter content (MoR-ratio), specific XRF elemental counts and ratios (log(MoR), log(Ti/Ca), log(Fe), log(Al + Si) ( Figure 5), as well as counts of the SC (mollusks (Moll), foraminifera (benthic + planktonic)). PCA, LDA, and HCL were conducted on these data matrices with slightly different variables (In the LDA the textural types are excluded and SC counts and log(Al + Si) are excluded for PCA and HCL). The results of the analyses are graphically presented in Figure 8. LDA results show that the pre-defined chronological layers can be explained by the data set, although the grouping of samples is clearer in core M107 than in core M106 (Figure 8a,b). Here, post-tsunami samples and the tsunami samples are overlapping, indicating a more similar sedimentary characteristic (>90% of correct classification), and sample 61 (depth 61 cm to 62 cm) is classified in the tsunami-group. However, higher Al + Si-content, more abundant mollusk fragments, lower Fe-content, and total foraminifera counts are separating the post-tsunami and tsunami samples. Mean grain size is increasing with core depth and is the main factor responsible for the distinction between the pre-tsunami samples and the two younger units. In core M107, the increase of mean grain size is the main factor for the differentiation of the tsunami layer. Higher Al + Si-content, lower abundance of foraminifera, and lower Ti/Ca ratio is contributing to a minor degree to the differentiation. Post tsunami samples show higher MS and carbonate-content compared to the lower layers. Pre-tsunami samples have higher organic matter (MoR-ratio) content and low MS values.
The PCA includes samples of both cores and depicts general differences between the cores besides information regarding the possible tsunami imprint (Figure 8d). Core M107 differs from core M106 with higher proportions of sand-sized particles, higher carbonate content, and higher mean grain size, while Fe-content, Ti/Ca-ratio, organic matter content, and MS are lower compared to core M106. Also, the variability between the samples is higher in core M107 as indicated by the higher variation in PCA loadings. In general, core M107 has a greater intravariability of the studied sediments in their sedimentary characteristic compared to core M106. Nevertheless, differences in the sedimentology between the three chronological units (post-tsunami, tsunami, and pre-tsunami) are visible in both cores using PCA. Results in respect of the distinction between the chronological For the multivariate statistical analyses a data matrix was created consisting of the supporting data set including magnetic susceptibility (log(MS)), mean grain size (log(MGS)), textural type percentages (sand, silt, clay) carbonate content (CaCO 3 ), organic matter content (MoR-ratio), specific XRF elemental counts and ratios (log(MoR), log(Ti/Ca), log(Fe), log(Al + Si) ( Figure 5), as well as counts of the SC (mollusks (Moll), foraminifera (benthic + planktonic)). PCA, LDA, and HCL were conducted on these data matrices with slightly different variables (In the LDA the textural types are excluded and SC counts and log(Al + Si) are excluded for PCA and HCL). The results of the analyses are graphically presented in Figure 8. LDA results show that the pre-defined chronological layers can be explained by the data set, although the grouping of samples is clearer in core M107 than in core M106 (Figure 8a,b). Here, post-tsunami samples and the tsunami samples are overlapping, indicating a more similar sedimentary characteristic (>90% of correct classification), and sample 61 (depth 61 cm to 62 cm) is classified in the tsunami-group. However, higher Al + Si-content, more abundant mollusk fragments, lower Fe-content, and total foraminifera counts are separating the post-tsunami and tsunami samples. Mean grain size is increasing with core depth and is the main factor responsible for the distinction between the pre-tsunami samples and the two younger units. In core M107, the increase of mean grain size is the main factor for the differentiation of the tsunami layer. Higher Al + Si-content, lower abundance of foraminifera, and lower Ti/Ca ratio is contributing to a minor degree to the differentiation. Post tsunami samples show higher MS and carbonate-content compared to the lower layers. Pre-tsunami samples have higher organic matter (MoR-ratio) content and low MS values.

Sand Composition
For the SC analysis, a total of 51 samples of the two cores were analyzed. By counting and identifying 100 grains per sample a total of 5100 random grains were processed, while 159 grains were classified as non-identified (ca. 3.1% of total). Overall, sand components show minor fluctuations in their percentages in the studied sections in both cores (52-80 cm in core M106 and 8-32 cm in core M107) (Figure 9). The PCA includes samples of both cores and depicts general differences between the cores besides information regarding the possible tsunami imprint (Figure 8d). Core M107 differs from core M106 with higher proportions of sand-sized particles, higher carbonate content, and higher mean grain size, while Fe-content, Ti/Ca-ratio, organic matter content, and MS are lower compared to core M106. Also, the variability between the samples is higher in core M107 as indicated by the higher variation in PCA loadings. In general, core M107 has a greater intravariability of the studied sediments in their sedimentary characteristic compared to core M106. Nevertheless, differences in the sedimentology between the three chronological units (post-tsunami, tsunami, and pre-tsunami) are visible in both cores using PCA. Results in respect of the distinction between the chronological units already shown by the two LDAs are mimicked by the PCA. The cluster analysis of the PCA results ( Figure 8c) with a solution of four clusters is separating the analyzed samples in two units, roughly before and after the tsunami event unit.

Sand Composition
For the SC analysis, a total of 51 samples of the two cores were analyzed. By counting and identifying 100 grains per sample a total of 5100 random grains were processed, while 159 grains were classified as non-identified (ca. 3.1% of total). Overall, sand components show minor fluctuations in their percentages in the studied sections in both cores (52-80 cm in core M106 and 8-32 cm in core M107) (Figure 9). A general difference in SC between core M107 and M106 is apparent: (1) In the shallower shelf core M107 higher percentages of terrigenous components (ca. 50%) are present compared to the deeper shelf core M106 (ca. 30%). (2) In core M106, percentages of benthic foraminifera are higher, and percentages of mollusks are lower compared to core M107. The terrigenous part of core M107 is made up of mainly quartz and mica, while opaques, other terrigenous, and aggregates are low in percentages (<8%) (Figure 9b). The biogenic components are mainly other biogenic, mollusk, and benthic foraminifera. Terrestrial biogenic components have very low percentages throughout the studied core section. Although core M106 shows general fluctuations in the SC there is no general trend in the studied section (Figure 9a). The terrigenous part consists mainly of quartz, mica, and aggregates, while opaques and other terrigenous components are underrepresented. A peak of aggregates at depth 73 cm in core M106 is accompanied by lower densities as depicted by the CTdata (Figure 7). The biogenic part is dominated by other biogenic components and benthic foraminifera (up to 36%). Mollusk and terrestrial biogenic components are low in this core with maximal percentages of ca 5%. In the possible tsunami layer, no remarkable peak of terrigenous components can be observed in either cores, that could serve as a backwash tracer. However, higher percentages of mollusk shells/fragments (only in core M106) and terrestrial biogenic components are present in the high-energy event layer with peak percentages within the total studied section, indicative of sediment transportation from beach or nearshore areas to the shelf.
For each core, a separate PCA of the SC was conducted, and results are visualized in Figure 10. A solution of 4 principal components was chosen, and in both cores selected principal components explained >77% of the total variance in the data set. For core M107, selected component assemblages are PC1_107 and PC2_107, explaining 49.9% of the total variance and PC1_106 and PC4_106 for core M106 that explain 50% of the total variance. Both cores show a distinction between the chronological groups by their SC, but differences of the background sedimentation and the possible tsunami units are small. Still, important results are depicted: (1) possible tsunami samples of both cores differ from the background sediments by higher abundance of quartz, mica, and terrestrial biogenic components and mainly in more abundant mollusk fragments (core M106). Also, in both cores, tsunami samples have lower amounts of foraminifera. (2) Core M107, representative of the inner shelf, shows a more A general difference in SC between core M107 and M106 is apparent: (1) In the shallower shelf core M107 higher percentages of terrigenous components (ca. 50%) are present compared to the deeper shelf core M106 (ca. 30%). (2) In core M106, percentages of benthic foraminifera are higher, and percentages of mollusks are lower compared to core M107. The terrigenous part of core M107 is made up of mainly quartz and mica, while opaques, other terrigenous, and aggregates are low in percentages (<8%) (Figure 9b). The biogenic components are mainly other biogenic, mollusk, and benthic foraminifera. Terrestrial biogenic components have very low percentages throughout the studied core section. Although core M106 shows general fluctuations in the SC there is no general trend in the studied section (Figure 9a). The terrigenous part consists mainly of quartz, mica, and aggregates, while opaques and other terrigenous components are underrepresented. A peak of aggregates at depth 73 cm in core M106 is accompanied by lower densities as depicted by the CT-data (Figure 7). The biogenic part is dominated by other biogenic components and benthic foraminifera (up to 36%). Mollusk and terrestrial biogenic components are low in this core with maximal percentages of ca 5%. In the possible tsunami layer, no remarkable peak of terrigenous components can be observed in either cores, that could serve as a backwash tracer. However, higher percentages of mollusk shells/fragments (only in core M106) and terrestrial biogenic components are present in the high-energy event layer with peak percentages within the total studied section, indicative of sediment transportation from beach or nearshore areas to the shelf.
For each core, a separate PCA of the SC was conducted, and results are visualized in Figure 10. A solution of 4 principal components was chosen, and in both cores selected principal components explained >77% of the total variance in the data set. For core M107, selected component assemblages are PC1_107 and PC2_107, explaining 49.9% of the total variance and PC1_106 and PC4_106 for core M106 that explain 50% of the total variance. Both cores show a distinction between the chronological groups by their SC, but differences of the background sedimentation and the possible tsunami units are small. Still, important results are depicted: (1) possible tsunami samples of both cores differ from the background sediments by higher abundance of quartz, mica, and terrestrial biogenic components and mainly in more abundant mollusk fragments (core M106). Also, in both cores, tsunami samples have lower amounts of foraminifera. (2) Core M107, representative of the inner shelf, shows a more distinct differentiation of the tsunami samples from the background sedimentation, while core M106, which is representative of the outer shelf, shows a weaker tsunami signature.
Minerals 2019, 9, x FOR PEER REVIEW 18 of 30 distinct differentiation of the tsunami samples from the background sedimentation, while core M106, which is representative of the outer shelf, shows a weaker tsunami signature.

Microtextural Analysis of Quartz Grains
In general, the quartz grain samples show minor differences between them. The predominant microtexture is dissolution with median values of 3 (25-50% of grain surface) up to the whole grain surfaces, in some cases ( Figure 11) while percussion marks are the second most abundant followed by fresh surfaces. Adhering particles occupy minor percentages of the grain surfaces. Roundness (R) and roughness factor (RF) show a small variation between samples and show similar trends. Grouping the samples in their respective chronological regime shows small and no significant differentiation between tsunami and non-tsunami samples. However, a tendency of more percussion

Microtextural Analysis of Quartz Grains
In general, the quartz grain samples show minor differences between them. The predominant microtexture is dissolution with median values of 3 (25-50% of grain surface) up to the whole grain surfaces, in some cases ( Figure 11) while percussion marks are the second most abundant followed by fresh surfaces. Adhering particles occupy minor percentages of the grain surfaces. Roundness (R) and roughness factor (RF) show a small variation between samples and show similar trends. Grouping the samples in their respective chronological regime shows small and no significant differentiation between tsunami and non-tsunami samples. However, a tendency of more percussion marks on possible tsunami samples can be seen and adhering particles are significantly more abundant on post tsunami grains.

Discussion
A multi-proxy approach was applied to investigate sedimentary signatures of the 1755 CE tsunami in the southern Portuguese continental outer shelf. The presumed 1755 CE tsunami layer in the southern Portuguese continental shelf off Faro is identifiable by the combination of several

Discussion
A multi-proxy approach was applied to investigate sedimentary signatures of the 1755 CE tsunami in the southern Portuguese continental outer shelf. The presumed 1755 CE tsunami layer in the southern Portuguese continental shelf off Faro is identifiable by the combination of several indicators frequently discussed in the literature. Before the interpretation of the presumed tsunami signatures, the correspondence of the studied sediments to the 1755 CE tsunami is discussed.

Correspondence to the 1755 CE Tsunami
Storms are capable of creating similar deposits as tsunamis and their discrimination is problematic [11,29,87]. To be certain that storms did not affect the tsunami deposits, the deposits have to be either in storm sheltered areas or beneath the storm wave base. After Weiss and Bahlburg [12] the water depth (d) where a storm wave affects the bottom sediments (d b ) is: where g is the gravitational acceleration (g ≈ 9.81 m·s −2 ) and A the amplitude (A = wave height (H) · 2 −1 ), assuming exactly the boundary between deep-water (0.5 · wavelength (L) < d) and intermediate-water (0.5 L > d > 0.05 L). Core M107 has the lower water depth of the studied cores with 58 m. By assuming d b = 58 m, the amplitude of the storm wave must be at least 7 m (or H = 14 m) to be capable of reworking the tsunami sediments. The storm of 1941 CE, assumed to be the biggest storm of the century [88], had (modeled) wave heights of max. 14 m and periods of 14 s [89]. Considering a SW direction of the storm, similar wave heights could have reached the present study area, but an erosion of the tsunami layers in water depths > 58 m can be precluded. Fluvial floods are other possible events leading to a terrestrial input to the offshore areas. Similar to tsunami backwash flows, flash floods can generate hyperpycnal flows because of their high suspension load. However, there are different sedimentological characteristics between tsunami and flood deposits. The latter show generally better sorting and are composed of less coarse material than tsunami backwash deposits, because of the higher energy involved in tsunami backwash events [90]. Regarding the lower energy in flood events, their sedimentary/hydrodynamic equilibrium is reached in closer proximity to the coast and thus unlikely to have a pronounced depositional signature in the outer shelf. Despite this, in the Algarve region, no extreme fluvial flood is reported in the study area, although maximal historical peak discharges of~11,000 m 3 ·s −1 in 1876 CE and 1603 CE were documented for the next biggest and closest river (Guadiana River) (80-100km) to the study area [91].
Recorded historical tsunamis have the advantage to allow chronological correlation between the year of the tsunami event and the potential event layer as a first basis. The established age models are based on a small number of 14 C-ages and 210 Pb ages in the topmost sections and some radiocarbon ages are taken from probable reworked material (Figure 4), which can lead to incorrect age estimations [92]. However, in the age model function of "rbacon", outliers are addressed by considering a student-t distribution instead of assuming an error of the radiocarbon age [93], minimizing the effect of outliers [76]. Thus, by using Bayesian age-depth models, the most realistic precision and robustness, compensating the small number of dated levels, can be expected compared to such classical methods as linear interpolation [94]. Results of the established age-depth relations suggest a highly probable correlation of the 1755 CE Lisbon event with the presumed tsunami deposit under study. The deposit in core M106 also likely correlates, on the basis of similar core depths (~75 cm to 40 cm), to possible 1755 tsunami deposits in proximity (~400 m) to core M106, as described by Quintela et al. [22]. Dating tsunami deposits can be challenging even when clear limits of the tsunami layer are established. The reason is the erosional action of tsunami events, incorporating material of different ages [95]. Thus, age errors of offshore tsunami deposits are expected to be relatively high, because confining the event unit is additionally complicated (see Section 4.2). Also, many assumptions have to be accepted, such as the marine reservoir effect, calibration curves, and constant rate of 210 Pb supply to the sediment (etc.), increasing the age error [82]. Concerning tsunami ages, it is suggested to have a precision within a few hundred years or less to allow the implication of tsunami frequency [82]. Here, errors indicated by the age model are 246 yrs and 468 yrs for the tsunami units of core M106 and M107, respectively ( Figure 4). The age model of core M107 should be thus supplemented with more radiocarbon ages to increase its precision. High age uncertainties were also reported in other offshore paleo-tsunami studies [3,96], stressing that age control of tsunami deposit is difficult and needs to be conducted carefully. Another issue of the age uncertainty is the distinguishing of the transatlantic tsunami in 1761 CE [97] from the 1755 CE event. The source area of the 1761 CE tsunami is proposed to be in the Southwest Iberian Margin similar to the 1755 CE tsunami [97]. Despite the effects of the earthquake in mainland Portugal, the related tsunami had about zero impact in respect to human and economic losses [47]. It is possible, that the 1761 CE event created tsunami related signatures in the studied sedimentary records. However, regarding the much lower severity of the 1761 CE tsunami compared to the 1755 CE event [47,98], it seems more unlikely to have signatures of the 1761 CE tsunami in the sedimentary record.

Identifying Outer Shelf Tsunami Deposits
Identifying tsunami event layers in outer shelf sedimentary records is not a straightforward task and several authors suggest a multidisciplinary approach [4] since tsunami deposits are notably site-specific and share only a few common characteristics [29]. Moreover, subtle sedimentological changes are expected instead of unequivocal indicators for the upper and lower boundary of an offshore tsunami deposit. Grain size increases in offshore tsunami deposits, as in core M107 (Figures 5b and 8b), is a common global feature [4,14,[22][23][24]96] and related to the higher energy involved during the tsunami event when compared to the background sedimentation. The absence of grain size variations [3,17], and decreased grain size [16] are also reported, especially in water depths > 50 m. In core M106, grain size variations are minor (Figures 5a and 8a). However, an increase in grain size was detected in proximity to core M106 for the correlative interpreted tsunami layer by Quintela et al. [22]. Considering a channelized backwash flow [28], grain size variation in offshore tsunami deposits might greatly depend on the local bathymetry. Some offshore tsunami deposits show stacked fining up sequences interpreted as different backwash phases [23,31,90], which is not well recorded in the here studied tsunami layer. However, both studied tsunami units show a fining upward sequence (Figure 5a,b) that might reflect the decreasing energy evolution of the backwash flow. SC analysis revealed subtle increases in quartz, mica, and terrestrial biogenic components (Figures 8 and 9). In the study area, quartz and mica (here Illite and Kaolinite) are abundant in onshore sediments/rocks [56] and the here used proxy, XRF elemental counts of Al + Si, increases in the tsunami deposit ( Figure 5). High Ti/Ca-ratios were used as indicator for increased terrestrial input in offshore tsunami sediments [3,4,15,19,30,90]. Although, in core M106 the highest value of Ti/Ca is observed in the tsunami layer ( Figure 5), this proxy shows no pronounced occurrence in the 1755 CE tsunami layer in the southern Portuguese shelf (present study) and the western Portuguese shelf [17]. In the southwestern Portuguese margin, a turbidite deposit related to the 1755 CE event shows, however, clear peaks in Ti/Ca (depicted as Ca/Ti) [96]. The weak peak of Ti/Ca can be linked to the different geology in the southwestern part of the Algarve compared to south coast in proximity to the studied deposits. Numerical modeling of the tsunami wave can aid in the detection of approximate sediment source areas. In the case of the deposit of this study, possible sediment provenance might be the coastal area between Faro and Quarteira ( Figure 1) [68,99]. This highlights the importance of the local geology from the event sediment sources, and proxies must be adapted accordingly. The accumulation of allochthonous mollusk fragments in offshore tsunami deposits is another frequently observed feature [4,8,18,25,[100][101][102], dragged from shallower water depths by backwash flows [8]. Abundant mollusk fragments were observed both, in the present study (Figures 8 and 9) and in the 1755 CE tsunami deposit southwest of the Tagus estuary (Portuguese west coast) from similar water depths [17,100]. The accumulation of mollusk fragments is further validated by a general high carbonate content in the studied tsunami layer (Figures 5 and 8) [17]. The possible provenance could be the inner shelf, where mollusk fragments are more abundant [103]. Besides the abundant occurrence of mollusk fragments in the coastal areas, their planar shape might increase their buoyancy, and thus their sedimentary equilibrium is reached farther offshore, possibly in the outer shelf. Similar behavior of planar shaped particles was proposed for mica, suggesting a higher abundance of those components in upper tsunami units [104]. In core M107, an increase of mica in the upper tsunami unit can be observed, which is subsequently decreasing again (Figure 9), possibly indicating the decrease in energy. MS increases accompanied by higher Fe and organic matter contents was interpreted as an input of fine material, consisting of fine-grained magnetic minerals [17]. The same can be found in the upper tsunami unit in core M106 (Figures 5a and 8a), suggesting a settling of finer material in the upper tsunami unit of sediment-rich plumes generated by the tsunami event. Also, in the present-day sedimentation, studies show that only~14% of fine land-derived material is deposited in the shelf with the tendency of fine-material transport beyond the outer shelf [57]. This could further explain (1) the minor grain size variation within the tsunami layer in core M106, and (2) the lack of the pronounced fine-material input in the upper tsunami unit of core M107 ( Figure 5) because of the more landward location. The MT analysis did not detect a clear tsunami-related signature, although percussion marks are slightly increased in tsunami samples (Figure 11), which coincide with results from onshore deposit of the 1755 CE tsunami [34] and in onshore tsunami samples in general [33,105]. An increase in percussion marks was also observed in laboratory experiments, where high sediment concentration and supercritical flow conditions carved percussion marks onto the grain surfaces [36]. Another signature is the wider range of observed microtextural types (i.e., adhering particles, dissolution, fresh surface, percussion marks, and angularity) on tsunami quartz grain surfaces which was not found here. Besides the identification of high energy conditions, MT analysis could be useful for provenance studies of tsunami sediments when the geological setting of the study area is well constrained [37] and samples of the potential source material are available for the analysis.
Detecting the lower and upper bounds of offshore paleo-tsunamis in the sedimentary record is difficult, especially for the upper limit, since biological and physical post-depositional alterations affect sedimentary structures and tsunami tracers within the tsunami layer [3,9,31,106]. The lower boundary is mostly indicated by an erosional feature caused by the incoming tsunami wave or by the turbiditic backwash current. In studies of recent and paleo offshore tsunami deposits in the outer shelf, erosional surfaces were detected and used for setting the lower boundary [3,20], while other studies lack this identification [2,22]. Smedile et al. [2] and Quintela et al. [22] explain this by the farther distance to the shore or a reworking of the sediments, rather than a strong erosional action. In the present study, the presumed tsunami layers lack a sharp basal contact at the onset of the tsunami unit. However, the CT-image of core M106 reveals a layer of less dense sediments in proximity to the base of the tsunami unit (Figure 7), coupled with a small peak of magnetic susceptibility and a peak in aggregates (Figures 5a, 6 and 9). The abrupt change in radio-density at ca. 70 cm core depth (Figure 7a) and the internal sedimentary structures (i.e., of lamination or ripples) below (Figure 7b) can be related to the action of erosion. Speculatively, this can be interpreted as a reworking of the sediments by the land-ward propagating tsunami wave. Amore significant piece of evidence is the low numbers of total foraminifera tests at the onset of the tsunami unit in core M106 (Figure 9). Interestingly, the same trend was observed in tsunami deposits close to the shelf break in Augusta Bay (Eastern Sicily, Italy) [3] and Sendai Bay (Japan) [20]. Smedile et al. [3] interpret the low foraminifera numbers as a major input of fine sediments that can be assumed in the present study as well, based on the increase in Fe-content, MoR, and magnetic susceptibility (Figures 5 and 6). In core M107, the base of the tsunami layer is furthermore identified through increased sand-sized particles and carbonate content (Figure 5b). The upper limits of the tsunami units are set by diminishing of the tsunami proxies and controlled with compatible age estimations, although the set limits remain speculative. Despite under such uncertainties, it is possible to verify that χ LF values within possible tsunami layer of both cores mostly agree with an interruption of the overall increasing trend to a steeper increase, relative to the underlying sediment ( Figure 6). This means that the general trend of χ. χLF values is marked by an increase at the uppermost 150 cm of core M106 and uppermost 30 cm of core M107 that are interrupted within the suggested tsunami layers for core M106 and the post-tsunami layer for core M107. The general trend is expected due to early sedimentary diagenesis [107] but the origin of such interruption is not clear. Reworking of the sedimentary column and/or a change of the background sedimentary regime could be a plausible explanation.
As seen above, the limits of the tsunami deposits are questionable and clear tsunami-related sedimentary patterns are not evident. However, the combination of many small pieces of evidence determined by different proxies, namely a subtle terrestrial fingerprint, grain size variation, more abundant mollusk fragments, a decrease in total foraminifera abundance, and fine material input backed up by matching age correlations, are consistent with a deposition by the 1755 CE tsunami of the two units under study. Here, multivariate statistical analysis is essential to detect small changes in the outer shelf sedimentary record and to aid interpretation of the significance in proxy-variations between the tsunami layer and the marine background sedimentation. Moreover, the small increase of percussion marks of tsunami samples ( Figure 11) could point to higher energy conditions during the tsunami event. Differences between the cores in the two tsunami units are present and can be ascribed to the differing core locations, and thus, to probable differing sedimentation processes. Based on the PCA and LDA results of both cores (Figure 8a,b,d), it might be concluded that it is easier to filter out differences between tsunami and background sediments in sand-rich deposits compared to finer deposits. However, the higher sedimentation rates of core M106 and the seawards decreasing energy of tsunami events [108] might be another reason for greater differences.
It is known that no remarkable terrestrial signatures are expected in the outer shelf environment [5,31]. Moreover, wave flume experiments showed, that offshore tsunami deposits are onlỹ 30% composed of onshore material [108]. Hence, wave motions and return flows might be a major process incorporating autochthonous material and creating tsunami-related sedimentary structures of offshore tsunami deposition [20,108]. This might additionally veil the tsunami signal in the sedimentary record and explains the subtle terrestrial fingerprint. Thus, tsunami-related signatures in the outer shelf sedimentary record are small discrepancies between periods of normal marine sedimentation. Nevertheless, the poor signatures of the severe 1755 CE tsunami event in offshore sedimentary records (present study and in others [17,22,109]) and the fact that, among several well-documented historical tsunamis, only the 1755 CE tsunami left geological evidence onshore [47], is contradicting. Recent (preliminary) offshore studies in the southwestern part of the south Portuguese shelf have also encountered relatively small signatures of the 1755 CE tsunami [109]. However, other possible older events show a very clear signature such as a sandy layer interbedded in silty normal marine background sediments with well supporting XRF-data and biomarker analysis [109]. The small signature of the 1755 CE tsunami could be explained by a change in sedimentation dynamics in general observed by PCA results including both cores (Figure 8c,d) (see also Figure 2), the mass normalized magnetic susceptibility variation ( Figure 6) and other studies of the southern Portuguese shelf sedimentary records [110][111][112]. Several hypotheses are discussed for this change in sedimentation, i.e., of a decreasing trend in grain size and higher sedimentation rates: (1) the change of coastal morphology, such as the opening/closing of the Ria Formosa barrier island system [113], (2) the intensification of fine-grained loaded river discharge associated with the wetter climatic condition since the Little Ice Age [114], and (3) the increased settlement density induced land-use change, that probably resulted in an increased input of terrigenous material into the shelf [111]. However, more studies regarding the provenance of the material need to be conducted to interpret this change and the contradicting small signature of the 1755 CE tsunami in the southern Portuguese continental shelf.
In summary, we presented a detailed analysis of two deposits related to the 1755 CE tsunami. SC analysis reinforced by XRF analysis seems to be an appropriate approach to detect tsunami deposits in similar environments. Improvements in sand grain analysis, such as automated recognition and counting [115], help to extend the analysis in the sedimentary record to establish possible provenances of the tsunami sediments. Also, changes in the background sedimentation could be established and linked to paleoenvironmental changes. However, the small terrestrial fingerprint of the tsunami deposits calls for methodologies that reveal sedimentary structures, such as CT on X-ray core images, produced by wave motions and backwash flows of the tsunami [1,116]. Although it is possible to detect tsunami deposits in outer shelf environments, recent works suggest the inner shelf and more protected bays as a more suitable study area [5,15]. Concerning the preservation potential of the event deposits, inner shelf environments are at risk of erosion by storm waves (see Section 4.1), while bays might be sheltered for some tsunamis. Open shelf settings have probably less developed tsunami deposits but possibly cover more tsunami events. Hence, investigating offshore tsunami sediments is important to complete the dataset of tsunami events and enable accurate estimations of recurrence intervals for tsunami hazard assessments [117]. The results of this study will aid other studies on outer shelf tsunami deposits and the findings of the detailed analysis of the presumed 1755 CE layer can be used on other possible tsunami layers already recognized in the same cores (see Section 1.2).

Conclusions
We presented a detailed analysis of the 1755 CE tsunami sediments in outer shelf settings applying a multi-proxy approach. From the results of our data, we can extract the following conclusions:

•
In the present study, there is no pronounced single-proxy evidence for tsunami sedimentation. Yet, taking the sum of multiple small pieces of evidence reveals evidence of a tsunami-related deposition. Relative to non-tsunami samples, presumed tsunami samples have a subtle terrestrial fingerprint (increases of Al + Si ratio, quartz, terrestrial biogenic components), grain size variation, more abundant mollusk fragments, a decrease in total foraminifera abundance, and fine material input with deposit age ranges that overlap the 1755 CE event. The subtle terrestrial fingerprint of the tsunami sediments found here, as well as results in the literature, implies that the usage of methodologies that can reveal sedimentary structures related to the action of the tsunami wave and the backwash flow might have a higher potential. Moreover, outer shelf tsunami deposits depend greatly on the local geology and bathymetry, and differences between the two tsunami units were encountered related to the differing core location at 58 m and 91 m bsl. The site-specific tsunami deposit characteristics call for adapting the methodology to identify these deposits in similar environments.

•
Microtextural analysis of quartz grain surfaces was performed for the first time in outer shelf tsunami samples. Although tsunami samples showed a minor increase in percussion marks, this increase was insignificant, and no other possible tsunami signature was detected, suggesting that microtextural analysis is not useful for identifying tsunami deposits in shelf settings, unlike for onshore high-energy deposits. • A possible change in the general sedimentation following the tsunami event might have influenced the signature of the 1755 CE tsunami in the outer shelf sedimentary record. The contradicting result, that the 1755 CE tsunami is so far the only historical tsunami event that left geological evidence onshore and its subtle signature on the shelf, needs further investigation.

•
The methods described here to study outer shelf tsunami deposits (CT scans, SC reinforced by XRF, magnetic parameters, and multivariate statistical analysis) can be applied for other potential tsunami deposits in the cores under study to ultimately improve the dataset of past tsunamis in Portugal, which lacks geological data and the extension to pre-historic times.