Skip to main content

Seismic moment and volume change of a spherical source

Abstract

A spherical source, one of the simplest seismic sources, has been represented in various ways in the literature. These representations include a spherical source with outward radial expansion (S1), a spherical crack source with outward and inward crack wall motions along the spherical surface (S2), an isotropic source represented by three mutually perpendicular vector dipoles (S3) or three mutually perpendicular tensile cracks (S4), and a spherical source undergoing a transformational expansion (S5). We systematically examined these sources and their static displacement fields to clarify how these representations are mutually related. We also considered the sources in a bimaterial medium, in which the source material is different from the surrounding medium, as a model of a magma or hydrothermal reservoir. Our examinations show that the source volume change of a spherical source (S1) (actual volume) can be uniquely determined from the seismic moment of an isotropic source (S3) regardless of our assumption of the source medium and that the actual volume of S1 is related to the seismic moment of S3 through the equivalence of the displacement fields due to these two sources. The seismic moment of S3 is also related through another equation to the source volume change of three tensile cracks (S4), which is equal to the source volume changes defined in S2 and S5. This relation has different forms, depending on the source medium and source process. This study provides a unified view for quantifying a spherical source using the seismic moment of an isotropic source determined from waveform inversion analysis.

Background

Spherical and cylindrical volume changes have been considered as sources of volcano seismic signals (e.g., Kanamori et al. 1980; Fukuyama and Takeo 1987; Uhira and Takeo 1994; Dreger et al. 2000; Legrand et al. 2000; Tameguri et al. 2002; Ohminato et al. 2006; Maeda and Takeo 2011; Kumagai et al. 2011). Although a spherical source is one of the simplest seismic sources, there are various ways to represent it. These representations include a spherical source with outward radial expansion (e.g., Mogi 1958) (S1, Figure 1a), a spherical crack source of Müller (2001) (S2, Figure 1b), an isotropic source represented by three mutually perpendicular vector dipoles (S3, Figure 1c) or three mutually perpendicular tensile cracks (S4, Figure 1d), and a spherical source undergoing a transformational expansion based on the concept of imaginary cutting, straining, and welding operations of Eshelby (1957) as introduced by Aki and richards (2002), which we call an Eshelby spherical source (S5, Figure 2a).

Figure 1
figure 1

Source representations. (a) A spherical source with purely outward radial expansion of Δs (S1), (b) a spherical crack source with outward and inward displacements (Δs and ds, respectively) along the spherical surface (S2), (c) an isotropic source with three mutually perpendicular vector dipoles (S3), and (d) a superposition of three mutually perpendicular tensile cracks (S4).

Figure 2
figure 2

Spherical source with stress-free strain and strain-free stress. (a) Schematic view of a spherical source undergoing a transformational expansion (S5) based on imaginary cutting, straining, and welding operations of Eshelby (1957). (b) A spherical source in which a strain-free stress or pressure (δ P) is applied to the removed source volume (S6) (see text for details).

Using these representations, the seismic moment of a spherical source has been defined in two different ways (Müller 2001; Richards and Kim 2005). The two definitions are

M 0 =(λ+2μ)Δ V s ,
(1)
M 0 =(λ+2μ/3)Δ V f ,
(2)

where λ and μ are Lamé’s constants and ΔVs and ΔVf are two different definitions of source volume changes (see below). The recent volcano seismology literature has used Equation 1 (e.g., Kawakatsu and Yamamoto 2007; Kumagai 2009; Chouet 2013), whereas earlier studies generally used Equation 2 (e.g., Chouet 1996; Julian et al. 1998; Dreger et al. 2000; Nishimura and Iguchi 2011).

According to Kawakatsu and Yamamoto (2007), the volume change in Equation 2Vf) corresponds to the stress-free volume, introduced by Eshelby (1957), whereas the actual volume in Equation 1Vs) may be smaller than the stress-free volume because of the confining pressure of the surrounding medium. Kawakatsu and Yamamoto (2007) also noted that the actual volume corresponds to the volume change used in the Mogi model (e.g., Mogi 1958).

