Abstract:We demonstrate the efficacy of a Bayesian statistical inversion framework for reconstructing the likely characteristics of large pre-instrumentation earthquakes from historical records of tsunami observations. Our framework is designed and implemented for the estimation of the location and magnitude of seismic events from anecdotal accounts of tsunamis including shoreline wave arrival times, heights, and inundation lengths over a variety of spatially separated observation locations. As an initial test case we use our framework to reconstruct the great 1852 earthquake and tsunami of eastern Indonesia. Relying on the assumption that these observations were produced by a subducting thrust event, the posterior distribution indicates that the observables were the result of a massive mega-thrust event with magnitude near 8.8 Mw and a likely rupture zone in the north-eastern Banda arc. The distribution of predicted epicentral locations overlaps with the largest major seismic gap in the region as indicated by instrumentally recorded seismic events. These results provide a geologic and seismic context for hazard risk assessment in coastal communities experiencing growing population and urbanization in Indonesia. In addition, the methodology demonstrated here highlights the potential for applying a Bayesian approach to enhance understanding of the seismic history of other subduction zones around the world.

What problem does this paper attempt to address?

The problem this paper attempts to address is reconstructing the characteristics of large earthquakes from tsunami observation data in historical records, particularly those that occurred before instrumental recordings. Specifically, the paper demonstrates how to use a Bayesian statistical inversion framework to extract key information about earthquake parameters such as location and magnitude from historical records. Using the example of the massive earthquake and tsunami that occurred in the Banda Arc region of Indonesia in 1852, the research team inferred from historical records of tsunami arrival times, wave heights, and inundation extents that the earthquake had a magnitude of approximately 8.8 Mw and that the rupture area was likely located in the northeastern Banda Arc. ### Main Issues Addressed in the Paper: 1. **How to reconstruct the characteristics of large earthquakes using tsunami observation data from historical records?** - The research team developed a Bayesian statistical inversion framework to extract key earthquake parameters from historical records. - This framework can handle uncertainties and imprecise historical data, thereby improving the accuracy of earthquake characteristic estimates. 2. **What are the specific characteristics of the 1852 Banda Arc earthquake?** - By using tsunami observation data from historical records, the research team inferred that the 1852 Banda Arc earthquake had a magnitude of approximately 8.8 Mw. - The rupture area of the earthquake was likely located in the northeastern Banda Arc, consistent with the region of seismic activity recorded by modern instruments. 3. **How to verify whether the tsunamis in historical records were caused by large megathrust earthquakes?** - The research team validated that the 1852 tsunami was more likely caused by a large megathrust earthquake rather than a submarine landslide by comparing model results under different hypotheses. - They used Bayesian methods to quantify the probability distributions under different hypotheses, leading to the most likely earthquake source. ### Solutions: - **Bayesian Statistical Inversion Framework**: This framework combines tsunami observation data from historical records with modern seismological knowledge to estimate earthquake parameters through probability distributions. - **Forward Modeling**: Using the Geoclaw software package to simulate tsunami wave propagation, combined with observation data from historical records, to validate the model's accuracy. - **Parameter Calibration**: By reducing the dimensionality of the parameter space, computational efficiency is improved, ensuring the model's operability. ### Significance: - **Improving the Accuracy of Earthquake Risk Assessment**: By utilizing data from historical records, a more comprehensive assessment of earthquake risks can be achieved, especially in regions lacking modern instrumental records. - **Promoting the Application of Bayesian Methods**: This study demonstrates the potential of Bayesian methods in seismology, providing a reference for similar studies in other regions. Overall, this paper successfully reconstructs the characteristics of the 1852 Banda Arc earthquake from historical records through systematic methods and technical means, providing important scientific evidence for earthquake risk assessment and disaster mitigation.

Methodological reconstruction of historical seismic events from anecdotal accounts of destructive tsunamis: a case study for the great 1852 Banda arc mega-thrust earthquake and tsunami

Methodological Reconstruction of Historical Landslide Tsunamis Using Bayesian Inference

A tale of two faults: Statistical reconstruction of the 1820 Flores Sea earthquake using tsunami observations alone

The 2016 Kaikōura Earthquake: Simultaneous Rupture of the Subduction Interface and Overlying Faults

Tsunami scenario triggered by a submarine landslide offshore of northern Sumatra Island and its hazard assessment

Testing Megathrust Rupture Models Using Tsunami Deposits

Bayesian Hierarchical Modeling for Probabilistic Estimation of Tsunami Amplitude From Far‐Field Earthquake Sources

Validation of Linearity Assumptions for Using Tsunami Waveforms in Joint Inversion of Kinematic Rupture Models: Application to the 2010 Mentawai M-W 7.8 Tsunami Earthquake

A Medium‐Sized Paleo‐Tsunami Reconstruction by a Deep Neural Network Processing Sedimentary Deposits

Cascading Geological Hazards and Risks of the 2018 Sulawesi Indonesia Earthquake and Sensitivity Analysis of Tsunami Inundation Simulations

Historical and geological evidence for the seventeenth-century tsunamis along Kuril and Japan trenches: implications for the origin of the AD 1611 Keicho earthquake and tsunami, and for the probable future risk potential

A Workflow to Integrate Numerical Simulation, Machine Learning Regression and Bayesian Inversion for Induced Seismicity Study: Principles and a Case Study

Modeling the 2024 Noto Peninsula earthquake tsunami: implications for tsunami sources in the eastern margin of the Japan Sea

Source Parameters of the Great Sumatran Megathrust Earthquakes of 1797 and 1833 Inferred from Coral Microatolls

Use of tsunami observations for resolving coseismic slip of recent and historic large submarine earthquakes

Probabilistic earthquake-tsunami financial risk evaluation for the District of Tofino, British Columbia, Canada

Understanding paleo-earthquakes in the Kuril Trench based on Late-Holocene tsunami deposits in the distal region from wave sources, northern Hidaka, Hokkaido, Japan

Shallow slip amplification and enhanced tsunami hazard unravelled by dynamic simulations of mega-thrust earthquakes

Fully Coupled Simulations of Megathrust Earthquakes and Tsunamis in the Japan Trench, Nankai Trough, and Cascadia Subduction Zone

Resolving Shallow Moment Releases of Megathrust Earthquakes by Empirical Green’s Function‐Based Inversions: A New Approach for Tsunami Warning