Surface Wave Phase Velocity Research Articles

Super-large diameter slurry shield tunnel encountering boulder formation: Detection, treatment and data analysis

Boulders pose substantial challenges in coastal urban tunnel engineering. This study, for the first time, utilizes extensive closely spaced borehole data from 3486 drilling points obtained during actual tunnelling through boulder strata to evaluate the effectiveness of microtremor detection, analyze variations in operational parameters, and assess the efficiency of blasting treatment. Additionally, leveraging borehole data, a 3D reconstruction of the soil-boulder-bedrock stratum was developed, and the stratum content ratio at the excavation face was derived. A multivariate linear regression model was established to explore the relationship between boulder geology and construction parameters. The results demonstrated that microtremor surveys effectively identified boulder distributions by analyzing horizontal-to-vertical spectral ratio curves and differences in surface wave phase velocity. Field drilling confirmed the high accuracy of microtremor predictions, particularly in areas with dense or thick boulder layers. Post-blasting secondary drilling revealed boulder fragments smaller than 30 cm, meeting the operational requirements for shield cutter and crusher systems. The presence of boulders was found to increase thrust, torque, and overturning moments, as well as cause abnormal cutter wear. Blasting significantly reduced the standard deviation of operational data, demonstrating its effectiveness in stabilizing tunnelling performance. Regression analysis showed a strong correlation between geological conditions and tunnelling parameters, with results indicating that blasting operations effectively mitigate the adverse impacts of boulders on shield tunnelling performance. This study aims to provide valuable insights and references for future engineering projects encountering similar complex geological conditions.

Dipolar Huygens-Kerker radiation for surface waves.

Exotic dipolar radiation with zero light emission in one direction but maximal light emission in the opposite direction was envisioned by Huygens in 1690, and it could emerge in vacuum if the ratio between the source's electric and magnetic dipole moments fulfills the Kerker condition as revealed by Kerker in 1983. Due to its intricate connection with both the Huygens principle and Kerker condition, this radiation phenomenon is suggested to be termed as dipolar Huygens-Kerker radiation, and at this moment, the ratio is termed as the Huygens-Kerker ratio. However, the dipolar Huygens-Kerker radiation remains underexplored in non-vacuum matters, inside which the source locates, especially for surface waves. Here we find that the dipolar Huygens-Kerker radiation of surface waves in principle could occur in non-vacuum matters and is essentially featured with the same normalized radiation pattern, which is closely related to the inclination factor that appears in the Fresnel-Kirchhoff diffraction theory. Moreover, the corresponding Huygens-Kerker ratio is intrinsically determined by the phase velocity of excited surface waves. To be specific, the Huygens-Kerker ratio is proportional to the phase velocity for transverse-magnetic surface waves but becomes inversely proportional to the phase velocity for transverse-electric surface waves.

On the temperature sensitivity of near-surface seismic wave speeds: application to the Groningen region, the Netherlands

SUMMARY Subsurface temperature measurements play a crucial role, for instance, in optimizing geothermal power plants and monitoring heat-storage systems. Previous studies have demonstrated that time-lapse variations in temperature can be correlated with variations in seismic wave speeds, offering the potential for temperature monitoring via seismic surveys. However, an apparent discrepancy has emerged between field and laboratory experiments. Field studies predominantly report positive correlations between temperature and seismic wave speeds, while laboratory experiments often show anticorrelations. This inconsistency underscores the need for a more comprehensive, physics-based understanding of temperature-induced wave speed changes. In this study, we strive to bridge the gap between field and laboratory findings by examining several mechanisms governing temperature-induced seismic wave speed changes, namely the intrinsic temperature dependency of elastic parameters and thermally induced elasticity. We present a physics-based modelling approach to identify the primary mechanisms responsible for temperature-induced seismic wave speed changes. By considering several end-member models, we find that intrinsic temperature dependency of elastic parameters (negative correlation) compete with thermal pressure effects (positive correlation). The precise initial and boundary conditions and physical parameters of the system under consideration will determine the weight of both effects. Temperature-related dilatation does not seem to play an important role. We apply our approach to loosely consolidated sediments in the shallow subsurface of the Groningen region, where subsurface temperature fluctuations are driven by seasonal atmospheric temperature fluctuations roughly between −5 and $30\, ^\circ$C. For these models, we predict seasonal temperature-induced changes in body-wave speeds of up to 8 per cent in the first few metres of the subsurface, high-frequency (above 2 Hz) surface wave phase velocity variations in the range of 1–2 per cent, and relative changes in site amplification on the order of 4 per cent. These findings contribute to a more comprehensive understanding of the intricate relationship between temperature and near-surface seismic properties, offering insights for applications as subsurface temperature monitoring systems.

