Surface wave inversion with unknown number of soil layers based on a hybrid learning procedure of deep learning and genetic algorithm

Surface wave inversion is a key step in the application of surface waves to soil velocity profiling.Currently,a common practice for the process of inversion is that the number of soil layers is assumed to be known before using heuristic search algorithms to compute the shear wave velocity profile or the number of soil layers is considered as an optimization variable.However,an improper selection of the number of layers may lead to an incorrect shear wave velocity profile.In this study,a deep learning and genetic algorithm hybrid learning procedure is proposed to perform the surface wave inversion without the need to assume the number of soil layers.First,a deep neural network is adapted to learn from a large number of synthetic dispersion curves for inferring the layer number.Then,the shear-wave velocity profile is determined by a genetic algorithm with the known layer number.By applying this procedure to both simulated and real-world cases,the results indicate that the proposed method is reliable and efficient for surface wave inversion.

Effect of trigger system on experimental dispersion characteristics of active surface wave testing

A trigger system is typically employed in active seismic testing to trigger and synchronize multichannel surface wave data acquisition.The effect of the trigger system on the dispersion image of surface waves is empirically known to be negligible,however,theoretical explanation regarding the effect of the trigger system is insufficient.This study systematically examines the theory for surface wave dispersion analysis and proves that the effect of the trigger system on a dispersion image is negligible via a solid theoretical explanation.Subsequently,based on the new theoretical explanation,an alternative method that uses only the relative phase difference between sensors to extract dispersion characteristics with better conceptual clarity is proposed.Two active surface wave testing cases are considered to validate the theory and method.The results indicate that(1)an accurate trigger system is not necessary for surface wave data acquisition,and(2)it is unnecessary to assume that the impact point is the generation point of the surface waves for the experimental dispersion analysis.

Crustal and uppermost mantle structure of the northeastern Qinghai-Xizang Plateau from joint inversion of surface wave dispersions and receiver functions with P velocity constraints

Lithospheric structure beneath the northeastern Qinghai-Xizang Plateau is of vital significance for studying the geodynamic processes of crustal thickening and expansion of the Qinghai-Xizang Plateau. We conducted a joint inversion of receiver functions and surface wave dispersions with P-wave velocity constraints using data from the Chin Array Ⅱ temporary stations deployed across the Qinghai-Xizang Plateau. Prior to joint inversion, we applied the H-κ-c method(Li JT et al., 2019) to the receiver function data in order to correct for the back-azimuthal variations in the arrival times of Ps phases and crustal multiples caused by crustal anisotropy and dipping interfaces. High-resolution images of vS, crustal thickness, and vP/vSstructures in the Qinghai-Xizang Plateau were simultaneously derived from the joint inversion. The seismic images reveal that crustal thickness decreases outward from the Qinghai-Xizang Plateau. The stable interiors of the Ordos and Alxa blocks exhibited higher velocities and lower crustal vP/vSratios. While, lower velocities and higher vP/vSratios were observed beneath the Qilian Orogen and Songpan-Ganzi terrane(SPGZ), which are geologically active and mechanically weak, especially in the mid-lower crust.Delamination or thermal erosion of the lithosphere triggered by hot asthenospheric flow contributes to the observed uppermost mantle low-velocity zones(LVZs) in the SPGZ. The crustal thickness, vS, and vP/vSratios suggest that whole lithospheric shortening is a plausible mechanism for crustal thickening in the Qinghai-Xizang Plateau, supporting the idea of coupled lithospheric-scale deformation in this region.

Theoretical analysis of surface waves in piezoelectric medium with periodic shunting circuits

