Group Velocity Inversion Research Articles

Crustal and uppermost mantle structure near the Gloria Fault, North Atlantic, from ocean bottom seismometer surface wave observations

SUMMARY In this work, we present both 1-D and 3-D shear wave velocity (Vs) models of the oceanic crust and uppermost mantle below the Deep OCean Test ARray area, located ∼ 70 km north of the central section of the Gloria Fault, in the eastern North-Atlantic Ocean. The velocity models are inferred from the dispersion of surface waves recorded on ocean bottom seismometers. Dispersion measurements are obtained from the analysis of ambient seismic noise at short periods (< 14 s) and teleseismic surface waves at long periods (> 14 s) using the two-station method. The 1-D Vs model is inferred from the joint inversion of Rayleigh wave phase velocities and Love wave group and phase velocities. The 3-D tomographic model is obtained by inversion of 2-D Love wave group velocity maps as a function of depth, further constrained by the average of Love wave phase velocities obtained from ambient noise (4–9 s) and the average Rayleigh and Love wave phase velocities calculated from teleseismic data (14–44 s). The 1-D Vs model shows a sediment layer with a low velocity of 1.05 km s−1, similar to previous studies in the region. Below the sediments, we find an oceanic crust with velocities ranging from 3.3 to 4.5 km s−1. The model reaches an unusually high velocity of 4.9 km s−1 in a 20 km thick layer at depths between 16 and 36 km. We interpret this fast velocity layer as indicative of the presence of harzburgite, a residue of enhanced melting that might have been formed by the proximity between the Mid-Atlantic Ridge and the Azores mantle plume. At greater depths the velocity decreases, forming a low-velocity zone that reaches a minimum at ∼ 70 km depth, which we interpret as the maximum depth for the lithosphere–asthenosphere boundary. The 3-D model shows a structure that is mostly horizontally layered, with Vs isocontours at 3.5–4.5 km s−1 highlighting oscillations of the crustal structure with wavelengths of ∼25–30 km. These oscillations may be due to changes in the rate of mantle upwelling and magma supply rate.

A New Method for Selecting the Phase and Group Velocity Dispersion Curves of Rayleigh and Love Surface Waves: Real Data Case of Central Anatolia, Turkey (Türkiye)

We propose a new method to select the dispersion curve of Rayleigh and Love surface waves obtained from local and regional earthquakes. This method is efficient for the crustal tomography studies and the Central Anatolia is chosen to test it. The single-station method utilizing both seismograms and accelerograms is used for the group velocities from local earthquakes while the two-station approach utilizing the cross correlograms is employed for the group and phase velocities using both local and regional earthquakes. The proposed method is made of two stages. The first stage based on the multiple filter technique – MFT is designed for an automatic selection of the desired dispersion curve. In case of the single-station approach, the MFT diagrams are computed for both Rayleigh and Love seismograms or accelerograms. An algorithm based on several selection criteria is employed to extract the corresponding group velocity dispersion curve from the MFT diagram. In case of the two-station method, the MFT diagrams are computed for both Rayleigh and Love cross correlograms. The latter algorithm developed for the MFT diagram is again applied to the cross correlogram to select the respective two-station group velocity dispersion curve. To select the corresponding two-station phase velocity dispersion curve, another algorithm is developed where the joint inversion of phase and group velocities is employed. The RMS value quantifying the misfit between the selected and inverted dispersion curves is found minimum when the selected phase velocity curve is compatible with the group velocity curve extracted via the MFT. The selection procedures defined up until this point constitute the first stage. In the second stage, phase and group velocity gathering around base pathways is performed. This way it is possible to compute the averages and standard deviations for the observed pathways traversing the similar geology. The outliers determined by using the averages and standard deviations are eliminated from the data set and this elimination is performed for phase and group velocities, separately.

Joint inversion of Rayleigh group and phase velocities for S-wave velocity structure of the 2021 MS6.0 Luxian earthquake source area, China

