Seismic Velocity Structure Research Articles

Substantial Global Radial Variations of Basalt Content Near the 660‐km Discontinuity

AbstractMid‐ocean ridges generate basalt and harzburgite, which are introduced into the mantle through subduction as a mechanical mixture, contributing to both lateral and radial compositional heterogeneity. The possible accumulation of basalt in the mantle transition zone has been examined, but details of the mantle composition below the 660‐km discontinuity (hereafter d660) remain poorly constrained. In this study, we utilize the subtle waveform details of S660S, the underside shear‐wave reflection off the d660, to interpret the seismic velocity, density, and compositional structure near, and particularly below, the d660. We identify a significant difference in S660S waveform shape in subduction zones compared to other regions. The inversion results reveal globally enriched basalt at the d660, with a notably higher content in subduction zones, consistent with the smaller impedance jump and S660S peak amplitude. The basalt fraction decreases significantly to less than 10% near 800‐km depth, forming a global harzburgite‐enriched layer and resulting in a steep seismic velocity gradient just below the d660, in agreement with 1D global reference models. The striking compositional radial variations near the d660 verify geodynamic predictions and challenge the applicability of homogeneous radial compositional models in the mantle. These variations may also affect the viscosity profile and, consequently, the dynamics at the boundary between the upper and lower mantle.

Development of a 3D Geologic Model Used in the Seismic Hazard and Liquefaction Hazard Analysis of Madison County, Tennessee

ABSTRACT Madison County of western Tennessee is approximately 100 km southeast of the New Madrid seismic zone and, thus, is subject to potentially dangerous seismic shaking and liquefaction events. Holocene floodplain alluvium and associated Pliocene/Pleistocene terraces of the Forked Deer and Hatchie Rivers cover a large part of Madison County. Underlying these unconsolidated sediments are approximately ∼650 m of poorly consolidated Paleogene and Cretaceous terrestrial and marine sediments, which are underlain by Paleozoic limestone bedrock. Uniform regional geology used in the 2014 United States Geological Survey seismic hazard maps assumes a homogenous seismic velocity structure above bedrock. Three-dimensional geologic mapping of the ∼700 m of Quaternary, Neogene, Paleogene, and Cretaceous sediments conducted in this study provides a more detailed velocity structure, resulting in a more accurate and detailed map of the expected ground motions in Madison County during future earthquakes. In addition, the surface geologic map reveals regions that are susceptible to liquefaction. This paper discusses the methods used to construct the surface geology map, subsurface structure contour maps, isopach maps, and 3D geologic model used in the seismic hazard and liquefaction hazard analyses of Madison County, Tennessee, and summaries of the results of these analyses are presented.

Impact of the Offshore Seismograph Network and 3‐D Seismic Velocity Structure Model on Centroid Moment Tensor Analysis for Offshore Earthquakes: Application to the Japan Trench Subduction Zone

Impact of the Offshore Seismograph Network and 3‐D Seismic Velocity Structure Model on Centroid Moment Tensor Analysis for Offshore Earthquakes: Application to the Japan Trench Subduction Zone

Vertical high-velocity structures and seismic activity in western Shandong Rise, China: Case study inspired by double-difference seismic tomography

Abstract In this article, we collected the seismic phase arrival data (115°–120°E, 34°–39°N) of 2,833 local natural earthquakes above magnitude 1.0 recorded by 128 seismograph stations provided by the China Earthquake Networks Center covering the period from January 2008 to October 2023. We extracted the first arrival P- and S-wave arrival time data and obtained 26,351 P- and 26,349 S-wave absolute arrival times and 99,627 P-differential and 99,625 S-differential arrival times. Then, we determined 3-D P and S wave velocity structures under the Luxi region by using the double-difference tomography method. The results show lateral heterogeneities under the Luxi region, and the distribution of surface basins and uplift zones is associated with seismic velocities within the crust. The vertical section shows the widespread presence of vertical high-velocity bodies within the crust of the Luxi region, indicating seismological evidence of mantle upwelling in the region. Meanwhile, our imaging results also clearly indicate the presence of a significant low-velocity anomaly at depths of 13–17 km, which corresponds to the presence of a high-conductivity layer at that depth obtained from electromagnetic bathymetry. The transition zone between high and low velocities exists at a depth of 7–13 km, which plays an important role in the transfer of lithospheric stresses from the deep to the shallow part of the lithosphere. Combined with the results of the previous deep seismic wave velocity structure, petrology, and geological investigations, the deep dynamical background of the formation of slip tectonics and seismic mechanism in the Luxi area can be attributed to the go-slip activity of the Tan-Lu fault Zone, the subduction of the Pacific Plate, and the magmatism generated by mantle upwelling.

