(back to Carl Tape home)

A geometric perspective on seismic moment tensors

GALLERY OF IMAGES

Go here to see some example figures related to the papers above.

----------------------------------------------------------------

A geometric perspective on seismic moment tensors

GALLERY OF IMAGES

A geometric setting for moment tensors (pdf) Geophysical Journal International, v. 190, p. 476-498, 2012a Walter Tape and Carl Tape We describe a parametrization of
moment tensors that is suitable for use in algorithms
for moment tensor inversion. The parameters are
conceptually natural and can be easily visualized. The
ingredients of the parametrization are present in the
literature; we have consolidated them into a concise
statement in a geometric setting. We treat several
familiar moment tensor topics in the same geometric
setting as well. These topics include moment tensor
decompositions, crack-plus-double-couple moment tensors,
and the parameter that measures the difference between
a deviatoric moment tensor and a double couple. The
geometric approach clarifies concepts that are sometimes
obscured by calculations.
A geometric comparison of source-type plots for moment tensors (pdf) Geophysical Journal International, v. 190, p. 499-510, 2012b Walter Tape and Carl Tape We describe a general geometric
framework for thinking about source-type plots for
moment tensors. We consider two fundamental examples,
one where the source-type plot is on the unit sphere,
and one where it is on the unit cube. The plot on the
sphere is preferable to the plot on the cube: it is
simpler, it embodies a more natural assumption about
eigenvalue probabilities and it is more consistent with
the conventional euclidian definition of scalar seismic
moment. We describe the source-type plots of Hudson,
Pearce and Rogers (1989) in our geometric context, and
we find that they are equivalent to a plot on the cube.
We therefore suggest the plot on the sphere as an
alternative.
Angle between principal axis triples (pdf)Geophysical Journal International, v. 191, p. 813-831, 2012c Walter Tape and Carl Tape The principal axis angle xi0, or
Kagan angle, is a measure of the difference between the
orientations of two seismic moment tensors. It is the
smallest angle needed to rotate the principal axes of
one moment tensor to the corresponding principal axes of
the other. We give a concise formula for calculating
xi0, but our main goal is to illustrate the behavior of
xi0 geometrically. When the first of two moment tensors
is fixed, the angle xi0 between them becomes a function
on the unit ball. The level surfaces of xi0 can then be
depicted in the unit ball, and they give insights into
xi0 that are not obvious from calculations alone. We
also include a derivation of the known probability
density of xi0. The density is proportional to the area
of a certain surface, whose variation explains the
rather peculiar shape of the probability density curve.
This shape needs to be considered when analyzing
distributions of empirical xi0 values. We recall an
example of Willemann which shows that xi0 may not always
be the most appropriate measure of separation for moment
tensor orientations, and we offer an alternative
measure.
The classical model for moment tensors (pdf, supplement)Geophysical Journal International, v. 195, p. 1701-1720, 2013 Walter Tape and Carl Tape A seismic moment tensor is a description of an earthquake source, but the description is indirect. The moment tensor describes seismic radiation rather than the actual physical process that initiates the radiation. A moment tensor ‘model’ then ties the physical process to the moment tensor. The model is not unique, and the physical process is therefore not unique. In the classical moment tensor model, an earthquake arises from slip along a planar fault, but with the slip not necessarily in the plane of the fault. The model specifies the resulting moment tensor in terms of the slip vector, the fault normal vector and the Lame elastic parameters, assuming isotropy. We review the classical model in the context of the fundamental lune. The lune is closely related to the space of moment tensors, and it provides a setting that is conceptually natural as well as pictorial. In addition to the classical model, we consider a crack plus double-couple model (CDC model) in which a moment tensor is regarded as the sum of a crack tensor and a double couple. A uniform parameterization of moment tensors (pdf, supplement)Geophysical Journal International, 2015 Walter Tape and Carl Tape A moment tensor is a 3 x 3 symmetric matrix that expresses an earthquake source. We construct a parameterization of the five-dimensional space of all moment tensors of unit norm. The coordinates associated with the parameterization are closely related to moment tensor orientations and source types. The parameterization is uniform, in the sense that equal volumes in the coordinate domain of the parameterization correspond to equal volumes of moment tensors. Uniformly distributed points in the coordinate domain therefore give uniformly distributed moment tensors. A cartesian grid in the coordinate domain can be used to search efficiently over moment tensors. We find that uniformly distributed moment tensors have uniformly distributed orientations (eigenframes), but that their source types (eigenvalue triples) are distributed so as to favor double couples. A confidence parameter for seismic moment tensors
(pdf, supplement)Geophysical Journal International, 2016 Walter Tape and Carl Tape Given a moment tensor m inferred from seismic data for an earthquake, we define Ps(V) to be the probability that the true moment tensor for the earthquake lies in the neighborhood of m that has fractional volume V. The average value of Ps(V) is then a measure of our confidence in m. The calculation of Ps(V) requires knowing both the probability P(ω) and the fractional volume V(ω) of the set of moment tensors within a given angular radius ω of m. We explain how to construct Ps(ω) from a misfit function derived from seismic data, and we show how to calculate V(ω), which depends on the set M of moment tensors under consideration. The two most important instances of M are where M is the set of all moment tensors of fixed norm, and where M is the set of all double couples of fixed norm. ------------------------------------------------------------------------------ Chapman-Leaney decompositions (pdf)This pdf is a set of notes we wrote regarding a 2012 GJI paper by Chapman and Leaney. Our notes were published as an on-line supplement of a 2014 GJI addendum paper by Chapman and Leaney. In our notes we use the lune in the context of general anisotropy to investigate the moment tensor decomposition proposed by Chapman and Leaney. |

I have posted online some
matlab scripts, GMT scripts, and mathematica scripts
relevant for interpreting and plotting moment tensors. These
files are available within the compearth
repository on github. They can be downloaded from github as a zipped file or using git from the command line. Please email me with corrections or suggestions. If you think you might use some of these concepts, here is a start-up guide:1. Check out this pdf. 2. Download the compearth repository from github (git clone https://github.com/carltape/compearth) 3. Follow the setup instructions in the file compearth/momenttensor/matlab/README (here) 4. If you have Matlab installed, try one of the scripts to reproduce results in one of our published papers: TT2012kaganAppE.m, TT2013AppA.m, TT2015AppA.m These scripts use several functions that you may find useful. 5. If you have GMT installed, try out lune.pl from compearth/momenttensor/plot/gmt/lune.pl (Please email me if things do not work for you.) Here are some notes about the use of percentages in moment tensor decompositions. |

Go here to see some example figures related to the papers above.

----------------------------------------------------------------