On September 16, 2021, a MS6.0 earthquake struck Luxian County, one of the shale gas blocks in the Southeastern Sichuan Basin, China. To understand the seismogenic environment and its mechanism, we inverted a fine three-dimensional S-wave velocity model from ambient noise tomography using data from a newly deployed dense seismic array around the epicenter, by extracting and jointly inverting the Rayleigh phase and group velocities in the period of 1.6–7.2 s. The results showed that the velocity model varied significantly beneath different geological units. The Yujiasi syncline is characterized by low velocity at depths of ∼ 3.0–4.0 km, corresponding to the stable sedimentary layer in the Sichuan Basin. The eastern and western branches of the Huayingshan fault belt generally exhibit high velocities in the NE-SW direction, with a few local low-velocity zones. The Luxian MS6.0 earthquake epicenter is located at the boundary between the high- and low-velocity zones, and the earthquake sequences expand eastward from the epicenter at depths of 3.0–5.0 km. Integrated with the velocity variations around the epicenter, distribution of aftershock sequences, and focal mechanism solution, it is speculated that the seismogenic mechanism of the main shock might be interpreted as the reactivation of pre-existing faults by hydraulic fracturing.

Along‐strike variation in the shear wave crustal structure of the NE Himalayan and Indo‐Burmese arc: Evidence based on surface wave dispersion analysis

Surface wave dispersion data is used to infer the shear wave velocity structure variation in the north‐east Himalaya and Indo‐Burmese arc regions. We have used 25 earthquakes data from four groups with a magnitude range of 5.0–6.7, epicentral distance range 368–800 km, and focal depth less than 50 km. Ray paths from the earthquake location to the seismic station are transversely passing different geotectonic units of the Himalayas, Indo‐Gangetic plains, and Indo‐Burma collision zones. The weighted average dispersion curve and the path averaged shear wave velocity models are computed for the four groups located at different azimuths around Shillong seismic station. Non‐linear least‐square inversion is performed to obtain the shear wave velocity structure of the crust and uppermost mantle from joint inversion of Rayleigh and Love waves group velocities. Subsequently, Rayleigh and Love waves group velocities are inverted separately to obtain vertical and horizontal components of S‐wave velocity. Next, using these components anisotropy coefficient values at different depths for each path are also estimated. A high variation of dispersion curves and the shear wave velocity models from one group to another indicates that the region is geotectonically very complex. An approximately 80 km‐thick zone beneath the study region has shear wave velocity as low as 1.7 km/s in the uppermost crust in the southern part and ~4.7 km in the uppermost mantle beneath the Eastern Himalayan Syntaxis (EHS). Inferred velocity is also lower than that of PREM and AK135 global models with a much thicker crust beneath the study region. Radial anisotropy varies even within the northern part from the Indo‐Eurasian collision zone to EHS and northern to southern Indo‐Burma ranges. Anisotropy is comparatively stronger in the deeper part below ~40 km for the three paths, except for the EHS, where the result is contrary.

The speed limit of optoelectronics

Light-field driven charge motion links semiconductor technology to electric fields with attosecond temporal control. Motivated by ultimate-speed electron-based signal processing, strong-field excitation has been identified viable for the ultrafast manipulation of a solid’s electronic properties but found to evoke perplexing post-excitation dynamics. Here, we report on single-photon-populating the conduction band of a wide-gap dielectric within approximately one femtosecond. We control the subsequent Bloch wavepacket motion with the electric field of visible light. The resulting current allows sampling optical fields and tracking charge motion driven by optical signals. Our approach utilizes a large fraction of the conduction-band bandwidth to maximize operating speed. We identify population transfer to adjacent bands and the associated group velocity inversion as the mechanism ultimately limiting how fast electric currents can be controlled in solids. Our results imply a fundamental limit for classical signal processing and suggest the feasibility of solid-state optoelectronics up to 1 PHz frequency.

