Seismic Coupling Research Articles

Systematic simulation of multicomponent receiver coupling to the seafloor using rheological models

This paper reports a comparison of three different rheological models used to characterize receiver coupling to the seafloor. We used a finite-element simulation tool to simulate the mechanical receiver coupling to the seafloor as a viscoelastic system with a combination of linear elastic springs and linear viscous dashpots, known as rheological models. Three models cover most of all mechanic coupling systems, the most commonly applied Kelvin-Voigt model (KVM), the Maxwell model (MM), and the standard linear solid (SLS) model. The models differ in behavior for different coupling aspects such as oscillation, creeping, stress relaxation, and their combinations. We tested these models’ ability and relevance for use in modeling seismic receiver coupling to the seafloor. For that purpose, we used an optimized mathematical approach to simulate coupling behavior under various coupling conditions. We found how receiver coupling will affect P- and S-waves for all three models and provided some insight into which model is most suitable to describe coupling under different circumstances. We found that the SLS model represents a general description of most of the coupling effects to the seafloor and should be used when the coupling acts as a viscoelastic system. The KVM and MM are applicable in extreme cases, such as for elastic waves in consolidated sediments (KVM) and dominant creeping effects, as in very soft biosediment (MM).

Along-strike variations in temperature and tectonic tremor activity along the Hikurangi subduction zone, New Zealand

In the Hikurangi subduction zone, situated along the east coast of the North Island, New Zealand, where the old oceanic Pacific Plate is subducting beneath the Australian Plate, several slow slip events and tectonic tremors have recently been documented. These observations are somewhat surprising because such slow seismic phenomena tend to be common in subduction zones where relatively young oceanic plate is subducting. The locations of tectonic tremors, down-dip limit of slow slip events and seismic coupling transition change along strike from greater depths in the south to shallower depths in the north, suggesting significant along-strike variations in the characteristics of the plate interface. Similar along-strike variations have been observed for other characteristic features of the Hikurangi subduction zone. Here, we demonstrate that along-strike variations observed for tectonic tremors, slow slip events, and seismic coupling can be explained by lateral differences in the thermal structure of the subduction zone, which are controlled mainly by variations in convergence rate and friction along the plate interface. To demonstrate this, we first confirm that tectonic tremors occur around the plate interface. Then, we calculate the thermal structure of the Hikurangi subduction zone using a two-dimensional finite difference code. To explain the along-strike variation in the heat flow observed in the forearc region, temperatures along the plate interface should be systematically higher in the northern region than in the southern region, which we interpret as a consequence of higher convergence rates and greater frictional heating in the northern region. We compare the along-strike variation of seismic characteristics with calculated thermal structure and highlight that this along-strike variation in temperature controls the depth of the brittle-ductile transition, which is consistent with the observed spatial variations in tectonic tremors, down-dip limit of slow slip events and seismic coupling. Our results suggest that tectonic tremors recorded within subduction zones reflect the transient rheology of the materials being subducted, which is controlled by variations in temperature along the plate interface.

Comparison of postseismic afterslip models with aftershock seismicity for three subduction-zone earthquakes: Nias 2005, Maule 2010 and Tohoku 2011

