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
Shallow Crustal Shear Velocity and Vp/Vs Across Southern California: Joint Inversion of Short‐Period Rayleigh Wave Ellipticity, Phase Velocity, and Teleseismic Receiver Functions
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
- Award ID(s):
- 1753362
- PAR ID:
- 10360013
- Publisher / Repository:
- DOI PREFIX: 10.1029
- Date Published:
- Journal Name:
- Geophysical Research Letters
- Volume:
- 48
- Issue:
- 15
- ISSN:
- 0094-8276
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
Abstract Through the Alaska Transportable Array deployment of over 200 stations, we create a 3‐D tomographic model of Alaska with sensitivity ranging from the near surface (<1 km) into the upper mantle (~140 km). We perform a Markov chain Monte Carlo joint inversion of Rayleigh wave ellipticity and phase velocities, from both ambient noise and earthquake measurements, along with receiver functions to create a shear wave velocity model. We also use a follow‐up phase velocity inversion to resolve interstation structure. By comparing our results to previous tomography, geology, and geophysical studies we are able to validate our findings and connect localized near‐surface studies with deeper, regional models. Specifically, we are able to resolve shallow basins, including the Copper River, Cook Inlet, Yukon Flats, Nenana, and a variety of other shallower basins. Additionally, we gain insight on the interaction between the upper mantle wedge, asthenosphere, and active and nonactive volcanism along the Aleutians and Denali volcanic gap, respectively. We observe thicker crust beneath the Brooks Range and south of the Denali fault within the Wrangellia Composite Terrane and thinner crust in the Yukon Composite Terrane in interior Alaska. We also gain new perspective on the Wrangell Volcanic Field and its interaction between surrounding asthenosphere and the Yakutat Terrane.more » « less
-
Abstract The Formosa array, with 137 broadband seismometers and ∼5 km station spacing, was deployed recently in Northern Taiwan. Here by using eight months of continuous ambient noise records, we construct the first high‐resolution three‐dimensional (3‐D) shear wave velocity model of the crust in the area. We first calculate multi‐component cross‐correlations to extract robust Rayleigh wave signals. We then determine phase velocity maps between 3 and 10 s periods using Eikonal tomography and measure Rayleigh wave ellipticity at each station location between 2 and 13 s periods. For each location, we jointly invert the two types of Rayleigh wave measurements with a Bayesian‐based inversion method for a one‐dimensional shear wave velocity model. All piecewise continuous one‐dimensional models are then used to construct the final 3‐D model. Our 3‐D model reveals upper crustal structures that correlate well with surface geological features. Near the surface, the model delineates the low‐velocity Taipei and Ilan Basins from the adjacent fast‐velocity mountainous areas, with basin geometries consistent with the results of previous geophysical exploration and geological studies. At a greater depth, low velocity anomalies are observed associated with the Linkou Tableland, Tatun Volcano Group, and a possible dyke intrusion beneath the Southern Ilan Basin. The model also provides new geometrical constraints on the major active fault systems in the area, which are important to understand the basin formation, orogeny dynamics, and regional seismic hazard. The new 3‐D shear wave velocity model allows a comprehensive investigation of shallow geologic structures in the Northern Taiwan.more » « less
-
SUMMARY The near-surface seismic structure (to a depth of about 1000 m), particularly the shear wave velocity (VS), can strongly affect the propagation of seismic waves and, therefore, must be accurately calibrated for ground motion simulations and seismic hazard assessment. The VS of the top (<300 m) crust is often well characterized from borehole studies, geotechnical measurements, and water and oil wells, while the velocities of the material deeper than about 1000 m are typically determined by tomography studies. However, in depth ranges lacking information on shallow lithological stratification, typically rock sites outside the sedimentary basins, the material parameters between these two regions are typically poorly characterized due to resolution limits of seismic tomography. When the alluded geological constraints are not available, models, such as the Southern California Earthquake Center (SCEC) Community Velocity Models (CVMs), default to regional tomographic estimates that do not resolve the uppermost VS values, and therefore deliver unrealistically high shallow VS estimates. The SCEC Unified Community Velocity Model (UCVM) software includes a method to incorporate the near-surface earth structure by applying a generic overlay based on measurements of time-averaged VS in top 30 m (VS30) to taper the upper part of the model to merge with tomography at a depth of 350 m, which can be applied to any of the velocity models accessible through UCVM. However, our 3-D simulations of the 2014 Mw 5.1 La Habra earthquake in the Los Angeles area using the CVM-S4.26.M01 model significantly underpredict low-frequency (<1 Hz) ground motions at sites where the material properties in the top 350 m are significantly modified by the generic overlay (‘taper’). On the other hand, extending the VS30-based taper of the shallow velocities down to a depth of about 1000 m improves the fit between our synthetics and seismic data at those sites, without compromising the fit at well-constrained sites. We explore various tapering depths, demonstrating increasing amplification as the tapering depth increases, and the model with 1000 m tapering depth yields overall favourable results. Effects of varying anelastic attenuation are small compared to effects of velocity tapering and do not significantly bias the estimated tapering depth. Although a uniform tapering depth is adopted in the models, we observe some spatial variabilities that may further improve our method.more » « less
-
Abstract We construct a 3D shear velocity model of the Salt Lake Valley using Rayleigh waves excited by the 31 March 2020 Mw 6.5 central Idaho earthquake recorded on a 168-station temporary nodal geophone network and the 49-station permanent regional network. The temporary array—deployed in response to the March 18 Mw 5.7 Magna earthquake—serendipitously recorded clear surface waves between 10 and 20 s period from the Idaho event at ∼500 km epicentral distance, from which we measure both Rayleigh wave phase velocity and ellipticity (H/V ratio). In addition, we employ multicomponent earthquake coda cross correlation to extend the measurements down to 5 s period. Because Rayleigh wave ellipticity features outstanding shallow sensitivity, we invert for a 3D upper crust VS model of the Salt Lake Valley. Our model shows basin structure in general agreement with and complements the current Community Velocity Model, which is mostly constrained by borehole and gravity measurements. Our model thus provides critical information for future earthquake hazard assessment studies, which require detailed shallow velocity structure.more » « less