Refraction of space-time wave packets: III. experiments at oblique incidence.

The refraction of space-time (ST) wave packets at planar interfaces between non-dispersive, homogeneous, isotropic dielectrics exhibits fascinating phenomena, even at normal incidence. Examples of such refractive phenomena include group-velocity invariance across the interface, anomalous refraction, and group-velocity inversion. Crucial differences emerge at oblique incidence with respect to the results established at normal incidence. For example, the group velocity of the refracted ST wave packet can be tuned simply by changing the angle of incidence. In the third paper, we present experimental verification of the refractive phenomena exhibited by ST wave packets at oblique incidence that were in the first paper of this sequence [J. Opt. Soc. Am. A38, 1409 (2021)10.1364/JOSAA.430105]. We also examine a proposal for "blind synchronization," whereby identical ST wave packets arrive simultaneously at different receivers without a priori knowledge of their locations except that they are all located at the same depth beyond an interface between two media. A first proof-of-principle experimental demonstration of this effect is provided.

Refraction of space-time wave packets: II. experiments at normal incidence.

The refraction of space-time (ST) wave packets offers many fascinating surprises with respect to conventional pulsed beams. In the first paper in this sequence [J. Opt. Soc. Am. A38, 1409 (2021)10.1364/JOSAA.430105], we theoretically described the refraction of all families of ST wave packets at normal and oblique incidence at a planar interface between two nondispersive, homogeneous, isotropic dielectrics. Here, in this second paper in the sequence, we present experimental verification of the refractive phenomena predicted for baseband ST wave packets upon normal incidence on a planar interface. Specifically, we observe group velocity invariance, normal and anomalous refraction, and group velocity inversion leading to group delay cancellation. These phenomena are verified in a set of optical materials with refractive indices ranging from 1.38 to 1.76, including MgF2, fused silica, BK7 glass, and sapphire. We also provide a geometrical representation of the physics associated with anomalous refraction in terms of the dynamics of the spectral support domain for ST wave packets on the surface of the light cone.

Refraction of space-time wave packets: I. theoretical principles.

Space-time (ST) wave packets are pulsed optical beams endowed with precise spatio-temporal structure by virtue of which they exhibit unique and useful characteristics such as propagation invariance and tunable group velocity. We study in detail here, and in two accompanying papers, the refraction of ST wave packets at planar interfaces between non-dispersive, homogeneous, and isotropic dielectrics. We formulate a law of refraction that determines the change in the ST wave-packet group velocity across such an interface as a consequence of a newly identified optical refractive invariant that we call the "spectral curvature". Because the spectral curvature vanishes in conventional optical fields where the spatial and temporal degrees of freedom are separable, these phenomena have not been observed to date. We derive the laws of refraction for baseband, X wave, and sideband ST wave packets that reveal fascinating refractive phenomena, especially for the former class of wave packets. We predict theoretically, and confirm experimentally in the accompanying papers, refractive phenomena such as group-velocity invariance (ST wave packets whose group velocity does not change across the interface), anomalous refraction (group-velocity increase in higher-index media), group-velocity inversion (change in the sign of the group velocity upon refraction but not its magnitude), and the dependence of the group velocity of the refracted ST wave packet on the angle of incidence.

Crustal model in eastern Qinghai-Tibet plateau and western Yangtze craton based on conditional variational autoencoder

