S-wave Velocity Models Research Articles

Sedimentary structure of the Sichuan Basin derived from seismic ambient noise tomography

SUMMARY Determining a detailed 3-D velocity model with high resolution for the sedimentary layer in the Sichuan Basin is potentially beneficial both to the industrial oil/gas exploration and earthquake hazards’ mitigation. In this study, we apply the ambient noise tomography method to construct a 3-D S-wave velocity model. This model focuses on the sedimentary layer of the Sichuan Basin, with a 0.3° × 0.3° grid precision. Dispersion curves of both group and phase velocities of Rayleigh wave at 4 to 40 s periods are utilized, which are extracted from 87 broad-band stations in the Sichuan Basin and the surrounding areas. The 3-D model reveals a thick sedimentary layer of the Sichuan Basin with S-wave velocity ranging from ∼2.0 to 3.4 km s−1. The sediment thickness in the margins of the Sichuan Basin is generally greater than the typical values of 6–10 km in the central areas due to surrounding orogenic activities, with a maximum depth of ∼13 km in the northwestern margin. Moreover, a prominent low S-wave velocity anomaly in the margins may be caused by the sediment accumulations from large-scale landslides and pronounced denudation of the surrounding orogenic belts. Major geologic units in the sedimentary layer are delineated in this study. The S-wave velocity values within each geologic unit and their bottom interfaces are obtained. Based on our model, we calculate synthetic ground motions for the 2013 Lushan earthquake and obtain the distribution of the peak ground acceleration from the earthquake epicentre to the western Sichuan Basin. The result clearly illustrates the basin amplification effect on the seismic waves.

Crustal seismogenic structures and deformation styles along the Longmen Shan Fault belt in the eastern Tibetan Plateau inferred from ambient noise tomography

A seismic array comprising 80 broadband stations with ~10–20 km inter-station distances was deployed along the Longmen Shan fault belt (LMSF), the eastern boundary of the Tibetan Plateau. The recorded ambient noise data provided densely distributed inter-station cross-correlated surface waves. A new 3-D crustal S-wave velocity model for the LMSF was constructed by carrying out ambient noise tomography. The inverted model strongly improved data fitting and decreased data misfit compared to the reference (initial) model. The model highlights several crustal structure features. The Baoxing and Pengguan Massifs on the mountain side of the southern-to-middle LMSF exhibit relatively high crustal velocities, probably indicating strong crust. Low crustal velocities that may reflect weak, deformable brittle crust, exist mainly beneath the middle-to-northern segment of the LMSF and partly around the periphery of the Baoxing and Pengguan Massifs in its southern-to-middle segment. Two SW dipping low-velocity (weak) belts approximately perpendicular to the LMSF are imaged respectively around the Wenchuan-earthquake hypocenter in the south and Beichuan in the north. The low velocities in the two belts may focus movement of the eastern Tibetan Plateau relative to the Sichuan Basin (Yangtze Craton), and the uplift of the Tibetan Plateau over long time periods. Based on the velocity and tectonic structures, the young, high topography and thickened crust but low GPS shortening rates around the southern-to-middle LMSF may be due to the dominant effect of vertical crustal deformation caused by the existence of the strong Baoxing and Pengguan Massifs. This would then be in contrast to the characteristic lateral movements due to ductile crustal flow or weak, deformable brittle crust typical of the middle-to-northern LMSF.

Retrieval of shallow S-wave profiles from seismic reflection surveying and traffic-induced noise

