High-velocity Layer Research Articles

Empirical H/V spectral ratios at the InSight landing site and implications for the martian subsurface structure

SUMMARY The horizontal-to-vertical (H/V) spectral ratio inversion is a traditional technique for deriving the local subsurface structure on Earth. We calculated the H/V from the ambient vibrations at different wind levels at the InSight landing site, on Mars, and also computed the H/V from the S-wave coda of the martian seismic events (marsquakes). Different H/V curves were obtained for different wind periods and from the marsquakes. From the ambient vibrations, the recordings during low-wind periods are close to the instrument self-noise level. During high-wind periods, the seismic recordings are highly contaminated by the interaction of the lander with the wind and the martian ground. Therefore, these recordings are less favourable for traditional H/V analysis. Instead, the recordings of the S-wave coda of marsquakes were preferred to derive the characteristic H/V curve of this site between 0.4 and 10 Hz. The final H/V curve presents a characteristic trough at 2.4 Hz and a strong peak at 8 Hz. Using a full diffuse wavefield approach as the forward computation and the Neighbourhood Algorithm as the sampling technique, we invert for the 1-D shear wave velocity structure at the InSight landing site. Based on our inversion results, we propose a strong site effect at the InSight site to be due to the presence of a shallow high-velocity layer (SHVL) over low-velocity units. The SHVL is likely placed below a layer of coarse blocky ejecta and can be associated with Early Amazonian basaltic lava flows. The units below the SHVL have lower velocities, possibly related to a Late Hesperian or Early Amazonian epoch with a different magmatic regime and/or a greater impact rate and more extensive weathering. An extremely weak buried low velocity layer (bLVL) between these lava flows explains the data around the 2.4 Hz trough, whereas a more competent bLVL would not generate this latter feature. These subsurface models are in good agreement with results from hammering experiment and compliance measurements at the InSight landing site. Finally, this site effect is revealed only by seismic events data and explains the larger horizontal than vertical ground motion recorded for certain type of marsquakes.

Geophysical insights on the crustal structure of Greenland's northern continental margin towards the Morris Jesup Spur

During expedition PS115/1, the German research vessel Polarstern acquired seismic refraction data along a 102-km-long profile crossing Greenland's northern continental margin and extending up to the southwestern limit of the Morris Jesup Spur. A P–wave velocity model is obtained and validated by gravity modelling. A nearby seismic reflection line provides insights on the structures within the sedimentary cover. Beneath a 2-km-thick sedimentary cover with velocities of 1.8 km/s to 3.4 km/s, an up to 1.5-km-thick layer is characterized by velocities of 4.2 km/s and is interpreted to consist of volcanic rocks. This is consistent with proposed volcanic activity on the Morris Jesup Spur and exposed volcano-sedimentary rocks of the nearby Kap Washington Group. Below the volcanic rocks, an up to 7 km-thick unit with velocities of 4.4 to 5.8 km/s is interpreted to consist of metasedimentary rocks still belonging to the deformed units of the Franklinian Basin. The seismic reflection data image a tectonic overprint of the acoustic basement and the lowermost sedimentary layers. Tectonic faults indicating the tectonic overprint do not reach to the uppermost sedimentary layers. This hints at an older fracture zone that is also observed in the velocity model and aligns well with observed anomalies in the regional gravity and magnetic fields. The seismic velocity model reveals a highly extended continental crust that shows a magmatic overprint at all crustal levels. The upper crust is 5 km thick with a velocity of 6.0 km/s. However, in three distinct zones its thickness is up to 8 km and the velocity increases to 6.3 km/s. These zones are interpreted as magmatic intrusions into the upper crust. The lower crust represents a 9-km-thick high-velocity layer (7.2 km/s) that is interpreted as magmatic underplating or lower crustal sill intrusions. Such a high-velocity lower crust is not present in models of adjacent Arctic margins.PACS: 0000, 11112000 MSC: 0000, 1111

A combined GA and DLS parallel algorithm for multimode Rayleigh wave inversion

