Earthquake Tsunami Research Articles

Data‐driven machine learning for multi‐hazard fragility surfaces in seismic resilience analysis

AbstractOffshore earthquakes and subsequent tsunamis pose significant risks to many coastal populations worldwide. This paper introduces a data‐driven machine learning model that synthesizes accurate 3D earthquake–tsunami fragility surfaces from randomly selected 2D fragility curves. The integration of physics‐based simulations enhances the model's reliability for these specific hazards, making it a valuable tool for multi‐hazard analysis in earthquake–tsunami contexts. Additionally, by shifting 2D fragility curves to represent retrofitted structural systems, the model can generate earthquake–tsunami fragility surfaces for community‐level mitigation studies. While the model is demonstrated for earthquake–tsunami scenarios, its methodology architecture has the potential to contribute to other multi‐hazard situations for the initial conditions in multi‐hazard community resilience analysis.

Not trench-parallel but trench-normal source fault of the 1994 Hokkaido Toho-oki earthquake as revealed by the aftershock relocation using HypoDD

In the southwestern area of the Kurile Islands subduction zone, three large earthquakes including the Mw8.2 1969 and the Mw8.3 1994 Hokkaido Toho-oki earthquakes and the Mw7.8 1975 tsunami earthquake occurred in almost the same location. The three main shocks and their aftershocks were re-determined simultaneously using a double-difference earthquake location method, HypoDD. As a result, the 1969 and the 1975 main shocks and their aftershocks were concentrated near the upper surface of the subducting Pacific plate, indicating that both events were plate boundary megathrust earthquakes as shown by previous studies. On the other hand, the 1994 main shock and its aftershocks were distributed within the Pacific plate indicating that the 1994 event was an intraslab event as reported by previous studies. The result that interplate earthquakes and intraslab earthquakes could be appropriately distinguished is strong evidence that the accuracy of hypocenter determination is sufficiently high. The centroid moment tensor solution of the 1994 main shock has two nodal planes which are parallel or normal to the Kurile Trench. According to the high-resolution hypocenter distribution obtained in the present study, the aftershocks appear to locate on the nodal plane normal to the trench, suggesting that the source fault ruptured by the 1994 main shock was not parallel but normal to the Kurile Trench. Moreover, we found that the Coulomb failure stress change (ΔCFS) produced by the trench-normal fault plane is more likely to trigger the largest aftershock than ΔCFS produced by the trench-parallel fault plane.Graphical abstract

Modeling ground motions and crustal deformation from tsunami earthquakes: Rupture parameter constraints from the 2010 Mentawai event

Modeling ground motions and crustal deformation from tsunami earthquakes: Rupture parameter constraints from the 2010 Mentawai event

Dataset of Post-Event Survey of the 2024 Noto Peninsula Earthquake Tsunami in Japan

An earthquake with a moment magnitude of 7.5 (Mw) struck the northern Noto Peninsula, Ishikawa Prefecture, Japan, at 16:10 local time on January 1, 2024. This earthquake triggered a tsunami that propagated along the coastline of Ishikawa, Toyama, and Niigata Prefectures facing the Sea of Japan and significantly damaged coastal communities and infrastructure. Approximately 70 researchers from 23 universities or other institutes throughout Japan formed a joint research group to conduct a post-tsunami survey along a 340 km stretch of the coast. Based on the watermarks and traces of the tsunami, the inundation and run-up heights were surveyed using total stations, automatic optical levels, laser range finders, and a real-time kinematic (RTK) Global Navigation Satellite System (GNSS). The tidal correction was adjusted using astronomical tidal tables. In total, 303 survey records have been compiled, generating the NP2024TS (Noto Peninsula 2024 Tsunami Survey) dataset. This dataset provides comprehensive information on the inundation and run-up heights of the tsunami, which is useful for understanding tsunami characteristics and validating numerical tsunami models.

The role of vehicular considerations in tsunami shelter location optimization

