Phase Velocity Dispersion Curves Research Articles

Low temperature coefficient of frequency AlN Lamb wave resonator using groove structure between interdigital transducers

The lithographically tunable and small size features of Lamb wave resonators (LWRs) take a bright future for their use in high frequency narrow band filters. Resonant frequency drift due to temperature variation has a large impact on the narrower passband and a temperature coefficient of frequency (TCF) closer to 0 can broaden the stable temperature range of the resonator. This paper proposes a method of etching grooves in the area of the piezoelectric layer not covered by the Interdigital transducer (IDT) electrodes to improve the temperature stability of the Lamb resonator. Through the finite element method and theoretical analysis, the phase velocity, group velocity, and TCF dispersion curve of different modes of the resonator with groove structure were determined. The reason for the change of TCF under different normalized thicknesses of AlN was explained through the conversion of the LWR from a contour mode resonator (CMR) to a Cross-Sectional Lamé Mode Resonator. Simulation and test have shown that etching grooves have the effect of reducing TCF by adjusting the electromechanical coupling coefficient. The test shows that 205 nm-depth grooves can lower the TCF of the S0 mode IDT-Open LWR from −13.3 to −5.1 ppm/°C and S1 mode from −36.8 to −9.0 ppm/°C, respectively. The TCF of S0 and S1 mode LWRs decreased by 61.7% and 75.5%, respectively, after 205 nm groove etching. The groove etching method greatly raises the temperature stability of the LWR, enabling Lamb wave filters to operate stably over a wider temperature range.

3-D Shear wave velocity and azimuthal anisotropy structure in the shallow crust of Binchuan Basin in Yunnan, SW China, from ambient noise tomography

The Binchuan Basin in northwest Yunnan, SW China, is a rift basin developed at the intersection of the Red River Fault and Chenghai Fault, where active faults developed and historical earthquakes occurred. Understanding the fine velocity structure of the shallow crust in this region is crucial for comprehending earthquake occurrences and mechanisms and for improving earthquake location accuracy. Using the continuous waveform data recorded by 381 dense array stations in 2017, we obtained 7915 Rayleigh-wave phase velocity dispersion curves in the period band of 0.2-6 s from ambient noise cross-correlation functions after rigorous data processing and quality control. We determined 3-D isotropic and azimuthally anisotropic shear wave velocity models at depths above 6 km in the shallow crust based on the direct surface wave azimuthal anisotropic tomography method. The isotropic model reveals a strong correspondence between the S-wave velocity structure at depths of 0-1 km and the regional topography and lithology. The Binchuan depocenter, Zhoucheng depocenter, Xiangyun Basin, and Xihai Rift Basin are primarily composed of Quaternary deposits, which show low-velocity anomalies, while the regions with the Paleozoic shale, limestone, and basalt exhibit high-velocity anomalies. The nearly N-S orientation of fast directions from azimuthal anisotropy models are mainly controlled by the active Binchuan Fault with N-S strike as well as the NNW oriented primary compressive stress.

A new guided mode so-called minimum group velocity in viscoelastic sandwich plates: A parametric numerical study

Zero-Group-Velocity (ZGV) guided mode has recently received particular attention because of its highly sensitive to structural changes, such as defects. Parallel to ZGV, in the present work, we introduce a new guided mode so-called Minimum-Group-Velocity (MGV) that is a remarkable phenomenon in defect detection applications. Strictly speaking, when any two phase-velocity dispersion curves of Lamb modes approach each other without overlapping, we observe this new mode in group-velocity dispersion curves. It is worth noting that this MGV mode has never been addressed in literature. We compute the Lamb-like modes in asymmetric viscoelastic sandwich plates using the Stiffness Matrix Method (SMM). We show that the change in properties of this sandwich can dramatically affect frequency ranges associated with ZGV and MGV modes. However, a linear correlation between phase velocities and this change on viscoelastic properties is discussed, where we observed an interconnection phenomenon in the dispersion curves. This interconnection happens between specific points on the dispersion curves, each corresponding to a ZGV- or MGV-mode. On the other hand, we verify that the detection of changes in level of viscoelastic properties or variations in the thickness of the buffer viscoelastic layer is possible within a limited segment of the damped mode curve.

Retrieval and precise phase-velocity estimation of Rayleigh waves by the spatial autocorrelation method between distributed acoustic sensing and seismometer data