Multichannel analysis of surface waves (MASW) is a vital method in near-surface exploration. Linear algorithms, such as damped least square method (DLS), are highly dependent on the initial model and are prone to trapping in local minima. While nonlinear algorithms such as genetic algorithm (GA) are globally optimal but inefficient, they are not very practical for complex calculations.In this paper, we propose a combined GA and DLS parallel algorithm (GA-DLS), which has a global search capability due to the GA framework, and has a quick local search capability due to the DLS parallel computing. To simulate the earth model in complex near-surface conditions, we set a low-velocity layer among a layered earth model (LVL) and a high-velocity layer among a layered earth model (HVL), we perform the inversion with DLS, GA, and GA-DLS, respectively. The comparison shows that the GA-DLS inversion gets an efficient and accurate result. Finally, we use the GA-DLS algorithm to invert the field data and achieve good results, verifying the effectiveness of the algorithm.

Breakup mechanism of the northern South China Sea: Evidence from the deep crustal structure across the continent-ocean transition

To improve constraints on breakup processes that eventually lead to seafloor spreading and the subsequent formation of new oceanic crust in the northern margin of the South China Sea (SCS), we collect a 172-km-long deep-penetration wide-angle seismic data OBS2018-L5. In the continental domain, the results present a gradually thinning crust with thickness ranging from ∼17 km to ∼8 km, with a slight thickening crust at the Outer Margin High. The attenuated continental crust is fairly rough and highly offset by distinctly different scale faults system. In the oceanic domain, the top of the basement is highly reflective with weak faulting and relatively constant crust thickness ranging from 5 km to 8 km (average 6 km). The continent-ocean boundary is characterized by an abrupt change in crustal thickness accompanied by a specific boundary between mostly positive gravity anomalies in the oceanic basin and grossly negative anomalies in the thinned continental crust, and between low amplitude, out-of-shape magnetic anomalies in the thinned continental domain and high amplitude, linear anomalies in the oceanic domain. A 0–5 km thick high-velocity layer (7.0–7.5 km/s) is present in the distal margin, which might represent mafic magma underplating. We also identify other magmatic features such as distinct buried volcanic edifices overprinting the original velocity structure. We propose that three mechanisms control the extension and breakup of the SCS: (1) the rapid transition from continental to oceanic crust, (2) elevated potential mantle temperature, and (3) the ductile lower crustal necking. Their interactions lead to the intense thinning of the continental lithosphere, decompression melting of the enriched asthenosphere followed by fast upwelling of mafic melt, relatively strong volcanic activities during the final stage of syn-rift and breakup, and finally generating a margin with intermediately rich magma, referred to as in intermediate-type margin.

Seismic Imaging of Complex Velocity Structures by 2D Pseudo-Viscoelastic Time-Domain Full-Waveform Inversion

In the presented study, multi-parameter inversion in the presence of attenuation is used for the reconstruction of the P- and the S- wave velocities and the density models of a synthetic shallow subsurface structure that contains a dipping high-velocity layer near the surface with varying thicknesses. The problem of high-velocity layers also complicates selection of an appropriate initial velocity model. The forward problem is solved with the finite difference, and the inverse problem is solved with the preconditioned conjugate gradient. We used also the adjoint wavefield approach for computing the gradient of the misfit function without explicitly build the sensitivity matrix. The proposed method is capable of either minimizing the least-squares norm of the data misfit or use the Born approximation for estimating partial derivative wavefields. It depends on which characteristics of the recorded data—such as amplitude, phase, logarithm of the complex-valued data, envelope in the misfit, or the linearization procedure of the inverse problem—are used. It showed that by a pseudo-viscoelastic time-domain full-waveform inversion, structures below the high-velocity layer can be imaged. However, by inverting attenuation of P- and S- waves simultaneously with the velocities and mass density, better results would be obtained.

Crust and upper mantle structure beneath Southwest China and its implications for Mesozoic multistage gold deposits

The Youjiang Basin (YB) on the southwestern margin of the Yangtze Craton (YC) contains >200 Triassic or Cretaceous gold deposits, but deep dynamic processes contributing to large-scale gold accumulation remain unknown. Based on a south–north trending linear broadband seismic array, we obtained the crust and mantle structure beneath the YB and southwestern YC by P-wave and S-wave receiver functions and ambient noise tomography. Our results show a remarkable lithospheric boundary beneath the northern margin of the YB, which separates the YB with felsic thin crust and thin mantle lithosphere from the YC with both a thick crust and thick mantle lithosphere. A new structural image shows that the southern end of the YC Moho dips southward and overlaps with the YB Moho. A 5-km-thick high-velocity layer is also found at a depth of 10 km beneath the northern YB. Combining this information with previous tectonic and geological results, we reveal the process of Mesozoic lithospheric evolution in the southwestern South China Block. We infer that both the pre-Mesozoic weak crust and mantle lithosphere and the intense intracontinental orogenesis in the Early-Middle Triassic contributed to the formation of first-stage Au deposits in the YB. Lithospheric delamination in the Early Cretaceous thinned the lithosphere of YB rather than YC, forming a steep lithosphere–asthenosphere boundary beneath the transition zone between YB and YC, which favored upwelling of Au-rich asthenosphere below the thick cratonic root to the bottom of the crust in the northern YB and contributed to the second-stage Au deposits in the YB.