Probing a planet from the subsurface to the atmosphere with infrasound data

The analysis of infrasound waves has a significant potential for retrieving a range of geophysical parameters across various scales, such as atmospheric structures, characteristics of surface and buried sources, and seismic velocity structures. This potential was recently showcased in infrasound studies that illustrated the ability to retrieve Earth's shallow velocities and Mars' near-surface atmospheric winds from remote acoustic observations. Consequently, infrasound is becoming an essential component of planetary missions, providing insights into the interiors and atmospheres of other terrestrial planets in the solar system. However, utilizing infrasound data requires efficient forward wave simulation techniques and an accurate description of waveform sensitivity to source and medium parameters. Even under ideal conditions, many inverse problems associated with infrasound-based probing are inherently ill-posed, necessitating regularization or other approaches to ensure reliable results. In this contribution, we highlight recent seismo-acoustic research findings, addressing atmospheric probing at both smaller and global scales, as well as innovative methods to accelerate infrasonic wave propagation modeling.

A 3D Seismotectonic Model and the Spatiotemporal Relationship of Two Historical Large Earthquakes in the Linfen Basin, North China

The Shanxi Graben is a transitional zone between the Ordos Block and North China Plain with complex structures and frequent earthquakes. Six earthquakes with M ≥ 7.0 have been recorded in the area, including the 1303 Hongtong M 8 and 1695 Linfen M 7.8 earthquakes in the Linfen Basin. Research on these two large earthquakes, closely related in time and space, is lacking. Our objective was to use deep seismic reflection profiles and 3D velocity structure data from previous research, along with seismological observation results, to interpret the geological structure near the source of the two earthquakes. A 3D geometric model of the seismogenic fault was constructed, and the relationships among the deep and shallow structures, deep seismogenic environment, and two large earthquakes were explored. Differences in seismogenic environment between the southern and northern Linfen Basin were identified. The distribution of small earthquakes in the southern Linfen Basin was scattered, and the overall distribution was at depths <25 km. The small earthquakes in the northern part of the basin were dense and concentrated at depths of 25–35 km. Low-velocity layers at an approximate depth of 15–20 km in the southern basin led to differences in seismogenesis between the two regions. Based on the area of the 3D geometric model of the Huoshan Fault, the maximum magnitude of an earthquake caused by fault rupture is Mw 7.7, so the magnitude of the 1303 Hongtong earthquake might be overestimated. Numerical simulation results of Coulomb stress showed that the 1303 Hongtong earthquake had a stress-loading effect on the 1695 Linfen earthquake. The change in Coulomb rupture stress was 1.008–2.543 bar, which is higher than the generally considered earthquake trigger threshold (0.1 bar). We created a new 3D source model of large earthquakes in the Linfen Basin, Shanxi Province, providing a reference and typical cases for risk assessment of large earthquakes in different regions of the Shanxi Graben.

The effect of antigorite dehydration on velocity structure and water migration in subduction zones