On the eigenvalues and eigendisplacement of the critical mode in horizontally layered media

Wave propagation in horizontally layered media is a classical problem in seismic-wave theory. In semi-infinite space, a nondispersive Rayleigh wave mode exists, and the eigendisplacement decays exponentially with depth. In a layered model with increasing layer velocity, the phase velocity of the Rayleigh wave varies between the S-wave velocity of the bottom half-space and that of the classical Rayleigh wave propagated in a supposed half-space formed by the parameters of the top layer. If the phase velocity is the same as the P- or S-wave velocity of the layer, which is called the critical mode or critical phase velocity of surface waves, the general solution of the wave equation is not a homogeneous (expressed by trigonometric functions) or inhomogeneous (expressed by exponential functions) plane wave, but one whose amplitude changes linearly with depth (expressed by a linear function). Theories based on a general solution containing only trigonometric or exponential functions do not apply to the critical mode, owing to the singularity at the critical phase velocity. In this study, based on the classical framework of generalized reflection and transmission coefficients, the propagation of surface waves in horizontally layered media was studied by introducing a solution for the linear function at the critical phase velocity. Therefore, the eigenvalues and eigenfunctions of the critical mode can be calculated by solving a singular problem. The eigendisplacement characteristics associated with the critical phase velocity were investigated for different layered models. In contrast to the normal mode, the eigendisplacement associated with the critical phase velocity exhibits different characteristics. If the phase velocity is equal to the S-wave velocity in the bottom half-space, the eigendisplacement remains constant with increasing depth.

Formation Model for Acoustic Characteristics of Solid Media with Ordered Fracturing

Introduction. The development of new structural materials and improvement of existing technologies for the production of new products on their basis lead to the emergence of new types of medium discontinuities. Therefore, the development of new models of discontinuities that take the previously ignored parameters into account seems to be relevant for the purposes of nondestructive testing and structural measurements. This concerns, e.g., the roughness of adjacent surfaces of microcrack ordered sets.Aim. Theoretical substantiation for the processes of elastic waves propagation through an elastic medium containing an ordered lattice of microcracks with boundary conditions in the linear slip approximation, modified by taking into account the parameters of micro-convexities of microcrack rough boundaries. Database formation for experimental studies aimed at determining the physical and mechanical characteristics of structural materials.Materials and methods. The acoustic characteristics of materials were determined based on the derivation and solutions of dispersion equations describing the formation and propagation of effective longitudinal, transverse, and Rayleigh surface elastic waves in elastic media with ordered cracking. Their values were also used to determine the effective speed of Rayleigh surface waves.Results. The conducted simulation of elastic wave formation processes showed that an increase in the concentration of microcracks leads to a decrease in the phase velocities of effective longitudinal, transverse, and surface waves, as well as to an increase in the attenuation coefficients at given ultrasound frequencies and material parameters.Conclusion. The radius of the microsphere that replaces the surface micro-convexity and the roughness parameter Rz have a significant impact on the formation of physical and mechanical characteristics of materials, which are determined by the results of ultrasonic measurements. The developed model can be recommended as a basis for interpreting the results of ultrasonic measurements.

