# FinDerS(+): Real-Time Earthquake Slip Profiles and Magnitudes Estimated from Backprojected Displacement with Consideration of Fault Source Maturity Gradient

Maren Böse, Allie Hutchison, Isabelle Manighetti, Jiawei Li, Frédérick Massin, John Francis Clinton

► To cite this version:

Maren Böse, Allie Hutchison, Isabelle Manighetti, Jiawei Li, Frédérick Massin, et al. FinDerS(+): Real-Time Earthquake Slip Profiles and Magnitudes Estimated from Backprojected Displacement with Consideration of Fault Source Maturity Gradient. Frontiers in Earth Science, 2021, 9, pp.685879. 10.3389/feart.2021.685879. hal-03578549

HAL Id: hal-03578549

https://hal.science/hal-03578549

Submitted on 17 Feb 2022

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
<figure><img src="image_2.png" /></figure>

Check for
updates

# FinDerS(+): Real-Time Earthquake Slip Profiles and Magnitudes Estimated from Backprojected Displacement with Consideration of Fault Source Maturity Gradient

Maren Böse¹*, Allie A. Hutchison¹,², Isabelle Manighetti², Jiawei Li³, Frédérick Massin¹ and John Francis Clinton¹

¹Eidgenössische Technische Hochschule (ETH), Swiss Seismological Service (SED), Zurich, Switzerland, ²Université Côte d'Azur, Observatoire de la Côte d'Azur, IRD, CNRS, Géoazur, Sophia Antipolis, Nice, France, ³Southern University of Science and Technology (SUSTech), Institute of Risk Analysis, Prediction and Management (Risks-X), Shenzhen, China

The Finite-Fault Rupture Detector (FinDer) algorithm computes rapid line-source rupture models from high-frequency seismic acceleration amplitudes (PGA). In this paper, we propose two extensions to FinDer, called FinDerS and FinDerS+, which have the advantage of taking into account a geological property of the source fault, its structural maturity, as well as its relation to the earthquake slip distribution. These two new algorithms calculate real-time earthquake slip profiles by backprojecting seismic and/or geodetic displacement amplitudes onto the FinDer line-source. This backprojection is based on a general empirical equation established in previous work that relates dynamic peak ground displacement (PGD) at the stations to on-fault coseismic slip. While FinDerS projects PGD onto the current FinDer line-source, FinDerS+ allows the rupture to grow beyond the current model extent to predict future rupture evolution. For an informed interpolation and smoothing of the estimated slip values, FinDerS and FinDerS+ both employ a generic empirical function that has been shown to relate the along-strike gradient of structural maturity of the ruptured fault, the earthquake slip distribution, and the rupture length. Therefore, while FinDer derives magnitudes from a relatively uncertain and general empirical rupture length-magnitude relations, FinDerS and FinDerS+ provide alternate and better informed magnitude estimates using the mean slip of the profiles derived from the integration of fault source maturity. The two new algorithms can incorporate both seismic strong-motion and geodetic displacement data. In order to recover PGD from strong-motion instruments, we double-integrate and high-pass filter (> 0.075 Hz) the seismic acceleration records. Together, the three algorithms exploit the full spectrum of ground-motions, including high frequencies to derive a source fault model (FinDer) and low frequencies to determine the static offsets along this model (FinDerS and FinDerS+). We test the three algorithms for the 2019 M<sub>W</sub> 7.1 Ridgecrest (California), 2016 M<sub>W</sub> 7.0 Kumamoto (Japan), and 2008 M<sub>W</sub> 7.9 Wenchuan (China) earthquakes. Conclusively, low-frequency PGD data and integration of the fault maturity gradient do not speed-up calculations for these events, but provide additional information on slip distribution and

OPEN ACCESS

**Edited by:**
Joanna Faure Walker,
University College London,
United Kingdom

**Reviewed by:**
Yuki Kodera,
Japan Meteorological Agency, Japan
Francesco Iezzi,
University of Studies G. d'Annunzio
Chieti and Pescara, Italy

*Correspondence:
Maren Böse
mboese@sed.ethz.ch

**Specialty section:**
This article was submitted to Geohazards and Georisks, a section of the journal Frontiers in Earth Science

**Received:** 26 March 2021
**Accepted:** 02 July 2021
**Published:** 04 August 2021

**Citation:**
Böse M, Hutchison AA, Manighetti I, Li J, Massin Fé and Clinton JF (2021) FinDerS(+): Real-Time Earthquake Slip Profiles and Magnitudes Estimated from Backprojected Displacement with Consideration of Fault Source Maturity Gradient. *Front. Earth Sci.* 9:685879
doi: 10.3389/feart.2021.685879
final rupture length, as well as alternative estimates of magnitudes that can be useful to check for consistency across the algorithm suite. The FinDer algorithms systematically outperform previously established real-time PGD-based magnitude estimates in terms of speed and accuracy. The resulting slip distributions can be useful for improved ground-motion prediction given the observed relationship between seismic radiation and fault maturity.

**Keywords:** earthquake early warning, seismology, earthquake, natural Hazard, earthquake magnitude, fault properties, fault maturity, rupture determinism

# 1 INTRODUCTION

Earthquake early warning (EEW) uses real-time data from an ongoing earthquake to provide seconds of warning to people and users prior to the arrival of strong ground motions (e.g. Allen et al., 2009; Clinton et al., 2016; Allen and Melgar, 2019). Using the earliest radiated energy, EEW systems attempt to rapidly characterize the final size of an earthquake and to predict seismic ground motions in potentially affected areas. The timing of EEW depends on multiple factors, including the speed of the earthquake fault source location and size characterization, the distance from the earthquake's source to the alert recipient, as well as delays for data transmission and processing. Different approaches exist to constrain the location and size of the ongoing earthquake, including for instance the Earthquake Point-Source Integrated Code (EPIC; Chung et al. (2019), Virtual Seismologist [VS; Cua et al. (2009)], PRobabilistic and Evolutionary early warning SysTem [PRESTo; Satriano et al. (2011)], Propagation of Local Undamped Motion [PLUM; Kodera et al. (2018)], and Finite-Fault Rupture Detector [FinDer; Böse et al. (2018)] algorithms. In addition to these regional, network-based EEW approaches, a number of faster, though less accurate, single-station onsite algorithms have been proposed (e.g. Wu et al., 2006; Böse et al., 2012a). Regional EEW algorithms differ in how they utilize the energy radiated from a growing earthquake to derive information about the source fault. EPIC, for instance, uses trilateration and a grid search to determine the earthquake location, while the magnitude is estimated from empirical scaling relations (Chung et al., 2019). FinDer generates line-source models (i.e., location, strike, and length of the fault source) from the spatial distribution of high-frequency ground motions (Böse et al., 2012b). PLUM, by contrast, does not determine the earthquake fault source properties, but simply extrapolates observed motions to larger distances (Kodera et al., 2018).

Large earthquakes provide the best opportunity to implement EEW: they are associated with long fault ruptures of tens to hundreds of kilometers in length. Since earthquake ruptures typically propagate at fairly low speed (about 2.8 km/s), warning times to affected areas can thus exceed several tens of seconds, allowing efficient EEW. In small and moderate-sized (M < 6.0) earthquakes, by contrast, the rupture length is short, such that the strongest shaking typically occurs in small areas around the epicenter only; to be effective, warnings would need to be issued within a few seconds or less, which is challenging and in many cases impossible. In large earthquakes, however, EEW

requires rapid determination of finite-source fault dimensions (in particular of rupture length) in order to predict ground motions and warning areas as those mainly depend on the distance to the fault rupture (Böse et al., 2012b). This is even more challenging, as it is unclear how large the rupture will eventually grow. For example, Meier et al. (2017) showed that the source time functions - the rates at which energy is released from the earthquake fault source - of subduction-zone earthquakes do not deviate until they are halfway over; this implies that the earthquake's final size is not implicit until the event is 50% over. Alternatively, other studies suggest that information is contained within the first seconds to tens of seconds following the earthquake origin time, suggesting some determinism in the rupture behavior (e.g., Colombelli et al. (2014); Goldberg et al. (2018); Denolle (2019); Melgar and Hayes (2019); Hutchison et al. (2020)).

A recent study by Hutchison et al. (2020) showed that the final rupture length of an earthquake can be predicted from 20%, and its magnitude from 15% of the way through the rupture length, if the earthquake slip is known accurately and some intrinsic long-term properties of the source fault, namely its structural maturity, were considered. Structural maturity relates to the longevity of fault slip over geological time; the longer the slip history, the more mature the fault is overall (Manighetti et al., 2007). Additionally, a fault extends laterally (i.e., propagates) as it grows over the long-term (commonly, millions of years), generating a gradient in structural maturity along its length: the more mature part of the fault is where it originally initiated, while the fault becomes increasingly more immature towards its propagating tip(s) (Manighetti et al., 2001). Interestingly, as a fault or fault section becomes more mature, some of its geometrical (i.e., segment connections) and mechanical (damage compliance and possibly friction) properties evolve, and these changes impact an earthquake's behavior (Wesnousky, 1988; Manighetti et al., 2007; Perrin et al., 2016a). In particular, coseismic earthquake slip is greatest on the most mature section of the ruptured fault, generating an asymmetry in the earthquake slip-length distribution (Manighetti et al., 2005; Perrin et al., 2016a). Hutchison et al. (2020) formalized this generic relation between the along-strike fault maturity gradient and the earthquake slip asymmetry, and demonstrated that this empirical equation can be used to anticipate the final earthquake rupture length from several slip values measured in the first stages of the rupture growth, provided these slip values are accurately determined.
<figure><img src="image_3.png" /></figure>

**FIGURE 1** | Schematic illustration of FinDer algorithms, including the original FinDer main algorithm that computes line-source models from seismic acceleration data, and the new FinDerS and FinDerS+ algorithms (from **top** to **bottom**). FinDerS determines slip profiles by backprojecting peak ground displacement (PGD) amplitudes onto the FinDer line-source. FinDerS+ does the same, but allows the line-source to grow beyond the current FinDer model. This enables FinDerS+ to predict future rupture growth, while FinDer and FinDerS both characterize the (close to) current rupture state. PGD amplitudes are derived from integrated and high-pass (> 0.075 Hz) filtered seismic strong-motion and/or geodetic displacement records. We use a generic empirical equation derived in Hutchison et al. (2020) for an informed interpolation and smoothing of slip values. The two new algorithms provide alternate magnitude estimates that are based on seismic slip, while FinDer magnitudes are computed from the length of the line-source estimate. Together, the three algorithms exploit the full spectrum of ground-motions, including high frequencies to derive a source model (FinDer) and low frequencies to determine the static offsets along this source (FinDerS and FinDerS+).

However, the study of Hutchison et al. (2020) was done on final static earthquake slip profiles, that is, the time required to reach and to calculate final displacements is neglected. Furthermore, the slip profiles are accurate as they were measured in the field or remote data right after the earthquakes. Here, we explore the application of the algorithm in a retrospective study by simulating the real-time streams of seismic and geodetic data in a regional network recording a growing earthquake. In order for this to operate, it is necessary that we can generate an earthquake slip profile in real-time. This requires two pieces of information: 1) a finite-source model that characterizes the spatial dimensions (i.e., location, length) and orientation (i.e., strike) of the fault source and that updates in real-time based on the growing available data, and 2) a real-time estimated slip distribution along this source model. For the first element, we utilize in this study the line-source models computed from the Finite-Fault Rupture Detector (FinDer) algorithm (Böse et al., 2012b, 2015, 2018). For the second element, we utilize an empirically-based algorithm (Aagaard et al., 2004; Yamada, 2007) that backprojects dynamic peak ground displacement (PGD) amplitudes from individual stations onto this source model to determine slip (FinDerS, FinDerS+).

We test the FinDer algorithms here on three earthquakes for which available data are dense enough: the 2019 M<sub>W</sub> 7.1 Ridgecrest, California, the 2016 M<sub>W</sub> 7.0

Kumamoto, Japan, and the 2008 M<sub>W</sub> 7.9 Wenchuan, China, earthquakes.

## 2 METHODS

The three FinDer algorithms presented in this paper exploit the full spectrum of seismic ground-motions (Figure 1): FinDer (Böse et al., 2012b; Böse et al., 2015; Böse et al., 2018) derives a line-source model from high-frequency amplitudes; FinDerS and FinDerS+ determine the static offsets along this model from low-frequency displacement. FinDerS and FinDerS+ both estimate slip profiles along the source fault by backprojecting dynamic displacement (PGD) amplitudes onto the FinDer source model. However, while FinDerS projects PGD onto the current line-source, FinDerS+ allows the rupture to grow beyond the current model to predict future rupture evolution. For an informed interpolation and smoothing of the estimated slip values, FinDerS and FinDerS+ both employ the generic empirical relationship developed by Hutchison et al. (2020) that relates the along-strike gradient of long-term structural maturity of the ruptured fault to the coseismic slip distribution along the rupture. With both PGA and PGD changing over time,
estimates from FinDer, FinDerS, and FinDerS+ are continuously updated. Each of these processes is described in more detail in the following.

