EVIDENCE OF MASS MORTALITY OF THE LONG-LIVED BIVALVE MERCENARIA STIMPSONI CAUSED BY A CATASTROPHIC TSUNAMI

ABSTRACT Tsunamis are huge disasters that can significantly damage benthic organisms and the sea-bottom environment in coastal areas. It is of great ecological importance to understand how benthic ecosystems respond to such destructive forces and how individual species are affected. Investigating the effect of such disasters on animals that are seldom caught alive is particularly difficult. Bivalve mollusks are especially suitable for investigating how a tsunami affects coastal benthic species because they preserve an environmental record in their shells that can be extended back in time by crossdating the records of multiple individuals. Here we studied dead shells of Mercenaria stimpsoni, a long-lived clam, and precisely determined the time of death by using nuclear bomb–induced radiocarbon (bomb-14C) and by counting annual growth increments. First, a quasi-continuous, regional bomb-14C record was created by analyzing the shells of 6 live-caught M. stimpsoni individuals. Then 27 dead shells collected from the seafloor of Funakoshi Bay were 14C-dated and analyzed. The results showed that the huge tsunami that struck northeastern Japan on 11 March 2011 caused mass mortality of this bivalve in Funakoshi Bay. Nine of the 27 clams died during the March 2011 tsunami, probably by starvation after burial by tsunami deposits or exposure above the seafloor as a result of sediment liquefaction during the earthquake. The dating method used in this study can help us understand how long-lived marine organisms with low population density are affected by huge natural disasters such as a tsunami.