Bayesian eikonal tomography using Gaussian processes


 Eikonal tomography has become a popular methodology for deriving phase velocity maps from surface wave phase delay measurements. Its high efficiency makes it popular for handling datasets deriving from large-N arrays, in particular in the ambient-noise tomography setting. However, the results of eikonal tomography are crucially dependent on the way in which phase delay measurements are predicted from data, a point which has not been thoroughly investigated. In this work, I provide a rigorous formulation for eikonal tomography using Gaussian processes (GPs) to smooth observed phase delay measurements, including uncertainties. GPs allow the posterior phase delay gradient to be analytically derived. From the phase delay gradient, an excellent approximate solution for phase velocities can be obtained using the saddlepoint method. The result is a fully Bayesian result for phase velocities of surface waves, incorporating the nonlinear wavefront bending inherent in eikonal tomography, with no sampling required. The results of this analysis imply that the uncertainties reported for eikonal tomography are often underestimated.

Surface-Wave Imaging with Nonrandom Traffic Noise Sources

We analyze the influence of nonrandom traffic noise sources on surface wave imaging. With nonrandom traffic noise sources in time (correlated sources), spurious signals are produced in the crosscorrelation functions. With nonrandom traffic noise sources in space (directional sources), phase velocities are overestimated from the retrieved surface waves. The retrieval of surface waves can be improved by stationary-phase segment selection, cross-coherence and stacking of the results from multiple traffic noise sources. Synthetic results demonstrate that cross-coherence and stacking of multiple traffic noise sources makes the source tend to be random distribution in time and suppresses the spurious signals in the crosscorrelation functions. Enhancing the stationary-phase contribution makes more accurate estimation of surface wave phase velocities.

Six‐Component Earthquake Synchronous Observations Across Taiwan Strait: Phase Velocity and Source Location

AbstractWith the development of rotational seismometers, especially the high‐sensitivity optical seismometers, rotational motions have shown increasing benefits for calculation of dispersion curves, estimation of back azimuth and many other applications. Here we present the six‐component observations of 2019 Hualien mww6.1 earthquake simultaneously on both sides of Taiwan Strait. Seismograms and time‐frequency spectrums show that the attenuations for rotations are more significant than translations. Calculations of the surface wave phase velocity are performed for both stations through single‐station records, indicating a distinct difference owning to the instrument installation conditions and site effect. Comparison of theoretical dispersion curves is consistently in trends with estimations, and reveals the influence of higher‐mode surface waves. For the first time, we locate the source of a earthquake through intersection of two stations' back azimuths estimated through rotational methods. Results of comparing to translational polarization verify that rotational observation can distinctly reduce the deviation of source location.

Influence of orientation of fibres on surface wave in fibre‐induced anisotropic elastic medium

AbstractAn isotropic elastic solid reinforced with a family of parallel fibres behaves anisotropic to any deformation. With orientation of all fibres along a coordinate axis in rectangular Cartesian system, this anisotropy reduces to transverse isotropy. Else with arbitrary orientation of fibre‐family, induced anisotropy lacks any symmetry. In this composite medium with arbitrary anisotropy, three‐dimensional wave motion decaying with depth propagates as a surface wave along all directions on plane boundary. Real phase velocity for this surface wave lies implicit in a system of three equations with complex coefficients. An appropriate transformation of this system yields a real characteristic equation, which is solved to get the anisotropic phase velocity of generalised surface wave. This phase velocity defines a complex slowness vector, which is used to calculate the particle motion at any point in the medium. Effects of fibre‐orientation are analysed, numerically, on the phase velocity and polarisation of surface wave. For example, in the reinforced fibres alters the contours of particle motion through increased ellipticity. However, in the presence of fibres, the decay of particle motion with depth slows down.

Determining subsurface temperature & lithospheric structure from joint geophysical-petrological inversion: A case study from Ireland