Müller (2001) noted that the discrepancy between Equations 1 and 2 could not be resolved, and he interpreted these two forms to correspond to the limits of the actual volume of the source. Wielandt (2003) discussed the discrepancy and noted that the source volume change depends on the source geometry. Richards and Kim (2005) argued that each relationship is based on a different definition of volume changes at the source and that Equation 1 is preferred for characterizing underground explosions.

It is not clear why a spherical source has been represented in different ways and how these representations, including the two seismic moment equations, are mutually related. In Eshelby’s operations and Müller’s spherical crack, the source material is assumed to be the same as the surrounding medium. In volcanic regions, sources may be filled with magmatic or hydrothermal fluids. No examinations have been made on a spherical source in such bimaterial media, although the source representation of fault slip on a bimaterial interface has been discussed in various studies (e.g., Ampuero and Dahlen 2005).

In this paper, we first review the spherical source representations mentioned above and then derive the static displacement fields due to these sources in an infinite medium and in a half-space. We compare the analytical forms of these displacements to clarify the relationships among these representations. We further examine a spherical source in a bimaterial medium, which is of fundamental importance as a model of a magma or hydrothermal reservoir.

Methods

Source representations

Müller (2001) proposed to define a spherical source in terms of the outward crack wall motion of a spherical crack. Let us first define a spherical crack source (S2). We consider the moment density tensor of a tensile crack with normal direction (θ,ϕ) (Figure 3a):

Figure 3
figure 3

Source coordinates for (a) a crack and (b) a spherical source in a half-space.

m (θ,ϕ)= D s λ + 2 μ sin 2 θ cos 2 ϕ 2 μ sin 2 θ sin ϕ cos ϕ 2 μ sin θ cos θ cos ϕ 2 μ sin 2 θ sin ϕ cos ϕ λ + 2 μ sin 2 θ sin 2 ϕ 2 μ sin θ cos θ sin ϕ 2 μ sin θ cos θ cos ϕ 2 μ sin θ cos θ sin ϕ λ + 2 μ cos 2 θ ,
(3)

where Ds is the displacement discontinuity of a tensile crack (e.g., Chouet 1996). We consider a spherical crack surface with radius R where a constant radial displacement discontinuity Ds=(dss) occurs (Figure 1b). Here, the inner wall of the crack moves inward by ds, and the outer wall moves outward by Δs. The moment tensor of the spherical crack source can be obtained by integrating the moment density tensor of the tensile crack over the surface with radius R (Figure 1b):

M = 0 2 π 0 π m ( θ , ϕ ) R 2 sin θ d θ d ϕ , = Δ V D λ + 2 μ / 3 0 0 0 λ + 2 μ / 3 0 0 0 λ + 2 μ / 3 ,
(4)

where ΔVD is the volume given as

Δ V D =4π R 2 D s =4π R 2 d s + Δ s .
(5)