Using vehicles during tsunami evacuation might exacerbate road congestion, potentially increasing casualties. However, since the 2011 Tohoku Earthquake Tsunami, vehicular evacuation during tsunamis has partially been permitted, especially to assist elderly and disabled individuals with mobility challenges. The present study aimed to elucidate the impact of vehicular use on reducing mortality rate by introducing a new tsunami evacuation shelter in a given coastal area. To do so, the authors applied an agent-based tsunami evacuation simulation model that can account for behaviors of both pedestrians and vehicular evacuees to three Japanese coastal cities (Mihama Town, Shingu City, and Minabe Town) considering future Nankai–Tonankai Earthquake Tsunami scenarios. The authors first determined mortality rates under scenarios both with and without vehicular evacuation using current shelter placements. Subsequently, mortality rates with the introduction of a new tsunami evacuation shelter were examined for both scenarios. The findings showed that effective shelter locations for minimizing mortality rates varied depending on whether vehicular evacuation was considered. It was also found that the effectiveness of vehicular use during evacuation and the strategic placement of new tsunami evacuation shelters are highly dependent on local factors such as population distribution, existing shelter locations, and road layouts. These findings underscore the necessity for area-specific tsunami evacuation simulations to construct robust evacuation strategies.

Adjoint Synthesis for Trans‐Oceanic Tsunami Waveforms and Simultaneous Inversion of Fault Geometry and Slip Distribution

AbstractTsunamis propagate over long distances and can cause widespread damage even after crossing ocean basins. Prediction of tsunamis in distant areas based on observations near their sources is critical to mitigating damage. In recent years, the accuracy of numerical models of trans‐oceanic tsunami propagation has improved significantly due to the incorporation of effects such as the solid earth response to tsunami loading and wave dispersion. However, these models are computationally expensive and have not been fully utilized for real‐time prediction. Here, we derive the adjoint operator for the linear set of equations describing deep‐ocean tsunami propagation and show how a pre‐computed database of adjoint states can achieve rapid synthesis of tsunami waveforms at target sites from nonpoint arbitrary tsunami sources. The adjoint synthesis method allows for an exhaustive parameter search for tsunami source estimation. A method for simultaneous inversion of fault geometry and slip distribution using adjoint synthesis with Sequential Monte Carlo method was proposed and applied to the 2012 Haida Gwaii earthquake tsunami. The influence of model accuracy and the amount of observed data on the estimation of tsunami sources and waveforms was examined. It was found that with a highly accurate propagation model, using only a limited amount of observed data produced source and waveform estimates very similar to the final models obtained with much larger data sets. The final inferred fault model involved megathrust slip distributed between the Haida Gwaii trench and the Queen Charlotte fault. The proposed method can also quantify the uncertainty of the waveform forecasts.

New source model for the 1771 Meiwa tsunami along the southern Ryukyu Trench inferred from high-resolution tsunami calculation

The 1771 Meiwa tsunami which struck the southern Ryukyu Islands (Sakishima Islands) had greater than 22 m run-up height, leaving about 12,000 casualties in its wake. At many places, the tsunami inundation or lack of inundation is well recorded in historical documents. Several tsunami source models have been proposed for this event using historical records as constraints of tsunami calculations. Nevertheless, the source model remains under discussion. This study re-evaluated the tsunami wave source model of the 1771 Meiwa tsunami using high-resolution (10 m mesh) bathymetric and topographical data for tsunami calculation, the latest historical record dataset, and seismological knowledge. Results demonstrated that a tsunami earthquake along the southern Ryukyu Trench was the likely cause of the 1771 event. However, it is noteworthy that assumption of a large slip with 30 m is necessary for a shallow and narrow region (fault depth = 5 km, fault width = 30 km, Mw = 8.49) of the plate boundary in the Ryukyu Trench, which is far larger than previously thought. This requirement of very large initial water level change at the source might involve not only the fault rupture along the plate boundary but also deformation by splay faults, inelastic deformation of unconsolidated sediments near the trench axis, and/or giant submarine landslides. Results also show that the effects of fault parameters on the run-up were quite different depending on the offshore coral reef width. This phenomenon strongly constrained the fault width to 30 km. Our tsunami ray tracing analysis further revealed the effects of bathymetry on tsunami propagation. It is noteworthy that meter-long huge tsunami boulders tend to be distributed along the specific coasts at which the tsunami was concentrated by bathymetric effects. This finding suggests that past tsunamis, including the 1771 event, might have affected the specific coral reefs on Sakishima Islands repeatedly, which is crucially important for understanding the heterogeneous distribution of tsunami boulders. This feature might also be useful to elucidate the effects of large tsunamis on the corals and reefs because a direct comparison of coral reefs that are damaged and not damaged by tsunami waves is testable in narrow areas in the case of the Sakishima Islands.

Assessing the Risk of Potential Tsunamigenic Earthquakes in the Mentawai Region by Seismic Imaging, Central Sumatra