High-Frequency and Spectrum-Clean Shear-Horizontal Acoustic Wave Resonators with AlN Overlay.

By bonding the sub-wavelength-thick lithium niobate (LiNbO3) layer to high-phase-velocity (vp) substrates, such as Si, the shear-horizontal (SH) modes no longer couple with the bulk modes leaking into substrates. As the propagation loss is no longer the major concern for these types of nonleaky SH wave devices, the YX-LiNbO3 with a low rotation angle providing ultra-large coupling coefficient (keff2) can be used. In addition, by overlaying a high-velocity layer such as AlN on top of LiNbO3/Si, the vp of the SH wave can be significantly enhanced at a small cost of keff2. By a careful design of the stack, both the wide-band spurious (Lamb wave) and near-band spurious (Rayleigh wave) are suppressed successfully. This paper focuses on the design of layered substrate not only to optimize its resonance characteristics—series frequency (fs), quality factor (Q), keff2, and temperature coefficient of frequency (TCF)—but also for eliminating the out-of-band spurious responses. The optimized substrate design demonstrates the minimal propagation loss, high fs of 3 GHz, large keff2 of 14.4% and a spurious-free response at 0–6 GHz. These novel nonleaky SH wave devices can potentially enable the low loss and wideband processing functions, which is promising for the 5G/6G radio frequency (RF) communication systems.

Determination of near-surface shear-velocity structure based on the joint inversion of Rayleigh-wave dispersion and ellipticity from multistation active-seismic records

The dispersive characteristics of high-frequency (≥2 Hz) Rayleigh waves have been widely applied to determine subsurface shear velocity ([Formula: see text]). However, the [Formula: see text] of the high-velocity surface layers (e.g., rigid roads in urban areas) and underlying layers is poorly constrained in the Rayleigh-wave dispersion inversion. Compared to the phase velocity, the Rayleigh-wave horizontal-to-vertical (H/V) ratios of particle motion (i.e., ellipticity) are particularly sensitive to shallow [Formula: see text] structures. Therefore, we have proposed a robust method for obtaining near-surface [Formula: see text] structures by inverting Rayleigh-wave dispersion and H/V data. Initially, we separate the fundamental-mode Rayleigh waves from the active-seismic records using a high-resolution linear Radon transformation method. Then, we extract the high-frequency H/V and interstation phase traveltimes from the separated fundamental mode using deconvolution interferometry and the frequency-time analysis method. Based on the phase traveltimes, we calculate the phase velocity dispersion using ray-based tomography. Finally, we perform the joint inversion of dispersion and H/V data to obtain the [Formula: see text] structure. According to the inversion results of synthetic and field data, the accuracy of the [Formula: see text] in the rigid surface layer and underlying layers is considerably improved. We find that the joint inversion of Rayleigh-wave phase velocity dispersion and H/V data from surface-wave surveys is a useful tool for accurately determining the complex shallow [Formula: see text] structures.

Chemical Stratification in a Long Gamma-Ray Burst Cocoon and Early-time Spectral Signatures of Supernovae Associated with Gamma-Ray Bursts

