An official website of the United States government
Abstract Seismic anisotropy provides essential information for characterizing the orientation of deformation and flow in the crust and mantle. The isotropic structure of the Antarctic crust and upper mantle has been determined by previous studies, but the azimuthal anisotropy structure has only been constrained by mantle core phase (SKS) splitting observations. This study determines the azimuthal anisotropic structure of the crust and mantle beneath the central and West Antarctica based on 8—55 s Rayleigh wave phase velocities from ambient noise cross‐correlation. An anisotropic Rayleigh wave phase velocity map was created using a ray—based tomography method. These data are inverted using a Bayesian Monte Carlo method to obtain an azimuthal anisotropy model with uncertainties. The azimuthal anisotropy structure in most of the study region can be fit by a two‐layer structure, with one layer at depths of 0–15 km in the shallow crust and the other layer in the uppermost mantle. The azimuthal anisotropic layer in the shallow crust of West Antarctica, where it coincides with strong positive radial anisotropy quantified by the previous study, shows a fast direction that is subparallel to the inferred extension direction of the West Antarctic Rift System. Fast directions of upper mantle azimuthal anisotropy generally align with teleseismic shear wave splitting fast directions, suggesting a thin lithosphere or similar lithosphere‐asthenosphere deformation. However, inconsistencies in this exist in Marie Byrd Land, indicating differing ancient deformation patterns in the shallow mantle lithosphere sampled by the surface waves and deformation in the deeper mantle and asthenosphere sampled more strongly by splitting measurements.
more »
« less
Spurious Rayleigh-wave apparent anisotropy near major structural boundaries: a numerical and theoretical investigation
SUMMARY The recent developments in array-based surface-wave tomography have made it possible to directly measure apparent phase velocities through wave front tracking. While directionally dependent measurements have been used to infer intrinsic $$2\psi $$ azimuthal anisotropy (with a 180° periodicity), a few studies have also demonstrated strong but spurious $$1\psi $$ azimuthal anisotropy (360° periodicity) near major structure boundaries particularly for long period surface waves. In such observations, Rayleigh waves propagating in the direction perpendicular to the boundary from the slow to the fast side persistently show a higher apparent velocity compared to waves propagating in the opposite direction. In this study, we conduct numerical and theoretical investigations to explore the effect of scattering on the apparent Rayleigh-wave phase velocity measurement. Using 2-D spectral-element numerical wavefield simulations, we first reproduce the observation that waves propagating in opposite directions show different apparent phase velocities when passing through a major velocity contrast. Based on mode coupling theory and the locked mode approximation, we then investigate the effect of the scattered fundamental-mode Rayleigh wave and body waves interfering with the incident Rayleigh wave separately. We show that scattered fundamental-mode Rayleigh waves, while dominating the scattered wavefield, mostly cause short wavelength apparent phase velocity variations that could only be studied if the station spacing is less than about one tenth of the surface wave wavelength. Scattered body waves, on the other hand, cause longer wavelength velocity variations that correspond to the existing real data observations. Because of the sensitivity of the $$1\psi $$ apparent anisotropy to velocity contrasts, incorporating such measurements in surface wave tomography could improve the resolution and sharpen the structural boundaries of the inverted model.
more »
« less
- Award ID(s):
- 1753362
- PAR ID:
- 10540859
- Publisher / Repository:
- Oxford University Press
- Date Published:
- Journal Name:
- Geophysical Journal International
- Volume:
- 239
- Issue:
- 2
- ISSN:
- 0956-540X
- Format(s):
- Medium: X Size: p. 901-913
- Size(s):
- p. 901-913
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
Abstract We use Eikonal tomography to derive phase and group velocities of surface waves for the plate boundary region in Southern California. Seismic noise data in the period range 2 and 20 s recorded in year 2014 by 346 stations with ~1‐ to 30‐km station spacing are analyzed. Rayleigh and Love wave phase travel times are measured using vertical‐vertical and transverse‐transverse noise cross correlations, and group travel times are derived from the phase measurements. Using the Eikonal equation for each location and period, isotropic phase and group velocities and 2‐psi azimuthal anisotropy are determined statistically with measurements from different virtual sources. Starting with the SCEC Community Velocity Model, the observed 2.5‐ to 16‐s isotropic phase and group dispersion curves are jointly inverted on a 0.05° × 0.05° grid to obtain local 1‐D piecewise shear wave velocity (Vs) models. Compared to the starting model, the final results have generally lowerVsin the shallow crust (top 3–10 km), particularly in areas such as basins and fault zones. The results also show clear velocity contrasts across the San Andreas, San Jacinto, Elsinore, and Garlock Faults and suggest that the San Andreas Fault southeast of San Gorgonio Pass is dipping to the northeast. Investigation of the nonuniqueness of the 1‐DVsinversion suggests that imaging the top 3‐kmVsstructure requires either shorter period (≤2 s) surface wave dispersion measurements or other types of data set such as Rayleigh wave ellipticity.more » « less
-
Summary For a weakly anisotropic medium, Rayleigh and Love wave phase speeds at angular frequency ω and propagation azimuth ψ are given approximately by V(ω, ψ) = A0 + A2ccos 2ψ + A2ssin 2ψ + A4ccos 4ψ + A4ssin 4ψ. Earlier theories of the propagation of surface waves in anisotropic media based on non-degenerate perturbation theory predict that the dominant components are expected to be 2ψ for Rayleigh waves and 4ψ for Love waves. This paper is motivated by recent observations of the the 2ψ component for Love waves and 4ψ for Rayleigh waves, referred to here as “unexpected anisotropy”. To explain these observations, we present a quasi-degenerate theory of Rayleigh-Love coupling in a weakly anisotropic medium based on Hamilton’s Principle in Cartesian coordinates, benchmarking this theory with numerical results based on SPECFEM3D. We show that unexpected anisotropy is expected to be present when Rayleigh-Love coupling is strong and recent observations of Rayleigh and Love wave 2ψ and 4ψ anisotropy can be fit successfully with physically plausible models of a depth-dependent tilted transversely isotropic (TTI) medium. In addition, when observations of the 2ψ and 4ψ components of Rayleigh and Love anisotropy are used in the inversion, the ellipticity parameter ηX, introduced here, is better constrained, we can constrain the absolute dip direction based on polarization measurements, and we provide evidence that the mantle should be modeled as a tilted orthorhombic medium rather than a TTI medium. Ignoring observations of unexpected anisotropy may bias the estimated seismic model significantly. We also provide information about the polarization of the quasi-Love waves and coupling between fundamental mode Love and overtone Rayleigh waves in both continental and oceanic settings. The theory of SV-SH coupling for horizontally propagating body waves is presented for comparison with the surface wave theory, with emphasis on results for a TTI medium.more » « less
-
null (Ed.)Two types of surface wave anisotropy are observed regularly by seismologists but are only rarely interpreted jointly: apparent radial anisotropy, which is the difference in propagation speed between horizontally and vertically polarized waves inferred from Love and Rayleigh waves, and apparent azimuthal anisotropy, which is the directional dependence of surface wave speeds (usually Rayleigh waves). We show that a new data set of Love and Rayleigh wave isotropic phase speeds and Rayleigh wave azimuthal anisotropy observed within and surrounding eastern Tibet can be explained simultaneously by modeling the crust as a depth-dependent tilted hexagonally symmetric (THS) medium. We specify the THS medium with depth-dependent hexagonally symmetric elastic tensors tilted and rotated through dip and strike angles and estimate these quantities using a Bayesian Monte Carlo inversion to produce a 3-D model of the crust and uppermost mantle on a 0.5° × 0.5° spatial grid. In the interior of eastern Tibet and in the Yunnan-Guizhou plateau, we infer a steeply dipping THS upper crustal medium overlying a shallowly dipping THS medium in the middle-to-lower crust. Such vertical stratification of anisotropy may reflect a brittle to ductile transition in which shallow fractures and faults control upper crustal anisotropy and the crystal-preferred orientation of anisotropic (perhaps micaceous) minerals governs the anisotropy of the deeper crust. In contrast, near the periphery of the Tibetan Plateau the anisotropic medium is steeply dipping throughout the entire crust, which may be caused by the reorientation of the symmetry axes of deeper crustal anisotropic minerals as crustal flows are rotated near the borders of Tibet.more » « less
-
SUMMARY Global variations in the propagation of fundamental-mode and overtone surface waves provide unique constraints on the low-frequency source properties and structure of the Earth’s upper mantle, transition zone and mid mantle. We construct a reference data set of multimode dispersion measurements by reconciling large and diverse catalogues of Love-wave (49.65 million) and Rayleigh-wave dispersion (177.66 million) from eight groups worldwide. The reference data set summarizes measurements of dispersion of fundamental-mode surface waves and up to six overtone branches from 44 871 earthquakes recorded on 12 222 globally distributed seismographic stations. Dispersion curves are specified at a set of reference periods between 25 and 250 s to determine propagation-phase anomalies with respect to a reference Earth model. Our procedures for reconciling data sets include: (1) controlling quality and salvaging missing metadata; (2) identifying discrepant measurements and reasons for discrepancies; (3) equalizing geographic coverage by constructing summary rays for travel-time observations and (4) constructing phase velocity maps at various wavelengths with combination of data types to evaluate inter-dataset consistency. We retrieved missing station and earthquake metadata in several legacy compilations and codified scalable formats to facilitate reproducibility, easy storage and fast input/output on high-performance-computing systems. Outliers can be attributed to cycle skipping, station polarity issues or overtone interference at specific epicentral distances. By assessing inter-dataset consistency across similar paths, we empirically quantified uncertainties in traveltime measurements. More than 95 per cent measurements of fundamental-mode dispersion are internally consistent, but agreement deteriorates for overtones especially branches 5 and 6. Systematic discrepancies between raw phase anomalies from various techniques can be attributed to discrepant theoretical approximations, reference Earth models and processing schemes. Phase-velocity variations yielded by the inversion of the summary data set are highly correlated (R ≥ 0.8) with those from the quality-controlled contributing data sets. Long-wavelength variations in fundamental-mode dispersion (50–100 s) are largely independent of the measurement technique with high correlations extending up to degree ∼25. Agreement degrades with increasing branch number and period; highly correlated structure is found only up to degree ∼10 at longer periods (T > 150 s) and up to degree ∼8 for overtones. Only 2ζ azimuthal variations in phase velocity of fundamental-mode Rayleigh waves were required by the reference data set; maps of 2ζ azimuthal variations are highly consistent between catalogues ( R = 0.6–0.8). Reference data with uncertainties are useful for improving existing measurement techniques, validating models of interior structure, calculating teleseismic data corrections in local or multiscale investigations and developing a 3-D reference Earth model.more » « less