High quality maps of the geothermal gradient and temperature are essential when assessing geothermal potential for a region. However, determining geothermal potential is a challenge as direct measurements of in situ temperature are sparse and individual geophysical methods are sensitive to a range of parameters, not solely temperature. Here, we develop a novel approach to determine the geothermal gradient using a new joint geophysical-petrological inversion where seismic velocities and density in the mantle are related to temperature and bulk composition within a thermodynamic framework. Large datasets of phase velocities of seismic surface-waves are now incorporated into the inversion, and provide essential constraints on the lithospheric thickness and temperature, which shape the crustal geotherms to a significant extent. We also include all available measurements of the surface heat flow, radiogenic heat production (RHP) and thermal conductivity within the crust, to further constrain the temperature and geothermal gradient, in particular in the top few kilometres of the crust. We use Ireland as a case study and show how our new methodology can reproduce the results of previous work but also improve on them, thanks to the complementary sensitivities of the full range of data. Lithospheric and crustal thicknesses have primary controls on the temperature gradient, with areas of thinner lithosphere showing higher geothermal gradients. In some locations, anomalously warm geotherms result from high RHP within crustal granitic rocks, mudstone and shales. RHP is above continental averages across all Ireland, likely due to a crust with mostly felsic lithology. The new methodology provides a robust workflow for determining the geothermal potential in areas with limited direct temperature measurements, facilitating knowledge creation for the transition to sustainable energy sources and energy self-sufficiency.

Simulation of ground motion emissions from wind turbines in low mountain ranges: implications for amplitude decay prediction

The characterization and prediction of wind turbine (WT) emissions are important steps in reducing their impact on humans or sensitive technologies such as seismic stations or physics experiments. Here, WT ground motion emissions are studied along two measurement lines set up at two wind farms on the Eastern Swabian Alb, southwest Germany. The main purpose of the data analysis is to estimate amplitude decay rates from vertical component data and surface wave phase velocities excited by the permanent motion of the WT towers. Phase velocities as well as geological information serve as input to build realistic subsurface models for numerical wave field simulations. Amplitude A decay rates are characterized by b-values through Asim 1/r^b depending on distance r and are derived from peaks in power spectral density (PSD). We find an increase of b_text {PSD} with frequency from 0.5 to 3.2 for field data. For low frequencies (1.2 Hz and 3.6 Hz), b_text {PSD} ranges from 0.5 to 1.1, hence close to the geometrical spreading factor of surface waves (b_text {PSD}=1). Anelastic damping and scattering seem not to be significant at these frequencies which also shows in numerical simulations for quality factors Q=50-200. We also find that the emitted wavefields from several WTs interfere, especially in the near-field, and produce strong local ground motion amplitudes. The inclusion of a steep topography present in low mountain ranges adds more wave field distortions which can further increase the amplitudes. This needs to be considered when predicting WT induced ground motions.

On the Accuracy of Surface-Wave Dispersion Measurements from Horizontal-Component Ambient Noise Cross Correlations

Abstract At present, both Rayleigh and Love waves can be extracted from cross-correlations of ambient noise in ambient noise tomography (ANT). It has been demonstrated that accurate Rayleigh-wave dispersion curves can be extracted from vertical-component noise when noise sources are evenly distributed in space. However, because the rotation of horizontal components alters the distributions of effective sources in horizontal–horizontal cross correlations between a pair of stations, it is still not completely sure if dispersion curves of Rayleigh and Love waves extracted from horizontal-component noise data are accurate for ANT. In this study, through numerical simulations, we demonstrate that surface-wave phase velocities extracted from the horizontal–horizontal cross-correlation functions (CCFs) systematically deviate from structural phase velocity even when noise sources are evenly distributed in space. The deviations of phase velocities quickly increase with decreasing station separations of CCFs when station separations are shorter than three wavelengths. Further analysis on field data recorded at seismic stations in the contiguous United States confirms the existence of the phase-velocity deviations of Rayleigh and Love waves when they are measured from horizontal–horizontal CCFs. Because the deviations become smaller with increasing station separations and are less than 0.25% at station separations longer than three wavelengths, we suggest that it is best to only select those CCFs with station separations longer than three wavelengths in ANT when measuring phase velocities of surface waves from horizontal-component CCFs, such as Love waves from radial–radial CCFs.

