Abstract
Slow slip transients on faults can last from seconds to months and stitch together the earthquake cycle. However, no single geophysical instrument is able to observe the full range of slow slip because of bandwidth limitations. Here, we connect seismic and geodetic data from the Mexican subduction zone to explore an instrumental blind spot. We establish a calibration of the daily median amplitude of the seismically recorded low-frequency earthquakes to the daily geodetically recorded moment rate of previously established slow slip events. This calibration allows us to use the precise evolution of low-frequency earthquake activity to quantitatively measure the moment of smaller, subdaily slip events that are unresolvable by geodesy alone. The resulting inferred slow slip moments scale with duration and inter-event time like ordinary earthquakes. These new quantifications help connect slow and fast events in a broad spectrum of transient slip and suggest that slow slip events behave much like ordinary earthquakes.
INTRODUCTION
Slow slip is a common feature of some of Earth’s most dangerous subduction zones (1). These gradual transients are distinct from ordinary earthquakes and demonstrate that plate motion between major earthquakes is not as steady as previously thought, but is rather composed of a rich spectrum of events. A complete continuum of slip modes has been speculated to exist, with event durations from a fraction of a second to many months (2). Geodetic observations of slow slip on plate boundaries span years with continuous GPS records, but episodic slip on shorter time scales is largely invisible with standard geodetic techniques. On time scales of several minutes to hours, tremor (3) accompanies slow slip and contains seconds-long repetitive earthquakes called low-frequency earthquakes (LFEs) that are indicative of record transient slip (4, 5). What is missing is a quantitative record of motion between the seismic and geodetically recorded events. Can slow slip occur on minute to hour time scales, and if so, how does its behavior compare to previously established relationships for better studied, ordinary earthquakes?
Here, we present a calibration between seismic and geodetic data recorded above a subduction zone and use it to demonstrate that slow slip occurs every day on this plate boundary. We analyze a 2.5-year time period in the Guerrero (Mexico) segment that contains millions of seismically detected LFEs (6) and eight geodetically constrained slow slip events of different sizes (7, 8). We first exploit LFE amplitudes to identify otherwise-undetectable slow transients and measure their seismic moment rates
OBSERVATIONS
Seismic moment rate of LFEs during slow slip
Earthquake displacement amplitudes in the farfield are proportional to the seismic moment rate
We then explore in Fig. 1 how LFE amplitudes vary with the slow slip that is geodetically resolvable and dominates the tectonic slip cycle in Guerrero. Every 4 years, the regional GPS network captures a 6-month-long slow slip event that releases as much built-up tectonic stress as an M7.5 earthquake (7, 10). Further downdip, smaller M6.4 slow slip events have previously been identified in the GPS record using recurring bursts of low-frequency seismicity as a guide (8, 11). We observe the strongest LFE amplitudes in the updip source region, where LFE activity is most strongly modulated by slow slip (6). Notably, LFE amplitudes systematically increase during major slow slip events, implying that the LFE seismic moment rate
(A) Median S-wave LFE displacement amplitude in 7-week time bins and 15-km spatial bins, normalized by the mean of all plotted median amplitudes; color scale is logarithmic. (B) North-south surface displacement observed at the MEZC GPS station during the 2.5-year LFE record. Colored patches mark geodetically observed slow slip events: the 2006 M7.5 slow slip events [red; (7)] and the M6.4 slow slip events [yellow; (8)]. The dark red patches indicate the intermittent slow transients during the 2006 slow slip event (12), highlighting the short slip duration compared to the 6-month event duration. (C) Schematic of the Guerrero subduction zone, with inset showing geographical location. The two LFE source regions are indicated by the blue and purple boxes. The red and orange patches correspond to areas where geodetically observed slow slip occurs.
We capitalize on this observation by focusing on the 58 LFE sources within the updip source region in Guerrero (fig. S2) to estimate the average LFE
Geodetic moment rate of slow slip during large slow slip events
With estimates of the average seismic moment rate
Because slow slip is most prominent in the geodetic record during LFE activity (12, 13), we suggest that the duration of LFE activity is a good proxy for the active slow slip period. We thus use the number of LFE per day NLFE as a proxy for slow slip duration. We estimate
Identification and measurement of additional slow slip events
The power of combining the LFEs with the geodetic data is in the ability to detect otherwise invisible slow slip events. We therefore gather together GPS data on days with similar median LFE amplitudes to find additional events. GPS position solutions are daily; therefore, higher temporal resolution is not possible for this dataset. We bin by amplitude as finely as possible while minimizing overlap and maintaining a similar number of data points in each bin (see Materials and Methods). Each amplitude bin corresponds to a slow transient that generates a given range of LFE moment rates (fig. S1). The surface displacement rate of each newly identified slow transient is measured as the velocity of the cumulative sum of the binned GPS displacements, following the method of Frank (13). The geodetic moment
Calibration of seismic to geodetic moment rates
The resulting moment rate estimates in Fig. 2 appear to follow a power-law trend of the form
Units of δT reflect the average duration of the aseismic pulses that drive LFE activity; we assume δT = 0.5 s (14, 15) for the right y axis. The squares indicate the moment rates of geodetically observed slow slip events. The blue circles reflect 1-day slow transients, representing the average moment rate for a given range of LFE amplitudes (fig. S1). See Materials and Methods for a discussion of error bars.
We note that the slow transients derived from the LFE analysis are critical to this calibration. If only the eight geodetically observed slow slip events in Fig. 2 (orange and red symbols) are considered, then the relationship between seismic and geodetic moment rate is not obvious. The largest M7.5 slow slip exhibits higher than average geodetic and seismic moment rates. The small slow slip event of February 2006 displays the lowest moment rates; this slow slip lasts five times longer than the other six M6.4 events and exhibits a flat GPS signature alongside small LFE amplitudes (Fig. 1). However, most geodetically observed slow slip events are tightly clustered around fairly similar moment rates, and LFEs are necessary to identify a fuller range of dynamics.
DISCUSSION
Daily record of slow slip inferred from LFEs
Now that we have a well-calibrated scaling of seismic and geodetic moment rates, we can use this relationship to estimate the slow transients that drive LFE activity during every day of the seismic record. We do this by converting the daily median LFE amplitude into geodetic moment rate and then by multiplying the individual event slip duration to find the total daily slip (T = NLFEδT). Figure 3B shows the resulting catalog of slow transients. We observe a M > 4.5 slow earthquake every day along the subduction interface beneath Guerrero; at least 6 × 1015 N-m of slow slip is released by the plate interface each day. The total moment released over the 2.5-year study period is equivalent to an M7.64 earthquake. Previous studies documented direct (13, 17) and indirect (18) evidence that slow slip occurs more often than is obvious in the geodetic record. Here, we demonstrate that slow slip occurs on subdaily time scales, and we are able to measure the moment released every day.
The yellow and red patches respectively highlight the geodetically observed small M6.4 and large M7.5 slow slip events (Fig. 1). The dark red patches indicate the intermittent slow transients during the 2006 slow slip event (12). (A) Waiting times between successive LFEs in the updip LFE source region highlight the clustered seismicity that accompanies slow slip (29). (B) Evolution of daily slow transient magnitudes, estimated with the moment rate scaling shown in Fig. 2. The small slow slip events are made up of a single group of slow transients, while the 2006 M7.5 slow slip event is composed of a complex cluster of these transients.
Subdaily slow transients have been sought before by using time series of repeating earthquakes as proxies for fault slip (19–21). The repeating earthquake slip proxy depends strongly on an assumed interaction between the seismic and aseismic fault processes. We do not need to rely on that assumption, because the seismic-geodetic calibration allows us to measure the daily slow slip moment for the first time.
The resulting calibration is similar in spirit to the calibration between tremor and geodetic moment rate suggested by Aguiar et al. (22). This previous study, however, used tremor duration rather than LFE amplitude and focused only on large, relatively long duration events. That work inferred a constant moment rate for geodetic slip, which is inconsistent with the empirical scaling of Fig. 2. With the help of LFEs, we are able to sample high moment rates during slow slip events that would otherwise be lost to averaging over a highly variable, intermittent slip process.
The large geodetically observed slow slip events in Fig. 3 appear to be made up of many shorter ~M6 transients. Most of the M6.4 slow slip events appear as single sequences of these large slow transients. The 2006 M7.5 slow slip event is an intermittent sequence of M6 slow transients (12). The characteristic size of these slow transients echoes the recently reported characteristic rates of tectonic tremor activity during both major and small slow slip events (23). The magnitude-frequency distribution of slow transients is not well represented by an earthquake-like power law and rather appears to have a characteristic size (fig. S4). However, the LFE amplitudes in isolation do follow a power law (fig. S1).
Moment scaling of daily slow transients
We now go further to examine the systematics of the newly discovered events. Past studies of slow earthquakes have suggested that their slip dynamics are fundamentally different than those of earthquakes (24). In previous work, a constant apparent moment rate implied that slow earthquake moment Mo scales linearly with duration T (25). However, Fig. 4 shows that slow transient moment Mo scales with the slip duration T cubed. This slow transient moment-duration scaling is the same scaling relationship observed for earthquakes (26), implying that slow slip and earthquakes are more similar than previously thought. We speculate that the previously reported Mo ∝ T scaling (24) reflects the interaction of two distinct processes: the slow rupture process and the clustering process that links disparate slow transients into a large slow slip event (12). GPS observations are unable to provide the necessary resolution to distinguish between these two mechanisms, and capture the combination of both processes. Once we strip away this clustering process from slow slip events to examine how moment scales with slip duration, slow slip is revealed to follow the same scaling as fast earthquake slip. We can only speculate on the physical significance of the exponent of 3 in Eq. 1, and it remains to be seen whether it varies from one plate boundary to the next.
(A) Slow transient moment Mo scales with the slip duration T3 (blue line), where duration is the number of LFEs multiplied by the average aseismic pulse duration of δT; we assume δT = 0.5 s (14, 15) for the right y axis. The moment-duration scaling is the same as for ordinary earthquakes. (B) Slow transient moment scales with the inter-event interval τ to the power of 5.5 (blue line). The moment–inter-event interval scaling is approximately the same as for repeating earthquakes (orange line).
At least one additional systematic is in the data. We define a proxy for the average slow transient inter-event interval τ based on the longest waiting times between LFEs each day; this quantity represents the amount of time per day that slow slip is not happening. We estimate τ as the sum of the 10% longest LFE waiting times. This τ is similar to a recurrence time, but we prefer the more general term “inter-event time” to avoid any implication of periodic behavior. Figure 4 shows that slow transient moment increases with the inter-event time and scales with the same exponent of 6 as previously observed for repeating earthquakes (19, 27). The observation suggests a commonality stemming from both processes involving seismic asperities surrounded by aseismic ruptures with varying degrees of partitioning between seismic and aseismic slip.
Conclusions
We have demonstrated here that slow transients occur every day by constraining the seismic-geodetic moment rate scaling of slow slip, linking the dynamics of low-frequency seismicity to the geodetic signature of the driving aseismic slip. The existence of these minutes-long daily transients fills in an observational gap, suggesting that the spectrum of slip is continuous (2, 28). Accounting for the clustered and intermittent rupture process of slow slip, we find that the moment-duration scaling of these daily slow transients is the same as ordinary earthquake scaling, and the moment-recurrence scaling follows the same relationship as repeating earthquakes. In aggregate, our observations suggest that slow slip on a plate interface is similar to ordinary earthquakes, and a broad continuum of transient slip governed by one set of dynamics may explain both phenomena.
MATERIALS AND METHODS
Estimating LFE displacement amplitudes and seismic moment rate
We measured the displacement amplitude A of all observed LFE S waves on each horizontal component of the 10 recording seismic stations (6). First, we measured the RMS amplitude in a 6-s window centered on the LFE S wave in the frequency band of 1 to 2 Hz. This window duration matches the average length of LFE S waves in Guerrero, and the frequency band contains the highest signal-to-noise ratios in these events (6). We corrected for geometric spreading of the form r−1, where r is the source-receiver distance and then averaged the resulting, corrected amplitudes over all stations and components.
Displacement amplitude is directly proportional to seismic moment rate (9)
LFE amplitude binning to derive slow transients
We selected the LFE amplitude bins to contain at least 30 data points per bin while minimizing overlap. The four bins chosen for Fig. 2, from lowest to highest moment rates, contained 52, 41, 47, and 35 days of data. Following the method of Frank (13), the GPS displacements in each bin were cumulatively summed, and the average displacement per day was estimated as the slope of the best-fit linear trend.
To investigate the sensitivity of the results to binning, we performed the same analysis as described in the main text with three rather than four amplitude bins. We used the same criteria described above to select three amplitude bins, except that we required a minimum of 40 data points in each bin. We observed a power-law relationship between seismic and geodetic moment rates in fig. S3 with scaling coefficients of α = 1013.89 ± 0.26 and β = 3.15 ± 0.40. This trend is similar to the relationship constrained in Fig. 2.
Error estimation in seismic-geodetic moment rate scaling
The geodetic moment rate errors in Fig. 2 for the geodetically observed slow slip events reflect the observational surface displacement errors. These errors reflect the uncertainty of the estimated static displacement offsets that control the inversion for the distribution of slow slip on the interface. The displacement error for the 2006 M7.5 slow slip event represents ~10% of the estimated moment (7). The average displacement error for the seven M6.4 slow slip events represents ~25% of the estimated moment (8). The uncertainties of the amplitude-derived slow transients were determined with a jackknife analysis. We performed a 10% jackknife of the LFEs in each amplitude bin, determined the average surface displacement per day following the method of Frank (13), and estimated the moment as the corresponding fault dislocation (16). We estimated the geodetic moment rate error of each slow transient as the 10th and 90th percentiles of the resulting jackknifed distribution of estimated moments.
The seismic moment rate errors were estimated with a similar jackknife analysis. We performed a 10% jackknife of the LFEs in each slow slip event or slow transient and recomputed the median LFE amplitude. We estimated the error bars as the 10th to 90th percentiles of the resulting jackknifed distribution of median LFE amplitudes.
The errors in the estimated power-law coefficients α and β were determined by performing a one-sample jackknife of the 12 moment rate data points (8 geodetically constrained slow slip events and 4 amplitude-derived slow transients) and then fitting a new power-law relationship to each jackknifed resampling. We estimated the errors of α and β as the standard deviation of the jackknifed distribution of coefficients.
SUPPLEMENTARY MATERIALS
Supplementary material for this article is available at http://advances.sciencemag.org/cgi/content/full/5/10/eaaw9386/DC1
Fig. S1. Distribution and evolution of LFE displacement amplitudes in Guerrero.
Fig. S2. Tectonic context of the subduction zone underneath Guerrero, Mexico.
Fig. S3. Alternative number of slow transients to constrain the seismic to geodetic moment rate relationship shown in Fig. 2.
Fig. S4. Distribution of slow transient magnitudes.
Fig. S5. Daily count of the 5% largest LFEs (2337 events); the plotted amplitudes are >12.4 nm.
Reference (30)
This is an open-access article distributed under the terms of the Creative Commons Attribution license, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.
REFERENCES AND NOTES
- Copyright © 2019 The Authors, some rights reserved; exclusive licensee American Association for the Advancement of Science. No claim to original U.S. Government Works. Distributed under a Creative Commons Attribution License 4.0 (CC BY).