## 2.1 FinDer: Current Rupture Line-Source Model and Magnitude, MFD

The Finite-Fault Rupture Detector (FinDer) algorithm (Böse et al. (2012b), Böse et al. (2015), Böse et al. (2018)) provides rapid earthquake line-source models, that describe the position, length (LFD) and strike of a fault rupture. These parameters are computed from the spatial distribution of high-frequency PGA amplitudes recorded across a dense seismic network. Using template matching, FinDer compares these motions with theoretical spatial PGA templates, which are calculated from an empirical ground-motion model (here: Cua and Heaton (2009)) derived from worldwide earthquake data with different magnitudes, rupture lengths, and source mechanisms. The best template and resulting FinDer line-source model is found from a combined grid-search and divide-and-conquer approach (Böse et al., 2018). The FinDer magnitude, MFD, is estimated from a selected empirical rupture length-magnitude relation [typically Wells and Coppersmith (1994)]. Comparing different relationships (Wells and Coppersmith, 1994; Blaser et al., 2010; Leonard, 2014) suggest that the resulting magnitude estimates have an uncertainty of about +/−0.3 magnitude units.

## 2.2 FinDerS: Slip Profile and Magnitude, MFD

FinDerS (“S” stands for slip) determines 1D slip profiles from the backprojection of PGD amplitudes onto the FinDer line-source. The backprojection is done by employing a general empirical equation developed by Yamada (2007) and Aagaard et al. (2004) from worldwide earthquake data and simulations. This equation relates the seismic slip (or static offset), D, to the dynamic PGD amplitudes measured at closest distance, R, of the surface projected rupture

$$D = \frac{PGD\sqrt{1 + 0.125R^{1.55}}}{0.7} \quad (1)$$

where D and PGD are in cm and R in km. According to Eq. 1 the dynamic PGD is approximately 2/3 of the static slip along the ruptured fault. As originally proposed by Böse et al. (2013), we apply Eq. 1 here to dynamically backproject PGD observations (i.e. backproject the time-dependent data) onto the FinDer line-source, which will also evolve and grow over time. Here, the FinDerS slip profile has, at each time step, the length of the FinDer line-source model. We only backproject PGD from stations within a certain distance to the fault (here: R < 30 km, R < 60 km, and R < 90 km, respectively).

Then, to interpolate and smooth the estimated backprojected slip values, FinDerS employs the empirical relationship developed in Hutchison et al. (2020) that relates the along-strike gradient of long-term structural maturity of the ruptured fault to the coseismic slip distribution along the rupture. This equation solves for D at any given point along the rupture, l, up to length L using a least squares fitting technique:

$$D(l) = \frac{3D_{peak}l^q}{L^{2q}}(L^{2q} - l^{2q})^{3/4} \quad (2)$$

where Dpeak and q represent amplitude and skewedness of the slip profile, respectively. These are interdependent variables that change with each other. See Hutchison et al. (2020) for details.

Finally, as described in Hutchison et al. (2020), we estimate the (moment) magnitude from the profile average slip value, Dmean, and the rupture length, L, at each time step:

$$M \approx \frac{2}{3}\log_{10}(2LD_{mean} * 10^7) + 1 \quad (3)$$

Since FinDerS is based upon the FinDer line-source, L corresponds in Eqs. 2, 3 to LFD. The resulting magnitude is named M = MFDS. Since both the FinDer line-source and PGD evolve over time, Eqs. 1–3) are dynamically recomputed.

## 2.3 FinDerS+: Predicted Final Rupture Length, Slip Profile and Magnitude, MFD

As opposed to FinDerS, which fits the slip profile to the current FinDer line-source, FinDerS+ does not restrict the final rupture length and allows the FinDer-determined fault rupture to continue growing towards both directions up to a maximum reasonable rupture length, which we here set to 500 km (Manighetti et al., 2007). The slip values are determined as in FinDerS, by backprojection of PGD amplitudes. Then, as for FinDerS, FinDerS+ interpolates the slip values by using Eq. 2, where L = LFD (same in Eq. 3). Again, FinDerS+ estimates are updated over time. From the slip profiles derived at each step, FinDerS+ predicts the final rupture length, the final slip distribution, and the final magnitude, M = MFDS.

## 2.4 PGD-Based Magnitude, MPGD

Throughout this paper, we will compare the magnitude results of the three FinDer algorithms to estimates obtained from more simple, previously established PGD-magnitude scaling relationships of the form

$$M_{PGD} = (\log(PGD) - A)/(B + \text{Clog}(R)) \quad (4)$$

where R is the distance in km between the earthquake hypocenter and the station at which the PGD is observed. We are testing here three sets of coefficients published by various authors: 1) A = −5.919, B = 1.009, C = −0.145 (Ruhl et al., 2019); 2) A = −4.434, B = 1.047, C = −0.138 (Melgar et al., 2015), and 3) A = −6.687, B = 1.5, C = −0.214 (Crowell et al., 2016), respectively. PGD, here defined as PGD = max√E(t)² + N(t)² + Z(t)² (with E, N, and Z the East, North and vertical displacement components), is in cm for Melgar et al. (2015) and Crowell et al. (2016) and in m for Ruhl et al. (2019). The regression for the three coefficients A, B, and C was done with different global GNSS datasets of continental and subduction-zone earthquakes (Ruhl et al., 2019). We name the three resulting PGD-based magnitudes M Melgar , M Crowell , and M Ruhl , respectively.
<table><caption>**TABLE 1** | Source parameters and final FinDer line-source parameters for the M<sub>w</sub> 7.1 Ridgecrest (California), M<sub>w</sub> 7.0 Kumamoto (Japan), and M<sub>w</sub> 7.9 Wenchuan (China) earthquakes.</caption><thead><tr><th>Name</th><th>Origin time [UTC]</th><th>Latitude [degrees]</th><th>Longitude [degrees]</th><th>Depth [km]</th><th>MW</th><th>Approx. rupture length and duration [km] [s]</th><th>Source mechanism</th><th>Final FinDer linesource parameters<br>• Time from origin [s]<br>• Length [km]<br>• Strike [degrees]<br>• lat1/lon1<br>• lat2/lon2</th></tr></thead><tbody><tr><td>Ridgecrest</td><td>2019-07-06 03:19:53</td><td>35.770</td><td>-117.599</td><td>8</td><td>7.1</td><td>50 20</td><td>Strike-slip</td><td>• 26<br>• 53<br>• 145<br>• 35.92/-117.73<br>• 35.53/-117.40</td></tr><tr><td>Kumamoto*</td><td>2016-04-15 16:25:06</td><td>32.791</td><td>130.754</td><td>10</td><td>7.0</td><td>40 (65) 15</td><td>Strike-slip and normal</td><td>• 36<br>• 85<br>• 30<br>• 32.61/130.81<br>• 33.27/131.26</td></tr><tr><td>Yufuin*</td><td>2016-04-15 16:25:39</td><td>33.266</td><td>131.340</td><td>5</td><td>5.7-6.5</td><td></td><td>Normal and strike-slip</td><td>• 124<br>• 290<br>• 50</td></tr><tr><td>Wenchuan</td><td>2008-05-12 06:28:01</td><td>31.002</td><td>103.322</td><td>19</td><td>7.9</td><td>300 100</td><td>Thrust and strike-slip</td><td>• 31.24/103.33<br>• 32.92/105.69</td></tr></tbody></table>

## 2.5 Magnitude from Seismic Moment-Rate Function, M<sub>mr</sub>

Finally, we will compare all magnitude estimates to the moment-rate function-derived magnitude, M<sub>mr</sub>. To determine M<sub>mr</sub> at time t relative to the rupture nucleation time, we use the moment magnitude definition of Hanks and Kanamori (1979), M<sub>mr</sub>(t) = (log<sub>10</sub>(m<sub>0</sub>(t)) - 9.1), where the seismic moment at time t is m<sub>0</sub>(t) = ∫<sub>0</sub><sup>t</sup> Ṁ<sub>0</sub>(t)dt and Ṁ<sub>0</sub> is the moment-rate function, which we here take from the US Geological Survey (USGS) finite-fault database.

## 3 DATA AND PREPROCESSING

### 3.1 Earthquakes

We will demonstrate our suite of FinDer algorithms for three continental earthquakes: the 2019 M<sub>w</sub> 7.1 Ridgecrest earthquake with a right-lateral slip, the 2016 M<sub>w</sub> 7.0 Kumamoto earthquake with a dominant right-lateral and additional normal slip, and finally the 2008 M<sub>w</sub> 7.9 Wenchuan earthquake that had a right lateral and reverse slip (Table 1). We select these earthquakes mainly because of their availability of seismic (and in the cases of Ridgecrest and Kumamoto of geodetic) data, and their large rupture sizes (M<sub>w</sub> 7.0+). As will be shown later, the slip distributions of the three earthquakes show the generic asymmetry encapsulated in the empirical equation from Hutchison et al. (2020) (Eq. 2), even though the surface slip data for the Kumamoto earthquake are few. This means that the gradient of maturity along the three ruptured faults can be

derived from the empirical Eq. 2. Furthermore, the performance of the FinDer algorithm, although for older code versions, has been demonstrated and documented for these three events in previous studies (Böse et al., 2018; Chung et al., 2020; Li et al., 2020). Finally, the empirical equation from Hutchison et al. (2020) has been previously applied to the static slip distribution of the Wenchuan earthquake. It demonstrated a very good performance on the magnitude prediction (from the rupture start), and modest results on the length prediction (from ~50% throughout the rupture), due to large displacements over the first 100 km of the rupture.

The July 5, 2019 M<sub>w</sub> 7.1 Ridgecrest earthquake ruptured a major NW-trending right-lateral fault in the eastern California shear zone, along with many secondary subparallel faults, and others of sub-perpendicular orientation (Milliner and Donnellan, 2020). The earthquake produced strong shaking with up to Modified Mercalli Intensity (MMI) IX in the epicentral area, and was felt in large parts of southern California and Nevada (e.g. Barnhart et al., 2019). The earthquake mainly produced a NW-SE surface rupture of ~50 km (Hudnut et al., 2020; Milliner and Donnellan, 2020) and a maximum coseismic slip of ~4.5 m at surface, located in the northernmost part of the fault, close to the epicenter (Wang and Bürgmann, 2020). Most slip occurred at shallow depths of less than 10 km (Barnhart et al., 2019). The slip-length profile revealed to be asymmetric, with greatest slip to the north (Milliner and Donnellan, 2020). While the overall maturity of the fault ruptured in the Ridgecrest earthquake has not been described in prior works, neighboring, similar faults in the

<footnote>*Secondary (triggered) event.</footnote>
eastern California shear zone have been shown to be immature (Perrin et al., 2016a). Furthermore, according to (Perrin et al., 2016b), the structure of the Ridgecrest fault zone, with a dense network of diverging splays at its southern tip, suggests that the fault is more mature to the north, in keeping with the asymmetry of the Ridgecrest earthquake slip profile.

The M<sub>w</sub> 7.0 Kumamoto earthquake of April 16, 2016, ruptured at least ~40 km (at surface) of the NE-trending Futugawa fault, including its southern Hinagu splay (Shirahama et al., 2016; Scott et al., 2018; Milliner et al., 2020). The rupture length at depth is likely longer, up to ~65 km (Uchide et al., 2016). The earthquake accommodated mainly right-lateral slip, yet with an additional small normal component (Shirahama et al., 2016). A maximum dextral slip of about 3.2 m was measured at the ground surface (Shirahama et al., 2016; Milliner et al., 2020), while the rupture also produced ~1.5 m of vertical slip (Milliner et al., 2020). The total largest slip at surface might thus be about 3.5 m. The largest slip occurred in the northeastern part of the ruptured fault (Shirahama et al., 2016; Milliner et al., 2020), reaching about 6-7 m at 4-5 km depth (Uchide et al., 2016; Milliner et al., 2020). The strong directivity of the Kumamoto event is thought to have dynamically triggered a second event close to Yufuin, approximately 80 km to the NE relative to the mainshock epicenter (Uchide et al., 2016; Yoshida, 2016). This event might have ruptured about 20 km of fault length, north of its epicenter, in the same overall fault zone as the mainshock (Uchide et al., 2016). Because of the difficulty to isolate its waveforms from those of the mainshock, the magnitude of this secondary event has various estimates: M5.7 (Japanese Meteorological Agency, JMA), M6.0 (Yoshida, 2016), and M6.5 (Uchide et al., 2016). Like the Kumamoto mainshock, the Yufuin event has a dominant right-lateral slip, and is thought to have occurred ~30 s after the onset of the Kumamoto event (Uchide et al., 2016; Yoshida, 2016). We include the location of this event in our analysis (Table 1), since the wave trains for the two events are nearly indistinguishable and affect our results. The overall structural maturity of the Futugawa fault is unknown, but its architecture with fan-splays at its western tip suggests that it might be more mature to the northeast (Perrin et al., 2016b), in keeping with the earthquake slip asymmetry.