Summary In distributed acoustic sensing (DAS), optical fibre is used as sensors, which enables us to observe strain over tens of kilometres at intervals of several metres. S-wave velocity (Vs) structures of shallow sediments of high resolution have been obtained from surface wave dispersion curves by applying seismic interferometry to DAS data both onshore and offshore. However, it is known that there is a disadvantage to DAS seismic interferometry. In addition to Rayleigh waves, Love waves are also included. Consequently, the accuracy of the estimated phase velocities for Rayleigh waves is reduced due to the contamination of Love waves. To address this shortcoming, we suggest a spatial autocorrelation (SPAC) method between DAS and the vertical component of seismometer data. The SPAC method is equivalent to seismic interferometry and is useful for obtaining phase velocity dispersion curves of surface waves from the cross-correlation functions (CCFs) between the records of two receivers. The CCFs obtained from a combination of DAS and vertical seismometer data should contain only Rayleigh waves because Love waves have no vertical component. CCFs between DAS and vertical seismometer data are therefore expected to give more accurate phase velocities of Rayleigh waves than CCFs with DAS data only. In this study, we first formulated analytical expressions of cross-spectra for DAS and three-component seismometer data because seismic observation is generally carried out using a three-component seismometer. A new SPAC method is presented in the form of analytical expressions. We showed that our formulation only includes Rayleigh and not Love waves in the cross-spectra with DAS and the vertical-component seismometer data. We applied our SPAC method to actual DAS and vertical seismometer data recorded on the seafloor. Then, we compared our new SPAC method for DAS and vertical seismometer data with a conventional SPAC method for only DAS data. The results reveal that our new SPAC method can estimate the phase velocities of Rayleigh waves more accurately than the conventional method. In addition, the analytical formulations of the cross-spectrum between DAS and three-component seismometer data, which we obtained in this study, are expected to be useful for the estimation of accurate three-dimensional structures in the future, although this is not available at the moment due to the lack of an applicable dataset.

Evaluation of Robustness of S-Transform Based Phase Velocity Estimation in Viscoelastic Phantoms and Renal Transplants.

Ultrasound shear wave elastography (SWE) methods are being used to differentiate healthy versus diseased tissue on the basis of their viscoelastic mechanical properties. Tissue viscoelasticity is often studied by analyzing shear wave phase velocity dispersion curves, which is the variation of phase velocity with frequency or wavelength. Recently, a unique approach using a generalized Stockwell transformation (GST-SFK) was proposed for the calculation of dispersion curves in viscoelastic media over expanded frequency band. In this work, the method's robustness was evaluated on data from five custom-made viscoelastic tissue-mimicking phantoms and sixty in vivo renal transplants. For each phantom, 15 shear wave motion data acquisitions were taken, while 10-13 acquisitions were acquired for renal transplants measured in the renal cortex. For each data-set mean and standard deviation (SD) of estimated phase velocity dispersion curves were studied. In addition, the viscoelastic parameters of the Zener model were examined, which were preceded by a convergence analysis. For viscoelastic phantoms scanned with a research ultrasound scanner, and for the in vivo renal transplants scanned with a clinical scanner, the decisive advantage of the GST-SFK method over the standard two-dimensional Fourier transform (2D-FT) method was shown. The GST-SFK method provided dispersion curve estimates with lower SD over a wider frequency band in comparison to the 2D-FT method. These advantages are relevant to the analysis of the mechanical properties of tissues in clinical practice to discriminate healthy from diseased tissue.

Microtremor Full-Wavefield Modeling of Effective Phase Velocity and Horizontal-to-Vertical Spectral Ratio at Kyoto Reference Borehole Site: Comparison with Surface-Wavefield Modeling Based on a Velocity Structure with a Cap Layer