The investigations of surface waves in the piezoelectric medium bring out great possibility in designing smart surface acoustic wave(SAW)devices.It is important to study the dispersion properties and manipulation mechanism of surface waves in the semi-infinite piezoelectric medium connected with periodic arrangement of shunting circuits.In this study,the extended Stroh formalism is developed to theoretically analyze the dispersion relations of surface waves under different external circuits.The band structures of both the Rayleigh wave and the Bleustein-Gulyaev(BG)wave can be determined and manipulated with proper electrical boundary conditions.Furthermore,the electromechanical coupling effects on the band structures of surface waves are discussed to figure out the manipulation mechanism of adjusting electric circuit.The results indicate that the proposed method can explain the propagation behaviors of surface waves under the periodic electrical boundary conditions,and can provide an important theoretical guidance for designing novel SAW devices and exploring extensive applications in practice.

A speckle noise suppression method based on surface waves investigation and monitoring data

The internal energy distribution of waves can be described using ocean-wave spectra.In many ways,obtaining wave spectra on a global scale is critical.Surface waves investigation and monitoring onboard the Chinese-French oceanography satellite is the first space-borne instrument for detecting wave spectra specially,which was launched on October 29,2018.It can avoid the shortage of synthetic aperture radar detection results while still having some problems,especially with the effects of speckle noise.In this study,a method to suppress the speckle noise is proposed.First,the empirical formula for background speckle noise is established.Second,many spatio-temporal representative fluctuation spectra are classified and averaged.Third,rational transfer function filtering is used to obtain speckle noise close to the along-track direction.Finally,a signal-to-noise ratio threshold is used to suppress the abnormal speckle noise.This method solves the problems existing in previous denoising methods,such as excessive denoising in the along-track direction and the inability of some abnormal noises to be denoised in the two-dimensional directional wave spectra.

An ISGW Filtering Antenna with Spurious Modes and Surface Wave Suppression for Millimeter Wave Communications

In this paper,an integrated substrate gap waveguide(ISGW)filtering antenna is proposed at millimeter wave band,whose surface wave and spurious modes are simultaneously suppressed.A secondorder filtering response is obtained through a coupling feeding scheme using one uniform impedance resonator(UIR)and two stepped-impedance resonators(SIRs).To increase the stopband width of the antenna,the spurious modes are suppressed by selecting the appropriate sizes of the ISGW unit cell.Furthermore,the ISGW is implemented to improve the radiation performance of the antenna by alleviating the propagation of surface wave.And an equivalent circuit is investigated to reveal the working principle of ISGW.To demonstrate this methodology,an ISGW filtering antenna operating at a center frequency of 25 GHz is designed,fabricated,and measured.The results show that the antenna achieves a stopband width of 1.6f0(center frequency),an out-of-band suppression level of 21 dB,and a peak realized gain of 8.5 dBi.

Method of measuring one-dimensional photonic crystal period-structure-film thickness based on Bloch surface wave enhanced Goos–H?nchen shift

This paper puts forward a novel method of measuring the thin period-structure-film thickness based on the Bloch surface wave(BSW) enhanced Goos–Hanchen(GH) shift in one-dimensional photonic crystal(1DPC). The BSW phenomenon appearing in 1DPC enhances the GH shift generated in the attenuated total internal reflection structure. The GH shift is closely related to the thickness of the film which is composed of layer-structure of 1DPC. The GH shifts under multiple different incident light conditions will be obtained by varying the wavelength and angle of the measured light, and the thickness distribution of the entire structure of 1DPC is calculated by the particle swarm optimization(PSO) algorithm.The relationship between the structure of a 1DPC film composed of TiO_(2) and SiO_(2) layers and the GH shift, is investigated.Under the specific photonic crystal structure and incident conditions, a giant GH shift, 5.1 × 10^(3) times the wavelength of incidence, can be obtained theoretically. Simulation and calculation results show that the thickness of termination layer and periodic structure bilayer of 1DPC film with 0.1-nm resolution can be obtained by measuring the GH shifts. The exact structure of a 1DPC film is innovatively measured by the BSW-enhanced GH shift.

Three dimensional shear wave velocity structure of the crust and upper mantle beneath China from ambient noise surface wave tomography