AbstractIn the marginal regions along subduction zones, oceanic plates subducting beneath continental plates produce the largest number of earthquakes on Earth and sometimes devastating tsunamis. The lateral segmentation of earthquakes along the Sumatra subduction zone is well documented. However, the entirely different seismic behaviors among the segments indicate that local structures are key elements controlling coseismic slip propagation; in particular, frontal accretionary prism structures are closely associated with tsunami generation. Offshore of Central Sumatra, in the Mentawai segment, large earthquakes nucleated in 2007 and 2010 caused many human casualties and a great deal of property loss. Using seismic reflection data, we show the subsurface structure of accretionary over a significant portion of the frontal wedge that did not rupture during the 2007 earthquake, and this area remains locked. The subsurface deformation structure at the wedge front, which is similar to that in the 2010 Mw7.8 tsunami earthquake rupture zone, suggests the potential for a large tsunami earthquake in the near future. On the other hand, the along‐strike variations of the effective basal friction at shallow depth may indicate that different coseismic behaviors are caused by a sudden failure of the deeper seismogenic zone.

Possible Maximum Earthquake Tsunami Along the South Coast of China Inferred From GPS‐Derived Surface Velocities

AbstractThis study aims to assess the tsunami hazard of the maximum possible earthquake along the south coast of China. GPS‐derived surface velocities around the Manila Trench have been used to clarify the active tectonics and estimate the seismic potential of Manila subduction, where catastrophic tsunamis may originate and cause devastating damage to the surrounding area of the South China Sea. Based on the maximum potential earthquake magnitude estimated from the accumulated seismic moment, more than 100,000 possible tsunami scenarios with varying epicenters have been simulated considering the slip heterogeneity. The impact of several uncertain parameters including the release period and release ratio of accumulated seismic moment are investigated, as well as some assumptions in the model of tsunami generation. Although different release conditions for accumulated seismic moment result in a significant impact on the value of tsunami height, similar trend in spatial distribution has been shown. The results show that the tsunami wave height along the south coast of China has an increasing trend from west to east, with the maximum tsunami wave height exceeding 10 m corresponding to the potential earthquake Mw = 9.1 with a 1,000‐year release period. The impact of fault slip heterogeneity on tsunami wave height is significant. Compared to uniform slip, considering stochastic slip will increase the average tsunami wave height by about 15%. Considering locking distribution will increase the average tsunami wave height by about 20%–30% and also increase the maximum tsunami wave height by about 60%–70%.

Submarine landslides and tsunami genesis in Sagami Bay, Japan, caused by the 1923 Great Kanto earthquake

The 1923 Great Kanto earthquake occurred on September 1, in Japan, and caused severe damage mainly in the Kanto region. Tsunamis were observed over wide regions from the east coast of the Izu Peninsula to the west coast of the Boso Peninsula, and particularly, the damage in Atami was devastating. Many earthquake fault models including those of Kanamori (1971) and Ando (1971) were proposed based on the records of land deformation. However, such fault models cannot sufficiently explain the tsunami elevation and its initial sea-level motion on the coasts. Hence, the detailed mechanisms remain elusive. This study examines the possibility that a leading mechanism of the tsunami in the 1923 Great Kanto earthquake was a large-scale submarine landslide that occurred at Sagami Bay and at the mouth of Tokyo Bay based on the records of depth data measured by the Imperial Japanese Navy (1924) before and after the earthquake. We first show that the tsunami calculated by each fault model was inconsistent with the waveform at Yokosuka, the coastal tsunami elevations, and initial sea-level motion. Then, based on statistical analysis of the depth changes at the 1923 Great Kanto earthquake, we found that the seafloor bathymetric changes represented large-scale submarine landslides that may correspond to long-runout submarine liquefied sediment flows. The seafloor gradient over a 40 km flow-out distance was equal to or less than 0.4∘\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$^\\circ$$\\end{document}. Through the identification of the submarine landslide source by tsunami backpropagation analysis and utilizing an analytical solution of a high-density gravity flow and a sensitivity analysis, we conducted a range of numerical simulations of the 1923 Great Kanto earthquake tsunamis using a fault model and a submarine landslide tsunami source model due to a high-density liquefied gravity flow. The results quantitatively accounted for the discrepancy between the observed tsunami records with maximum tsunami elevations over 12 m and the fault-model–based simulations with maximum tsunami elevations of 2 to 5 m and explained consistently the maximum tsunami elevation distributions as well as the time-series tsunami waveforms. These results may thus facilitate and deepen our understanding of the earthquake-induced submarine landslide tsunami risk as cascading multi-geohazards.