We focus on the relation between seismic and total postseismic afterslip following the Maule Mw 8.8 earthquake on 2010 February 27 in central Chile. First, we calculate the cumulative slip released by aftershock seismicity. We do this by summing up the aftershock regions and slip estimated from scaling relations. Comparing the cumulative seismic slip with afterslip models we show that seismic slip of individual aftershocks exceeds locally the inverted afterslip model from geodetic constraints. As the afterslip model implicitly contains the displacements from the aftershocks, this reflects the tendency of afterslip models to smear out the actual slip pattern. However, it also suggests that locally slip for a number of the larger aftershocks exceeds the aseismic slip in spite of the fact that the total equivalent moment of the afterslip exceeds the cumulative moment of aftershocks by a large factor. This effect, seen weakly for the Maule 2010 and also for the Tohoku 2011 earthquake, can be explained by taking into account the uncertainties of the seismicity and afterslip models. In spite of uncertainties, the hypocentral region of the Nias 2005 earthquake is suggested to release a large fraction of moment almost purely seismically. Therefore, these aftershocks are not driven solely by the afterslip but instead their slip areas have probably been stressed by interseismic loading and the mainshock rupture. In a second step, we divide the megathrust of the Maule 2010 rupture into discrete cells and count the number of aftershocks that occur within 50 km of the centre of each cell as a function of time. We then compare this number to a time-dependent afterslip model by defining the ‘afterslip to aftershock ratio’ (ASAR) for each cell as the slope of the best fitting line when the afterslip at time t is plotted against aftershock count. Although we find a linear relation between afterslip and aftershocks for most cells, there is significant variability in ASAR in both the downdip and along-strike directions of the megathrust. We compare the spatial distribution of ASAR with the spatial distribution of seismic coupling, coseismic slip and Bouguer gravity anomaly, and in each case we find no significant correlation.

Gradual unlocking of plate boundary controlled initiation of the 2014 Iquique earthquake.

On 1April 2014, Northern Chile was struck by a magnitude 8.1 earthquake following a protracted series of foreshocks. The Integrated Plate Boundary Observatory Chile monitored the entire sequence of events, providing unprecedented resolution of the build-up to the main event and its rupture evolution. Here we show that the Iquique earthquake broke a central fraction of the so-called northern Chile seismic gap, the last major segment of the South American plate boundary that had not ruptured in the past century. Since July 2013 three seismic clusters, each lasting a few weeks, hit this part of the plate boundary with earthquakes of increasing peak magnitudes. Starting with the second cluster, geodetic observations show surface displacements that can be associated with slip on the plate interface. These seismic clusters and their slip transients occupied a part of the plate interface that was transitional between a fully locked and a creeping portion. Leading up to this earthquake, the bvalue of the foreshocks gradually decreased during the years before the earthquake, reversing its trend a few days before the Iquique earthquake. The mainshock finally nucleated at the northern end of the foreshock area, which skirted a locked patch, and ruptured mainly downdip towards higher locking. Peak slip was attained immediately downdip of the foreshock region and at the margin of the locked patch. We conclude that gradual weakening of the central part of the seismic gap accentuated by the foreshock activity in a zone of intermediate seismic coupling was instrumental in causing final failure, distinguishing the Iquique earthquake from most great earthquakes. Finally, only one-third of the gap was broken and the remaining locked segments now pose a significant, increased seismic hazard with the potential to host an earthquake with a magnitude of >8.5.

Wet Fault or Dry Fault? A Laboratory Approach to Remotely Monitor the Hydro-Mechanical State of a Discontinuity Using Controlled-Source Seismics

Stress variation and fluid migration occur in deformation zones, which are expected to affect seismic waves reflected off or propagating across such structures. We developed a basic experimental approach to monitor the mechanical coupling with respect to seismic coupling across a single discontinuity between a granite sample in contact with a steel platen. Piezoceramics located on the platen were used to both generate and record the P and S wave fields reflected off the discontinuity at normal incidence. This way, normal (B n ) and tangential (B t ) compliances were calculated using Schoenberg’s linear slip theory (Schoenberg, J Acoust Soc Am 68:1516–1521, 1980) when the roughness, the effective pressure (P eff, up to 200 MPa), and the nature of the filling (gas or water) vary. We observe that increasing the effective pressure decreases B n and B t , which is interpreted as the effect of the closure of the voids at the interface, permitting more seismic energy to be transmitted across the interface. Values of B n are significantly higher than those of B t at low P eff (<60–80 MPa) in dry conditions, and significantly drop under water-saturated conditions. The water filling the voids therefore helps to transmit the seismic energy of compressional waves across the interface. These results show that the assumption B n ≈ B t commonly found in some theoretical approaches does not always stand. The ratio B n /B t actually reflects the type of saturating fluids and the effective pressure, in agreement with other experimental studies. However, we illustrate that only the relative variations of this ratio seem to be relevant, not its absolute value as suggested in previous studies. Consequently, the use of B n against B t plots may allow effective pressure variation and the nature of the pore fluid to be inferred. In this respect, this experimental approach at sample scale helps to pave the way for remotely monitoring in the field the hydro-mechanical state of deformation zones, such as seismogenic faults, fractured reservoirs, or lava conduits.

