Abstract We present two new seismic velocity models for Alaska from joint inversions of body-wave and ambient-noise-derived surface-wave data, using two different methods. Our work takes advantage of data from many recent temporary seismic networks, including the Incorporated Research Institutions for Seismology Alaska Transportable Array, Southern Alaska Lithosphere and Mantle Observation Network, and onshore stations of the Alaska Amphibious Community Seismic Experiment. The first model primarily covers south-central Alaska and uses body-wave arrival times with Rayleigh-wave group-velocity maps accounting for their period-dependent lateral sensitivity. The second model results from direct inversion of body-wave arrival times and surface-wave phase travel times, and covers the entire state of Alaska. The two models provide 3D compressional- (VP) and shear-wave velocity (VS) information at depths ∼0–100 km. There are many similarities as well as differences between the two models. The first model provides a clear image of the high-velocity subducting plate and the low-velocity mantle wedge, in terms of the seismic velocities and the VP/VS ratio. The statewide model provides clearer images of many features such as sedimentary basins, a high-velocity anomaly in the mantle wedge under the Denali volcanic gap, low VP in the lower crust under Brooks Range, and low velocities at the eastern edge of Yakutat terrane under the Wrangell volcanic field. From simultaneously relocated earthquakes, we also find that the depth to the subducting Pacific plate beneath southern Alaska appears to be deeper than previous models.
more »
« less
Crustal Architecture of Puerto Rico Using Body-Wave Seismic Tomography and High-Resolution Earthquake Relocation
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
- PAR ID:
- 10345481
- Date Published:
- Journal Name:
- Seismological Research Letters
- Volume:
- 93
- Issue:
- 2A
- ISSN:
- 0895-0695
- Page Range / eLocation ID:
- 555 to 566
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
SUMMARY We present our estimations and comparisons of the in situ Vp/Vs ratios and seismicity characteristics for the Parkfield segment of the San Andreas fault in northern California and the San Jacinto Fault Zone and its adjacent regions in southern California. Our results show that the high-resolution in situ Vp/Vs ratios are much more complex than the tomographic Vp/Vs models. They show similar variation patterns to those in the tomographic Vp models, indicating that Vp/Vs ratios are controlled by material properties but are also strongly influenced by fluid contents. In Parkfield, we observe velocity contrasts between the creeping and locked sections. In southern California, we see small-scale anomalous Vp/Vs variation patterns, especially where fault segments intersect, terminate and change orientations. In addition, our investigation confirms that the seismicity in Parkfield is more repeatable than in southern California. However, the earthquakes in the southernmost portion of the San Andreas fault, the trifurcation area of the San Jacinto Fault Zone and the Imperial fault are as much likely falling into clusters as those in Parkfield. The correlation of highly similar events with anomalous in situ Vp/Vs ratios supports the important role of fluids in the occurrence of repeating earthquakes. The high-resolution Vp/Vs ratio estimation method and the corresponding results are helpful for revealing roles of fluids in driving earthquake, fault interaction and stress distribution in fault zones.more » « less
-
Abstract The influence of Al substitution on the elastic properties of stishovite and its transition to post-stishovite is of great importance for interpreting the seismic wave velocities of subducted mid-ocean ridge basalt (MORB) within the mantle transition zone and the lower mantle. However, atomistic mechanisms of Al substitution effects on the transition and its associated elasticity remain debated. Here synchrotron single-crystal X-ray diffraction measurements have been performed at room temperature on Al1.3-SiO2 (1.3 mol% Al in the chemical formula of Si0.965(3)Al0.041(1)O2H0.017(4)) and Al2.1-SiO2 (2.1 mol% Al in Si0.948(2)Al0.064(1)O2H 0.018(3)) crystals in diamond anvil cells with Boehler-Almax designed anvils up to 38.0 GPa and 28.5 GPa, respectively. Refinements of the diffraction patterns show that a transformation from stishovite (space group P42/mnm; No. 136) to CaCl2-typed post-stishovite (space group Pnnm; No. 58) is accompanied by splitting of O coordinates. The Al substitution in stishovite results in a faster decrease in the O coordinate, softer apical (Si,Al)-O bonds, and a softer and less distorted (Si,Al)O6 octahedron under compression. This leads to reduced adiabatic bulk modulus (KS), shear modulus (G), shear wave velocity (VS), and compressional wave velocity (VP) in the stishovite phase, explaining seismic wave perturbations in the mantle transition zone. Together with Raman data, Landau theory modeling shows that Al substitution increases the order parameter and excess free energy, stabilizing the post-stishovite phase at lower pressures. Correlation between elasticity and octahedral distortion index (D) reveals that at certain D, the Al substitution reduces KS, G, VS, and VP of the stishovite phase while increasing G, VS, and VP of the post-stishovite phase. Importantly, the maximum shear reduction is slightly enhanced at D = 0.00620(9) at the transition point. Our results help explain the seismically observed small-scale VS anomalies beneath subduction regions in the shallow lower mantle where Al,H-bearing stishovite undergoes the post-stishovite transition.more » « less
-
SUMMARY Following reanalysis of data from eight seismic networks that operated in the region surrounding the Pampean flat slab during the past several decades, we generated 3-D images of Vp, Vs and $$V_{\rm p}/V_{\rm s}$$ from a combination of arrival times of P and S waves from local earthquakes, and Rayleigh wave dispersion curves from both ambient noise and existing shear wave models. Among the robust features in these images is a low velocity, root-like structure that extends beneath the high Andes to a deflection in the flat slab, which suggests the presence of an overthickened Andean crust rather than a hypothesized continental lithospheric root. Most of the larger scale features observed in both the subducted Nazca plate and the overriding continental lithosphere are related to the intense seismic activity in and around the Juan Fernandez Ridge Seismic Zone (JFRSZ). $$V_{\rm p}/V_{\rm s}$$ ratios beneath, within and above the JFRSZ are generally lower (~1.65–1.68) than those in the surrounding Nazca and continental lithosphere (~1.74–1.80). While the higher continental lithosphere ratios are due to reduced Vs and likely a result of hydration, the lower JFRSZ related ratios are due to reduced Vp and can be explained by increased silica and CO2 originating from beneath the slab, perhaps in concert with supercritical fluid located within the fracture and fault networks associated with the JFR. These and related features such as a region of high Vp and Vs observed at the leading edge of the JFRSZ are consistent with a basal displacement model previously proposed for the Laramide flat-slab event, in which the eroded base of the continental lithosphere accumulates as a keel at the front end of the flat slab while compressional horizontal stresses cause it to buckle. An initial concave up bend in the slab facilitates the infiltration of silica and CO2-rich melts from beneath the slab in a manner analogous to petit spot volcanism, while a second, concave down bend, releases CO2 and supercritical fluid into the overlying continental lithosphere.more » « less
-
Abstract Near‐surface seismic velocity structure plays a critical role in ground motion amplification during large earthquakes. In particular, the local Vp/Vs ratio strongly influences the amplitude of Rayleigh waves. Previous studies have separately imaged 3D seismic velocity and Vp/Vs ratio at seismogenic depth, but lack regional coverage and/or fail to constrain the shallowest structure. Here, we combine three datasets with complementary sensitivity in a Bayesian joint inversion for shallow crustal shear velocity and near‐surface Vp/Vs ratio across Southern California. Receiver functions–including with an apparent delayed initial peak in sedimentary basins, and long considered a nuisance in receiver function imaging studies–highly correlate with short‐period Rayleigh wave ellipticity measurements and require the inclusion of a Vp/Vs parameter. The updated model includes near‐surface low shear velocity more in line with geotechnical layer estimates, and generally lower than expected Vp/Vs outside the basins suggesting widespread shallow fracturing and/or groundwater undersaturation.more » « less
An official website of the United States government

