Seismic Discontinuity Research Articles

Double-diffusive layering during the evolution of planetary mantles

Summary The degree of layering during planetary mantle evolution has been a key issue and is still heavily debated, especially since the observation of seismic discontinuities within the Earth mantle and their correlation to phase transitions within the mantle mineral assemblage and the detection of geochemical differences in mantle derived basalts (Mid Ocean Ridge Basalts and Ocean Island Basalts). We present the phenomenon of double-diffusive layering in a fluid with strongly temperature dependent viscosity, as relevant for planetary mantle material, to study the influence of an initial compositional gradient with regard to the thermal and chemical evolution of a planet. The numerical experiments show that in a wide parameter range distinct layers are formed self-organized from a continuously stratified state by dynamical fractionation and are thus likely to appear as a generic. Considering this as a plausible model for planetary mantle evolution it provides a dynamical explanation for the existence of distinct chemical reservoirs within the history of a planet's mantle.

Imaging of seismic discontinuities using an adjoint method

Summary For imaging of seismic discontinuities at depth, reverse time migration (RTM) is a powerful method to apply to recordings of seismic events. It is especially powerful when an extensive receiver array, numerous seismic sources, or both, permit adequate reconstruction of incident and scattered wavefields at depth. Reconstructing either the incident or scattered wavefield at depth becomes less accurate when relatively few recordings of seismic events are available. Here we explore an inverse scattering approach to imaging discontinuities based on an adjoint method, employing sensitivity kernels (Fréchet derivatives) that represent jumps in material properties across seismic-discontinuity surfaces. When combined with ray-based requirements on scattering geometry, it constitutes a powerful approach to determining the locations and amplitudes of the discontinuities, recovering only those properties that can be resolved by a spatially limited source and/or receiver distribution. This is illustrated by synthetic examples with local sources followed by a field example in a subduction zone setting.

A shallow mantle seismic discontinuity beneath northeast India: Evidence from receiver function analyses

The crust and shallow upper mantle structure beneath the Upper Brahmaputra Valley, Indo-Burma Ranges, and Bengal Basin of Northeast India have been investigated based on receiver function (RF) analysis of teleseismic earthquakes recorded by 11 seismological stations. The study reveals a thin crust (∼35 km) beneath the Brahmaputra Valley (at JORH station) with a surface sedimentary layer of ∼4 km thick. The crustal thickness is observed to increase towards the north in the Himalaya (∼40 km at ZIRO and ITAN) and to the south (up to ∼46 km at KOHI). The crustal thickness near the Tripura fold-belt and Bengal Basin varies within ∼36–40 km. The study reveals the existence of a shallow mantle discontinuity (Hales discontinuity) at a variable depth range of ∼54–78 km characterized by a step increase (∼7.5–11 %) in shear wave velocity observed in the inverted models. The mineralogical phase transformation from spinel to garnet is considered as the origin of this discontinuity. The shallow depth of the discontinuity indicates an increase in upper mantle temperature which conforms to the high geothermal gradient reported in the region. The variation of depth of the discontinuity can be interpreted in terms of the addition of Cr+3 that shifts the spinel-garnet stability field to higher depths whereas Fe+2 shifts it to lower depths. Despite the high temperature in the upper mantle, the observed low Vp/Vs ratio (1.65–1.75) below the Hales discontinuity can be explained by the presence of a high fraction of orthopyroxene.

Seismic discontinuity in the Martian crust possibly caused by water-filled cracks

Recent seismic data acquired by the InSight lander have revealed seismic discontinuities in the Martian crust that have been interpreted as sharp transitions in porosity or chemical composition. Here we propose an alternative model in which the transition from dry cracks to water-filled cracks could explain the observed seismic discontinuity in the Martian crust. Our model can explain sharp increases in seismic velocity and Vp/Vs at ∼10 km depth with no associated changes in porosity or chemical composition. The present model suggests the local existence of liquid water in the Martian crust, which could potentially serve as a subsurface habitat for life.

Topography of the 410 and 660 km discontinuities beneath the Tibetan plateau and adjacent areas