ABSTRACT Two types of data commonly used for microtremor exploration are phase-velocity dispersion curves obtained through an array measurement and horizontal-to-vertical spectral ratios (HVSRs) obtainable by a single-station measurement. Phase-velocity dispersion curves obtained by applying the spatial autocorrelation method to the array waveforms have a characteristic peaked shape in some cases. This dispersion curve shape has traditionally been explained as a consequence of the predominance of higher modes over fundamental mode in the Rayleigh waves. In this study, the effects of body waves on phase velocities and HVSRs were investigated based on both field measurements and theoretical calculations of microtremors. We used vertical-component array waveforms and single-station three-component waveforms of microtremors, obtained at and around a site where combined P-wave–S-wave (PS) and density loggings were conducted in the Kyoto basin, Japan (site KD-1), to identify phase velocities and HVSRs at frequencies in the range 0.2–2 Hz. The corresponding theoretical phase velocities and HVSRs were identified using full-wavefield synthetic data, which were generated assuming excitation points randomly distributed over the surface of a horizontally stratified velocity structure model created based on the logging data. The following key results were obtained. The measured phase-velocity dispersion curve exhibits a peaked shape with the value exceeding the S-wave velocity of the Tamba Group (Tb-Group), which is the bedrock (half-space) of the velocity structure model. Theoretical calculations based on the surface-wavefield theory were unable to reproduce this peaked shape; however, theoretical calculations based on the full-wavefield theory reproduced it with extraordinary accuracy. To reproduce the peaked shape based on the surface-wavefield theory, it was necessary to construct a model containing a cap (i.e., high-velocity layer) connected under the Tb-Group. The theoretical calculation based on the full wavefield also accurately reproduced the peak value and peak frequency of the measured HVSRs.

Method of automated extracting dispersion curves based on time-frequency distribution of seismic data

This paper discusses examples of testing a new implementation of the method of multi-channel surface wave analysis on synthetic data computed for elastic media with complex boundary geometry. The new implementation of the method includes developed algorithms for noise-resistant spectral analysis based on time-frequency domain filtering of seismograms and inversion of dispersion curves of phase velocities based on determination of ranges of possible transverse wave velocity models and application of artificial neural networks. Based on the results of synthetic data processing, the accuracy, lateral resolution limitations and applicability limits of the method under consideration are evaluated.

Potential geothermal reservoir systems in the Kenyan Great Rift Valley and volcanic region assessed by ambient noise analysis

To identify potential geothermal reservoirs in the Kenyan Great Rift, we analyzed ambient noise data to reveal seismic velocities and map azimuth-based Rayleigh-wave phase velocities. We extracted phase velocity dispersion curves from the ambient noise cross-correlation using the zero-crossing method and constructed S-wave velocity models using 3D surface-wave tomography. S-wave velocities were lower under Silali and Paka volcanoes (i.e., in magmatic regions) and in fractured areas (i.e., geothermal reservoirs) normal to the rift axis. We detected magmas at depths of >3 km beneath the volcanoes. Low-velocity anomalies may be related to fluid migration pathways. Most potential geothermal reservoirs are located above magma bodies and below the rift axis and are associated with intensive fracturing and porous rocks. Our approach used a limited number of seismometers, making it useful for application in areas with restricted access and challenging deployment; therefore, it could be used for the assessment of preliminary results in other magmatically active regions.

Detection of Low-Velocity Layer Using Passive Surface-Wave Data and Deep Learning: An Application in The Urban Area of Hangzhou, China

Passive surface-wave methods using dense seismic arrays have gained growing attention in near-surface high-resolution imaging in urban environments. Deep learning (DL) can release a tremendous workload brought by dense seismic arrays. We presented a case study of shear-wave velocity (Vs) structure imaging in the Hangzhou urban area (eastern China) using DL inversion. Noise data were recorded by dense linear arrays with approximately 5 m spacing deployed along two crossing roads for investigating the top 80 m of the subsurface. Phase-velocity dispersion curves are extracted from virtual shot gathers using multichannel analysis of surface waves. We divided the area where the low-velocity layer (LVL) may exist into three layers with a thickness of 5 m. We gave the four layers weak constraints to generate training dataset and adopted a convolutional neural network to directly invert fundamental-mode Rayleigh-wave phase velocity for 1D Vs models. To improve the accuracy, we further applied the sensitivities to weight the loss function in DL inversion. The obtained pseudo-2D Vs profiles correspond to the velocities estimated from logging data and previous survey. The well-trained neural network successfully identified that the LVL is located at 50-60 m deep. And this network was also achieved accurately the inversion of a dense seismic network nearby. The results of this survey demonstrate the accuracy and efficiency of delineating near-surface structures from traffic-induced noise using the DL technique, which has great potential for monitoring subsurface changes in urban areas.

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

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

Shear-wave velocity determination by combining data from passive and active source field investigations in Kumamoto city, Japan