Abstract We present the results of 3D hydrodynamic simulations of a gamma-ray burst (GRB) jet emanating from a massive star with a particular focus on the formation of high-velocity quasi-spherical ejecta and the jet-induced chemical mixing. Recent early-time optical observations of supernovae associated with GRBs (e.g., GRB 171205A/SN 2017iuk) indicate a considerable amount of heavy metals in the high-velocity outer layers of the ejecta. Using our jet simulations, we show that the density and chemical structure of the outer ejecta implied by observations can be naturally reproduced by a powerful jet penetrating the progenitor star. We consider three representative jet models with a stripped massive star, a standard jet, a weak jet, and a jet choked by an extended circumstellar medium, to clarify the differences in the dynamical evolution and the chemical properties of the ejected materials. The standard jet successfully penetrates the progenitor star and creates a quasi-spherical ejecta component (cocoon). The jet-induced mixing significantly contaminates the cocoon with heavy elements that have been otherwise embedded in the inner layer of the ejecta. The weak and choked jet models fail to produce an ultrarelativistic jet but produce a quasi-spherical cocoon with different chemical properties. We discuss the impact of the different jet−star interactions on the expected early-time electromagnetic signatures of long GRBs and how to probe the jet dynamics from observations.

Walkaway Vertical Seismic Profiling (WVSP) Modeling and Imaging Study along FaultedCoal Seams over a High Velocity Limestone Model: A Synthetic Study

The problem of collecting reflection data from the layers below a high velocity layer (HVL) and imaging those layers accurately is solved using walkaway vertical seismic profiling (WVSP) technique. The procedure is shown on a dataset from a model of highly faulted, thin-bedded coal layers separated by a high velocity layer of limestone. By locating the receivers above, in, and below the HVL, the reflections are expected to be recorded locally in situ in the borehole below the HVL before they travel back to the surface as in the case of surface seismic. The imaging results obtained from the model (or synthetic) WVSP data lead to the conclusion that promising results may be obtained via the acquisition of field WVSP data such that the layers below the HVL are well-imaged. Analysis of the test results showed that the layers below the HVL can be imaged accurately and robustly when the receivers are situated below the HVL only, while the strong images below the HVL tend to fade away and lose their strength when the receivers are moved above the HVL. However, both these cases of the receiver location geometry are needed to obtain images that are vertically and horizontally stable. When a full instrumentation of the borehole is provided, good images as far as 500 m away from the VSP borehole are obtained from the layers with various dips and faults that are nearly 800 m in depth from the surface. Two cases of velocities (4000 m/s and 5000 m/s) are tested for the high velocity limestone layer. Both cases showed good images below the HVL. Unfortunately, some pulled-up effects are seen in the images right below the HVL where the layers that are immediately below the HVL are imaged above their correct locations, the pulled-up effects are corrected in the deeper sections below the HVL however. Comparison of a depth converted corridor stack from the zero offset VSP data and the depth migrated image show good agreement.

Walkaway Vertical Seismic Profiling (WVSP) Modeling and Imaging Study along FaultedCoal Seams over a High Velocity Limestone Model: A Synthetic Study

The problem of collecting reflection data from the layers below a high velocity layer (HVL) and imaging those layers accurately is solved using walkaway vertical seismic profiling (WVSP) technique. The procedure is shown on a dataset from a model of highly faulted, thin-bedded coal layers separated by a high velocity layer of limestone. By locating the receivers above, in, and below the HVL, the reflections are expected to be recorded locally in situ in the borehole below the HVL before they travel back to the surface as in the case of surface seismic. The imaging results obtained from the model (or synthetic) WVSP data lead to the conclusion that promising results may be obtained via the acquisition of field WVSP data such that the layers below the HVL are well-imaged. Analysis of the test results showed that the layers below the HVL can be imaged accurately and robustly when the receivers are situated below the HVL only, while the strong images below the HVL tend to fade away and lose their strength when the receivers are moved above the HVL. However, both these cases of the receiver location geometry are needed to obtain images that are vertically and horizontally stable. When a full instrumentation of the borehole is provided, good images as far as 500 m away from the VSP borehole are obtained from the layers with various dips and faults that are nearly 800 m in depth from the surface. Two cases of velocities (4000 m/s and 5000 m/s) are tested for the high velocity limestone layer. Both cases showed good images below the HVL. Unfortunately, some pulled-up effects are seen in the images right below the HVL where the layers that are immediately below the HVL are imaged above their correct locations, the pulled-up effects are corrected in the deeper sections below the HVL however. Comparison of a depth converted corridor stack from the zero offset VSP data and the depth migrated image show good agreement.

Application of active-source surface waves in urban underground space detection: a case study of Rongcheng County, Hebei, China