In urban subsurface exploration, seismic surveys are mostly conducted along roads where seismic vibrators can be extensively used to generate strong seismic energy due to economic and environmental constraints. Generally, Rayleigh waves also are excited by the compressional wave profiling process. Shear-wave (S-wave) velocities can be inferred using Rayleigh waves to complement near-surface characterization. Most vibrators cannot excite seismic energy at lower frequencies (<5 Hz) to map greater depths during surface-wave analysis in areas with low S-wave velocities, but low-frequency surface waves ([Formula: see text]) can be extracted from traffic-induced noise, which can be easily obtained at marginal additional cost. We have implemented synthetic tests to evaluate the velocity deviation caused by offline sources, finding a reasonably small relative bias of surface-wave dispersion curves due to vehicle sources on roads. Using a 2D reflection survey and traffic-induced noise from the central North China Plain, we apply seismic interferometry to a series of 10.0 s segments of passive data. Then, each segment is selectively stacked on the acausal-to-causal ratio of the mean signal-to-noise ratio to generate virtual shot gathers with better dispersion energy images. We next use the dispersion curves derived by combining controlled source surveying with vehicle noise to retrieve the shallow S-wave velocity structure. A maximum exploration depth of 90 m is achieved, and the inverted S-wave profile and interval S-wave velocity model obtained from reflection processing appear consistent. The data set demonstrates that using surface waves derived from seismic reflection surveying and traffic-induced noise provides an efficient supplementary technique for delineating shallow structures in areas featuring thick Quaternary overburden. Additionally, the field test indicates that traffic noise can be created using vehicles or vibrators to capture surface waves within a reliable frequency band of 2–25 Hz if no vehicles are moving along the survey line.

Converted-wave model building and imaging based on common-focus-point methodology

Seismic images can be viewed as photographs for underground rocks. These images can be generated from different reflections of elastic waves with different rock properties. Although the dominant seismic data processing is still based on the acoustic wave assumption, elastic wave processing and imaging have become increasingly popular in recent years. A major challenge in elastic wave processing is shear-wave (S-wave) velocity model building. For this reason, we have developed a sequence of procedures for estimating seismic S-wave velocities and the subsequent generation of seismic images using converted waves. We have two main essential new supporting techniques. The first technique is the decoupling of the S-wave information by generating common-focus-point gathers via application of the compressional-wave (P-wave) velocity on the converted seismic data. The second technique is to assume one common VP/ VS ratio to approximate two types of ratios, namely, the ratio of the average earth layer velocity and the ratio of the stacking velocity. The benefit is that we reduce two unknown ratios into one, so it can be easily scanned and picked in practice. The PS-wave images produced by this technology could be aligned with the PP-wave images such that both can be produced in the same coordinate system. The registration between the PP and PS images provides cross-validation of the migrated structures and a better estimation of underground rock and fluid properties. The S-wave velocity, computed from the picked optimal ratio, can be used not only for generating the PS-wave images, but also to ensure well registration between the converted-wave and P-wave images.

RomUkrSeis: Seismic model of the crust and upper mantle across the Eastern Carpathians – From the Apuseni Mountains to the Ukrainian Shield

RomUkrSeis is a controlled source wide-angle reflection and refraction (WARR) profile acquired in August 2014. It is 675 km long, running roughly SW-NE from the Apuseni Mountains in Romania and the Transylvanian Basin, crossing the arc of the Eastern Carpathian orogen and terminating in the East European Craton (EEC) in SW Ukraine. Well-constrained 2-D ray-tracing P- and partly S-wave velocity models have been constructed along the profile from 348 single-component seismic recorders and eleven shot points. The Eastern Carpathian arc formed in the Cenozoic and have obscured the pre-existing Teisseyre-Tornquist Zone (TTZ), which is a transition zone between the Precambrian EEC and continental terranes accreted to it from the southwest in the Palaeozoic. The TTZ is characterised by low-velocity through its entire crust (6.0–6.3 km/s) and a considerable width (~140 km). It is interpreted as EEC crust stretched during rifting and continental margin formation in the Neoproterozoic and early Palaeozoic. The crust of the TTZ has a “trough in trough” structure wherein an upper body of ~40 km width comprising Outer Carpathian (Vp 4.9 km/s) and Late Palaeozoic-Mesozoic (Vp 5.4 km/s) units to 15 km depth lies above a wider, deeper one of inferred Neoproterozoic-early Palaeozoic strata. The crust of the Transylvanian Basin and Apuseni Mountains is relatively thin (~32 km). A high-velocity body at 4–12 km depth in this area is interpreted as a rootless fragment of an ophiolite complex exposed at the surface in this area. The lower crust beneath the Transylvanian Basin displays higher velocities than adjacent segments. Moho topography is strongly differentiated along the profile, varying from 32 to 50 km. The Moho shape, especially in the area between the Inner and Outer Carpathians, suggests a NE dip and, hence, thrusting of the Tisza-Dacia lowermost crustal and upper mantle units under the TTZ domain which, in turn, could be thrust under the cratonic (EEC) block.