INTRODUCTION
Natural disasters such as tsunamis, earthquakes, typhoons, and volcanic ash falls can substantially disturb benthic ecosystems in coastal areas (López et al. 2008;Jaramillo et al. 2012;Seike et al. 2013;Harris 2014). In particular, tsunamis can cause either rapid sedimentation, thus burying benthic organisms deep under the seafloor, or rapid erosion, thus washing them out of the upper layer of substrate and exposing them to predation (Kranz 1974;Bromley 1996;López et al. 2008). On 11 March 2011, the Tohoku-Oki Earthquake, a M9.0 megathrust earthquake, ruptured the plate boundary off the Pacific coast of northeastern Japan (Ide et al. 2011) and generated a huge tsunami that affected 2000 km of the Pacific coast of northeastern Japan (Mori et al. 2011). In Funakoshi Bay, Iwate prefecture, the tsunami run-up height was 29.4 m (Ishimura and Yamada 2019), and largescale seafloor erosion and sedimentation occurred, followed by the re-establishment of bioturbation (Seike et al. 2016(Seike et al. , 2019. Furthermore, population densities of sand dollar Scaphechinus mirabilis (an echinoid) and Sakhalin surf clam Pseudocardium sachalinense (a mollusk) in the bay were significantly decreased by the tsunami disturbance (Seike et al. 2013). Although investigating the population dynamics of animals with a large population density before (and after) the tsunami is possible if the data from before the event are available, investigating those with a small population density and rare catch opportunity is extremely difficult. In addition, in most coastal areas, environmental or ecosystem survey records from before an event are unlikely to be available; as a result, it is usually very difficult or even impossible to assess the effect of a catastrophic event on coastal ecosystems. Some marine animals such as mollusks have a countinuously grown calcium carbonate shell. As a result, mollusk shells continuously record the ambient environment, making them excellent archives of past environments and environmental disturbances (e.g., Murakami-Sugihara et al. 2019). The cold-water venerid bivalve Mercenaria stimpsoni (Stimpson's hard clam), which lives in Funakoshi Bay, is a long-lived mollusk with a shell that displays distinct annual growth increments (Kubota et al. , 2018aShirai et al. 2018). It is thus possible to test the response of this marine mollusk to environmental perturbation, despite its low population density in the bay (Seike et al. 2013), because gaps in the shell record can be filled by using crossdating to combine the records of multiple individuals.
Mercenaria stimpsoni is the most long-lived among modern bivalves living near the Japanese coast , and its shell growth pattern and the paleoenvironmental implications of the pattern have been documented previously Kubota et al. 2017Kubota et al. , 2018aShirai et al. 2018). This species inhabits fine sandy substrata in the shallow coastal zone (water depth 5-30 m) along mid-to high-latitude coasts of the northwest Pacific, the Sea of Japan, and the Sea of Okhotsk (Habe 1977). It is an infaunal filter feeder that consumes phytoplankton (predominantly diatoms) and suspended organic matter Seike et al. 2017). At ontogenetically younger stages (<15 years old), the shells grow from early spring to late summer at a relatively high rate (∼10 mm/yr), but shell growth ceases or is substantially reduced during winter, when the water temperature is below ∼10°C . At ontogenetically older stages (>15 years), shell growth occurs only in summer and the growth rate is very low (0.05-0.1 mm/yr). Growth cessation lines are visible in shell cross sections, and the calendar year of past growth increments of a live-caught clam can be obtained by counting the annual growth lines backward from the most recent (Kubota et al. 2018a;Shirai et al. 2018).
Many dead M. stimpsoni shells can be found on the seafloor of Funakoshi Bay. While a modern specimen provides a record that extends back to when the specimen first developed a shell, by combining, through crossdating (e.g., Black et al. 2016), the annual growth increment records of modern specimens with those of specimens that lived during an earlier period, the paleoenvironmental record can be extended further back in time to produce a master chronology. However, to obtain a continuous record with annual resolution, we need to determine the dates of death of dead shells very precisely, with a precision of one year. However, crossdating of dead shells is more complicated than crossdating of live-caught specimens because the record of a dead shell, especially if it is short, may match several calendar year intervals in the master chronology. Conventional radiocarbon dating of marine samples has a large age determination error of more than ±50 years, which stems from the analytical precision and uncertainties of the seawater 14 C calibration curve (e.g., Marine20, Heaton et al. 2020) and the marine 14 C reservoir age (Alves et al. 2018). One way to overcome this problem is to use 14 C produced by atomic bomb testing (Hua et al. 2013). Intense atmospheric testing of thermonuclear weapons conducted in the 1950s and 1960s artificially produced 14 C through reactions of high-energy neutrons with atmospheric nitrogen (bomb-14 C). The bomb-14 C entered the Earth surface system and has already been widely used as a tracer for studying biogeochemical cycles and to estimate the longevity of both terrestrial and marine organisms (Kalish 1993;Campana et al. 2002;Hamady et al. 2014;Nielsen et al. 2016).
Although mixing occurs rapidly in the atmosphere, in seawater mixing occurs more slowly; as a result, bomb-14 C is heterogeneously distributed in the surface of the ocean (Mahadevan 2001;Druffel 2002;Grottoli and Eakin 2007;Scourse et al. 2012;Toggweiler et al. 2019). Therefore, to use the bomb-14 C record as an age-determination tool, a continuous regional bomb-14 C record must first be established. Kubota et al. (2018a) established a bomb-14 C record for Funakoshi Bay by crossdating the 6 live-caught M. stimpsoni individuals mentioned above (this study improved the curve by adding 4 new data points; see Online Supplementary Material).
In the present study, we used crossdating and the bomb-14 C record to determine the precise (annually resolved) ages of dead shells of M. stimpsoni collected from the seafloor of Funakoshi Bay. We then compared our analysis results with the timing and amplitude of past tsunamis to determine how they might have affected the mortality of this species in the bay. We also measured the 14 C content of an articulated M. stimpsoni shell found in a marine sediment core obtained from tsunami deposits in Funakoshi Bay in 2014. The goal of this study is to assess how past tsunamis, especially one occurred at 11 March 2011, affected the mortality of M. stimpsoni in Funakoshi Bay.

METHODS
Funakoshi Bay is located on a ria coast, a part of the Sanriku coast of northeastern Japan ( Figure 1). The coastal plain along this section of coast is narrow, and the bay, which faces the open ocean, is directly affected by ocean wave activity. The floor of the bay is covered by poorly sorted sand and inclines gradually toward the open ocean (Seike et al. 2013). Otsuchi Bay is a rectangular bay, located immediately south of Funakoshi Bay, that measures 8 km from east to west and 2 km from north to south ( Figure 1). Water depth increases gradually toward the bay mouth, where it reaches 70 m. Otsuchi Bay is fed by three main rivers (the Otsuchi, Kozuchi, and Unosumai Rivers). Each river discharges freshwater in the bay at a rate of 3-10 m 3 /s (Kubota et al. 2018b).
In this study, we used 5 live M. stimpsoni specimens collected from the seafloor of Funakoshi Bay (site F-1, 20 m water depth; Figure 1) and 1 from Otsuchi Bay (site O-3, 5 m water depth) during a dive or by using a Smith-McIntyre mud sampler (see Online Supplementary Material) (Kubota et al. 2018a). These specimens were frozen after recovery and then transported to the laboratory for the removal of their soft tissues. In the laboratory, the shells were cleaned and dried at room temperature. We also collected both the right and left valves of dead shells (in rare cases articulate shells, too) of M. stimpsoni from the seafloor at site F-1 in Funakoshi Bay ( Figure 1) during SCUBA dives. We selected 27 right valves of the dead shells for further analysis to avoid duplication of the analysis in the same individuals.
For the analysis we also used one articulated M. stimpsoni specimen (see Online Supplementary Material) that was removed from a sediment core (6 cm in diameter, 100 cm long) collected from the seafloor of Funakoshi Bay (F-1) in September 2014 ). The removed M. stimpsoni specimen was stored in a plastic bag until analysis. Sedimentary structures preserved in the core have already been reported in Seike et al. (2017). Briefly, the structures were examined by using a computed tomography (CT) scanner (LightSpeed Ultra16, GE Healthcare Japan Corp.) at the Center for Advanced Marine Core Research, Kochi University, Kochi, Japan. To observe both physical and biogenic sedimentary structures in the core, the CT scans were visualized by using OsiriX imaging software (version 4.1.2). After the core was scanned, it was sliced into 5-cm-long subsamples for analysis of the grain-size distribution with a laser granulometer (Mastersizer 2000, Malvern Instruments) at the Center for Advanced Marine Core Research.
The preparation of the live specimens for sclerochronological analysis was described by Kubota et al. (2018a), and the dead shells were prepared following the same method. First, we carefully removed materials adhering to the dead shell surfaces using a brush and then dried the shells at room temperature. From the 27 individuals, we selected 22 with relatively thick shells (i.e., with greater longevity), because short-lived individuals are not suitable for sclerochronological analysis. Then, each valve was wrapped in epoxy resin and sectioned from the umbo to the ventral margin along the maximum growth axis with a low-speed saw (Buehler Isomet). A cross-section was prepared and polished, and then images of the section were taken at 100× magnification with a KEYENCE digital microscope.
On the photographic image of each cross section, we counted the number of annual increments and measured their width by using ImageJ freeware (http://imagej.nih.gov/ij/). Growth in M. stimpsoni is similar to growth in trees in that growth increments vary in width and show strong ontogenetic trends Kubota et al. 2017;Shirai et al. 2018). We used the ARSTAN software package (https://www.ldeo.columbia.edu/tree-ring-laboratory/resources/ software), following Butler et al. (2009aButler et al. ( , 2009b, to produce a chronology by detrending and standardizing the time series data on annual growth increment widths. A negative exponential function was fitted to the time series of annual growth increment widths and then the residuals from the fitted function were calculated. We used the following equation to calculate the growth index (GI) for each year for each shell: where m t is the measured annual increment width, and p t is the annual increment predicted from the fitted curve, at age t. Annual increment width records of modern shells has already been reported in Kubota et al. (2018a), but those of dead shell are reported in this study for the first time. Then we calculated the standardized shell growth index (SGI; Jones et al. 1989) from the GI: where x is the average and σ is the standard deviation of all GI values of the time series. The SGI is a dimensionless measure of shell growth that can be used for direct comparisons of shell growth between different specimens and different ontogenetic stages. According to Black et al. (2016), we visually crossdated both the live-caught and dead shells to produce a master chronology for the study area. We also calculated the expressed population signal (EPS) to evaluate the strength of the correlations among the individual time series and between those time series and the master chronology (Wigley et al. 1984). Window length was set to 10 years in this exercise.
In this study, we used the time series data of annual growth increment widths for the shells of the 6 live-caught clams reported by Kubota et al. (2018a), together with data from 15 dead shells (among the 27 14 C-dated dead shells), reported here for the first time. In 22 dead shells that were sclerochronologicaly analyzed 7 showed little annual variation and was not able to obtain dates, thus we did not include them in the master chronology. A cross section of the articulated M. stimpsoni specimen was not made because the shell was broken during coring (see Online Supplementary Material).
To obtain the most recent record from each shell, we traced the surface of the last-precipitated portion of the shell (last increment; see Online Supplementary Material) with a grinder to remove surface contamination and then collected calcium carbonate powder from the fresh surface. We also sampled the articulated M. stimpsoni from the sediment core in this manner.
We first photographed each cross section. Then we used a high-precision micromilling system (Geomill326, Izumo Web, Japan) and a triangular tungsten carbide drill bit (Shofu) and drilled into the outer shell layer to obtain samples of the annual growth layers for 14 C analysis (e.g., Kubota et al. 2015aKubota et al. , 2015bKubota et al. , 2017Sakamoto et al. 2017 Most of the radiocarbon analysis results for the live-caught shells (20 of 24 14 C data points) have already been reported by Kubota et al. (2018a), but all radiocarbon data of dead shells are reported here for the first time. Five calcium carbonate powder samples were submitted to Paleo Labo Co. Ltd., Gifu, Japan, for routine 14 C analysis with a compact accelerator mass spectrometer (AMS; NEC, USA). Other calcium carbonate powder samples were reacted with phosphoric acid, and the carbon dioxide produced was converted to graphite by a Fe catalyst at 620°C for 6-12 hr. The target graphite samples were measured with a Tandetron AMS (High Voltage Engineering Europa, the Netherlands) installed at the Institute for Space-Earth Environmental Research, Nagoya University.
The 14 C measurements are expressed as the fraction of modern 14 C, F 14 C, calculated with the following equation: where A SN and A ON are the 14 C/ 12 C ratios of the sample and of the NIST oxalic acid standard (HOxII), respectively, normalized to δ 13 C PDB = −25‰. The error of the F 14 C analysis by Paleo Labo was better than 0.0026 (1σ; see Online Supplementary Material). The error of the F 14 C analysis at Nagoya University was better than 0.0040 (1σ), except for two samples (0.0070 and 0.0068) with low beam intensities (see Online Supplementary Material).
The time averaging of the samples drilled from the live-caught specimens was 1-7 years, but in half of the samples, high-resolution sampling of <2 years was achieved (see Online Supplementary Material). We previously confirmed that 14 C of the last increment of livecaught individuals represents that of dissolved inorganic carbon in the seawater (for details, see Kubota et al. 2018a). For samples that died in the pre-bomb era, conventional 14 C ages were converted to calendar ages before 1950 CE by using the Marine20 calibration curve (Heaton et al. 2020) and the local marine 14 C reservoir age of −49 ± 20 14 C yr, which was calculated by using 7 14 C measurements of the pre-bomb era portion of live-caught individuals (Alves et al. 2018) (see Online Supplementary Material).

RESULTS
The present study added 4 new 14 C data points to the bomb-14 C record for Funakoshi Bay established by Kubota et al. (2018a). The expanded record improves the time resolution of the period from 1953 to 1963. We also merged radiocarbon data (6 data points) from the live-caught M. stimpsoni shell from Otsuchi Bay (specimen 120901OTW_alive), which is immediately south of Funakoshi Bay (Figure 1), with the bomb-14 C record for Funakoshi Bay, because there is no discernible difference in seawater and shell F 14 C values between the two bays (for more details, see Kubota et al. 2018a). One likely reason of similar 14 C content irrespective of river water input is that freshwater run away at the very surface (0-3 m) of the Otsuchi Bay. (e.g., Tanaka et al. 2017). In the refined bomb-14 C record, the last appearance of pre-bomb 14 C is in 1955 (in the previous record, it is 1952).
In the refined bomb-14 C record peak atmospheric value (e.g., Konishi et al. 1981;Fallon and Guilderson 2008;Glynn et al. 2013;Andrews et al. 2016;Hirabayashi et al. 2017). Given that the estimation error of F 14 C is about 0.0030, the bomb-14 C peak is strong enough to be used for age determination of the dead shells.
We found that the correlation between an average of standardized shell growth indexes (SGI) of the 6 live-caught individuals and the SGIs of 10 dead shells was statistically significant in all combinations (r = 0.57-0.92, p < 0.01; see Online Supplementary Material). We calculated these correlations for the records covering the years 1960-1985, rather than for the entire period of the sclerochronological record, because only during these years were statistically significant correlations found among the 6 live-caught individuals by Kubota et al. (2018a).  Mass Mortality of Long-Lived Bivalve 1637 The widths of the last annual increment can be used to constrain the timing of death with seasonal resolution. The last annual increments of 9 of the 10 individuals described above were as wide as the annual increments of the previous year. Therefore, the possibility that they died (and shell growth ceased) during the summer is low; that is, they died not in summer 2010, but later, after shell growth had ceased. Because the growth period of M. stimpsoni in Funakoshi Bay begins in June and continues until the following January, it is likely that these individuals precipitated their last shell increments in 2010, but died in early 2011. Thus, the March 2011 tsunami likely led to the deaths of these individuals. However, radiocarbon dating can only suggest that the specimens died around 2010 because, given the various uncertainties, it is quite difficult to obtain ages with the necessary precision of one year to identify the cause of death. These uncertainties arise from the analytical uncertainty of the radiocarbon measurement, the slight differences in growth period (month) among individuals, the heterogeneity of seawater 14 C, and the relatively large period of averaging of inner part of the shell required by the sampling technique.
However, by combining bomb-14 C dating with a sclerochronological analysis, age determination with a one-year precision is possible (we selected only relatively long-lived individuals for sclerochronological analysis; thus, we have not been able to determine ages of all of the post-bomb specimens). To sum up the sclerochronological analysis, we found that among 14 C-dated 10 individuals 9 died in 2010, but 1 (specimen 110909FKW_dead_ A02R) died in 2015.
The F 14 C value of the last increment of 1 individual caught alive in September 2010 (100906FKW_alive) was 1.0344, which is in good agreement with the F 14 C values of 7 of 9 dead shells (1.0304-1.0411) (Figure 2). The F 14 C values of the other 2 dead shells, however, were only 1.0157 and 1.0163 (Figure 2). If only the F 14 C data are taken into account, it would seem that these 2 dead shells precipitated their last annual increment after 2010, because they were collected in 2016 from the seafloor and because there F 14 C values are closer to those of specimens caught alive in 2012 and 2016 (120901OTW_alive and 160608FKW_alive_01R, respectively). However, our sclerochronological analysis revealed that these 2 dead shells also precipitated their last annual increment in 2010 (r = 0.81-0.82, p < 0.0001; see Online Supplementary Material). The F 14 C values of the tenth individual inferred by radiocarbon dating alone to have died in 2010 (160608FKW_dead_14R) were 1.0195 for the last annual increment and 1.1035 for increments precipitated when it was ontogenetically younger. However, sclerochronological analysis of this specimen revealed that the last annual increment was precipitated in 2015 (r = 0.82, p < 0.0001) (see Online Supplementary Material). It is noteworthy that, different from the other 9 individuals inferred to have died in 2010, the last increment of this specimen was only half as wide as the increment of the previous year. Therefore, the shell growth in this individual was interrupted in summer 2015, when the individual died for an unknown reason (the summertime temperature in 2015 was not unusual). In summary, on the basis of 14 C data alone, 10 individuals appeared to have died in 2010, but by combining 14 C data with a sclerochronological analysis, it was possible to show that 1 of the 10 died not in 2010 but in 2015.
From these results, we inferred that the deaths of the 9 individuals (dead shells) that precipitated their last annual increments in 2010 were likely related to the March 2011 tsunami for the following reasons: (1) the F 14 C values of the last increments and of the inner parts of the dead shells coincided well with the bomb-14 C record in the study area; (2) correlations between the average SGI of the live-caught individuals and the SGIs of the dead shells in all combinations were statistically significant; and (3) the width of the last increments showed that these animals died, not during the shell growth season, but later, after growth had ceased. We showed that mass mortality in M. stimpsoni (9 of 20 individuals that died in the post-bomb era) was likely related to the March 2011 tsunami, but other possible reasons for this mass mortality need to be carefully addressed. In particular, during winter 2010 to spring 2011, water temperature may have become unsuitable or predation by other animals may have increased. With regard to the former, the mean bottom water temperature in March 2011 was 7.05°C, 0.97°C higher than the 30-year climatological mean; therefore, unsuitably high water temperatures probably did not cause death in M. stimpsoni, because this clam species can survive and grow (calcify) in water temperatures of up to ∼20°C in summer. This species is considered to have a wide temperature tolerance because it is distributed over a wide latitudinal range in the northwest Pacific, from temperate to subpolar regions (e.g., off South Korea, Khim et al. 1998;off Hokkaido, northern Japan, Tanabe et al. 2017). This wide temperature tolerance suggests that it is unlikely that unsuitable temperatures were a major cause of death of this clam. With regard to the latter, it is also unlikely that predation by other animals caused mass mortality in this clam for three reasons: (1) M. stimpsoni lives beneath the seafloor and is thus inaccessible to many predators; (2) the thick shells of adult individuals protect them from predation; and (3) the shells did not show any damage, such as holes drilled by carnivorous gastropods, or other signs of predation (Nakaoka 2000;Chiba and Sato 2013).
The distribution of tsunami deposits in Funakoshi Bay suggests that the seafloor was first greatly eroded by sub-bottom flows caused by the tsunami and subsequently covered by laminated tsunami deposits . Clams that became buried beneath these tsunami deposits might have suffocated because they were unable to escape from deep burial. This possibility is supported by the finding of an articulated dead shell within the tsunami deposits; furthermore, the F 14 C value of the last annual increment of this shell (1.0392) was close to those of specimens caught live during 2010-2016 and also to those of the clams that died in March 2011. However, sclerochronological analysis of this shell was not performed because the shell became broken during the coring (Figure 2). It is also possible that sediment liquefaction caused by ground shaking during the earthquake caused M. stimpsoni individuals to become exposed above the seafloor, because density of this animal is lower than that of the ambient substrate (Seike et al. 2019). Such exposure might lead to the death of individuals that were not able to re-enter the sediment (e.g., by predation by octopus). Although it is clear that the mass mortality of M. stimpsoni was related to March 2011 earthquake and tsunami, the specific mechanism of death is difficult to determine. Because we collected dead shells only from the seafloor surface, we do not know whether many dead individuals might remain buried in the tsunami deposits, although wave action associated with, for example, an atmospheric low pressure system might bring some once-buried individuals to the seafloor.

Relationship between M. stimpsoni Deaths and Past Tsunamis
Although the deaths of at least 9 of the 27 studied individuals collected as dead shells were related to the March 2011 earthquake and tsunami, the other deaths did not coincide with Mass Mortality of Long-Lived Bivalve 1639 the occurrence of a tsunami. Among those clams that died during the post-bomb era, 14 C dating and sclerochronological analysis revealed that 2 individuals died in 1966and 4 others died in 1968, 1971, 1986, and 2015). Because neither water temperature (high or low) nor the shell growth period (longer or shorter) were unusual in these years, no environmental factor seems to have directly influenced their deaths. Furthermore, it is noteworthy that we found no evidence that the tsunami generated by the 1960 Chilean earthquake that struck northeast Japan in 1960 caused mass mortality of M. stimpsoni (Figure 4). One possible reason is that the amplitude of the 1960 tsunami (maximum runup height along the Japanese coast was ∼6 m) was not large enough to disturb the seafloor environment and damage the benthic community. It would be worthwhile to test this hypothesis by collecting and dating additional dead shells.
Huge tsunamis strike coastal northeastern Japan at intervals of 50-100 years, so it is possible that past tsunamis similar in scale to the March 2011 tsunami might have caused mass mortality of M. stimpsoni. One huge tsunami occurred in 1933 (Showa Sanriku earthquake: run-up height, ∼29 m) and another in 1896 (Meiji Sanriku earthquake: run-up height, ∼38 m). However, it is not possible to link a tsunami occurring during the pre-bomb era with the time of death of M. stimpsoni individuals by using only conventional 14 C dating. The time of death estimated by conventional 14 C dating of the 7 individuals that died during the pre-bomb era overlap with the occurrence times of these two tsunamis (Figure 4), but to confirm this association, it will be necessary to extend the master chronology further back in time.
The oldest clam among the live-caught individuals was 92 years old ), but one of the dead shells was from an individual that had lived for 135 years (110909FKW_dead_B01R). Therefore, there is a chance that additional long-lived individuals will be found by future field surveys. Intriguingly, the longer-lived individual (and many others) must have survived two huge tsunamis, the Showa Sanriku earthquake tsunami in 1933 and the Meiji Sanriku earthquake tsunami in 1896, only to be killed by the March 2011 tsunami. Without the huge environmental disturbance caused by the tsunami, this animal might have lived even longer. Our findings imply that tsunamis act as a constraint on the longevity of M. stimpsoni clams living in this region. It would be worthwhile to test this idea by comparing the longevity of M. stimpsoni living in this region with that of this species living in other regions that rarely experience tsunamis.