Discrete characteristics of the aftershock sequence of the 2011 Van earthquake

An intraplate earthquake of magnitude Mw 7.2 occurred on a NE–SW trending blind oblique thrust fault in accretionary orogen, the Van region of Eastern Anatolia on October 23, 2011. The aftershock seismicity in the Van earthquake was not continuous but, rather, highly discrete. This shed light on the chaotic nonuniformity of the event distribution and played key roles in determining the seismic coupling between the rupturing process and seismogeneity.I analyzed the discrete statistical mechanics of the 2011 Van mainshock–aftershock sequence with an estimation of the non-dimensional tuning parameters consisting of; temporal clusters (C) and the random (RN) distribution of aftershocks, range of size scales (ROSS), strength change (εD), temperature (T), P-value of temporal decay, material parameter R-value, seismic coupling χ, and Q-value of aftershock distribution. I also investigated the frequency–size (FS), temporal (T) statistics and the sequential characteristics of aftershock dynamics using discrete approach and examined the discrete evolutionary periods of the Van earthquake Gutenberg–Richter (GR) distribution. My study revealed that the FS and T statistical properties of aftershock sequence represent the Gutenberg–Richter (GR) distribution, clustered (C) in time and random (RN) Poisson distribution, respectively.The overall statistical behavior of the aftershock sequence shows that the Van earthquake originated in a discrete structural framework with high seismic coupling under highly variable faulting conditions. My analyses relate this larger dip–slip event to a discrete seismogenesis with two main components of complex fracturing and branching framework of the ruptured fault and dynamic strengthening and hardening behavior of the earthquake. The results indicate two dynamic cases. The first is associated with aperiodic nature of aftershock distribution, indicating a time-independent Poissonian event. The second is associated with variable slip model, varying loading rates and high seismic coupling along the ruptured area. This study concludes that the aperiodic statistical behavior of the Van aftershock sequence is one of the results of the time-independent Poissonian seismic deformation in the 2011 Van earthquake.

A trial derivation of seismic plate coupling by focusing on the activity of shallow slow earthquakes

To understand the effect of plate coupling on very low-frequency event (VLFE) activity resulting from megathrust earthquakes, we performed long-term multiscale earthquake cycle simulations (including a megathrust earthquake and slow earthquakes) on a 3-D subduction plate boundary model, based on a rate- and state-dependent friction law. Our simulation suggests that quiescence of shallow VLFEs off Miyagi may be explained by the location in the shallow central part of the 2011 Tohoku earthquake because of the locally strong coupling, while observed activation of VLFEs off Iwate (northern part of Tohoku district), Fukushima (southern part of Tohoku district), and Ibaraki (northern part of Kanto district) is explained by the location on the outer rim. The area and duration of the quiescence off Miyagi may be a new clue to evaluate the potential for plate coupling strong enough to cause the next megathrust earthquake.

Seismic attenuation tomography of the Northeast Japan arc: Insight into the 2011 Tohoku earthquake (Mw9.0) and subduction dynamics

