Seismic evidence for an 850 km thick low-velocity structure in the Earth's lowermost mantle beneath Kamchatka

2014 ◽  
Vol 41 (20) ◽  
pp. 7073-7079 ◽  
Author(s):  
Yumei He ◽  
Lianxing Wen ◽  
Tianyu Zheng
2021 ◽  
Author(s):  
Sima Mousavi ◽  
Hrvoje Tkalčić ◽  
Rhys Hawkins ◽  
Malcolm Sambridge

The core-mantle boundary (CMB) is the most extreme boundary within the Earth where the liquid, iron-rich outer core interacts with the rocky, silicate mantle. The nature of the lowermost mantle atop the CMB, and its role in mantle dynamics, is not completely understood. Various regional studies have documented significant heterogeneities at different spatial scales. While there is a consensus on the long scale-length structure of the inferred S-wave speed tomograms, there are also notable differences stemming from different imaging methods and datasets. Here we aim to overcome over-smoothing and avoid over-fitting data for the case where the spatial coverage is sparse and the inverse problem ill-posed. Here we present an S-wave tomography model at global scale for the Lowermost Mantle (LM) using the Hierarchical Trans-dimensional Bayesian Inversion (HTDBI) framework, LM-HTDBI. Our HTDBI analysis of ScS-S travel times includes uncertainty, and the complexity of the model is deduced from the data itself through an implicit parameterization of the model space. Our comprehensive resolution estimates indicate that short-scale anomalies are significant and resolvable features of the lowermost mantle regardless of the chosen mantle-model reference to correct the travel times above the D’’ layer. The recovered morphology of the Large-Low-Shear-wave Velocity Provinces (LLSVPs) is complex, featuring small high-velocity patches among low-velocity domains. Instead of two large, unified, and smooth LLSVPs, the newly obtained images suggest that their margins are not uniformly flat.


1992 ◽  
Vol 29 (7) ◽  
pp. 1509-1529 ◽  
Author(s):  
Tianson Yuan ◽  
G. D. Spence ◽  
R. D. Hyndman

A combined multichannel seismic reflection and refraction survey was carried out in July 1988 to study the Tertiary sedimentary basin architecture and formation and to define the crustal structure and associated plate interactions in the Queen Charlotte Islands region. Simultaneously with the collection of the multichannel reflection data, refractions and wide-angle reflections from the airgun array shots were recorded on single-channel seismographs distributed on land around Hecate Strait and Queen Charlotte Sound. For this paper a subset of the resulting data set was chosen to study the crustal structure in Queen Charlotte Sound and the nearby subduction zone.Two-dimensional ray tracing and synthetic seismogram modelling produced a velocity structure model in Queen Charlotte Sound. On a margin-parallel line, Moho depth was modelled at 27 km off southern Moresby Island but only 23 km north of Vancouver Island. Excluding the approximately 5 km of the Tertiary sediments, the crust in the latter area is only about 18 km thick, suggesting substantial crustal thinning in Queen Charlotte Sound. Such thinning of the crust supports an extensional mechanism for the origin of the sedimentary basin. Deep crustal layers with velocities of more than 7 km/s were interpreted in the southern portion of Queen Charlotte Sound and beneath the continental margin. They could represent high-velocity material emplaced in the crust from earlier subduction episodes or mafic intrusion associated with the Tertiary volcanics.Seismic velocities of both sediment and upper crust layers are lower in the southern part of Queen Charlotte Sound than in the region near Moresby Island. Well velocity logs indicate a similar velocity variation. Gravity modelling along the survey line parallel to the margin provides additional constraints on the structure. The data require lower densities in the sediment and upper crust of southern Queen Charlotte Sound. The low-velocity, low-density sediments in the south correspond to high-porosity marine sediments found in wells in that region and contrast with lower porosity nonmarine sediments in wells farther north.


1997 ◽  
Vol 40 (1) ◽  
Author(s):  
G. Drakatos ◽  
G. Karantonis ◽  
G. N. Stavrakakis