The water migration in subduction zones, primarily driven by the phase transition in hydrous minerals, can give rise to hydrated regions with reduced velocity. A fundamental element in comprehending and deciphering these low-velocity zones revolves around acquiring insights into the stability and elasticity of relevant hydrous minerals. As one of the main water carriers in shallow areas, antigorite can dehydrate to form talc, forsterite, and fluid (talc–bearing peridotites) in deep areas of subduction zones, and then the talc thus serves as one of the minerals that can bring water to the deep Earth. Here, the elasticity of talc up to 24 GPa and forsterite up to 12 GPa are calculated by using the first principles method. The result supposes that the talc structure transforming from talc I to talc II is at a pressure between 6 GPa and 8 GPa, impacting the trend of elastic wave velocity in response to pressure. Furthermore, the elastic wave velocity of forsterite can be significantly affected by iron concentration. Meanwhile, a variation velocity model with antigorite consumption and talc content is set up for talc-bearing serpentinized peridotite based on the elastic properties of talc and forsterite in this study, and antigorite in Wang et al. (2022). The results of our model demonstrate a decrease in the low-velocity anomaly in subduction zones, particularly in deep regions or areas with higher initial serpentinization degrees. The results also suggest that the mode of antigorite dehydration can diminish the estimation of water content transported to depths of subduction zones, such as the Mariana Trench and Northern Japan subduction zones. The mode of antigorite dehydration thus provides a useful tool for constraining the composition, seismic velocity structure, and water migration in subduction zones.

Extracting Multimodal Surface-Wave Dispersion Curves from Ambient Seismic Noise Using High-Resolution Linear Radon Transform

Abstract In the past two decades or so, ambient noise tomography (ANT) has emerged as an established method for imaging subsurface seismic velocity structures. One of the key steps in ANT is to extract surface-wave dispersion curves. The predominant approach for subsurface shear-wave velocity structure inversion involves utilizing fundamental-mode surface waves in ANT. Nevertheless, a notable challenge encountered is the issue of nonuniqueness when employing the dispersion information of fundamental-mode surface waves to invert for shear-wave velocity models. The inclusion of higher-mode dispersion curves in the inversion offers several benefits, including the reduction of nonuniqueness, enhancement of inversion stability, and decreased dependence on the initial model. In this study, we illustrate the applicability of the high-resolution linear radon transform method (HRLRT) for extracting multimodal surface-wave dispersion energy from ambient seismic noise data. We apply HRLRT to both the synthetic noise data and real data recorded by USArray. Our results of applications show that the HRLRT method can extract multimodal surface-wave dispersion information. Compared with established methods such as the frequency–Bessel transform and multicomponent frequency–Bessel transform, the HRLRT exhibits an advantage in suppressing “crossed” artifacts and the second-/third-type artifacts caused by sparse spatial sampling, and the resulting dispersion energy from HRLRT has narrower peaks, meaning high resolution of dispersion curves based on the HRLRT method.

On the Footsteps of Active Faults from the Saronic Gulf to the Eastern Corinth Gulf: Application of Tomographic Inversion Using Recent Seismic Activity

This study examines the body-wave velocity structure of Attica, Greece. The region is located between two major rifts, the Gulf of Corinth and the Euboekos Gulf, and has experienced significant earthquakes throughout history. The distribution of seismic activity in the area necessitates a thorough investigation of geophysical properties, such as seismic velocities, to reveal the extent of significant fault zones or the presence of potential hidden faults. This case study utilized over 3000 revised events to conduct a local earthquake tomography (LET). P- and S-wave travel-time data were analyzed to explore small- to medium-scale (~10 km) anomalies that could be linked to local neotectonic structures. The study presents a detailed 3-D seismic velocity structure for Attica and its adjacent regions. The results of the study revealed strong lateral body-wave velocity anomalies in the upper crust were related to activated faults and that a significant portion of the observed seismicity is concentrated near the sites of the 1999 and 2019 events.

Global variability of the composition and temperature at the 410-km discontinuity from receiver function analysis of dense arrays