We determine the three-dimensional shear wave velocity structure of the crust and upper mantle in China using Green＇s functions obtained from seismic ambient noise cross-correlation. The data we use are from the China National Seismic Network, global and regional networks and PASSCAL stations in the region. We first acquire cross-correlation seismograms between all possible station pairs. We then measure the Rayleigh wave group and phase dispersion curves using a frequency-time analysis method from 8 s to 60 s. After that, Rayleigh wave group and phase velocity dispersion maps on 1°by 1°spatial grids are obtained at different periods. Finally, we invert these maps for the 3-D shear wave velocity structure of the crust and upper mantle beneath China at each grid node. The inversion results show large-scale structures that correlate well with surface geology. Near the surface, velocities in major basins are anomalously slow, consistent with the thick sediments. East-west contrasts are striking in Moho depth. There is also a fast mid-to-lower crust and mantle lithosphere beneath the major basins surrounding the Tibetan plateau （TP） and Tianshan （Junggar, Tarim, Ordos, and Sichuan）. These strong blocks, therefore, appear to play an important role in confining the deformation of the TP and constraining its geometry to form its current triangular shape. In northwest TP in Qiangtang, slow anomalies extend from the crust to the mantle lithosphere. Meanwhile, widespread, a prominent low-velocity zone is observed in the middle crust beneath most of the central, eastern and southeastern Tibetan plateau, consistent with a weak （and perhaps mobile） middle crust.

Lithospheric structure and deformation in SE Tibet revealed by ambient noise and earthquake surface wave tomography: Recent advances and perspectives

High-resolution lithospheric structure is essential for understanding the tectonic evolution and deformation patterns of the southeastern Tibetan plateau. This is now possible due to recent advances in ambient noise and earthquake surface wave tomography, and great improvements in data coverage from dense portable array stations deployed in SE Tibet. In this review paper, I first give a brief overview of the tomographic methods from ambient noise and earthquake surface waves, and then summarize the major findings about the lithospheric structure and deformation in SE Tibet revealed by ambient noise and earthquake surface wave tomography as well as by other seismic and geophysical observations. These findings mainly include the 3-D distribution of mechanically weak zones in the mid-lower crust, lateral and vertical variations in radial and azimuthal anisotropy, possible interplay of some fault zones with crustal weak zones, and importance of strike-slip faulting on upper crustal deformation. These results suggest that integration of block extrusion in the more rigid upper-middle crust and channel flow in the more ductile mid-lower crust will be more compatible with the current geophysical observations. Finally I discuss some future perspective researches in SE Tibet, including array-based tomography, joint inversion using multiple seismic data, and integration of geodynamic modeling and seismic observations.

Shear velocity structure of crust and uppermost mantle in China from surface wave tomography using ambient noise and earthquake data

We present a 3D model of shear velocity of crust and upper mantle in China and surrounding regions from surface wave tomography. We combine dispersion measurements from ambient noise correlation and traditional earthquake data. The stations include the China National Seismic Network, global networks, and all the available PASSCAL stations in the region over the years. The combined data sets provide excellent data coverage of the region for surface wave measurements from 8 to 120 s, which are used to invert for 3D shear wave velocity structure of the crust and upper mantle down to about 150 kin. We also derive new models of the study region for crustal thickness and averaged S velocities for upper, mid, and lower crust and the uppermost mantle. The models provide a fundamental data set for understanding continental dynamics and evolution. The tomography results reveal significant features of crust and upper mantle structure, including major basins, Moho depth variation, mantle velocity contrast between eastern and western North China Craton, widespread low-velocity zone in mid- crust in much of the Tibetan Plateau, and clear velocity contrasts of the mantle lithosphere between north and southern Tibet with significant E-W variations. The low velocity structure in the upper mantle under north and eastern TP correlates with surface geological boundaries. A patch of high velocity anomaly is found under the eastern part of the TP, which may indicate intact mantle lithosphere. Mantle lithosphere change from the western to The Tanlu Fault appears boundary. shows striking systematic eastern North China Craton. to be a major lithosphere

Wave height measurement in the Taiwan Strait with a portable high frequency surface wave radar