Computational Analysis of Tsunami Wave Behaviour for Three Historical Tsunami Events using T-Impulse Model

Natural catastrophes pose a serious threat to both human life and the environment because they are unpredictable. One of the most devastating natural disasters is a tsunami, and forecasting models are essential to preventing catastrophic damage to the environment and people along the coast. In the Impulse model, the generation of a tsunami depends on the impulse force generated during the event. Understanding tsunamis begins with simulating the tsunami generation process. This process involves simulating both the motion of the seafloor and the subsequent motion of the water above for tsunamis caused by underwater earthquakes. This modeling strategy can mimic all three stages of a tsunami: generation, propagation, and run-up. Three separate earthquake tsunami events—the 1755 Lisbon earthquake, the 1964 Alaska earthquake, the 2004 Sumatra earthquake are each investigated in this research. To demonstrate its relevance to current events and various ocean locations, the results of these events are compared and confirmed with the observed data. Analyzing the parameters used in this modeling study and identifying the parameter that has the most influence will demonstrate their significance in tsunami generation. The seabed displacement profile, seawater deformation, changes in tsunami characteristics during propagation, the tsunami’s travel time, earliest arrival time, the tsunami wave height at the coast, and inundation distance are the anticipated findings from this study. The major objective of this study is to obtain the maximum and most accurate result possible using the fewest parameters possible.

Do Subducted Seamounts Act as Weak Asperities?

AbstractThe additional work of ploughing makes seamounts more resistant to subduction and more strongly coupled than smoother areas. Nevertheless, the idea that subducted seamounts are weakly coupled and slip aseismically has become dominant in the last decade. This idea is primarily based on the claim that a seamount being subducted in the southern Japan Trench behaves this way. The key element in this assertion is that large M ∼ 7 earthquakes that abut the leading edge of the seamount require that the seamount be aseismically sliding to initiate them. More recent observations show instead that the surrounding region is aseismically sliding while the seamount acts as a stationary buttress. Here we re‐examine this case and model it with both weak and strong asperity assumptions. Our modeling results show that only a strong asperity model can produce this type of earthquake. Strong asperities also rupture the seamount in great earthquakes with long recurrence times. This provides the previously unknown source for a series of great tsunami earthquakes that have occurred along the southern Japan Trench, the most recent being the 1677 M8.3–8.6 Enpo Boso‐oki tsunami earthquake. The “weak asperity” hypothesis is thus found to be false in this foundational example.

Multihazard capacity optimization of an NPP using a multi-objective genetic algorithm and sampling-based PSA

After the Tohoku earthquake and tsunami (Japan, 2011), regulatory efforts to mitigate external hazards have increased both the safety requirements and the total capital cost of nuclear power plants (NPPs). In these circumstances, identifying not only disaster robustness but also cost-effective capacity setting of NPPs has become one of the most important tasks for the nuclear power industry. A few studies have been performed to relocate the seismic capacity of NPPs, yet the effects of multiple hazards have not been accounted for in NPP capacity optimization. The major challenges in extending this problem to the multihazard dimension are (1) the high computational costs for both multihazard risk quantification and system-level optimization and (2) the lack of capital cost databases of NPPs. To resolve these issues, this paper proposes an effective method that identifies the optimal multihazard capacity of NPPs using a multi-objective genetic algorithm and the two-stage direct quantification of fault trees using Monte Carlo simulation method, called the two-stage DQFM. Also, a capacity-based indirect capital cost measure is proposed. Such a proposed method enables NPP to achieve safety and cost-effectiveness against multi-hazard simultaneously within the computationally efficient platform. The proposed multihazard capacity optimization framework is demonstrated and tested with an earthquake–tsunami example.

Seamount and Ridge Subduction at the Java Margin, Indonesia: Effects on Structural Geology and Seismogenesis