We present the 1D subsoil structure and local site effects at KUMA strong ground motion station in Kumamoto City, Japan. We analyze data from a field campaign conducted in the framework of the Blind Prediction BP1 test of the 6th IASPEI/IAEE International Symposium: Effects of Surface Geology on Seismic Motion. In parallel with other participants of the BP1 test, we process data from passive and active source measurements aiming to determine the shear-wave velocity, Vs, structure and the site response at KUMA station. Passive measurements are associated to five microtremor arrays. In each array, seven seismometers have been deployed in a common-center triangle shape, recording microtremors simultaneously for 1 to 2 h. The vertical component of microtremors was analyzed using the spatial autocorrelation (SPAC) method. Cross-correlation coefficients were computed for all station pairs available for each array. By fitting the average SPAC’s coefficients to the first-kind zero-order Bessel function, J0, and assuming that microtremors primarily comprise fundamental mode Rayleigh waves, phase velocity dispersion curves were determined. Phase velocity values for frequencies > 15 Hz were obtained from data of a close-by active source geophone profile. We integrated the data with those of the passive measurements and obtained an experimental phase velocity dispersion curve. The resulting curve shows low velocity variation, from 150 to 200 m/s, in the surface layers, whereas significant dispersion appears in frequencies below 2.5 Hz. By inverting this curve, we achieved to determine the 1D shear-wave velocity structure at KUMA station. Site response characteristics were determined by applying the Horizontal-to-Vertical-Spectral-Ratios method. Significantly amplified peaks in the frequency range between 0.3 to 1.5 Hz dominate HVSR spectral ratios. These peaks correspond to resonant frequencies of soils and originate from different impedance contrasts within the substratum of the site.Graphical

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

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

Seismic structure of the Eastern European crust and upper mantle from probabilistic ambient noise tomography

The Eastern European lithosphere is a natural laboratory to study continental formation and evolution through time, comprising Archean continental remnants, Proterozoic rifts and belts, and younger accreted terranes. We investigate the seismic structure of the East European Craton (EEC) crust and uppermost mantle, and the transition from Precambrian to Phanerozoic Europe across the Trans European Suture Zone (TESZ) using probabilistic transdimensional ambient noise tomography. We cross-correlate noise recorded at broadband seismic stations from Eastern, Northern, and Central Europe, remove earthquake signals using continuous wavelet transform, and extract Rayleigh wave phase velocity dispersion curves. We invert these for the highest resolution shear wave velocity model of the Eastern European lithosphere to date, using Markov chain Monte Carlo Bayesian inversion. Our shear wave velocity model exhibits spatial correlation with major tectonic units and bears similarities with active seismic survey profiles in terms of seismic velocity patterns and main discontinuities. The crust thickens across the TESZ boundary and the mantle is seismically faster than beneath younger terranes, consistent with a less dense Precambrian lithosphere in the EEC. The crust and lithosphere beneath the Pannonian region is hyper-extended but the adjacent Transylvanian basin crust shows significant heterogeneity. The Precambrian building blocks of the EEC exhibit contrasting seismic fabrics. The Baltic orogens of Fennoscandia are underlain by uniform crust with a flat Moho, while Sarmatia shows alternating high and low velocity layers and a regional south-dipping crustal boundary from beneath the Ukrainian Shield towards the Crimean Peninsula. The last observation supports a geodynamic style driven by horizontal rather than vertical tectonics, with fundamental implications for the formation and evolution of early continents.

Determination of the contrast boundaries of the low velocity zone by analyzing the features of the dispersion curves of the Rayleigh wave velocities based on empirical dependencies

An important application of the multichannel surface wave (MASW) method is seismic safety assessment. The key parameters in the calculation of increments in the input of seismic microzoning are the thickness and average velocity of the soil mass lying on the base of more rigid rocks. To determine these parameters, we propose a new method for constructing horizontally layered models of the upper part of the geological section using the features (positions of extrema of the second derivative) of the dispersion curves of phase velocities of the Rayleigh surface wave, which, as shown by numerical experiments, are associated with the position of contrasting boundaries in the medium under study. (e.g., the boundary between soils and rocks). This approach is much simpler than the problem of recovering a horizontally layered model traditionally solved in the MASW method from a set of phase velocities for a sequence of frequencies and does not require an initial approximation and/or any restrictions on the possible values of the model parameters. In the case of two-layer and three-layer media, our approach is reduced to a simple and fast application of explicit formulas.

A comprehensive analysis of guided-wave propagation in 3D-printed PLA plates with different infill densities – Experimental study