The three-dimensional velocity structure of the crust in the Aegean sea and the surrounding regions (34.0º-42.OºN, 19.0ºE-29.0ºE) is investigated by inversion of about 10000 residuals of arrival times of P-wave from local events. The resulting velocity structure shows strong horizontal variations due to the complicated crustal structure and the variations of crustal thickness. The northern part of the region generally shows high velocities. In the inner part of the volcanic arc (Southern Aegean area), relatively low velocities are observed, suggesting a large-scale absorption of seismic energy as confirmed by the low seismicity of the region. A low velocity zone was observed along the subduction zone of the region, up to a depth of 4 km. The existence of such a zone could be due to granitic or other intrusions in the crust during the uplift of the region during Alpidic orogenesis.


1996 ◽  
Vol 39 (6) ◽  
Author(s):  
C. Chiarabba ◽  
A. Amato

In this paper we provide P-wave velocity images of the crust underneath the Apennines (Italy), focusing on the lower crustal structure and the Moho topography. We inverted P-wave arrival times of earthquakes which occurred from 1986 to 1993 within the Apenninic area. To overcome inversion instabilities due to noisy data (we used bulletin data) we decided to resolve a minimum number of velocity parameters, inverting for only two layers in the crust and one in the uppermost mantle underneath the Moho. A partial inversion of only 55% of the overall dataset yields velocity images similar to those obtained with the whole data set, indicating that the depicted tomograms are stable and fairly insensitive to the number of data used. We find a low-velocity anomaly in the lower crust extending underneath the whole Apenninic belt. This feature is segmented by a relative high-velocity zone in correspondence with the Ortona-Roccamonfina line, that separates the northern from the southern Apenninic arcs. The Moho has a variable depth in the study area, and is deeper (more than 37 km) in the Adriatic side of the Northern Apennines with respect to the Tyrrhenian side, where it is found in the depth interval 22-34 km.


2019 ◽  
Vol 219 (Supplement_1) ◽  
pp. S152-S166 ◽  
Author(s):  
Jonathan Wolf ◽  
Neala Creasy ◽  
Angelo Pisconti ◽  
Maureen D Long ◽  
Christine Thomas

SUMMARY Iceland represents one of the most well-known examples of hotspot volcanism, but the details of how surface volcanism connects to geodynamic processes in the deep mantle remain poorly understood. Recent work has identified evidence for an ultra-low velocity zone in the lowermost mantle beneath Iceland and argued for a cylindrically symmetric upwelling at the base of a deep mantle plume. This scenario makes a specific prediction about flow and deformation in the lowermost mantle, which can potentially be tested with observations of seismic anisotropy. Here we present an investigation of seismic anisotropy in the lowermost mantle beneath Iceland, using differential shear wave splitting measurements of S–ScS and SKS–SKKS phases. We apply our techniques to waves propagating at multiple azimuths, with the goal of gaining good geographical and azimuthal coverage of the region. Practical limitations imposed by the suboptimal distribution of global seismicity at the relevant distance ranges resulted in a relatively small data set, particularly for S–ScS. Despite this, however, our measurements of ScS splitting due to lowermost mantle anisotropy clearly show a rotation of the fast splitting direction from nearly horizontal for two sets of paths that sample away from the low velocity region (implying VSH > VSV) to nearly vertical for a set of paths that sample directly beneath Iceland (implying VSV > VSH). We also find evidence for sporadic SKS–SKKS discrepancies beneath our study region; while the geographic distribution of discrepant pairs is scattered, those pairs that sample closest to the base of the Iceland plume tend to be discrepant. Our measurements do not uniquely constrain the pattern of mantle flow. However, we carried out simple ray-theoretical forward modelling for a suite of plausible anisotropy mechanisms, including those based on single-crystal elastic tensors, those obtained via effective medium modelling for partial melt scenarios, and those derived from global or regional models of flow and texture development in the deep mantle. These simplified models do not take into account details such as possible transitions in anisotropy mechanism or deformation regime, and test a simplified flow field (vertical flow beneath the plume and horizontal flow outside it) rather than more detailed flow scenarios. Nevertheless, our modelling results demonstrate that our ScS splitting observations are generally consistent with a flow scenario that invokes nearly vertical flow directly beneath the Iceland hotspot, with horizontal flow just outside this region.


Geophysics ◽  
1995 ◽  
Vol 60 (6) ◽  
pp. 1917-1929 ◽  
Author(s):  
Joseph P. Stefani

