Search for: All records

Abstract Seismic risk estimates are greatly improved with an increased understanding of historical (and pre‐historical) seismic events. Although Bayesian inference has been shown to provide reasonable estimates of the location and magnitude of historical earthquakes from anecdotal tsunamigenic evidence, the validity and robustness of such an approach has yet to be definitively demonstrated. Thus, in this article we present a careful analysis of the uncertainty inherent to this statistical recreation of historical seismic events. Using a priori estimates on the posterior and numerical approximations of the Hessian, we demonstrate that the 1852 Banda Sea earthquake and tsunami is well‐understood given certain explicit hypotheses. Using the same techniques we also find that the 1820 south Sulawesi event may best be explained by a dual fault rupture, best attributed to the Kalatoa fault potentially conjoining the Flores thrust and Walanae/Selayar fault. 

Abstract This article studies two particular algorithms, a relaxation least squares algorithm and a relaxation Newton iteration scheme, for reconstructing unknown parameters in dissipative dynamical systems. Both algorithms are based on a continuous data assimilation (CDA) algorithm for state reconstruction of Azouaniet al(2014J. Nonlinear Sci.24277–304). Due to the CDA origins of these parameter recovery algorithms, these schemes provide on-the-fly reconstruction, that is, as data is collected, of unknown state and parameters simultaneously. It is shown how both algorithms give way to a robust general framework for simultaneous state and parameter estimation. In particular, we develop a general theory, applicable to a large class of dissipative dynamical systems, which identifies structural and algorithmic conditions under which the proposed algorithms achieve reconstruction of the true parameters. The algorithms are implemented on a high-dimensional two-layer Lorenz 96 model, where the theoretical conditions of the general framework are explicitly verifiable. They are also implemented on the two-dimensional Rayleigh–Bénard convection system to demonstrate the applicability of the algorithms beyond the finite-dimensional setting. In each case, systematic numerical experiments are carried out probing the efficacy of the proposed algorithms, in addition to the apparent benefits and drawbacks between them.