To investigate the three-dimensional spatial distribution of subducted oceanic slab segments and their consequential effect on the thermal and hydrous composition of the mantle transition zone (MTZ) beneath the entire Tibetan Plateau, we conducted a comprehensive analysis of approximately 655,000 P-to-s receiver functions (RFs), obtained from 735 broadband seismic stations. These RFs were utilized to delineate the 410 km (d410) and 660 km (d660) seismic discontinuities, which represent the upper and lower boundaries of the MTZ, respectively. The RFs were grouped within circular bins with a radius of 1 degree and stacked to image the discontinuities. The mean apparent depths derived based on the 1-D IASP91 Earth model for the d410 and d660 across the entire study area are 412.2±8.3 km and 668.9±8.5 km respectively, and the MTZ thickness is 256.5±6.9 km. The observed apparent depths underwent subsequent correction utilizing multiple velocity models. Several areas in central Tibet exhibit a normal d410 and an anomalously deep d660, which can be attributed to the combined effect of the negative thermal anomaly and dehydration associated with subducted slab segments that have penetrated at least to the d660 depth. The anomalous thickening of the MTZ beneath the southeastern Tibetan Plateau surrounding the Tengchong volcanic field can be explained by the dehydration of the subducted Indian Slab. Significant thinning of the MTZ associated with the deepening of the d410 beneath the western Tian Shan may indicate active thermal upwelling originating from the MTZ.

Assessing foundation characteristics at the war dam site, lake tana basin, Ethiopia: A geophysical and geotechnical perspective

An integrated geophysical and geotechnical study evaluated the foundation conditions at the War dam site in northwest Ethiopia. This investigation included the classification of rock quality, shallow seismic refraction, and magnetic approaches. The dam's location comprises quaternary soil deposits and rhyolite rock units that have undergone varied weathering and fracturing. The shallow seismic refraction method distinguishes three layers of p-wave velocities that are less than 1.5 km per second with a depth range of 2–6 m, 1.5–2.5 km per second at a depth range of 15–20 m, and 2.5–3.5 km per second ranging from 20 to 40 m, respectively. Magnetic data were used to identify lineaments, and the RQD value acquired from boreholes ranged from extremely poor to excellent. Lineaments were recognized using the tilt angle approach. The results of the permeability tests demonstrated that the rock mass that serves as the dam's foundation had characteristics that are resistant to low permeability. The maximum and minimum lugeon values obtained from the testing were 9Lu and 0.81Lu, respectively. There are weak zones at and below the surface of the dam site, according to the overall findings acquired from seismic refraction, magnetic, and discontinuity surveying. These results were obtained from monitoring the dam site. These significant structures are directed towards a SW-NE, NE-SW, NNW-SSE, and SSW-NNE orientation. The study assessed the geological suitability of a proposed dam site using seismic refraction and magnetic survey methods. Significant geological variations were observed, particularly in the right abutment and valley floor, indicating the need for targeted grouting. The findings suggest that while the site is generally suitable for dam construction, specific areas require further ground improvement to ensure stability.

Stabilizing inverse ringwoodite with defects, and a possible origin for the 560-km seismic discontinuity

Ringwoodite is an important mineral in the mantle transition zone, and its cationic disorder can profoundly affect its physicochemical properties, but there is currently much controversy about this disorder. In this study, we investigate the cation disorder states of pure Mg2SiO4-ringwoodite and defective ringwoodite under mantle transition zone conditions through DFT calculations and thermodynamic models. Two stable endmembers are seen, one with normal ringwoodite structure and the other with inverted structure (its Si atoms and half of its Mg atoms have swapped sites). Our results indicate that pure ringwoodite does not invert (swap Mg and Si cations) under normal mantle temperatures but the introduction of a Si-excess, Mg-deficient defect induces a swap at normal mantle temperatures and this swap is likely induced by a wide range of defects including water. Thus, in the presence of such a defect or similar defects the olivine phase transition sequence may then go from olivine to wadsleyite to inverse ringwoodite, and then normal ringwoodite. We calculate the seismic properties of normal and inverse ringwoodite and find significantly slower wave speeds in inverted ringwoodite. Due to this difference the presence of inverse ringwoodite may provide a potential explanation for the discontinuous interface of seismic waves at the depth of ∼560 km.