Seismic boundaries caused by phase transitions between olivine polymorphs in Earth's mantle provide thermal and compositional markers that inform mantle dynamics. Seismic studies of the mantle transition zone often use either global averaging with sparse arrays or regional sampling from a single dense array. The intermediate approach of this study utilizes many densely spaced seismic arrays distributed around the globe. We systematically compute teleseismic P-to-S receiver functions for each seismic array and invert for the 1-D seismic velocity structure of the mantle transition zone beneath each array to facilitate a comparison between densely sampled regions. We stack 3,600 receiver functions on average at 67 arrays in total. The stack is used in a probabilistic inversion to estimate the mantle transition zone interface depths and velocities beneath each array. We focus on the 410-km discontinuity (410) because it is a prominent seismic interface that is clearly linked to a single mineral phase transition between olivine and wadsleyite. The depths and velocity contrasts of the 410 are mapped to temperatures and compositions using mineral physics constraints. The depth of the 410 ranges from ∼405–440 km, which is consistent with a ∼360 K temperature range in a dry mantle and a ∼260 K temperature range in a wet mantle (2 wt. % water). The Vs contrast across the 410 ranges from ∼2.5–8 %, which is consistent with ∼20–70 vol. % olivine composition in a dry mantle and ∼25–80 vol. % in a wet mantle. The bulk composition of the upper mantle near the 410-km discontinuity is typically considered to be well-mixed because there is no thermodynamic impediment to convection at the olivine to wadsleyite phase transition. However, the wide range of inferred olivine content from our study suggests that there are large lateral variations in the bulk composition of the upper mantle near the 410-km discontinuity.

Seismic Imaging of Shallow Depths using Ambient Noise Tomography in the Hyderabad Region, Eastern Dharwar Craton, India

ABSTRACT Hyderabad and its surrounding region are situated within Eastern Dharwar Craton, primarily composed of Archean granite. We use ambient noise tomography to obtain the 3-D shallow subsurface shear seismic velocity structure of 150 x 200 km2 area around Hyderabad. We used continuous waveforms that were recorded by a digital seismic network with 10 three-component broadband seismographs from November 2020 to December 2021. Using the vertical component of waveforms, the Empirical Green’s functions (EGFs) have been estimated, and the group velocity dispersion curves in the time period range of 1.0–10 s have been extracted from these. We inverted these dispersion curves to estimate 3-D variation of shear wave velocity by using direct surface wave tomographic method. The results suggest that the subsurface velocity structure is heterogeneous. The shear wave velocity varies from 3.4 to 3.7 km/s in the shallow subsurface. The limited seismicity in this area is mostly confined to the boundary between high and low-velocity changes.

Exploring a planet with infrasound: Challenges in probing the subsurface and the atmosphere

Infrasound waves have shown great potential to retrieve a range of geophysical parameters across scales, such as atmospheric structures, surface and buried source characteristics, and seismic velocity structures. In particular, infrasound generated by seismic waves coupling to the atmosphere can provide unique insight into a planet's interior. However, utilizing infrasound data require efficient forward wave simulation techniques and an accurate description of waveform sensitivity to source and medium parameters. Even under idealized conditions, many of the inverse problems associated with infrasound-based probing are inherently ill-posed, requiring regularization or other approaches to yield reliable output. In this contribution, we highlight recent research results from our group and collaborators, tackling atmospheric probing at smaller and global scales, as well as innovative approaches to speed up infrasonic wave propagation modeling. We also call attention to the objectives and early results from a newly launched project focusing on the feasibility of exploring the interior of Venus using infrasound recorded from balloon platforms.

Upper crustal seismic velocity structure of the Hayward fault zone, San Francisco Bay, California, USA: Results from the 2016 East Bay Seismic Experiment (EBSI-16)

