Stress drops and radiated energies of aftershocks of the 1994 Northridge, California, earthquake

We study stress levels and radiated energy to infer the rupture characteristics and scaling relationships of aftershocks and other southern California earthquakes. We use empirical Green functions to obtain source time functions for 47 of the larger (M ≥ 4.0) aftershocks of the 1994 Northridge, California earthquake (M6.7). We estimate static and dynamic stress drops from the source time functions and compare them to well-calibrated estimates of the radiated energy. Our measurements of radiated energy are relatively low compared to the static stress drops, indicating that the static and dynamic stress drops are of similar magnitude. This is confirmed by our direct estimates of the dynamic stress drops. Combining our results for the Northridge aftershocks with data from other southern California earthquakes appears to show an increase in the ratio of radiated energy to moment, with increasing moment. There is no corresponding increase in the static stress drop. This systematic change in earthquake scaling from smaller to larger (M3 to M7) earthquakes suggests differences in rupture properties that may be attributed to differences of dynamic friction or stress levels on the faults

Earthquake rupture below the brittle-ductile transition in continental lithospheric mantle

Earthquakes deep in the continental lithosphere are rare and hard to interpret in our current understanding of temperature control on brittle failure. The recent lithospheric mantle earthquake with a moment magnitude of 4.8 at a depth of ~75 km in the Wyoming Craton was exceptionally well recorded and thus enabled us to probe the cause of these unusual earthquakes. On the basis of complete earthquake energy balance estimates using broadband waveforms and temperature estimates using surface heat flow and shear wave velocities, we argue that this earthquake occurred in response to ductile deformation at temperatures above 750°C. The high stress drop, low rupture velocity, and low radiation efficiency are all consistent with a dissipative mechanism. Our results imply that earthquake nucleation in the lithospheric mantle is not exclusively limited to the brittle regime; weakening mechanisms in the ductile regime can allow earthquakes to initiate and propagate. This finding has significant implications for understanding deep earthquake rupture mechanics and rheology of the continental lithosphere

SMatStack to enhance noisy teleseismic seismic phases: validation and application to resolving depths of oceanic transform earthquakes

The depths of most moderate‐sized oceanic earthquakes are poorly constrained because of the lack of local recording stations and noisy teleseismic recordings. This hampers our understanding of slip behaviors along oceanic faults and the mechanical properties of the oceanic lithosphere. In this study, we develop a new method to enhance the weak body‐wave signals, particularly the depth phases, associated with earthquakes on oceanic transform faults using large‐aperture arrays in the teleseismic range. We simulate the enhanced teleseismic signals to refine the centroid depths of moderate‐sized earthquakes. We validate the new approach using synthetic waveforms and show it outperforms conventional methods when dealing with noisy signals. We obtain the depth estimates for three moderate‐sized earthquakes on the Chain transform fault in the equatorial Atlantic Ocean and find two of them are consistent with a local catalog derived using oceanic bottom seismometers. Application of the new method to the past decades of teleseismic recordings of moderate‐sized earthquakes on the large and slow‐slipping transform faults will provide significantly improved constraints on the width of the seismogenic zone, thus advancing our understanding of the rheology of oceanic lithosphere and earthquake processes in oceanic settings and, by comparison, their more dangerous continental counterparts. The new method is not limited to oceanic transform earthquakes, but can be easily adapted to other seismological studies in which noisy but coherent signals could be revisited for better usage.https://doi.org/10.1029/2023GC011109Published versio

Investigating spectral estimates of stress drop for small to moderate earthquakes with heterogeneous slip distribution: examples from the 2016–2017 Amatrice earthquake sequence

Estimates of spectral stress drop are fundamental to understanding the factors controlling earthquake rupture and high frequency ground motion, but are known to include large, poorly‐constrained uncertainties. We use earthquakes from the 2016–2017 sequence in the Italian Appenines (largest event at Norcia, M_w 6.3) to investigate these uncertainties and their causes. The similarly‐sized events near Amatrice (M_w 6.0) and Visso (M_w 5.9) enable better constrained relative analysis. We calculate S wave source spectra, corner frequencies, and spectral stress drop for 30 of the larger events. We compare both empirical and modeling approaches to isolate the source spectra and calculate source parameters; we also compare our results with those from published studies. Both random and systematic inter‐study variations are larger than the standard errors reported by any individual study. The reported magnitude dependence of stress drop varies between studies, being largest for generalized inversions and smallest for more individual event based approaches. The relative spectral estimates of inter‐event stress drop are more consistent; all approaches estimated higher stress drop in the Amatrice earthquake than the similar‐sized Visso earthquake. In contrast, finite fault inversions of these two earthquakes found that the Visso earthquake had the larger region of concentrated, higher slip, whereas the Amatrice earthquake had multiple, lower slip, subevents. The Amatrice spectra contain more high frequency energy than those of the Visso earthquake. This comparison suggests that consistent measurement of a higher spectral stress drop indicates greater high‐frequency ground motion but may correspond to greater rupture complexity rather than higher stress drop.EAR-2043281 - National Science Foundationhttps://doi.org/10.1029/2022JB025022Published versio