Generating 1D depth profile of shear-wave velocity from ambient noise cross-correlation: application to Jakarta array

In the past decade, cross-correlations of ambient seismic noise have been exploited in various applications to model the shallow-to-deep structure of Earth’s interior through tomographic inversions. The stack of cross-correlations between a 2-station pair represents empirical Green’s function and comprises the information of the subsurface structure between those stations. In practice, noise correlation function (NCF) is analyzed to reconstruct surface wave group or phase velocity dispersion; then, the dispersion data is used to model shear-wave velocity (Vs). This study presents a case for temporary seismic networks deployed in the Jakarta Basin; we applied a two-step routine to obtain a representative 1D Vs profile beneath an array. First, we extracted our array’s average phase velocity dispersion based on the relationship between NCF’s spectra and the Bessel function. Then, we invert for the 1D depth profile of Vs using a transdimensional Bayesian inversion to allow for exploring a number of layers in parameterizations. We successfully generate a 1D Vs profile up to 5 km depth reflecting the regional stratigraphy of the Jakarta Basin. In general, a sedimentary basin fill covers the area reaching a depth of 650 m. We suggest that this simple routine can be undertaken for other ambient noise cross-correlation cases; such a 1D depth profile would be beneficial to be used as a reference model.

Early Results of Surface Tomography Using Automated Surface-Wave Phase-Velocity Measuring System (ASWMS) in Hawaii

The Hawaiian Islands are a prominent place to conduct seismicity research. Geologically, occurring earthquakes in this region are unique since it is caused by the presence of a hot spot rather than the common fault activity. The research was conducted using a tomography method called Surface Wave Tomography. The goals of this study are to produce tomographic images of the Hawaiian Islands and to test the accuracy of the Automated Surface-Wave Phase-Velocity Measuring System (ASWMS). ASWMS employs a cross-correlation process to calculate the phase delay between stations, then Eikonal equation to invert the slowness vector to obtain apparent phase velocities, and lastly Helmholtz equation to correct the amplitude to obtain structural phase velocities. The teleseismic data used in this study were collected from stations located throughout the Hawaiian Islands between 2004 and 2009. This study discovered a low-velocity anomaly beneath the Hawaiian Islands with a value of 3.8 – 3.9 km/s, indicating the presence of a low-velocity body. The decrease in velocity with increasing depth suggests an increase in temperature associated with the path that magma takes from the Earth’s mantle to the Earth’s crust beneath the Hawaiian Islands.

Data Processing of Surface Wave Tomography in Australia using Automated Surface Wave Phase-Velocity Measuring System (ASWMS)

The Australian continent is divided into tectonic blocks, including Archaean and Proterozoic units in the western and central part, and Phanerozoic units in the eastern part. Surface wave tomography is an effective method of getting information about velocity variations in a particular region. Surface wave tomography was used in this project to map the subsurface of the Australian continent down to a depth of 100 km with the ASWMS (Automated Surface-Wave Phase-Velocity Measuring System). From 2010 to 2015, we gathered waveform data from 767 events and 219 permanent and temporal seismic stations dispersed across Australia. We selected a set of events that has a minimum magnitude of 6. Rayleigh wave extraction and cross-correlation between two nearby stations are among the first processing processes. Following that, the apparent phase-velocity value is calculated using the Eikonal and tomographic inversion. To generate a tomographic map with a true structural phase velocity, each phase velocity data is stacked with amplitude correction. According to current tomographic imaging studies, low-velocity values under Australia tend to be dispersed easterly, with Rayleigh velocity ranging from 3.7 to 3.9 km/s. The tectonic history and evolution of eastern Australia are related to a sequence of orogenic events that are pushing toward the continent’s eastern boundary.

Wave Dispersion Behavior in Quasi-Solid State Concrete Hydration.