Evidence of a Low‐Velocity Zone in the Upper Mantle Beneath Cumbre Vieja Volcano (Canary Islands) Through Receiver Functions Analysis

AbstractOceanic volcanic islands have a complex internal structure due to their geological evolution. La Palma is the second youngest in the Canary Islands. Historical eruptions occurred in the Cumbre Vieja (CV) volcanic system, including the last 2021 Tajogaite eruption. The seismicity before, during, and after the Tajogaite eruption provides valuable information for a better understanding of deep volcanic processes. We applied the receiver function technique to image the crustal and upper mantle structures beneath the CV volcanic complex. We investigated the first 50 km of the lithosphere, identifying various seismic discontinuities. After a thin layer of recent volcanic rocks, we found a layer corresponding to the old oceanic crust. Beneath the Moho, we found a significant decrease in the seismic velocities until about 37 km of depth. This low‐velocity layer corresponds to a zone of partial melting with interconnected magmatic chambers that feed the volcanic activity in CV.

Unveiling the Distinct Structure of the Upper Mantle Beneath the Canary and Madeira Hotspots, as Depicted by the 660, 410, and X Discontinuities

AbstractThe Canary and Madeira Islands are two distinct hotspots in the Central‐East Atlantic that are close to each other. Their volcanism is generally attributed to underlying mantle plumes, but the detailed structure of these plumes is still not well understood. The thermal and compositional structure of the plume introduces complexities in the phase transitions of the mantle, which impact the depth and magnitude of seismic discontinuities. We use 1,268 high‐quality receiver functions from stations located at the two hotspots to detect P‐to‐s converted phases through a common‐conversion point stacking approach and conduct a detailed analysis of mantle seismic discontinuities. The results show that both hotspots are characterized by a thin mantle transition zone (MTZ), with sharp 410 and 660 discontinuities at depths of 428–421 km and 647–664 km, beneath the Canaries and Madeira respectively. The results indicate that the Canary plume crosses the MTZ, whereas the Madeira plume mainly influences the upper portion of the MTZ. Furthermore, we find reliable detections of a sharp X discontinuity beneath the Canaries at 287 km. Its presence suggests the accumulation of silica‐rich recycled eclogite at these depths. We also use the amplitudes of P410s and PXs to derive velocity jumps at corresponding discontinuities. Based on these measurements, we estimate that the basalt proportion is 60%–80%, with accumulation being more significant in the Canaries than in Madeira. The MTZ thickness, the presence of the X discontinuity, and the high basalt proportion provide compelling evidence for a deep‐rooted thermochemical plume beneath the study area.

A global SS precursor method for imaging discontinuities: the Moho and beyond

SUMMARY Imaging seismic velocity discontinuities within the Earth's interior offers important insight into our understanding of the tectonic plate, associated mantle dynamics, and the evolution of the planet. However, imaging velocity discontinuities in locations where station coverage is sparse, is sometimes challenging. Here we demonstrate the effectiveness of a new imaging approach using deconvolved SS precursor phases. We demonstrate its effectiveness by applying it to synthetic seismograms. We also apply it to ∼1.6 M SS precursor waveforms from the global seismic database (1990–2018) for comparison with CRUST1.0. We migrate to depth and stack the data in circular 6° bins. The synthetic tests demonstrate that we can recover Moho depths as shallow as 20 km. Globally, the Moho is resolved at 21–67 km depth beneath continental regions. The Moho increases in depth from 21 km ± 4 km beneath the continental shelf to 45–50 km beneath the continental interiors and is as deep as 67 ± 4 km beneath Tibet. We resolve the Moho in 77 percent of all continental bins, within 10 km of CRUST1.0, with all outliers located in coastal regions. We also demonstrate the feasibility of using this method to image discontinuities associated with the mantle transition zone with both synthetic and real data. Overall, the approach shows broad promise for imaging mantle discontinuities.

Ridgelet Transform Based on Optimal Basic Wavelet and Its Application in Seismic Discontinuity Detection

Ridgelet Transform Based on Optimal Basic Wavelet and Its Application in Seismic Discontinuity Detection