Thanks to its many advantages, additive manufacturing (AM) is witnessing a gigantic shift from small applications to the production of complex industrial components, including safety–criticalapplications like aerospace and civil construction. This arises the need to develop accurate and robust technologies for evaluating and monitoring the structural integrity of AM components. The objective of this study is to comprehend how ultrasonic guided waves (GWs) propagate in 3D-printed poly(lactic acid) (PLA) plates of different infill densities. The experimental investigation involves five intact PLA plates printed with infill densities ranging from 20 to 100% at a step of 20%. A scanning laser Doppler vibrometer was used to visualize GW propagation by performing full-wavefield area-scans. Spatial-time, frequency-wavenumber, and wavenumber-wavenumber analyses were performed to scrutinize the availableGW modes and their slowness profiles within the anisotropic structure. In addition, the tuning of the fundamental Lamb-wave (LW) modes at different excitation frequencies was examined, and the amplitude attenuation of the generally most pronounced mode (A0 mode) was investigated for different infill densities at an excitation of 50 kHz. A reference methodology is presented to experimentallyidentify more than five propagating modes in the frequency-wavenumber domain then obtain their phase-velocity dispersion curves. Significant variations are observed in the dispersion characteristics when lowering the infill density leading to the appearance of more modes at lower cutoff frequencies and drops in both the group and phase velocities in most of the cases. A complex behavior of the guided waves was observed in the dispersion curves of the 20%-infill plate, with possible mode interruptions/disturbance at an excitation frequency beyond 70 kHz. The appearance of shear-horizontal modes is noticed/strengthened at oblique propagation directions which may be attributed to mode conversion of the LW modes after obliquely interacting with the internal structure of the plates. Though high wave attenuation is a characteristic of the used polymer material, this attenuation is further increased at lower infill densities. The tuning behavior of the fundamental LW modes is also shown to be significantly affected by the infill density. The presented results can serve as reference propagation and dispersion characteristics for researchers working on similar materials and geometries. Such knowledge and understanding can be particularly useful for GW-based structural healthy monitoring of AM structures.

Dispersion of Scholte wave under horizontally layered viscoelastic seabed

SUMMARY Dispersion inversion of Scholte wave is an effective method for constructing the shear wave velocity models of seabed sediments, but it is usually conducted based on the elastic layered medium theory, which ignores the viscoelasticity of sediments. In this work, we use the transitive matrix method to establish the dispersion equation for Scholte wave under horizontally layered viscoelastic seabed. This equation integrates the kinematic property of seismic wave in a viscoelastic media and the fluid–solid coupling mechanism. The phase velocity and attenuation coefficient dispersion curves of Scholte wave are presented by the real and imaginary parts of the complex-valued roots of the dispersion equation at different frequencies solved by Muller iteration algorithm, respectively. We perform numerical comparisons and analyses on the dispersion curves of Scholte waves for three typical seabed models under both elastic and viscoelastic conditions. Results demonstrate that the seabed viscoelasticity could greatly affect the propagation and dispersion characteristics of Scholte wave. Moreover, the dispersion curves of Scholte wave are sensitive to the variations in S-wave velocity and quality factor of seabed sediments.

Dispersion Curves of Love Waves in Elastic Waveguides Loaded with a Newtonian Liquid Layer of Finite Thickness

In this paper, the authors analyse the propagation of surface Love waves in an elastic layered waveguide (elastic guiding layer deposited on an elastic substrate) covered on its surface with a Newtonian liquid layer of finite thickness. By solving the equations of motion in the constituent regions (elastic substrate, elastic surface layer and Newtonian liquid) and imposing the appropriate boundary conditions, the authors established an analytical form of the complex dispersion equation for Love surface waves. Further, decomposition of the complex dispersion equation into its real and imaginary part, enabled for evaluation of the phase velocity and attenuation dispersion curves of the Love wave. Subsequently, the influence of the finite thickness of a Newtonian liquid on the dispersion curves was evaluated. Theoretical (numerical) analysis shows that when the thickness of the Newtonian liquid layer exceeds approximately four penetration depths 4δ of the wave in a Newtonian liquid, then this Newtonian liquid layer can be regarded as a semi-infinite half-space. The results obtained in this paper can be important in the design and optimization of ultrasonic Love wave sensors such as: biosensors, chemosensors and viscosity sensors. Love wave viscosity sensors can be used to assess the viscosity of various liquids, e.g. liquid polymers.