AbstractThe Java ‐ Lesser Sunda margin, which features multi‐scale subducting oceanic basement relief, is classified as neutral (Lombok and Sumbawa) to erosional (Central Java to Bali) in comparison to its accretionary counterpart offshore Sumatra. However, a comprehensive analysis of how plate boundary and upper plate structure across the neutral to erosional transition are modulated by the subduction of oceanic basement relief is lacking to date. To shed light on the tectonic parameters that push the margin into the neutral or erosional domain, we combine multi‐channel reflection seismic images derived through a grid‐based P‐wave velocity inversion, and multibeam bathymetric maps. The data document how different scales of subducting topography modify seafloor morphology, upper plate structure, and décollement position. Large‐scale subducting features cause a landward shift of the deformation front, shortening of the accretionary wedge, and seafloor steepening at the relief's trailing edge. Small‐scale subducting ridges primarily impact the frontal prism resulting in over‐steepening at the trench and localized slope failure. Ahead of subducting relief, deformation of the accretionary wedge encompasses enhanced compression and a reduction in seafloor slope but appears independent of the size of the relief. Ridge and seamount subduction induce frontal erosion and basal erosion offshore Lombok and Bali, respectively. Our P‐wave velocity models indicate that the rigidity of the upper plate's base along the eastern Sunda margin is significantly lower than the worldwide trend. We conclude that this favors the genesis of tsunami earthquakes that have occurred on the Java margin.

Study on the dynamic response of offshore bridge under earthquake action and tsunami impact

This paper takes an offshore continuous rigid bridge as the background bridge to reveal the dynamic response law of an offshore bridge under the action of earthquake and tsunami. A tsunami wave is simulated by a dam-break wave. A numerical model of a dam-break impact of a bridge pier is established in Fluent software, and k−ε turbulence model is chosen to solve the Navier–Stokes equation. Meanwhile, ABAQUS software is used to establish a three-dimensional nonlinear full bridge model. The dam-break wave working condition is proposed according to the Stoker theory. The bridge dynamic response index under a single earthquake action and a continuous earthquake tsunami action is calculated, and the dynamic response rules under the two actions are summarized. The results show that compared with the tsunami action alone, the sequential action of the earthquake and tsunami will cause a greater dynamic response of the bridge structure. With the increase in the height of the incoming tsunami wave, the waveform trend difference of the response index will be more significant. Moreover, the influence coefficient of the dynamic response index of the bridge tends to increase as the probability level of the earthquake occurrence decreases.

Prediction of shelter needs based on simulation results of evacuees choosing evacuation facilities due to a potential tsunami earthquake in Padang, Indonesia

Tsunami disasters are predicted to occur in the West Sumatra region due to earthquakes. To deal with the tsunami hazard that will hit Padang City, the government has built several shelters and evacuation routes. At present, the number of evacuation facilities is still limited, this will cause community vulnerability in facing tsunami hazards. To overcome this vulnerable condition, it is necessary to determine the number of shelters needed to accommodate evacuees. In this study, the determination of the required number of shelters is done by simulating the number of evacuees who choose shelters during evacuation. The model of evacuee generation from 5 tsunami-prone zones in Padang City was created using a linear regression method. The model for selecting the type of evacuation facility in each zone uses a binary logit model. Prediction of shelter needs is carried out for minimum, normal and maximum conditions that may occur. From the simulation results, the number of evacuees who choose the shelter as a means of evacuation is obtained. The number of shelters needed in Padang City is a minimum of 6 units, normal 23 units and maximum 50 units.

The natural warning signs of tsunami earthquake in Indonesia: case of the 2006 Cilacap event

ABSTRACT Natural warning signs and self-evacuation are essential to save lives in the event of a tsunami caused by an earthquake, especially if there are no clear early signs for example ground shaking. This paper explores the participants’ experiences focusing on how participants’ perceptions and responses to what they felt, saw or heard before seeing the first wave of 2006 tsunami in Cilacap, Indonesia. Using a phenomenological approach, this paper analyzes the stories of ten participants in three sub-districts, namely Cilacap Selatan, covering islands of Nusakambangan, Adipala, and Binangun. The tsunami was caused by a slow-motion earthquake, so people did not feel the ground shaking. The participants did not recognise natural signs of a tsunami, so they ignored them. Before seeing the first wave, the witness heard sound of explosions, extreme low tide, unusual animal behavior around, the appearance of a large number of marine products – lobsters (Crustacea) and sea shells (Mytilus spp.), and unusual color and shape of the waves. These signs should be socialised to the community in tsunami-prone areas, thereby increasing community’s disaster mitigation capacity. This study presents some natural warning signs that have not been reported by previous studies.

Reappraisal of tsunami hazard for the Northern Coastal of Egypt considering sea level rise and delta subsidence scenarios