Turning‐ray tomography is useful for estimating near‐surface velocity structure in areas where conventional refraction statics techniques fail because of poor data or lack of smooth refractor/velocity structure. This paper explores the accuracy and inherent smoothing of turning‐ray tomography in its capacity to estimate absolute near‐surface velocity and the statics times derived from these velocities, and the fidelity with which wavefields collapse to point diffractors when migrated through these estimated velocities. The method comprises nonlinear iterations of forward ray tracing through triangular cells linear in slowness squared, coupled with the LSQR linear inversion algorithm. It is applied to two synthetic finite‐ difference data sets of types that usually foil conventional refraction statics techniques. These models represent a complex hard‐rock overthrust structure with a low‐velocity zone and pinchouts, and a contemporaneous near‐shore marine trench filled with low‐ velocity unconsolidated deposits exhibiting no seismically apparent internal structure. In both cases velocities are estimated accurately to a depth of one‐ fifth the maximum offset, as are the associated statics times. Of equal importance, the velocities are sufficiently accurate to correctly focus synthetic wavefields back to their initial point sources, so migration/datuming applications can also use these velocities. The method is applied to a real data example from the Timbalier Trench in the Gulf of Mexico, which exhibits the same essential features as the marine trench synthetic model. The Timbalier velocity inversion is geologically reasonable and yields long and short wavelength statics that improve the CMP gathers and stack and that correctly align reflections to known well markers. Turning‐ray tomography estimates near‐surface velocities accurately enough for the three purposes of lithology interpretation, statics calculations, and wavefield focusing for shallow migration and datuming.


2020 ◽  
Author(s):  
Yan Cai ◽  
Jianping Wu

<p>North China Craton is the oldest craton in the world. It contains the eastern, central and western part. Shanxi rift and Taihang mountain contribute the central part. With strong tectonic deformation and intense seismic activity, its crust-mantle deformation and deep structure have always been highly concerned. In recent years, China Earthquake Administration has deployed a dense temporary seismic array in North China. With the permanent and temporary stations, we obtained the crust-mantle S-wave velocity structure in the central North China Craton by using the joint inversion of receiver function and surface wave dispersion. The results show that the crustal thickness is thick in the north of the Shanxi rift (42km) and thin in the south (35km). Datong basin, located in the north of the rift, exhibits large-scale low-velocity anomalies in the middle-lower crust and upper mantle; the Taiyuan basin and Linfen basin, located in the central part, have high velocities in the lower crust and upper mantle; the Yuncheng basin, in the southern part, has low velocities in the lower crust and upper mantle velocities, but has a high-velocity layer below 80 km. We speculate that an upwelling channel beneath the west of the Datong basin caused the low velocity anomalies there. In the central part of the Shanxi rift, magmatic bottom intrusion occurred before the tension rifting, so that the heated lithosphere has enough time to cool down to form high velocity. Its current lithosphere with high temperature may indicate the future deformation and damage. There may be a hot lithospheric uplift in the south of the Shanxi rift, heating the crust and the lithospheric mantle. The high-velocity layer in its upper mantle suggests that the bottom of the lithosphere after the intrusion of the magma began to cool down.</p>


Geophysics ◽  
1994 ◽  
Vol 59 (8) ◽  
pp. 1278-1289 ◽  
Author(s):  
William J. Lutter ◽  
Rufus D. Catchings ◽  
Craig M. Jarchow

We use a method of traveltime inversion of high‐resolution seismic data to provide the first reliable images of internal details of the Columbia River Basalt Group (CRBG), the subsurface basalt/sediment interface, and the deeper sediment/basement interface. Velocity structure within the basalts, delineated on the order of 1 km horizontally and 0.2 km vertically, is constrained to within ±0.1 km/s for most of the seismic profile. Over 5000 observed traveltimes fit our model with an rms error of 0.018 s. The maximum depth of penetration of the basalt diving waves (truncated by underlying low‐velocity sediments) provides a reliable estimate of the depth to the base of the basalt, which agrees with well‐log measurements to within 0.05 km (165 ft). We use image blurring, calculated from the resolution matrix, to estimate the aspect ratio of imaged velocity anomaly widths to true widths for velocity features within the basalt. From our calculations of image blurring, we interpret low velocity zones (LVZ) within the basalts at Boylston Mountain and the Whiskey Dick anticline to have widths of 4.5 and 3 km, respectively, within the upper 1.5 km of the model. At greater depth, the widths of these imaged LVZs thin to approximately 2 km or less. We interpret these linear, subparallel, low‐velocity zones imaged adjacent to anticlines of the Yakima Fold Belt to be brecciated fault zones. These fault zones dip to the south at angles between 15 to 45 degrees.