Lithospheric structures of and tectonic implications for the central–east Tibetan plateau inferred from joint tomography of receiver functions and surface waves

SUMMARYWe present an updated joint tomographic method to simultaneously invert receiver function waveforms and surface wave dispersions for a 3-D S-wave velocity (Vs) model. By applying this method to observations from ∼900 seismic stations and with a priori Moho constraints from previous studies, we construct a 3-D lithospheric S-wave velocity model and crustal-thickness map for the central–east Tibetan plateau. Data misfit/fitting shows that the inverted model can fit the receiver functions and surface wave dispersions reasonably well, and checkerboard tests show the model can retrieve major structural information. The results highlight several features. Within the plateau crustal thickness is >60 km and outwith the plateau it is ∼40 km. Obvious Moho offsets and lateral variations of crustal velocities exist beneath the eastern (Longmen Shan Fault), northern (central–east Kunlun Fault) and northeastern (east Kunlun Fault) boundaries of the plateau, but with decreasing intensity. Segmented high upper-mantle velocities have varied occurrences and depth extents from south/southwest to north/northeast in the plateau. A Z-shaped upper-mantle low-velocity channel, which was taken as Tibetan lithospheric mantle, reflecting deformable material lies along the northern and eastern periphery of the Tibetan plateau, seemingly separating two large high-velocity mantle areas that, respectively, correspond to the Indian and Asian lithospheres. Other small high-velocity mantle segments overlain by the Z-shaped channel are possibly remnants of cold microplates/slabs associated with subductions/collisions prior to the Indian–Eurasian collision during the accretion of the Tibetan region. By integrating the Vs structures with known tectonic information, we derive that the Indian slab generally underlies the plateau south of the Bangong–Nujiang suture in central Tibet and the Jinsha River suture in eastern Tibet and west of the Lanchangjiang suture in southeastern Tibet. The eastern, northern, northeastern and southeastern boundaries of the Tibetan plateau have undergone deformation with decreasing intensity. The weakly resisting northeast and southeast margins, bounded by a wider softer channel of uppermost mantle material, are two potential regions for plateau expansion in the future.

P- and S-wave velocity structure of central West Antarctica: Implications for the tectonic evolution of the West Antarctic Rift System

New P- and S-wave velocity models of the upper mantle from 100 to 400 km depth beneath the central portions of West Antarctica, obtained by inverting relative travel-times from teleseismic earthquakes recorded on Polar Earth Observing Network (POLENET/ANET) and UK Antarctic Network (UKANET) seismic stations between 2007 and 2017, reveal a heterogeneous upper mantle. A low velocity anomaly (−1.0% Vp; −2.0% Vs) imaged beneath Marie Byrd Land is attributed to thermally perturbed upper mantle of possible plume origin, and a low velocity anomaly imaged beneath the Pine Island Glacier and the mouth of Thwaites Glacier is interpreted as a rift-related thermal structure that may include warm mantle flowing from Marie Byrd Land. High velocity anomalies (≤0.8% Vp; 1.5% Vs) imaged in the central portion of the West Antarctic Rift System indicate the presence of lithosphere unmodified by tectonic activity since the Late Cretaceous formation of the rift system. Within the region of high velocities, localized low velocity anomalies beneath parts of the Bentley Subglacial Trench are suggestive of focused Cenozoic rifting. The models also show variable velocity structure beneath the Haag-Ellsworth Whitmore crustal block and low velocities beneath the Thurston Island-Eights Coast crustal block. The heterogenous upper mantle structure of central West Antarctica indicates that upper mantle temperatures could vary by 100 K or more over distances of less than 100 km, which may add complexity to solid earth-ice interactions and influence basal ice sheet conditions.

Random objective waveform inversion of surface waves