Seismic attribute-assisted seismic fault interpretation

Fault mapping is one of the main tasks of 3D seismic interpretation, and seismic discontinuity attributes are often used to support fault mapping in vertical sections and time slices. The most commonly used fault mapping procedure involves three passes/generations. First, the approximate positions of faults (generation I) on some vertical sections are interpreted manually. Second, the generation I faults are projected on some time slices, and the individual fault traces (generation II) are picked on time slices. Finally, the generation II faults are projected onto the vertical sections and the fault sticks (generation III) are picked in the vertical sections. To speed up the fault mapping process, we use the seismic discontinuity attribute as input to automatically extract the generation III faults in the vertical sections under the generation I and II fault mapping conditions. A quality control (QC) step is essential to improve the accuracy of the automatically generated fault sticks. To reduce the time required for QC, we propose to output the fault sticks on the seismic vertical sections defined by the interpreters. The workflow is applied to a real seismic survey to demonstrate its effectiveness.

Structure of the crust-uppermost mantle beneath the Ethiopian volcanic province using ambient seismic noise and teleseismic P wave coda autocorrelation

The Ethiopian-Yemen plateaus have been affected by plume-related magmatism for at least 40 My, yet the nature and distribution of intrusive magmas and their relation to subsequent East Africa rifting remain poorly characterized. This paper presents a 3D shear (S) - wave velocity structure of the crust and uppermost mantle beneath the northern sectors of the East African rift and the dynamically-supported Ethiopian Plateau. We apply a Bayesian inversion of group and phase velocity dispersion curves for periods 5–40 s extracted from the ambient noise tomography of data from over 125 broadband seismic stations, and we identify seismic discontinuities using autocorrelation of P-wave coda. The resulting crustal structure affecting much of the Main Ethiopian Rift (MER) indicates a pronounced low S-wave velocity lower crust (∼3.6 km/s) that becomes thinner towards the northeast rising to almost the same depth level (∼20 km) as the base crust beneath the Afar triple junction zone. We also find high velocity features in the crust and uppermost mantle beneath the Afar rift axes, the current locus of strain and volcanism, that may indicate voluminous new basaltic crust. Our Moho-depth map indicates that the entire Afar Depression is underlain by crust thinner than 25 km, whereas the MER and central Ethiopian plateau are underlain by low S-wave velocity crust thicker than 35 km. These patterns suggest active underplating beneath the largely unfaulted plateau, in contrast to the efficient magmatic feeder system supplying a narrow zone of incipient seafloor spreading in the Afar depression. Within the uppermost mantle we report evidence of a localized asthenospheric upwelling beneath the western part of the plateau southeast of Lake Tana, overlying a cooled magma body within the upper crust. A broader low-velocity uppermost mantle underlies the crust beneath the Western Plateau. These observations attest to active mantle upwelling, heating, and melt intrusion at relatively different length scales beneath the plateau in agreement with the ongoing active strains reported away from the rift zones. At 45–60 km depth beneath the rift and its shoulders, we observe a fast lid (4.2–4.4 km/s). Beneath the plateau, upper mantle velocity highs that are at least 5 km below the Moho are observed. These areas have most likely experienced prolonged magma extraction and may mark highly depleted areas.

Extraction of Mantle Discontinuities From Teleseismic Body‐Wave Microseisms

AbstractOcean swell activities excite body‐wave microseisms that contain information on the Earth's internal structure. Although seismic interferometry is feasible for exploring structures, it faces the problem of spurious phases stemming from an inhomogeneous source distribution. This paper proposes a new method for inferring seismic discontinuity structures beneath receivers using body‐wave microseisms. This method considers the excitation sources of body‐wave microseisms to be spatially localized and persistent over time. To detect the P‐s conversion beneath the receivers, we generalize the receiver function analysis for earthquakes to body‐wave microseisms. The resultant receiver functions are migrated to the depth section. The detected 410‐ and 660‐km mantle discontinuities are consistent with the results obtained using earthquakes, thereby demonstrating the feasibility of our method for exploring deep‐earth interiors. This study is a significant step toward body‐wave exploration considering the sources of P‐wave microseisms to be isolated events.