Eastern margin of the Qinghai-Tibet plateau and western margin of the Yangtze craton has a complex crust-mantle structure and is an important region for studying the crust-mantle deformation mechanism. Because the Rayleigh surface wave group velocity has a strong constraint on the structure of the crust and upper mantle, a lot of researches have used the Rayleigh surface wave group velocity to invert the crust-mantle structure in this region, but the layered multi-parameter joint inversion and the uncertainty of the inversion result evaluation is still a research hotspot. Motived by Conditional Variational Autoencoder (CVAE) which is composed of a decoding and encoding process that uses a deep neural network to predict the variational distribution of parameters, we propose a crustal model by CVAE to attain the forward and inversion relationship between the group velocity of Rayleigh wave and the crustal model. In which the CVAE decoding corresponds to the forward process and the encoding corresponds to the inverse process. Moreover, the inversion results of this method are constrained from three aspects as follows: The network inputs and outputs have the same group velocities; The decoded outputs are normal distribution of hidden variables representing the crustal model; The mean values of the decoded output during CVAE training are the sampled crustal models. The proposed method can invert multi-parameters and evaluate uncertainty of the inversion efficiently. Based on the latest group velocity model, we invert for mean and variance of thickness in the upper, middle, and lower sedimentary cover, and of thickness, P-wave velocity, S-wave velocity and density in the upper, middle, and lower crystalline crust for eastern Qinghai-Tibet plateau and western Yangtze craton. Compared with the existing crustal models obtained by different methods, the results of this study are in good agreement. We conclude that CVAE-based group velocity inversion of the crustal model is a feasible and reliable method.

Group velocity maps using subspace and Trans-dimensional inversions: ambient noise tomography in the Western part of Java, Indonesia

Summary In this paper, we compare two different methods for group velocity inversion: iterative, least-squares subspace optimization, and probabilistic sampling based on the Trans-dimensional Bayesian method with tree-based wavelet parameterization. The wavelet parameterization used a hierarchical prior for wavelet coefficients which could adapt to the data. We applied these inversion methods for ambient noise tomography of the western part of Java, Indonesia. This area is an area prone to multiple geological hazards due to its proximity to the subduction of the Australia Plate beneath Eurasia. It is therefore important to have a better understanding of upper crustal structure to support seismic hazard and disaster mitigation efforts in this area. We utilized a new waveform dataset collected from 85 temporary seismometers deployed during 2016–2018. Cross-correlation of the waveform data was applied to retrieve empirical Rayleigh wave Green's functions between station pairs, and the spatial distribution of group velocity was obtained by inverting dispersion curves. Our results show that, although computationally expensive, the Trans-dimensional Bayesian approach offered important advantages over optimization, including more effective explorative of the model space and more robust characterization of the spatial pattern of Rayleigh wave group velocity. Meanwhile, the iterative, least-square subspace optimization suffered from the subjectivity of choice for reference velocity model and regularization parameter values. Our Rayleigh wave group velocity results show that for short (1–10 s) periods group velocity correlates well with surface geology, and for longer periods (13–25 s) it correlates with centers of volcanic activity.

3D joint inversion of gravity data and Rayleigh wave group velocities to resolve shear-wave velocity and density structure in the Makran subduction zone, south-east Iran

In this study, we developed a method to invert jointly Rayleigh wave group velocities and gravity anomalies for velocity and density structure of the lithosphere. We applied the method to the Makran accretionary prism, SE Iran. The reason for using different data sets is that each of these data sets is sensitive to different parameters. Surface wave group velocities are sensitive mainly to shear wave velocity distribution in depth but do not well resolve density variations. Therefore, joint inversion with gravity data increases the resolution of density distribution. Our approach differs from others mainly in the model parameterization: Instead of subdividing the model into a large number of thin layers, we invert for the properties of only four layers: thickness, P- and S-wave velocities and densities and their vertical gradients in sediments, upper-crust, lower-crust and upper mantle. The method is applied first to synthetic models in order to demonstrate its usefulness. We then applied the method to real data to investigate the lithosphere structure beneath the Makran. The resulting model shows that Moho depth increases from Oman Sea (18–33 km) and Makran fore-arc (33–37 km) to the volcanic-arc (44–46 km). The crustal density is high in the Oman Sea as should be expected for the oceanic crust. We also find a high-velocity anomaly in the upper mantle under the Oman Sea corresponding to the subducting slab. The crust under the fore-arc, volcanic-arc and back-arc settings of Makran subduction zone is characterized by low-velocity zones.