The M<sub>w</sub> 7.9 Wenchuan earthquake on May 12, 2008 produced seismic intensities of up to XI on the MMI scale (Chen and Booth, 2011). The event primarily ruptured the NE-trending Beichuan fault with a dextral and reverse slip, along with the adjacent Pengguan reverse fault, both within the Longmenshen fault zone. The total rupture length was about 280 km, while the earthquake produced 12-13 m of slip at the ground surface. The largest slip occurred in the southwest, close to the epicenter (Shen et al., 2009; Perrin et al., 2016a). The earthquake slip profile is markedly asymmetric with slip tapering from the SW to the NE (Shen et al. (2009); Perrin et al. (2016a)). The fault has been described as of intermediate maturity, with a gradient of maturity decreasing towards NE (Perrin et al., 2016a).

## 3.2 Data Preprocessing

For the Ridgecrest earthquake, we use strong-motion data from the Southern California Seismic Network (SCSN), as well as preprocessed 1-Hz displacement time series from 10 GNSS stations from the Geodetic Facility for the Advancement of Geoscience (GAGE) Network of the Americas (NOTA) that were obtained through UNAVCO (Mattioli et al., 2020). For the Kumamoto earthquake, we use strong-motion records from both KiK-net and K-NET stations, as well as preprocessed GEONET Global Positioning System time series downloaded from Ruhl et al. (2019). For the Wenchuan earthquake, we use strong-motion records obtained through the China Strong Motion Networks Center (CSMNC) at the Institute of Engineering Mechanics, China Earthquake Administration. These records start 20 s before the P-wave arrival. As described in Li et al. (2020) we reconstruct absolute times from the event location and origin time (Table 1) and align in each record the theoretical and observed P-wave arrivals assuming vp = 6.1 km/s. Figure 2 shows the distributions of seismic and geodetic sensors used in this study along with the final FinDer line-source models (Table 1) for the Ridgecrest, Kumamoto, and Wenchuan earthquakes at 26 s, 36 s, and 124 s from origin times, respectively.

While the FinDer line-source models are recovered directly from the high-frequency strong-motion data (Böse et al., 2018), the slip values are computed from the backprojected dynamic displacement amplitudes using Eq. 1. Following Yamada (2007; page 45), we determine the maximum peak ground displacement at each seismic/geodetic station as $PGD = max(PGD_z, \sqrt{PGD_E^2 + PGD_N^2})$, where PGD<sub>Z</sub>, PGD<sub>N</sub>, and PGD<sub>E</sub> are the peak ground displacement amplitudes along the vertical, north and east components, respectively.

Recovering displacement from strong-motion recordings is generally challenging and will be discussed further in later sections. We are doing the following: first, we cut all strong-motion waveforms so that they begin at the origin time, $t_0$. Then we calculate and remove the background median noise of each waveform by determining the median amplitude between the start time of the waveform and the arrival of the P-wave and subtracting this median value from the entire waveform. At this stage, we also detrend the data, removing the change in the mean as it changes over time. We do not apply a taper, which would artificially decrease the amplitudes at the edges of the time window. Next, we double-integrate the waveforms to obtain displacements and apply a 4th order high-pass causal butterworth filter with a cutoff-frequency of 0.075 Hz (Yamada, 2007). This filter is made to reduce long period noise artifacts that are typically induced when integrating accelerograms due to the lack of sensitivity of an accelerometer to longer periods. High-pass filtered PGDs are typically smaller than true PGDs (Yamada, 2007). We only use PGD data from seismic stations at perpendicular fault rupture distances of up to a certain distance (here 30 km, 60 km, and 90 km relative to the FinDer line-source). We remove stations that backproject off the line-source. This means the list of stations used is updated with every new solution, as the line-source geometry is constantly changing over time.

The GNSS data for Ridgecrest and Kumamoto is already preprocessed with a precise-point-positioning (PPP) algorithm
A 2019 Ridgecrest, M 7.1

<figure><img src="image_4.png" /></figure>

B 2016 Kumamoto, M 7.0

<figure><img src="image_5.png" /></figure>

C 2008 Wenchuan, M 7.9

<figure><img src="image_6.png" /></figure>

**FIGURE 2** | Maps showing earthquake epicenters (stars), seismic strong-motion and geodetic GNSS stations (triangles and squares, respectively), and FinDer line-source models (black lines) for the a) M<sub>W</sub> 7.1 Ridgecrest, California, b) M<sub>W</sub> 7.0 Kumamoto (and Yufuin), Japan, and c) M<sub>W</sub> 7.9 Wenchuan, China, earthquakes. The line-source models were calculated from seismic waveform playback and correspond to the final FinDer solutions at 26, 36, and 124 s, respectively (**Table 1**). Color of station markers shows their respective distance range relative to the FinDer line-source.
A 2019 Ridgecrest, M 7.1

<figure><img src="image_7.png" /></figure>

C 2008 Wenchuan, M 7.9

<figure><img src="image_8.png" /></figure>

B 2016 Kumamoto, M 7.0

<figure><img src="image_9.png" /></figure>

**FIGURE 3** | Peak ground displacement (PGD) amplitudes recovered from high-pass filtered seismic and geodetic recordings of the Ridgecrest, Kumamoto, and Wenchuan earthquakes, sorted as a function of the closest distance to the final FinDer line-source. Values are listed in **Supplementary Table 1 (Supplementary Material)**. Seismic (triangles) and geodetic (squares) PGD values generally agree well. In the new FinDerS and FinDerS+ algorithms we backproject the time-varying PGD amplitudes onto the evolving FinDer line-source to estimate seismic slip profiles at any given time through the rupture.
<figure><img src="image_10.png" /></figure>

**FIGURE 4** | Fault slip, D, estimated from the backprojection of final PGD amplitudes onto the final FinDer line-source (Figure 2) for the a) Ridgecrest, b) Kumamoto, and c) Wenchuan earthquakes. Gray dots show surface slip profiles published by various authors as given in the legends. For Ridgecrest and Kumamoto, the published values show lateral slip, while for Wenchuan, they show net slip. Note that the backprojected slip values and epicenters are relative to the respective FinDer line-source model, oriented as indicated on the plots.

and in displacement (Ruhl et al., 2019; Mattioli et al., 2020), so we do not apply any filters. We simply detrend the data and, like the seismic data, we remove the pre-event median of the data by calculating the median amplitude between the start time of the waveform and the arrival of the P-wave, assuming vp = 6.1 km/s, and subtracting this value from the entire waveform. We use all geodetic stations within 100 km from the FinDer line-source regardless of the seismic threshold we are testing (30 km, 60 km, 90 km). We remove again stations that backproject onto the tips of the FinDer line-source.

The resulting displacement waveforms are shown in the Supplementary Material of this paper (**Supplementary Figures 2A–4**). The extracted final PGD amplitudes from the preprocessed seismic and geodetic waveforms are summarized in **Supplementary Table 1** (*Supplementary Material*) and plotted in **Figure 3**. As expected the amplitudes decay as a function of distance from the (final) FinDer line-source with values of about 1 m close to the rupture and of a few cm at about 100 km distance. PGD amplitudes extracted from seismic and GNSS recordings generally agree well.

# 4 RESULTS

## 4.1 Static Application: Final Rupture Lengths and Slips

Figure 4 shows the backprojection results of final dynamic PGD amplitudes onto the final FinDer line-source models (Table 1)

compared to measured surface slip profiles available in the literature. Even though Eq. 1 has been calibrated with the slip at depth, for a steep-dipping fault (as examined here) the results should provide a first-order approximation of surface slip. As shown in (Manighetti et al., 2005), the general (oblique-triangular) pattern of the slip profiles, which we take advantage of in Eq. 2, persists at depth.

For Ridgecrest (Figure 4A), we compare backprojected values with surface slip values derived from subpixel correlation of high-resolution optical imagery from two different satellites (Milliner and Donnellan, 2020). The backprojected seismic and geodetic slip values follow the overall distribution of surface slip along strike, though, the maximum slip value identified by Milliner and Donnellan (2020) of ~400 cm, is not matched neither by the backprojected geodetic nor seismic data, where the maximum slip value is ~270 cm. Consistent with the inference of Milliner and Donnellan (2020), the closest stations which best constrain the slip suggest that largest displacements occurred in the northern part of the rupture.

For Kumamoto (Figure 4B), the rupture trace observed at the ground surface was only ~40 km long (Shirahama et al., 2016; Scott et al., 2018; Milliner et al., 2020), while the rupture at depth was likely longer, ~65 km (Uchide et al., 2016). The final FinDer line-source is estimated as ~84 km long, that is 20 km longer than the likely rupture length. This is because the model integrates both the mainshock and the subsequent triggered Yufuin event towards the NE, and also because it slightly overshoots the rupture towards the SW relative to the
epicenter, possibly because of FinDer's simplified line-source approximation. We can thus only compare backprojected slip values with the first ~ 40 km of the FinDer line-source. Here, we compare our values with horizontal surface slip values measured at the surface by Shirahama et al. (2016). Like for the Ridgecrest earthquake, the slip values compare quite well, particularly for stations closer than 60 km to the rupture. Overall, as observed in the field, they suggest a slip increase towards the NE over the considered 40 km length.

For Wenchuan (Figure 4C) backprojected slip values also compare quite well with the surface slip profile measured in the field. Like in the other events, however, the maximum backprojected slip values (~ 8 m) never reach the surface slip values reported in the literature (~ 12 m). At the beginning of the slip profile, that is near the hypocenter, the backprojected slip values are notably lower than the measured surface slip. This may in part be due to limited station coverage, as there are only two stations within the first 50 km of the rupture. However, overall, the largest slip is well located to the SW of the rupture, as observed in the field.

In all three earthquakes, we note that the backprojected slip tends to under-estimate the actual surface slip. In general, there is no clear distance range of stations (here: 30 km, 60 km and 90 km from the line-source model) that works best for the backprojection, but stations less than 60 km from the line-source provide slip values in closer agreement with those in the literature. Therefore, in the subsequent analyses we will select the 60 km station cutoff for further demonstration of our approach. Results for the other two cutoffs are shown in the **Supplementary Material (Supplementary Figures 5, 6).**

## 4.2 Dynamic Application: Evolving Ruptures and Slips

For the dynamic application of the three FinDer algorithms (Figure 1) we run playbacks for the three earthquakes using the preprocessed seismic and geodetic waveform records (see Section 3.2) to emulate the evolution of possible output of FinDer, FinDerS, and FinDerS+ over time using the respectively available information (even if final peak values have not yet been reached). Table 2 summarizes the results (here with parameters being updated every 5 s even though a higher resolution is possible), while Figures 5-7 illustrate the results at some selected time steps. Figure 8 compares the evolution of estimated magnitudes for the various approaches.

For simplicity we neglect data latencies, which typically are on the order of 1–2 s for a fast seismic network (e.g. Behr et al., 2015).

Details on the FinDer results for the three earthquakes are given in Böse et al. (2018), Chung et al. (2020), and Li et al. (2020). Due to the FinDer trigger settings applied in this study, which requires the P-wave to have reached at least 4 neighbouring stations, FinDer triggers in this retrospective study 8 s after $t_0$ for Ridgecrest, 4 s for Kumamoto, and 13 s for Wenchuan. FinDerS and FinDerS+ require PGD data from at least three stations (and the FinDer line-source model), and results from these two algorithms are typically given a bit later (Table 2).

### 4.2.1 Playback Results for $M_W$ 7.1 Ridgecrest

For the $M_W$ 7.1 Ridgecrest earthquake (Figure 5 and Table 2) FinDer triggers 8 s after $t_0$ (Chung et al., 2020) and gives a rupture length of 29 km at 15 s, and of 45 km at 20 s. The final rupture length and magnitude at 30 s are estimated as $L_{FD}$ ~ 53 km and $M_{FD}$ 6.9, respectively. These values are in excellent agreement in terms of both length and orientation with the 54 km-long rupture determined through high-resolution optical imagery (Milliner and Donnellan, 2020). The performance of FinDerS closely follows that of FinDer, but provides an independent and more informed estimate of the magnitude as the latter is derived from the average slip resulting from backprojection and interpolation with the maturity gradient-based empirical relation (Eq. 2). We note, however, that the backprojected slip is significantly under-estimated. Furthermore, the slip profile fit of FinDerS is somewhat symmetric (Figure 5), likely due to the rupture length being fixed to the FinDer line-source determination. As a matter of fact, FinDerS+ best reproduces the asymmetry of the slip profile and produces higher performance in terms of magnitude. However, FinDerS+ over-estimates the rupture length, probably as a result of the under-estimated slip.

