Category Archives: Applied Mathematics

Seismic wave attenuation: geometrical spreading, anelasticity, multipathing and scattering

As we discussed in class, seismic waves can lose energy through reflection, geometrical spreading and intrinsic attenuation, also referred as anelasticity.

Geometrical spreading depends on the distance r the wave has propagated from the source. In a uniform material, seismic waves propagate away from their source as spherical wave front of increasing area. Because of the conservation of energy, the energy per unit area of wave front decreases as the distance from the source increases. For surface waves, in the case of a homogeneous flat earth, the energy per unit area of wave front decreases as 1/r and hence the amplitude, which is proportional to the root square of the energy, decreases as 1/√r. For body waves, the energy per unit area of wave front decreases as 1/r2 and hence the amplitude decreases as 1/r.

On the other hand, anelasticity reduces seismic wave amplitudes by converting part of their kinetic energy to frictional heat by permanent deformation of the medium. Anelasticity is characterized by the frequency-dependent quality factor Q, which is a measure of the the energy lost per oscillation of the seismic wave : Q = 2πE/∆E . So the smaller Q, the larger the energy loss.The loss of energy will lead to exponential decay of the seismic wave amplitude : A(t) = A0*e-πft/Q . The smaller Q and the larger the frequency (i.e. more oscillations per second), the larger the attenuation and the seismic wave amplitude decay.

Two other processes can also reduce seismic wave amplitudes: multipathing and scattering. Multipathing and scattering can be thought of as elastic processes. They conserve energy and decrease or increase the amplitude of an incoming wave by shifting its energy to an earlier or later arrival.

Seismic wave multipathing is caused by velocity variations within the medium of propagation. According to Fermat’s principle, seismic waves follow the least-time path of propagation between two points in a medium. Lateral velocity variations in the medium will then cause seismic waves to focus in high velocity regions and defocus in low velocity regions.The spacing between seismic rays in a region represents the energy density in this region. The further apart the rays are, the lower the amplitudes of the recorded wave. By contrast, the closer the rays are, the larger the wave amplitudes. So the seismic waves arriving at a station have usually followed different ray paths in addition to the ideal, direct path and the region of the earth they sampled forms a volume called Fresnel zone. Multipathing can be a significant attenuation effect because most seismic activity occur at plate boundaries and velocity heterogeneities are important in these regions.

Likewise, heterogeneities within the propagation medium cause a propagating wave field to be scattered. These heterogeneities can be velocity anomalies but also material heterogeneities such as mineral boundaries, pore edges, cracks… Scattering will cause part of the energy released by an earthquake to arrive later at a receiver (i.e. after the initial pulse) as a coda (i.e. tail of incoherent energy that decays over a few seconds to a few minutes). Whether a seismic wave will be scattered or not when encountering a heterogeneity depends on the ratio of the heterogeneity size to the wavelength and the propagation distance in the heterogeneous medium. If the heterogeneity is large compared to the wavelength, the seismic energy will follow a different ray path (i.e. multipathing effect). However if the heterogeneity and the wavelength have the same order of magnitude, the seismic energy will be scattered. Heterogeneities much smaller than the wavelength will just change the medium’s “bulk” properties. Scattering can be significant in the continental crust because of the presence of many small-scale geologic structures that can significantly affect short wavelength waves (i.e. tens of kilometers or smaller).

Reference : Stein, S., & Wysession, M. (2009). An introduction to seismology, earthquakes, and earth structure. John Wiley & Sons.

Frontier beneath our feet: Seismic study aims to map Earth’s interior in 3-D

Frontier beneath our feet: Seismic study aims to map Earth’s interior in 3-D

Seismic waves carry information about the Earth’s structure. Thus seismologists combine seismology and computer science to map the Earth’s interior. Here, the Princeton University attempts to map the deep structures on 3D-map. The project will use M>5 worldwide earthquakes recorded on thousands of seismic stations through NSF and research institutions for seismology.

Scientist from Princeton University are interested on map the mantle up welling and plumes, so it will be great make some cross correlations with the Africa rift system currently studied by Andy Nyblade’s group.

http://www.princeton.edu/main/news/archive/S42/59/33Q27/index.xml?section

Using Cross-Correlation to Detect Small Earthquakes: Part I

This post will be the first in a small series. The purpose of the small series is to describe aspects of my Master’s research project (Detecting Small Earthquakes on Remote Oceanic Transform Faults) in several steps. In this post I introduce our main analysis approach, cross-correlation, and briefly describe our analysis process.

I’m hoping these posts will help me with writing my thesis, so I would greatly appreciate any suggestions to help improve the material.

You can access the content for Part I here. Please let me know if you have trouble viewing the file.

Thanks!
Rhiannon

Moment-Tensor Decomposition

While the moment tensor makes relating seismograms to seismic sources more convenient, it adds a non-unique interpretational step once you have the Mij values. You have to decompose the moment tensor into double-couples, CLVD sources, isotropic sources, etc.  Here is a link to a script I wrote for Matlab to do the decomposition.

mtensor_decomp.m

“All models are wrong, but some are useful”

Scientists use observations as input to construct models to explain the nature, but a model’s value depends on the quality of data and the assumptions behind it. Below is a link to an opinion piece by Edward H. Field of the USGS about models related to earthquakes forecasts for California. The author notes that some models account for segmented fault ruptures while others exclude multi-fault ruptures and that although both use different assumptions, both are useful. The author concludes that we have to quantify relatively which model will be useful under which conditions.

http://srl.geoscienceworld.org/content/86/2A/291.full.pdf+html?sid=df3d027d-8c38-45b6-bf81-f49685385f89

Understanding Einstein Summation Notation

In class we are using a continuum mechanics approach to study stresses and strains within the earth. Many of the equations we will use to describe the motion of the earth are written using Einstein summation notation. Therefore, it is important to understand how this convention works.

When I first encountered the Einstein summation notation I had some difficulty understanding what the equations actually meant. I’ve included some of the external resources that helped me to better understand this convention and thought they may be helpful to others who are less familiar with this material.

The MOOC Alternative (Youtube Videos):

This is a series of short youtube videos related to the Einstein summation convention. The first video starts with basic concepts of scalars, vectors, and tensors. With those building blocks, the subsequent videos go on to discuss summation over indices, dummy indices, and the Kronecker delta (all of which are important to the derivations we do in class).

Intro to indicial notation by Theo Hopman:

For a basic introduction, I found the first two sections of this document to be most helpful. The first section discusses free indices, summation over an index, and dummy indices. The second part discusses the Kronecker delta and Levi-Civita functions. I found the other sections interesting to read although I am not sure how much of that material will be necessary for this course. Nevertheless, I would recommend having a quick look over them particularly the examples in section 4.