Detailed three-dimensional (3-D) P and S wave attenuation (Qp and Qs) models of the crust and upper mantle under the entire Northeast Japan (Tohoku) arc from the Japan Trench to the Japan Sea coast are determined, for the first time, using a large number of high-quality t* data measured precisely from P and S wave spectra of local earthquakes. The suboceanic earthquakes used in this work are relocated precisely using sP depth phases. Our results reveal a prominent landward dipping high-Q zone representing the subducting Pacific slab, a landward dipping intermediate- to high-Q zone in the mantle wedge between the Pacific coast and the volcanic front, and significant low-Q anomalies in the crust and mantle wedge between the volcanic front and the Japan Sea coast. Prominent high-Q patches surrounded by low-Q anomalies are revealed in the interplate megathrust zone under the Tohoku fore arc where the great 2011 Tohoku-oki earthquake (Mw 9.0) occurred. The high-Q patches in the megathrust zone generally exhibit large coseismic slips of megathrust earthquakes and large slip deficit on the plate interface. We think that these high-Q patches represent asperities in the megathrust zone, whereas the low-Q anomalies reflect weakly coupled areas. We also find that the hypocenters of the 2011 Tohoku-oki interplate earthquakes (Mw > 7.0) are located in areas where Qp, Qs, and Qp/Qs change abruptly. These results suggest that structural heterogeneities in the megathrust zone control the interplate seismic coupling and the nucleation of megathrust earthquakes.

A high-resolution, time-variable afterslip model for the 2010 Maule Mw = 8.8, Chile megathrust earthquake

The excellent spatial coverage of continuous GPS stations in the region affected by the Maule Mw = 8.8 2010 earthquake, combined with the proximity of the coast to the seismogenic zone, allows us to model megathrust afterslip on the plate interface with unprecedented detail. We invert post-seismic observations from continuous GPS sites to derive a time-variable model of the first 420 d of afterslip. We also invert co-seismic GPS displacements to create a new co-seismic slip model. The afterslip pattern appears to be transient and non-stationary, with the cumulative afterslip pattern being formed from afterslip pulses. Changes in static stress on the plate interface from the co- and post-seismic slip cannot solely explain the aftershock patterns, suggesting that another process – perhaps fluid related – is controlling the lower magnitude aftershocks. We use aftershock data to quantify the seismic coupling distribution during the post-seismic phase. Comparison of the post-seismic behaviour to interseismic locking suggests that highly locked regions do not necessarily behave as rate-weakening in the post-seismic period. By comparing the inter-seismic locking, co-seismic slip, afterslip, and aftershocks we attempt to classify the heterogeneous frictional behaviour of the plate interface.

Megathrust splay faults at the focus of the Prince William Sound asperity, Alaska

AbstractHigh‐resolution sparker and crustal‐scale air gun seismic reflection data, coupled with repeat bathymetric surveys, document a region of repeated coseismic uplift on the portion of the Alaska subduction zone that ruptured in 1964. This area defines the western limit of Prince William Sound. Differencing of vintage and modern bathymetric surveys shows that the region of greatest uplift related to the 1964 Great Alaska earthquake was focused along a series of subparallel faults beneath Prince William Sound and the adjacent Gulf of Alaska shelf. Bathymetric differencing indicates that 12 m of coseismic uplift occurred along two faults that reached the seafloor as submarine terraces on the Cape Cleare bank southwest of Montague Island. Sparker seismic reflection data provide cumulative Holocene slip estimates as high as 9 mm/yr along a series of splay thrust faults within both the inner wedge and transition zone of the accretionary prism. Crustal seismic data show that these megathrust splay faults root separately into the subduction zone décollement. Splay fault divergence from this megathrust correlates with changes in midcrustal seismic velocity and magnetic susceptibility values, best explained by duplexing of the subducted Yakutat terrane rocks above Pacific plate rocks along the trailing edge of the Yakutat terrane. Although each splay fault is capable of independent motion, we conclude that the identified splay faults rupture in a similar pattern during successive megathrust earthquakes and that the region of greatest seismic coupling has remained consistent throughout the Holocene.

Large earthquake rupture process variations on the Middle America megathrust