Full-waveform inversion (FWI) of surface waves is becoming increasingly popular among shallow-seismic methods. Due to a huge amount of data and the high nonlinearity of the objective function, FWI usually requires heavy computational costs and may converge toward a local minimum. To mitigate these problems, we have reformulated FWI under a multiobjective framework and adopted a random objective waveform inversion (ROWI) method for surface-wave characterization. Three different measure functions were used, whereas the combination of one measure function with one shot independently provided one of the [Formula: see text] objective functions ([Formula: see text] is the total number of shots). We have randomly chose and optimized one objective function at each iteration. We performed a synthetic test to compare the performance of the ROWI and conventional FWI approaches, which showed that the convergence of ROWI is faster and more robust compared with conventional FWI approaches. We also applied ROWI to a field data set acquired in Rheinstetten, Germany. ROWI successfully reconstructed the main geologic feature, a refilled trench, in the final result. The comparison between the ROWI result and a migrated ground-penetrating radar profile further proved the effectiveness of ROWI in reconstructing the near-surface S-wave velocity model. We also ran the same field example by using a poor initial model. In this case, conventional FWI failed whereas ROWI still reconstructed the subsurface model to a fairly good level, which highlighted the relatively low dependency of ROWI on the initial model.

Elastic Full Waveform Inversion With Angle Decomposition and Wavefield Decoupling

Full waveform inversion (FWI) is a powerful tool to understand the real complicated earth model. As FWI is a highly nonlinear problem and depends strongly on the initial model, how to effectively retrieve the large-scale background model is critical for the success of FWI. For elastic FWI (EFWI), the inversion challenge increases because the P-wave and S-wave are coupled together if no mode separation technologies are applied. In this article, we develop a new EFWI strategy, where we simultaneously implement the angle decomposition and mode separation for the wavefield. Based on the analysis of radiation patterns of different parameters and the fact that small scattering angles correspond to large-scale model perturbations, we can retrieve the large-scale background model of the P-wave velocity with pure small scattering angle P-P mode wavefield. On the other hand, the pure small scattering angle S-S, S-P, and P-S mode wavefields are used to estimate the large-scale background model of the S-wave velocity. The correctly retrieved large-scale background models further guarantee the success of subsequent fine structure retrieving for the P- and S-wave velocity models by using different wave modes. The proposed method is able to reduce the cycle-skipping problem and the multiparameter crosstalk problem simultaneously. Numerical examples show that the proposed method provides much improved inversion results than the conventional EFWI, which demonstrates the validity of the proposed method.

3D wave-equation dispersion inversion of surface waves recorded on irregular topography

Irregular topography can cause strong scattering and defocusing of propagating surface waves, so it is important to account for such effects when inverting surface waves for shallow S-wave velocity structures. We have developed a 3D surface-wave dispersion inversion method that takes into account the topographic effects modeled by a 3D spectral element solver. The objective function is the frequency summation of the squared wavenumber differences [Formula: see text] along each azimuthal angle of the fundamental mode or higher-order modes of Rayleigh waves in each shot gather. The wavenumbers [Formula: see text] associated with the dispersion curves are calculated using the data recorded along the irregular free surface. Numerical tests on synthetic and field data demonstrate that 3D topographic wave equation dispersion inversion (TWD) can accurately invert for the S-wave velocity model from surface-wave data recorded on irregular topography. Field data tests for data recorded across an Arizona fault demonstrate that, for this example, the 2D TWD model can be as accurate as the 3D tomographic model. This suggests that in some cases, the 2D TWD inversion is preferred over 3D TWD because of its significant reduction in computational costs. Compared to the 3D P-wave velocity tomogram, the 3D S-wave tomogram agrees much more closely with the geologic model taken from the trench log. The agreement with the trench log is even better when the [Formula: see text] tomogram is computed, which reveals a sharp change in velocity across the fault. The localized velocity anomaly in the [Formula: see text] tomogram is in very good agreement with the well log. Our results suggest that integrating the [Formula: see text] and [Formula: see text] tomograms can sometimes give the most accurate estimates of the subsurface geology across normal faults.

2-D multiparameter viscoelastic shallow-seismic full-waveform inversion: reconstruction tests and first field-data application