### 4.2.2 Playback Results for $M_W$ 7.0 Kumamoto

For the $M_W$ 7.0 Kumamoto earthquake (Figure 6 and Table 2) FinDer triggers 4 s after $t_0$ (Böse et al., 2018). We show in Figure 6 only the results until 25 s, because just after 30 s the dynamically triggered M ~ 6 Yufuin earthquake at ~ 80 km distance starts and may contaminate the results from FinDer, FinDerS, and FinDerS+. However, in Table 2 we show results of all the algorithms until they stabilize at ~ 40 s.

FinDerS+ performs very similarly to FinDerS, which indicates a well-formed slip profile resulting from a good station coverage, which does not leave much room for fitting a longer rupture in FinDerS+. Notably, despite the over-estimated rupture length, the maximum slip values are fairly consistent with observed surface slip values (Shirahama et al., 2016), particularly in the earliest time windows.

### 4.2.3 Playback Results for $M_W$ 7.9 Wenchuan

FinDer triggers 13 s after $t_0$ for the MW 7.9 Wenchuan earthquake (Li et al., 2020). FinDer performs quite well (Figure 7 and Table 2), nearly predicting the magnitude ($M_{FD}$ 7.8) correctly 50% of the way through the event duration (Figure 7C). The final magnitude prediction ($M_{FD}$ 8.0) is only slightly larger than the true magnitude.

The backprojected slip values are significantly under-estimated, by a factor of 3–6. With such poor slip estimates, we cannot expect that FinDerS and FinDerS+ perfom well. FinDerS, using the FinDer rupture length, performs less well than FinDer for magnitude, never arriving at the cataloged final value, $M_W$ 7.9. FinDerS+, however, performs somewhat better in terms of magnitude, at least in the middle time steps, with $M_{FDS+}$ 7.7 at 45 s, compared to $M_{FDS}$ 7.3 (Figure 7B). This, however, results from an over-prediction of the final rupture length, where at 45–85 s, the predicted final rupture length for FinDerS+ is 500 km (our theoretically allowed maximum value), showing that
<table><caption>**TABLE 2** | Results from FinDer (rupture length, L<sub>FD</sub>, and magnitude, M<sub>FD</sub>), FinDerS (mean slip, Dmean, and magnitude, M<sub>FDS</sub>, both from fitted slip profile) and FinDerS+ (predicted final rupture length, L<sub>FDS+</sub>, predicted final mean slip, mean, and predicted final magnitude, M<sub>FDS+</sub>) for the a) Ridgecrest, b) Kumamoto, and c) Wenchuan earthquake as a function of time from origin for a 60 km station distance cutoff. Magnitudes M<sub>mr</sub> and M<sub>PGD</sub> are estimated from the USGS moment rate function and empirical PGD-distance relationships (Melgar et al., 2015; Crowell et al., 2016; Ruhl et al., 2019), respectively. Table shows updates in increments of 5 s.</caption><caption>Time from origin [s]</caption><thead><tr><th rowspan="2"></th><th rowspan="2">Moment Rate</th><th colspan="3">FinDer (current line-source)</th><th colspan="2">FinDerS (current slip profile)</th><th colspan="3">FinDerS+ (predicted slip profile)</th><th colspan="7">PGD scaling (M<sub>PGD</sub>)</th></tr><tr><th>M<sub>mr</sub></th><th>L<sub>FD</sub> [km]</th><th>M<sub>FD</sub></th><th>D<sub>mean</sub> [m]</th><th>M<sub>FDS</sub></th><th>L<sub>FDS+</sub> [km]</th><th>D<sub>mean</sub> [m]</th><th>M<sub>FDS+</sub></th><th>M<sub>PGD</sub><sup>Crowell</sup></th><th>M<sub>PGD</sub><sup>Crowell</sup><sub>GNSS</sub></th><th>M<sub>PGD</sub><sup>Melgar</sup></th><th>M<sub>PGD</sub><sup>Melgar</sup><sub>GNSS</sub></th><th>M<sub>PGD</sub><sup>Ruhl</sup></th><th>M<sub>PGD</sub><sup>Ruhl</sup><sub>GNSS</sub></th></tr></thead><tbody><tr><td colspan="16"><strong>a) Ridgecrest earthquake</strong></td></tr><tr><td>5</td><td>6.1</td><td></td><td></td><td></td><td>—</td><td>—</td><td>—</td><td>—</td><td>—</td><td>5.4</td><td>—</td><td>5.4</td><td>—</td><td>5.0</td><td>—</td></tr><tr><td>10</td><td>6.8</td><td></td><td></td><td></td><td>—</td><td>—</td><td>—</td><td>—</td><td>—</td><td>6.0</td><td>6.2</td><td>5.9</td><td>6.1</td><td>5.7</td><td>5.8</td></tr><tr><td>15</td><td>6.9</td><td></td><td></td><td></td><td>—</td><td>—</td><td>—</td><td>—</td><td>—</td><td>6.5</td><td>6.7</td><td>6.5</td><td>6.7</td><td>6.3</td><td>6.5</td></tr><tr><td>20</td><td>7.0</td><td></td><td></td><td></td><td>—</td><td>—</td><td>—</td><td>—</td><td>—</td><td>6.6</td><td>6.8</td><td>6.6</td><td>6.8</td><td>6.4</td><td>6.6</td></tr><tr><td>25</td><td>7.0</td><td></td><td></td><td></td><td>—</td><td>—</td><td>—</td><td>—</td><td>—</td><td>6.7</td><td>6.8</td><td>6.6</td><td>6.8</td><td>6.4</td><td>6.6</td></tr><tr><td>30</td><td>7.0</td><td></td><td></td><td></td><td>—</td><td>—</td><td>—</td><td>—</td><td>—</td><td>6.8</td><td>6.9</td><td>6.7</td><td>6.9</td><td>6.4</td><td>6.7</td></tr><tr><td colspan="16"><strong>b) Kumamoto earthquake</strong></td></tr><tr><td>5</td><td>6.6</td><td></td><td></td><td></td><td>—</td><td>—</td><td>—</td><td>—</td><td>—</td><td>5.4</td><td>4.9</td><td>5.4</td><td>4.5</td><td>5.0</td><td>4.2</td></tr><tr><td>10</td><td>6.9</td><td></td><td></td><td></td><td>—</td><td>—</td><td>—</td><td>—</td><td>—</td><td>6.0</td><td>6.0</td><td>5.9</td><td>5.8</td><td>5.7</td><td>5.5</td></tr><tr><td>15</td><td>7.0</td><td></td><td></td><td></td><td>—</td><td>—</td><td>—</td><td>—</td><td>—</td><td>6.5</td><td>6.6</td><td>6.5</td><td>6.6</td><td>6.3</td><td>6.4</td></tr><tr><td>20</td><td>7.0</td><td></td><td></td><td></td><td>—</td><td>—</td><td>—</td><td>—</td><td>—</td><td>6.6</td><td>6.8</td><td>6.6</td><td>6.7</td><td>6.4</td><td>6.5</td></tr><tr><td>25</td><td>7.1</td><td></td><td></td><td></td><td>—</td><td>—</td><td>—</td><td>—</td><td>—</td><td>6.7</td><td>6.8</td><td>6.6</td><td>6.8</td><td>6.4</td><td>6.5</td></tr><tr><td>30</td><td>7.1</td><td></td><td></td><td></td><td>—</td><td>—</td><td>—</td><td>—</td><td>—</td><td>6.8</td><td>6.9</td><td>6.7</td><td>6.8</td><td>6.4</td><td>6.6</td></tr><tr><td>35</td><td>7.1</td><td></td><td></td><td></td><td>—</td><td>—</td><td>—</td><td>—</td><td>—</td><td>6.9</td><td>6.9</td><td>6.7</td><td>6.8</td><td>6.5</td><td>6.6</td></tr><tr><td>40</td><td>7.1</td><td></td><td></td><td></td><td>—</td><td>—</td><td>—</td><td>—</td><td>—</td><td>6.8</td><td>6.9</td><td>6.7</td><td>6.8</td><td>6.5</td><td>6.6</td></tr><tr><td colspan="16"><strong>c) Wenchuan earthquake</strong></td></tr><tr><td>5</td><td>6.8</td><td></td><td></td><td></td><td>—</td><td>—</td><td>—</td><td>—</td><td>—</td><td>5.9</td><td>—</td><td>5.9</td><td>—</td><td>5.6</td><td>—</td></tr><tr><td>10</td><td>6.9</td><td></td><td></td><td></td><td>—</td><td>—</td><td>—</td><td>—</td><td>—</td><td>6.1</td><td>—</td><td>6.2</td><td>—</td><td>5.9</td><td>—</td></tr><tr><td>15</td><td>7.1</td><td></td><td></td><td></td><td>—</td><td>—</td><td>—</td><td>—</td><td>—</td><td>6.3</td><td>—</td><td>6.4</td><td>—</td><td>6.1</td><td>—</td></tr><tr><td>20</td><td>7.3</td><td></td><td></td><td></td><td>—</td><td>—</td><td>—</td><td>—</td><td>—</td><td>6.7</td><td>—</td><td>6.6</td><td>—</td><td>6.3</td><td>—</td></tr><tr><td>25</td><td>7.5</td><td></td><td></td><td></td><td>—</td><td>—</td><td>—</td><td>—</td><td>—</td><td>6.8</td><td>—</td><td>6.7</td><td>—</td><td>6.7</td><td>—</td></tr><tr><td>30</td><td>7.6</td><td></td><td></td><td></td><td>—</td><td>—</td><td>—</td><td>—</td><td>—</td><td>6.9</td><td>—</td><td>6.8</td><td>—</td><td>6.7</td><td>—</td></tr><tr><td>35</td><td>7.7</td><td></td><td></td><td></td><td>—</td><td>—</td><td>—</td><td>—</td><td>—</td><td>6.9</td><td>—</td><td>6.8</td><td>—</td><td>6.6</td><td>—</td></tr><tr><td>40</td><td>7.7</td><td></td><td></td><td></td><td>—</td><td>—</td><td>—</td><td>—</td><td>—</td><td>7.0</td><td>—</td><td>7.0</td><td>—</td><td>6.8</td><td>—</td></tr><tr><td>45</td><td>7.7</td><td></td><td></td><td></td><td>—</td><td>—</td><td>—</td><td>—</td><td>—</td><td>7.4</td><td>—</td><td>7.7</td><td>—</td><td>6.9</td><td>—</td></tr><tr><td>50</td><td>7.7</td><td></td><td></td><td></td><td>—</td><td>—</td><td>—</td><td>—</td><td>—</td><td>7.7</td><td>—</td><td>7.2</td><td>—</td><td>7.0</td><td>—</td></tr><tr><td>55</td><td>7.8</td><td></td><td></td><td></td><td>—</td><td>—</td><td>—</td><td>—</td><td>—</td><td>7.7</td><td>—</td><td>7.2</td><td>—</td><td>7.1</td><td>—</td></tr><tr><td>60</td><td>7.8</td><td></td><td></td><td></td><td>—</td><td>—</td><td>—</td><td>—</td><td>—</td><td>7.7</td><td>—</td><td>7.2</td><td>—</td><td>7.1</td><td>—</td></tr><tr><td>65</td><td>7.8</td><td></td><td></td><td></td><td>—</td><td>—</td><td>—</td><td>—</td><td>—</td><td>7.7</td><td>—</td><td>7.2</td><td>—</td><td>7.1</td><td>—</td></tr><tr><td>70</td><td>7.8</td><td></td><td></td><td></td><td>—</td><td>—</td><td>—</td><td>—</td><td>—</td><td>7.7</td><td>—</td><td>7.2</td><td>—</td><td>7.1</td><td>—</td></tr><tr><td>75</td><td>7.8</td><td></td><td></td><td></td><td>—</td><td>—</td><td>—</td><td>—</td><td>—</td><td>7.7</td><td>—</td><td>7.2</td><td>—</td><td>7.1</td><td>—</td></tr><tr><td>80</td><td>7.9</td><td></td><td></td><td></td><td>—</td><td>—</td><td>—</td><td>—</td><td>—</td><td>7.7</td><td>—</td><td>7.2</td><td>—</td><td>7.1</td><td>—</td></tr><tr><td>85</td><td>7.9</td><td></td><td></td><td></td><td>—</td><td>—</td><td>—</td><td>—</td><td>—</td><td>7.7</td><td>—</td><td>7.2</td><td>—</td><td>7.1</td><td>—</td></tr></tbody></table>