We developed Vp, Vs, Vp/Vs ratio, and Poisson’s ratio models of the uppermost crust (<4 km depth) from the eastern San Francisco (SF) Bay (California, USA) to near the Calaveras fault along a 15-km-long, linear profile. Upper crustal velocities are highly variable beneath, west, and well east of the Hayward fault. We observe eight notable features, from west to east: (1) Near San Francisco Bay, there is an ∼2-km-wide structure with high Vp/Vs ratios (up to 5) and Poisson’s ratios (up to 0.48) extending from the surface to the base of our model, which we suggest the structure is a near-vertical fault that lies along a straight-line projection between the Silver Creek fault to the south and the Point Richmond fault to the north. The structure may be part of an ∼90-km-long fault along the eastern SF Bay. (2) The western East Bay Plain, the lower lying area between the bay and the hills, includes up to 800 m of low-velocity sediments (Vp ∼1600−3000 m/s, Vs ∼500 m/s to ∼1000 m/s), underlain by higher velocity basement rocks (Vp ∼3000−5800 m/s; Vs ∼1000−1500 m/s). (3) Between ∼1 km and 3 km east of the Bay shoreline, sediments thin in a series of steps (likely faults) toward the Hayward fault. (4) Between ∼3 km west and ∼1 km east of the Hayward fault (at the East Chabot fault) at depths greater than 1 km, basement Vp (up to 6000 m/s) and Vs (up to 2800 m/s) are high, and Vp/Vs ratios (<2) and Poisson’s ratios (<0.3) are low, suggesting crystalline rocks. Furthermore, a near-vertical zone of low Vp/Vs ratios and Poisson’s ratios is between near-surface traces of the Hayward and East Chabot faults, likely corresponding to the San Leandro Gabbro of Ponce et al. (2003). (5) Eastward of the East Chabot fault in the upper 1.5 km, basement Vp (∼3000 m/s to ∼4200 m/s) and Vs (∼1200−2000 m/s) are lower than those west of the fault. (6) In the eastern Hayward/Oakland Hills, there are zones of laterally varying, high- and low-velocity (Vp ∼2500−3000 m/s) Jurassic−Cretaceous and Tertiary sediments in the shallow subsurface that likely extend much deeper than imaged. (7) Seismic energy that propagates westward from sources east of the Hayward fault (HF) appear weaker than energy that propagates eastward from sources west of the HF, suggesting that the HF acts as a partial barrier to shallow seismic energy propagation into the more populated eastern SF Bay area. (8) Unlike many fault zones, it appears that the active trace of the Hayward fault (in our study area) is not cored by a prominent, low-velocity zone relative to rocks to the east and west of the active trace. However, the active trace does mark a prominent change from relatively higher velocities to the west and lower velocities to the east.

Multi-scale anisotropy in NE China: Evidence for localized mantle upwelling

It is commonly proposed that the subduction of the Pacific plate has been responsible for widespread Holocene intraplate volcanism across NE China and the Korean Peninsula. Yet, how this process drives volcanism and even if it plays a critical role remains a topic of vigorous debate. In this study, using seismic data from four networks across NE China and northern Democratic People's Republic of Korea (DPRK), we analyze shear wave splitting in converted P to S-waves at the Moho (Pms), S-waves from the subducted slab interface (local S), and SKS phases. The Pms phases show a relatively weak crustal anisotropy (<0.25 s), with fast polarization directions aligned sub-parallel to major tectonic features. For the local S and SKS phases, fast polarization directions show significant lateral variations. We further perform a quantitative inversion to show that the depth of the anisotropy is ∼150 km, thus driven by flow within the asthenosphere associated with Pacific subduction. However, the presence of many null SKS splitting phases, together with scattered local S anisotropy across a wide range of incidence angles suggests a localized region of vertical flow directly beneath Changbaishan volcano. Such patterns correspond well to regional upper-mantle seismic velocity structure, and suggest that a localized upwelling with a relatively deep origin drives volcanism in the Changbaishan region. Furthermore, we infer that this mantle upwelling is deflected to the SW beneath Changbaishan and spreads asymmetrically at the base of the lithosphere, possibly because of the long history of volcanism in the region.

Structural Heterogeneity Controlled Rupture Process of the 2021 Mw 7.1 Fukushima, Japan, Earthquake Revealed by Joint Inversion of Seismic and Geodetic Data