SUMMARY 2-D full-waveform inversion (FWI) of shallow-seismic wavefields has recently become a novel way to reconstruct S-wave velocity models of the shallow subsurface with high vertical and lateral resolution. In most applications, seismic wave attenuation is ignored or considered as a passive modelling parameter only. In this study, we explore the feasibility and performance of multiparameter viscoelastic 2-D FWI in which seismic velocities and attenuation of P and S waves, respectively, and mass density are inverted simultaneously. Synthetic reconstruction experiments reveal that multiple crosstalks between all viscoelastic material parameters may occur. The reconstruction of S-wave velocity is always robust and of high quality. The parameters P-wave velocity and density exhibit weaker sensitivity and can be reconstructed more reliably by multiparameter viscoelastic FWI. Anomalies in S-wave attenuation can be recovered but with limited resolution. In a field-data application, a small-scale refilled trench is nicely delineated as a low P- and S-wave velocity anomaly. The reconstruction of P-wave velocity is improved by the simultaneous inversion of attenuation. The reconstructed S-wave attenuation reveals higher attenuation in the shallow weathering zone and weaker attenuation below. The variations in the reconstructed P- and S-wave velocity models are consistent with the reflectivity observed in a ground penetrating radar (GPR) profile.

Parameterization analysis and field validation of VTI-elastic full-waveform inversion in a walk-away vertical seismic profile configuration

Elastic full-waveform inversion (FWI) in transversely isotropic media with a vertical symmetry axis (VTI) is applied to field walk-away vertical seismic profile (W-VSP) data acquired in Western Canada. The performance of VTI-elastic FWI is significantly influenced by the model parameterization choice. Synthetic analysis based on specific field survey configuration is carried out to evaluate three different VTI-elastic model parameterizations. Interparameter trade-offs are quantified using the recently introduced interparameter contamination sensitivity kernel approach. Synthetic results suggest that neglecting anisotropy leads to inaccurate velocity estimations. For the conventional vertical velocity-Thomsen’s parameter parameterization (i.e., vertical P-wave velocity, vertical S-wave velocity, Thomsen’s parameters [Formula: see text] and [Formula: see text]), a sequential inversion strategy is designed to reduce strong natural interparameter trade-offs. The model parameterizations of elastic-constant ([Formula: see text], [Formula: see text], [Formula: see text], and [Formula: see text]) and velocity-based (vertical, horizontal, and normal move-out P-wave velocities and vertical S-wave velocity) models appear to suffer from fewer interparameter trade-offs, providing more reliable velocity and anisotropy models. Results derived from application of VTI-elastic FWI to the field W-VSP data set tend to support the synthetic conclusions. Multiparameter point spread functions are calculated to quantify the local interparameter trade-offs of the inverted models. The output inversion results are interpreted to provide valuable references regarding the target hydrocarbon reservoir.

Identification of a nascent tectonic boundary in the San-in area, southwest Japan, using a 3D S-wave velocity structure obtained by ambient noise surface wave tomography

We derived a three-dimensional S-wave velocity model for the San-in area of southwest Japan to examine heterogeneous structures such as tectonic faults. Many earthquakes occur in this area, but much of the activity has been relatively recent, so the fault distribution has yet to be fully clarified. Here, we used continuous ambient noise data from a dense seismic network, deployed from November 2009 to extract Rayleigh and Love wave dispersion data between station pairs, and then applied a direct surface wave inversion to the phase velocities of each station pair to determine a three-dimensional S-wave velocity model. In the resulting model, faults and a previously unrecognized tectonic boundary appeared as low-velocity anomalies or velocity boundaries, and the velocity anomalies were also associated with many past earthquake hypocenters. These results contribute to our understanding of heterogeneous structures caused by recent tectonic motion and of possible future tectonic activity, such as intraplate earthquakes. Surface wave tomography using ambient noise recorded in dense seismic networks could also be applied in other parts of the world to reveal new heterogeneous geological structures (i.e., unrevealed tectonic faults) and could contribute to disaster mitigation.

Multiple source downwellings beneath eastern North China revealed by 3-D CCP migration of receiver function data

We collected and processed a total of 97,700 teleseismic receiver functions recorded by 580 broadband seismic stations from 686 earthquakes to study the mantle transition zone beneath eastern North China. We employed a recently developed 3-D P- and S-wave velocity model (the East Asia Radially Anisotropic Model, EARA2014) to compute Pds moveouts to migrate receiver function data. Amplitudes of the stacked P660s and P410s using the 3-D Pds moveout table are approximately 30% higher than those calculated from the 1-D iasp91 model. The 3-D volume of CCP (common-conversion-point) stacked images that covers the area of 110–125°E and 36–43°N revealed two remarkable mantle transition zone anomalies with significant depression of the 660-km discontinuity. Together with other seismic observations, we speculated that these two anomalies are resulted from two different dynamic processes beneath eastern North China.