In this research, we have compiled a catalogue of the worst tsunami earthquakes that caused destruction and losses in the northern coastal region of Egypt. Variability in seismic source parameters reflects the level of cognitive uncertainty that exists because most of these destructive tsunamis, if not all, occurred in the pre-instrumental period. We also synthesized a set of hypothetical tsunamis constrained by seismotectonic knowledge of the area of influence to embrace the uncertainty. We also collected the topo-bathymetry information available to make a simulation on the scale of the Egyptian north coast. The multiple scenarios were computed considering references and hypothetical scenarios, and results were presented in maps and tables for wave heights, inundation depth, and estimated arrivals. Comparisons were made between the results calculated in this study and those from previous studies. These were done to identify and measure the reasons for the difference and similarities. The results may contribute to developing risk reduction strategies and a national early warning system. Variability in tsunami intensity measures is prepared for the hypothetical scenarios due to the difference in the fault’s magnitude, location, and geometry. Also, for the first time, temporally based scenarios the sea-level rise and delta sinking were considered in order to highlight the importance of multi-risk analysis for the region of interest. The insight on incorporating other hazards of the temporal resolution indicates a significant increase in inundated areas and wave height. Therefore, proper and detailed estimation of multi-hazard maps for the coastal region of the Nile Delta is highly required to protect the ongoing sustainable development and protect the people and their assets.

Do upper-plate material properties or fault frictional properties play more important roles in tsunami earthquake characteristics?

Tsunami earthquakes are a type of shallow subduction zone events that rupture slowly (<1.5 km/s) with exceptionally long duration and depleted high frequency radiation, resulting in a large discrepancy of Mw and Ms magnitudes and abnormally large tsunami along coastal areas. Heterogeneous fault frictional properties at shallow depth have been thought to dominate tsunami earthquake generation. Some recent studies propose heterogeneous upper-plate material properties determine rupture behavior of megathrust earthquakes, including characteristics of tsunami earthquakes. In this study, we use a recently developed 3D dynamic earthquake simulator to explore tsunami earthquake generation and systematically examine roles of upper-plate material properties and fault frictional properties in tsunami earthquake characteristics in a physics-based framework. For heterogeneous fault friction, we consider isolated asperities with strongly velocity-weakening properties embedded in a conditionally stable zone with weakly velocity-weakening properties. For heterogeneous upper-plate properties, we consider a generic depth profile of seismic velocity and rigidity constrained from seismic surveys. We design a set of models to explore their effects on tsunami earthquake generation and characteristics. We find that the conditionally stable zone can significantly slow down rupture speeds of earthquakes that nucleate on asperities to be <1.5 km/s over a large depth range (1–20 km), while heterogeneous upper-plate properties can only reduce rupture speeds to be ∼1.5–2.0 km/s over a narrow depth range (1-3 km). Nevertheless, heterogeneous upper-plate properties promote cascading rupture over multiple isolated asperities on the shallow subduction plane, contributing to large tsunami earthquake generation. We also find that fault frictional properties play much more important roles than upper-plate material properties in long normalized duration, high-frequency depletion and low moment-scaled radiated energy in tsunami earthquakes. In addition, the effective normal stress on the subduction plane, which affects fault frictional strength, also influences the characteristics of tsunami earthquakes, including duration, stress drop and moment-scaled radiated energy.

Cover Image, Volume 38, Issue 2

Computer-Aided Civil and Infrastructure EngineeringVolume 38, Issue 3 COVERFree Access Cover Image, Volume 38, Issue 2 First published: 17 January 2023 https://doi.org/10.1111/mice.12975AboutPDF ToolsRequest permissionExport citationAdd to favoritesTrack citation ShareShare Give accessShare full text accessShare full-text accessPlease review our Terms and Conditions of Use and check box below to share full-text version of article.I have read and accept the Wiley Online Library Terms and Conditions of UseShareable LinkUse the link below to share a full-text version of this article with your friends and colleagues. Learn more.Copy URL Share a linkShare onFacebookTwitterLinkedInRedditWechat Graphical Abstract On the cover: The cover image is based on the Research Article A computational framework for fluid–structure interaction with applications on stability evaluation of breakwater under combined tsunami—earthquake activity by Shuai Huang and Chuanzheng Liu, https://doi.org/10.1111/mice.12880. Volume38, Issue3February 2023 RelatedInformation