Washover deposits related to tsunami and storm surge along the north coast of the Shimokita Peninsula in northern Japan

A decade after the 2011 Tohoku-oki earthquake (Mw 9.0), geological surveys were conducted at multiple sites along the Pacific Coast of the tsunami-inundated Tohoku region in Japan, providing thousands of years of tsunami history. However, the challenges of correlation between historical records and geological tsunami deposits and identifying sources of historical and paleotsunamis have newly surfaced. Particularly the simultaneity and source of the 1611 Keicho tsunami in the Tohoku region and the seventeenth-century tsunami in the Hokkaido region are problematic. To solve such major issues, we conducted a tsunami-deposit survey at Sekinehama on the north coast of Shimokita Peninsula, near the junction of the Japan and Kuril trenches. We performed nondestructive analyses (X-ray computed tomography and micro-X-ray-fluorescence core scanning), grain-size analysis, tephra analysis, and radiocarbon dating of sediments from two coastal outcrops and inland drill cores. We identified five tsunami deposits (TD1–TD5) during the last 6 kyr and correlated them at a 200–400 m distance from the coast. They also correlate with previously identified tsunami deposits around the Shimokita Peninsula. From our study on tsunami deposits, we found other washover deposits in the coastal outcrops that are not represented in the inland cores. These indicate minor washover events related to small tsunamis and infrequent storm surges. The modeled age of the latest tsunami deposit is 500–300 cal yr BP (1450–1650 cal CE). This either correlates with two known tsunamis (the 1611 Keicho tsunami and another seventeenth-century tsunami) or is a previously unknown tsunami that occurred in the fifteenth–seventeenth centuries. If the latest tsunami deposit is to be accurately correlated with tsunami deposits previously identified within a 50-km distance from the study site, we need to consider an unknown fifteenth-century tsunami. Our investigation yields insights regarding the tsunami source in the vicinity of the junction of the Japan and Kuril trenches.


