Localized double-array stacking analysis of PcP: D″ and ULVZ structure beneath the Cocos plate, Mexico, central Pacific, and north Pacific

Alexander R. Hutko; Thorne Lay; Justin Revenaugh

doi:10.1016/j.pepi.2008.11.003

Home/Archive/USGS

USGS · 70034980

Localized double-array stacking analysis of PcP: D″ and ULVZ structure beneath the Cocos plate, Mexico, central Pacific, and north Pacific

Alexander R. Hutko·Thorne Lay·Justin Revenaugh

View original source ↗Export metadata ↓

Abstract

A large, high quality P-wave data set comprising short-period and broadband signals sampling four separate regions in the lowermost mantle beneath the Cocos plate, Mexico, the central Pacific, and the north Pacific is analyzed using regional one-dimensional double-array stacking and modelling with reflectivity synthetics. A data-screening criterion retains only events with stable PcP energy in the final data stacks used for modelling and interpretation. This significantly improves the signal stacks relative to including unscreened observations, allows confident alignment on the PcP arrival and allows tight bounds to be placed on P-wave velocity structure above the core–mantle boundary (CMB). The PcP reflections under the Cocos plate are well modelled without any ultra-low velocity zone from 5 to 20°N. At latitudes from 15 to 20°N, we find evidence for two P-wave velocity discontinuities in the D ″ region. The first is &sim;182 km above the CMB with a δln V p of +1.5%, near the same depth as a weaker discontinuity (<+0.5%) observed from 5 to 15°N in prior work. The other reflector is &sim;454 km above the CMB, with a δln V p of +0.4%; this appears to be a shallower continuation of the joint P- and S-wave discontinuity previously detected south of 15° N, which is presumed to be the perovskite to post-perovskite phase transition. The data stacks for paths bottoming below Mexico have PcP images that are well matched with the simple IASP91 structure, contradicting previous inferences of ULVZ presence in this region. These particular data are not very sensitive to any D ″ discontinuities, and simply bound them to be <&sim;2%, if present. Data sampling the lowermost mantle beneath the central Pacific confirm the presence of a &sim;15-km thick ultra-low velocity zone (ULVZ) just above the CMB, with δln V p and δln V s of around −3 to −4% and −4 to −8%, respectively. The ULVZ models predict previous S-wave data stacks well. The data for this region indicate laterally varying V p discontinuities in D ″, with one subregion having a δln V p of 0.5% 140 km above the CMB. Beneath the north Pacific, the PcP arrivals are compatible with only weak ULVZ (δln V p &sim; 0 to −3%), and there is a weak D ″ reflector with δln V p = 0.5%, near 314 km above the CMB. These results indicate localized occurrence of detectable ULVZ structures rather than ubiquitous ULVZ structure and emphasize the distinctiveness between the large low shear velocity province under the central Pacific and circum-Pacific regions.

Explore related subjects

Physics of the Earth and Planetary Interiors Article

Source-reported bounding extent: -42.0329743324414° to 48.574789910928864° latitude; -175.60546875° to -68.37890625° longitude. This indicates report coverage, not an exact sampling location. View area on OpenStreetMap.

Keep this discovery

Explore connections, maps & timelines

BibTeX RIS

Alexander R. Hutko, Thorne Lay, Justin Revenaugh. 2009. Localized double-array stacking analysis of PcP: D″ and ULVZ structure beneath the Cocos plate, Mexico, central Pacific, and north Pacific. https://doi.org/10.1016/j.pepi.2008.11.003

Cite the original work for its findings. Save a collection to share your selection of sources.

KEEP EXPLORING

Related USGS reports

USGS2021

Electrical conductivity of the lithosphere-asthenosphere system

Electromagnetic geophysical methods image the electrical conductivity of the subsurface. Electrical conductivity is an intrinsic material property that is sensitive to temperature, composition, porosity, volatile and/or melt content, and other physical properties relevant to the solid Earth. Therefore, imaging the electrical structure of the crust and mantle yields valuable information on the physical and chemical state of the lithosphere-asthenosphere system. Here we explore the viability of the passive magnetotelluric (MT) method for constraining upper mantle properties. We approach this problem in four successive steps: 1) review the electrical conductivity behavior of relevant materials; 2) predict the bulk electrical conductivity structure of oceanic and continental lithosphere for a suite of representative physical states; 3) generate synthetic MT data from the conductivity predictions; 4) compare and discuss the conductivity predictions and the synthetic data with select case studies from oceanic and continental settings. Our aim is to clarify the uncertainties associated with drawing inferences from electrical conductivity observations and ultimately to provide a basis for assigning confidence levels to interpretations.