Abstract We determined the rupture model of the 2021 Mw 7.1 Fukushima earthquake near northeastern Japan in this study and adopted this model to investigate the cause of this earthquake and its aftershocks. The rupture model was obtained through joint inversion of teleseismic, strong-motion and geodetic data. It is shown that the slips were predominantly distributed on the southwest side of the earthquake epicenter, indicating a unilateral rupture event. We observed that the seismic moment was released in three time periods, producing four slip patches on the fault plane. Through comparison, we demonstrated that our joint inversion model was more reliable in describing the rupture process of the Fukushima earthquake than the automatic inversion models determined using only strong-motion data. By jointly analyzing the slip distribution and seismic velocity structure, we found a good correlation between the slip patches and VP/VS anomalies, suggesting that structural heterogeneities along the fault zone played a critical role in controlling the rupture process of the Fukushima earthquake. In addition, most aftershocks were located in the region characterized by small slips and high VP/VS, and we demonstrated that they were caused by stress changes due to the presence of fluids and the rupture of the mainshock.

Seismic velocity structure and seismicity at the Kairei hydrothermal vent field near the Rodriguez Triple Junction in the Indian Ocean

The Kairei hydrothermal vent field (KHF) is on the northeastern side of the first segment of the Central Indian Ridge and discharges H2-rich hydrothermal fluids. To identify the source of these fluids based on the crustal and upper mantle structures and tectonic stress fields, we conducted a seismic refraction survey and observed natural seismicity in and around the KHF. We estimated 3D P-wave and S-wave velocity (VP and VS, respectively) structure models from airgun and natural seismic sources and determined 1861 hypocenters during a period of about two months. Beneath Yokoniwa Rise, VP in a convex zone exceeded 6.0 km/s at ~2 km depth below the seafloor and exceeded 7.0 km/s at 3–4 km depth; we thus interpret the rise to be a non-transform offset massif formed by the uplift of deep rocks. We detected a seismic cluster in a normal fault zone 3–5 km northwest of the KHF. Many events in this cluster occurred at 2–5 km depth below the seafloor in an area of high VP/VS ratios (~2.0) that are consistent with serpentinized peridotite, suggesting that seawater seeping downward through the fault zone has serpentinized peridotite at depth, thus becoming H2-rich hydrothermal fluid. This high-VP/VS area extends laterally at depth, and other such areas are observed at <1 km below the seafloor near the KHF. This structural configuration suggests that the high-VP/VS areas might provide conduits for the hydrothermal fluids that issue from the vents of the KHF. At the first segment of the ridge, we detected distinct shallow and deep seismic clusters. Seismicity in the shallow cluster was intermittent whereas deeper seismicity was nearly continuous, possibly indicating distinct magmatic activities at each depth.

High-resolution seismic velocity structure using induced inter-event interferometry at the geothermal sites, Geysers, Eastern California

The possibility of retrieving Rayleigh wave empirical Green's functions (EGFs) using virtual seismometers provides an opportunity to calculate high-resolution velocity models without a set of dense stations and/or expensive seismic imaging/monitoring. The high-resolution tomographic map can give an overview of the fracturing within the saturated rock. Theoretically, the cross-correlations of earthquakes contain significant information about the EGFs between event-pairs. In this work, we developed an application of virtual seismometers on a very small scale in superficial layers where anthropogenic seismicity occurred. To test this method's capabilities, we applied it to investigate the NW-part of The Geysers (TG) geothermal field. We used 1276 micro-earthquakes, occurred in a cuboid with the horizontal side of ∼1 × ∼2 km and the vertical edge at a depth from 0.9 to 2.2 km, to obtain the high-resolution (100 × 100 m) velocity structure. A group of strict criteria was imposed separately for single event waveform and event pairs to stabilize the tomographic results. After retrieving the inter-event EGF signal and calculating the Rayleigh wave group velocity dispersion curves, we obtained the group velocity maps for each horizontal layer. The tomographic maps as a function of depth indicated a low-velocity anomaly related to existing inactive faults, topography variation of the top of the steam zone, and low-velocity anomalies possibly stimulated by fluid injection. The topography of the top of the steam varies from ∼1.0 km to 1.3 km of maps, and also anomalies connected with fluid migration appear at the depth of injection wells (1.5–1.9 km).

Solidified magma reservoir derived from active source seismic experiments in the Aira caldera, southern Kyushu, Japan

