skip to main content


Title: Seismic Traveltime Tomography of Southern California Using Poisson‐Voronoi Cells and 20 Years of Data
Abstract

We derive new, 3D, isotropic models of seismic compressional and shear wavespeeds, Vp and Vs, respectively, their ratio, Vp/Vs, and a catalog of relocated earthquakes for Southern California from more than 10 million P‐ and S‐wave arrivals associated with over 0.3 million earthquakes that occurred between 2000 and 2020. We augment high‐quality analyst‐reviewed phase arrival picks from the Southern California Earthquake Data Center with S‐wave arrival picks obtained with an automated algorithm, and we derive new wavespeed models via traveltime tomography formulated using Poisson‐Voronoi cells (Fang et al., 2020,https://doi.org/10.1785/0220190141). The results contribute to improved regional wavespeed models, particularly the Vp/Vs model, and absolute event locations. The obtained models correlate well with regional geological features and yield more accurate synthetic waveforms than other regional models do for waves with periods shorter than 5 s in much of the modeled region. The derived event catalog exhibits tighter spatial clustering than the standard regional catalog, thereby helping to characterize subsurface features of major faults. The regional 1D averaged Vp/Vs ratio shows high values at shallow depths, decreases to a minimum at about 10 km, then increases again at greater depths below 15 km. Deep seismicity correlates well with regions of Vp/Vs ratio lower than 1.75, which may indicate an increased brittle‐to‐ductile transition depth with an important influence on crustal mechanics. The new wavespeed models and seismic catalog can be useful for various studies including analyses of seismicity patterns and simulations of crustal deformation and ground motion.

 
more » « less
NSF-PAR ID:
10444650
Author(s) / Creator(s):
 ;  ;  ;  
Publisher / Repository:
DOI PREFIX: 10.1029
Date Published:
Journal Name:
Journal of Geophysical Research: Solid Earth
Volume:
127
Issue:
5
ISSN:
2169-9313
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract

    Earthquake focal mechanisms, determined with P‐wave polarities and S/P amplitude ratios, are primary data for analyzing fault zone geometry, sense of slip, and the crustal stress field. Solving for the focal mechanisms of small earthquakes is often challenging because phase arrivals and first‐motion polarities are hard to be separated from noise. To overcome this challenge, we implement convolutional‐neural‐network algorithms (Ross, Meier, & Hauksson, 2018, Ross, Meier, Hauksson, & Heaton, 2018,https://doi.org/10.1029/2017jb015251,https://doi.org/10.1785/0120180080) to detect additional phases and polarities. Using both existing and these new data, we build a high‐quality focal mechanism catalog of 297,478 events that occurred from 1981 to 2021 in southern California with the HASH method of Hardebeck and Shearer (2002),https://doi.org/10.1785/0120010200, Hardebeck and Shearer (2003),https://doi.org/10.1785/0120020236. The new focal mechanism catalog is overall consistent with the standard catalog (Yang et al., 2012,https://doi.org/10.1785/0120110311) but includes 40% more focal mechanisms, and is more consistent with moment tensor solutions derived using waveform‐fitting methods. We apply the new catalog to identify changes in focal mechanism properties caused by the occurrences of large mainshocks such as the 2010Mw7.2 El Mayor‐Cucapah and 2019Mw7.1 Ridgecrest earthquakes. Such changes may be associated with co‐seismic stress drops, post‐seismic deformation processes, and static stress changes on a regional scale. The new high‐resolution catalog will contribute to improved understanding of the crustal stress field, earthquake triggering mechanisms, fault zone geometry, and sense of slip on the faults in southern California.

     
    more » « less
  2. Abstract

    Hundreds of earthquakes were recorded during a nine‐month ocean bottom seismometer deployment surrounding Lō'ihi submarine volcano, Hawai'i. The 12‐station ocean bottom seismometer network widened the aperture of earthquake detection around the Big Island, allowing better constraints on the location of seismicity offshore Hawai'i. Although this deployment occurred during a time of volcanic quiescence for Lō'ihi, it establishes an important basis for background seismicity of the volcano. Offshore seismicity during this study was dominated by events located in the mantle fault zone at depths of 25–40 km. These events reflect rupture on preexisting faults in the lower lithosphere caused by stresses induced by volcano loading and flexure of the Pacific Plate (Pritchard et al., 2007,https://doi.org/10.1111/j.1365‐246X.2006.03169.x; Wolfe et al., 2004,https://doi.org/10.1029/2003GC000618). Tomography was performed using double‐difference seismic tomography and showed shallow velocities to be slower than the regional velocity model (HG50; Klein, 1981,https://pubs.geoscienceworld.org/ssa/bssa/article/71/5/1503/118231/A‐linear‐gradient‐crustal‐model‐for‐south‐Hawaii). A broad, low‐velocity anomaly was observed from 20–40‐km depth, and is suggestive of the central plume conduit that supplies magma to Lō'ihi and the active volcanoes of the Big Island. A localized high‐velocity body is observed 4–6‐km depth beneath Lō'ihi's summit, extending 10 km to the north and south. Following Lō'ihi's active rift zones and crossing the summit, this high‐velocity body is characteristic of intrusive material. Two low‐velocity anomalies are observed below the oceanic crust, interpreted as melt accumulation beneath Lō'ihi and magmatic underplating beneath Hawai'i Island.

     
    more » « less
  3. Abstract

    This study provides a global assessment of the abundance of the major oxides in the deep continental crust. The combination of geochemistry and seismology better constrains the composition of the middle and lower continental crust better than either discipline can achieve alone. The inaccessible nature of the deep crust (typically >15 km) forces reliance on analog samples and modeling results to interpret its bulk composition, evolution, and physical properties. A common practice relates major oxide compositions of small‐ to medium‐scale samples (e.g., medium to high metamorphic grade terrains and xenoliths) to large scale measurements of seismic velocities (Vp, Vs, Vp/Vs) to determine the composition of the deep crust. We provide a framework for building crustal models with multidisciplinary constraints on composition. We present a global deep crustal model that documents compositional changes with depth and accounts for uncertainties in Moho depth, temperature, and physical and chemical properties. Our 3D compositional model of the deep crust uses the USGS Global Seismic Structure Catalog (Mooney, 2015) and a compilation of geochemical analyses on amphibolite and granulite facies lithologies (Sammon & McDonough, 2021,https://doi.org/10.1029/2021JB022791). We find a SiO2gradient from 61.2 ± 7.3 to 53.3 ± 4.8 wt.% from the middle to the base of the crust, with the equivalent lithological gradient ranging from quartz monzonite to gabbronorite. In addition, we calculate trace element abundances as a function of depth from their correlations with major oxides. From here, other lithospheric properties, such as Moho heat flux ( mW/m2), are derived.

     
    more » « less
  4. Abstract Puerto Rico is a highly seismically active island, where several damaging historical earthquakes have occurred and frequent small events persist. It situates at the boundary between the Caribbean and North American plates, featuring a complex fault system. Here, we investigate the seismotectonic crustal structure of the island by interpreting the 3D compressional-wave velocity VP and compressional- to shear-wave velocity ratio VP/VS models and by analyzing the distribution of the relocated earthquakes. The 3D velocity models are obtained by applying the simul2000 tomographic inversion algorithm based on the phase arrivals recorded by the Puerto Rico seismic network. We find high-VP and low-VP/VS anomalies in the eastern and central province between the Great Northern Puerto Rico fault zone and the Great Southern Puerto Rico fault zone, correlating with the Utuado pluton. Further, there are low-VP anomalies beneath both the Great Southern Puerto Rico fault zone and the South Lajas fault, indicating northerly dipping structures from the southwest to the northwest of the island. We relocate 19,095 earthquakes from May 2017 to April 2021 using the new 3D velocity model and waveform cross-correlation data. The relocated seismicity shows trends along the Investigator fault, the Ponce faults, the Guayanilla rift, and the Punta Montalva fault. The majority of the 2019–2021 Southwestern Puerto Rico earthquakes are associated with the Punta Montalva fault. Earthquakes forming 17° northward-dipping structures at various depths possibly manifest continuation of the Muertos trough, along which the Caribbean plate is being subducted beneath the Puerto Rico microplate. Our results show complex fault geometries of a diffuse fault network, suggesting possible subduction process accommodated by faults within a low-velocity zone. 
    more » « less
  5. Abstract

    Seven three‐component ocean bottom seismometers (OBS) of the Ocean Observatories Initiative (OOI) Cabled Array on top of Axial Seamount are continuously streaming data in real time to the Incorporated Research Institutions for Seismology (IRIS). The OBS array records earthquakes from the submarine volcano which last erupted on 24 April 2015, about 4 months after the array came online. The OBS data have proven crucial in providing insight into the volcano structure and dynamics (Wilcock et al., 2016,https://doi.org/10.1126/science.aah5563). We implemented a real‐time double‐difference (RT‐DD) monitoring system that automatically computes high‐precision (tens of meters) locations of new earthquakes. The system's underlying double‐difference base catalog includes nearly 100,000 earthquakes and was computed using kurtosis phase onset picks, cross‐correlation phase delay times, and 3‐DPandSvelocity models to predict the data. The relocations reveal the fine‐scale structures of long‐lived, narrow (<200 m wide), outward dipping, convex faults on the east and west walls of the caldera that appear to form a figure 8‐shaped ring fault system. These faults accommodate stresses caused by the inflation of magma prior to and deflation during eruptions. The east fault is segmented and pulled apart in east‐west direction due to its interaction with the Juan de Fuca Ridge, which at this location forms an overlapping spreading center. The RT‐DD system enables the monitoring and rapid analysis of variations in fine‐scale seismic and fault properties and has the potential to improve prediction of timing and location of the next Axial eruption expected to occur in the 2022–2023 time frame.

     
    more » « less