This paper aims to investigate wave dispersion behavior in the quasi-solid state of concrete to better understand microstructure hydration interactions. The quasi-solid state refers to the consistency of the mixture between the initial liquid-solid stage and the hardened stage, where the concrete has not yet fully solidified but still exhibits viscous behavior. The study seeks to enable a more accurate evaluation of the optimal time for the quasi-liquid product of concrete using both contact and noncontact sensors, as current set time measurement approaches based on group velocity may not provide a comprehensive understanding of the hydration phenomenon. To achieve this goal, the wave dispersion behavior of P-wave and surface wave with transducers and sensors is studied. The dispersion behavior with different concrete mixtures and the phase velocity comparison of dispersion behavior are investigated. The analytical solutions are used to validate the measured data. The laboratory test specimen with w/c = 0.5 was subjected to an impulse in a frequency range of 40 kHz to 150 kHz. The results demonstrate that the P-wave results exhibit well-fitted waveform trends with analytical solutions, showing a maximum phase velocity when the impulse frequency is at 50 kHz. The surface wave phase velocity shows distinct patterns at different scanning times, which is attributed to the effect of the microstructure on the wave dispersion behavior. This investigation delivers profound knowledge of hydration and quality control in the quasi-solid state of concrete with wave dispersion behavior, providing a new approach for determining the optimal time of the quasi-liquid product. The criteria and methods developed in this paper can be applied to optimal timing for additive manufacturing of concrete material for 3D printers by utilizing sensors.

Structure and evolution of the Australian plate and underlying upper mantle from waveform tomography with massive data sets

SUMMARY We present a new S-wave velocity tomographic model of the upper mantle beneath the Australian Plate and its boundaries that we call Aus22. It includes azimuthal anisotropy and was constrained by waveforms from 0.9 million vertical-component seismograms, with the densest data sampling in the hemisphere centred on the Australian continent, using all available data covering this hemisphere. Waveform inversion extracted structural information from surface waves, S- and multiple S-waves and constrained S- and P-wave speeds and S-wave azimuthal anisotropy of the crust and upper mantle, down to the 660-km discontinuity. The model was validated by resolution tests and, for particular locations in Australia with notable differences from previous models, independent inter-station measurements of surface-wave phase velocities. Aus22 can be used to constrain the structure and evolution of the Australian Plate and its boundaries in fine detail at the regional scale. Thick, high-velocity (and, by inference, cold) cratonic lithosphere occupies nearly all of western and central Australia but shows substantial lateral heterogeneity. It extends up to the northern edge of the plate, where it collides with island arcs, without subducting. Diamondiferous kimberlites and lamproite deposits are underlain by cratonic lithosphere, except for the most recent diamondiferous lamproites in the King Leopold Orogen. The rugged eastern boundary of the cratonic lithosphere resolved by the model provides a lithospheric definition of the Tasman Line. Just east of the Tasman Line, an area of intermediate-thick lithosphere is observed in the southern part of the continent. The eastern part of Australia is underlain by thin, warm lithosphere, evidenced by low seismic velocities. All the sites of Cenozoic intraplate volcanism in eastern Australia are located on thin lithosphere. A low-velocity anomaly is present in the mantle transition zone (410–660 km depths) beneath the Lord Howe and Tasmanid hotspots, indicative of anomalously high temperature and consistent with a deep mantle upwelling feeding these hotspots and, possibly, also the East Australia hotspot. High seismic velocities at 200–410 km depth below New Guinea indicate the presence of slab fragments, probably linked to the subduction of the Australian Plate. High seismic velocities are observed in the transition zone below northeast Australia and indicate the presence of subducted lithospheric fragments trapped in the transition zone, possibly parts of the former northern continental margin of Australia.

Seismic source mapping by surface wave time reversal: application to the great 2004 Sumatra earthquake