(Continued on following page)
<table><caption>**TABLE 2** | (Continued) Results from FinDer (rupture length, $L_{FD}$, and magnitude, $M_{FD}$), FinDerS (mean slip, Dmean, and magnitude, $M_{FDS}$, both from fitted slip profile) and FinDerS+ (predicted final rupture length, $L_{FDS+}$, predicted final mean slip, mean, and predicted final magnitude, $M_{FDS+}$) for the a) Ridgecrest, b) Kumamoto, and c) Wenchuan earthquake as a function of time from origin for a 60 km station distance cutoff. Magnitudes Mnr and $M_{FDS}$ are estimated from the USGS moment rate function and empirical PGD-distance relationships (Melgar et al., 2015; Crowell et al., 2016; Ruhl et al., 2019), respectively. Table shows updates in increments of 5 s.</caption><caption>(c) Wenchuan earthquake</caption><thead><tr><th></th><th>M<sub>mr</sub></th><th>L<sub>FD</sub> [km]</th><th>M<sub>FD</sub></th><th>D<sub>mean</sub> [m]</th><th>M<sub>FDS</sub></th><th>L<sub>FDS+</sub> [km]</th><th>D<sub>mean</sub> [m]</th><th>M<sub>FDS+</sub></th><th>M<sub>PGD</sub> Crowell</th><th>M<sub>PGD</sub> Melgar</th><th>M<sub>PGD</sub> Ruhl</th></tr></thead><tbody><tr><td>90</td><td>7.9</td><td>249</td><td>7.9</td><td>1.4</td><td>7.5</td><td>263</td><td>1.3</td><td>7.6</td><td>7.3</td><td>7.1</td><td>6.9</td></tr><tr><td>95</td><td>7.9</td><td>249</td><td>7.9</td><td>1.3</td><td>7.5</td><td>276</td><td>1.2</td><td>7.6</td><td>7.3</td><td>7.1</td><td>6.9</td></tr><tr><td>100</td><td>7.9</td><td>290</td><td>8.0</td><td>1.3</td><td>7.6</td><td>423</td><td>1.1</td><td>7.6</td><td>7.3</td><td>7.2</td><td>7.0</td></tr><tr><td>105</td><td>7.9</td><td>290</td><td>8.0</td><td>1.4</td><td>7.6</td><td>446</td><td>1.1</td><td>7.7</td><td>7.3</td><td>7.2</td><td>7.0</td></tr><tr><td>110</td><td>7.9</td><td>290</td><td>8.0</td><td>1.5</td><td>7.6</td><td>323</td><td>1.3</td><td>7.6</td><td>7.3</td><td>7.2</td><td>7.0</td></tr><tr><td>115</td><td>7.9</td><td>290</td><td>8.0</td><td>1.5</td><td>7.6</td><td>323</td><td>1.3</td><td>7.6</td><td>7.3</td><td>7.2</td><td>7.0</td></tr><tr><td>120</td><td>7.9</td><td>290</td><td>8.0</td><td>1.5</td><td>7.6</td><td>323</td><td>1.3</td><td>7.6</td><td>7.3</td><td>7.2</td><td>7.0</td></tr></tbody></table>

the fitting does not converge. The fitting resumes converging in the later stages, predicting the final rupture length fairly well.

### 4.2.4 Magnitudes

**Figure 8** and **Table 2** compare the magnitude estimates from the FinDer algorithms with magnitudes derived from the USGS moment rate function, $M_{mr}$, as well as from simple PGD-magnitude scaling relationships (Eq. 4). While these PGD-magnitude relations (Melgar et al., 2015; Crowell et al., 2016; Ruhl et al., 2019) were regressed originally from GNSS displacement data, we are applying them here to both seismic+GNSS and GNSS-only amplitudes. We compute $M_{PGD}$ for each of the three earthquakes by taking the median magnitude estimate taken over all stations at which the seismic waves should have arrived at a given time assuming a moveout velocity of 3 km/s.

For the $M_W$ 7.1 Ridgecrest earthquake (Figure 8A and Table 2), $M_{PGD}$ for both GNSS data and the combined seismic and GNSS datasets are the same for each of the three empirical scaling relationships, so we discuss them together. On the whole, $M_{PGD}$ under-estimates the magnitude of the Ridgecrest earthquake, and is similar to what is predicted by the FinDer family of algorithms. At the earliest time step when each FinDer algorithm is giving magnitudes, 15 s after $t_0$, the FinDer algorithms give magnitudes between $M_{FD}$ 6.5 and $M_{FDS+}$ 6.6, whereas the PGD-based magnitudes range between $M_{PGD}$ 6.3 (seismic+GNSS) and $M_{PGD}$ 6.7 (GNSS-only). The final estimate for magnitudes based on PGD ranges from $M_{PGD}$ 6.4 (seismic+GNSS) to $M_{PGD}$ 6.9 (GNSS-only), whereas the FinDer algorithms give magnitudes between $M_{FD}$ 6.9 to $M_{FDS+}$ 7.3. The range of FinDer algorithms give magnitudes closer to the moment rates, $M_{mr}$, cataloged by the USGS (Table 2).

PGD-based magnitude estimates for the $M_W$ 7.0 Kumamoto earthquake (Figure 8B and Table 2) continuously underestimate the true magnitude. For the first calculation of magnitude at 5 s, using Eq. 4, the values range from $M_{PGD}$ 4.2 (GNSS-only) to $M_{PGD}$ 5.4 (seismic+GNSS). For the same time step, the FinDer family of algorithms give $M_{FD}$ 5.4. At 10 s, the PGD based magnitudes range from $M_{PGD}$ 5.5 (GNSS-only) to $M_{PGD}$ 6.0 (seismic+GNSS), and FinDer magnitudes range from $M_{FDS}$ 6.2 to $M_{FD}$ 6.7. At 20 s, halfway through the time period, we calculate $M_{PGD}$ 6.4 (seismic+GNSS) to $M_{PGD}$ 6.8 (GNSS-only), however, the FinDer algorithms give larger values that are closer to the cataloged magnitude of $M_W$ 7.0, $M_{FDS}$ 6.9 to $M_{FD}$ 7.1. Finally, at 40 s, the PGD based magnitudes give a range of $M_{PGD}$ 6.5 (seismic+GNSS) to $M_{PGD}$ 6.9 (GNSS-only), and the FinDer algorithms give a range of $M_{FDS}$ 7.1 to $M_{FD}$ 7.2. It is notable that at first the combined geodetic and seismic PGD-based magnitude estimates are higher, and thus closer to the cataloged values. However, in the final time step, the geodetic (only) dataset gives higher PGD-based magnitudes.

For the $M_W$ 7.9 Wenchuan earthquake (Figure 8C and Table 2), we have only seismic records, so our PGD-based magnitude estimates are solely based on seismic data. Throughout the duration of the event, there is a relatively large discrepancy between $M_{PGD}$ and the true magnitude; the FinDer algorithms arrive at magnitude estimates much closer to
<figure><img src="image_11.png" /></figure>

A

<figure><img src="image_12.png" /></figure>

$M_{\text{FDS}}^{6.6}$

<figure><img src="image_13.png" /></figure>

$M_{\text{FDS+}}^{6.6}$

<figure><img src="image_14.png" /></figure>

B

<figure><img src="image_15.png" /></figure>

$M_{\text{FDS}}^{6.8}$

<figure><img src="image_16.png" /></figure>

$M_{\text{FDS+}}^{6.9}$

<figure><img src="image_17.png" /></figure>

C

<figure><img src="image_18.png" /></figure>

$M_{\text{FDS}}^{6.9}$

<figure><img src="image_19.png" /></figure>

$M_{\text{FDS+}}^{7}$

<figure><img src="image_20.png" /></figure>

D

<figure><img src="image_21.png" /></figure>

$M_{\text{FDS}}^{6.9}$

<figure><img src="image_22.png" /></figure>

$M_{\text{FDS+}}^{7.3}$

<figure><img src="image_23.png" /></figure>

**FIGURE 5** | Playback results for the $M_W$ 7.1 Ridgecrest earthquake (L ~ 50 km, ~ 20 s rupture duration) at a) 15 s, b) 20 s, c) 25 s, and d) 30 s from event origin. Plots on the left show FinDer line-source (black line) and PGD measurements at seismic (triangles) and geodetic (squares) sensors. Plots in the middle show interpolated slip profiles from FinDerS. Plots on the right show predicted slip profiles from FinDerS+. Epicenters are relative to the respective FinDer line-source model, oriented as indicated on the plots. See **Table 2** for details.
<figure><img src="image_24.png" /></figure>

<figure><img src="image_25.png" /></figure>

<figure><img src="image_26.png" /></figure>

<figure><img src="image_27.png" /></figure>

<figure><img src="image_28.png" /></figure>

<figure><img src="image_29.png" /></figure>

<figure><img src="image_30.png" /></figure>

<figure><img src="image_31.png" /></figure>

<figure><img src="image_32.png" /></figure>

<figure><img src="image_33.png" /></figure>

<figure><img src="image_34.png" /></figure>

<figure><img src="image_35.png" /></figure>

<figure><img src="image_36.png" /></figure>

**FIGURE 6** | Playback results for the M<sub>W</sub> 7.0 Kumamoto earthquake (L ~ 40 km, ~ 15 s rupture duration) at a) 10 s, b) 15 s, c) 20 s, and d) 25 s from event origin. Follows **Figure 5**. See **Table 2** for details.

that of the event. To begin with, at 20 s, once more than three
stations can be used to calculate the entire FinDer suite of
algorithms, the PGD based magnitude estimates range from

MP<sub>GD</sub> 6.3 to M<sub>PGD</sub> 6.7. Meanwhile, the FinDer algorithms
give a range of M<sub>FDS+</sub> 6.8 to M<sub>FD</sub> 7.0. At 60 s, the FinDer
suite of algorithms gives a range of M<sub>FDS</sub> 7.4 to M<sub>FD</sub> 7.7,
<figure><img src="image_37.png" /><figcaption>**FIGURE 7** | Playback results for the M<sub>W</sub> 7.9 Wenchuan earthquake (L ~ 300 km, ~ 100 s rupture duration) at a) 30 s, b) 40 s, c) 65 s, d) 90 s and e) 110 s from event origin. Follows **Figure 5**. See **Table 2** for details.</figcaption></figure>

while the PGD based magnitudes give a range between M<sub>PGD</sub> 6.9 and M<sub>PGD</sub> 7.2. The final spread of magnitude estimates at 120 s from PGD based estimates is M<sub>PGD</sub> 7.0 to M<sub>PGD</sub> 7.3, while the FinDer family of algorithms gives magnitude estimates much closer to the cataloged final magnitude M<sub>W</sub> 7.9, M<sub>FDS(+)</sub> 7.6 to M<sub>FD</sub> 8.0.

Generally, in terms of magnitude, the FinDer suite of algorithms performs quite well with respect to Mmr. FinDer continues to perform consistently for all three events, matching the M<sub>mr</sub> values for all three events, though for Kumamoto and Ridgecrest, FinDerS and FinDerS+ perform similarly well. The FinDerS+ magnitude estimates are

slightly better aside from the Wenchuan earthquake, where FinDer has the best results. This is mainly due to FinDerS+ under-estimating the Wenchuan slip significantly. The algorithms consistently outperform the PGD-based magnitude estimates.

# 5 DISCUSSION

Ideally, FinDerS+ is expected to provide the best performance because, on the one hand, it builds on the FinDer results which determine the fault location and strike, and on the other hand, it takes into account an important property of the source fault, its structural maturity gradient, in the form of the empirical equation Eq. 2. The latter allows integration of the location and strike of the fault delivered by FinDer with the slip estimates independently obtained from PGD amplitudes, and derives the best-fitting slip-length profile at every stage of the growing rupture. FinDerS provides a more basic approach as, at each stage of the rupture growth, it adopts the rupture length delivered by FinDer.

## 5.1 Performance for M<sub>W</sub> 7.1 Ridgecrest

Since March 2018, the U.S. West Coast ShakeAlert system (Given et al., 2018) is comprised of two independent algorithms - FinDer and the EPIC point-source algorithm (Chung et al., 2019). Both algorithms detected the 2019 M<sub>W</sub> 7.1 Ridgecrest shortly after its nucleation and sent an alert within 6.9 s from t<sub>0</sub>, even though they both under-estimated the magnitude of the mainshock by 0.8 units (Chung et al., 2020). In the case of FinDer this magnitude under-estimation was caused mostly by strongly increasing data latencies ( > 20–30 s) of seismic data telemetered from the stations to the main data center in Pasadena due to inefficient data compression and limited bandwidth (Stubailo et al., 2021). After the Ridgecrest earthquake, the FinDer code was updated to better handle latent data, as well to allow for a faster magnitude convergence (Böse et al., in prep.). We are presenting results from this new FinDer code here.

In our playback, FinDer determines the final rupture length of the M<sub>W</sub> 7.1 Ridgecrest earthquake within 25–30 s after t<sub>0</sub> as 53 km, with a close prediction of 45 km at 15–20 s (Table 2). FinDerS+, however, mildly over-estimates the final rupture length at 20–25 s (Figures 5B,C; Table 2), and grossly over-estimates the final rupture length at a value of 239 km at 30 s, when FinDer (and thus also FinDerS) has an accurate grip on the final rupture length. This length over-estimation with FinDerS+ is due to its inability to recover the earthquake slip asymmetry at 30 s.