Active-source surface wave exploration is advantageous because it has high imaging accuracy, is not affected by high-speed layers and has a low cost. It has unique advantages in the investigation of shallow surface structures. For the development and utilization of urban underground space, there are two important parameters in the shallow surface, namely, shear wave velocity (vs.) and site predominant period, which determine the elevation and aseismic grade of the building design. The traditional method mainly obtains the above two parameters through the testing and measurement of drilling samples. However, this method is extremely expensive and time consuming. Therefore, this paper uses the multichannel surface wave acquisition method to extract the fundamental dispersion curve of the single shot data using the phase shift method and obtains the vs. characteristics in the uppermost 40 m by inversion. According to the vs. profile, the following two conclusions were obtained. First, the study area can be roughly divided into 5 layers, among which the layers 0-8 m, 14-20 m, and 20-30 m are low-velocity layers, corresponding to miscellaneous fill, a water-bearing sand layer and a sand layer. Therefore, the vs. is relatively low, and the layers at 8-14 m and 30-40 m are high-velocity layers that are mainly composed of clay, with a relatively better compactness and relatively high vs. values. In addition, the low-speed anomaly suddenly appears in the high-speed area at 20-40 m. Combining this with geological data suggests that it is an ancient river channel. Second, according to the vs. value, the V_se(equivalent shear wave velocity) was calculated. The construction site soil is classified as Category III, with good engineering geological conditions. In addition, according to vs., we calculate the site predominant period to be 0.56-0.77 s. Therefore, in the overall structural design of foundation engineering, the natural vibration period of the structure should be strictly controlled to avoid the predominant period of the site.

From Subduction to LLSVP: The Core‐Mantle Boundary Heterogeneities Across North Atlantic

AbstractThe distribution of Ultra‐low Velocity Zones (ULVZs) is generally associated with the Large Low Shear Velocity Zones (LLSVPs). However, subducted slabs may be important for the characteristics of the ULVZs as well. Here we use ScS precursor (SdS) and postcursor (ScscS) phases, as well as Scd phases recorded by USArray to investigate the lowermost mantle structure across the north Atlantic, which spans the African LLSVP and the subducted Farallon slab. Our results show a ULVZ located near the edge of the African LLSVP, with thicknesses ranging from 11–16 km and S wave velocity reductions of 12%–24%. At the same time, a thin high‐velocity layer (HVL) with thicknesses of 20–30 km, Vs increases of 3%–4%, is also detected outside the LLSVP in the northwest region. Further west there is a clear D" layer at 250 km above the core mantle boundary (CMB), overlaid by a low velocity layer with Vs reduction of 3%. This transition profile from D" to HVL to ULVZ provides an excellent opportunity for us to understand the heterogeneities near the CMB region. We propose that interactions between the subducted Farallon slab with ambient mantle and LLSVP material, together with thermal effects, may result in this transition between different seismic structures along the CMB.

An Investigation of Seismic Velocity Variation through a Tectonic Boundaries-Case Study in Central Iraq

A 3D velocity model was created by using stacking velocity of 9 seismic lines and average velocity of 6 wells drilled in Iraq. The model was achieved by creating a time model to 25 surfaces with an interval time between each two successive surfaces of about 100 msec. The summation time of all surfaces reached about 2400 msec, that was adopted according to West Kifl-1 well, which penetrated to a depth of 6000 m, representing the deepest well in the study area. The seismic lines and well data were converted to build a 3D cube time model and the velocity was spread on the model. The seismic inversion modeling of the elastic properties of the horizon and well data was applied to achieve a corrected velocity cube. Then, the velocity cube was converted to a time model and, finally, a corrected 3D depth model was obtained. This model shows that the western side of the study area, which is a part of the stable shelf, is characterized by relatively low thickness and high velocity layers. While the eastern side of the study area, which is a part of the Mesopotamian, is characterized by high thickness and low velocity of the Cretaceous succession. The Abu Jir fault is considered as a boundary between the stable and unstable shelves in Iraq, situated at the extreme west part of the study area. The area of relatively high velocity gradient is considered as the limit of the western side of the Mesopotamian basin. This area extends from Najaf-Karbala axis in the west to the Euphrates River in the east. It is found that the 3D stacking velocity model can be used to obtain good results concerning the tectonic boundary.

Tailoring microstructure and mechanical property of laser powder deposited Inconel 718 superalloy via hybrid processing strategies