SUMMARY Different approaches to map seismic rupture in space and time often lead to incoherent results for the same event. Building on earlier work by our team, we ‘time-reverse’ and ‘backpropagate’ seismic surface wave recordings to study the focusing of the time-reversed field at the seismic source. Currently used source-imaging methods relying on seismic recordings neglect the information carried by surface waves, and mostly focus on the P-wave arrival alone. Our new method combines seismic time reversal approach with a surface wave ray-tracing algorithm based on a generalized spherical-harmonic parametrization of surface wave phase velocity, accounting for azimuthal anisotropy. It is applied to surface wave signal filtered within narrow-frequency bands, so that the inherently 3-D problem of simulating surface wave propagation is separated into a suite of 2-D problems, each of relatively limited computational cost. We validate our method through a number of synthetic tests, then apply it to the great 2004 Sumatra–Andaman earthquake, characterized by the extremely large extent of the ruptured fault. Many studies have estimated its rupture characteristics from seismological data (e.g. Lomax, Ni et al., Guilbert et al., Ishii et al., Krüger & Ohrnberger, Jaffe et al.) and geodetic data (e.g. Banerjee et al., Catherine et al., Vigny et al., Hashimoto et al., Bletery et al.). Applying our technique to recordings from only 89 stations of the Global Seismographic Network (GSN) and bandpass filtering the corresponding surface wave signal around 80-to-120, 50-to-110 and 40-to-90 s, we reproduce the findings of earlier studies, including in particular the northward direction of rupture propagation, its approximate spatial extent and duration, and the locations of the areas where most energy appears to be released.

Study on the Elimination Method of Wind Field Influence in Retrieving a Sea Surface Current Field.

An along-the-track interferometric synthetic aperture radar (ATI-SAR) system can estimate the radial velocity of a moving target on the ground and on a sea surface current. This acquires the interference phase by combining two composite SAR images obtained by two antennas spatially separated along the direction of movement of the platform. The key to retrieving the sea surface current is to remove the interference of sea surface waves, wind-generated current, and Bragg phase velocity in the interference Doppler velocity. Previous methods removed the surface waves, Bragg phase velocity, and other interferences based on externally-assisted wind fields (e.g., ECMWF), using the M4S or other models. However, the wind fields obtained from ECMWF and other external information are often average results of a large temporal and spatial scale, while the images obtained from SAR are high-resolution images of sea surface transients, which are quite different in time and space. This paper takes the SAR image data of the Gaofen-3 satellite as the research object and employs an SAR-based wind field retrieval method to obtain an SAR-observed transient wind field. Combined with the CDOP model, the interference of Doppler velocities, such as the sea surface wave, wind-generated current, and Bragg wave phase velocity, was calculated and subtracted from the Doppler velocity, to obtain the sea surface velocity result. Then, the current field measured by the shore-based HF radar was compared with that obtained by correcting the ATI Doppler velocity based on the SAR retrieved wind field and the ECMWF wind field. The comparison of results indicated that the wind field correction result based on the SAR retrieved wind field was closer to the current field measured by the shore-based HF radar than the wind field correction result based on the ECMWF wind field.

Application of six-component ambient seismic noise data for high-resolution imaging of lateral heterogeneities

SUMMARY We develop a novel approach for imaging subsurface lateral heterogeneities using six-component (6C) ambient seismic noise data, consisting of three translational components and three rotational components. We first derive the 6C cross-correlation functions (CCFs) from ambient seismic noise data for surface waves and then apply the modified frequency Bessel (MF-J) transform to the 6C CCFs for surface wave dispersion curves. The phase velocities of Rayleigh and Love waves can be directly calculated from the 6C CCFs of a single pair of stations. Traditional array-based surface wave methods derive the surface wave phase velocity based on the presumption of a horizontally layered model. When this assumption breaks down due to the presence of, for example, dipping layer or heterogeneities, the resulting phase velocity can be severely smeared. By contrast, the proposed approach is not limited by the layered model assumption because it relies on single-point measurements to calculate the localized dispersion relations of the formation right beneath a receiver. Our numerical modelling results demonstrate that this approach is applicable to heterogeneous models and can image small-scale subsurface anomalies with very high lateral resolution.