An official website of the United States government
An elastic map T associates stress with strain in some material. A symmetry of T is a rotation of the material that leaves T unchanged, and the symmetry group of T consists of all such rotations. The symmetry class of T describes the symmetry group but without the orientation information. With an eye toward geophysical applications, Browaeys & Chevrot developed a theory which, for any elastic map T and for each of six symmetry classes Σ, computes the "Σ-percentage" of T. The theory also finds a "hexagonal approximation"—an approximation to T whose symmetry class is at least transverse isotropic. We reexamine their theory and recommend that the Σ-percentages be abandoned. We also recommend that the hexagonal approximations to T be replaced with the closest transverse isotropic maps to T.
more »
« less
Navigating the space of seismic anisotropy for crystal and whole-Earth scales
SUMMARY Evidence of seismic anisotropy is widespread within the Earth, including from individual crystals, rocks, borehole measurements, active-source seismic data, and global seismic data. The seismic anisotropy of a material determines how wave speeds vary as a function of propagation direction and polarization, and it is characterized by density and the elastic map, which relates strain and stress in the material. Associated with the elastic map is a symmetric $$6 \times 6$$ matrix, which therefore has 21 parameters. The 21-D space of elastic maps is vast and poses challenges for both theoretical analysis and typical inverse problems. Most estimation approaches using a given set of directional wave speed measurements assume a high-symmetry approximation, typically either in the form of isotropy (2 parameters), vertical transverse isotropy (radial anisotropy: 5 parameters), or horizontal transverse isotropy (azimuthal anisotropy: 6 parameters). We offer a general approach to explore the space of elastic maps by starting with a given elastic map $$\mathbf {T}$$. Using a combined minimization and projection procedure, we calculate the closest $$\Sigma$$-maps to $$\mathbf {T}$$, where $$\Sigma$$ is one of the eight elastic symmetry classes: isotropic, cubic, transverse isotropic, trigonal, tetragonal, orthorhombic, monoclinic and trivial. We apply this approach to 21-parameter elastic maps derived from laboratory measurements of minerals; the measurements include dependencies on pressure, temperature, and composition. We also examine global elasticity models derived from subduction flow modelling. Our approach offers a different perspective on seismic anisotropy and motivates new interpretations, such as for why elasticity varies as a function of pressure, temperature, and composition. The two primary advances of this study are (1) to provide visualization of elastic maps, including along specific pathways through the space of model parameters, and (2) to offer distinct options for reducing the complexity of a given elastic map by providing a higher-symmetry approximation or a lower-anisotropic version. This could contribute to improved imaging and interpretation of Earth structure and dynamics from seismic anisotropy.
more »
« less
- PAR ID:
- 10589990
- Publisher / Repository:
- Oxford University Press
- Date Published:
- Journal Name:
- Geophysical Journal International
- Volume:
- 242
- Issue:
- 1
- ISSN:
- 0956-540X
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
Universal displacements are those displacements that can be maintained, in the absence of body forces, by applying only boundary tractions for any material in a given class of materials. Therefore, equilibrium equations must be satisfied for arbitrary elastic moduli for a given anisotropy class. These conditions can be expressed as a set of partial differential equations for the displacement field that we call universality constraints. The classification of universal displacements in homogeneous linear elasticity has been completed for all the eight anisotropy classes. Here, we extend our previous work by studying universal displacements in inhomogeneous anisotropic linear elasticity assuming that the directions of anisotropy are known. We show that universality constraints of inhomogeneous linear elasticity include those of homogeneous linear elasticity. For each class and for its known universal displacements, we find the most general inhomogeneous elastic moduli that are consistent with the universality constrains. It is known that the larger the symmetry group, the larger the space of universal displacements. We show that the larger the symmetry group, the more severe the universality constraints are on the inhomogeneities of the elastic moduli. In particular, we show that inhomogeneous isotropic and inhomogeneous cubic linear elastic solids do not admit universal displacements and we completely characterize the universal inhomogeneities for the other six anisotropy classes.more » « less
-
Abstract Comprehensive observations of surface wave anisotropy across Alaska and the Aleutian subduction zone would help to improve understanding of its tectonics, mantle dynamics, and earthquake risk. We produce such observations, using stations from the USArray Transportable Array, regional networks across Alaska, and the Alaska Amphibious Community Seismic Experiment in the Alaska‐Aleutian subduction zone both onshore and offshore. Our data include Rayleigh and Love wave phase dispersion from earthquakes (28–85 s) and ambient noise two‐ and three‐station interferometry (8–50 s). Compared with using two‐station interferometry alone, three‐station interferometry significantly improves the signal‐to‐noise ratio and approximately doubles the number of measurements retained. Average differences between both isotropic and anisotropic tomographic maps constructed from different methods lie within their uncertainties, which is justification for combining the measurements. The composite tomographic maps include Rayleigh wave isotropy and azimuthal anisotropy from 8 to 85 s both onshore and offshore, and onshore Love wave isotropy from 8 to 80 s. In the Alaska‐Aleutian subduction zone, Rayleigh wave fast directions vary from trench parallel to perpendicular and back to parallel with increasing periods, apparently reflecting the effect of the subducted Pacific Plate. The tomographic maps provide a basis for inferring the 3‐D anisotropic crustal and uppermost mantle structure across Alaska and the Aleutian subduction zone.more » « less
-
Abstract An elastic map $$\mathbf {T}$$ T describes the strain-stress relation at a particular point $$\mathbf {p}$$ p in some material. A symmetry of $$\mathbf {T}$$ T is a rotation of the material, about $$\mathbf {p}$$ p , that does not change $$\mathbf {T}$$ T . We describe two ways of inferring the group $$\mathcal {S} _{ \mathbf {T} }$$ S T of symmetries of any elastic map $$\mathbf {T}$$ T ; one way is qualitative and visual, the other is quantitative. In the first method, we associate to each $$\mathbf {T}$$ T its “monoclinic distance function” "Equation missing" on the unit sphere. The function "Equation missing" is invariant under all of the symmetries of $$\mathbf {T}$$ T , so the group $$\mathcal {S} _{ \mathbf {T} }$$ S T is seen, approximately, in a contour plot of "Equation missing" . The second method is harder to summarize, but it complements the first by providing an algorithm to compute the symmetry group $$\mathcal {S} _{ \mathbf {T} }$$ S T . In addition to $$\mathcal {S} _{ \mathbf {T} }$$ S T , the algorithm gives a quantitative description of the overall approximate symmetry of $$\mathbf {T}$$ T . Mathematica codes are provided for implementing both the visual and the quantitative approaches.more » « less
-
A theoretical approach was recently introduced [Guidetti and Royston, J. Acoust. Soc. Am. 144, 2312–2323 (2018)] for the radially converging slow shear wave pattern in transverse isotropic materials subjected to axisymmetric excitation normal to the axis of isotropy at the outer boundary of the material. This approach is enabled via transformation to an elliptic coordinate system with isotropic properties. The approach is extended to converging fast shear waves driven by axisymmetric torsional motion polarized in a plane containing the axis of isotropy. The approach involves transformation to a super-elliptic shape with isotropic properties and use of a numerically efficient boundary value approximation.more » « less