Geophysics ◽  
2001 ◽  
Vol 66 (6) ◽  
pp. 1925-1936 ◽  
Author(s):  
Moritz M. Fliedner ◽  
Robert S. White

We use the wide‐angle wavefield to constrain estimates of the seismic velocity and thickness of basalt flows overlying sediments. Wide angle means the seismic wavefield recorded at offsets beyond the emergence of the direct wave. This wide‐angle wavefield contains arrivals that are returned from within and below the basalt flows, including the diving wave through the basalts as the first arrival and P‐wave reflections from the base of the basalts and from subbasalt structures. The velocity structure of basalt flows can be determined to first order from traveltime information by ray tracing the basalt turning rays and the wide‐angle base‐basalt reflection. This can be refined by using the amplitude variation with offset (AVO) of the basalt diving wave. Synthetic seismogram models with varying flow thicknesses and velocity gradients demonstrate the sensitivity to the velocity structure of the basalt diving wave and of reflections from the base of the basalt layer and below. The diving‐wave amplitudes of the models containing velocity gradients show a local amplitude minimum followed by a maximum at a greater range if the basalt thickness exceeds one wavelength and beyond that an exponential amplitude decay. The offset at which the maximum occurs can be used to determine the basalt thickness. The velocity gradient within the basalt can be determined from the slope of the exponential amplitude decay. The amplitudes of subbasalt reflections can be used to determine seismic velocities of the overburden and the impedance contrast at the reflector. Combining wide‐angle traveltimes and amplitudes of the basalt diving wave and subbasalt reflections enables us to obtain a more detailed velocity profile than is possible with the NMO velocities of small‐offset reflections. This paper concentrates on the subbasalt problem, but the results are more generally applicable to situations where high‐velocity bodies overlie a low‐velocity target, such as subsalt structures.


Geophysics ◽  
2006 ◽  
Vol 71 (6) ◽  
pp. U77-U84 ◽  
Author(s):  
B. Bergman ◽  
A. Tryggvason ◽  
C. Juhlin

Reflection seismic imaging of the uppermost kilometer of crystalline bedrock is an important component in site surveys for locating potential storage sites for nuclear waste in Sweden. To obtain high-quality images, refraction statics are calculated using first-break traveltimes. These first-break picks may also be used to produce tomographic velocity images of the uppermost bedrock. In an earlier study, we presented a method applicable to data sets where the vast majority of shots are located in the bedrock below the glacial deposits, or cover, typical for northern latitudes. A by-product of this method was an estimate of the cover thickness from the receiver static that was introduced to sharpen the image. We now present a modified version of this method that is applicable for sources located in or on the cover, the general situation for nuclear waste site surveys. This modified methodalso solves for 3D velocity structure and static correctionssimultaneously in the inversion process. The static corrections can then be used to estimate the cover thickness. First, we test our tomography method on synthetic data withthe shot points in the bedrock below the cover. Next, we developa strategy for the case when the sources are within the cover. Themethod is then applied to field data from five crooked-line,high-resolution reflection seismic profiles ranging in lengthfrom 2 to [Formula: see text]. The crooked-line profiles make the study 2.5dimensional regarding bedrock velocities. The cover thicknessalong the profiles varies from 0 to [Formula: see text]. Estimated thickness ofthe cover agrees well with data from boreholes drilled near theprofiles. Low-velocity zones in the uppermost bedrock generallycorrelate with locations where reflections from the stackedsections project to the surface. Thus, the method is functional,both for imaging the uppermost bedrock velocities as well as for estimating the cover thickness.


Sign in / Sign up

Export Citation Format

Share Document