skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.


Title: Sharpen data-driven prediction rules of individual large earthquakes with aid of Fourier and Gauss
Abstract Predicting individual large earthquakes (EQs)’ locations, magnitudes, and timing remains unreachable. The author’s prior study shows that individual large EQs have unique signatures obtained from multi-layered data transformations. Via spatio-temporal convolutions, decades-long EQ catalog data are transformed into pseudo-physics quantities (e.g., energy, power, vorticity, and Laplacian), which turn into surface-like information via Gauss curvatures. Using these new features, a rule-learning machine learning approach unravels promising prediction rules. This paper suggests further data transformation via Fourier transformation (FT). Results show that FT-based new feature can help sharpen the prediction rules. Feasibility tests of large EQs ($$M\ge$$ M 6.5) over the past 40 years in the western U.S. show promise, shedding light on data-driven prediction of individual large EQs. The handshake among ML methods, Fourier, and Gauss may help answer the long-standing enigma of seismogenesis.  more » « less
Award ID(s):
1931380 2129796
PAR ID:
10491691
Author(s) / Creator(s):
Publisher / Repository:
Nature
Date Published:
Journal Name:
Scientific Reports
Volume:
13
Issue:
1
ISSN:
2045-2322
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract Statistical descriptions of earthquakes offer important probabilistic information, and newly emerging technologies of high-precision observations and machine learning collectively advance our knowledge regarding complex earthquake behaviors. Still, there remains a formidable knowledge gap for predicting individual large earthquakes’ locations and magnitudes. Here, this study shows that the individual large earthquakes may have unique signatures that can be represented by new high-dimensional features—Gauss curvature-based coordinates. Particularly, the observed earthquake catalog data are transformed into a number of pseudo physics quantities (i.e., energy, power, vorticity, and Laplacian) which turn into smooth surface-like information via spatio-temporal convolution, giving rise to the new high-dimensional coordinates. Validations with 40-year earthquakes in the West U.S. region show that the new coordinates appear to hold uniqueness for individual large earthquakes ($$M_w \ge 7.0$$ M w 7.0 ), and the pseudo physics quantities help identify a customized data-driven prediction model. A Bayesian evolutionary algorithm in conjunction with flexible bases can identify a data-driven model, demonstrating its promising reproduction of individual large earthquake’s location and magnitude. Results imply that an individual large earthquake can be distinguished and remembered while its best-so-far model can be customized by machine learning. This study paves a new way to data-driven automated evolution of individual earthquake prediction. 
    more » « less
  2. Abstract We introduce a family of Finsler metrics, called the$$L^p$$ L p -Fisher–Rao metrics$$F_p$$ F p , for$$p\in (1,\infty )$$ p ( 1 , ) , which generalizes the classical Fisher–Rao metric$$F_2$$ F 2 , both on the space of densities$${\text {Dens}}_+(M)$$ Dens + ( M ) and probability densities$${\text {Prob}}(M)$$ Prob ( M ) . We then study their relations to the Amari–C̆encov$$\alpha $$ α -connections$$\nabla ^{(\alpha )}$$ ( α ) from information geometry: on$${\text {Dens}}_+(M)$$ Dens + ( M ) , the geodesic equations of$$F_p$$ F p and$$\nabla ^{(\alpha )}$$ ( α ) coincide, for$$p = 2/(1-\alpha )$$ p = 2 / ( 1 - α ) . Both are pullbacks of canonical constructions on$$L^p(M)$$ L p ( M ) , in which geodesics are simply straight lines. In particular, this gives a new variational interpretation of$$\alpha $$ α -geodesics as being energy minimizing curves. On$${\text {Prob}}(M)$$ Prob ( M ) , the$$F_p$$ F p and$$\nabla ^{(\alpha )}$$ ( α ) geodesics can still be thought as pullbacks of natural operations on the unit sphere in$$L^p(M)$$ L p ( M ) , but in this case they no longer coincide unless$$p=2$$ p = 2 . Using this transformation, we solve the geodesic equation of the$$\alpha $$ α -connection by showing that the geodesic are pullbacks of projections of straight lines onto the unit sphere, and they always cease to exists after finite time when they leave the positive part of the sphere. This unveils the geometric structure of solutions to the generalized Proudman–Johnson equations, and generalizes them to higher dimensions. In addition, we calculate the associate tensors of$$F_p$$ F p , and study their relation to$$\nabla ^{(\alpha )}$$ ( α )
    more » « less
  3. Abstract Let$$\phi $$ ϕ be a positive map from the$$n\times n$$ n × n matrices$$\mathcal {M}_n$$ M n to the$$m\times m$$ m × m matrices$$\mathcal {M}_m$$ M m . It is known that$$\phi $$ ϕ is 2-positive if and only if for all$$K\in \mathcal {M}_n$$ K M n and all strictly positive$$X\in \mathcal {M}_n$$ X M n ,$$\phi (K^*X^{-1}K) \geqslant \phi (K)^*\phi (X)^{-1}\phi (K)$$ ϕ ( K X - 1 K ) ϕ ( K ) ϕ ( X ) - 1 ϕ ( K ) . This inequality is not generally true if$$\phi $$ ϕ is merely a Schwarz map. We show that the corresponding tracial inequality$${{\,\textrm{Tr}\,}}[\phi (K^*X^{-1}K)] \geqslant {{\,\textrm{Tr}\,}}[\phi (K)^*\phi (X)^{-1}\phi (K)]$$ Tr [ ϕ ( K X - 1 K ) ] Tr [ ϕ ( K ) ϕ ( X ) - 1 ϕ ( K ) ] holds for a wider class of positive maps that is specified here. We also comment on the connections of this inequality with various monotonicity statements that have found wide use in mathematical physics, and apply it, and a close relative, to obtain some new, definitive results. 
    more » « less
  4. Abstract Let us fix a primepand a homogeneous system ofmlinear equations$$a_{j,1}x_1+\dots +a_{j,k}x_k=0$$ a j , 1 x 1 + + a j , k x k = 0 for$$j=1,\dots ,m$$ j = 1 , , m with coefficients$$a_{j,i}\in \mathbb {F}_p$$ a j , i F p . Suppose that$$k\ge 3m$$ k 3 m , that$$a_{j,1}+\dots +a_{j,k}=0$$ a j , 1 + + a j , k = 0 for$$j=1,\dots ,m$$ j = 1 , , m and that every$$m\times m$$ m × m minor of the$$m\times k$$ m × k matrix$$(a_{j,i})_{j,i}$$ ( a j , i ) j , i is non-singular. Then we prove that for any (large)n, any subset$$A\subseteq \mathbb {F}_p^n$$ A F p n of size$$|A|> C\cdot \Gamma ^n$$ | A | > C · Γ n contains a solution$$(x_1,\dots ,x_k)\in A^k$$ ( x 1 , , x k ) A k to the given system of equations such that the vectors$$x_1,\dots ,x_k\in A$$ x 1 , , x k A are all distinct. Here,Cand$$\Gamma $$ Γ are constants only depending onp,mandksuch that$$\Gamma Γ < p . The crucial point here is the condition for the vectors$$x_1,\dots ,x_k$$ x 1 , , x k in the solution$$(x_1,\dots ,x_k)\in A^k$$ ( x 1 , , x k ) A k to be distinct. If we relax this condition and only demand that$$x_1,\dots ,x_k$$ x 1 , , x k are not all equal, then the statement would follow easily from Tao’s slice rank polynomial method. However, handling the distinctness condition is much harder, and requires a new approach. While all previous combinatorial applications of the slice rank polynomial method have relied on the slice rank of diagonal tensors, we use a slice rank argument for a non-diagonal tensor in combination with combinatorial and probabilistic arguments. 
    more » « less
  5. We derive BM-like continuum models for the bands of superlattice heterostructures formed out of Fe-chalcogenide monolayers: (I) a single monolayer experiencing an external periodic potential, and (II) twisted bilayers with long-range moire tunneling. A symmetry derivation for the inter-layer moire tunnelling is provided for both the\Gamma Γ andM M high-symmetry points. In this paper, we focus on moire bands formed from hole-band maxima centered on\Gamma Γ , and show the possibility of moire bands withC=0 C = 0 or±1 ± 1 topological quantum numbers without breaking time-reversal symmetry. In theC=0 C = 0 region for\theta→0 θ 0 (and similarly in the limit of large superlattice period for I), the system becomes a square lattice of 2D harmonic oscillators. We fit our model to FeSe and argue that it is a viable platform for the simulation of the square Hubbard model with tunable interaction strength. 
    more » « less