Search for: All records

Earthquakes are rupture-like processes that propagate along tectonic faults and cause seismic waves. The propagation speed and final area of the rupture, which determine an earthquake’s potential impact, are directly related to the nature and quantity of the energy dissipation involved in the rupture process. Here, we present the challenges associated with defining and measuring the energy dissipation in laboratory and natural earthquakes across many scales. We discuss the importance and implications of distinguishing between energy dissipation that occurs close to and far behind the rupture tip, and we identify open scientific questions related to a consistent modeling framework for earthquake physics that extends beyond classical Linear Elastic Fracture Mechanics.

Low-frequency earthquakes, atypical seismic events distinct from regular earthquakes, occur downdip of the seismogenic megathrust where an aseismic rheology dominates the subduction plate boundary. Well situated to provide clues on the slip regime of this unique faulting environment, their distinctive waveforms reflect either an unusual rupture process or unusually strong attenuation in their source zone. We take advantage of the unique geometry of seismicity in the Nankai Trough to isolate the spectral signature of low-frequency earthquakes after correcting for empirically derived attenuation. We observe that low-frequency earthquake spectra are consistent with the classical earthquake model, yet their rupture duration and stress drop are orders of magnitude different from ordinary earthquakes. We conclude their low-frequency nature primarily results from an atypical seismic rupture process rather than near-source attenuation.

It is well known that large earthquakes often exhibit significant rupture complexity such as well separated subevents. With improved recording and data processing techniques, small earthquakes have been found to exhibit rupture complexity as well. Studying these small earthquakes offers the opportunity to better understand the possible causes of rupture complexities. Specifically, if they are random or are related to fault properties. We examine microearthquakes (M < 3) in the Parkfield, California, area that are recorded by a high-resolution borehole network. We quantify earthquake complexity by the deviation of source time functions and source spectra from simple circular (omega-square) source models. We establish thresholds to declare complexity, and find that it can be detected in earthquakes larger than magnitude 2, with the best resolution above M2.5. Comparison between the two approaches reveals good agreement (>90 per cent), implying both methods are characterizing the same source complexity. For the two methods, 60–80 per cent (M 2.6–3) of the resolved events are complex depending on the method. The complex events we observe tend to cluster in areas of previously identified structural complexity; a larger fraction of the earthquakes exhibit complexity in the days following the Mw 6 2004 Parkfield earthquake. Ignoring the complexity of these small events can introduce artefacts or add uncertainty to stress drop measurements. Focusing only on simple events however could lead to systematic bias, scaling artefacts and the lack of measurements of stress in structurally complex regions.

Seismicity along transform faults provides important constraints for our understanding of the factors that control earthquake ruptures. Oceanic transform faults are particularly informative due to their relatively simple structure in comparison to their continental counterparts. The seismicity of several fast-moving transform faults has been investigated by local networks, but as of today there been few studies of transform faults in slow spreading ridges. Here, we present the first local seismicity catalogue based on event data recorded by a temporary broad-band network of 39 ocean–bottom seismometers located around the slow-moving Chain Transform Fault (CTF) along the Mid-Atlantic Ridge (MAR) from 2016 to 2017 March. We locate 972 events in the area by simultaneously inverting for a 1-D velocity model informed by the event P- and S-arrival times. We refine the depths and focal mechanisms of the larger events using deviatoric moment tensor inversion. Most of the earthquakes are located along the CTF (700) and Romanche transform fault (94) and the MAR (155); a smaller number (23) can be observed on the continuing fracture zones or in intraplate locations. The ridge events are characterized by normal faulting and most of the transform events are characterized by strike-slip faulting, but with several reverse mechanisms that are likely related to transpressional stresses in the region. CTF events range in magnitude from 1.1 to 5.6 with a magnitude of completeness around 2.3. Along the CTF we calculate a b-value of 0.81 ± 0.09. The event depths are mostly shallower than 15 km below sea level (523), but a small number of high-quality earthquakes (16) are located deeper, with some (8) located deeper than the brittle-ductile transition as predicted by the 600 °C-isotherm from a simple thermal model. The deeper events could be explained by the control of sea water infiltration on the brittle failure limit.