While, for the most part, the slip profiles predicted through FinDerS and especially FinDerS+ match the overall pattern of the observed slip profile (Milliner and Donnellan, 2020), the backprojected slip values are under-estimated, with a maximum slip that is about half as large as measured at surface shortly after the earthquake (thus, these measures are not expected to include significant post-seismic slip, Barnhart et al., 2019; Wang and Bürgmann, 2020) (Figures 4A, 5A–D). This under-estimation may be a product of our approach using a very simple equation to estimate slip (Eq. 1), combined with
<figure><img src="image_38.png" /><figcaption>B 2016 Kumamoto, Mw 7.0</figcaption></figure>

<figure><img src="image_39.png" /><figcaption>C 2008 Wenchuan, Mw 7.9</figcaption></figure>

<figure><img src="image_40.png" /></figure>

<figure><img src="image_41.png" /></figure>

**FIGURE 8** | Comparison of magnitudes from FinDer, FinDerS, FinDerS+, and PGD-based magnitude estimates as function of time from event origin for a) M<sub>W</sub> 7.1 Ridgecrest, b) M<sub>W</sub> 7.0 Kumamoto, and c) M<sub>W</sub> 7.9 Wenchuan. Triangles mark estimates derived from seismic and geodetic data, squares mark estimates from geodetic data only. M<sub>mr</sub> is determined from the seismic moment rate and shown for reference. See **Table 2** for details.

using filtered PGD data, which is known to artificially reduce PGD amplitudes (Yamada, 2007). We come back to this point further below. The under-estimated slip values are likely the source of low magnitude predictions in FinDerS, and of the difficulty of FinDerS+ to recover the asymmetry of the final slip profile and hence, its actual (shorter) final length.

The main output of the FinDer algorithm is a line-source model, not a magnitude estimate (Böse et al., 2018). As a matter of

fact, while the FinDer rupture length estimate of 53 km is excellent, the final FinDer magnitude M<sub>FD</sub> 6.9 under-estimates the cataloged magnitude by 0.2 units. This is mainly due to the use of a simple empirical rupture length-magnitude relation (Wells and Coppersmith, 1994) in FinDer; although valuable, this relation has significant uncertainties. Using the slip information in addition to rupture length is expected to compensate for limitations in the empirical error-prone
relation used in FinDer. In the Ridgecrest earthquake, FinDerS, however, also never arrives at the cataloged magnitude, $M_W$ 7.1 (Figures 5A–D). FinDerS+, while over-predicting the final rupture length, does provide a better magnitude prediction. This, however, comes at the cost of a largely over-estimated rupture length, in turn likely related to slip under-estimation. This raises a fundamental question about trade-offs in algorithms, where accuracy in one area (magnitude) may mean a sacrifice in other metrics (final rupture length). In the future, it may be useful to combine the family of FinDer algorithms to find a magnitude and a length prediction that considers all three values.

## 5.2 Performance for $M_W$ 7.0 Kumamoto

The Japanese EEW system operated by Japan Meteorological Agency (JMA) issued a first alert 8 s after $t_0$ of the $M_W$ 7.0 Kumamoto earthquake with an initial magnitude estimate of $M_{JMA}$ 5.9. About 5 s later, the magnitude was updated to $M_{JMA}$ 6.9 (Kodera et al., 2016). In comparison, 5 s after $t_0$ (neglecting data latencies, which typically are on the order of 1–2 s), FinDer gives $M_{FD}$ 5.4. Then, 5 s later, FinDer updates its estimate to $M_{FD}$ 6.7, and to $M_{FD}$ 7.0 after 15 s from $t_0$.

For the Kumamoto earthquake, the family of FinDer algorithms predicts the final rupture length at surface from 10 s after $t_0$ (Table 2), and the final rupture length at depth at 15 s. Subsequently, FinDer, as well as FinDerS and FinDerS+ (due to their reliance on FinDer for the rupture length), over-estimate and over-predict the rupture length of the Kumamoto earthquake by ~20 km. However, the final rupture length of the event determined by FinDer(S) is 84 km, which is approximately the distance from the hypocenter to the triggered event at Yufuin. Since this secondary event did not occur until > 30 s after $t_0$ of the original Kumamoto earthquake, the results do not reflect the combined energy from both events until after ~35 s. While we do not have an explanation, we suggest it is possible that this finding (which was already observed in (Böse et al., 2018)) is not simply fortuitous; the mainshock originally had enough energy to rupture, in one way or another (i.e., in one or two earthquakes), the total fault length of 84 km. This is supported by the anomalously strong ground motions towards NE direction, which might have been a product of strong directivity, the simultaneous slipping of two oblique faults with strong pulse like ruptures, and high slip rates (Kobayashi et al., 2017; Somei et al., 2019). These motions may have been strong enough to break the fault entirely, although in two parts somewhat delayed (dynamic triggering of the Yufuin event Uchide et al., 2016; Yoshida, 2016).

Looking at slip distributions up to ~40 km along strike (the final rupture length at the ground surface), the maximum slip values calculated through backprojection at 25 s are consistent with the largest surface slip measures. However, overall, the backprojected slips are about half that measured at surface (Shirahama et al., 2016) (Figure 4B), while the asymmetry of the slip profile is recovered at 20 s only.

## 5.3 Performance for $M_W$ 7.9 Wenchuan

At the time of the $M_W$ 7.9 Wenchuan earthquake, which is by far the largest event in this study and hence has the longest rupture

duration, China had no operational EEW system. Thus, we cannot compare the performance of the FinDer suite of algorithms to another early warning system. However, in a previous study using the FinDer algorithm and the same dataset, Li et al. (2020) reported a slightly slower magnitude convergence for the Wenchuan earthquake, which can be attributed to their use of an older version of the FinDer code.

The FinDer algorithms all anticipate from the very start (25–30 s after $t_0$) that the rupture will be long, at least 100–150 km, and the magnitude large, at least 7.2 (Table 2). FinDer and FinDerS then predict an increasing length, up to the accurate estimate of the final rupture length by 95–100 s after $t_0$ (Table 2). While FinDerS+ over-predicts the length between about 45 and 85 s, it predicts it well from 90 s on, even though its final estimate is slightly greater than the actual rupture length. As for the Ridgecrest and the Kumamoto earthquakes, the backprojected slip values are under-estimated, in this case significantly with the maximum slip about three times lower than the largest displacements measured at surface, and the mean backprojected slip about 5–6 times smaller than what was observed. The actual asymmetry of the slip profile is not well recovered either, even though FinDerS+ anticipates a westward asymmetry from 45 s.

In terms of magnitude, FinDer predicts a MFD 7.8 by 60–65 s after $t_0$, which is approximately 50% of the rupture duration. At 100 s, FinDer reaches a final magnitude of $M_{FD}$ 8.0, which is a reasonable estimate. Compared to Ridgecrest and Kumamoto, seismic recordings of the Wenchuan earthquake are sparse, which partially explains the longer duration needed for source characterization, along with the generally expected longer duration of such a large earthquake (Table 1).

## 5.4 Seismic/Geodetic Displacement Versus Fault Slip

Long-period PGD amplitudes provide information about fault slip and magnitudes (e.g. Melgar et al., 2015; Crowell et al., 2016; Ruhl et al., 2019). As demonstrated in this study, both GNSS and seismic data can produce PGD datasets that can help to quickly constrain the metrics of an ongoing earthquake. Although more reliable at long periods and in particular at the static offset, GNSS data is generally more noisy compared to strong-motion data. The GNSS noise level is typically 1–2 cm and data is thus useful only for moderate and large earthquakes at local to regional distances (Ruhl et al., 2019).

It is well known that recovering seismic displacement from double-integrated strong-motion data can present issues in the presence of strong nonphysical drifts and saturation of the sensor (Bock et al., 2011). These drifts typically emerge from small baseline offsets in the acceleration data that are accentuated when double-integrated to displacement (Emore et al., 2007). Since this issue mainly affects long-period motions (and we are mainly interested in dynamic rather than static displacement), we apply in this study a high-pass filter with a cutoff-frequency at 0.075 Hz. A drawback of this approach, however, is that PGD amplitudes are reduced (Yamada, 2007), which will be reflected in smaller slip values resulting from the backprojection of PGD
using **Eq. 1**. Nonetheless, first-order differences between seismic and geodetic dynamic PGD amplitudes are small (Figure 3) and the utility of the filtered amplitudes is clearly demonstrated in this paper.

Contrary to the results of the same algorithm using accurate static slip profiles, where final rupture lengths and magnitudes could be determined within 10–20% of the way through the signal’s duration (Hutchison et al., 2020), the results from the dynamic application of FinDerS+ in this paper suggest that the true convergence might actually take longer. Instead, as shown by the examples of the earthquakes presented here, FinDerS and FinDerS+ rarely reach the cataloged magnitude until 40–50% of the rupture duration. For Wenchuan, the actual magnitude is never reached, though both algorithms do converge at a reasonable solution (within 0.2–0.5 magnitude units) about 30% of the way through the event. This slower ability to predict the final length and magnitude of the earthquake arises from two issues: first, it takes some time until final slip values in an earthquake are reached. Second, the (final) slip tends to be significantly under-estimated in our backprojection algorithm. This difficulty to reproduce the slip values accurately leads in some cases to over-predict the rupture length (see Figures 5D, 7C). While this over-prediction might not be too much an issue in large earthquakes as those analyzed here, it would become more critical for smaller magnitude earthquakes. The approach developed here is thus best suited for earthquakes with $M_W$ greater than 7. However, more systematic studies (and adjustments of the algorithms) are required to determine the best range of applicability.

The systematic under-estimation of slip from the backprojected PGD amplitudes largely relates to the function **Eq. 1** we are applying. This function, taken from Yamada (2007), was derived from near-source ground motion simulations by Aagaard et al. (2004), made to relate dynamic ground displacements at some distance from the fault rupture to on-rupture slip. These simulations were done using reference rupture source models, having variable dips and slip modes, but similar length of 100 km and similar largest slips of about 6 m. Would these ruptures be real, available earthquake slip-length scaling relations (Manighetti et al., 2007) would suggest that these reference earthquakes occurred on fairly mature faults. In contrast, the three earthquakes we analyzed here ruptured fairly immature faults, as described in **Section 3.1**. Their slip-to-length ratios actually confirm this inference (compare their slip-length values to **Figure 5A** in (Manighetti et al., 2007)). Earthquakes on immature faults have been shown to produce larger displacements and stronger ground accelerations than earthquakes on more mature faults (Radiguet et al., 2009). As the function from Aagaard et al. (2004) does not take these source differences into account, it likely smooths them out, in effect under-estimating the slips produced by earthquakes on fairly immature faults, as those analyzed here.

To address the problem of slip under-estimation in the future, we might consider an updated backprojection function that takes the overall maturity of the ruptured fault into account, provided that the latter can be known in real-time (Böse and Heaton, 2010). Alternatively, if the ruptured fault is not identified or if its overall maturity is unknown, in future we might be able to derive from

analysis of the increasing number of well instrumented earthquakes, a “generic” coefficient of correction that would compensate for the systematic slip under-estimation for immature faults. As a first estimate, for the three earthquakes analyzed here, a factor of two is suggested. By using simultaneously the “raw” and the corrected slip data, the FinDer algorithms might be able to bracket more efficiently, i.e., faster, the actual final rupture length and magnitude.

## 5.5 FinDer-Versus PGD-Based Magnitudes

For the three earthquakes studied here, the FinDer family of algorithms systematically outperforms the simple PGD-based magnitude estimates, both using exclusively GNSS data or combined GNSS and seismic datasets. While the PGD-hypocentral distance scaling relationship is simple, the fact that the coefficients in **Eq. 4** change repeatedly based on the introduction of new earthquakes (Melgar et al., 2015; Crowell et al., 2016; Ruhl et al., 2019), and that station selectivity can have a significant impact on the accuracy of the magnitude estimate (e.g. Dahmen et al., 2020), make this scaling relationship imprecise. Furthermore, these PGD-based magnitude relations ignore the earthquake fault source dimensions and properties, while those have been shown to impact both the ground accelerations (Radiguet et al., 2009) and displacements (Manighetti et al., 2007). The FinDer group of algorithms consistently does a better job of calculating magnitude accurately, without the introduction of complex station selection parameters or subsequently extending empirical datasets. FinDerS and FinDerS+, which take some of the earthquake fault source properties into account, are generally more efficient at predicting the final magnitudes (Figures 8A,B).