Estimation of Shallow-to-Deep Shear Wave Velocity Structure from Joint Inversion of Surface-wave Phase and Group Velocities

Estimation of Shallow-to-Deep Shear Wave Velocity Structure from Joint Inversion of Surface-wave Phase and Group Velocities

Lateral variation of crust and upper mantle structures in NW Iran derived from surface wave analysis

To obtain the shear velocity structure across North-West of Iran and surrounding areas to a depth of 160 km, we performed a namely Hedgehog nonlinear inversion on Rayleigh wave group velocity dispersion curves in the period range from 7 to 60 s. The distributed dispersion curves are the results of our surface wave dispersion tomography using the data of 280 local and regional seismic events, recorded by the medium- and broad-band seismic stations in the region. We outline different crust and upper mantle structures for the study area based on calculated group and shear velocities. Our results reveal relatively low velocities at the shorter periods (7–10 s) in the presence of sedimentary basins (e.g., South Caspian Basin) and for eastern Anatolia and relatively high velocities along the Sanandaj–Sirjan Metamorphic zone, Alborz, Talesh, and the Lesser Caucasus Mountains. By depth inversion of group velocities, we observed 14-km-thick sediments in South Caspian Basin and Kura Depression. Based on our maps at 20 s, we outline different crustal models for the region and highlight the differences between South Caspian Basin and NW Iran, on one side, and the similarities between the South Caspian Basin and Kura Depression that extend beneath Talesh, Alborz, and Lesser Caucasus, on the other. Comparing the shear velocity of lower crust in South Caspian Basin and Kura Depression with that of NW Iran proves different origination of lower crust in the basin, probably oceanic source, because of its significant higher shear velocity rather than NW Iran. In Talesh, we observe indications of an under-thrusting of the lower crust of SCB beneath NW Iran while the middle crust is locked. The analysis of group velocities at longer periods (≥ 35 s) and obtained shear velocity models allows us to outline different lithospheric structures and crustal depth in the region. The high group velocities in Talesh, South Caspian Sea, and Lesser Caucasus on one side and Zagros Folding and Thrust Belt on the other, beside the result of shear velocity models, suggest the presence of a stable and thick mantle lid that seems to be thin or absent in the eastern Anatolia and much of NW Iran. The shallowest Moho and Lithosphere Asthenosphere boundary depths of 37 and 63 km were observed in Eastern Anatolian Accretionary Complex. The thin mantle lid in this region has affected the whole crust in such a way that we observed the lowest shear velocities inside the crust in this region. We observed a significant thickening of both crust and lithosphere in Sanandaj–Sirjan Metamorphic zone comparing to Urmieh Dokhtar Magmatic Arc and Zagros Folding and Thrust Belt on its two sides.

Shear wave velocity model beneath CBJI station West Java, Indonesia from joint inversion of teleseismic receiver functions and surface wave dispersion

Earthquake signal observations around the world allow seismologists to obtain the information of internal structure of the Earth especially the Earth’s crust. In this study, we used joint inversion of receiver functions and surface wave group velocities to investigate crustal structure beneath CBJI station in West Java, Indonesia. Receiver function were calculated from earthquakes with magnitude more than 5 and at distance 30°-90°. Surface wave group velocities were calculated using frequency time analysis from earthquakes at distance of 30°- 40°. We inverted shear wave velocity model beneath the station by conducting joint inversion from receiver functions and surface wave dispersions. We suggest that the crustal thickness beneath CBJI station, West Java, Indonesia is about 35 km.

Crustal structure of Nigeria and Southern Ghana, West Africa from P-wave receiver functions