The megathrust fault between the underthrusting Cocos plate and overriding Caribbean plate recently experienced three large ruptures: the August 27, 2012 (Mw 7.3) El Salvador; September 5, 2012 (Mw 7.6) Costa Rica; and November 7, 2012 (Mw 7.4) Guatemala earthquakes. All three events involve shallow-dipping thrust faulting on the plate boundary, but they had variable rupture processes. The El Salvador earthquake ruptured from about 4 to 20 km depth, with a relatively large centroid time of ∼19 s, low seismic moment-scaled energy release, and a depleted teleseismic short-period source spectrum similar to that of the September 2, 1992 (Mw 7.6) Nicaragua tsunami earthquake that ruptured the adjacent shallow portion of the plate boundary. The Costa Rica and Guatemala earthquakes had large slip in the depth range 15 to 30 km, and more typical teleseismic source spectra. Regional seismic recordings have higher short-period energy levels for the Costa Rica event relative to the El Salvador event, consistent with the teleseismic observations. A broadband regional waveform template correlation analysis is applied to categorize the focal mechanisms for larger aftershocks of the three events. Modeling of regional wave spectral ratios for clustered events with similar mechanisms indicates that interplate thrust events have corner frequencies, normalized by a reference model, that increase down-dip from anomalously low values near the Middle America trench. Relatively high corner frequencies are found for thrust events near Costa Rica; thus, variations along strike of the trench may also be important. Geodetic observations indicate trench-parallel motion of a forearc sliver extending from Costa Rica to Guatemala, and low seismic coupling on the megathrust has been inferred from a lack of boundary-perpendicular strain accumulation. The slip distributions and seismic radiation from the large regional thrust events indicate relatively strong seismic coupling near Nicoya, Costa Rica, patchy zones of strong seismic coupling in the shallowest megathrust region along Nicaragua and El Salvador, and small deeper patchy zones of strong seismic coupling near Guatemala, which can be reconciled with the geodetic observations as long as the strong coupling is limited to a small fraction of the megathrust area.

Along-trench structural variation and seismic coupling in the northern Japan subduction zone

Large destructive interplate earthquakes, such as the 2011 Mw 9.0 Tohoku-oki earthquake, have occurred repeatedly in the northern Japan subduction zone. The spatial distribution of large interplate earthquakes shows distinct along-trench variations, implying regional variations in interplate coupling. We conducted an extensive wide-angle seismic survey to elucidate the along-trench variation in the seismic structure of the forearc and to examine structural factors affecting the interplate coupling beneath the forearc mantle wedge. Seismic structure models derived from wide-angle traveltimes showed significant along-trench variation within the overlying plate. In a weakly coupled segment, (i) the sediment layer was thick and flat, (ii) the forearc upper crust was extremely thin, (iii) the forearc Moho was remarkably shallow (about 5 km), and (iv) the P-wave velocity within the forearc mantle wedge was low, whereas in the strongly coupled segments, opposite conditions were found. The good correlation between the seismic structure and the segmentation of the interplate coupling implies that variations in the forearc structure are closely related to those in the interplate coupling.

Seismicity patterns along the Ecuadorian subduction zone: new constraints from earthquake location in a 3-D a priori velocity model