Our PGD-based magnitudes for the Kumamoto earthquake (Figure 8) differ from the results of a previous GNSS data-based study (Dahmen et al., 2020), in which magnitude predictions are much closer to the cataloged magnitude. For example, Dahmen et al. (2020) estimate $M_{PGD}$ 6.3 14 s after the event detection, $M_{PGD}$ 6.8 after 20 s, and finally $M_{PGD}$ 7.0 after 100 s. However, there are a couple of differences compared to our study. First, at least five stations that were included in the Dahmen et al. (2020) study, were not included in the GNSS dataset available through Ruhl et al. (2019) that we used here. Second, Dahmen et al. (2020) utilized a noise-based threshold criteria to select stations. Because of this, Dahmen et al. (2020) used only 27 of the 39 available GNSS stations that we used (Ruhl et al., 2019). Without applying this selection criteria, the PGD-based magnitude estimate of Dahmen et al. (2020) 14 s after the event detection is $M_{PGD}$ 6.7, which is similar to our estimates. It is notable that - despite not applying station selection criteria - the FinDer algorithm suite magnitude estimates are consistent with the PGD-based magnitude estimates from Dahmen et al. (2020) and the true magnitude evolution reflected in $M_{mr}$ (Table 2). The same applies for the $M_W$ 7.1 Ridgecrest earthquake, where the $M_{PGD}$ values of Hodgkinson et al. (2020) differ slightly from our results.

## 5.6 Implications and Benefits for EEW

The main output of the FinDer algorithm is a line-source model, which is determined from the spatial distribution of
high-frequency PGA amplitudes, FinDer magnitudes, M<sub>FD</sub>, are a secondary product only, that are estimated from the application of general empirical rupture length-magnitudes relationships. These relations are known to have significant uncertainties. Our two new algorithms, FinDerS and FinDerS+, can help to improve magnitude estimates by incorporating the additional information on the earthquake slip-fault maturity relation and on long-period motions, which are closely related to the static fault offset and seismic moment (Hanks and Kanamori, 1979).

Furthermore, the new addition of a fault slip profile estimated by FinDerS and FinDerS+ can improve the spatial prediction of the ground motions, the ultimate goal of EEW. Slip distributions with large and smooth slip patches affect mainly long-period ground motions along the fault rupture and thus matter most for high-rise buildings at close distance as well as tsunami generation (e.g. Scholz, 2019). Furthermore, by associating the greatest coseismic slip with the most mature parts of the rupture, we can extrapolate the locations of the strongest radiation of high-frequency ground motions. It has indeed been shown that the amplitude of the ground motions is partly controlled by the structural maturity of the source fault, with immature faults or fault sections producing the largest ground motions (Radiguet et al., 2009), furthermore, likely, of highest frequencies. This implies that the output of FinDerS and FinDerS+ does not help only with magnitude estimates (at least for redundancy checks), but potentially improves the seismic ground-motion predictions for EEW.

As proposed by Hutchison et al. (2020), we could even go further by deriving and using two extreme slip-length-maturity empirical functions (Eq. 1 in Hutchison et al. (2020)), one describing the slip behavior of earthquakes on mature faults, and one that of earthquakes on immature faults (see Manighetti et al., 2007). By using simultaneously these two extreme equations, the FinDer algorithms would better bracket the actual range of lengths and magnitudes and hence probably converge faster towards the actual final length and magnitude of the earthquake.

# 6 CONCLUSION

We recover earthquake line-source models from high-frequency seismic acceleration data using FinDer (Böse et al., 2012b, 2015, 2018) and give 1D slip profile estimates from the backprojection of long-period displacement amplitudes onto these models using two new algorithms, FinDerS and FinDerS+. While FinDerS backprojects displacement onto the current FinDer line-source, FinDerS+ allows for additional rupture growth along this model. Both FinDerS and FinDerS+ make use of a general relationship that we established in an earlier study (Hutchison et al., 2020) that relates the along-strike gradient of long-term structural maturity of the ruptured fault to the coseismic slip distribution along the rupture. As such, FinDerS and FinDerS+ incorporate a part of the earthquake physics described through this empirical relation, which we use to produce an informed interpolation and smoothing of the backprojected slip values. The two new algorithms provide independent magnitude estimates from

FinDer based on slip which could be useful to speed-up magnitude convergence and for redundancy checks with FinDer estimated magnitudes.

FinDerS+ can over-predict the rupture length in two cases: 1) when the slip data are not well determined, as is the case here, in particular showing sparse measures with abrupt fluctuations as in Ridgecrest; 2) when the ongoing slip is gradually increasing away from the hypocenter [see Hutchison et al. (2020)]. This is because the empirical relation between earthquake slip, rupture length, and gradient in fault structural maturity relies on the envelope shape of the along-fault slip distribution and the generic asymmetry of the earthquake slip profiles [see Manighetti et al. (2005)]. Therefore, in the cases where slip fluctuations are too large, or the major slip drop occurs at the end of the rupture growth, the model is unable to recover the envelope shape of the slip distribution. To overcome the problem of slip fluctuations (in the presence of multiple large asperities), dense and accurate slip data are needed. To overcome the problem of the few earthquake cases where slip gradually increasing away from the hypocenter (those cases are fairly rare), in future work, we may incorporate the overall maturity of the source faults: the earthquake slip-length curves have indeed an amplitude that decreases with the overall maturity of the ruptured fault (Manighetti et al., 2007). These amplitude relations will provide slip maximums and force tighter constraints on upper slip bounds of the evolving slip curves, especially for these very asymmetric earthquake profiles. Even though our approach over-predicts the final length of the Ridgecrest and Wenchuan earthquakes at specific stages of their growth, it does predict their magnitude well, from fairly early stages of the rupture growth. In playback, these early large magnitudes indicate that these earthquakes will end up being large, even during the initial stages of the rupture.

For the earthquakes shown here (Kumamoto, Ridgecrest, Wenchuan) the magnitudes computed by the FinDer algorithms converge faster and reach values closer to the cataloged magnitudes than compared to those computed directly from PGD. FinDerS and FinDerS+ benefit from real-time GNSS data streams, but could also use seismic data only. This study has been a preliminary work to understand the basic utility of two new FinDer-based algorithms, but further tests are needed to better constrain the capabilities and limitations of these new EEW algorithms. Further work is also needed to improve our capacity to estimate ongoing slip in real-time. The accuracy of slip estimation is critical to make FinDerS and FinDerS+ efficient. Furthermore, larger sets of earthquakes need to be analyzed, including subduction earthquakes for which information about slip is especially critical in a tsunami context, and rupture length of vital importance to warn the populations along the coastline.

# DATA AVAILABILITY STATEMENT

The original contributions presented in the study are included in the article/Supplementary Material, further inquiries can be directed to the corresponding author.
## AUTHOR CONTRIBUTIONS

MB: gave idea and guided research, computed FinDer models, wrote manuscript; AH: conducted analysis, co-developed FinDerS(+) algorithm, wrote manuscript; IM: provided fault background/interpretation, wrote manuscript, leads the project in which this work was done; FM: helped with discussions and writing manuscript; JL: provided the strong ground motion records of the Wenchuan earthquake, assisted in analysis of Wenchuan earthquake; JC: helped with discussions and writing manuscript.

## FUNDING

This project is funded by the ANR Grant FAULTS_R_GEMS #ANR-17-CE31-0008.

## REFERENCES

Aagaard, B. T., Hall, J. F., and Heaton, T. H. (2004). Effects of Fault Dip and Slip Rake Angles on Near-Source Ground Motions: Why Rupture Directivity Was Minimal in the 1999 Chi-Chi, Taiwan, Earthquake. *Bull. Seismological Soc. America* 94, 155–170. doi:10.1785/0120030053

Allen, R. M., Gasparini, P., Kamigaichi, O., and Bose, M. (2009). The Status of Earthquake Early Warning Around the World: An Introductory Overview. *Seismological Res. Lett.* 80, 682–693. doi:10.1785/gssrl.80.5.682

Allen, R. M., and Melgar, D. (2019). Earthquake Early Warning: Advances, Scientific Challenges, and Societal Needs. *Annu. Rev. Earth Planet. Sci.* 47, 361–388. doi:10.1146/annurev-earth-053018-060457

Barnhart, W. D., Hayes, G. P., and Gold, R. D. (2019). The July 2019 Ridgecrest, California, Earthquake Sequence: Kinematics of Slip and Stressing in Cross-Fault Ruptures. *Geophys. Res. Lett.* 46, 11859–11867. doi:10.1029/2019gl084741

Behr, Y., Clinton, J., Kästli, P., Cauzzi, C., Racine, R., and Meier, M.-A. (2015). Anatomy of an Earthquake Early Warning (EEW) Alert: Predicting Time Delays for an End-To-End EEW System. *Seismological Res. Lett.* 86, 830–840. doi:10.1785/0220140179

Blaser, L., Krüger, F., Ohrnberger, M., and Scherbaum, F. (2010). Scaling Relations of Earthquake Source Parameter Estimates with Special Focus on Subduction Environment. *Bull. Seismological Soc. America* 100, 2914–2926. doi:10.1785/0120100111

Bock, Y., Melgar, D., and Crowell, B. W. (2011). Real-time strong-motion Broadband Displacements from Collocated GPS and Accelerometers. *Bull. Seismological Soc. America* 101, 2904–2925. doi:10.1785/0120110007

Böse, M., Felizardo, C., and Heaton, T. H. (2015). Finite-Fault Rupture Detector (FinDer): Going Real-Time in CalifornianShakeAlertWarning System. *Seismological Res. Lett.* 86, 1692–1704. doi:10.1785/0220150154

Böse, M., Heaton, T., and Hauksson, E. (2012a). Rapid Estimation of Earthquake Source and Ground-Motion Parameters for Earthquake Early Warning Using Data from a Single Three-Component Broadband or strong-motion Sensor. *Bull. Seismological Soc. America* 102, 738–750. doi:10.1785/0120110152

Böse, M., Heaton, T. H., and Hauksson, E. (2012b). Real-time Finite Fault Rupture Detector (FinDer) for Large Earthquakes. *Geophys. J. Int.* 191, 803–812. doi:10.1111/j.1365-246x.2012.05657.x

Böse, M., and Heaton, T. H. (2010). Probabilistic Prediction of Rupture Length, Slip and Seismic Ground Motions for an Ongoing Rupture: Implications for Early Warning for Large Earthquakes. *Geophys. J. Int.* 183, 1014–1030. doi:10.1111/j.1365-246x.2010.04774.x

Böse, M., Heaton, T., and Hudnut, K. (2013). Combining Real-Time Seismic and GPS Data for Earthquake Early Warning. *AGU Fall Meet. Abstr.* 2013, G51B–G05.

Böse, M., Smith, D. E., Felizardo, C., Meier, M.-A., Heaton, T. H., and Clinton, J. F. (2018). FinDer v.2: Improved Real-Time Ground-Motion Predictions for M2-M9 with Seismic Finite-Source Characterization. *Geophys. J. Int.* 212, 725–742. doi:10.1093/gji/ggx430

## ACKNOWLEDGMENTS

This material is based on services provided by UNAVCO, Inc., the National Research Institute for Earth Science and Disaster Resilience (NIED), the GEONET Global Positioning System, the China Strong MotionNetworks Center (CSMNC), and the Southern California Earthquake Data Center (SCEDC) and Southern California Seismic Network (SCSN). The SCEDC and SCSN are funded through U.S. Geological Survey Grant G10AP00091.

## SUPPLEMENTARY MATERIAL

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/feart.2021.685879/full#supplementary-material

Chen, Y., and Booth, D. C. (2011). *The Wenchuan Earthquake of 2008: Anatomy of a Disaster*. Springer Science & Business Media.

Chung, A. I., Henson, I., and Allen, R. M. (2019). Optimizing Earthquake Early Warning Performance: ElarmS-3. *Seismological Res. Lett.* 90, 727–743. doi:10.1785/0220180192

Chung, A. I., Meier, M.-A., Andrews, J., Böse, M., Crowell, B. W., McGuire, J. J., et al. (2020). Shakealert Earthquake Early Warning System Performance during the 2019 Ridgecrest Earthquake Sequence. *Bull. Seismological Soc. America* 110, 1904–1923. doi:10.1785/0120200032

Clinton, J., Zollo, A., Marmureanu, A., Zulfikar, C., and Parolai, S. (2016). State-of-the Art and Future of Earthquake Early Warning in the European Region. *Bull. Earthquake Eng.* 14, 2441–2458. doi:10.1007/s10518-016-9922-7

Colombelli, S., Zollo, A., Festa, G., and Picozzi, M. (2014). Evidence for a Difference in Rupture Initiation between Small and Large Earthquakes. *Nat. Commun.* 5, 3958. doi:10.1038/ncomms4958

Crowell, B. W., Schmidt, D. A., Bodin, P., Vidale, J. E., Gomberg, J., Renate Hartog, J., et al. (2016). Demonstration of the Cascadia G-FAST Geodetic Earthquake Early Warning System for the Nisqually, Washington, Earthquake. *Seismological Res. Lett.* 87, 930–943. doi:10.1785/0220150255