Introduction
In the 10 years since the 2011 Tohoku-oki earthquake (Mw 9.0), tsunami-deposit research has shown dramatic progress not only in Japan but also across the world (Goto et al. 2014(Goto et al. , 2021Costa and Andrade 2020;Sawai 2020;Sugawara 2021). In addition to the discovery of new tsunami deposits, research on numerical modeling of tsunami and sediment transport (Sugawara 2021) and on tsunami size estimation (Namegaya and Satake 2014; Ishimura and Yamada 2019; Naruse and Abe 2017; Mitra et al. 2020Mitra et al. , 2021 has been conducted. In particular, geological surveys were performed at multiple sites along the Pacific Coast of the Tohoku region where the 2011 tsunami inundated provided thousands of years of tsunami history (e.g., Goto et al. 2015;Miyauchi 2015, 2017;Takada et al. 2016;Inoue et al. 2017). Although long-term tsunami-deposit history was obtained at multiple sites, correlation between tsunami deposits from site to site and between historical records and geological surveys and estimation of the probable sources of historical and paleotsunamis have been challenging (Takashimizu 2014;Goto et al. 2021;Ishizawa et al. 2022).
Particularly for the fifteenth-seventeenth-century tsunami deposits whose ages cannot be sufficiently constrained by radiocarbon dating owing to wiggles in the calibration curve, the identification of the 1611 Keicho (e.g., Ebina and Imai 2014) and the 1454 Kyotoku (Namegaya and Yata 2014) tsunami deposits and the simultaneity and source of the 1611 Keicho tsunami in the Tohoku region and the seventeenth-century tsunami in the Hokkaido region remain major foci for investigation. Ishizawa et al. (2022) tackled the long-standing challenge of distinguishing the 1611 Keicho and 1454 Kyotoku tsunami deposits by applying high-resolution radiocarbon dating and concluded that only the 1611 Keicho tsunami deposit is distributed along the northern Sanriku Coast (Fig. 1a). Takashimizu (2014) discussed the simultaneity and wave source of the seventeenth-century tsunami in Hokkaido and the 1611 Keicho tsunami in Tohoku and the distribution and age of tsunami deposits in western Hokkaido, such as along the Hidaka Coast (Fig. 1a), to determine how widely distributed the seventeenth-century tsunami deposits are in Hokkaido. The northern part of the Pacific Coast of the Tohoku region, such as the Shimokita Peninsula, has high coastal dunes that limit appropriate sites for tsunami-deposit surveys in the north of Hachinohe (Fig. 1a, b) (Minoura et al. 2013;Tanigawa et al. 2014).
In our study, we conducted a tsunami-deposit survey at Sekinehama on the north coast of the Shimokita Peninsula in the northernmost part of the Pacific Coast of the Tohoku region ( Fig. 1) to provide continuous and long-term tsunami history and fill in blank areas of tsunami-deposit information near the junction of the Japan and Kuril trenches. The multiple event deposits that we identified in the Holocene sediments and dated as from 6 ka or later are considered to be transported by tsunamis and storm surges. This study is in coordination with Sato et al. (in press in this special issue), who investigated the wave source of the latest tsunami deposit discussed in our study.

Settings
The geologic units along the north coast of the Shimokita Peninsula (Fig. 1b, c) consist of the Pliocene Sunagomata Formation, the Pleistocene Tanabu Formation, and superposed marine-terrace deposits of multiple ages within the Pleistocene (Tsushima and Takizawa 1977). Marine Isotope Stage 5e marine terraces are distributed across the study area ( Fig. 1c) (Koike and Machida 2001), and Holocene sediments are distributed in the lowlands and incised valleys. Tephra and drift pumice layers were already identified in the Holocene sediments at two outcrops along the coast (SKN-E and SKN-W in Fig. 1c). In SKN-W, Kobayashi (2019) found fallout ash of B-Tm tephra (946 CE; Oppenheimer et al. 2017) from Changbaishan. In SKN-E, fallout ash of To-Cu tephra (5986-5899 cal yr BP; Mclean et al. 2018) from Towada volcano was identified by Ishimura and Hiramine (2020), and drift pumice of To-Cu and U-2 tephras (5681-5619 cal yr BP; Mclean et al., 2018) from Ulleung Island, Korea, was identified by Hiramine et al. (in press).
Historical and observational records of tsunamis hitting this region (Table 1) show that they are of two types: seismic and volcanic tsunamis. Most of the seismic tsunamis in Hokkaido and Tohoku regions have their sources along the Pacific Coast. Records show that past tsunamis were < 3 m high and may not have inundated inland in the study area because the topographic height behind the present-day beach is 4-5 m (Fig. 1c). The recorded tsunami height of the 2011 Tohoku-oki earthquake 200 m west of the study site was 2.2-2.7 m (Mori et al. 2012), suggesting the absence of a large inundation of inland area in Sekinehama. Volcanic tsunamis are thought to be generated by sector collapses of volcanoes due to volcanic activity at Komagatake and Oshima-Oshima volcanoes in Hokkaido (Fig. 1a) (Nishimura and Miyaji 1998;Satake 2007). The tsunamis' inundation heights were documented and estimated near their source areas (Nishimura and Miyaji 1998;Satake 2007), but no such information is of record for on the Shimokita Peninsula.
Based on prior research or observational records, the following summarizes information known or proposed about widely recognized tsunami deposits along the Pacific Coast in the Hokkaido and Tohoku regions. The 1611 Keicho tsunami deposits have been identified from the Sendai Plain to the northern Sanriku Coast (e.g., Sawai et al. 2012;Ishimura and Miyauchi 2015;Takada et al. 2016;Ishizawa et al. 2022). The highest tsunami run-up height along the Sanriku Coast was estimated to be more than 25 m (Ebina and Imai 2014;Ishimura and Ebina 2021). Recently, using high-resolution radiocarbon dating, Ishizawa et al. (2022) confirmed that 1611 Keicho tsunami deposits are distributed along the northern Sanriku Coast (Fig. 1a). The seventeenth-century tsunami deposits are widely recognized in eastern Hokkaido (e.g., Nanayama et al. 2003;Sawai et al. 2009;Ishizawa et al. 2017). Tsunami source models for the seventeenth-century tsunami were proposed by Nanayama et al. (2003), Satake et al. (2008), and Ioki and Tanioka (2016). The latest model proposed by Ioki and Tanioka (2016) setting the fault off eastern Hokkaido calculated the total seismic moment magnitude as Mw 8.8. The simultaneity of this Page 3 of 27 Ishimura et al. Progress in Earth and Planetary Science (2022) Ishimura et al. Progress in Earth and Planetary Science (2022) 9:69 tsunami with the 1611 Keicho tsunami has been debated but is not yet fully understood (e.g., Tetsuka et al. 2020). From research in the vicinity of the study site, two studies, Minoura et al. (2013) and Tanigawa et al. (2014), identified tsunami deposits along the east coast of the Shimokita Peninsula (Fig. 1b). Minoura et al. (2013) identified one tsunami deposit as having occurred between the late thirteenth and fifteenth centuries. Tanigawa et al. (2014) identified five sandy event deposits interbedded within marsh deposits during the last 6000 years, three of which were attributed to inundation from the coast on the basis of diatom assemblages. The latest one is estimated to have occurred 470-180 cal yr BP. Along the Sanriku Coast, many tsunami-deposit surveys were conducted after the 2011 event. They show the relatively high frequency of paleotsunamis from the north to the central parts of this coast (Fig. 1), where their average recurrence interval is ~ 300-500 years (Ishimura 2017;Ishimura and Yamada. 2021;Ishizawa et al. 2022;Velasco-Reyes et al. 2022). On the Hokkaido side, many drilling surveys were also conducted. In eastern Hokkaido, a short average recurrence interval (close to 400 years) was estimated from long-term tsunami history based on the tsunami deposits (e.g., Sawai et al. 2009), whereas prior tsunamideposit researchers found a longer average recurrence interval to the west of the Hidaka Coast than eastern Hokkaido (Takashimizu et al. 2007;Takashimizu 2014;Nakanishi et al., 2020Nakanishi et al., , 2022. To the west of the Hidaka Coast (Fig. 1a), only one tsunami deposit was identified since ca. 3 ka and was correlated with the seventeenthcentury tsunami (Takashimizu et al. 2007). Additionally, the lateral correlation and sources of the tsunami deposits formed during the last 500 years in the west of the Hidaka Coast (Fig. 1a) are not well understood (Takashimizu 2014). Recently, Kase et al. (2022) reported tsunami deposits around Hakodate on the opposite side of the Tsugaru Strait from the study area (Fig. 1a). They drilled at four sites and identified tsunami deposits formed since ~ 3000 cal yr BP.

Methods
We conducted a field survey to log two outcrops (SKN-E and SKN-W) along the Sekinehama Coast and cores drilled along the incised valley, which is a subsidiary of the Bitsuke River (Fig. 1c). An additional coring survey was conducted in the south of the SKN-W outcrop to observe the younger sediments. A handy-Geoslicer (GS) (Takada et al. 2002) was used for the coring survey, by which we obtained a total of 11 GS cores.
For nondestructive analyses, we extracted block samples from outcrops and GS cores with a 50-cm-long, 4-cm-wide, and 3-cm-thick case. These block samples were used for X-ray computed tomography (XCT) and high-resolution micro-X-ray-fluorescence (µXRF) scanning at the Center for Advanced Marine Core Research, Kochi University. The XCT images were taken by using an Aquilion PRIME/Focus Edition scanner (Canon Medical Systems Corporation). For this analysis, a sample was imaged according the X-ray absorption so that high-density materials could be identified by their relatively light color (Falvard and Paris 2017). A vertical slice of the block sample was taken every 0.5 mm, and these vertical slices were then processed using Horos software (a free and open-source code software  Ishimura et al. Progress in Earth and Planetary Science (2022) 9:69 under the GNU Lesser General Public license, Version 3 (LGPL-3.0)) to build a 3D image of the block sample. Using the ITRAX µXRF core scanner (Croudace et al. 2006), the block sample was measured at 1-mm intervals. The scanner, which irradiated the block sample surface by applying a molybdenum X-ray source, operated at 60 kV and 50 mA and allowed us to obtain data on 30 elements. We used 12 element intensities (Si, S, K, Ca, Ti, Mn, Fe, Ni, Rb, Sr, Y, and Zr) after considering previous studies (Chagué -Goff et al. 2016;Martin-Puertas et al. 2017;Ishizawa et al. 2019). Because the elemental data were influenced by the sample's condition (Löwemark et al., 2011;Chagué-Goff et al. 2016), we normalized each elemental datum over counts per second for comparison between cores (e.g., Davies et al. 2015;Chagué-Goff et al. 2016;Ishizawa et al. 2019). We determined major element compositions of volcanic glass shards for correlating known tephra. We sampled sediments from GS cores for tephra analysis, while considering the already-identified B-Tm and To-Cu tephras in the coastal outcrops. For the cores with no visible tephra, we conducted continuous sampling and calculated the content of volcanic glass shards. First, we used an ultrasonic generator to rinse the samples in water and then sieved them with water through a 62-μm nylon mesh. Next, we dried the samples at 50 °C and sieved the dried samples through a 120-μm nylon mesh. Then, using a polarizing microscope, we counted volcanic glass shards in the 62-120-μm fraction. Finally, we determined the major element composition of the volcanic glass shards in that size fraction by energy-dispersive X-ray spectrometry (EDAX Genesis APEX2) and scanning electron microscopy (JEOL JSM-6390) following the method and analysis conditions described by Suzuki et al. (2014). Volcanic glass shards of the AT tephra sampled at Chigaki, Toyama Prefecture (Machida and Arai 2013), were used as a working standard to check data reproducibility and instrument stability.
For grain-size analysis, we collected all event-deposit samples from the outcrops and cores, and present-day sediment samples surrounding the study site (Fig. 1c). As for the present beach sediments, we confirmed that there has been no artificial beach nourishment by the interview to the city office. First, each sample was pretreated with 3% H 2 O 2 to eliminate organic matter and then rinsed with water using an ultrasonic generator. Then, we sieved the samples with water through a 62-μm nylon mesh and then used a 2-mm sieve. Finally, we measured the weight of sand and gravel fractions to calculate mud, sand, and gravel contents. Grain-size analysis of sand fractions was performed using Morphologi G3 (Malvern Instruments Ltd.) according to Ishizawa et al. (2019), but in our study, we used the 5X magnification setting.
Radiocarbon dating using an accelerator mass spectrometer was performed at the Institute of Accelerator Analysis Ltd., Japan. We obtained 42 radiocarbon dates: 32 for samples from GS cores and 10 for samples from outcrops. Radiocarbon dates were calibrated based on the IntCal20 calibration curve (Reimer et al. 2020) using the OxCal 4.4 program (Bronk Ramsey 2009). We used the sequence model (Bronk Ramsey 2008) for calculating depositional ages of tsunami deposits (Additional file 1).

SKN-E outcrop
The SKN-E outcrop (Fig. 2) is ~ 300 m east of the mouth of the Bitsuke River (Fig. 1c) and is on a sea cliff that is 2.5 m high and 10 m wide. The sediments exposed in this outcrop apparently belong to two different sedimentary facies. The lower part consists of blue-gray silt (possibly marine sediments), and the upper part includes marsh deposits consisting of organic silt and peat. The uppermost part is artificially modified sediments that contain concrete blocks and fishnets. In this study, we described and sampled sandy event deposits in the marsh sediments.
We recognized 12 sandy event deposits in the SKN-E marsh sediments (SKN-E-E1 to SKN-E-E12, from youngest to oldest) ( Fig. 2; Table 2). Except for SKN-E-E7, they show good horizontal continuity. Most consist of well-sorted fine to medium sand. SKN-E-E5, SKN-E-E6, SKN-E-E8, and SKN-E-E9 contain coarser grains such as granules. SKN-E-E5 and SKN-E-E6 have an erosive contact at the base and contain rip-up clasts in the lower part. SKN-E-E8 also contains rip-up clasts. Some of the sandy event deposits show normal grading. The bottom shape of SKN-E-E7 is convex downward, suggesting a channel, and the deposit has parallel laminations.
From the grain-size analysis (Fig. 3), SKN-E-E1 to SKN-E-E4 show grain-size distributions concentrated in the range 1-2 phi. Meanwhile, SKN-E-E5 to SKN-E-E10 show modes of 1-2 phi, but exhibit a wider grainsize distribution. In the event deposits below SKN-E-E11, the mode changes from 1-2 phi to 1.5-2.5 phi with finer grains. This is apparent from the mud content (Fig. 3). The gravel contents also support the above descriptions indicating that SKN-E-E5, SKN-E-E6, and SKN-E-E8 include coarser grains.
Three tephra-related layers in the marsh sediments were already identified (Fig. 2) and described in detail in Ishimura and Hiramine (2020) and Hiramine et al. (in press). To-Cu (5986-5899 cal yr BP) tephra (finegrained volcanic ash) is interbedded between SKN-E-E9 and SKN-E-E10 (Ishimura and Hiramine 2020), and drift pumice clasts that originated from To-Cu tephra Page 6 of 27 Ishimura et al. Progress in Earth and Planetary Science (2022) 9:69 are concentrated above the ash layer (Hiramine et al. in press). Other drift pumice clasts, concentrated between SKN-E-E8 and SKN-E-E9, were correlated with the U-2 tephra (5681-5619 cal yr BP) (Hiramine et al. in press). Furthermore, they suggested that the depositional timings of these pumice concentrated layers are shortly after the eruptions, and thus, their tephra dates can be used as ages of the horizons.

SKN-W outcrop
The SKN-W outcrop ( Fig. 4) is ~ 150 m east of the mouth of the Bitsuke River and is exposed over a 2 m height and a 1 m width. There are 17 sandy event deposits (SKN-W-E1 to SKN-W-E17, from youngest to oldest) interbedded within organic silt and peat ( Fig. 3; Table 2). For observing the uppermost sandy event deposit, two GS cores (SKN-W-GS1 and SKN-W-GS2) were obtained behind the outcrop. Except for SKN-W-E1 and SKN-W-E8, which are notably thick (> 15 cm), the sandy event deposits are ~ 10 cm thick or less. Most of the sandy event deposits consist of well-sorted fine-to-medium sand. SKN-W-E1, SKN-W-E2, SKN-W-E11, and SKN-W-E14 include coarse sand. Our detailed observation revealed that the 1-2-cm-thick black layer between SKN-W-E3 and SKN-W-E4 was a sand layer with a little organic material rather than soil. Therefore, SKN-W-E3 and SKN-W-E4 were determined to be a single sandy event deposit (hereinafter SKN-W-E3/E4). The bases of SKN-W-E3/E4, SKN-W-E11, and SKN-W-E14 are erosive contacts, but the upper and lower boundaries   of other event deposits are unclear. We recognized normal grading in SKN-W-E11, SKN-W-E14, and SKN-W-E16, whereas SKN-W-E1 and SKN-W-E8 characteristically have multiple parallel laminae. The GS cores behind the outcrop also show alternating layers of light and heavy minerals in the lower part of SKN-W-E1, while the upper part is massive owing to bioturbation (Fig. 4c).
In the grain-size analysis (Fig. 3), almost all event deposits show a mode at 1-2.5 phi, but their distribution ranges are variable. Below SKN-W-E14, event deposits contained more fine sand fractions compared to those above it. The event deposits in SKN-W-GS2 show a very sharp mode from SKN-W-E1 to SKN-W-E5. The upper part of SKN-W-E1 has more fine-grained material compared to its lower part, possibly influenced by bioturbation and aeolian deposits.
We recognized the volcanic ash layer between SKN-W-E5 and SKN-W-E6 that was already correlated with B-Tm tephra (946 CE; Oppenheimer et al., 2017) by Kobayashi (2019), who also detected volcanic glass shards of To-a tephra (915 CE; Machida and Arai 2013) from the lower part of B-Tm tephra, suggesting that the Vertical axis indicates the probability of particles. Pie charts show gravel, sand, and mud contents. N means the number of particle samples measured. Mean, median, mode, skewness, and kurtosis of each histogram are shown in Additional file 2: Table S1 Page 10 of 27 Ishimura et al. Progress in Earth and Planetary Science (2022) 9:69 age of sediments above and below this tephra is post-946 CE and pre-915 CE, respectively.

Cores drilled in the incised valley
Nine GS cores (SKN-GS1 to SKN-GS9) were obtained 200-400 m inland from the shore (Figs. 1c, 5). Coring depths were 1.3-2.8 m. We cut 50-cm-length block samples from each core for nondestructive analysis. Except for SKN-GS2, the sediments mainly consisted of organicrich sediments (marsh deposits) with interbedded sandy event deposits. Only SKN-GS2 contains a massive grayish-white clay at deeper than ~ 2.6 m depth. We recognized sandy event deposits consisting of a single unit or multiple units. Many of the single-unit deposits were well sorted but lacked any sedimentary structures. SKN-GS1-E3, SKN-GS6-E3, and SKN-GS6-E-6 characteristically have laminae that tilt inland. We interpreted that these laminae may be cross lamination by considering the relationship between the bedding plane and laminae and by taking into account the some deformations (e.g., dragging) due to GS coring, although the observation window of GS cores is narrow compared to outcrops. Therefore, they indicate that the paleocurrent direction was from seaside to inland (Fig. 6). The multiunit deposits were further divided into thin (< 10-20 cm) and thick (> 35 cm) deposits. Among the multiunit event deposits, thin layers were indistinct to the naked eye but could be identified clearly on XCT images. The uppermost event deposit (E1 in each core) consists of a lower visible sand layer and one or two unclear sand layers that can be recognized in XCT images taken at the seaward sites but that show as a single sand layer at the inland sites (Figs. 5,6).   Ishimura et al. Progress in Earth and Planetary Science (2022) 9:69 From the tephra analysis, we identified visible tephra layers and volcanic glass concentrated horizons (Table 3). Visible tephra layers were identified below E1 deposits in all SKN-GS cores except for SKN-GS4 and SKN-GS5, for which continuous sampling revealed volcanic glass concentrated horizons (Fig. 5, Table 3). The volcanic glass shards in the visible tephra layers and volcanic glass concentrated horizons below E1 deposits in all cores show mainly platy-type morphology. However, the horizon at 170-172 cm depth in SKN-GS5 showed the highest volcanic glass content within 163-186 cm depth (Table 3). The volcanic glass shards at that horizon have a pumiceous (highly vesiculated) morphology.
The results of the grain-size analyses of the SKN-GS samples (Fig. 7) show more variety than those of the coastal outcrops SKN-E and SKN-W. Three types were observed: sandy event deposits with a 1-2-phi mode and a narrow grain-size distribution, those with a 3-4-phi mode and a wide grain-size distribution, and those with a bimodal mode (1-2 phi and 3-4 phi). Some SKN-GS samples contained a large amount of mud, but this may be owing to the bioturbation and indistinct sand layers.
The uppermost E1 deposit in all cores with the same stratigraphic position (above the visible tephra or volcanic glass concentrated horizon) has the same grain-size distributions.
According to the ITRAX analysis results ( Fig. 8 and Additional file 2: Fig. S1) for the 12 elements that we investigated, the values of Si, K, Ca, Ti, Mn, Fe, and Sr are higher in the sand layers. In the visible tephra layers and volcanic glass concentrated horizons below E1 deposits, K and trace elements (Rb, Y, and Zr) characteristically increase. Data for Ni are somewhat noisy but are macroscopically stable. S is slightly elevated in a few event deposits (SKN-GS1-E3, SKN-GS1-E6, and SKN-GS6-E6 in Fig. 8).

Grain-size analysis of present-day sediments
We did grain-size analysis of present-day sediments (Figs. 1c,9) surrounding the study site to compare with event deposits in the outcrops and SKN-GS cores. We included foreshore and backshore beach sediments, aeolian sediments, Pleistocene marine-terrace deposits, and fluvial sediments in the incised valley where we obtained SKN-GS cores. We first observed the grain compositions and particle shapes of present-day sediments, which were composed mainly of quartz, pyroxenes, and rounded rock fragments with weathering rind. The volcanic rocks and tephras from the volcanoes located at the west of the study site contain quartz and pyroxenes (Tsushima and Takizawa 1977); thus, sources of the sediments are primarily the ejecta of volcanoes. Furthermore, our findings suggest no substantial differences in the type and shape of the particles that comprise any deposits; therefore, the grain-size distribution is considered important for determining the source of the sediments. The foreshore sediments show a mode in the range 1.5-2.5 phi and are characterized by a narrow grainsize distribution. The backshore sediments contain finer-grained particles (3-4 phi), which may reflect a contribution of aeolian sand. Aeolian sediments on the Pleistocene marine terrace have multiple modes: 2-phi, 2.5-phi, and 3-4-phi particles. These finer-grained particles are thought to be the result of wind-driven transport. The mode of the Pleistocene marine-terrace deposits is similar to that of present-day beach sediments but has a wider grain-size distribution and contains coarser particles. Fluvial sediments from the creek in the incised valley show a similar size distribution to that of the Pleistocene marine-terrace deposits, which is consistent with the creek's having incised the terrace.

Tephra correlation
Using the major element compositions of volcanic glass shards (Table 4; Additional file 2: Fig. S2), the visible tephra at the 88-89 cm depth in SKN-GS1 was correlated with B-Tm tephra in SKN-W, and the volcanic glass concentration horizon at the 170-172 cm depth in SKN-GS5 included volcanic glass shards of To-Cu tephra. The former result means that the visible tephra layers and volcanic glass concentrated horizons below E1 deposits in all GS cores are B-Tm tephra based on the stratigraphic position and similarity of the morphology of volcanic glass shards. Additionally, Nara et al. (2021) reported that the B-Tm tephra characteristically included K 2 O and trace elements, which is consistent with our tephra correlations and ITRAX results ( Fig. 8 and Additional file 2: Fig. S1). The latter means that the horizon at the 170-172 cm depth in SKN-GS5 is interpreted to be deposited at least after the To-Cu tephra (5986-5899 cal yr BP) as there is a possibility of its secondary deposition. Although two visibly separate volcanic ash layers were identified at the 108-111 cm and 115-117 cm depths in SKN-GS8, their volcanic glass morphology and grain  Ishimura et al. Progress in Earth and Planetary Science (2022) 9:69 compositions are similar, and therefore, we interpreted that both are B-Tm tephra and the lower tephra might have migrated downward by bioturbation.

Cross-checking radiocarbon dates of organic bulk sediments
Because appropriate measurable organic materials (e.g., seed and leaf ) were not available, we used organic bulk sediments for radiocarbon dating (Table 5). When dating organic bulk sediments, if the origin of the organic material is uncertain, it must be cross-checked with known ages. Therefore, in this study, the organic bulk sediments were dated in the To-Cu horizon of SKN-GS5 and horizons above and below the B-Tm tephra in SKN-GS2. As a result, the dates of the upper and lower parts of the B-Tm tephra in SKN-GS2 were 1003 ± 21 yr BP (960-800 cal yr BP) and 1239 ± 21 yr BP (1270-1070 cal yr BP), respectively, which are consistent with the tephra age. However, the To-Cu volcanic glass concentrated horizon of SKN-GS5 has an age of 6673 ± 29 yr BP (7610-7470 cal yr BP), which is older than the To-Cu tephra (5986-5899 cal yr BP). Because we expected that the horizon contained secondary deposits of To-Cu tephra, we assumed that a younger age than the tephra would not be a problem, but the obtained age is older than the tephra age. Therefore, we considered the possibility that old carbon was included in the horizon. Although by visual inspection we thought the sedimentary facies consisted of clay-rich sediments (Fig. 6), the sedimentary structure from the XCT indicated that it was a poorly sorted clay with sand, indicating that the sediments were transported by fluvial processes and contained old carbon. The horizons above and below B-Tm tephra are organic-rich sediments containing few sand grains (Fig. 6). Because organic-rich sediments deposited in a stable environment provide a more accurate age, we chose the horizon with little old organic carbon for radiocarbon dating of organic bulk sediments. We determined that horizons with less flow-in by fluvial processes (i.e., with fewer clastic particles) were those with less old organic carbon and used XCT images to select such horizons. Therefore, some radiocarbon dates were obtained at some distance from the event deposits (Fig. 5).

Paleogeography of Sekinehama
Before the identification of the tsunami deposits, we looked at the topographic change around the study site since ~ 6 ka. Using an aerial photograph taken in 1952 and a corresponding geomorphic classification (Fig. 10a,  b), we reconstructed the paleogeography of Sekinehama. The 1952 aerial photograph shows that dunes used to be distributed along the coast where currently the front of the coastal outcrops are located, while the other aerial photographs show that the topography in 2001 was almost the same as today. This suggests that the outcrops that we studied appeared at the surface owing to coastal erosion in the last 50 years. Thus, until more than 50 years ago, these outcrops had been in a lowland behind the coastal dune and unexposed. Hiramine et al. (in press) found two layers composed of abundant drift pumice in SKN-E and concluded that SKN-E was under a lagoon-like wetland environment connected to the sea via a small channel at ~ 6 ka (Fig. 10c). The presence of drift pumice and underlying marine sediments at SKN-E roughly indicates a paleosea level at 2 m above the present sea level (Hiramine et al. in press). Okuno et al. (2014) showed that theoretical relative sea level is 0-2 m in 6 ka at Ominato (15 km southwest from the study site), indicating that study cite is stable or slightly uplifted during the last 6 ka. This tendency is consistent with the long-term uplift rate of 0.1-0.2 mm/yr based on the MIS 5 marine-terrace surfaces (Koike and Machida 2001). Even if this uplift is taken into account, the sediments in the GS cores are considered to be higher in elevation than the sea level at that time. From the above-estimated paleogeography, the coastline at ~ 6 ka was estimated to be more seaward than at present, indicating that the present coastline position is the nearest position to the GS core sites during the last 6 kyrs. Additionally, the layers of the drift pumice indicate the paleoshoreline position of ~ 6 ka in the inland lagoonal environment (SKN-E), meaning that there is a distance of at least 200 m from the paleoshoreline of the lagoon behind the spit or dune (Figs. 10c, 11) to the GS core site and that there is a greater distance from the coastline facing the sea (Tsugaru Strait). Although we could not collect additional data on the paleogeography after 6 ka, we interpreted that the study area has tended to be land from 6 ka; therefore, the marine influence has decreased. Additionally, the estimated tidal inlet was limited, and the GS core site was far from it (Fig. 10c). Therefore, a storm is not likely to reach the GS core site after 6 ka.

Interpretation of event deposits in GS cores
We compared source sediments surrounding the study site with event deposits in the GS cores. From the grain-size analysis of source sediments, only the beach sediments show a mode of 1.5-2.5 phi with a narrow grain-size distribution (Fig. 9). This indicates that the sandy event deposits with a grain-size distribution similar to that of present-day beach sediments might have been transported from the beach (pink histograms in Fig. 7). However, coarse-grained sandy event deposits (including -1 to 1 phi particles) with a wider grain-size distribution likely originated from the Pleistocene marine-terrace deposits and are similar to creek sediments. Therefore, they are estimated to be deposited by fluvial processes. Conversely, fine-grained particles (3-4 phi) could be derived from fine sediments above and below event deposits due to bioturbation. If they are mixed by bioturbation to event deposits, especially thin event deposits, Pie charts show gravel, sand, and mud contents. N means the number of particle samples measured. In the histogram, blue means that the size distribution is different from that of present-day beach sediments (presented in Fig. 9); pink means that it is similar to that of present-day beach sediments; and orange means that it is somewhat similar. Mean, median, mode, skewness, and kurtosis of each histogram are shown in Additional file 2: Table S1 Page 15 of 27 Ishimura et al. Progress in Earth and Planetary Science (2022)  Grain-size analysis of present-day sediments. Histograms show phi-scale grain-size distribution. Vertical axis indicates the probability of particles. N means the number of particle samples measured. Mean, median, mode, skewness, and kurtosis of each histogram are shown in Additional file 2: Table S1 Page 17 of 27 Ishimura et al. Progress in Earth and Planetary Science (2022) 9:69 grain-size distribution becomes a little wider than the present beach sediments, such as SKN-GS1-E6, or shows bimodal distribution, such as SKN-GS8-E5. If the contribution of fine-grained particles is minor, this possibility should be considered in the determination of source sediments. Therefore, for thin event deposits, it is possible that the fine-grained component, if present, originally had a grain-size distribution similar to that of the present beach sand (orange histograms in Fig. 7). These two event-deposit types (pink and orange histograms in Fig. 7) are candidates for tsunami deposits, as discussed below. Some event deposits whose grain-size distributions are similar to present-day beach sediments show cross lamination tilting inland (SKN-GS1-E3, SKN-GS6-E3, and SKN-GS6-E6), indicating a paleocurrent directed inland. Furthermore, SKN-GS1-E3, SKN-GS6-E3, and SKN-GS6-E6 show sharp basal contacts, indicating erosional surfaces. These features indicate that they were deposited by the strong flow from the coast. Considering the distance from the coast and estimated paleogeography since ~ 6 ka (Fig. 10c) as discussed in Sect. 5.3, they were deposited by tsunamis rather than storm surges. The E1 deposits above the B-Tm tephra in each core also show almost the same grain-size distribution as the present-day beach sediments, thin inland, and have multiple sand layers. These characteristics also suggest that they are a tsunami deposit and were formed by multiple waves. Because other event deposits do not show any relevant features other than grain-size distributions, we interpret them based on their distribution in the incised valley. We discuss their correlations between GS cores using radiocarbon dating later in this section.
SKN-GS1-E4, SKN-GS4-E6, SKN-GS7-E3, and SKN-GS8-E4 are characterized by cross laminations showing a seaward paleocurrent direction and have thicknesses of > 35 cm. Almost of their grain-size distributions are similar to the present creek sediments (Fig. 9). From these characteristics of sediments and the present creek topography (~ 1 m wide and 30-50 cm deep), they are considered to be buried channel fills. However, considering the tsunami behavior, such channel-fill deposits could have been formed during a tsunami return-wave period and thus could be either fluvial or tsunami deposits. The grain-size distributions of lower part of SKN-GS4-E6 (Fig. 7) showed similar distribution of present-day beach sediments, indicating the possibility of tsunami return-wave sediments. Other event deposits with a wide grain-size distribution and large amounts of particles both coarser and finer grain sizes than 1-2 phi are more likely to be fluvial instead of tsunami deposits.
Finally, we conducted radiocarbon dating to determine the correlation between tsunami deposits as inferred above and candidates in the GS cores in the incised valley. As a result, we identified five tsunami deposits (TD1-TD5). The modeled ages of each tsunami deposit obtained using OxCal are shown in Table 6 and Additional file 2: Fig. S3. TD1 (500-300 cal yr   Page 20 of 27 Ishimura et al. Progress in Earth and Planetary Science (2022) 9:69 BP (1450-1650 cal CE)) is the latest tsunami deposit (E1 deposit in all cores) in this study and is distributed above B-Tm tephra. TD1 in many of the cores shows a grain-size distribution similar to that of the beach sediments, with multiple units and thinning inland. TD2 (4490-4300 cal yr BP; SKN-GS1-E2, SKN-GS8-E2, and SKN-GS9-E3), found up to SKN-GS1 and not farther inland, shows a grain-size distribution similar to that of the beach sediments, with an erosive contact at the base in SKN-GS1. TD3 (4790-4580 cal yr BP; SKN-GS1-E3, SKN-GS6-E2, SKN-GS8-E3, and SKN-GS9-E4) shows a grain-size distribution similar to that of the beach sediments, with an erosive contact at the base in SKN-GS1. The sedimentary structure of SKN-GS1-E3 shows a paleocurrent direction inland. SKN-GS7-E3 was a channel-fill deposit that may have been formed by the return wave during the TD3 deposition. TD4 (5420-5030 cal yr BP; SKN-GS2-E2 and SKN-GS6-E3) has a slightly finergrain-size distribution in the beach sediments. In SKN-GS6-E3, the paleocurrent direction inland was estimated from the sedimentary structure. SKN-GS1-E4 and SKN-GS8-E4 were channel-fill deposits and may have been formed by the return wave during TD4 deposition. TD5 (6030-5790 cal yr BP; SKN-GS1-E5, SKN-GS2-E5, SKN-GS6-E6, and SKN-GS8-E5) shows a grain-size distribution similar to that of the beach sediments, with an erosive contact at the base. The sedimentary structure of SKN-GS6-E6 shows a paleocurrent direction inland. SKN-GS4-E6 was channel-fill deposit and may have been formed by the return wave during TD5 deposition. In the  Page 21 of 27 Ishimura et al. Progress in Earth and Planetary Science (2022) 9:69 ITRAX results (Fig. 8), some sand layers (e.g., SKN-GS1-E3 and SKN-GS6-E6) show an increase in S, which is an indicator of marine intrusion (e.g., Chagué-Goff et al. 2016), supporting the above identifications of tsunami deposits.
In addition, SKN-GS1-E6, SKN-GS6-E7, and SKN-GS8-E6 may be correlated. SKN-GS6-E7 shows a finer-grain-size distribution than that of the beach sediments, while SKN-GS1-E6 and SKN-GS8-E6 show grain-size distributions similar to those of the beach sediments. The ITRAX results show a slight increase in S at SKN-GS1-E6 (Fig. 8). However, because they are near the base of the cores, and the underlying sediments were not fully observable, they are not identified here as tsunami deposits.

Interpretation of event deposits in SKN-E and SKN-W outcrops
Candidates for tsunami deposits in SKN-W are SKN-W-E3/E4, SKN-W-E10, SKN-W-E11, and SKN-W-E14 because all of them have erosional contacts at the base and are distinct sand layers compared to the other event deposits (Table 2). In contrast, SKN-W-E1 and SKN-W-E8 characteristically exhibit thicker layers and parallel laminae. These features are commonly observed in washover deposits caused by storm surges behind the coast (e.g., Morton et al. 2007). Additionally, SKN-W-E1 is distributed up to the ground surface (Fig. 4). Considering the most recent topographic changes (Fig. 10a, b), it could be a washover deposit formed by a storm surge that eroded the coastal sand dune in the past 50 years. In SKN-E, SKN-E-E5, SKN-E-E6, and SKN-E-E8 are candidates for tsunami deposits because all of them have erosive contact at the base with rip-up clasts. The rip-up clast is often used as one of the indicators for identifying tsunami deposits (e.g., Goff et al. 2012) and for distinguishing tsunami and storm-surge deposits (e.g., Kortekaas and Dawson 2007). SKN-E-E5 in particular includes a few centimeters of gravels and pumice clasts, suggesting a stronger flow compared to other event deposits in the outcrop. SKN-E-E7 is a channel-fill deposit with parallel laminae. Taking into account the paleogeography (Fig. 10c), it could be fluvial sediments that buried the creek flowing through the wetland or tidal-channel sediments that filled the channel connecting to the sea.
To correlate the sediments between the two outcrops, we used radiocarbon dating on deposits below the tsunami-deposit candidates mentioned above, which are expected to be widely distributed in the study area. Dates   Ishimura et al. Progress in Earth and Planetary Science (2022) 9:69 for samples collected below SKN-E-E5 and SKN-W-E14 are 3970 ± 27 yr BP and 3992 ± 26 yr BP, respectively. This coincidence means that they are the sandy deposits related to the same event. These dates added to the tephra ages and other radiocarbon dates confirmed having continuous sediments from the present to 6300 cal yr BP exposed in the outcrops. Grain-size analyses of the outcrop samples show that almost all the event deposits have a mode at 1-2 phi, similar to the present-day beach (Figs. 3,9). In the SKN-E outcrop, SKN-E-E1 to SKN-E-E4 (younger than ~ 4500 cal yr BP) show a narrow grain-size distribution with a 1-2-phi mode. In contrast, SKN-E-E5 to SKN-E-E10 (4500-6000 cal yr BP) exhibit 1-2-phi modes but wider grain-size distributions. Below SKN-E-E11 (older than ~ 6000 cal yr BP), the mode shifts to a somewhat-smaller size (2 phi), and the proportion of fine-grained particles increases. A similar trend of grain-size variation was observed in the SKN-W outcrop. From SKN-W-E1 to SKN-W-E13 (present to 4500 cal yr BP), they show narrow grain-size distributions with 1-2-phi modes. SKN-W-E14 to SKN-W-E16 (older than 4500 cal yr BP) show 1-2-phi modes but with wider grain-size distributions. SKN-W-E17 has a smaller mode (2 phi) and contains a higher proportion of fine-grained particles compared to the event deposits above it. This similarity in the grain-size distribution of the event deposits at the two sites according to age suggests the existence of common event deposits at the two sites. However, the numbers of event deposits between the two sites do not necessarily match, and correlations with radiocarbon dates are difficult because of the close distributions of event deposits. In Sect. 5.6, we correlate the tsunami deposits (TD1-TD5) identified in the SKN-GS cores to the event deposits in these outcrops.

Tsunami deposits distribution in Sekinehama
From the dating results, we correlated tsunami-deposit candidates among the coastal outcrops to TD1-TD5 in the GS cores (Figs. 2, 4, 11). TD1 correlates with SKN-W-E3/E4, which has an erosional base, and was traced in the SKN-W-GS cores drilled behind the outcrop. In outcrop, they present as two separate units, indicating the multiple waves. This feature is consistent with the sedimentary facies of E1 deposits of the SKN-GS cores. Accordingly, TD1 was deposited by a tsunami with multiple waves during 500-300 cal yr BP (1450-1650 cal CE).
TD2, TD3, and TD4 were correlated with SKN-E-E5, SKN-E-E6, and SKN-E-E8, respectively. Because these three layers have coarse grains, an erosional base, and rip-up clasts, we consider them to have been deposited by strong wave currents. Considering the radiocarbon dates, TD5 may correlate with SKN-E-E9 or SKN-E-E10.
In SKN-W, we correlated SKN-W-E14 to TD2. SKN-W-E16 may correspond to TD3, but the age could not be measured owing to the low organic carbon content. Considering the age of the sediments, it is unlikely that TD4 and TD5 are represented in SKN-W.
These results allowed us to correlate tsunami deposits from the coastal outcrops to deposits 200-400 m inland (Fig. 11). These distributions also indicate that TD1-TD5 were deposited by large tsunamis compared to historical tsunamis (Table 1). TD1, TD4, and TD5 are distributed the farthest inland among the five tsunami deposits. However, the fact that the sediments in the coastal outcrops were deposited in lowlands behind dunes and that the other event deposits not correlated with TD1-TD5 show similar grain-size distributions with those of the beach sediments suggests that they may have been deposited by smaller tsunamis or by low-frequency stormsurge events. This is consistent with our interpretation that the thick event deposits with multiple laminations and limited distributions (SKN-W-E1 and SKN-W-E8) are washover deposits associated with storm surges.
Between TD2 and B-Tm tephra, multiple event deposits are commonly recognized in the coastal outcrops, whereas no distinct event deposits are recognized between them in the GS cores. The GS cores yield dates of ~ 3000 cal yr BP ~ 50 cm below the B-Tm tephra. The sedimentary facies show no unconformity below B-Tm, but this indicates that the sedimentation rate was low or it represents a hiatus. Given the absence of event deposits between TD2 and B-Tm tephra in the GS core, apparently event deposits either are distributed only near the coast or are not well preserved in the GS cores.

Tsunami deposits correlation in Tohoku and Hokkaido regions
Considering the correlation to more-distant tsunami deposits, TD1 (500-300 cal yr BP (1450-1650 cal CE)) can be correlated with the 1611 Keicho tsunami deposits on the Sanriku Coast, the seventeenth-century tsunami deposits in western Hokkaido, and the 1640 Komagatake tsunami deposits. The 1454 Kyotoku tsunami (Namegaya and Yata 2014) also appeared to be a possible correlative event linked to TD1, but Ishizawa et al. (2022) showed that there is no deposit from the 1454 Kyotoku tsunami on the northern Sanriku Coast (Fig. 1a). As for the 1640 Komagatake tsunami, numerical simulations conducted by Nakanishi and Okamura (2019) to explain the tsunami-deposit distributions show that with distance from the source the wave heights decreased. As Sekinehama is ~ 100 km from the Komagatake, the calculated tsunami height was < 5 m, which is below the present sea-cliff height of ~ 5 m (Fig. 11), meaning that the 1640 Komagatake tsunami likely did not inundate the locations of the Page 23 of 27 Ishimura et al. Progress in Earth and Planetary Science (2022) 9:69 studied GS cores. Also, Sato et al. (in press in this special issue) confirmed that the simulated 1640 Komagatake tsunami would not have inundated the GS core sites in our study.
To explain the tsunami deposits near the study site (Minoura et al. 2013;Tanigawa et al. 2014;Kase et al. 2022; Fig. 1), it was necessary to assume the occurrence of an unknown tsunami in the fifteenth century, considering that only one layer of tsunami deposits at each site has been identified during the last 1500 years (Fig. 12), the topographical settings (distance from the coastline and elevation) of each site, and the uncertainties of ages due to differences in dating materials. We determined the depositional ages of the latest tsunami deposits within a 50-km distance from Sekinehama (Figs. 1, 12) to be as follows: Ta1 deposit (589-516 cal yr BP (1361-1434 cal CE)) at the Takadai-cho site in Hokkaido (Elevation is 6.4 m, and distance from the coast is 200 m; Kase et al. 2022), TD1 (500-300 cal yr BP (1450-1650 cal CE)) in this study (Elevation is ~ 5-7 m, and distance from the coast is ~ 200-400 m), a tsunami deposit (681-490 cal yr BP (1269-1460 cal CE)) on the east coast of Shimokita Peninsula (Elevation is ~ 21-22 m, and distance from the coast is ~ 1350 m; Minoura et al. 2013), and S1 deposit (470-180 cal yr BP (1480-1770 cal CE)) on the east coast of Shimokita Peninsula (Elevation is ~ 3-5.5 m, and distance from the coast is ~ 700-1050 m; Tanigawa et al. 2014.). Although not completely overlapping, they are very close to the mid-fifteenth century (Fig. 12). If they were all deposited by the same tsunami event, we  Machida and Arai (2013). For study-site locations, see Fig. 1 Page 24 of 27 Ishimura et al. Progress in Earth and Planetary Science (2022) 9:69 assume they were deposited in the mid-fifteenth century. Accordingly, and given the seventeenth-century tsunami in Hokkaido and the 1611 Keicho tsunami on the Sanriku Coast, the latest tsunami deposits at Sekinehama could have been deposited by an indeterminate fifteenth-century tsunami that has yet to be recognized in the historical record. We correlated TD2-TD5 to tsunami deposits identified near the Sekinehama Coast. On the other side of the Tsugaru Strait, in addition to Ta1 mentioned above, Kase et al. (2022) identified two tsunami deposits (Ta2 and Ta3) at the Takadai-cho site (Fig. 1a), which is the closest to our study site among their study sites. The ages of Ta2 and Ta3 are 1745-1639 cal yr BP and 2401-2265 cal yr BP, respectively (Fig. 12). Although they do not correlate with TD2-TD5, radiocarbon dates indicate that they may correlate with SKN-W-E7 to SKN-W-E10. However, Kase et al. (2022) stated that the areal distributions of Ta1-Ta3 were not confirmed and were limited to a narrow area of 10 m or less. Therefore, more excavation and drilling surveys similar to our current study are needed to confirm the inland distribution and correlate laterally.
In addition to the already-described S1, Tanigawa et al. (2014) identified four sandy event deposits (S2-S5) that formed during the last 6000 years and concluded that S1, S3, and S4 were deposited by seawater inundation on the basis of diatom analysis. The ages of the four deposits are as follows (Fig. 12): S2 (4800-4500 cal yr BP), S3 (5300-4870 cal yr BP), S4 (5540-5320 cal yr BP), and S5 (5890-5490 cal yr BP). We correlated TD4 (5420-5030 cal yr BP) to S3 or S4. If S2 and S5 are also tsunami deposits, we would correlate TD3 (4790-4580 cal yr BP) to S2 and TD5 (6030-5790 cal yr BP) to S5. Tanigawa et al. (2014) also found no seawater inundation over a long period between S1 and S2 or S3. Thus, the depositional ages and intervals that they determined for the tsunami deposits are similar to those in our study. In addition, S1 and S3 are distributed farthest inland in the drilling area, which is also consistent with the distribution of TD1 and TD4 in our study (Fig. 11).
We also investigated correlation to deposits in more-distant areas. West of the Hidaka Coast in Hokkaido, tsunami deposit during 4000-6000 cal yr BP is unknown (Takashimizu, 2014). In the southern Hidaka Coast area, dates for tsunami deposits that formed during 5.5-2.2 ka were obtained by Nakanishi et al. (2020) (Fig. 12). Although tsunami deposits that correlate with our results are distributed, Nakanishi et al. (2021) considered that their tsunami sources were likely to be along the Kuril Trench (Fig. 1a). Therefore, to confirm correlation, tsunami sources need to be addressed. Similarly, although TD2-TD5 could be correlated with tsunami deposits identified along the Sanriku Coast (Fig. 12;Ishimura 2017;Ishimura and Yamada 2021;Ishizawa et al. 2022;Velasco-Reyes et al., 2022), we need to assess whether tsunamis off the Sanriku Coast could generate tsunamis inundating Sekinehama. Thus, although there is a possibility that the tsunami deposits in Sekinehama correlate with those in the Hokkaido and the Sanriku Coast, further investigation involving tsunami sources is needed to correlate the Sekinehama deposits to more-distant tsunami deposits.

Conclusions
In this study, we performed nondestructive analyses (XCT and ITRAX µXRF core scanning), grain-size analysis, tephra analysis, and radiocarbon dating of sediments exposed in the two coastal outcrops and GS cores at Sekinehama on the north coast of the Shimokita Peninsula, which is near the junction of the Japan and Kuril trenches. By comparing grain-size distribution between sandy event deposits and present-day source sediments and by investigating detailed sedimentary structures observed in XCT images, we found that some event deposits show transportation inward from the coast. Finally, we identified five tsunami deposits (TD1-TD5) that formed since ~ 6 ka and correlated them from coastal outcrops to the incised valley based on radiocarbon dating and tephrochronology. We also found that other washover deposits were among the coastal outcrops but not in the GS cores. This limited distribution indicates minor washover events related to small tsunamis and low-frequency storm surges.
The modeled age of the latest tsunami deposit is 500-300 cal yr BP (1450-1650 cal CE). This may be correlative with two known tsunamis (seventeenthcentury tsunami in Hokkaido and 1611 Keicho tsunami in Tohoku) or may be an unknown tsunami that possibly occurred during the fifteenth-seventeenth centuries. If the latest tsunami deposit is correlated with either or both of the seventeenth-century and 1611 Keicho tsunamis, the tsunami inundation area would be expanded, resulting in the need for a complex tsunami source model. However, if the latest tsunami deposit is correlated with tsunami deposits previously identified within a 50-km distance from the study site, we need to consider an unknown fifteenth-century tsunami and use a different source model set off the Shimokita Peninsula. Nonetheless, our determination of the latest tsunami deposit is one constraint with high reliability for constraining the tsunami source models in the future even if not correlated with any specific tsunamis.