To improve earthquake location, we create a 3-D a priori P-wave velocity model (3-DVM) that approximates the large velocity variations of the Ecuadorian subduction system. The 3-DVM is constructed from the integration of geophysical and geological data that depend on the structural geometry and velocity properties of the crust and the upper mantle. In addition, specific station selection is carried out to compensate for the high station density on the Andean Chain. 3-D synthetic experiments are then designed to evaluate the network capacity to recover the event position using only P arrivals and the MAXI technique. Three synthetic earthquake location experiments are proposed: (1) noise-free and (2) noisy arrivals used in the 3-DVM, and (3) noise-free arrivals used in a 1-DVM. Synthetic results indicate that, under the best conditions (exact arrival data set and 3-DVM), the spatiotemporal configuration of the Ecuadorian network can accurately locate 70 per cent of events in the frontal part of the subduction zone (average azimuthal gap is 289° ± 44°). Noisy P arrivals (up to ± 0.3 s) can accurately located 50 per cent of earthquakes. Processing earthquake location within a 1-DVM almost never allows accurate hypocentre position for offshore earthquakes (15 per cent), which highlights the role of using a 3-DVM in subduction zone. For the application to real data, the seismicity distribution from the 3-D-MAXI catalogue is also compared to the determinations obtained in a 1-D-layered VM. In addition to good-quality location uncertainties, the clustering and the depth distribution confirm the 3-D-MAXI catalogue reliability. The pattern of the seismicity distribution (a 13 yr record during the inter-seismic period of the seismic cycle) is compared to the pattern of rupture zone and asperity of the Mw = 7.9 1942 and the Mw = 7.7 1958 events (the Mw = 8.8 1906 asperity patch is not defined). We observe that the nucleation of 1942, 1958 and 1906 events coincides with areas of positive Simple Bouguer anomalies and areas where marine terraces are still preserved on the coastal morphology. From north to south: (1) the 1958 rupture zone is almost aseismic and is attributed to a zone of high coupling; (2) south of the Galera alignment (perpendicular to the trench), the 1942 rupture zone presents moderate seismicity, deeper on the seismogenic interplate zone, and abutting on the Jama cluster (to the south). This cluster is facing the Cabo Pasado cap and positive Bouguer anomalies on the overriding margin. We suspect that this cluster reflects a zone of local asperity (partial coupling). South of the Jama cluster, the spherical aseismic zone in the Bahia area is interpreted as having a low seismic coupling (steady creep motion or slow slip events). We suspect that the site that generated the three M > 7 events (1896, 1956 and 1998) correspond to a small patch of strong coupling. To the south, in the Manta-Puerto Lopez zone, the seismicity is mainly organized in earthquake swarms (1998, 2002, 2005). Although slow slip events have been observed in the area (Vallée et al. submitted), we infer from the coastline shape, the marine terraces and the high positive Bouguer anomalies that the seismicity here might reveal a significant amount of seismic coupling.

Seismic velocity constraints on the material properties that control earthquake behavior at the Quebrada‐Discovery‐Gofar transform faults, East Pacific Rise

Mid‐ocean ridge transform faults (RTFs) vary strongly along strike in their ability to generate large earthquakes. This general observation suggests that local variations in material properties along RTFs exert a primary control on earthquake rupture dynamics. We explore these relationships by examining the seismic structure of two RTFs that have distinctly different seismic coupling. We determine the seismic velocity structure at the Gofar and Quebrada faults on the East Pacific Rise (EPR) using P wave traveltime tomography with data from two active‐source wide‐angle refraction lines crossing the faults. We image low‐velocity zones (LVZs) at both faults, where P wave velocities are reduced by as much as 0.5–1.0 km/s (∼10–20%) within a several kilometer wide region. At the Gofar fault, the LVZ extends through the entire crust, into the seismogenic zone. We rule out widespread serpentinization as an explanation for the low velocities, owing to the lack of a corresponding signal in the locally measured gravity field. The reduced velocities can be explained if the plate boundary region is composed of fault material with enhanced fluid‐filled porosity (1.5–8%). Local seismic observations indicate that the high‐porosity region lies within a ∼10 km long portion of the fault that fails in large swarms of microearthquakes and acts as a barrier to the propagation of large (M ∼ 6.0) earthquakes. Tomographic images of fault structure combined with observed earthquake behavior suggest that EPR transform segments capable of generating large earthquakes have relatively intact gabbro within the seismogenic zone, whereas segments that slip aseismically or via earthquake swarms are composed of highly fractured, ≥2 km wide damage zones that extend throughout the crust.

The seismic coupling of subduction zones revisited

The nature of seismic coupling for many of the world's subduction zones has been reevaluated. Geodetic estimates of seismic coupling obtained from GPS measurements of upper plate deformation during the interseismic period are summarized. We compared those with new estimates of seismic coupling obtained from seismological data. The results show that with a few notable exceptions the two methods agree to within about 10%. The seismological estimates have been greatly improved over those made 20–30 years ago because of an abundance of paleoseismological data that greatly extend the temporal record of great subduction earthquakes and by the occurrence, in the intervening years, of an unusual number of great and giant earthquakes that have filled in some of the most critical holes in the seismic record. The data also, again with a few notable exceptions, support the frictional instability theory of seismic coupling, and in particular, the test of that theory made by Scholz and Campos (1995). Overall, the results support their prediction that high coupling occurs for subduction zones subjected to high normal forces with a switch to low coupling occurring fairly abruptly as the normal force decreases below a critical value. There is also considerable variation of coupling within individual subduction zones. Earthquake asperities correlate with areas of high coupling and hence have a semblance of permanence, but the rupture zones and asperity distributions of great earthquakes may differ greatly between seismic cycles because of differences in the phase of seismic flux accumulation.