By solving the elastostatic boundary-value problem (see Subsection 'Displacement fields’), we obtain

Δ V s Δ V D = λ + 2 μ / 3 λ + 2 μ ,
(6)

where ΔVs = 4π R2Δs.

An isotropic source is defined as three mutually perpendicular vector dipoles (S3, Figure 1c), of which the moment tensor is given as

M = M 0 Iso 1 0 0 0 1 0 0 0 1 ,
(7)

where M 0 Iso is the seismic moment for an isotropic source. The three mutually perpendicular vector dipoles can be represented as the superposition of three mutually perpendicular tensile cracks (S4, Figure 1d). In such a case, M 0 Iso can be defined as

M 0 Iso =(λ+2μ/3)Δ V i .
(8)

Here, ΔVi = 3ΔVc, where ΔVc is the incremental change in the volume of each crack.

Finally, we consider an Eshelby spherical source following Aki and richards (2002) (S5, Figure 2a). A spherical volume with radius R removed from its surroundings expands to Ra with stress-free volumetric strain. We apply the surface traction that will restore the volume source to its original radius R, put the volume back in its hole, and weld the material across the cut. We release the applied traction, which results in the source volume expanding to Rs. The volume defined as ΔVf = 4π R2Δa is the stress-free volume, and the volume defined as ΔVs = 4π R2Δs is the actual volume or Mogi volume. We define the seismic moment using the stress-free volume as

M 0 =(λ+2μ/3)Δ V f ,
(9)

and the seismic moment using the Mogi volume may be given as

M 0 =(λ+2μ)Δ V s .
(10)

Displacement fields

Let us consider the static displacement fields due to each of the described sources. For the spherical source (S1), spherical crack source (S2), and Eshelby spherical source (S5), the displacement fields in an infinite medium are described by the following static equilibrium equation:

2 u r r 2 + 2 r u r ∂r - 2 r 2 u r =0,
(11)

where ur is the radial displacement. We denote the regions inside and outside the sphere with radius R as regions 1 and 2, respectively. Equation 11 has two solutions: ur = A r and ur = B/r2, where A and B are constants. The former is the interior solution for region 1 (r ≤ R), and the latter is the exterior solution for region 2 (r ≥ R).

Spherical source

We first consider the spherical source (S1). If the pressure in the spherical source increases by ΔP, the spherical source surface moves outward by Δs at r = R. Then, the boundary conditions at r = R are given as

u r e = Δ s
(12)
σ rr e = - ΔP ,
(13)

where u r e and σ rr e are the radial displacement and radial stress in region 2, respectively. Using these boundary conditions and the exterior solution ur = B/r2, we find that

ΔP=4μ Δ s R ,
(14)

and

u r e = Δ V s 4 π r 2 ,
(15)

where ΔVs = 4π R2Δs.

Following Sezawa (1931), we consider the spherical source in a half-space (Figure 3b). The radial displacement on the free surface at distance r is given as

u r e = ( λ + 2 μ ) 2 μ ( λ + μ ) R 3 ΔP r 2 ,
(16)
= ( λ + 2 μ ) 2 π ( λ + μ ) Δ V s r 2 .
(17)

If we assume λ = μ, the vertical and radial components for u r e in cylindrical coordinates are given as

u r e cos θ = 3 R 3 ΔP 4 μ c d 2 + c 2 3 / 2 ,
(18)
u r e sin θ = 3 R 3 ΔP 4 μ d d 2 + c 2 3 / 2 ,
(19)

respectively, where c is the source depth and d is the horizontal source distance (Figure 3b). Equations 18 and 19 represent the displacement field of the Mogi model (Mogi 1958).

Comparing Equation 17 for a half-space with Equation 15 for an infinite medium, we find that these two equations differ by a factor of 2(λ + 2μ)/(λ + μ) (Figure 4a). For this consideration, we used ΔVs = 4π R2Δs, in which neither the stress-free volume nor Eshelby’s operations are introduced.

Figure 4
figure 4

Static displacement fields. These fields are due to (a) a spherical source (S1) and (b) an isotropic source (S4) in an infinite medium and in a half-space.

Spherical crack source

Next, we consider the spherical crack source (S2) of Müller (2001) in an infinite medium. The boundary conditions at r = R are

u r e - u r i = 1 R 2 B - RA = D s ,
(20)
σ rr e - σ rr i = - 4 μ R 3 B - ( 3 λ + 2 μ ) A = 0 ,
(21)

where u r i and σ rr i are the radial displacement and radial stress in region 1, respectively. Then, Müller (2001) obtained

d s = 4 μ D s 3 ( λ + 2 μ ) ,
(22)
Δ s = ( λ + 2 μ / 3 ) D s ( λ + 2 μ ) .
(23)

Using volume ΔVs = 4π R2Δs, the radial displacement is written as

u r e = Δ V s 4 π r 2 ,
(24)

where

Δ V s = ( λ + 2 μ / 3 ) ( λ + 2 μ ) Δ V D .
(25)

Here, ΔVD = 4π R2Ds.

Isotropic source

We next consider the displacement field due to an isotropic source represented by three vector dipoles (S3). The static displacement field due to the moment tensor M p q is given as

u i = 1 8 πμ r 2 λ + μ λ + 2 μ 3 γ i γ p γ q - γ i δ pq + γ p δ iq + γ q δ ip + 2 γ q δ ip M pq ,
(26)

where γ i is the directional cosine and δ i j is the Kronecker delta (Aki and richards 2002). Using the moment tensor for the isotropic source in Equation 7, we obtain the radial displacement due to the isotropic source in spherical coordinates as follows:

u r = M 0 Iso 4 π ( λ + 2 μ ) r 2 .
(27)

For an isotropic source represented by three tensile cracks (S4), Equations 8 and 27 yield

u r = ( λ + 2 μ / 3 ) 4 π ( λ + 2 μ ) Δ V i r 2 .
(28)

Following Okada (1992), we find that the radial displacement on the free surface at the distance r from the isotropic source in a half-space is given as

u r = ( λ + 2 μ / 3 ) 2 π ( λ + μ ) Δ V i r 2 .
(29)

Comparing Equation 28 with Equation 29, we find, as in the spherical source or Mogi model (Equations ?? and 17), that they differ by a factor of 2(λ+2μ)/(λ+μ) (Figure 4b).

Eshelby spherical source

We finally consider the displacement field due to an Eshelby spherical source (S5). Before the release of the applied traction, the radial stress in region 1 is given as

σ rr i =-(λ+2μ/3) 3 Δa R .
(30)

In the equilibrium state after the traction release, the displacements and stresses at r=R are given as

u r i = u r e = Δ s ,
(31)
σ rr i = σ rr e = - ΔP.
(32)

Using these boundary conditions, we find that

u r e = Δ V s 4 π r 2 ,
(33)

where

Δ V s = ( λ + 2 μ / 3 ) ( λ + 2 μ ) Δ V f .
(34)

Here, ΔVf = 4π R2Δa. Comparing Equation 34 with Equation 25, we see that Ds = Δa.

Instead of applying the stress-free strain, we may apply a stress or pressure (δ P) to the removed source volume without deforming it (strain-free stress) (S6, Figure 2b). We put the stressed volume back in the hole, weld the material across the cut, and release the stress, which results in the volume source expanding to R + Δs. Before the release of the stress, the radial stress in region 1 is given as

σ rr i =-δP.
(35)

The boundary conditions after the stress release are the same as those for the stress-free strain (Equations 31 and 32), and so we find that the displacement field due to this source is

u r e = Δ V s 4 π r 2 ,
(36)

where

Δ V s = 4 π R 3 3 ( λ + 2 μ ) δP.
(37)

We note that δ P is the applied pressure, which may correspond to the stress glut (Backus and Mulcahy 1976), and ΔP is the pressure after the equilibrium state.

Sources in a bimaterial medium

We consider a bimaterial medium, in which Lamé’s constants are λ and μ in region 1 and λ and μ in region 2. This medium may be viewed as a model of a magma or hydrothermal reservoir, and its source representations are critically important to understand volcanic processes.

For a spherical source (S1), no assumption is made on region 1, and the boundary conditions at r = R in a bimaterial medium are given as those in a homogeneous medium (Equations 12 and 13). Therefore, the displacement field is the same as in a homogeneous medium given in Equation 15.

For an isotropic source, three vector dipoles (S3) can be represented as three tensile cracks (S4). For the moment tensor of each tensile crack, no assumption is made for the material inside a tensile crack, and we may not able to define an isotropic source in a bimaterial medium.

For a spherical crack source (S2), the boundary conditions at r = R in a bimaterial medium may be given as

u r e - u r i = 1 R 2 B - RA = D s ,
(38)
σ rr e - σ rr i = - 4 μ R 3 B - ( 3 λ + 2 μ ) A = 0 ,
(39)

where D s = d s + Δ s . Here, d s is the inward crack wall displacement in region 1. We then obtain

u r e = Δ V s 4 π r 2 ,
(40)

and

Δ V s = ( λ + 2 μ / 3 ) ( λ + 2 μ / 3 + 4 μ / 3 ) Δ V D ,
(41)

where V D =4π R 2 D s .

For an Eshelby spherical source with stress-free strain (S5), the source volume of region 1 expands to R + Δa, and the radial stress in region 1 before the release of the applied traction may be given as

σ rr i =-( λ +2 μ /3) 3 Δ a R .
(42)

The boundary conditions at r = R are the same as those in a homogeneous medium (Equations 31 and 32), and so we have

u r e = Δ V s 4 π r 2 ,
(43)

and

Δ V s = ( λ + 2 μ / 3 ) ( λ + 2 μ / 3 + 4 μ / 3 ) Δ V f ,
(44)

where Δ V f =4π R 2 Δ a .

For an Eshelby spherical source with strain-free stress δ P (S6), we obtain

u r e = Δ V s 4 π r 2 ,
(45)

and

Δ V s = 4 π R 3 3 ( λ + 2 μ / 3 + 4 μ / 3 ) δ P .
(46)

Results and discussion

We summarize the radial displacement fields due to the sources examined above in Figure 5. Except for an isotropic source (Figure 5c,d), the displacement fields are given as ΔVs/(4π r2). Let us first consider the relationship between a spherical source (S1, Figure 5a) and an isotropic source represented by three vector dipoles (S3, Figure 5c). The displacement fields of the spherical source and isotropic source are given as

Figure 5
figure 5

Summary of displacement fields and source parameters in homogeneous and bimaterial media for various sources. (a) A spherical source (S1), (b) a spherical crack source (S2), (c) an isotropic source (three mutually perpendicular vector dipoles) (S3), (d) an isotropic source (three mutually perpendicular tensile cracks) (S4), (e) an Eshelby spherical source with stress-free strain (S5), and (f) an Eshelby spherical source with strain-free stress (S6).

u r e = Δ V s 4 π r 2 ,
(47)
u r = M 0 Iso 4 π ( λ + 2 μ ) r 2 ,
(48)

respectively. The equivalence of these displacements at r > R provides the following relation:

Δ V s = M 0 Iso ( λ + 2 μ ) .
(49)

For an isotropic source represented by three tensile cracks (S4), we obtain Equation 8, i.e.,

M 0 Iso =(λ+2μ/3)Δ V i .
(50)

The relation of Equation 49 is also held between a spherical crack source (S2) and an isotropic source (S3) and between an Eshelby spherical source (S5 and S6) and an isotropic source (S3) (see Figure 5). Then, using Equations 25 and 34, we find that the following relation is held for a spherical crack source (S2) and an Eshelby spherical source with stress-free strain (S5):

M 0 Iso =(λ+2μ/3)Δ V f ,
(51)

where ΔVf = 4π R2Δa = 4π R2Ds. We see from Equations 50 and 51 that ΔVi = ΔVf: the sum of three tensile crack volumes is equal to a spherical crack volume or stress-free volume. Therefore, the two seismic moment equations (Equations 1 and 2) can be understood as follows: the first equation from the equivalence of the displacement fields due to a spherical source (S1) and an isotropic source (S3) (Equation 49) and the second equation from the equivalence of three tensile crack volumes and a spherical crack volume or stress-free volume (Equations 50 and 51). The first equation can be used to estimate the actual volume from the seismic moment of an isotropic source (S3).

As shown in Subsection 'Sources in a bimaterial medium’, no assumption is made on region 1 for a spherical source (S1), and its source parameter is given as ΔVs in both homogeneous and bimaterial media (Figure 5a). Thus, we may use Equation 49 with M 0 Iso in a homogeneous medium to estimate ΔVs, regardless of the assumption on the medium.

It is interesting to note that we may be able to estimate a spherical crack volume or stress-free volume (Equation 51) that cannot be observed. This estimation, however, is based on our interpretation of the source process in region 1. For example, if we consider an Eshelby spherical source with strain-free stress (S6, Figure 5f), we obtain

M 0 Iso = 4 π R 3 3 δP,
(52)

which gives an estimation of δ P in region 1. If we assume a bimaterial medium for a spherical crack source (S2) and consider an Eshelby spherical source with stress-free strain (S5), we have

M 0 Iso = ( λ + 2 μ ) ( λ + 2 μ / 3 ) ( λ + 2 μ / 3 + 4 μ / 3 ) Δ V f ,
(53)

where Δ V f =4π R 2 D s =4π R 2 Δ a . For an Eshelby spherical source with strain-free stress (S6) in a bimaterial medium, we obtain

M 0 Iso = 4 π ( λ + 2 μ ) R 3 3 ( λ + 2 μ / 3 + 4 μ / 3 ) δ P .
(54)

These results suggest that the source estimation in region 1 depends on our interpretation, whereas we can uniquely estimate the source parameter of ΔVs in region 2 from M 0 Iso .

In volcano seismological studies, waveform inversions have been performed to determine moment tensor solutions. If one interprets the isotropic source determined by waveform inversion to be spherical, then Equation 49 can be used to estimate the actual volume ΔVs. In volcano geodetic studies, the Mogi model has been used to interpret volcano deformations. Because the Mogi model is based on a spherical source, ΔVs is estimated directly. On the other hand, Okada (1992) model uses an isotropic source. Thus, if one estimates the seismic moment of an isotropic source by using Okada’s model, the relation given in Equation 49 should be used to estimate the volume change due to the corresponding spherical source or Mogi model.

We see that the displacement fields due to the described sources in an infinite medium are all expressed as a function of r-2 (Figure 5). The divergence of such a function is zero, which means that a source volume change must be conserved everywhere in an infinite medium. For all the sources except for an isotropic source, we obtain the following relation in an infinite medium:

4π r 2 u r =Δ V s .
(55)

The left-hand side represents the volume change of a sphere with an arbitrary radius r centered at the source, which is equal to ΔVs. The source volume change in this displacement field is, thus, conserved in an arbitrary subspace enclosing the source. For an isotropic source represented by three tensile cracks (S4), we obtain the following relation:

4 π r 2 u r = ( λ + 2 μ / 3 ) ( λ + 2 μ ) Δ V i .
(56)

This relation indicates that the source volume change is not conserved. The non-conserved volume may be attributed to compressed volumes in regions close to the cracks (Figure 6). The non-conserved volume is equal to the spherical crack volume or stress-free volume, although the reason why these volumes are equal is not clear. We note that the conservation of a source volume does not hold for displacement fields in a half-space because of the effect of the free surface.

Figure 6
figure 6

Force systems for three mutually perpendicular tensile cracks.

Let us further consider a spherical crack source in a homogeneous medium with the following boundary conditions at r = R:

u r e = Δ s ,
(57)
σ rr e = - ΔP ,
(58)
u r i = - d s ,
(59)
σ rr i = Δ P ̄ ,
(60)
u r e - u r i = 1 R 2 B - RA = Δ s + d s ,
(61)
σ rr e - σ rr i = - 4 μ R 3 B - ( 3 λ + 2 μ ) A = - ( ΔP + Δ P ̄ ) .
(62)

These boundary conditions represent the displacement discontinuity Δs + ds and the stress or traction discontinuity -(ΔP+Δ P ̄ ) along the spherical surface. Using Equations 61 and 62, we obtain

u r e = R 3 ( λ + 2 μ ) ( ΔP + Δ P ̄ ) + ( λ + 2 μ / 3 ) ( λ + 2 μ ) ( Δ s + d s ) R r 2 .
(63)

This equation represents the displacement field due to a generalized spherical crack source having both displacement and traction discontinuities. Using Equations 59 and 60, we obtain

Δ P ̄ =-(3λ+2μ) d s R .
(64)

Substituting Equation 64 into Equation 63, we obtain

u r e = R 3 ( λ + 2 μ ) ΔP + ( λ + 2 μ / 3 ) ( λ + 2 μ ) Δ s R r 2 .
(65)

From Equations 57 and 58, we see that Equation 65 becomes Equation 15 for the displacement field due to a spherical source. Therefore, Equation 65 may be regarded as the general form of the equation for the displacement field due to a spherical source represented by ΔP and Δs.

Conclusions

We systematically examined the following seismic sources and their static displacement fields: a spherical source (radius R) with outward radial expansion (R + Δs) (S1), a spherical crack source with outward and inward crack wall motions along the spherical surface with displacement discontinuity Ds (S2), an isotropic source represented by three mutually perpendicular vector dipoles (S3) or three tensile cracks (S4), and an Eshelby spherical source undergoing a transformational expansion (R + Δa) (S5) or strain-free stress (δ P) (S6). Our examinations show that the static radial displacement fields due to these sources except for an isotropic source are given as ΔVs/(4π r2), where ΔVs = 4π R2Δs is the actual volume. The equivalence of these fields and that due to an isotropic source (S3) gives Δ V s = M 0 Iso /(λ+2μ), where M 0 Iso is the seismic moment of an isotropic source. The actual volume ΔVs may be uniquely determined from M 0 Iso in homogeneous and bimaterial media. We see that the source volume change for an isotropic source defined by three tensile cracks (S4) is equal to the spherical crack volume (4π R2Ds) and stress-free volume (4π R2Δa). We then obtain M 0 Iso =(λ+2μ/3)Δ V f , where ΔVf is the spherical crack or stress-free volume. This relation depends on our interpretation of the source process. If we consider an Eshelby spherical source with strain-free stress (S6), M 0 Iso is related to δ P. We obtain a different estimation of ΔVf or δ P in a bimaterial medium. This study provides a unified view for quantifying a spherical source using the seismic moment of an isotropic source determined from waveform inversion analysis in volcano seismological studies.

References

  • Aki K, Richards PG: Quantitative seismology, 2nd edn. Sausalito: University Science Books; 2002.

    Google Scholar 

  • Ampuero JP, FA Dahlen: Ambiguity of the moment tensor. Bull Seism Soc Am 95 2005, 390–400.

    Google Scholar 

  • Backus G, Mulcahy M: Moment tensors and other phenomenological descriptions of seismic sources–I. Continuous displacements. Geophys J R Astr Soc 1976, 46: 341–361. 10.1111/j.1365-246X.1976.tb04162.x

    Article  Google Scholar 

  • Chouet B: A seismic model for the source of long-period events and harmonic tremor. In Volcanic seismology. Edited by: Gasparini P, Scarpa R, Aki K. Berlin: Springer; 1996:133–156.

    Google Scholar 

  • Chouet B: Introduction to quantitative volcano seismology: fluid-driven sources. In Modeling volcanic processes the physics and mathematics of volcanism. Edited by: RMC, Lopes, TKP, Gregg, Fragents SA, Gregg , TKP , Lopes , RMC . Cambridge: Cambridge University Press; 2013:331–358.

    Chapter  Google Scholar 

  • Dreger D, Johnston M, Tkalčić H: Dilational processes accompanying earthquakes in the Long Valley Caldera. Science 2000, 288: 122–125. 10.1126/science.288.5463.122

    Article  Google Scholar 

  • Eshelby JD: The determination of the elastic field of an ellipsoidal inclusion, and related problems. Proc Roy Soc Ser A 1957, 241: 376–396. 10.1098/rspa.1957.0133

    Article  Google Scholar 

  • Fukuyama E, Takeo M: Analysis of the near-field seismogram observed during the eruption of Izu-Oshima Volcano on November 16. Bull Volcanol Soc Jpn 1987, 35: 283–297.

    Google Scholar 

  • Julian BR, Miller AD, Foulger GR: Non-double-couple earthquakes 1. Theory. Rev Geophys 1998, 36: 525–549. 10.1029/98RG00716

    Article  Google Scholar 

  • Kanamori H, Given J, Lay T: Analysis of seismic body waves excited by the Mount St. Helens eruption of May 18,. J Geophys Res 1980, 89: 1856–1866.

    Article  Google Scholar 

  • Kawakatsu H, Yamamoto M: Volcano seismology. In Kamanori H (ed) Treatise on geophysics volume 4—earthquake seismology. Amsterdam: Elsevier; 2007:389–420.

    Google Scholar 

  • Kumagai H: Source quantification of volcano seismic signals. In Encyclopedia of complexity and systems science. Edited by: Mayers RA, Mayers RA . New York: Springer; 2009.

    Google Scholar 

  • Kumagai H, Placios P, Ruiz M, Yepes H, Kozono T: Ascending seismic source during an explosive eruption at Tungurahua volcano. Ecuador Geophys Res Lett 2011., 38: L01306 doi:10.1029/2010GL045944 L01306 doi:10.1029/2010GL045944

    Google Scholar 

  • Legrand D, Kaneshima S, Kawakatsu H: Moment tensor analysis of near-field broadband waveforms observed at Aso volcano, Japan. J Volcanol Geotherm Res 2000, 101: 155–169. 10.1016/S0377-0273(00)00167-0

    Article  Google Scholar 

  • Maeda M, Takeo M: Take very-long-period pulses at Asama volcano, central Japan, inferred from dense seismic observations. Geophys J Int 2011, 185: 265–282. 10.1111/j.1365-246X.2011.04938.x

    Article  Google Scholar 

  • Mogi K: Relation between the eruptions of various volcanoes and the deformations of the ground surfaces around them. Bull Earthq Res Inst Univ Tokyo 1958, 36: 99–134.

    Google Scholar 

  • Müller G: Volume change of seismic. Bull Seism Soc Am 2001, 91: 880–884. 10.1785/0120000261

    Article  Google Scholar 

  • Nishimura T, Iguchi M: Volcanic earthquakes and tremor in Japan. Kyoto University Press, Kyoto 2011.

    Google Scholar 

  • Ohminato T, Takeo M, Kumagai H, Yamashina T, Oikawa J, Koyama E, Tsuji H, Urabe T: Vulcanian eruptions with dominant single force components observed during the Asama 2004. volcanic activity in Japan Earth Planets Space 2006, 58: 583–593.

    Article  Google Scholar 

  • Okada Y: Internal deformation due to shear and tensile faults in a half-space. Bull Seism Soc Am 1992, 82: 1018–1040.

    Google Scholar 

  • Richards PG, Kim WY: Equivalent volume sources for explosions at depth: theory and observations. 2005, 95: 401–407.

    Google Scholar 

  • Sezawa K: The plastico-elastic deformation of a semi-infinite solid body due to an internal force. 1931, 9: 398–406.

    Google Scholar 

  • Tameguri T, Iguchi M, Ishihara K: Mechanism of explosive eruptions from moment tensor analyses of explosion earthquakes at Sakurajima volcano, Japan. Bull Volcanol Soc Jpn 2002, 47: 197–216.

    Google Scholar 

  • Uhira K, Takeo M: The source of explosive eruptions of Sakurajima volcano, Japan. J Geophys Res 1994, 99: 17775–17789. 10.1029/94JB00990

    Article  Google Scholar 

  • Wielandt E: On the relationship between seismic moment and source volume. 2003.http://www.geophys.uni-stuttgart.de/~erhard/skripte/ew/isomoment/ Accessed May 30 2003

    Google Scholar 

Download references

Acknowledgements

We thank Mare Yamamoto, Minoru Takeo, and Shiro Hirano for their useful discussions. This work was supported by the Japan Science and Technology Agency (JST).

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to Hiroyuki Kumagai.

Additional information

Competing interests

The authors declare that they have no competing interests.

Authors’ contributions

HK developed the main framework of the study and drafted the manuscript. All authors carried out the theoretical considerations and read and approved the final manuscript.

Authors’ original submitted files for images

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 2.0 International License (https://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Reprints and permissions

About this article

Cite this article

Kumagai, H., Maeda, Y., Ichihara, M. et al. Seismic moment and volume change of a spherical source. Earth Planet Sp 66, 7 (2014). https://doi.org/10.1186/1880-5981-66-7

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/1880-5981-66-7

Keywords