Elastic SH- and Love-wave Full-Waveform Inversion for shallow shear wave velocity with a preconditioned technique

In contrast to the Rayleigh waves, the system of SH- and Love-wave Full-Waveform Inversion (FWI) is simpler (Only S-wave velocity and density is inverted). However, the gradient is not scaled with the increasing depth in gradient-based FWI. The main reasons for the phenomenon are the limited frequency band of the wavelet, non-uniform coverage of the excitation-receiving points and the double scattering. The illumination can be enhanced by adding a preconditioned operator to the inversion. To compensate the illumination, we presented a preconditioned technique abiding by the theory of inverse scattering in SH- and Love-wave FWI. The synthetic tests demonstrate that the preconditioned technique can expedite the convergence rate and enhance the illumination intensity to the deeper medium compared with the conventional gradient-based FWI. The results also indicate that the preconditioned conjugate gradient (PCG) method has the best parameter reconstruction performance among the local optimization algorithms. Subsequently, we adopt the PCG method to inverse the field data acquired in the village of Čachtice, northwestern Slovakia. In addition, a multiscale strategy is employed and the Wiechert-Herglotz method is followed to estimate the initial S-wave velocity model in inversion. Eventually, this paper illustrates the credibility of the inverted result from the perspectives of waveform fitting, common-offset profile, signal similarity, relative error of channel energy and normalized misfit function values. This study primarily seeks to delineate the foundations of historical buildings under the survey line. According to the range of S-wave velocity, we assert that there are no foundation remains in the subsurface.

Two-station continuous wavelet transform cross-coherence analysis for surface-wave tomography using active-source seismic data

ABSTRACT Surface-wave tomography has great potential to improve the lateral resolution of near-surface characterization compared to 2D surface-wave analysis with multichannel analysis of surface waves (MASW). Surface-wave tomography has been widely applied to obtain high-resolution maps of phase or group velocity from dispersion curves between pairs of stations in seismological studies. However, very few studies have done surface-wave tomography with active-source (exploration) seismic data, probably because extracting surface-wave dispersion curves between two stations is difficult due to the complex wave propagation in heterogeneous near-surface structures. Here, we describe a method to estimate reliable phase-velocity dispersion curves between two stations from exploration seismic data. In our approach, we compute cross coherences between pairs of stations to extract phase information, stacking the cross coherences from different shot gathers to improve the signal-to-noise ratio. To further distinguish surface-wave signals from noise in the time domain, we perform a time-frequency analysis using the continuous wavelet transform (CWT) on the stacked cross coherences. We used modeling of the wavelet transform between station pairs to extract phase-velocity dispersion curves from the stacked cross coherences. We apply this two-station CWT cross-coherence method to synthetic and field data sets. Both applications demonstrate that our method can extract stable phase-velocity dispersion curves between two stations better than two-station or multistation analysis without time-domain filtering. In phase-velocity distributions constructed by surface-wave tomography from the dispersion curves between two stations, the horizontal resolution is improved over MASW-based analyses. Improvement of the horizontal resolution is also achieved in S-wave velocity structures derived by inversion of the phase-velocity distributions. Our method is effective in estimating reliable phase-velocity dispersion curves and may contribute to constructing high-resolution S-wave velocity models located with a laterally heterogeneous structure, by subsequent surface-wave tomography and S-wave velocity inversion.

Constraints on S-wave velocity structures of the lithosphere in mainland China from broadband ambient noise tomography