Frictional behavior of oceanic transform faults and its influence on earthquake characteristics

We use a three‐dimensional strike‐slip fault model in the framework of rate and state‐dependent friction to investigate earthquake behavior and scaling relations on oceanic transform faults (OTFs). Gabbro friction data under hydrothermal conditions are mapped onto OTFs using temperatures from (1) a half‐space cooling model, and (2) a thermal model that incorporates a visco‐plastic rheology, non‐Newtonian viscous flow and the effects of shear heating and hydrothermal circulation. Without introducing small‐scale frictional heterogeneities on the fault, our model predicts that an OTF segment can transition between seismic and aseismic slip over many earthquake cycles, consistent with the multimode hypothesis for OTF ruptures. The average seismic coupling coefficientχ is strongly dependent on the ratio of seismogenic zone width W to earthquake nucleation size h*; χ increases by four orders of magnitude as W/h* increases from ∼1 to 2. Specifically, the average χ = 0.15 ± 0.05 derived from global OTF earthquake catalogs can be reached at W/h* ≈ 1.2–1.7. Further, in all simulations the area of the largest earthquake rupture is less than the total seismogenic area and we predict a deficiency of large earthquakes on long transforms, which is also consistent with observations. To match these observations over this narrow range of W/h* requires an increase in the characteristic slip distance dcas the seismogenic zone becomes wider and normal stress is higher on long transforms. Earthquake magnitude and distribution on the Gofar and Romanche transforms are better predicted by simulations using the visco‐plastic model than the half‐space cooling model.

Reducing the effect of seismic noise in LIGO searches by targeted veto generation

One of the major obstacles to the detection and study of gravitational waves using ground-based laser interferometers is the effect of seismic noise on instrument sensitivity. Environmental disturbances cause motion of the interferometer optics, coupling as noise in the gravitational wave data output whose magnitude can be much greater than that of an astrophysical signal. We present an improved method of identifying times of high seismic noise coupling by tuning a gravitational-wave burst detection algorithm to the low-frequency signature of these events and testing for coincidence with a low-latency compact binary coalescence detection algorithm. This method has been proven highly effective in removing transients of seismic origin, with 60% of all compact binary coalescence candidate events correlated with seismic noise in just 6% of analysis time

The Sanriku‐Oki low‐seismicity region on the northern margin of the great 2011 Tohoku‐Oki earthquake rupture

We examine a region of the megathrust fault offshore of northeastern Honshu (38.75°–40.25°N, 141.5°–143.25°E) that we designate as the Sanriku‐Oki low‐seismicity region (SLSR). The SLSR, located near the northern termination of the 2011 Tohoku‐Oki (Mw 9.0) rupture, lacks historical great earthquake ruptures and has relatively low levels of moderate‐size (Mj ≥ 5.0) earthquakes, with subregions having many small events (Mj 2.5–5.0) in the Japan Meteorological Agency unified catalog. The SLSR is located downdip along the megathrust from the rupture zone of the great 1896 Sanriku tsunami earthquake and the great 1933 Sanriku outer trench slope normal‐faulting event; weak seismic coupling of the SLSR had been deduced based on the occurrences of those unusual events. Relatively low slip deficit on the SLSR megathrust was estimated based on GPS deformations prior to 2011 compared with adjacent areas with strong inferred locking to the south and north. The southern portion of the SLSR appears to have had, at most, modest levels (&lt;5 m) of coseismic slip during the 2011 event. Some thrust‐faulting aftershocks did occur in the SLSR, primarily at depths near 40 km where there had previously been small (Mj ∼5.0) repeating earthquakes (e.g., the Kamaishi repeater). An Mw ∼ 7.4 underthrusting aftershock occurred near the northeastern edge of the SLSR ∼22.5 min after the great 2011 event. Postseismic convergence along the megathrust is peaked in the SLSR. The collective observations indicate that the SLSR is primarily undergoing quasi‐static aseismic convergence, and the lack of regional strain accumulation likely delimited the northern extent of the great 2011 rupture as well as the downdip extent of the 1896 rupture. The triggering of the Mw 7.4 aftershock and heightened activity in the downdip repeater regions suggest that high postseismic strain rates drove the region to have ephemerally increased seismic failure, but it appears unlikely that a great earthquake will nucleate or rupture through this region. Similar properties may exist on the megathrust near the southern end of the 2011 rupture.