Physics of the Earth and Planetary Interiors↗

USGS2021

Slip distribution and rupture history of the August 11, 2012, double earthquakes in Ahar – Varzaghan, Iran, using joint inversion of teleseismic broadband and local strong motion data

We use combined teleseismic and strong motion data sets to investigate finite-fault slip models for a double of earthquakes that occurred on August 11, 2012, in northwestern Iran near the cities of Ahar and Varzaghan. The data include teleseismic P-waveforms retrieved from broadband seismic stations located between 30°–94° from the earthquakes and local strong motion data recorded by the Iran Strong Motion Network, installed and operated by the Building and Housing Research Centre. We first invert teleseismic P-waveforms and local strong motion data separately. For the first event (12:23 UTC), the teleseismic broadband inversion yields a somewhat deeper and simpler distribution of slip than the local strong motion inversion. The strong motion inversion results in a more complex distribution because of higher frequency content but can also be influenced by complexities in the propagation path. For the second event (12:34 UTC), the slip distribution from strong motion data is more similar to the teleseismic result and shows a simple slip area with a small relative movement to the west. To resolve the differences between the results of these two data sets and obtain a better constrained slip model, we perform a joint inversion of teleseismic broadband and local strong motion data. The joint inversion for the first event shows two asperities with a maximum slip of 3.9 m up- dip from the hypocenter and extending to the west between depths of 1 and 5 km. A second narrower high-slip area is seen just above the hypocenter from 6 to 10 km depth. The total moment for this earthquake is calculated to be M o = 3.8 × 10 25 dyn-cm (3.8 × 10 18 N.m) (M w 6.4). For the second event, the results of the joint inversion show a simple slip distribution that is mainly confined in a single patch around the hypocenter with a depth range from about 10 to 13 km and maximum slip of 1.9 m. We compute a total seismic moment of M o = 1.6 × 10 25 dyn-cm (1.6 × 10 18 N.m) (M w 6.1) for the second event. The largest stress drops for the first event occur above the hypocenter with an average stress drop over the rupture area of 120 bar (12 Mpa). For the second event, the maximum stress drop occurs at the reported focal depth with an average stress drop over the rupture area of 80 bar (8 Mpa).

East Anatolian Fault↗

USGS2018

The 1952 Kern County, California earthquake: A case study of issues in the analysis of historical intensity data for estimation of source parameters

Seismic intensity data based on first-hand accounts of shaking give valuable insight into historical and early instrumental earthquakes. Comparing an observed intensity distribution to intensity-prediction models based on modern calibration events allows the magnitude to be estimated for many historic earthquakes. Magnitude estimates can also potentially be refined for earthquakes for which limited instrumental data are available. However, the complicated nature of macroseismic data and the methods used to collect and interpret the data introduce significant uncertainties. In this paper, we illustrate these challenges and possible solutions using the 1952 Kern County, California, earthquake as a case study. Published estimates of its magnitude vary from M W 7.2–7.5, making it possibly the second largest in California during the 20th century. We considered over 1100 first-hand reports of shaking, supplemented with other data, and inferred the magnitude in several ways using intensity prediction equations, yielding a preferred intensity magnitude M I 7.2 ± 0.2, where the uncertainty reflects our judgement. The revised intensity distribution reveals stronger shaking on the hanging wall, south of the surface expression of the White Wolf fault, than on the footwall. Characterizing the magnitude and shaking distribution of this early instrumental earthquake can help improve estimation of the seismic hazard of the region. Such reinterpreted intensities for historic earthquakes, combined with U.S. Geological Survey (USGS) Did You Feel It? data for more recent events, can be used to produce a uniform shaking dataset with which earthquake hazard map performance can be assessed.

California↗