The Aira caldera, located in southern Kyushu, Japan, originally formed 100 ka, and its current shape reflects the more recent 30 ka caldera-forming eruptions (hereafter, called the AT eruptions). This study aimed to delineate the detailed two-dimensional (2D) seismic velocity structure of the Aira caldera down to approximately 15 km, by means of the travel-time tomography analysis of the seismic profile across the caldera acquired in 2017 and 2018. A substantial structural difference in thickness in the subsurface low-velocity areas in the Aira caldera between the eastern and western sides, suggest that the Aira caldera comprises at least two calderas, identified as the AT and Wakamiko calderas. The most interesting feature of the caldera structure is the existence of a substantial high-velocity zone (HVZ) with a velocity of more than 6.8 km/s at depths of about 6–11 km beneath the central area of the AT caldera. Because no high ratio of P- to S-wave velocity zones in the depth range were detected from the previous three-dimensional velocity model beneath the AT caldera region, we infer that the HVZ is not an active magma reservoir but comprises a solidified and cool remnant. In addition, a poorly resolved low-velocity zone around 15 km in depth suggests the existence of a deep active magma reservoir. By superimposing the distribution of the known pressure sources derived from the observed ground inflation and the volcanic earthquake distribution onto the 2D velocity model, the magma transportation path in the crust was imaged. This image suggested that the HVZ plays an important role in magma transportation in the upper crust. Moreover, we estimated that the AT magma reservoir in the 30 ka Aira caldera-forming eruptions has the total volume of 490 km3 DRE and is distributed in a depth range of 4–11 km.Graphical

Broad fault zones enable deep fluid transport and limit earthquake magnitudes

Constraining the controlling factors of fault rupture is fundamentally important. Fluids influence earthquake locations and magnitudes, although the exact pathways through the lithosphere are not well-known. Ocean transform faults are ideal for studying faults and fluid pathways given their relative simplicity. We analyse seismicity recorded by the Passive Imaging of the Lithosphere-Asthenosphere Boundary (PI-LAB) experiment, centred around the Chain Fracture Zone. We find earthquakes beneath morphological transpressional features occur deeper than the brittle-ductile transition predicted by simple thermal models, but elsewhere occur shallower. These features are characterised by multiple parallel fault segments and step overs, higher proportions of smaller events, gaps in large historical earthquakes, and seismic velocity structures consistent with hydrothermal alteration. Therefore, broader fault damage zones preferentially facilitate fluid transport. This cools the mantle and reduces the potential for large earthquakes at localized barriers that divide the transform into shorter asperity regions, limiting earthquake magnitudes on the transform.

PREM-like velocity structure in the outermost core from global SKS and ScS waveform modeling

Seismic structure in the Earth's outermost core is important for our understanding of thermochemical stratification in the outer core, and is yet heavily debated. Here we study the compressional velocity structure in the Earth's outermost core based on waveform modeling of a unique untapped SKS and ScS dataset near bifurcation distances, collected from global seismic arrays for earthquakes occurring from 2000 to 2020. Using the SKS-ScS array dataset minimizes the effects of many uncertainties associated with earthquake source parameters and seismic heterogeneities in the mantle, and affords an opportunity to study and assess the seismic structure in the outermost core. We study outer core structure by testing two end-member models: 1) the D″ model that attributes any anomalous seismic observations to the effects of the lowermost mantle structure and is paired with a PREM (the Preliminary Reference Earth's Model) structure in the outermost core and 2) the outer core model that is paired with either a PREM or a tomographic structure in the lowermost mantle and attributes other unexplained seismic signals to an outer core structure. The results of the outer core models present unreasonable large lateral variations of >3.1% in the outermost core, while the inferred D″ models exhibit large-scale seismic anomalies that are consistent with the tomographic models and small-scale anomalies that are confirmed by further analysis of the seismic array data. Our analyses suggest a PREM-like seismic velocity structure and a lack of strong thermochemical anomalies in the topmost ∼200 km of the outer core, placing bounds on possible thermal and compositional conditions in the region of the Earth's outermost core. Our study also identifies the existence of small-scale seismic anomalies and sharp velocity variations in the lowermost mantle beneath the south coast of Alaska, northwestern Atlantic and the middle of Central America.