Spatially consistent small-scale stress heterogeneity revealed by the 2008 Mogul, Nevada, earthquakes

We compute and analyze stress drops for 4175 earthquakes (M_L 0–5) in the 2008 Mogul, Nevada, swarm–mainshock sequence using a spectral decomposition approach that uses depth-dependent path corrections. We find that the highest stress-drop foreshocks occur within the fault zone of the M_w 4.9 mainshock, nucleating at the edges of seismicity voids and concentrating near complexities in the fault geometry, confirming and extending inferences from prior work based on empirical Green’s functions for ∼150 of the larger Mogul earthquakes. The region of the highest stress-drop foreshocks is not reruptured by aftershocks, whereas low-stress-drop areas are consistently low during both the foreshock and aftershock periods, implying that stress drop depends on inherent individual fault properties rather than timing within the sequence. These results have implications for swarm evolution and fault activation within complex 3D structures.EAR-2043281 - National Science Foundationhttps://doi.org/10.1785/0320230026Published versio

Local seismicity around the Chain Transform Fault at the Mid-Atlantic Ridge from OBS observations

Summary 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 broadband network of 39 ocean bottom seismometers located around the slow-moving Chain Transform Fault (CTF) along the Mid-Atlantic Ridge (MAR) from March 2016 to March 2017. 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 characterised by normal faulting and most of the transform events are characterised 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 seawater infiltration on the brittle failure limit

Planetary Interiors

This report identifies two main themes to guide planetary science in the next two decades: understanding planetary origins, and understanding the constitution and fundamental processes of the planets themselves. Within the latter theme, four specific goals related to interior measurements addressing the theme. These are: (1) Understanding the internal structure and dynamics of at least one solid body, other than the Earth or Moon, that is actively convecting, (2) Determine the characteristics of the magnetic fields of Mercury and the outer planets to provide insight into the generation of planetary magnetic fields, (3) Specify the nature and sources of stress that are responsible for the global tectonics of Mars, Venus, and several icy satellites of the outer planets, and (4) Advance significantly our understanding of crust-mantle structure for all the solid planets. These goals can be addressed almost exclusively by measurements made on the surfaces of planetary bodies

Quantifying rupture characteristics of microearthquakes in the parkfield area using a high-resolution borehole network

[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 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 M_w 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.]EAR-2043281 - National Science Foundationhttps://doi.org/10.1093/gji/ggad023.Published versio

Quantifying Site Effects and Their Influence on Earthquake Source Parameter Estimations Using a Dense Array in Oklahoma

We investigate the effects of site response on source parameter estimates using earthquakes recorded by the LArge-n Seismic Survey in Oklahoma (LASSO). While it is well known that near-surface unconsolidated sediments can cause an apparent breakdown of earthquake self-similarity, the influence of laterally varying site conditions remains unclear. We analyze site conditions across the 1825-station array on a river plain within an area of 40 km by 23 km using vertical ground motions from 14 regional earthquakes. While the source radiation pattern controls P-wave ground motions below 8 Hz, the surface geology correlates with P-wave ground motions above 8 Hz and S-wave ground motions at 2–21 Hz. Stations installed in alluvial sediments have vertical ground motions that can exceed three times the array median. We use the variation of ground motion of regional earthquakes across the array as a proxy for site effects. The corner frequencies and stress drops of local earthquakes (ML = 0.01–3) estimated using a standard single-spectra approach show negative correlations with the site-effect proxy, while the seismic moments show positive correlations. In contrast, the spectral-ratio approach effectively shows no correlation. The overall bias is small as expected for this relatively homogeneous structure; accurate estimation of site-related biases requires at least 30 stations. Correcting for site-related biases reduces the standard deviations of the source parameters by less than 13% of the total variations. Remaining variations are partially associated with source directivity and model misfits— as small earthquakes can have complex ruptures

Energy dissipation in earthquakes

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.Comment: 8 pages, 2 figure