To see the other types of publications on this topic, follow the link: Rayleigh`s phase function.
Journal articles on the topic 'Rayleigh`s phase function'
Create a spot-on reference in APA, MLA, Chicago, Harvard, and other styles
Consult the top 50 journal articles for your research on the topic 'Rayleigh`s phase function.'
Next to every source in the list of references, there is an 'Add to bibliography' button. Press on it, and we will generate automatically the bibliographic reference to the chosen work in the citation style you need: APA, MLA, Harvard, Chicago, Vancouver, etc.
You can also download the full text of the academic publication as pdf and read online its abstract whenever available in the metadata.
Browse journal articles on a wide variety of disciplines and organise your bibliography correctly.
1
Koroleva, Tat’iana. "Rayleigh wave velocity maps beneath the Caucasus from the Caucasus seismic network (CNET)." Russian Journal of Seismology 2, no. 3 (September 30, 2020): 70–77. http://dx.doi.org/10.35540/2686-7907.2020.3.06.
Full textAbstract:
Ambient noise surface wave tomography is a widely used method for determining the velocity structure of the upper layers of the Earth. It is based on the fact that the cross-correlation function (CCF) of noise at two stations, averaged over a long time interval, determines the Green's function of the surface wave. This allows us to estimate the group and phase velocities of surface waves on the paths between stations. The method was applied to the records of the vertical components of 67 seismic stations of the Caucasian network CNET network, which were obtained during 2018. The cross-correlation functions for all interstation paths were calculated. The dispersion curves of group and phase Rayleigh velocities for periods from 5 to 30 s were obtained by means of frequency-time analysis. The lateral distribution of the velocities was received for periods from 7 to 22 s, which correlate with velocity structure at depths of 5-25 km. The group and phase velocity maps for Rayleigh wave for periods 7, 10, 12, 14, 17, 22 s are presented.
2
Forbriger, Thomas, Lingli Gao, Peter Malischewsky, Matthias Ohrnberger, and Yudi Pan. "A single Rayleigh mode may exist with multiple values of phase-velocity at one frequency." Geophysical Journal International 222, no. 1 (March 17, 2020): 582–94. http://dx.doi.org/10.1093/gji/ggaa123.
Full textAbstract:
SUMMARY Other than commonly assumed in seismology, the phase velocity of Rayleigh waves is not necessarily a single-valued function of frequency. In fact, a single Rayleigh mode can exist with three different values of phase velocity at one frequency. We demonstrate this for the first higher mode on a realistic shallow seismic structure of a homogeneous layer of unconsolidated sediments on top of a half-space of solid rock (LOH). In the case of LOH a significant contrast to the half-space is required to produce the phenomenon. In a simpler structure of a homogeneous layer with fixed (rigid) bottom (LFB) the phenomenon exists for values of Poisson’s ratio between 0.19 and 0.5 and is most pronounced for P-wave velocity being three times S-wave velocity (Poisson’s ratio of 0.4375). A pavement-like structure (PAV) of two layers on top of a half-space produces the multivaluedness for the fundamental mode. Programs for the computation of synthetic dispersion curves are prone to trouble in such cases. Many of them use mode-follower algorithms which loose track of the dispersion curve and miss the multivalued section. We show results for well established programs. Their inability to properly handle these cases might be one reason why the phenomenon of multivaluedness went unnoticed in seismological Rayleigh wave research for so long. For the very same reason methods of dispersion analysis must fail if they imply wave number kl(ω) for the lth Rayleigh mode to be a single-valued function of frequency ω. This applies in particular to deconvolution methods like phase-matched filters. We demonstrate that a slant-stack analysis fails in the multivalued section, while a Fourier–Bessel transformation captures the complete Rayleigh-wave signal. Waves of finite bandwidth in the multivalued section propagate with positive group-velocity and negative phase-velocity. Their eigenfunctions appear conventional and contain no conspicuous feature.
3
Burns, P., and E. Meiburg. "Sediment-laden fresh water above salt water: nonlinear simulations." Journal of Fluid Mechanics 762 (November 27, 2014): 156–95. http://dx.doi.org/10.1017/jfm.2014.645.
Full textAbstract:
AbstractWhen a layer of particle-laden fresh water is placed above clear, saline water, both double-diffusive and Rayleigh–Taylor instabilities may arise. The present investigation extends the linear stability analysis of Burns & Meiburg (J. Fluid Mech., vol. 691, 2012, pp. 279–314) into the nonlinear regime, by means of two- and three-dimensional direct numerical simulations (DNS). The initial instability growth in the DNS is seen to be consistent with the dominant modes predicted by the linear stability analysis. The subsequent vigorous growth of individual fingers gives rise to a secondary instability, and eventually to the formation of intense plumes that become detached from the interfacial region. The simulations show that the presence of particles with a Stokes settling velocity modifies the traditional double-diffusive fingering by creating an unstable ‘nose region’ in the horizontally averaged profiles, located between the upward-moving salinity and the downward-moving sediment interface. The effective thickness $l_{s}$ ($l_{c}$) of the salinity (sediment) interface grows diffusively, as does the height $H$ of the nose region. The ratio $H/l_{s}$ initially grows and then plateaus, at a value that is determined by the balance between the flux of sediment into the rose region from above, the double-diffusive/Rayleigh–Taylor flux out of the nose region below, and the rate of sediment accumulation within the nose region. For small values of $H/l_{s}\leqslant O(0.1)$, double-diffusive fingering dominates, while for larger values $H/l_{s}\geqslant O(0.1)$ the sediment and salinity interfaces become increasingly separated in space and the dominant instability mode becomes Rayleigh–Taylor like. A scaling analysis based on the results of a parametric study indicates that $H/l_{s}$ is a linear function of a single dimensionless grouping that can be interpreted as the ratio of inflow and outflow of sediment into the nose region. The simulation results furthermore indicate that double-diffusive and Rayleigh–Taylor instability mechanisms cause the effective settling velocity of the sediment to scale with the overall buoyancy velocity of the system, which can be orders of magnitude larger than the Stokes settling velocity. While the power spectra of double-diffusive and Rayleigh–Taylor-dominated flows are qualitatively similar, the difference between flows dominated by fingering and leaking is clearly seen when analysing the spectral phase shift. For leaking-dominated flows a phase-locking mechanism is observed, which intensifies with time. Hence, the leaking mode can be interpreted as a fingering mode which has become phase-locked due to large-scale overturning events in the nose region, as a result of a Rayleigh–Taylor instability.
4
Xia, Jianghai, Richard D. Miller, and Choon B. Park. "Estimation of near‐surface shear‐wave velocity by inversion of Rayleigh waves." GEOPHYSICS 64, no. 3 (May 1999): 691–700. http://dx.doi.org/10.1190/1.1444578.
Full textAbstract:
The shear‐wave (S-wave) velocity of near‐surface materials (soil, rocks, pavement) and its effect on seismic‐wave propagation are of fundamental interest in many groundwater, engineering, and environmental studies. Rayleigh‐wave phase velocity of a layered‐earth model is a function of frequency and four groups of earth properties: P-wave velocity, S-wave velocity, density, and thickness of layers. Analysis of the Jacobian matrix provides a measure of dispersion‐curve sensitivity to earth properties. S-wave velocities are the dominant influence on a dispersion curve in a high‐frequency range (>5 Hz) followed by layer thickness. An iterative solution technique to the weighted equation proved very effective in the high‐frequency range when using the Levenberg‐Marquardt and singular‐value decomposition techniques. Convergence of the weighted solution is guaranteed through selection of the damping factor using the Levenberg‐Marquardt method. Synthetic examples demonstrated calculation efficiency and stability of inverse procedures. We verify our method using borehole S-wave velocity measurements.
5
Zhang, Zhen-Dong, and Tariq Alkhalifah. "Wave-equation Rayleigh-wave dispersion inversion using fundamental and higher modes." GEOPHYSICS 84, no. 4 (July 1, 2019): EN57—EN65. http://dx.doi.org/10.1190/geo2018-0506.1.
Full textAbstract:
Recorded surface waves often provide reasonable estimates of the S-wave velocity in the near surface. However, existing algorithms are mainly based on the 1D layered-model assumption and require picking the dispersion curves either automatically or manually. We have developed a wave-equation-based inversion algorithm that inverts for S-wave velocities using fundamental and higher mode Rayleigh waves without picking an explicit dispersion curve. Our method aims to maximize the similarity of the phase velocity spectrum ([Formula: see text]) of the observed and predicted surface waves with all Rayleigh-wave modes (if they exist) included in the inversion. The [Formula: see text] spectrum is calculated using the linear Radon transform applied to a local similarity-based objective function; thus, we do not need to pick velocities in spectrum plots. As a result, the best match between the predicted and observed [Formula: see text] spectrum provides the optimal estimation of the S-wave velocity. We derive S-wave velocity updates using the adjoint-state method and solve the optimization problem using a limited-memory Broyden-Fletcher-Goldfarb-Shanno algorithm. Our method excels in cases in which the S-wave velocity has vertical reversals and lateral variations because we used all-modes dispersion, and it can suppress the local minimum problem often associated with full-waveform inversion applications. Synthetic and field examples are used to verify the effectiveness of our method.
6
Hsu, Kai, and Cengiz Esmersoy. "Parametric estimation of phase and group slownesses from sonic logging waveforms." GEOPHYSICS 57, no. 8 (August 1992): 978–85. http://dx.doi.org/10.1190/1.1443323.
Full textAbstract:
Sonic logging waveforms consist of a mixture of nondispersive waves, such as the P‐ and S‐headwaves, and dispersive waves, such as the Stoneley and pseudo‐Rayleigh waves in monopole logging and the flexural wave in dipole logging. Conventionally, slowness dispersion curves of various waves are estimated at each frequency, independent of data at other frequencies. This approach does not account for the fact that slowness dispersion functions in sonic logging are continuous and, in most cases, smooth functions of frequency. We describe a parametric slowness estimation method that uses this property by locally approximating the wavenumber of each wave as a linear function of frequency. This provides a parametric model for the phase and group slownesses of the waves propagating across the receiver array. The estimation of phase and group slownesses is then carried out by minimizing the squared difference between the predicted and observed waveforms. The minimization problem is nonlinear and is solved by an iterative algorithm. Examples using synthetic and field data are shown and the results are compared with those obtained by the conventional Prony method. Based on the comparison, we conclude that the parametric method is better than the conventional Prony method in providing robust and stable slowness estimates.
7
Mokhtar, T. A., R. B. Herrmann, and D. R. Russell. "Seismic velocity and Q model for the shallow structure of the Arabian shield from short‐period Rayleigh waves." GEOPHYSICS 53, no. 11 (November 1988): 1379–87. http://dx.doi.org/10.1190/1.1442417.
Full textAbstract:
The shear velocity and [Formula: see text] structure for the upper 1 km in different tectonic regions of the Arabian shield have been investigated using high‐frequency vertical component records of Rayleigh waves (1–20 Hz), which were recorded at source‐to‐receiver separations 55–80 km during a deep seismic refraction survey. Group and phase velocities of the fundamental and first higher modes were inverted for the shear‐wave velocity structure; Rayleigh‐wave attenuation coefficients were determined from the decay of the amplitude spectrum of the fundamental mode and used to invert for the [Formula: see text] structure. Models derived from the data were tested by calculating synthetic seismograms for the fundamental and the first higher modes from surface‐wave theory with a center of compression used to represent the source; both band‐pass filtered step and Dirac delta source time functions were tested. Modeling indicates that the shear‐wave velocity of the shield increases from 2.6 km/s to 3.4 km/s in the upper 400 m of the crust. [Formula: see text] increases from 30 in the upper 50 m to 150 at 500 m depth. The underlying material has a [Formula: see text] of 400–500 for the out‐cropping igneous rocks such as granite and may reach values higher than 700 for the metamorphic green schist rock. A band‐pass filtered Dirac delta source time function produces the synthetic that is the best fit with observations.
8
Hobiger, Manuel, Paolo Bergamo, Walter Imperatori, Francesco Panzera, Agostiny Marrios Lontsi, Vincent Perron, Clotaire Michel, Jan Burjánek, and Donat Fäh. "Site Characterization of Swiss Strong-Motion Stations: The Benefit of Advanced Processing Algorithms." Bulletin of the Seismological Society of America 111, no. 4 (June 8, 2021): 1713–39. http://dx.doi.org/10.1785/0120200316.
Full textAbstract:
ABSTRACT Since 2009, 91 new strong-motion stations were built for the renewal of the Swiss Strong Motion Network. Another nine stations will be installed until 2022. For each new station, an extensive site characterization study is performed to model the 1D seismic-velocity profile and, for some sites, the liquefaction potential. Geophysical (passive and active surface-wave methods) and geotechnical cone penetration test (CPT) with additional pore-pressure measurement (CPTu) and seismic CPT (SCPT) methods are used. Analyzing the passive and active recordings with a variety of established and advanced methods, the fundamental frequency of the site, the polarization of the wavefield, the Love- and Rayleigh-wave phase-velocity dispersion curves, and the Rayleigh-wave ellipticity function are retrieved. The liquefaction potential is assessed using CPTu. SCPT measurements are sometimes used to determine the shallow underground structure. The benefits of the combination of different appropriate methods are shown for two examples—the borehole station SBUS in Buochs and the upcoming borehole station SCME in Collombey-Muraz. At both the sites, the CPTu measurements show an elevated liquefaction potential. Combining the passive and active data, the dispersion curves for Love and Rayleigh waves and Rayleigh-wave ellipticity curves are retrieved over a wide-frequency range and inverted for the S-wave velocity profile, in which the shallow part is constrained by the active or SCPT data, the intermediate part by the dispersion curves of the passive methods, and the deepest part by the ellipticity information. For Buochs, the 1D SH-wave amplification functions modeled for the velocity profiles are compared with the empirical amplification for earthquake recordings. Finally, an overview of the site characterization results for 52 of the newly installed seismic stations is given.
9
MIURA, Hiroyuki, Atsuko MATSUO, Tatsuo KANNO, Michiko SHIGEFUJI, and Tetsuo ABIRU. "Estimation of S-Wave Velocity Structure Model by Joint Inversion of Site Amplification, Receiver Function and Phase Velocity of Rayleigh Wave." Journal of JAEE 17, no. 5 (2017): 5_78–5_95. http://dx.doi.org/10.5610/jaee.17.5_78.
Full text
10
Asten, Michael W., William J. Stephenson, and Stephen Hartzell. "Spatially averaged coherencies (krSPAC) and Rayleigh effective-mode modeling of microtremor data from asymmetric arrays." GEOPHYSICS 84, no. 3 (May 1, 2019): EN47—EN56. http://dx.doi.org/10.1190/geo2018-0524.1.
Full textAbstract:
The spatial autocorrelation and spatially averaged coherency (SPAC) methods of processing Rayleigh-wave microtremor noise observations for estimation of S-wave velocity profiles traditionally require a circular or triangular array symmetry to allow spatial (azimuthal) averaging of interstation coherencies over a constant station separation. Common processing methods allow for station separations to vary by typically [Formula: see text] in the azimuthal averaging before degradation of the SPAC spectrum is excessive. Transformation of a set of frequency-coherency spectra to wavenumber-coherency spectra (kr spectra) allows spatial averaging of spectra from multiple pairs of sensors irrespective of differences in spatial separation of the pairs. The method is called krSPAC and is implemented by iterative direct fitting of observed and model kr spectra to determine an optimal layered-earth S-wave velocity model. The observed kr spectra are updated with each iteration of the velocity model because the wavenumber is a function of model phase velocity that varies with each iteration of the modeling process. The method proves applicable when modeling either with the assumption of fundamental mode Rayleigh-wave propagation or with a summation of fundamental and higher modes. The method proves robust when compared with alternative methodologies using symmetric and asymmetric arrays on a sample of synthetic data and on field data in which station spacings vary from 70 to 800 m side lengths.
11
Hansen, Samantha E., Andrew A. Nyblade, David S. Heeszel, Douglas A. Wiens, Patrick Shore, and Masaki Kanao. "Crustal structure of the Gamburtsev Mountains, East Antarctica, from S-wave receiver functions and Rayleigh wave phase velocities." Earth and Planetary Science Letters 300, no. 3-4 (December 2010): 395–401. http://dx.doi.org/10.1016/j.epsl.2010.10.022.
Full text
12
Al-Amri, Abdullah M. "Lithospheric structure of the Arabian Shield from joint inversion of P- and S-wave receiver functions and dispersion velocities." Acta Geologica Polonica 65, no. 2 (June 1, 2015): 239–55. http://dx.doi.org/10.1515/agp-2015-0009.
Full textAbstract:
Abstract New velocity models of lithospheric thickness and velocity structure have been developed for the Arabian Shield by three tasks: 1) Computing P-Wave Receiver Functions (PRFs) and S-Wave Receiver Functions (SRFs) for all the broadband stations within the Saudi seismic networks. The number of receiver function waveforms depends on the recording time window and quality of the broadband station. 2) Computing ambient noise correlation Green’s functions for all available station pairs within the Saudi seismic networks to image the shear velocity in the crust and uppermost mantle beneath the Arabian Peninsula. Together they provided hundreds of additional, unique paths exclusively sampling the region of interest. Both phase and group velocities for all the resulting empirical Green’s functions have been measured and to be used in the joint inversion. 3) Jointly inverted the PRFs and SRFs obtained in task 1 with dispersion velocities measured on the Green’s functions obtained in task 2 and with fundamental-mode, Rayleigh-wave, group and phase velocities borrowed from the tomographic studies to precisely determine 1D crustal velocity structure and upper mantle. The analysis of the PRFs revealed values of 25-45 km for crustal thickness, with the thin crust next to the Red Sea and Gulf of Aqaba and the thicker crust under the platform, and Vp/Vs ratios in the 1.70-1.80 range, suggesting a range of compositions (felsic to mafic) for the shield’s crust. The migrated SRFs suggest lithospheric thicknesses in the 80-100 km range for portions of the shield close to the Red Sea and Gulf of Aqaba and near the Arabian Gulf. Generally, the novelty of the velocity models developed under this paper has consisted in the addition of SRF data to extend the velocity models down to lithospheric and sub-lithospheric depths.
13
Kassaras, I., F. Louis, A. Magganas, K. Makropoulos, and G. Kaviris. "Anelasticity beneath the Aegean inferred from Rayleigh wave attenuation." Bulletin of the Geological Society of Greece 40, no. 3 (June 5, 2018): 1103. http://dx.doi.org/10.12681/bgsg.16829.
Full textAbstract:
Anelasticity of the Earth crucially affects the propagation of seismic waves especially, in the long period range. However, even though the elastic properties of the Aegean deep lithosphère and upper mantle have been thoroughly investigated, their quantitative anelastic properties that influence the long period wavefield are still largely unknown. This work is towards contributing to the better knowledge of the deep structure of the Aegean by introducing experimental anelastic parameters via the study of long period Rayleigh waves attenuation. For this scope, fundamental mode attenuation coefficients (γ%) have been obtained for different two-station great-circle paths across the Aegean. The data used were provided by a broadband array installed in the area for 6 months in 1997. More than 1100 seismograms were analyzed in the 10-100 s range to obtain 17 sets of path average γR(T) functions. The attenuation coefficients are in the range 2.5*10~3 — 0.15 x 10' km' and correlate sufficiently with both experimental measurements in active tectonic regions elsewhere and synthetics generated with the use of an attenuation reference model inferred from other sources. By applying a stochastic uncoupled causal inversion method an average joint Qß'1 and shear velocity model representative of the under study area was obtained. Furthermore, path average JR(T) functions were combined in a continuous regionalization tomographic scheme to obtain local yR(T) and tomograms were constructed in the range 10-60 s. The most prominent feature in the tomograms is a high attenuation region in the central and north Aegean. This region is located south of the North Anatolian Trough and correlates well with a low shear velocity zone inferred from surface wave phase velocities. Moreover, it is associated with observed intense extensional deformation rates, mantle olivine anisotropy, recent volcanism and high heat flow.
14
Bhuiyan, Mohammad Zahidul H., and Elena Simona Lohan. "Advanced Multipath Mitigation Techniques for Satellite-Based Positioning Applications." International Journal of Navigation and Observation 2010 (December 9, 2010): 1–15. http://dx.doi.org/10.1155/2010/412393.
Full textAbstract:
Multipath remains a dominant source of ranging errors in Global Navigation Satellite Systems (GNSS), such as the Global Positioning System (GPS) or the future European satellite navigation system Galileo. Multipath is generally considered undesirable in the context of GNSS, since the reception of multipath can make significant distortion to the shape of the correlation function used for time delay estimation. However, some wireless communications techniques exploit multipath in order to provide signal diversity though in GNSS, the major challenge is to effectively mitigate the multipath, since we are interested only in the satellite-receiver transit time offset of the Line-Of-Sight (LOS) signal for the receiver's position estimate. Therefore, the multipath problem has been approached from several directions in order to mitigate the impact of multipath on navigation receivers, including the development of novel signal processing techniques. In this paper, we propose a maximum likelihood-based technique, namely, the Reduced Search Space Maximum Likelihood (RSSML) delay estimator, which is capable of mitigating the multipath effects reasonably well at the expense of increased complexity. The proposed RSSML attempts to compensate the multipath error contribution by performing a nonlinear curve fit on the input correlation function, which finds a perfect match from a set of ideal reference correlation functions with certain amplitude(s), phase(s), and delay(s) of the multipath signal. It also incorporates a threshold-based peak detection method, which eventually reduces the code-delay search space significantly. However, the downfall of RSSML is the memory requirement which it uses to store the reference correlation functions. The multipath performance of other delay-tracking methods previously studied for Binary Phase Shift Keying-(BPSK-) and Sine Binary Offset Carrier- (SinBOC-) modulated signals is also analyzed in closed loop model with the new Composite BOC (CBOC) modulation chosen for Galileo E1 signal. The simulation results show that the RSSML achieves the best multipath mitigation performance in a uniformly distributed two-to-four paths Rayleigh fading channel model for all three modulated signals.
15
Movaghari, R., and G. Javan Doloei. "3-D crustal structure of the Iran plateau using phase velocity ambient noise tomography." Geophysical Journal International 220, no. 3 (December 17, 2019): 1555–68. http://dx.doi.org/10.1093/gji/ggz537.
Full textAbstract:
SUMMARY More accurate crustal structure models will help us to better understand the tectonic convergence between Arabian and Eurasian plates in the Iran plateau. In this study, the crustal and uppermost mantle velocity structure of the Iran plateau is investigated using ambient noise tomography. Three years of continuous data are correlated to retrieve Rayleigh wave empirical Green's functions, and phase velocity dispersion curves are extracted using the spectral method. High-resolution Rayleigh wave phase velocity maps are presented at periods of 8–60 s. The tomographic maps show a clear consistency with geological structures such as sedimentary basins and seismotectonic zones, especially at short periods. A quasi-3-D shear wave velocity model is determined from the surface down to 100 km beneath the Iran plateau. A transect of the shear wave velocity model has been considered along with a profile extending across the southern Zagros, the Sanandaj-Sirjan Zone (SSZ), the Urumieh-Dokhtar Magmatic Arc (UDMA) and Central Iran and Kopeh-Dagh (KD). Obvious crustal thinning and thickening are observable along the transect of the shear wave velocity model beneath Central Iran and the SSZ, respectively. The observed shear wave velocities beneath the Iran plateau, specifically Central Iran, support the slab break-off idea in which low density asthenospheric materials drive towards the upper layers, replacing materials in the subcrustal lithosphere.
16
Lu, Zhiqu. "An acoustic near surface soil profiler using surface wave method." Journal of the Acoustical Society of America 151, no. 4 (April 2022): A58. http://dx.doi.org/10.1121/10.0010649.
Full textAbstract:
An acoustic soil profiler, using a so-called the high-frequency multi-channel analysis of surface waves (HF-MASW) method, has been developed, which uses surface (Rayleigh) waves to measure soil profile in terms of the shear (S) wave velocity as a function of depth, up to a 2.5 m deep below the surface. Several practical techniques have been developed to enhance the HF-MASW method, including (1) a variable sensor spacing configuration, (2) the self-adaptive method, and (3) the phase-only signal processing. Fundamentally, the S-wave velocity is related to soil mechanical and hydrological properties through the principle of effective stress. Therefore, the measured two-dimensional S-wave velocity images reflect the temporal and spatial variations of soils due to weather effects, geological anomalies, and anthropologic activities. Several HF-MASW applications will be reported, including (1) near surface soil profiling, (2) a long-term-survey for studying weather and seasonal effects, (3) short-term monitoring rain fall events, (4) detecting fraigpan layers, and (5) a farmland compaction study. This acoustic soil profiler can be used for agricultural, environmental, civil engineering, and military applications.
17
Dolan, Brenda, and Steven A. Rutledge. "A Theory-Based Hydrometeor Identification Algorithm for X-Band Polarimetric Radars." Journal of Atmospheric and Oceanic Technology 26, no. 10 (October 1, 2009): 2071–88. http://dx.doi.org/10.1175/2009jtecha1208.1.
Full textAbstract:
Abstract Although much work has been done at S band to automatically identify hydrometeors by using polarimetric radar, several challenges are presented when adapting such algorithms to X band. At X band, attenuation and non-Rayleigh scattering can pose significant problems. This study seeks to develop a hydrometeor identification (HID) algorithm for X band based on theoretical simulations using the T-matrix scattering model of seven different hydrometeor types: rain, drizzle, aggregates, pristine ice crystals, low-density graupel, high-density graupel, and vertical ice. Hail and mixed-phase hydrometeors are excluded for the purposes of this study. Non-Rayleigh scattering effects are explored by comparison with S-band simulations. Variable ranges based on the theoretical simulations are used to create one-dimensional fuzzy-logic membership beta functions that form the basis of the new X-band HID. The theory-based X-band HID is applied to a case from the Collaborative Adaptive Sensing of the Atmosphere (CASA) Integrated Project 1 (IP1) network of X-band radars, and comparisons are made with similar S-band hydrometeor identification algorithms applied to data from the S-band polarimetric Next Generation Weather Radar (NEXRAD) prototype radar, KOUN. The X-band HID shows promise for illustrating bulk hydrometeor types and qualitatively agrees with analysis from KOUN. A simple reflectivity- and temperature-only HID is also applied to both KOUN and CASA IP1 data to reveal benefits of the polarimetric-based HID algorithms, especially in the classification of ice hydrometeors and oriented ice crystals.
18
Dunham, C. K., J. P. O’Donnell, G. W. Stuart, A. M. Brisbourne, S. Rost, T. A. Jordan, A. A. Nyblade, D. A. Wiens, and R. C. Aster. "A joint inversion of receiver function and Rayleigh wave phase velocity dispersion data to estimate crustal structure in West Antarctica." Geophysical Journal International 223, no. 3 (August 22, 2020): 1644–57. http://dx.doi.org/10.1093/gji/ggaa398.
Full textAbstract:
SUMMARY We determine crustal shear wave velocity structure and crustal thickness at recently deployed seismic stations across West Antarctica, using a joint inversion of receiver functions and fundamental mode Rayleigh wave phase velocity dispersion. The stations are from both the UK Antarctic Network (UKANET) and Polar Earth Observing Network/Antarctic Network (POLENET/ANET). The former include, for the first time, four stations along the spine of the Antarctic Peninsula, three in the Ellsworth Land and five stations in the vicinity of the Pine Island Rift. Within the West Antarctic Rift System (WARS) we model a crustal thickness range of 18–28 km, and show that the thinnest crust (∼18 km) is in the vicinity of the Byrd Subglacial Basin and Bentley Subglacial Trench. In these regions we also find the highest ratio of fast (Vs = 4.0–4.3 km s–1, likely mafic) lower crust to felsic/intermediate upper crust. The thickest mafic lower crust we model is in Ellsworth Land, a critical area for constraining the eastern limits of the WARS. Although we find thinner crust in this region (∼30 km) than in the neighbouring Antarctic Peninsula and Haag-Ellsworth Whitmore block (HEW), the Ellsworth Land crust has not undergone as much extension as the central WARS. This suggests that the WARS does not link with the Weddell Sea Rift System through Ellsworth Land, and instead has progressed during its formation towards the Bellingshausen and Amundsen Sea Embayments. We also find that the thin WARS crust extends towards the Pine Island Rift, suggesting that the boundary between the WARS and the Thurston Island block lies in this region, ∼200 km north of its previously accepted position. The thickest crust (38–40 km) we model in this study is in the Ellsworth Mountain section of the HEW block. We find thinner crust (30–33 km) in the Whitmore Mountains and Haag Nunatak sectors of the HEW, consistent with the composite nature of the block. In the Antarctic Peninsula we find a crustal thickness range of 30–38 km and a likely dominantly felsic/intermediate crustal composition. By forward modelling high frequency receiver functions we also assess if any thick, low velocity subglacial sediment accumulations are present, and find a 0.1–0.8-km-thick layer at 10 stations within the WARS, Thurston Island and Ellsworth Land. We suggest that these units of subglacial sediment could provide a source region for the soft basal till layers found beneath numerous outlet glaciers, and may act to accelerate ice flow.
19
Borcherdt, Roger D., and Leif Wennerberg. "General P, type-I S, and type-II S waves in anelastic solids; inhomogeneous wave fields in low-loss solids." Bulletin of the Seismological Society of America 75, no. 6 (December 1, 1985): 1729–63. http://dx.doi.org/10.1785/bssa0750061729.
Full textAbstract:
Abstract The physical characteristics for general plane-wave radiation fields in an arbitrary linear viscoelastic solid are derived. Expressions for the characteristics of inhomogeneous wave fields, derived in terms of those for homogeneous fields, are utilized to specify the characteristics and a set of reference curves for general P and S wave fields in arbitrary viscoelastic solids as a function of wave inhomogeneity and intrinsic material absorption. The expressions show that an increase in inhomogeneity of the wave fields causes the velocity to decrease, the fractional-energy loss (Q−1) to increase, the deviation of maximum energy flow with respect to phase propagation to increase, and the elliptical particle motions for P and type-I S waves to approach circularity. Q−1 for inhomogeneous type-I S waves is shown to be greater than that for type-II S waves, with the deviation first increasing then decreasing with inhomogeneity. The mean energy densities (kinetic, potential, and total), the mean rate of energy dissipation, the mean energy flux, and Q−1 for inhomogeneous waves are shown to be greater than corresponding characteristics for homogeneous waves, with the deviations increasing as the inhomogeneity is increased for waves of fixed maximum displacement amplitude. For inhomogeneous wave fields in low-loss solids, only the tilt of the particle motion ellipse for P and type-I S waves is independent to first order of the degree of inhomogeneity. Quantitative estimates for the characteristics of inhomogeneous plane body waves in layered low-loss solids are derived and guidelines established for estimating the effect of inhomogeneity on seismic body waves and a Rayleigh-type surface wave in low-loss media.
20
Yamanaka, Hiroaki, and Kentaro Motoki. "Joint inversion of Rayleigh wave phase velocity and receiver function for estimation of S-wave velocity of deep sedimentary layers in the Kashiwazaki city, Japan." BUTSURI-TANSA(Geophysical Exploration) 62, no. 2 (2009): 237–47. http://dx.doi.org/10.3124/segj.62.237.
Full text
21
Acevedo, Jorge, Gabriela Fernández-Viejo, Sergio Llana-Fúnez, Carlos López-Fernández, and Javier Olona. "Ambient noise tomography of the southern sector of the Cantabrian Mountains, NW Spain." Geophysical Journal International 219, no. 1 (July 8, 2019): 479–95. http://dx.doi.org/10.1093/gji/ggz308.
Full textAbstract:
SUMMARY This study presents the first detailed analysis of ambient noise tomography in an area of the continental upper crust in the Cantabrian Mountains (NW Spain), where a confluence of crustal scale faults occurs at depth. Ambient noise data from two different seismic networks have been analysed. In one side, a 10-short-period station network was set recording continuously for 19 months. A second set of data from 13 broad-band stations was used to extend at depth the models. The phase cross-correlation processing technique was used to compute in total more than 34 000 cross-correlations from 123 station pairs. The empirical Green's functions were obtained by applying the time–frequency, phase-weighted stacking methodology and provided the emergence of Rayleigh waves. After measuring group velocities, Rayleigh-wave group velocity tomographic maps were computed at different periods and then they were inverted in order to calculate S-wave velocities as a function of depth, reaching the first 12 km of the crust. The results show that shallow velocity patterns are dominated by geological features that can be observed at the surface, particularly bedding and/or lithology and fracturing associated with faults. In contrast, velocity patterns below 4 km depth seem to be segmented by large structures, which show a velocity reduction along fault zones. The best example is the visualization in the tomography of the frontal thrust of the Cantabrian Mountains at depth, which places higher velocity Palaeozoic rocks over Cenozoic sediments of the foreland Duero basin. One of the major findings in the tomographic images is the reduction of seismic velocities above the area in the crust where one seismicity cluster is nucleated within the otherwise quiet seismic area of the range. The noise tomography reveals itself as a valuable technique to identify shear zones associated with crustal scale fractures and hence, lower strain areas favourable to seismicity.
22
Chrapkiewicz, Kajetan, Monika Wilde-Piórko, Marcin Polkowski, and Marek Grad. "Reliable workflow for inversion of seismic receiver function and surface wave dispersion data: a “13 BB Star” case study." Journal of Seismology 24, no. 1 (December 16, 2019): 101–20. http://dx.doi.org/10.1007/s10950-019-09888-1.
Full textAbstract:
AbstractNon-linear inverse problems arising in seismology are usually addressed either by linearization or by Monte Carlo methods. Neither approach is flawless. The former needs an accurate starting model; the latter is computationally intensive. Both require careful tuning of inversion parameters. An additional challenge is posed by joint inversion of data of different sensitivities and noise levels such as receiver functions and surface wave dispersion curves. We propose a generic workflow that combines advantages of both methods by endowing the linearized approach with an ensemble of homogeneous starting models. It successfully addresses several fundamental issues inherent in a wide range of inverse problems, such as trapping by local minima, exploitation of a priori knowledge, choice of a model depth, proper weighting of data sets characterized by different uncertainties, and credibility of final models. Some of them are tackled with the aid of novel 1D checkerboard tests—an intuitive and feasible addition to the resolution matrix. We applied our workflow to study the south-western margin of the East European Craton. Rayleigh wave phase velocity dispersion and P-wave receiver function data were gathered in the passive seismic experiment “13 BB Star” (2013–2016) in the area of the crust recognized by previous borehole and refraction surveys. Final models of S-wave velocity down to 300 km depth beneath the array are characterized by proximity in the parameter space and very good data fit. The maximum value in the mantle is higher by 0.1–0.2 km/s than reported for other cratons.
23
Koroleva, Tat’iana, and Evgeniia Lyskova. "Causes of artifacts in ambient noise surface wave tomography in mantle investigations and ways for their elimination." Russian Journal of Seismology 2, no. 2 (June 23, 2020): 58–65. http://dx.doi.org/10.35540/2686-7907.2020.2.05.
Full textAbstract:
Ambient noise surface wave tomography is a widely used method for determining the velocity structure of the upper layers of the Earth. It is based on the fact that the cross-correlation function (CCF) of noise at two stations, averaged over a long time interval, determines the Green's function of the surface wave. This allows us to estimate the group and phase velocities of surface waves on the paths between stations, which are used in surface-wave tomography. This makes it possible to ultimately estimate the spatial distribution of the S-wave velocities. The method is well-grounded on the assumption that the “noise” is a result of the superposition of surface waves propagating from sources uniformly distributed over the surface. Therefore, the initial data, which are long-period seismic records, are subjected to preliminary processing, an important stage of which is normalization, which allows reducing the effect of earthquakes and averaging the resulting CCFs over a long time interval. At the same time, we have shown that earthquakes mainly contribute to noise at periods above 30-40 s, whose sources are distributed unevenly. Therefore, in cases of clustering of foci in a certain limited area, for example, because of aftershocks after a strong earthquake, the CCF maxima, which determines the dispersion curve of the surface wave, are shifted to shorter times, and the group velocities are correspondingly overestimated. In determining the dispersion of Love waves from the CCF transversal (T-T) noise component, the presence of clusters leads to an additional underestimation of the group velocity due to the superposition on the T component (perpendicular to the inter-station path) of the radial component of the Rayleigh wave having a velocity less than the Love wave velocity. Therefore, the anisotropy coefficient, determined from the noise, is underestimated as compared to that obtained from the records of earthquakes along nearby paths. Obviously, to obtain more correct dispersion curves of both Rayleigh and Love waves, it is necessary, for summing the CCFs, to use time intervals in which earthquake clusters would be absent as far as possible.
24
Liu, Xin, and Dapeng Zhao. "Seismic evidence for a plume-modified oceanic lithosphere–asthenosphere system beneath Cape Verde." Geophysical Journal International 225, no. 2 (January 11, 2021): 872–86. http://dx.doi.org/10.1093/gji/ggab012.
Full textAbstract:
SUMMARY We determine a new 3-D shear wave velocity (Vs) model down to 400 km depth beneath the Cape Verde hotspot that is far from plate boundaries. This Vs model is obtained by using a new method of jointly inverting P- and S-wave receiver functions, Rayleigh-wave phase-velocity data and S-wave arrival times of teleseismic events. Two Vs discontinuities at ∼15 and ∼60 km depths are revealed beneath volcanic islands, which are interpreted as the Moho discontinuity and the Gutenberg (G) discontinuity. Between the north and south islands, obvious high-Vs anomalies exist in the uppermost mantle down to a depth of ∼100–150 km beneath the Atlantic Ocean, whereas obvious low-Vs anomalies exist in the uppermost mantle beneath the volcanic islands including the active Fogo volcano. These low-Vs anomalies merge into a significant column-like low-Vs zone at depths of ∼150–400 km beneath the Cape Verde swell. We propose that these features in the upper mantle reflect a plume-modified oceanic lithosphere–asthenosphere system beneath the Cape Verde hotspot.
25
Cheng, Ningya, and Chuen Hon Cheng. "Estimations of formation velocity, permeability, and shear‐wave anisotropy using acoustic logs." GEOPHYSICS 61, no. 2 (March 1996): 437–43. http://dx.doi.org/10.1190/1.1443971.
Full textAbstract:
Field data sets collected by an array monopole acoustic logging tool and a shear wave logging tool are processed and interpreted. The P‐ and S‐wave velocities of the formation are determined by threshold detection with cross‐correlation correction from the full waveform and the shear‐wave log, respectively. The array monopole acoustic logging data are also processed using the extended Prony’s method to estimate the borehole Stoneley wave phase velocity and attenuation as a function of frequency. The well formation between depths of 2950 and 3150 ft (899 and 960 m) can be described as an isotropic elastic medium. The inverted [Formula: see text] from the Stoneley wave phase velocity is in excellent agreement with the shear‐wave log results in this section. The well formation between the depths of 3715 and 3780 ft (1132 and 1152 m) can be described as a porous medium with shear‐wave velocity anisotropy about 10% to 20% and with the symmetry axis perpendicular to the borehole axis. The disagreement between the shear‐wave velocity from the Stoneley wave inversion and the direct shear‐wave log velocity in this section is beyond the errors in the measurements. Estimated permeabilities from low‐frequency Stoneley wave velocity and attenuation data are in good agreement with the core measurements. Also it is proven that the formation permeability is not the cause of the discrepancy. From the estimated “shear/pseudo‐Rayleigh” phase velocities in the array monopole log and the 3-D finite‐difference synthetics in the anisotropic formation, the discrepancy can best be explained as shear‐wave anisotropy.
26
Kleiss, Jessica M., and W. Kendall Melville. "Observations of Wave Breaking Kinematics in Fetch-Limited Seas." Journal of Physical Oceanography 40, no. 12 (December 1, 2010): 2575–604. http://dx.doi.org/10.1175/2010jpo4383.1.
Full textAbstract:
Abstract Breaking waves play an important role in air–sea interaction, enhancing momentum flux from the atmosphere to the ocean, dissipating wave energy that is then available for turbulent mixing, injecting aerosols and sea spray into the atmosphere, and affecting air–sea gas transfer due to air entrainment. In this paper observations are presented of the occurrence of breaking waves under conditions of strong winds (10–25 m s−1) and fetch-limited seas (0–500 km) in the Gulf of Tehuantepec Experiment (GOTEX) in 2004. An airborne nadir-looking video camera, along with a global positioning system (GPS) and inertial motion unit (IMU), provided digital videos of the breaking sea surface and position in an earth frame. In particular, the authors present observations of Λ(c), which is the distribution of breaking wave crest lengths per unit sea surface area, per unit increment in velocity c or scalar speed c, first introduced by O. M. Phillips. In another paper, the authors discuss the effect of processing methodology on the resulting shape of the Λ(c) distribution. In this paper, the elemental method of measuring breaking crests is used to investigate the Λ(c) distributions under a variety of wind and wave conditions. The integral and the first two moments of the Λ(c) distributions are highly correlated with the active breaking rate and the active whitecap coverage. The computation of whitecap coverage yields a larger observational dataset from which the variability of whitecap coverage with wind speed, friction velocity, wave age, and wave slope is presented and compared to previous observations. The dependence of the active breaking rate on the spectral peak steepness is in agreement with previous studies. Dimensional analysis of Λ(c) indicates that scaling with friction velocity and gravity, as in the classical fetch relations, collapses the breaking distributions more effectively than scaling with dominant wave parameters. Significant wave breaking is observed at speeds near the spectral peak in young seas only, consistent with previous studies. The fourth and fifth moments of Λ(c) are related to the flux of momentum transferred by breaking waves to the underlying water and the rate of wave energy dissipation, respectively. The maximum in the fourth moment occurs at breaking speeds of 5–5.5 m s−1, and the maximum in the fifth moment occurs at 5.8–6.8 m s−1, apparently independent of wave age. However, when nondimensionalized by the phase speed at the peak of the local wave spectrum cp, the maxima in the nondimensionalized fourth and fifth moments show a decreasing trend with wave age, obtaining the maxima at dimensionless speeds c/cp near unity at smaller wave ages and moving to lower dimensionless speeds c/cp ≪ 1 at larger wave ages. The angular dependence of Λ(c) is predominantly unimodal and better aligned with the wind direction than the dominant wave direction. However, the directional distribution of Λ(c) is broadest for small c and often exhibits a bimodal structure for slow breaking speeds under developing seas. An asymmetry in the directional distribution is also observed for moderately developed seas. Observations are compared to the Phillips model for Λ(c) in the equilibrium range of the wave spectrum. Although the ensemble of Λ(c) distributions appears consistent with a c−6 function, the distributions are not described by a constant power-law exponent. However, the Λ(c) observations are described well by the Rayleigh distribution for slow and intermediate speeds, yet fall above the Rayleigh distribution for the fastest breaking speeds. From the Rayleigh description, it is found that the dimensionless width of the Λ(c) distribution increases weakly with dimensionless fetch, s/u*e = 1.69χ0.06, where s is the Rayleigh parameter, u*e is the effective friction velocity, and the dimensionless fetch is a function of the fetch X and gravitational acceleration g. The nondimensionalized total length of breaking per unit sea surface area is found to decrease with dimensionless fetch for intermediate to fully developed seas, , where A is the total length of breaking crests per unit sea surface area.
27
Li, Xuantao, Jinli Huang, and Zhikun Liu. "Ambient-Noise Tomography of the Baiyun Gold Deposit in Liaoning, China." Seismological Research Letters 91, no. 5 (July 15, 2020): 2791–802. http://dx.doi.org/10.1785/0220190393.
Full textAbstract:
Abstract Ambient-noise tomography (ANT) has become an effective method for determining the fine velocity structure of the shallow crust. However, studies on metal mines using this method are rarely reported. To investigate the tectonic background and prospecting of the deep mine in the Baiyun gold deposit (BYGD) of eastern Liaoning Province, China, we use ANT to determine a 3D S-wave velocity structure model of the BYGD. A total of 21 broadband seismic stations were installed in an area of 15×14 km, centered at the BYGD. Continuous observations for approximately three months were made. After single-station preprocessing, cross correlation of ambient noise, and phase-weighted stacking, the empirical Green’s function for the Rayleigh waves between stations was recovered. Next, group-velocity dispersion with 0.8–3 s periods was measured. A direct inversion method of surface-wave dispersion based on raytracing was then adopted to determine a 3D S-wave velocity structure of the BYGD from the ground surface to a depth of 1.8 km. The results show that the distribution of S-wave velocities in the study area well reflected the geological characteristics of the surface. The velocities were significantly low within the “ore field” and the regional ore-controlling Jianshanzi fault. Combining this with the fact that a large number of magmatic veins were visible inside both structures, it was deduced that both structures had experienced large-scale magmatic intrusion activities, thus confirming that BYGD was a magmatic hydrothermal deposit. The significantly low S-wave velocities beneath the gold deposit extended to a depth of 1.8 km. This might imply the occurrence of blind ore bodies at that depth. The fine velocity structure of the BYGD reconstructed by this study provided a direction for subsequent prospecting of deep regions and demonstrated that ANT has good potential in metal mine exploration.
28
WU, XUESONG, P. A. STEWART, and S. J. COWLEY. "On the catalytic role of the phase-locked interaction of Tollmien–Schlichting waves in boundary-layer transition." Journal of Fluid Mechanics 590 (October 15, 2007): 265–94. http://dx.doi.org/10.1017/s002211200700804x.
Full textAbstract:
This paper is concerned with the nonlinear interaction between a planar and a pair of oblique Tollmien–Schlichting (T-S) waves which are phase-locked in that they travel with (nearly) the same phase speed. The evolution of such a disturbance is described using a high-Reynolds-number asymptotic approach in the so-called ‘upper--branch’ scaling regime. It follows that there exists a well-defined common critical layer (i.e. a thin region surrounding the level at which the basic flow velocity equals the phase speed of the waves to leading order) and the dominant interactions take place there. The disturbance is shown to evolve through several distinctive stages. In the first of these, the critical layer is in equilibrium and viscosity dominated. If a small mismatching exists in the phase speeds, the interaction between the planar and oblique waves leads directly to super-exponential growth/decay of the oblique modes. However, if the modes are perfectly phase-locked, the interaction in the first instance affects only the phase of the amplitude function of the oblique modes (so causing rapid wavelength shortening), while the modulus of the amplitude still evolves exponentially until the wavelength shortening produces a back reaction on the modulus (which then induces a super-exponential growth). Whether or not there is a small mismatch or a perfect match in the phase speeds, once the growth rate of the oblique modes becomes sufficiently large, the disturbance enters a second stage, in which the critical layer becomes both non-equilibrium and viscous in nature. The oblique modes continue to experience super-exponential growth, albeit of a different form from that in the previous stages, until the self-interaction between them, as well as their back effect on the planar mode, becomes important. At that point, the disturbance enters a third, fully interactive stage, during which the development of the disturbance is governed by the amplitude equations with the same nonlinear terms as previously derived for the phase-locked interaction of Rayleigh instability waves. The solution develops a singularity, leading to the final stage where the flow is governed by fully nonlinear three-dimensional inviscid triple-deck equations. The present work indicates that seeding a planar T-S wave can enhance the amplification of all oblique modes which share approximately its phase speed.
29
Pang, Mengqiang, Jing Ba, Li-Yun Fu, José M. Carcione, Uti I. Markus, and Lin Zhang. "Estimation of microfracture porosity in deep carbonate reservoirs based on 3D rock-physics templates." Interpretation 8, no. 4 (July 23, 2020): SP43—SP52. http://dx.doi.org/10.1190/int-2019-0258.1.
Full textAbstract:
Carbonate reservoirs in the S area of the Tarim Basin (China) are ultradeep hydrocarbon resources, with low porosity, complex fracture systems, and dissolved pores. Microfracturing is a key factor of reservoir connectivity and storage space. We have performed measurements on limestone samples, under different confining pressures, and we used the self-consistent approximation model and the Biot-Rayleigh theory of double porosity to study the microfractures. We have computed the fluid properties (mainly oil) as a function of temperature and pressure. Using the dependence of seismic [Formula: see text] on the microfractures, a multiscale 3D rock-physics template (RPT) is built, based on the attenuation, P-wave impedance, and phase velocity ratio. We estimate the ultrasonic and seismic attenuation with the spectral-ratio method and the improved frequency-shift method, respectively. Then, calibration of the RPTs is performed at ultrasonic and seismic frequencies. We use the RPTs to estimate the total and microfracture porosities. The results indicate that the total porosity is low and the microfracture porosity is relatively high, which is consistent with the well log data and actual oil production reports. This work presents a method for identification of deep carbonate reservoirs by using the microfracture porosity estimated from the 3D RPT, which could be exploited in oil and gas exploration.
30
Dhar, J., P. Meunier, F. Nadal, and Y. Méheust. "Convective dissolution of carbon dioxide in two- and three-dimensional porous media: The impact of hydrodynamic dispersion." Physics of Fluids 34, no. 6 (June 2022): 064114. http://dx.doi.org/10.1063/5.0086370.
Full textAbstract:
Convective dissolution is the process by which CO2 injected in geological formations dissolves into the aqueous phase and thus remains stored perennially by gravity. It can be modeled by buoyancy-coupled Darcy flow and solute transport. The transport equation should include a diffusive term accounting for hydrodynamic dispersion, wherein the effective diffusion coefficient is proportional to the local interstitial velocity. We investigate the impact of the hydrodynamic dispersion tensor on convective dissolution in two-dimensional (2D) and three-dimensional (3D) homogeneous porous media. Using a novel numerical model, we systematically analyze, among other observables, the time evolution of the fingers' structure, dissolution flux in the quasi-constant flux regime, and mean concentration of the dissolved CO2; we also determine the onset time of convection, [Formula: see text]. For a given Rayleigh number Ra, the efficiency of convective dissolution over long times is controlled by [Formula: see text]. For porous media with a dispersion anisotropy commonly found in the subsurface, [Formula: see text] increases as a function of the longitudinal dispersion's strength ( S), in agreement with previous experimental findings and in contrast to previous numerical findings, a discrepancy that we explain. More generally, for a given strength of transverse dispersion, longitudinal dispersion always slows down convective dissolution, while for a given strength of longitudinal dispersion, transverse dispersion always accelerates it. Furthermore, a systematic comparison between 2D and 3D results shows that they are consistent on all accounts, except for a slight difference in [Formula: see text] and a significant impact of Ra on the dependence of the finger number density on S in 3D.
31
Zhang, Shane, Lili Feng, and Michael H. Ritzwoller. "Three-station interferometry and tomography: coda versus direct waves." Geophysical Journal International 221, no. 1 (January 28, 2020): 521–41. http://dx.doi.org/10.1093/gji/ggaa046.
Full textAbstract:
SUMMARY Traditional two-station ambient noise interferometry estimates the Green’s function between a pair of synchronously deployed seismic stations. Three-station interferometry considers records observed three stations at a time, where two of the stations are considered receiver–stations and the third is a source–station. Cross-correlations between records at the source–station with each of the receiver–stations are correlated or convolved again to estimate the Green’s function between the receiver–stations, which may be deployed asynchronously. We use data from the EarthScope USArray in the western United States to compare Rayleigh wave dispersion obtained from two-station and three-station interferometry. Three three-station interferometric methods are distinguished by the data segment utilized (coda-wave or direct-wave) and whether the source–stations are constrained to lie in stationary phase zones approximately inline with the receiver–stations. The primary finding is that the three-station direct wave methods perform considerably better than the three-station coda-wave method and two-station ambient noise interferometry for obtaining surface wave dispersion measurements in terms of signal-to-noise ratio, bandwidth, and the number of measurements obtained, but possess small biases relative to two-station interferometry. We present a ray-theoretic correction method that largely removes the bias below 40 s period and reduces it at longer periods. Three-station direct-wave interferometry provides substantial value for imaging the crust and uppermost mantle, and its ability to bridge asynchronously deployed stations may impact the design of seismic networks in the future.
32
Vinnik, L. P., G. D. Georgieva, S. I. Oreshin, L. I. Makeyeva, D. N. Dragomirov, V. D. Buchakchiev, and L. D. Dimitrova. "Deep Structure and Dynamics of the Central Balkan Peninsula from Seismic Data." Izvestiya, Physics of the Solid Earth 57, no. 6 (November 2021): 849–63. http://dx.doi.org/10.1134/s1069351321060124.
Full textAbstract:
Abstract—Analysis of P- and S-receiver functions for 19 seismic stations on the Balkan Peninsula has been performed. Half of the stations are in Bulgaria. The crustal thickness varies from 28–30 to 50 km. The ratio of longitudinal and shear wave velocities in the upper crust reaches 2.0 in some places. In the southwest of the study area, the 410-km seismic boundary is uplifted by 10 km relative to nominal depth. The elevation may be caused by hydration and/or cooling of the mantle transition zone under the influence of the Hellenic subduction zone. A low S-wave velocity layer related to the 410-km boundary may be located atop this boundary. In the northwestern part of the study area this layer is present in spite of the absence of the 410-km boundary. A similar paradox has been previously noted in central Anatolia. Indications of a low-velocity layer are present at a depth exceeding 410 km. The simultaneous inversion of the receiver functions of the two types (P and S) and the Rayleigh wave phase velocities reveals a large (7–9%) decrease in the S-wave velocity in the upper mantle of southern Bulgaria and northern Greece. The thickness of the low-velocity layer (asthenosphere) is about 50 km. The lithosphere-asthenosphere boundary (LAB) is at depths of 40 to 60 km. In terms of tectonics, this zone is characterized as the South Balkan extension system. To the north of 43° N, the S-wave velocity in the upper mantle is usually at least 4.4 km/s and the LAB is not detected or is detected at a depth of over 80 km. The SKS analysis of azimuthal anisotropy reveals lateral zoning in the upper mantle that is correlated to velocity zoning. Probably, the mechanically weak low-velocity mantle of the South Balkan system is easily deformed, and the azimuth of the fast direction of anisotropy (20°) indicates the direction of extension. At the northern stations, the fast direction (about –30°) may be a reflection of an older process.
33
Jung, Youngsun, Guifu Zhang, and Ming Xue. "Assimilation of Simulated Polarimetric Radar Data for a Convective Storm Using the Ensemble Kalman Filter. Part I: Observation Operators for Reflectivity and Polarimetric Variables." Monthly Weather Review 136, no. 6 (June 1, 2008): 2228–45. http://dx.doi.org/10.1175/2007mwr2083.1.
Full textAbstract:
Abstract A radar simulator for polarimetric radar variables, including reflectivities at horizontal and vertical polarizations, the differential reflectivity, and the specific differential phase, has been developed. This simulator serves as a test bed for developing and testing forward observation operators of polarimetric radar variables that are needed when directly assimilating these variables into storm-scale numerical weather prediction (NWP) models, using either variational or ensemble-based assimilation methods. The simulator takes as input the results of high-resolution NWP model simulations with ice microphysics and produces simulated polarimetric radar data that may also contain simulated errors. It is developed based on calculations of electromagnetic wave propagation and scattering at the S band of wavelength 10.7 cm in a hydrometeor-containing atmosphere. The T-matrix method is used for the scattering calculation of raindrops and the Rayleigh scattering approximation is applied to snow and hail particles. The polarimetric variables are expressed as functions of the hydrometeor mixing ratios as well as their corresponding drop size distribution parameters and densities. The presence of wet snow and wet hail in the melting layer is accounted for by using a new, relatively simple melting model that defines the water fraction in the melting snow or hail. The effect of varying density due to the melting snow or hail is also included. Vertical cross sections and profiles of the polarimetric variables for a simulated mature multicellular squall-line system and a supercell storm show that polarimetric signatures of the bright band in the stratiform region and those associated with deep convection are well captured by the simulator.
34
Frisvad, Jeppe Revall. "Importance sampling the Rayleigh phase function." Journal of the Optical Society of America A 28, no. 12 (November 10, 2011): 2436. http://dx.doi.org/10.1364/josaa.28.002436.
Full text
35
Liu, Quanhua, and Fuzhong Weng. "Combined Henyey-Greenstein and Rayleigh phase function." Applied Optics 45, no. 28 (October 1, 2006): 7475. http://dx.doi.org/10.1364/ao.45.007475.
Full text
36
Rindraharisaona, E. J., F. Tilmann, X. Yuan, J. Dreiling, J. Giese, K. Priestley, and G. Rümpker. "Velocity structure and radial anisotropy of the lithosphere in southern Madagascar from surface wave dispersion." Geophysical Journal International 224, no. 3 (November 17, 2020): 1930–44. http://dx.doi.org/10.1093/gji/ggaa550.
Full textAbstract:
SUMMARY We investigate the upper mantle seismic structure beneath southern Madagascar and infer the imprint of geodynamic events since Madagascar’s break-up from Africa and India and earlier rifting episodes. Rayleigh and Love wave phase velocities along a profile across southern Madagascar were determined by application of the two-station method to teleseismic earthquake data. For shorter periods (<20 s), these data were supplemented by previously published dispersion curves determined from ambient noise correlation. First, tomographic models of the phase velocities were determined. In a second step, 1-D models of SV and SH wave velocities were inverted based on the dispersion curves extracted from the tomographic models. As the lithospheric mantle is represented by high velocities we identify the lithosphere–asthenosphere boundary by the strongest negative velocity gradient. Finally, the radial anisotropy (RA) is derived from the difference between the SV and SH velocity models. An additional constraint on the lithospheric thickness is provided by the presence of a negative conversion seen in S receiver functions, which results in comparable estimates under most of Madagascar. We infer a lithospheric thickness of 110−150 km beneath southern Madagascar, significantly thinner than beneath the mobile belts in East Africa (150−200 km), where the crust is of comparable age and which were located close to Madagascar in Gondwanaland. The lithospheric thickness is correlated with the geological domains. The thinnest lithosphere (∼110 km) is found beneath the Morondava basin. The pre-breakup Karoo failed rifting, the rifting and breakup of Gondwanaland have likely thinned the lithosphere there. The thickness of the lithosphere in the Proterozoic terranes (Androyen and Anosyen domains) ranges from 125 to 140 km, which is still ∼30 km thinner than in the Mozambique belt in Tanzania. The lithosphere is the thickest beneath Ikalamavony domain (Proterozoic) and the west part of the Antananarivo domain (Archean) with a thickness of ∼150 km. Below the eastern part of Archean domain the lithosphere thickness reduces to ∼130 km. The lithosphere below the entire profile is characterized by positive RA. The strongest RA is observed in the uppermost mantle beneath the Morondava basin (maximum value of ∼9 per cent), which is understandable from the strong stretching that the basin was exposed to during the Karoo and subsequent rifting episode. Anisotropy is still significantly positive below the Proterozoic (maximum value of ∼5 per cent) and Archean (maximum value of ∼6 per cent) domains, which may result from lithospheric extension during the Mesozoic and/or thereafter. In the asthenosphere, a positive RA is observed beneath the eastern part Morondava sedimentary basin and the Proterozoic domain, indicating a horizontal asthenospheric flow pattern. Negative RA is found beneath the Archean in the east, suggesting a small-scale asthenospheric upwelling, consistent with previous studies. Alternatively, the relatively high shear wave velocity in the asthenosphere in this region indicate that the negative RA could be associated to the Réunion mantle plume, at least beneath the volcanic formation, along the eastern coast.
37
Pramatadie, Andi Muhamad, Hiroaki Yamanaka, and Afnimar Afnimar. "Shallow S-wave Velocity Profile Estimation using Surface Velocity and Microtremor HVSR with a Linear Velocity Increase Approach." Journal of Mathematical and Fundamental Sciences 54, no. 3 (July 26, 2023): 330–58. http://dx.doi.org/10.5614/j.math.fund.sci.2023.54.3.4.
Full textAbstract:
We propose a simple method for 1D S-wave velocity (Vs) profile estimation using a measured surface S-wave velocity (V1) and peak frequency of the observed microtremor horizontal-to-vertical spectral ratio (HVSR). In this method, the S-wave velocity profile is presented as linear velocity increase with depth in sediments over a bedrock layer that has a given constant S-wave velocity. Thus, the profile can be parameterized with the measured surface S-wave velocity and the velocity gradient. The gradient can be estimated based on the agreement of the peak frequencies of the observed microtremor HVSR and the theoretical ellipticity of the fundamental mode of the Rayleigh wave. We examined the applicability of the proposed method using numerical tests as well as application to actual data at five sites in the Bandung Basin, Indonesia, where observed Rayleigh wave phase velocities from microtremor array surveys were available. The applicability was confirmed in numerical tests using sample models of soil profiles in the basin. Actual application indicated the appropriateness of the estimated S-wave velocity profiles due to the similarity of their theoretical Rayleigh wave phase velocities with the observed Rayleigh wave phase velocities. Since the proposed method needs prior confirmation of the linear increase of the S-wave velocity, it is suitable for use in spatial interpolation of shallow S-wave velocity profiles with simple data acquisition.
38
Zhang, Jinyun, Zhifu Zhang, Yinjie Zhang, Xuhan Lu, Xianhong Jiang, Peng Li, and Kang Li. "Rayleigh Wave Dispersion and Inversion for Shallow Surface with a High-velocity Rigid Pavement." Journal of Physics: Conference Series 2651, no. 1 (December 1, 2023): 012027. http://dx.doi.org/10.1088/1742-6596/2651/1/012027.
Full textAbstract:
Abstract We investigate the rigid pavement influence on phase velocity dispersion of the Rayleigh wave by using numerical modeling. An overestimation of Rayleigh wave phase velocity is observed more than 10% deviation caused by high-velocity rigid pavement. Based on the generalized reflection and transmission (R/T) coefficient method, the linear eigenvalue problem is solved for vertically heterogeneous media with a rigid pavement structure. The secular function demonstrates the multi-valuedness which is different phase velocity at the same frequency for the single Rayleigh mode. To address mode misidentification problem, we adopt a new misfit function for multimodal dispersion inversion. A Markov chain Monte Carlo (MCMC) sampling algorithm is used to solve dispersion inversion. The synthetic test shows that Rayleigh wave dispersion inversion incorporating a pavement layer can improve estimation accuracy of shear wave velocity.
39
Cui, Qing Yi, Ya Wei Wang, and Min Bu. "Study of Phase Function of the Biological Cell." Applied Mechanics and Materials 433-435 (October 2013): 795–98. http://dx.doi.org/10.4028/www.scientific.net/amm.433-435.795.
Full textAbstract:
In order to settle the influence of measuring results of cells with its shape-variety in the examination of biological cell, MCEM (modify co-central ellipse model) is adopted based on the theories of Rayleigh-Debye-Gans, and the changing tendency of scattering phase function of nucleated cells with different body factor, nuclear size and shooting angle are studied. These provide a useful theory foundation for improving the measurement and distinguishing of biological cell.
40
Kattawar, George W. "Irradiance invariance for scattering according to a Rayleigh phase function compared to a Rayleigh phase matrix for a plane-parallel medium." Applied Optics 29, no. 16 (June 1, 1990): 2365. http://dx.doi.org/10.1364/ao.29.002365.
Full text
41
Krim, M. S. Abdel. "Stochastic radiative transfer in finite plane for Rayleigh scattering phase function." Journal of Quantitative Spectroscopy and Radiative Transfer 69, no. 6 (June 2001): 745–59. http://dx.doi.org/10.1016/s0022-4073(00)00110-2.
Full text
42
SCHUMACHER, JÖRG, and OLIVIER PAULUIS. "Buoyancy statistics in moist turbulent Rayleigh–Bénard convection." Journal of Fluid Mechanics 648 (April 7, 2010): 509–19. http://dx.doi.org/10.1017/s0022112010000030.
Full textAbstract:
We study shallow moist Rayleigh–Bénard convection in the Boussinesq approximation in three-dimensional direct numerical simulations. The thermodynamics of phase changes is approximated by a piecewise linear equation of state close to the phase boundary. The impact of phase changes on the turbulent fluctuations and the transfer of buoyancy through the layer is discussed as a function of the Rayleigh number and the ability to form liquid water. The enhanced buoyancy flux due to phase changes is compared with dry convection reference cases and related to the cloud cover in the convection layer. This study indicates that the moist Rayleigh–Bénard problem offers a practical framework for the development and evaluation of parameterizations for atmospheric convection.
43
Vita, A. N., Zulfakriza, A. A. Martha, S. Rohadi, N. Heryandoko, and C. Milkerreit. "Preliminary Result of Rayleigh Wave Tomography beneath Jailolo Volcanic Complex, North Moluccas, Indonesia using Ambient Noise." Journal of Physics: Conference Series 2243, no. 1 (June 1, 2022): 012024. http://dx.doi.org/10.1088/1742-6596/2243/1/012024.
Full textAbstract:
Abstract From November 2015 to February 2016, seismic swarm activity occurred around Jailolo volcano, a long-dormant and few-studied volcano located on Halmahera Island, North Moluccas, Indonesia. We perform seismic ambient noise tomography method to delineate the sub-surface seismic structure beneath Jailolo volcanic complex. We use the vertical seismogram data from 29 temporary short-period seismic stations of 7G Seismic Network from January to March 2017 observation. The 7G Seismic Network was collaborative research between the Indonesian Agency of Meteorology, Climatology, Geophysics (BMKG), Indonesia, and the German Research Centre for Geoscience (GFZ). We cross-correlated vertical components of pair stations to estimate Green’s Rayleigh wave function. More than 350 group velocities were estimated. We measure the dispersion characteristic of Rayleigh wave group velocity by applying the frequency-time analysis (FTAN) technique. We compare Green Function and the dispersion curve by applying different band-pass filters. The Green Function of Rayleigh wave is clearly seen in causal and/or acausal part by applied 0.01-0.5 Hz band-pass filter and then dispersion curve is picked within 2-10 s. The Rayleigh group velocities average is 0.6-0.9 km/s.
44
Kattawar, George W. "Irradiance invariance for scattering according to a Rayleigh phase function compared to a Rayleigh phase matrix for a plane-parallel medium: erratum." Applied Optics 30, no. 30 (October 20, 1991): 4288. http://dx.doi.org/10.1364/ao.30.004288.
Full text
45
Yang, Xinyu, Haijiang He, Jun Xu, Yikun Wei, and Hua Zhang. "Entropy Generation Rates in Two-Dimensional Rayleigh–Taylor Turbulence Mixing." Entropy 20, no. 10 (September 26, 2018): 738. http://dx.doi.org/10.3390/e20100738.
Full textAbstract:
Entropy generation rates in two-dimensional Rayleigh–Taylor (RT) turbulence mixing are investigated by numerical calculation. We mainly focus on the behavior of thermal entropy generation and viscous entropy generation of global quantities with time evolution in Rayleigh–Taylor turbulence mixing. Our results mainly indicate that, with time evolution, the intense viscous entropy generation rate s u and the intense thermal entropy generation rate S θ occur in the large gradient of velocity and interfaces between hot and cold fluids in the RT mixing process. Furthermore, it is also noted that the mixed changing gradient of two quantities from the center of the region to both sides decrease as time evolves, and that the viscous entropy generation rate ⟨ S u ⟩ V and thermal entropy generation rate ⟨ S θ ⟩ V constantly increase with time evolution; the thermal entropy generation rate ⟨ S θ ⟩ V with time evolution always dominates in the entropy generation of the RT mixing region. It is further found that a “smooth” function ⟨ S u ⟩ V ∼ t 1 / 2 and a linear function ⟨ S θ ⟩ V ∼ t are achieved in the spatial averaging entropy generation of RT mixing process, respectively.
46
Nayak, Avinash, and Clifford H. Thurber. "Using multicomponent ambient seismic noise cross-correlations to identify higher mode Rayleigh waves and improve dispersion measurements." Geophysical Journal International 222, no. 3 (June 1, 2020): 1590–605. http://dx.doi.org/10.1093/gji/ggaa270.
Full textAbstract:
SUMMARY Ambient seismic noise cross-correlation with three-component sensors yields a nine-component empirical Green's tensor, in which four components of the radial–vertical plane contain Rayleigh waves. We exploit the retrograde elliptical nature of particle motion of the fundamental mode Rayleigh wave to correct the phase of the four radial–vertical components and stack them to obtain an average fundamental mode Rayleigh-wave time-series. This technique can suppress incoherent noise and wave packets that do not follow the targeted elliptical particle motion. The same technique can be used to isolate the first higher mode Rayleigh wave that follows prograde elliptical particle motion. We first demonstrate the effectiveness of the method on synthetic waveforms and then apply it on noise cross-correlations computed in Central California. Using this method, we isolate 1st higher mode Rayleigh waves on noise cross-correlations in the Great Valley, California, which provides new phase velocity constraints for estimating velocity structure in the sedimentary basin. We also obtain improved estimates of fundamental mode Rayleigh-wave dispersion for surface-wave tomography. The waveforms stacked assuming retrograde particle motion return at least ∼20 per cent more group velocity dispersion measurements satisfying a minimum signal-to-noise ratio (SNR) criterion than the individual components for periods ∼4–18 s. For equivalent group velocity measurements, SNR for the stacked estimate of the fundamental mode Rayleigh wave is on average 40 per cent greater than that measured on the individual components at periods less than 10 s. The technique also provides an easy way to detect large errors in sensor orientation.
47
Dai, Chuanshan, and Hideo Inaba. "Neutral Instability and Optimum Convective Mode in a Fluid Layer with PCM Particles." Journal of Heat Transfer 127, no. 12 (June 10, 2005): 1289–95. http://dx.doi.org/10.1115/1.2060728.
Full textAbstract:
Linear stability analysis is performed to determine the critical Rayleigh number for the onset of convection in a fluid layer with phase-change-material particles. Sine and Gaussian functions are used for describing the large variation of apparent specific heat in a narrow phase changing temperature range. The critical conditions are numerically obtained using the fourth order Runge-Kutta-Gill finite difference method with Newton-Raphson iteration. The critical eigenfunctions of temperature and velocity perturbations are obtained. The results show that the critical Rayleigh number decreases monotonically with the amplitude of Sine or Gaussian function. There is a minimum critical Rayleigh number while the phase angle is between π∕2 and π, which corresponds to the optimum experimental convective mode.
48
Rozi Kurniawan, Muhammad Fachrul, Shindy Rosalia, Andri Dian Nugraha, Zulfakriza, David P. Sahara, Abdul Muhari, Andi Azhar Nurdin, et al. "Ambient Seismic Noise Cross - correlation of Ambon Island and Surrounding Area, Eastern Indonesia: Preliminary Result." IOP Conference Series: Earth and Environmental Science 873, no. 1 (October 1, 2021): 012023. http://dx.doi.org/10.1088/1755-1315/873/1/012023.
Full textAbstract:
Abstract The island of Ambon lies on complex tectonics, part of Banda Arc which is driven by the Australia – Eurasia collision. Historical earthquake data show that an earthquake resulting the greatest tsunami in Indonesia had occurred at Ambon Island. On 26 September 2019, Ambon was shaken by an M 6.5 earthquake at a depth of 10 km (BMKG). In this study, we use ambient noise data from 11 temporary stations deployed by ITB and 4 permanent stations owned BMKG which are recorded from October until December 2019. Here, we purely use the vertical component of seismogram to retrieve the Empirical Green’s Function of Rayleigh waves. Cross-correlations were obtained from the daily data series and stacked the day-by-day cross-correlation data into one inter-station cross-correlation. The Empirical Green’s Function is seen at the band period 1-15 s. As a part of our study, we analyze the Green’s Function with frequency-time analysis (FTAN) to get Rayleigh wave group velocity. The group velocity of Rayleigh waves varies from 1.04 km/s – 3.75 km/s. Low group velocity might be indicated the presence of sediment or volcanic deposits and high group velocity might be indicated metamorphic rocks. The result of this study might give a finer velocity model of the shallow crustal beneath Ambon Island and the surrounding area.
49
Pei, Donghong, John N. Louie, and Satish K. Pullammanappallil. "Application of simulated annealing inversion on high-frequency fundamental-mode Rayleigh wave dispersion curves." GEOPHYSICS 72, no. 5 (September 2007): R77—R85. http://dx.doi.org/10.1190/1.2752529.
Full textAbstract:
The simulated annealing (SA) inversion technique has been successfully applied for solving various nonlinear geophysical problems. Following previous developments, we modified the SA inversion, yielding 1D shallow S-wave velocity profiles from high frequency fundamental-mode Rayleigh dispersion curves, and validated the inversion with blind tests. Unlike previous applications of SA, this study draws random numbers from a standard Gaussian distribution. The numbers simultaneously perturb both S-wave velocities and the layer thickness of models. The annealing temperature is gradually decreased following a polynomial-time cooling schedule. Phase velocities are calculated using the reflectivity-transmission coefficient method. The reliability of the model resulting from our implementation is evaluated by statistically calculating the expected values of model parameters and their covariance matrices. Blind tests on two field and 12 synthetic Rayleigh dispersion data sets show that our SA implementation works well for S-wave velocity inversion of dispersion curves from high-frequency fundamental-mode Rayleigh waves. Blind estimates of layer S-wave velocities fall within one standard deviation of the velocities of the original synthetic models in 78% of cases.
50
Li, Lun, and Yuanyuan V. Fu. "Surface-Wave Tomography of Eastern and Central Tibet from Two-Plane-Wave Inversion: Rayleigh-Wave and Love-Wave Phase Velocity Maps." Bulletin of the Seismological Society of America 110, no. 3 (March 17, 2020): 1359–71. http://dx.doi.org/10.1785/0120190199.
Full textAbstract:
ABSTRACT An understanding of mantle dynamics occurring beneath the Tibetan plateau requires a detailed image of its seismic velocity and anisotropic structure. Surface waves at long periods (>50 s) could provide such critical information. Though Rayleigh-wave phase velocity maps have been constructed in the Tibetan regions using ambient-noise tomography (ANT) and regional earthquake surface-wave tomography, Love-wave phase velocity maps, especially those at longer periods (>50 s), are rare. In this study, two-plane-wave teleseismic surface-wave tomography is applied to develop 2D Rayleigh-wave and Love-wave phase velocity maps at periods between 20 and 143 s across eastern and central Tibet and its surroundings using four temporary broadband seismic experiments. These phase velocity maps share similar patterns and show high consistency with those previously obtained from ANT at overlapping periods (20–50 s), whereas our phase velocity maps carry useful information at longer periods (50–143 s). Prominent slow velocity is imaged at periods of 20–143 s beneath the interior of the Tibetan plateau (i.e., the Songpan–Ganzi terrane, the Qiangtang terrane, and the Lhasa terrane), implying the existence of thick Tibetan crust along with warm and weak Tibetan lithosphere. In contrast, the dispersal of fast velocity anomalies coincides with mechanically strong, cold tectonic blocks, such as the Sichuan basin and the Qaidam basin. These phase velocity maps could be used to construct 3D shear-wave velocity and radial seismic anisotropy models of the crust and upper mantle down to 250 km across the eastern and central Tibetan plateau.