CONCLUSIONS
By combining 14 C and sclerochronological analyses, we showed that many individuals of the long-lived bivalve M. stimpsoni in Funakoshi Bay died during the huge tsunami that took place in March 2011. Further extension of crossdated annual growth increment records should make it possible to identify past mass mortality events related to tsunamis. In particular, the association of M. stimpsoni deaths with the tsunamis that occurred in 1933 and 1896 should be investigated by the method used here for the March 2011 tsunami. The 1896 tsunami occurred in June, which is during the M. stimpsoni shell growth season; thus, evidence for the tsunami may be recorded in the shell by trace element variations, as Murakami-Sugihara et al. (2019) demonstrated in mussel shells. Trace elements in the shell of an individual that experienced the 1896 tsunami as a juvenile, when its annual shell growth rate was large, might serve as a distinct time marker and chronological constraint.
This dating method based on bomb-14 C and sclerochronological analyses of the shells of a long-lived bivalve presented in this study can be applied to the shells of other mollusks that inhabit the Pacific coastal waters of Japan or other oceanic regions affected by tsunamis. However, application of this method to bivalves living in regions where the bomb-14 C peak is dampened (e.g., upwelling and high-latitude ocean regions) would not be effective, nor is it applicable to bivalves that live at abyssal depths, where the amount of bomb-induced 14 C Figure 4 (A) Bomb-14 C record of the study area as in Figure 2 (black line with circles) (Kubota et al. 2018a) and the F 14 C values of the last increment of post-bomb era dead shells collected from the seafloor of Funakoshi Bay (red squares). Black arrows show the dates of major tsunamis that have struck northeastern Japan during the postbomb era. (B) Block plot of the 1σ (68%, thick bars) and 2σ (95%, thin bars) calibrated age ranges for pre-bomb era dead shells collected from Funakoshi Bay, as determined by conventional radiocarbon dating. Black arrows at the top of the panel indicate the dates of major tsunamis that struck northeastern Japan during the pre-bomb era. It is difficult to associate any tsunami with the deaths of M. stimpsoni that died during the pre-bomb era because of the relatively large dating errors.
Mass Mortality of Long-Lived Bivalve 1641 is too small. Nevertheless, this method may prove effective for understanding the responses of shallow-water bivalve mollusks to environmental disturbances, especially in countries and remote regions where observations of past tsunamis are scarce. Further understanding of such responses may lead to better protection and sustainable use of fishery resources, especially commercially important shellfish.