Category Archives: Seismic Wave Propagation

What’s all this talk about earthquakes? Part II

Moonquakes – a review of lunar seismicity

Beginning in 1969, a mere 66 years after the first ever engine-powered flight, the crew of the Apollo 12 mission deployed a set of seismic instruments on the surface of the Moon. Data was collected and radioed to Earth for 8 years until the machines were switched off in 1977.

A11setup5-2

Buzz Aldrin deploys a seismometer in the Sea of Tranquillity. Image taken from science.nasa.gov

The Moon is as close of a partner to Earth as it is a satellite, but there is still one large difference between the two bodies, active tectonics. The moon is considered to be tectonically dead; no new crust has been made since mare basalt formations roughly 3 billion years ago. But that doesn’t mean that we didn’t see seismicity!

Seismic instruments on the Moon recorded four fundamental sources of seismic energy.

  • deep moonquakes ~700km below the surface
  • impact vibrations
  • thermal quakes caused by the expansion of the crust due to exposure to light after a two weeklong, deep freeze night
  • shallow 20-40km depth seismicity

Seismic energy release on the moon is on the order of a billion times less than on Earth and is almost all a result of tidal forces from the earth.

quakes

Moon cartoon showing depth distribution of moonquakes. Image taken from the-moon.wikispaces.com

All of this activity helped lunar seismologists to explore the interior of the planet. They found a solid inner core, fluid outer core, a layer of partial melt atop the fluid outer core and of course a large mantle. They also found that the moon as a whole is much more rigid than Earth. The largest of the 28 shallow moonquakes recorded between 1972 and 1977 had a magnitude of 5.5 and energy was recorded for over 10 minutes (the Moon exhibits very little attenuation).

First person to guess why the moon is so much more rigid than the earth gets a high five!

Origin of the Earth’s hum

Earthquakes can cause the Earth to vibrate over period of days to months. However even in the absence of earthquakes, the Earth keeps vibrating at very low frequencies. This continuous vibration generated by very slow seismic waves with periods greater than 50 s was first discovered in the late 1990s and since then several explanations involving ocean wave propagation have been proposed to explain the origin of this phenomenon. One of the theories suggests that this continuous hum is generated by the constructive interference of microseismic waves that are created during the collision of ocean waves moving in opposite directions. However the period of these microseismic waves is lower than 13 s. So the origin of this Earth’s hum remained unexplained until recently. In February 2015, a group of french researchers proposed the interaction of long ocean waves with the seafloor as a possible explanation and by integrating this new theory into their models, were able to generate microseismic signals with periods ranging from 13 to 300 s. They concluded that both the collision of opposing ocean waves and in a larger extent the movement of long ocean waves over the ocean bottom are responsible for the Earth’s hum.

To learn more about this study, please see the following paper: Ardhuin, F., Gualtieri, L., & Stutzmann, E. (2015). How ocean waves rock the Earth: two mechanisms explain microseisms with periods 3 to 300 s. Geophysical Research Letters.

Reference : http://blogs.agu.org/geospace/2015/04/07/new-study-explains-source-of-earths-mysterious-ringing/

Counteracting the effects of an earthquake with a seismic metamaterial

In 2013, a team of scientists at a french construction firm published a paper describing a new way of counteracting the destructive effects of an earthquake using a seismic metamaterial. In this paper, they explain how they tried to attenuate the amplitude of seismic waves at the free surface by modifying the energy distribution. Their experiment consisted in simulating an earthquake and using a metamaterial made of a grid of vertical and empty cylindrical columns bored into soil near the earthquake source to attenuate the energy released by the surface waves. They were able to reflect the energy of the incoming surface waves and hence significantly dampen their energy.

To learn more about the details and the results of this experiment, please see their paper:http://arxiv.org/pdf/1301.7642v1.pdf

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

Pennsylvania P-Wave Travel Times

For this blog post I figured I would use some of my prior work and construct a P wave travel time curve.  These travel times are from small seismic events recorded across Pennsylvania. The travel time versus distance observations are shown below.

Observed P-wave travel times for small earthquakes across Pennsylvania (Kyle Homann, 2015).

Observed P-wave travel times for small earthquakes across Pennsylvania (Kyle Homann, 2015).

This curve is about what I expected, with some picks that could be questionable, but the majority are consistent.  I think that around 150 km distance, some S wave picks managed to sneak into data.  The slope is slightly less at greater distances along with a “bulge” around 200 km.  This may be a result of Pn overtaking Pg.  If you recall, Pg is the P wave that travels within the crust and Pn is the refracted wave from the crust-mantle boundary. As expected, the Pg wave speed indicated by the slope of the graph is ~6 km/s.

Books & Videos – Tsunami, Observational Seismology, & Programming

Penn State has good libraries and learning resources that you should take advantage of and immerse yourself in your studies. Here are links to two books available from PSU that you may find interesting (or that you should, since you are in this class…). Browse them and read what you can.

Tsunami, The Underrated Hazard (Second Edition) by Edward Bryant

Routine Data Processing in Earthquake Seismology With Sample Data, Exercises and Software by Jens Havskov and Lars Ottemoller

And if you are new to programming, but interested enough to invest some time to learn, try this video from http://psu.lynda.com, which requires about 5 hrs of investment, but my help clarify some of the fundamentals for you. I don’t think you need to learn to program in javascript, but the ideas of loops, conditionals, etc are common to all languages.