Cua, G., Fischer, M., Heaton, T., and Wiemer, S. (2009). Real-time Performance of the Virtual Seismologist Earthquake Early Warning Algorithm in Southern California. *Seismological Res. Lett.* 80, 740–747. doi:10.1785/gssrl.80.5.740

Cua, G., and Heaton, T. H. (2009). Characterizing Average Properties of Southern California Ground Motion Amplitudes and Envelopes.

Dahmen, N., Hohensinn, R., and Clinton, J. (2020). Comparison and Combination of GNSS and Strong-Motion Observations: A Case Study of the 2016 Mw 7.0 Kumamoto Earthquake. *Bull. Seismological Soc. America* 110, 2647–2660. doi:10.1785/0120200135

Denolle, M. A. (2019). Energetic Onset of Earthquakes. *Geophys. Res. Lett.* 46, 2458–2466. doi:10.1029/2018gl080687

Emore, G. L., Haase, J. S., Choi, K., Larson, K. M., and Yamagiwa, A. (2007). Recovering Seismic Displacements through Combined Use of 1-Hz GPS and strong-motion Accelerometers. *Bull. Seismological Soc. America* 97, 357–378. doi:10.1785/0120060153

Given, D. D., Allen, R. M., Baltay, A. S., Bodin, P., Cochran, E. S., Creager, K., et al. (2018). Revised Technical Implementation Plan for the ShakeAlert System–An Earthquake Early Warning System for the West Coast of the United States. *U.S. Geological Survey Open-File Report 2018–1155* [Supersedes USGS Open-File Report 2014–1097], 42. doi:10.3133/ofr20181155

Goldberg, D. E., Melgar, D., Bock, Y., and Allen, R. M. (2018). Geodetic Observations of Weak Determinism in Rupture Evolution of Large Earthquakes. *J. Geophys. Res. Solid Earth* 123, 9950–9962. doi:10.1029/2018jb015962

Hanks, T. C., and Kanamori, H. (1979). A Moment Magnitude Scale. *J. Geophys. Res.* 84, 2348–2350. doi:10.1029/jb084ib05p02348

Hodgkinson, K. M., Mencin, D. J., Feaux, K., Sievers, C., and Mattioli, G. S. (2020). Evaluation of Earthquake Magnitude Estimation and Event Detection
Thresholds for Real-Time GNSS Networks: Examples from Recent Events Captured by the Network of the Americas. *Seismological Res. Lett.* 91, 1628–1645. doi:10.1785/0220190269

Hudnut, K. W., Brooks, B. A., Scharer, K., Hernandez, J. L., Dawson, T. E., Oskin, M. E., et al. (2020). Airborne Lidar and Electro-Optical Imagery along Surface Ruptures of the 2019 Ridgecrest Earthquake Sequence, Southern California. *Seismological Res. Lett.* 91, 2096–2107. doi:10.1785/0220190338

Hutchison, A. A., Böse, M., and Manighetti, I. (2020). Improving Early Estimates of Large Earthquake's Final Fault Lengths and Magnitudes Leveraging Source Fault Structural Maturity Information. *Geophys. Res. Lett.* 47, e2020GL087539. doi:10.1029/2020gl087539

Kobayashi, H., Koketsu, K., and Miyake, H. (2017). Rupture Processes of the 2016 Kumamoto Earthquake Sequence: Causes for Extreme Ground Motions. *Geophys. Res. Lett.* 44, 6002–6010. doi:10.1002/2017gl073857

Kodera, Y., Saitou, J., Hayashimoto, N., Adachi, S., Morimoto, M., Nishimae, Y., et al. (2016). Earthquake Early Warning for the 2016 Kumamoto Earthquake: Performance Evaluation of the Current System and the Next-Generation Methods of the Japan Meteorological Agency. *Earth, Planets and Space* 68, 1–14. doi:10.1186/s40623-016-0567-1

Kodera, Y., Yamada, Y., Hirano, K., Tamaribuchi, K., Adachi, S., Hayashimoto, N., et al. (2018). The Propagation of Local Undamped Motion (PLUM) Method: A Simple and Robust Seismic Wavefield Estimation Approach for Earthquake Early Warning. *Bull. Seismological Soc. America* 108, 983–1003. doi:10.1785/0120170085

Leonard, M. (2014). Self-consistent Earthquake Fault-Scaling Relations: Update and Extension to Stable continental Strike-Slip Faults. *Bull. Seismological Soc. America* 104, 2953–2965. doi:10.1785/0120140087

Li, J., Böse, M., Wyss, M., Wald, D. J., Hutchison, A., Clinton, J. F., et al. (2020). Estimating Rupture Dimensions of Three Major Earthquakes in Sichuan, China, for Early Warning and Rapid Loss Estimates. *Bull. Seismological Soc. America* 110, 920–936. doi:10.1785/0120190117

Manighetti, I., Campillo, M., Bouley, S., and Cotton, F. (2007). Earthquake Scaling, Fault Segmentation, and Structural Maturity. *Earth Planet. Sci. Lett.* 253, 429–438. doi:10.1016/j.epsl.2006.11.004

Manighetti, I., Campillo, M., Sammis, C., Mai, P., and King, G. (2005). Evidence for Self-Similar, Triangular Slip Distributions on Earthquakes: Implications for Earthquake and Fault Mechanics. *J. Geophys. Res. Solid Earth* 110. doi:10.1029/2004jb003174

Manighetti, I., King, G. C. P., Gaudemer, Y., Scholz, C. H., and Doubre, C. (2001). Slip Accumulation and Lateral Propagation of Active normal Faults in Afar. *J. Geophys. Res.* 106, 13667–13696. doi:10.1029/2000jb900471

Mattioli, G. S., Phillips, D. A., Hodgkinson, K. M., Walls, C., Mencin, D. J., Bartel, B. A., et al. (2020). The GAGE Data and Field Response to the 2019 Ridgecrest Earthquake Sequence. *Seismological Res. Lett.* 91, 2075–2086. doi:10.1785/0220190283

Meier, M.-A., Ampuero, J. P., and Heaton, T. H. (2017). The Hidden Simplicity of Subduction Megahrust Earthquakes. *Science* 357, 1277–1281. doi:10.1126/science.aan5643

Melgar, D., Crowell, B. W., Geng, J., Allen, R. M., Bock, Y., Riquelme, S., et al. (2015). Earthquake Magnitude Calculation without Saturation from the Scaling of Peak Ground Displacement. *Geophys. Res. Lett.* 42, 5197–5205. doi:10.1002/2015gl064278

Melgar, D., and Hayes, G. P. (2019). Characterizing Large Earthquakes before Rupture Is Complete. *Sci. Adv.* 5, eaav2032. doi:10.1126/sciadv.aav2032

Milliner, C., Bürgmann, R., Inbal, A., Wang, T., and Liang, C. (2020). Resolving the Kinematics and Moment Release of Early Afterslip within the First Hours Following the 2016 M<sub>w</sub> 7.1 Kumamoto Earthquake: Implications for the Shallow Slip Deficit and Frictional Behavior of Aseismic Creep. *J. Geophys. Res. Solid Earth* 125, e2019JB018928. doi:10.1029/2019jb018928

Milliner, C., and Donnellan, A. (2020). Using Daily Observations from Planet Labs Satellite Imagery to Separate the Surface Deformation between the 4 July Mw 6.4 Foreshock and 5 July Mw 7.1 Mainshock during the 2019 Ridgecrest Earthquake Sequence. *Seismological Res. Lett.* 91, 1986–1997. doi:10.1785/0220190271

Perrin, C., Manighetti, I., Ampuero, J.-P., Cappa, F., and Gaudemer, Y. (2016a). Location of Largest Earthquake Slip and Fast Rupture Controlled by Along-Strike Change in Fault Structural Maturity Due to Fault Growth. *J. Geophys. Res. Solid Earth* 121, 3666–3685. doi:10.1002/2015jb012671

Perrin, C., Manighetti, I., and Gaudemer, Y. (2016b). Off-fault Tip Splay Networks: A Genetic and Generic Property of Faults Indicative of Their Long-Term Propagation. *Comptes Rendus Geosci.* 348, 52–60. doi:10.1016/j.crte.2015.05.002

Radiguet, M., Cotton, F., Manighetti, I., Campillo, M., and Douglas, J. (2009). Dependency of Near-Field Ground Motions on the Structural Maturity of the Ruptured Faults. *Bull. Seismological Soc. America* 99, 2572–2581. doi:10.1785/0120080340

Ruhl, C. J., Melgar, D., Geng, J., Goldberg, D. E., Crowell, B. W., Allen, R. M., et al. (2019). A Global Database of Strong-Motion Displacement GNSS Recordings and an Example Application to PGD Scaling. *Seismological Res. Lett.* 90, 271–279. doi:10.1785/0220180177

Satriano, C., Elia, L., Martino, C., Lancieri, M., Zollo, A., and Iannaccone, G. (2011). PRESTo, the Earthquake Early Warning System for Southern Italy: Concepts, Capabilities and Future Perspectives. *Soil Dyn. Earthquake Eng.* 31, 137–153. doi:10.1016/j.soildyn.2010.06.008

Scholz, C. H. (2019). *The Mechanics of Earthquakes and Faulting*. Cambridge University Press.

Scott, C. P., Arrowsmith, J. R., Nissen, E., Lajoie, L., Maruyama, T., and Chiba, T. (2018). The M 7 2016 Kumamoto, Japan, Earthquake: 3-D Deformation along the Fault and within the Damage Zone Constrained from Differential Lidar Topography. *J. Geophys. Res. Solid Earth* 123, 6138–6155. doi:10.1029/2018jb015581

Shen, Z.-K., Sun, J., Zhang, P., Wan, Y., Wang, M., Bürgmann, R., et al. (2009). Slip Maxima at Fault Junctions and Rupturing of Barriers during the 2008 Wenchuan Earthquake. *Nat. Geosci* 2, 718–724. doi:10.1038/ngeo0636

Shirahama, Y., Yoshimi, M., Awata, Y., Maruyama, T., Azuma, T., Miyashita, Y., et al. (2016). Characteristics of the Surface Ruptures Associated with the 2016 Kumamoto Earthquake Sequence, central Kyushu, Japan. *Earth, Planets and Space* 68, 1–12. doi:10.1186/s40623-016-0559-1

Somei, K., Miyakoshi, K., Yoshida, K., Kurahashi, S., and Irikura, K. (2020). Near-Source Strong Pulses During Two Large MJMA 6.5 and MJMA 7.3 Events in the 2016 Kumamoto, Japan, Earthquakes. *Pure Appl. Geophys.*, 177. doi:10.1007/s00024-019-02095-6

Subailo, I., Alvarez, M., Biasi, G., Bhadha, R., and Hauksson, E. (2021). Latency of Waveform Data Delivery from the Southern California Seismic Network during the 2019 Ridgecrest Earthquake Sequence and its Effect on ShakeAlert. *Seismological Soc. America* 92, 170–186. doi:10.1785/0220200211

Uchide, T., Horikawa, H., Nakai, M., Matsushita, R., Shigematsu, N., Ando, R., et al. (2016). The 2016 Kumamoto-Oita Earthquake Sequence: Aftershock Seismicity gap and Dynamic Triggering in Volcanic Areas. *Earth, Planets and Space* 68, 1–10. doi:10.1186/s40623-016-0556-4

Wang, K., and Bürgmann, R. (2020). Co- and Early Postseismic Deformation Due to the 2019 Ridgecrest Earthquake Sequence Constrained by Sentinel-1 and COSMO-SkyMed SAR Data. *Seismological Res. Lett.* 91, 1998–2009. doi:10.1785/0220190299

Wells, D. L., and Coppersmith, K. J. (1994). New Empirical Relationships Among Magnitude, Rupture Length, Rupture Width, Rupture Area, and Surface Displacement. *Bull. seismological Soc. America* 84, 974–1002.

Wesnousky, S. G. (1988). Seismological and Structural Evolution of Strike-Slip Faults. *Nature* 335, 340–343. doi:10.1038/335340a0

Wu, Y.-M., Yen, H.-Y., Zhao, L., Huang, B.-S., and Liang, W.-T. (2006). Magnitude Determination Using Initial P Waves: A Single-Station Approach. *Geophys. Res. Lett.* 33. doi:10.1029/2005gl025395

Yamada, M. (2007). *Early Warning for Earthquakes with Large Rupture Dimension*.

Yoshida, S. (2016). Earthquakes in Oita Triggered by the 2016 M 7.3 Kumamoto Earthquake. *Earth, Planets and Space* 68, 1–13. doi:10.1186/s40623-016-0552-8

**Conflict of Interest:** The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

**Publisher's Note:** All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Copyright © 2021 Böse, Hutchison, Manighetti, Li, Massin and Clinton. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.