Customized tailoring of the microstructure and mechanical properties is a longstanding and challenging requirement for the high-performance implementations of laser additive manufacturing technology in the aircraft industry. In this study, two sets of processing parameters were alternated layer by layer, forming hybrid processing strategies. Thin-wall builds of Inconel 718 superalloy were produced by laser powder deposition process to preliminarily study the effect of hybrid processing strategies on tailoring the microstructure and mechanical properties. The results indicate the hybrid processing strategies have a distinct influence on the crystal growth behavior, grain orientations and sizes, precipitations and associated tensile strengths. Under the bidirectional scanning pattern, the hybrid processing strategy could effectively tailor the macroscale orientation of coarse zigzag microstructure through adjusting and disturbing the primary-secondary dendritic arms alternate growth behavior of columnar crystals. The columnar crystal trunk size exhibits a banded feature with large value in low-velocity deposition layer and small value in high-velocity layer. Under the unidirectional scanning pattern, the hybrid processing strategy can distinctly balance the size and aspect ratio of oriented fibrous grains. The strong texture of microstructure acts as an important contributor to determine the fracture behavior and associated ultimate tensile strength. Comparing to the single processing strategies, the hybrid processing strategies show little effect on the ultimate tensile strength with a slightly increased anisotropy under bidirectional scanning pattern, while distinctly improve the ultimate tensile strength with a reduced anisotropy under unidirectional scanning pattern, showing the potential to tailor the microstructure distribution and mechanical properties with customized high-performance requirements by purposely adjusting the hybridization of processing parameters.

Using ambient noise tomography and MAPS for high resolution stratigraphic identification in Hangzhou urban area

Near-surface (<200 m) stratigraphic identification is of scientific and engineering significance to urban sustainable development and expansion. Complex urbanized environments, however, limit the application of traditional active-source seismic surveys. Passive seismic surveys provide an alternative for investigating near-surface structures by utilizing abundant traffic noise. We presented a case study of fine-scale shear-wave velocity (Vs) structure imaging in the Hangzhou urban area (eastern China) and compared inversion results of ambient noise tomography and multichannel analysis of passive surface wave (MAPS). Noise data were recorded within one day by a dense seismic network (130 m × 120 m). The average adjacent station interval is about 10 m. Our goal is to investigate the top 100 m of the subsurface at the survey area. The two inverted Vs models by tomography and MAPS are similar and they both tally with the logging data. The subsurface geology structures are analyzed based on tomography inversion results and borehole data. The bedrock surface is determined at about 60 m depth and a high-velocity layer is found at around 45– 50 m depth. The sedimentary formation and bedrock in the survey area are further stratified under the control of borehole data.

Large-offset P-wave traveltime in layered transversely isotropic media

Large-offset seismic data processing, imaging, and velocity estimation require an accurate traveltime approximation over a wide range of offsets. In layered transversely isotropic media with a vertical symmetry axis, the accuracy of traditional traveltime approximations is limited to near offsets. We have developed a new traveltime approximation that maintains the accuracy of classic equations around the zero offset and exhibits the correct curvilinear asymptote at infinitely large offsets. Our approximation is based on the conventional acoustic assumption. Its equation incorporates six parameters. To define them, we use the Taylor series expansion of the exact traveltime around the zero offset and a new asymptotic series for the infinite offset. Our asymptotic equation indicates that the traveltime behavior at infinitely large offsets is dominated by the properties of the layer with the maximum horizontal velocity in the sequence. The parameters of our approximation depend on the effective zero-offset traveltime, the normal moveout velocity, the anellipticity, a new large-offset heterogeneity parameter, and the properties of the layer with the maximum horizontal velocity in the sequence. We have applied our traveltime approximation (1) to directly calculate traveltime and ray parameter at given offsets, as analytical forward modeling, and (2) to estimate the first four of the aforementioned parameters for the layers beneath a known high-velocity layer. Our large-offset heterogeneity parameter includes the layering effect on the reflections’ traveltime.

Rock probes and seismic methods to constrain the structure, composition and evolution of the deep crust beneath North China Block