We report new estimates of crustal thickness (Moho depth), Poisson's ratio and shear-wave velocities for eleven broadband seismological stations in Nigeria and Ghana. Data used for this study came from teleseismic earthquakes recorded at epicentral distances between 30° and 95° and with moment magnitudes greater than or equal to 5.5. P-wave receiver functions were modeled using the Moho Ps arrival times, H–k stacking, and joint inversion of receiver functions and Rayleigh wave group velocities. The average crustal thickness of the stations in the Neoproterozoic basement complex of Nigeria is 36km, and 23km for the stations in the Cretaceous Benue Trough. The crustal structure of the Paleoproterozoic Birimian Terrain, and Neoproterozoic Dahomeyan Terrain and Togo Structural Unit in southern Ghana is similar, with an average Moho depth of 44km. Poisson's ratios for all the stations range from 0.24 to 0.26, indicating a bulk felsic to intermediate crustal composition. The crustal structure of the basement complex in Nigeria is similar to the average crustal structure of Neoproterozoic terrains in other parts of Africa, but the two Neoproterozoic terrains in southern Ghana have a thicker crust with a thick mafic lower crust, ranging in thickness from 12 to 17km. Both the thicker crust and thick mafic lower crustal section are consistent with many Precambrian suture zones, and thus we suggest that both features are relict from the collisional event during the formation of Gondwana.

Crustal S-wave structure beneath Eastern Black Sea Region revealed by Rayleigh-wave group velocities

In this study, the crustal S-wave structure beneath the Eastern Black Sea Region (including the Eastern Black Sea Basin (EBSB) and Eastern Pontides (EP)) has been revealed using inversion of single-station, fundamental-mode Rayleigh-wave group velocities in the period range of 4–40seconds. We used digital broadband recordings of 13 regional earthquakes that recently occurred in the easternmost EBSB recorded at stations of the Kandilli Observatory and Earthquake Research Institute (KOERI). The average group-velocity-dispersion curves were generated from 26 paths for the EBSB, and 16 paths for the EP, and they were inverted to determine the average 1-D shear-wave structure of the region. We have created a pseudo-section, roughly depicting the crustal structure of the region based on the group velocity inversion results of all station-earthquake paths. The thickness of the sedimentary layer reaches 12km in the center of EBSB (Vs=2.5–3.1km/s) and decreases 4km in the EP. There is a thin sedimentary layer in the EP (Vs=2.7km/s). A consolidated thin crust that exists in the EBSB possesses a high seismic velocity (Vs=3.8km/s). While a thin (∼26km) and transitional crust exists beneath the EBSB, a thick (about 42km) continental crust exists beneath the EP where the Conrad is clearly seen at about a 24km depth. Thick continental crust in the EP region is clearly distinguished from a gradational velocity change (Vs=3.4–3.8km/s). The Moho dips approximately southwards, and the Vs velocity (4.25–4.15km/s) beneath the Moho discontinuity decreases from the EBSB to the EP in the N–S direction. This may be an indication of a southward subduction.

Crustal velocity structure of the Deccan Volcanic Province, Indian Peninsula, from observed surface wave dispersion

<p>Through inversion of fundamental mode group velocities of Love and Rayleigh waves, we study the crustal and subcrustal structure across the central Deccan Volcanic Province (DVP), which is one of the world’s largest terrestrial flood basalts. Our analysis is based on broadband seismograms recorded at seismological station Bhopal (BHPL) in the central India from earthquakes located near west coast of India, with an average epicentral distance about 768 km. The recording station and epicentral zone are situated respectively on the northern and southern edges of DVP with wave paths across central DVP. The period of group velocity data ranges from 5 to 60 s for Rayleigh waves and 5 to 45 s for Love waves. Using the genetic algorithm, the observed data have been inverted to obtain the crust and subcrustal velocity structure along the wavepaths. Using this procedure, a similar velocity structure was also obtained earlier for the northwestern DVP, which is in the west of the present study region. Comparison of results show that the crustal thickness decreases westward from central DVP (39.6 km) to northwestern DVP (37.8 km) along with the decrease of thickness of upper crust; while the thickness of lower crust remains nearly same. From east to west S-wave velocity in the upper crust decreases by 2 to 3 per cent, while P-wave velocity in the whole crust and subcrust decreases by 3 to 6 per cent. The P- and S-wave velocities are positively correlated with crustal thickness and negatively correlated with earth’s heat flow. It appears that the elevated crustal and subcrustal temperature in the western side is the main factor for low velocities on this side.</p>