In this study, a 3-D S-wave velocity model of the lithosphere in mainland China down to 150 km depth was determined from ambient noise data. First, we collected 2-year continuous waveforms recorded by 1031 broadband stations of the China Regional Seismic Network and NECESSArray. Then, by applying the procedures of noise cross-correlation and time-frequency domain phase-weighted stacking, we obtained the interstation empirical Green's functions of Rayleigh waves. We measured the group and phase velocity dispersions at periods of 5–125 s from the EGFs and inverted both group and phase velocities for the S-wave velocity model. Our model revealed that the lithospheric thickness in the northern Songliao basin is ~70–110 km. A weak high-velocity anomaly that is discontinuous in the vertical direction was detected beneath the southern Songliao basin, which may indicate that the lithosphere is delaminated in this area. The model also revealed that deep (>150 km) high-velocity lithospheric roots exist beneath the Ordos and Sichuan basins. The lithospheric thicknesses in eastern parts of the Northeast China Block, the North China Craton, and the South China Block are only 60–100 km, which may be the result of a series of deep processes caused by the subduction of the Pacific plate to the Eurasian plate that produced strong lithospheric thinning in eastern China. The lithospheric thickness within the Tanlu fault belt is generally thinner than that in the adjacent areas of the fault belt, and the S-wave velocity of the mantle lithosphere beneath the fault belt is also lower than that beneath the adjacent areas, indicating that the Tanlu fault belt is a deep fault belt cutting through the lithosphere.

The formation of the Dabashan orocline, central China: Insights from high-resolution 3D crustal shear-wave velocity structure

The Dabashan orocline is a large thrust belt in central China. Determination of its high-resolution three-dimensional (3D) crustal shear-wave (S-wave) velocity structure can enhance our understanding of the intracontinental evolution of the South China plate and North China plate. In this study, we estimated the Rayleigh wave phase and group velocities at periods of 5–50 s, P-wave receiver functions, and Ps and SsPmp travel times from 63 permanent stations and 29 portable stations; we also constructed a 3D S-wave velocity model and crustal thickness and Vp/Vs ratio maps of the Dabashan orocline and adjacent areas. The crustal thickness and Vp/Vs ratio maps show prominent thick crust (50–55 km) and high Vp/Vs ratios (1.8–1.85), suggesting a mafic lower crust in the Dabashan orocline. The 3D S-wave velocity structure shows northeast-dipping low-velocity anomalies in the upper crust of the Dabashan orocline, and high velocities extended from near the surface to the lower crust in the Hannan-Micang and Shennong-Huangling domes, suggesting deep roots of the two domes. We propose that the two domes acted as rigid indenters and rotated clockwise together with the South China plate, penetrating into the Qinling-Dabie orogenic belt. Weak sediments in the upper crust of the Dabashan orocline were compressed and blocked by the two domes and the rigid Qinling-Dabie orogenic belt, resulting in vertical accretion of the sediments and crust and subsequent formation of the northeast-dipping low-velocity anomalies and the thick crust beneath the Dabashan orocline.

Seismic exploration survey northeast of the Tres Virgenes Geothermal Field, Baja California Sur, Mexico: A new Geothermal prospect

Geophysical prospecting has become essential to identify low enthalpy areas with geothermal potential, mainly, for direct uses. The present study uses two seismic techniques and lithological information from shallow boreholes located northeast of the Tres Virgenes Geothermal Field, B.C.S., Mexico, to identify geological structures of geothermal interesting. Using the Multichannel Analysis of Surface Waves method (MASW), it was possible to correlate the S-wave velocity model with the lithology from two wells drilled over the seismic profile. Additionally, we generated a seismic reflection model to describe the geological distribution of deeper structures. The reflection profile, which reached depths of about 500 m, was complemented with superficial information from a MASW velocity model. We found a relationship between the velocity distribution of P and S waves from both models. The integration of seismic and lithological information allows us to identify a very fractured subsurface and a deep zone with possible geothermal interest.

3D wave-equation dispersion inversion of Rayleigh waves

The 2D wave-equation dispersion (WD) inversion method is extended to 3D wave-equation dispersion inversion of surface waves for the shear-velocity distribution. The objective function of 3D WD is the frequency summation of the squared wavenumber [Formula: see text] differences along each azimuth angle of the fundamental or higher modes of Rayleigh waves in each shot gather. The S-wave velocity model is updated by the weighted zero-lag crosscorrelation between the weighted source-side wavefield and the back-projected receiver-side wavefield for each azimuth angle. A multiscale 3D WD strategy is provided, which starts from the pseudo-1D S-velocity model, which is then used to get the 2D WD tomogram, which in turn is used as the starting model for 3D WD. The synthetic and field data examples demonstrate that 3D WD can accurately reconstruct the 3D S-wave velocity model of a laterally heterogeneous medium and has much less of a tendency to getting stuck in a local minimum compared with full-waveform inversion.