<p indent="0mm">Deep continental crust archives critical information on the bulk crustal composition as well as its formation. It is commonly the inherent zone for complex crust-mantle interaction. Unlike the shallow crust, the deep crust is generally inaccessible for direct geological sampling. Instead, knowledge of the deep crust can be obtained from a combination study of geophysical techniques and “lithoprobe” (deep crustal fragments entrained in magmatic outcrops). Wide-angle reflection/refraction seismic profiles have the highest resolution on crustal velocity structure and provide information on the averaged present-day nature of large deep crust domains; “lithoprobe”, instead, give detailed information of locality-specific deep crustal sections through determining the physical and chemical properties, dating the formation and alteration ages and constraining the melting and equilibration depths of the xenoliths. Here, the two methodologies are integrated to constrain the crustal structure on the North China Block. The seismic crustal structure model of the block is compiled from all available seismic wide-angle reflection/refraction profiles and receiver function results. The model parameters include Moho depth, <italic>V</italic>p velocity and thicknesses for four crustal layers (sedimentary cover, and the upper, middle and lower crust). The crust as a whole is thin (&lt; <sc>36 km)</sc> in the eastern part relative to the western part <sc>(38−44 km).</sc> In addition, the crust is shown to have a thick upper and middle crust and a relatively thin (typically <sc>&lt;5 km)</sc> lower mafic crust; no spatial variations in average crustal <italic>V</italic>p velocity and densities are observed. The average crustal <italic>V</italic>p and density for the whole North China Block are <sc>6.33−6.37 km/s</sc> and <sc>2.72−2.78 g/cm³,</sc> respectively. The overall absence of high P-wave velocity layer and the relatively low average <italic>V</italic>p and density indicate that the block has an intermediate to felsic bulk crustal composition. In contrast, regional crustal xenoliths include mafic to felsic granulites and amphibolites and rare eclogite xenoliths (e.g., Xinyang and Xuhuai), indicating the presence of thick lower crust (e.g., <sc>&gt;45 km)</sc> at that time. Meanwhile, mafic xenoliths in Hannuoba and Nüshan show strong mafic underplating at the lower crust in the Meso-Cenozoic while no high-velocity layer is shown in the seismic crustal model. These inconsistences can be explained by the fact that the two methodologies tap the crustal information at different timeslices and dimensions and be complementary to understand the chemical and physical structure of the deep crust. Therefore, the combination study of petrological and geochemical of the crustal xenoliths and geophysical observations can provide essential information on making and destroying the deep continental crust through time.

Study on reflection spherical wave at the near surface

In real seismic surveys, seismic wave propagates in the form of spherical wave and plane wave is the approximations of spherical wave in the far-field and for high frequencies. The traditional seismic processing and interpretation methods based on the plane wave theory are not exact. Furthermore the difference between the plane wave and spherical wave is strong, especially in the near surface. Studying the characteristics of reflection spherical wave near the surface can develop more accurate wave theory. In this paper, a two-layer model with a high-velocity upper layer is built to analyze the spherical reflection wave characteristics. Tests show when the incident angle reaches a certain level, there will be an irregular wave, the pseudo-spherical wave, following the main reflection spherical wave, which has a great influence on the main reflection wave. This discovery has enriched our understanding of wave propagation theory and it has important guiding significance for accurate study on wave phenomena.

Passive seismic recording of cryoseisms in Adventdalen, Svalbard

Abstract. A series of transient seismic events were discovered in passive seismic recordings from 2-D geophone arrays deployed at a frost polygon site in Adventdalen, Svalbard. These events contain a high proportion of surface wave energy and produce high-quality dispersion images using an apparent offset re-sorting and inter-trace delay minimisation technique to locate the seismic source, followed by cross-correlation beamforming dispersion imaging. The dispersion images are highly analogous to surface wave studies of pavements and display a complex multimodal dispersion pattern. Supported by theoretical modelling based on a highly simplified arrangement of horizontal layers, we infer that a ∼3.5–4.5 m thick, stiff, high-velocity layer overlies a ∼30 m thick layer that is significantly softer and slower at our study site. Based on previous studies we link the upper layer with syngenetic ground ice formed in aeolian sediments, while the underlying layer is linked to epigenetic permafrost in marine-deltaic sediments containing unfrozen saline pore water. Comparing events from spring and autumn indicates that temporal variation can be resolved via passive seismic monitoring. The transient seismic events that we record occur during periods of rapidly changing air temperature. This correlation, along with the spatial clustering along the elevated river terrace in a known frost polygon, ice-wedge area and the high proportion of surface wave energy, constitutes the primary evidence for us to interpret these events as frost quakes, a class of cryoseism. In this study we have proved the concept of passive seismic monitoring of permafrost in Adventdalen, Svalbard.