Earthquake stress drop is an important source parameter that directly links to strong ground motion and fundamental questions in earthquake physics. Stress drop estimations may contain significant uncertainties due to such factors as variations in material properties and data limitations, which limit the applications of stress drop interpretations. Using a high‐resolution borehole network, we estimate stress drop for 4551 (M0‐4) earthquakes on the San Andreas Fault at Parkfield, California, between 2001 and 2016 using spectral decomposition and an improved stacking method. To evaluate the influence of spatiotemporal variations of material properties on stress drop estimations, we apply different strategies to account for spatial variations of velocity and attenuation changes, and divide earthquakes into three separate time periods to correct temporal variations of attenuation. These results show that appropriate corrections can significantly reduce the scatter in stress drop estimates, and decrease apparent depth and magnitude dependence. We find that insufficient bandwidth can cause systematic underestimation of stress drop estimates and increased scatter. The stress drop measurements from the high‐frequency borehole recordings exhibit complex stable spatial patterns with no clear correlation with the nature of fault slip, or the slip distribution of the 2004 M6 earthquake. Temporal variations are significantly smaller, less well resolved and varying spatially. They do not affect the long‐term stress drop spatial variations, suggesting local material properties may control the spatial heterogeneity of stress drop.

We calculate rupture directivity and velocity for earthquakes in three well‐recorded repeating sequences (2001–2016) on the San Andreas Fault at Parkfield usingPwaves from borehole recordings and the empirical Green's function method. The individual events in each sequence all show the same directivity; the largest magnitude sequence (M ~ 2.7, 8 events) ruptures unilaterally NW (at ~0.8Vs), the second sequence (M ~ 2.3, 9 events) ruptures unilaterally SE, and the smallest magnitude sequence (M ~ 2, 11 events) is less well resolved. The highly repetitive rupture suggests that geometry or material properties might control nucleation of small locked patches. The source spectra of theM ~ 2.7 sequence exhibit no detectable temporal variation. The smallerMsequences both exhibit a decrease in high‐frequency energy following theM6 earthquake that recovers with time. This could indicate a decrease in stress drop, an increase in attenuation, or a combination of the two, followed by gradual healing.

We compare source parameter estimates for earthquakes in the 2011 Prague Mw 5.7, Oklahoma, sequence to investigate random uncertainty and systematic bias, and resolve reliable relative variations in stress drop. Source parameters provide insight into the earthquake rupture processes but large variations between studies occur. The Prague earthquake sequence is a prime example of this, with different studies reaching contrasting interpretations of the effects of injection on source parameters. We examine the Prague earthquake sequence using a single coherent catalog for all the events detected by the Oklahoma Geological Survey (OGS) and McMahon et al. (2017). We use three principal approaches to estimate stress drop in order to understand the biases of each: a spectral decomposition method based on stacking, individual event spectral modeling, and a spectral ratio method based on highly correlated events. We also compare our results with previous studies for the Prague sequences aftershocks, as well as past results for the Mw 4.8 foreshock and Mw 4.8 aftershock and Mw 5.7 mainshock. The absolute values of stress drop vary significantly between methods, but the relative patterns remain consistent, except when low quality or low bandwidth data are included. The consistent relative patterns reveal that the stress drops of aftershocks are dependent on the fault orientation and the proximity of the events to the mainshocks slip. These results indicate that fault structure as well as past events play an important role in stress drop patterns.

We image seismic attenuation near the Hikurangi trench offshore New Zealand, using ocean bottom and land‐based seismometers, revealing high attenuation above a recurring shallow slow‐slip event and within the subducting Hikurangi Plateau. The Hikurangi subduction margin east of the North Island, New Zealand is the site of frequent shallow slow slip events. Overpressured fluids are hypothesized to lead to slow slip at shallow depths close to the oceanic trench. Seismic attenuation, energy loss of seismic waves, can be used to detect high temperatures, melt, the presence of fluids, and fractures. We use local earthquake P‐ and S‐waves from 180 earthquakes to invert fort*, and subsequently invert for Qp and Qs, offshore the North Island directly above the area of slow slip. We image Qp and Qs to ∼25 km depth, increasing resolution of previously identified coastal lowQ(100–300), and finding a new region of even higher attenuation (Qp and Qs < 50–100) directly above the shallow slow slip event of 2014–2015, beneath the offshore seismic array. This highest attenuation is downdip of a subducting seamount, and is spatially correlated with a high seismic reflectivity zone and Vp/Vs > 1.85, all of which provide evidence for the presence of fluids. The Qp and Qs is low at the trench (<50–100) and in the subducting plate (100–200), suggesting that seismic wave scattering due to faults, fractures, and the inherent heterogeneous composition of the Hikurangi Plateau, a large igneous province, plays a role in seismic attenuation.