Temperature gradients in a low-shear-velocity province in the lowermost mantle (D'' region) beneath the central Pacific Ocean were inferred from the observation of a rapid S-wave velocity increase overlying a rapid decrease. These paired seismic discontinuities are attributed to a phase change from perovskite to post-perovskite and then back to perovskite as the temperature increases with depth. Iron enrichment could explain the occurrence of post-perovskite several hundred kilometers above the core-mantle boundary in this warm, chemically distinct province. The double phase-boundary crossing directly constrains the lowermost mantle temperature gradients. Assuming a standard but unconstrained choice of thermal conductivity, the regional core-mantle boundary heat flux (approximately 85 +/- 25 milliwatts per square meter), comparable to the average at Earth's surface, was estimated, along with a lower bound on global core-mantle boundary heat flow in the range of 13 +/- 4 terawatts. Mapped velocity-contrast variations indicate that the lens of post-perovskite minerals thins and vanishes over 1000 kilometers laterally toward the margin of the chemical distinct region as a result of a approximately 500-kelvin temperature increase.

[1] Anomalous boundary layer structure at the core-mantle boundary (CMB) is investigated using a global set of broadband SKS and SPdKS waves from permanent and portable broadband seismometer arrays. SPdKS is an SKS wave that intersects the CMB at the critical angle for ScP, thus initiating a diffracted P wave (P diff ) along the CMB at the core entry and exit locations. The waveshape and timing of SPdKS data are analyzed relative to SKS, with some SPdKS data showing significant delays and broadening. Broadband data from several hundred deep focus earthquakes were analyzed; retaining data with simple sources and high signal-to-noise ratios resulted in 53 high-quality earthquakes. For each earthquake an empirical source was constructed by stacking pre-SPdKS distance range SKS pulses ($90°-100°). These were utilized in our synthetic modeling process, whereby reflectivity synthetic seismograms are produced for three classes of models: (1) mantle-side ultralow-velocity zones (UVLZ), (2) underside CMB core rigidity zones, and (3) core-mantle transition zones. For ULVZ structures, ratios of P-to-S velocity reductions of 1:1 and 1:3 are explored, where 1:3 is appropriate for a partial melt origin of ULVZ. Over 330 unique CMB boundary layer models have been constructed and tested, corroborating previous work suggesting strong trade-offs between the three model spaces. We produce maps of inferred boundary layer structure from the global data and find evidence for extremely fine-scale heterogeneity where our wave path sampling is the densest. While uncertainties are present relating to the source versus receiver sides of the SPdKS wave path geometry, our data are consistent with the hypothesis that ULVZ presence (or absence) correlates with reduced (or average) heterogeneity in the overlying mantle.