An official website of the United States government
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 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
-
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
-
The Santorini arc volcano in the Hellenic subduction zone has a history of caldera-forming Plinian eruptions, most recently in the Late Bronze Age 3.4 kya, and it remains volcanically active. To inform volcanic hazard assessments, it is crucial to understand where melt is distributed. The PROTEUS experiment in 2015 recorded >14,000 controlled marine sound sources on 165 land and seafloor seismic stations. Tomographic inversion of this data revealed low P-wave velocities in the upper 4 kilometers beneath the caldera and nearby Kolumbo seamount interpreted as the magma system (McVey et al., 2020; Chrapkewiecz et al., 2022). However, structure of the magma system was only determined in the upper (<4-6km) crust and melt content is only weakly constrained. In this study we improve constraints on the deeper magma system and subsurface melt content with a tomographic P and S wave velocity structure. To do so, we add to the inverse problem arrival times from ~1500 local earthquakes with magnitudes from 0.5 to 3.0 that occurred between 5 and 20 km depth. The events were recorded on 142 3-component ocean bottom and island seismic stations that span the seafloor ~60 km west and east of the island and the nearby islands. Results beneath Santorini and Kolumbo suggest that the upper crustal magma reservoirs extend deeper than previously found, and we image a high Vp layer (~5-8 km) under the magma reservoir at Kolumbo. We identify this layer as strong, cooled, intruded magma and correlate it to the location of earthquakes, within which, swarms of rapidly upward propagating seismicity support prior inferences of melt conduits traversing a rheologically strong layer (Schmid et al, 2022). We give values for melt content of the upper crustal reservoirs using a scaled Vp/Vs model. Since the number of arrivals, apriori assigned uncertainty, and differences in ray geometry can result in P and S waves with different resolving power, we use measured resolution to scale the Vs perturbations and create a more realistic Vp/Vs model. The addition of earthquake arrivals allows us to map the magma reservoirs beneath the Santorini-Kolumbo magma system to 8 km depth and identify regions of elevated melt content.more » « less
- Free Publicly Accessible Full Text
-
Accepted Manuscript1.0
- Journal Article:
- https://doi.org/10.1785/0220210223