Variation of b-Value of Normal, Reverse and Strike-Slip Faulting Earthquakes with Focal Depth

A crucial factor in evaluating seismic activity, seismotectonics, and seismic hazard assessment is the distribution of earthquake sizes (b-value). The goal of this research is to explore the relation between depth and the b-constant values of normal, inverse, and strike-slip faulting earthquakes. The earthquake catalog used in this work was extracted from the Global Centroid Moment Tensor Catalog (GCMTC) and covers the period from 2008 to 2017. The focal depth of the selected tremors ranges from zero to 700 km. The magnitude of completeness for normal and reverse faulting events is 5.2 and 5.1 for strike-slip faulting events. The value of the b- constant was estimated using the maximum likelihood technique. The findings of regression analysis showed that there is an insignificant negative poor correlation relation between the b- value of the normal and strike-slip faulting events and depth, while that relation is an insignificant positive poor correlation for reverse faulting. The depth of the turning point for the b-constant value of normal and reverse faulting earthquakes is near the mantle seismic discontinuities 410 and 520 km. The depth of inflection points of b-constant value strike-slip faulting earthquakes is near to lithosphere-asthenosphere boundary (LAB) at depth of 250km and the mantle seismic discontinuity of 520km.

Seismic Architecture of the Lithosphere-Asthenosphere System in the Western United States from a Joint Inversion of Body- and Surface-wave Observations: Distribution of Partial Melt in the Upper Mantle

Quantitative evaluation of the physical state of the upper mantle, including mapping temperature variations and the possible distribution of partial melt, requires accurately characterizing absolute seismic velocities near seismic discontinuities. We present a joint inversion for absolute but discontinuous models of shear-wave velocity (Vs) using 4 types of data: Rayleigh wave phases velocities, P-to-s receiver functions, S-to-p receiver functions, and Pn velocities. Application to the western United States clarifies where upper mantle discontinuities are lithosphere-asthenosphere boundaries (LAB) or mid-lithospheric discontinuities (MLD). Values of Vs below 4 km/s are observed below the LAB over much of the Basin and Range and below the edges of the Colorado Plateau; the current generation of experimentally based models for shear-wave velocity in the mantle cannot explain such low Vs without invoking the presence of melt. Large gradients of Vs below the LAB also require a gradient in melt-fraction. Nearly all volcanism of Pleistocene or younger age occurred where we infer the presence of melt below the LAB. Only the ultrapotassic Leucite Hills in the Wyoming Craton lie above an MLD. Here, the seismic constraints allow for the melting of phlogopite below the MLD.

Geophysical Subsoil Characterization and Modeling Using Cluster Analysis for Seismic Microzonation Purposes

In the municipality of Enna, 80 HVSR measurements were performed, and some of these were combined with MASW seismic measurements, which made it possible to constrain the data inversion and obtain significant shear wave velocity models. A reconstruction of the depth of the seismic bedrock was performed for the whole territory, showing different depths for the higher and lower areas, as evidenced also by the Vseq parameter map. The frequency peaks identified in the H/V curve were analyzed through a cluster analysis algorithm to evaluate similarities that allow these peaks to be divided according to their stratigraphic origin. A non-hierarchical analysis algorithm modified in such a way as to avoid any a priori choice that could influence the partition has been used. The cluster analysis made it possible to divide the frequency peaks into five groupings, each of which was then associated with a seismic discontinuity, according to the geological contacts expected in the subsoil. Finally, the inversion of the data made it possible to reconstruct the geometries of these geological contact surfaces and to reconstruct a 3D model of the subsoil, which agrees well with the surface geology of the area.

Converted-wave reverse time migration imaging in subduction zone settings