An ultrasonic scanning system for the inspection of composite stiffened panels from elastic constant identification via inversion of guided waves

This paper describes an ultrasonic scanning system for the inspection of stiffened composite panels used in modern aircraft construction based on elastic constant identification. Conventional ultrasonic scanning systems track individual features of the propagating waves (e.g. amplitude, arrival time, etc.) to detect damage by tracking the associated wave scattering. An ability to map the elastic constants of the composite part can greatly enhance the sensitivity of the scan to structural damage that can be detected and quantified by a corresponding reduction in the stiffness components. Elastic constants can be also correlated to residual strength. The proposed scanning system utilizes a “single-input-dual-output“ (SIDO) scheme whereby ultrasonic guided waves are excited by an impact and detected by two air-coupled ultrasonic sensors. The phase velocity dispersion curves of the waves are experimentally measured under a single excitation at each scanning line using a phase-spectrum technique. These curves are then inverted using a Semi-Analytical Finite Element (SAFE) method as the wave propagation forward model and Simulated Annealing as the optimization routine. This approach yields estimates of several elastic effective constants of the composite panel at each scanning line. These constants are then related to impact damage in stiffened stringer-skin composite panels.

Isotropic and radially anisotropic S-velocity structure beneath the Arabian plate inferred from surface wave tomography

We estimated Rayleigh- and Love-wave group-velocity dispersion curves using dense seismic networks in Saudi Arabia and nearby countries. We merged our data with group- and phase-velocity dispersion curves from previous studies to produce a new isotropic and radially anisotropic S-velocity model for the Arabian plate and nearby regions. Strong negative radial anisotropy (VSH < VSV) is observed in the crustal depth beneath the southern Red Sea, but weak positive radial anisotropy (VSH > VSV) and relatively thick lithosphere are found beneath the northern Red Sea, implying different tectonics between the two regions. We found low-velocity anomalies beneath Cenozoic volcanism in the Arabian Peninsula, but corresponding radial anisotropic perturbations are observed differently in spatiotemporal regimes depending on source origin. From the strong positive radial anisotropy at asthenospheric depth, which correlates with a circular pattern of azimuthal anisotropy from shear-wave splitting, we infer toroidal mantle flow around regions of thick lithosphere including an eastern part of the Arabian platform, the Mesopotamian Plain, and the Zagros belt. On the contrary, weak radial anisotropy is found beneath the thin lithosphere of the Arabian shield and in the thick lithosphere of the Mesopotamian, which may imply vertical flow due to asthenospheric upwelling and cratonic growth, respectively.

Shallow seismic structure around the Vrancea Seismic Zone from joint inversion of ambient noise H/V ratios and surface wave dispersion

The South-East Carpathians in Romania host the Vrancea Seismic Zone (VSZ), an intraplate seismic nest of frequent high magnitude earthquakes that cause unusually high damage over extended areas. Estimating the amplification of seismic signals strongly depends on the depth to the geophysical bedrock and the complex near surface geological structures. To improve seismic hazard estimates and ultimately mitigate seismic risk, we determine the frequency of resonance at 81 broadband seismic stations in southeast Romania using the horizontal-to-vertical spectral ratio (HVSR) of ambient seismic noise. We jointly invert HVSR functions with Rayleigh wave phase velocity dispersion curves using a Markov-chain Monte Carlo algorithm under the assumption of a diffuse field. Our results show excellent correlation between the fundamental frequency of resonance and local surface geology: stations located on foreland sedimentary basins show HVSR peaks in the range 0–1 Hz, while those located in orogenic nappes reach values as high as 10 Hz. Seismic bedrock topography mirrors the crystalline bedrock depth inferred from gravity modeling, but is consistently ∼2000 m shallower, corresponding to the interface between Paleozoic-Triassic to Jurassic-Cretaceous sedimentary levels on the Moesian Platform constrained by borehole data. The near surface seismic model shows anomalously low velocities in the Focsani forearc basin, adjacent to the VSZ, consistent with extremely low resonance frequencies and deep sedimentary layer successions (>5 km). In the back-arc of the Carpathian Orogen, at the contact with the Transylvanian Basin, average resonance frequencies are <2 Hz and inferred seismic bedrock depths vary laterally, likely affected by buried volcanic structures and diapiric salt domes. Our study presents the largest endeavor to study anomalous intensity patterns and local site effects in and around the VSZ and has fundamental implications for seismic risk evaluation and future civil engineering design planning.