As an important equipment for sea state remote sensing, high frequency surface wave radar （HFSWR） has received more and more attention. The conventional method for wave height inversion is based on the ratio of the integration of the second-order spectral continuum to that of the first-order region, where the strong external noise and the incorrect delineation of the first- and second-order Doppler spectral regions due to spectral aliasing are two major sources of errors in the wave height. To account for these factors, two more indices are introduced to the wave height estimation, i.e., the ratio of the maximum power of the second-or- der continuum to that of the Bragg spectral region （RSCB） and the ratio of the power of the second harmonic peak to that of the Bragg peak （RSHB）. Both indices also have a strong correlation with the underlying wave height. On the basis of all these indices an empirical model is proposed to estimate the wave height. This method has been used in a three-months long experiment of the ocean state measuring and analyzing ra- dar, type S （OSMAR-S）, which is a portable HFSWR with compact cross-loop/monopole receive antennas developed by Wuhan University since 2006. During the experiment in the Taiwan Strait, the significant wave height varied from 0 to 5 m. The significant wave heights estimated by the OSMAR-S correlate well with the data provided by the Oceanweather Inc. for comparison, with a correlation coefficient of 0.74 and a root mean square error （RMSE） of 0.77 m. The proposed method has made an effective improvement to the wave height estimation and thus a further step toward operational use of the OSMAR-S in the wave height extraction.

On the Interaction of Surface Waves with An Elastic Plate of Finite Length in Head Seas

An eigen-function expansion method based on a new orthogonal inner product is proposed by Sahoo et al. (2000) for the study of the hydroelastic response of mat-type VLFS in head seas. However, their main emphasis is on the effect of edge conditions and they assume that the plate is of a semi-infinite length. In reality, the plate is of finite length. For consideration of the finite length effect, the reflection and transmission from the other end must be considered. The effect of this reflection and transmission on the hydroelastic response of VLFS is of interest for practical application. Furthermore, the physical meaning of the new inner product was not given in their paper. In this paper, it is shown that the new inner product can he derived from the governing equation and the bottom boundary conditions. Then the same eigen-function expansion method is adopted for the study of the hydroelastic response of an elastic plate of finite length in surface waves. Detailed comparisons are made between the present finite length model and the semi-infinite model and between the present model predictions and the experimental results. It is found that that the finite length effect is significant and the accuracy of present model is higher than the semi-infinite model. Furthermore, a new phenomenon, which is not mentioned in Sahoo et al. (2000), is found. Taht is, for larger L/h ratios, the reflection and transmission coefficients will oscillate with the non-dimensional parameter k(0) h. Further study is needed for full understanding of this phenomenon.

Imaging shallow structure with active-source surface wave signal recorded by distributed acoustic sensing arrays

Distributed acoustic sensing(DAS) is one recently developed seismic acquisition technique that is based on fiber-optic sensing. DAS provides dense spatial spacing that is useful to image shallow structure with surface waves.To test the feasibility of DAS in shallow structure imaging,the PoroTomo team conducted a DAS experiment with the vibroseis truck T-Rex in Brady’s Hot Springs, Nevada, USA.The Rayleigh waves excited by the vertical mode of the vibroseis truck were analyzed with the Multichannel Analysis of Surface Waves(MASW) method. Phase velocities between5 and 20 Hz were successfully extracted for one segment of cable and were employed to build a shear-wave velocity model for the top 50 meters. The dispersion curves obtained with DAS agree well with the ones extracted from co-located geophones data and from the passive source Noise Correlation Functions(NCF). Comparing to the co-located geophone array, the higher sensor density that DAS arrays provides help reducing aliasing in dispersion analysis, and separating different surface wave modes. This study demonstrates the feasibility and advantage of DAS in imaging shallow structure with surface waves.

A numerical simulation of surface wave excitation in a rectangular planar-type plasma source

