skip to main content


Title: 3D Seismic Velocity Models for Alaska from Joint Tomographic Inversion of Body-Wave and Surface-Wave Data
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
Award ID(s):
1829440 1251971
NSF-PAR ID:
10196405
Author(s) / Creator(s):
; ; ; ; ; ; ; ;
Date Published:
Journal Name:
Seismological Research Letters
ISSN:
0895-0695
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract

    The along‐strike variations of the velocity, thickness, and dip of subducting slabs and the volcano distribution have been observed globally. It is, however, unclear what controls the distribution of volcanoes and the associated magma generation. With the presence of nonuniform volcanism, the Aleutian‐Alaska subduction zone (AASZ) is an ideal place to investigate subduction segmentation and its relationship with volcanism. Using full‐wave ambient noise tomography, we present a high‐resolution 3‐D shear wave velocity model of the AASZ for the depths of 15–110 km. The velocity model reveals the distinct high‐velocity Pacific slab, the thicker, flatter, and more heterogeneous Yakutat slab, and the northeasterly dipping Wrangell slab. We observe low velocities within the uppermost mantle (at depth <60 km) below the Aleutian arc volcanoes, representing partial melt accumulation. The large crustal low‐velocity anomaly beneath the Wrangell volcanic field suggests a large magma reservoir, likely responsible for the clustering of volcanoes. The Denali volcanic gap is above an average‐velocity crust but an extremely fast mantle wedge, suggesting the lack of subsurface melt. This is in contrast with the lower‐velocity back‐arc mantle beneath the adjacent Buzzard Creek‐Jumbo Dome volcanoes to the east. The back‐arc low velocities associated with the Pacific, the eastern Yakutat, and the Wrangell slabs may reflect subduction‐driven mantle upwelling. The structural variation of the downgoing slabs and the overriding plate explains the change of volcanic activity along the AASZ. Our findings demonstrate the combined role of the subducting slab and the overriding plate in controlling the characteristics of arc magmatism.

     
    more » « less
  2. Abstract

    We report newPandSwave velocity models of the upper mantle beneath southern Africa using data recorded on seismic stations spanning the entire subcontinent. Beneath most of the Damara Belt, including the Okavango Rift, our models show lower than average velocities (−0.8% Vp; −1.2% Vs) with an abrupt increase in velocities along the terrane's southern margin. We attribute the lower than average velocities to thinner lithosphere (~130 km thick) compared to thicker lithosphere (~200 km thick) immediately to the south under the Kalahari Craton. Beneath the Etendeka Flood Basalt Province, higher than average velocities (0.25% Vp; 0.75% Vs) indicate thicker and/or compositionally distinct lithosphere compared to other parts of the Damara Belt. In the Rehoboth Province, higher than average velocities (0.3% Vp; 0.5% Vs) suggest the presence of a microcraton, as do higher than average velocities (1.0% Vp; 1.5% Vs) under the Southern Irumide Belt. Lower than average velocities (−0.4% Vp; −0.7% Vs) beneath the Bushveld Complex and parts of the Mgondi and Okwa terranes are consistent with previous studies, which attributed them to compositionally modified lithosphere resulting from Precambrian magmatic events. There is little evidence for thermally modified upper mantle beneath any of these terranes which could provide a source of uplift for the Southern African Plateau. In contrast, beneath parts of the Irumide Belt in southern and central Zambia and the Mozambique Belt in central Mozambique, deep‐seated low velocity anomalies (−0.7% Vp; −0.8% Vs) can be attributed to upper mantle extensions of the African superplume structure.

     
    more » « less
  3. SUMMARY EarthScope's USArray seismic component provided unprecedented coverage of the contiguous United States and has therefore spurred significant advances in tomographic imaging and geodynamic modelling. Here, we present a new global, radially anisotropic shear wave velocity tomography model to investigate upper mantle structure and North American Plate dynamics, with a focus on the contiguous United States. The model uses a data-adaptive mesh and traveltimes of both surface waves and body waves to constrain structure in the crust and mantle in order to arrive at a more consistent representation of the subsurface compared to what is provided by existing models. The resulting model is broadly consistent with previous global models at the largest scales, but there are substantial differences under the contiguous United States where we can achieve higher resolution. On these regional scales, the new model contains short wavelength anomalies consistent with regional models derived from USArray data alone. We use the model to explore the geometry of the subducting Farallon Slab, the presence of upper mantle high velocity anomalies, low velocity zones in the central and eastern United States and evaluate models of dynamic topography in the Cordillera. Our models indicate a single, shallowly dipping, discontinuous slab associated with the Farallon Plate, but there are remaining imaging challenges. Inferring dynamic topography from the new model captures both the long-wavelength anomalies common in global models and the short-wavelength anomalies apparent in regional models. Our model thus bridges the gap between high-resolution regional models within the proper uppermost mantle context provided by global models, which is crucial for understanding many of the fundamental questions in continental dynamics. 
    more » « less
  4. Abstract

    We investigated the seismic velocity structure of the Hikurangi margin in New Zealand to uncover the physical features of the subduction zone and explore the relationships between microearthquake seismicity, seismic velocity structure, and slow slip events. Using local earthquake tomography with data collected from both temporary ocean bottom seismometers and on‐land permanent seismic stations, we used the tomography code TomoFD to iteratively perform a damped least squared inversion of absolute P and S arrival times to obtain relocated hypocenters and generate 3D velocity models for Vp and Vp/Vs. The seismic tomography images show two high Vp/Vs anomalies, one offshore and adjacent to a subducted seamount and the other beneath the North Island of New Zealand. The ∼50‐km wide offshore anomaly extends ∼10 km beneath the plate interface and lies directly beneath the area that slipped at least 50 mm during the 2 week‐long 2014 slow slip event. High Vp/Vs values may be related to high pore fluid pressures from subducted sediments, and such increases in pore fluid pressures have been suggested to trigger the occurrence of slow slip events in active subduction zones. The second onshore high Vp/Vs anomaly is located in the overlying plate and subducting slab and correlates with areas suggested by other geophysical techniques to be rich in fluids. Our seismic imaging supports interpretations that subduction processes in the Hikurangi margin are highly dependent on physical features such as subducted seamounts and fluid‐rich sediments.

     
    more » « less
  5. Abstract

    Shear‐wave splitting observations can provide insight into mantle flow, due to the link between the deformation of mantle rocks and their direction‐dependent seismic wave velocities. We identify anisotropy in the Cook Inlet segment of the Alaska subduction zone by analyzing splitting parameters of S waves from local intraslab earthquakes between 50 and 200 km depths, recorded from 2015–2017 and emphasizing stations from the Southern Alaska Lithosphere and Mantle Observation Network experiment. We classify 678 high‐quality local shear‐wave splitting observations into four regions, from northwest to southeast: (L1b) splitting measurements parallel to Pacific plate motion, (L1a) arc‐perpendicular splitting pattern, (L2) sharp transition to arc‐parallel splitting, and (L3) splitting parallel to Pacific plate motion. Forward modeling of splitting from various mantle fabrics shows that no one simple model fully explains the observed splitting patterns. An A‐type olivine fabric with fast direction dipping 45° to the northwest (300°)—aligned with the dipping slab—predicts fast directions that fit L1a observations well, but not L2. The inability of the forward model fabrics to fit all the observed splitting patterns suggests that the anisotropy variations are not due to variable ray angles, but require distinct differences in the anisotropy regime below the arc, forearc, and subducting plate.

     
    more » « less