Earth structure and instrumental seismicity of Madagascar: Implications on the seismotectonics

The aim of this study is to improve the knowledge of the seismotectonics of Madagascar. We first investigate the structure of the Earth beneath Madagascar through the joint inversion of receiver functions and Rayleigh wave group velocities. Then we use the obtained velocity models to relocate local earthquakes in order to analyse the distribution of seismicity. Finally, we use structural models and earthquake coordinates to compute focal mechanisms.Our retrieved Earth structure models confirm a thin lithosphere beneath Madagascar when compared to the nearby East African Rift. The High Plateau in the Central region coincides with the thinnest lithosphere over the slowest asthenosphere. Our results are in good agreement with the gravity anomalies and likely confirm a localised asthenospheric upwelling beneath the central part of Madagascar. The surface expression of the asthenospheric upwelling consists in a horst–graben structure. The moderate seismicity is localised along pre-existing structures reflecting an E-W extension that is mostly accommodated in the lower crust.

Crustal structure of the Khartoum Basin, Sudan

The crustal structure of the northern part of the Khartoum Basin has been investigated using data from 3 permanent seismic stations within 40km of Khartoum and two modeling methods, H–k stacking of receiver functions and a joint inversion of receiver functions and Rayleigh wave group velocities. The Khartoum Basin is one of several Mesozoic rift basins in Sudan associated with the Central African Rift System. Results from the H–k-stacking indicate that crustal thickness beneath the Khartoum Basin ranges between 33 and 37km, with an average of 35km, and that the crustal Vp/Vs ratio ranges from 1.74 to 1.81, with an average of 1.78. From the joint inversion of receiver functions and Rayleigh wave group velocities, we obtained similar results for Moho depth, as well as an average shear wave velocity of 3.7km/s for the crust. These results provide the first seismic estimates of Moho depth for a basin in Sudan. When compared to average crustal thickness for unrifted Proterozoic crust in eastern Africa, our results indicate that at most only a few km of crustal thinning may have occurred beneath the Khartoum Basin. This finding is consistent with estimates of effective elastic plate thickness, which indicate little modification of the Proterozoic lithosphere beneath the basin, and suggests that there may be insufficient topography on the lithosphere–asthenosphere boundary beneath the Sudanese basins to channel plume material westward from Ethiopia.

S-wave velocity structure of northeastern China from joint inversion of Rayleigh wave phase and group velocities

SUMMARY We imaged the crust and uppermost mantle structure beneath northeastern (NE) China with fundamental mode Rayleigh waves recorded by 125 broad-band stations deployed in the region. Rayleigh wave phase and group velocities along more than 700 interstation paths were estimated using the wavelet transformation method and then these data were utilized to construct 2-D phase and group velocity maps in the period range of 15–60 s. Owing to the dense ray coverage and short ray path, our results provide better lateral resolution in the NE China region compared with previous phase and group velocity studies. The regularized Rayleigh wave phase and group velocity dispersions at each cell were jointly inverted to determine 1-D shear wave velocity structure using the linear inversion method and then assembled into a 3-D model. Our results show that obvious low velocities exist in the uppermost mantle beneath the Changbaishan volcanic region, which may be due to asthenospheric upwelling. The thin lithosphere with fast S-wave velocity in the lower crust of the Songliao Basin implies that the lithospheric mantle beneath NE China is partly removed.