The principle of surface wave plasma discharge in a rectangular cavity is introduced simply based on surface plasmon polariton theory. The distribution of surface-wave electric field at the interface of the plasma-dielectric slab is investigated by using the three-dimensional finite-difference time-domain method （3D-FDTD） with different slotantenna structures. And the experimental image of discharge with a novel slot antenna array and the simulation of the electric field with this slot antenna array are both displayed. Combined with the distribution of surface wave excitation and experimental results, the numerical simulation performed by using 3D-FDTD is shown to be a useful tool in the computer-aided antenna design for large area planar-type surface-wave plasma sources.

Effects of Sea-Surface Waves and Ocean Spray on Air-Sea Momentum Fluxes

The effects of sea-surface waves and ocean spray on the marine atmospheric boundary layer （MABL） at different wind speeds and wave ages were investigated. An MABL model was developed that introduces a wave-induced component and spray force to the total surface stress. The theoretical model solution was determined assuming the eddy viscosity coefficient varied linearly with height above the sea surface. The wave-induced component was evaluated using a directional wave spectrum and growth rate. Spray force was described using interactions between ocean-spray droplets and wind-velocity shear. Wind profiles and sea-surface drag coefficients were calculated for low to high wind speeds for wind-generated sea at different wave ages to examine surface-wave and ocean-spray effects on MABL momentum distribution. The theoretical solutions were compared with model solutions neglecting wave-induced stress and/or spray stress. Surface waves strongly affected near-surface wind profiles and sea-surface drag coefficients at low to moderate wind speeds. Drag coefficients and near-surface wind speeds were lower for young than for old waves. At high wind speeds, ocean-spray droplets produced by wind-tearing breaking-wave crests affected the MABL strongly in comparison with surface waves, implying that wave age affects the MABL only negligibly. Low drag coefficients at high wind caused by ocean-spray production increased turbulent stress in the sea-spray generation layer, accelerating near-sea-surface wind. Comparing the analytical drag coefficient values with laboratory measurements and field observations indicated that surface waves and ocean spray significantly affect the MABL at different wind speeds and wave ages.

Smooth constraint inversion technique in genetic algorithms and its application to surface wave study in the Tibetan Plateau

Smooth constraint is important in linear inversion, but it is difficult to apply directly to model parameters in genetic algorithms. If the model parameters are smoothed in iteration, the diversity of models will be greatly suppressed and all the models in population will tend to equal in a few iterations, so the optimal solution meeting requirement can not be obtained. In this paper, an indirect smooth constraint technique is introduced to genetic inversion. In this method, the new models produced in iteration are smoothed, then used as theoretical models in calculation of misfit function, but in process of iteration only the original models are used in order to keep the diversity of models. The technique is effective in inversion of surface wave and receiver function. Using this technique, we invert the phase velocity of Raleigh wave in the Tibetan Plateau, revealing the horizontal variation of S wave velocity structure near the center of the Tibetan Plateau. The results show that the S wave velocity in the north is relatively lower than that in the south. For most paths there is a lower velocity zone with 12-25 km thick at the depth of 15-40 km. The lower velocity zone in upper mantle is located below the depth of 100 km, and the thickness is usually 40-80 km, but for a few paths reach to 100 km thick. Among the area of Ando, Maqi and Ushu stations, there is an obvious lower velocity zone with the lowest velocity of 4.2-4.3 km/s at the depth of 90-230 km. Based on the S wave velocity structures of different paths and former data, we infer that the subduction of the Indian Plate is delimited nearby the Yarlung Zangbo suture zone.

Vessel fusion tracking with a dual-frequency high-frequency surface wave radar and calibrated by an automatic identification system