Dependence of earthquake stress drop on critical slip‐weakening distance

Numerical simulations are conducted to investigate the dependence of the static stress drop of earthquakes on the critical slip‐weakening distance. A fault model in a two‐dimensional elastic medium is used, in which a locked zone (asperity) is assumed to exist between creeping zones. The presence of such a zone is commonly assumed for plate boundaries that have relatively low seismic coupling. Shear stress is concentrated at the edges of the locked zone and seismic rupture occurs when the strain energy release by rupture extension becomes sufficiently large to overcome the fracture energy. The static stress drop is calculated by taking the ratio of the average seismic slip to the length of seismic slip zone for the simulated earthquakes in the model, using various values for the normal stress applied to the fault and the characteristic slip distance of a rate‐ and state‐dependent friction law. The simulation results indicate that the static stress drop is proportional to the square root of the product of the critical slip‐weakening distance and the normal stress, when the seismic rupture starts near an edge of the locked zone. In contrast, when the seismic rupture starts near the center of the locked zone, the static stress drop increases linearly with normal stress and is independent of the critical slip‐weakening distance. These simulation results may explain the relatively small depth dependence, as well as the large scatter, in observations of the static stress drop.

Earthquakes along Eltanin transform system, SE Pacific Ocean: fault segments characterized by strong and poor seismic coupling and implications for long-term earthquake prediction

SUMMARY Centroid moment tensor solutions are recomputed for 190 earthquakes from 1976 to 2010 along the Heezen, Tharp and Hollister transform faults of the Eltanin system using a 3-D seismic velocity model. The total length of the three en echelon faults is nearly 1000 km; each is characterized by fast long-term rates of displacement of about 80 mm yr–1. Strike-slip faulting with moment magnitudes Mw up to 6.4 characterizes most of these events. The few involving normal faulting are located up to 40 km on either side of the transforms and involve extension nearly normal to the transforms. This partitioning of slip likely results from changes during the last few million years in the Euler pole for relative motion between the Antarctic and Pacific plates. Some parts of the Heezen and Tharp transforms exhibit strong seismic coupling but others were aseismic at the resolution of our study, Mw > 5.0–5.5. Earthquakes were not found along nearby fast spreading ridges at that resolution. We calculate downdip widths of seismic coupling of about 5 km for four strongly coupled segments from observed moment rates and lengths along strike assuming earthquake activity accounts for the entire plate motion. Major differences in seismic coupling along strike are not in accord with common thermal models of plate cooling but instead are attributed to varying degrees of metamorphism, rock type and effective normal stress and possibly to the presence of short intratransform spreading centres. One 30–42-km-long segment of the Heezen transform that appears to be an isolated well-coupled asperity has ruptured in eight earthquakes of Mw 5.9–6.1 quasi-periodically with a coefficient of variation of 0.26 every 4.0 ± 1.0 yr. Other well-coupled fault segments, which were sites with earthquakes up to Mw 6.39 and fewer events since 1976, have average repeat times of about 7–24 yr. The fast rate of plate motion, maximum size of events and relatively short repeat times make these fault segments a good laboratory for research on quasi-periodic behaviour and earthquake prediction.