5:15 PM - 6:30 PM
[SSS32-P11] Theoretical calculation of internal stress/strain changes caused by earthquakes: the effectiveness of the reciprocity theorem in a spherical earth
Keywords:co-seismic deformation, internal deformation, spherically symmetric earth, reciprocity theorem
After Takeuchi & Saito (1972), the equations governing co-seismic deformation field result in first-order inhomogeneous differential equations by expanding displacement and stress by the spherical harmonics. To solve these differential equations, in the conventional method, (i) we obtain the complementary solutions, (ii) find the particular solution, and (iii) add them so that the final solution satisfies the surface boundary condition. The magnitude of the final solution is the order of (rs/rp)n, where the rs and rp (>rs) are the radii where the source is located and the deformation is evaluated, respectively, and n is the degree of the spherical harmonics. On the other hand, the magnitude of the solutions obtained by the process (i) and (ii) is the order of (rp/rs)n. This means that, in process (iii), LSD cannot be avoided at a large degree n because we need to add the numbers whose magnitude is (rp/rs)n/(rs/rp)n=(rp/rs)2n times larger than that of the final solution. For example, the ratio (rp/rs)2n becomes 1012 at n=8,000 when the deformation at a depth of 10 km (rp=6361 km) due to a source at a depth of 20 km (rs=6351 km) is considered. We confirmed that LSD occurs around n=8,000 in actual computation.
In the method using the reciprocity theorem, we obtain (a) the solutions x1 at the radius rs caused by external sources such as tide and (b) a solution x2 that has a unit jump at the radius rp. They are easily calculated by numerical integration. (c) Finally, the final solution is obtained by multiplying x1 and x2 together. This means that we can avoid the LSD that occurs in doing addition.