SUMMARYWe use a newly developed 2-D elastic reverse time migration (RTM) imaging algorithm based on the Helmholtz decomposition to test approaches for imaging the descending slab in subduction zone regions using local earthquake sources. Our elastic RTM method is designed to reconstruct incident and scattered wavefields at depth, isolate constituent P- and S-wave components via Helmholtz decomposition, and evaluate normalized imaging functions that leverage dominant P and S signals. This method allows us to target particular converted-wave scattering geometries, for example incident S to scattered P, which may be expected to have dominant signals in any given data set. The method is intended to be applied to dense seismic array observations that adequately capture both incident and converted wavefields. We draw a direct connection between our imaging functions and the first-order contrasts in shear wave material properties across seismic discontinuities. Through tests on synthetic data using either S → P or P → S conversions, we find that our technique can successfully recover the structure of a subducting slab using data from a dense wide-angle array of surface stations. We also calculate images with a small-aperture array to test the impact of array geometry on image resolution and interpretability. Our results show that our imaging technique is capable of imaging multiple seismic discontinuities at depth, even with a small number of earthquakes, but that limitations arise when a small aperture array is considered. In this case, the presence of artefacts makes it more difficult to determine the location of seismic discontinuities.

Earth's evolving geodynamic regime recorded by titanium isotopes.

Earth's mantle has a two-layered structure, with the upper and lower mantle domains separated by a seismic discontinuity at about 660 km (refs. 1,2). The extent of mass transfer between these mantle domains throughout Earth's history is, however, poorly understood. Continental crust extraction results in Ti-stable isotopic fractionation, producing isotopically light melting residues3-7. Mantle recycling of these components can impart Ti isotope variability that is trackable in deep time. We report ultrahigh-precision 49Ti/47Ti ratios for chondrites, ancient terrestrial mantle-derived lavas ranging from 3.8 to 2.0 billion years ago (Ga) and modern ocean island basalts (OIBs). Our new Ti bulk silicate Earth (BSE) estimate based on chondrites is 0.052 ± 0.006‰ heavier than the modern upper mantle sampled by normal mid-ocean ridge basalts (N-MORBs). The 49Ti/47Ti ratio of Earth's upper mantle was chondritic before 3.5 Ga and evolved to a N-MORB-like composition between approximately 3.5 and 2.7 Ga, establishing that more continental crust was extracted during this epoch. The +0.052 ± 0.006‰ offset between BSE and N-MORBs requires that <30% of Earth's mantle equilibrated with recycled crustal material, implying limited mass exchange between the upper and lower mantle and, therefore, preservation of a primordial lower-mantle reservoir for most of Earth's geologic history. Modern OIBs record variable 49Ti/47Ti ratios ranging from chondritic to N-MORBs compositions, indicating continuing disruption of Earth's primordial mantle. Thus, modern-style plate tectonics with high mass transfer between the upper and lower mantle only represents a recent feature of Earth's history.

Regional stratigraphically discordant dolomite mapping based on a dolomitization mechanism on the Arabian carbonate shelf

Dolomite mapping, as a first step of carbonate diagenetic modeling, is critical to understanding dolomitization mechanisms and predicting reservoir quality. Stratigraphically discordant dolomite bodies within the Upper Jurassic intervals have long been studied on the Arabian shelf. Our dolomitization mechanisms find that faults/fractures serve as migration conduits and determine dolomite distribution by controlling fluid flow. We establish an innovative approach based on the diagenetic mechanisms and comprehensive characterization of fracture systems at multiple scales to delineate dolomite in the study area. The complex fracture networks in the subsurface are categorized into macroscale, mesoscale, and microscale fractures. Macroscale fractures, i.e., faults, are much greater than the seismic wavelength and are easily observed in seismic sections due to their obvious seismic discontinuities. Mesoscale fractures are slightly greater than or equal to the seismic wavelength and are recognized using seismic attributes. Microscale fractures are significantly smaller than the seismic wavelength and are observed mainly on core samples and thin sections. The dolomite mapping workflow includes three steps: (1) calibrate the seismic response of multiscale fractures with borehole image logs and core interpretations, (2) implement multiscale fracture characterization using multiple seismic attributes, and (3) interpret dolomite bodies based on fracture characterization and log identification. The mapping results show that massive dolomite bodies are heterogeneously developed and distributed mainly in the northeastern part of the study area, with a southward decrease in dolomite content. The application demonstrates that multiscale fracture systems play critical roles in massive dolomitization, and multiscale fracture prediction provides an innovative solution to delineate complex dolomitization systems and the combination of reservoir and seal in diagenetic traps.