High-frequency surface wave radar（HFSWR） and automatic identification system（AIS） are the two most important sensors used for vessel tracking.The HFSWR can be applied to tracking all vessels in a detection area,while the AIS is usually used to verify the information of cooperative vessels.Because of interference from sea clutter,employing single-frequency HFSWR for vessel tracking may obscure vessels located in the blind zones of Bragg peaks.Analyzing changes in the detection frequencies constitutes an effective method for addressing this deficiency.A solution consisting of vessel fusion tracking is proposed using dual-frequency HFSWR data calibrated by the AIS.Since different systematic biases exist between HFSWR frequency measurements and AIS measurements,AIS information is used to estimate and correct the HFSWR systematic biases at each frequency.First,AIS point measurements for cooperative vessels are associated with the HFSWR measurements using a JVC assignment algorithm.From the association results of the cooperative vessels,the systematic biases in the dualfrequency HFSWR data are estimated and corrected.Then,based on the corrected dual-frequency HFSWR data,the vessels are tracked using a dual-frequency fusion joint probabilistic data association（JPDA）-unscented Kalman filter（UKF） algorithm.Experimental results using real-life detection data show that the proposed method is efficient at tracking vessels in real time and can improve the tracking capability and accuracy compared with tracking processes involving single-frequency data.

Surface wave tomography on a large-scale seismic array combining ambient noise and teleseismic earthquake data

We discuss two array-based tomography methods, ambient noise tomography （ANT） and two-plane- wave earthquake tomography （TPWT）, which are capable of taking advantage of emerging large-scale broadband seismic arrays to generate high resolution phase velocity maps, but in complementary period band： ANT at 8-40 s and TPWT at 25-100 s period. Combining these two methods generates surface wave dispersion maps from 8 to 100 s periods, which can be used to construct a 3D vs model from the surface to -200 km depth. As an illustration, we apply the two methods to the USArray/Transportable Array. We process seismic noise data from over 1 500 stations obtained from 2005 through 2009 to produce Rayleigh wave phase velocity maps from 8 to 40 s period, and also perform TPWT using -450 teleseismic earthquakes to obtain phase velocity maps between 25 and 100 s period. Combining dispersion maps from ANT and TPWT, we construct a 3D vs model from the surface to a depth of 160 km in the western and central USA. These surface wave tomography methods can also be applied to other rapidly growing seismic networks such as those in China.

Surface tidal currents in the open sea area to the east of the Zhoushan Islands measured with high frequency surface wave radar

Based on the quasi-harmonic analysis of 11 d vector ocean currents obtained from two high frequency sur- face wave radars located at Zhujiajian Island and Shengshan Island, the spatial distribution characteristics of surface tidal currents in the open sea area to the east of the Zhoushan Islands of Zhejiang Province, China are studied. The following conclusions are drawn from the analysis： the tidal current pattern in the open sea area to the east of Zhoushan Islands is primarily regular semidiurnal, which is significantly affected by the shallow water constituents. The directions of the major axes of tidal current ellipses of M2 lie approx- imately in the NW-SE direction. With the increasing of distance away from the coast, the directions of the tidal current ellipses gradually shift toward the E-W direction. The tidal currents are mainly reversing cur- rents. The spatial distribution of probable maximum current velocities decreases gradually from northeast to southwest which is basically in accordance with the spatial distribution of measured maximum current velocities. The residual currents near the coast are larger than those far away from the coast. The directions of the residual currents are basically north by east, and the angle to the due north increases gradually with the increasing distance away from the coast. The topography shows a certain impact on the spatial distri- bution of shallow water constituents, the rotation of tidal currents, the probable maximum currents and the residual currents.

Noise Reduction Technique Applied to the Multichannel Analysis of Surface Waves

This paper shows the presence of noises and technique to reduce these noises during the surface wave analysis. The frequency-dependent properties of Rayleigh-type surface waves can be used for imaging and characterizing the shallow subsurface. Interference by coherent source-generated noise inhibits the reliability of shear-wave velocities determined through inversion of the phase velocities of Rayleigh waves. Among these interferences by non-planar, non-fundamental mode Rayleigh waves （noise） are body waves, scattered and non-source-generated surface waves, and highermode surface waves. For the reduction of noise, the filtering technique is implemented in this paper for the multichannel analysis of surface wave method （MASW）. With the de-noising technique during the MASW method, more robust and reliable outcome is achieved. The significance of this paper is to obtain pre-awareness about noises during surface wave analysis and take better outcomes with denoising performance in near surface soil investigations.