Receiver function

Last updated

The receiver function technique is a way to image the structure of the Earth and its internal boundaries by using the information from teleseismic earthquakes recorded at a three-component seismograph.

Contents

A teleseismic P-wave will generate P-to-S conversions at boundaries, such as the Moho (crust-mantle boundary), beneath the seismograph. The difference in travel time between the generated S-wave and P-wave contains information about the depth to the boundary and about the P- and S-wave velocities. If further reverberations are included, more detailed structure can be resolved. [1] [2] This is done by deconvolution of the incoming vertical and longitudinal components of the seismogram, which removes the common part of the components - namely, the source and travel path information. [3] The resulting waveform is the receiver function.

Similarly, a teleseismic S-wave will generate an S-to-P conversion beneath the seismic station.

Method

Incident P and S wave phases for a receiver function. Receiver Function Incident Wave.svg
Incident P and S wave phases for a receiver function.

As a P wave in the mantle passes upwards through the Moho, it is partially converted into an S wave. Both the P-wave and S-wave (known as Ps) are picked up by the seismometer on the Earth's surface and can be used to analyze discontinuities within the Earth. In addition to these P and Ps waves, additional phases are created by multiple reflections. These phases include: PpPmp PpSmp, PpPms, and PpSms, as well as PsPmp, PsSmp, PsPms and PsSms. (See seismic phase notation for more info). The phases for which the last "leg" is a P-wave (P, Ps, Smp, etc.) are dominantly recorded on the vertical component of the seismograph, whereas the phases for which the last "leg" is an S-wave (Ps, Pp, Sms, etc.) are dominantly recorded on the horizontal component. [4]

The primary method for creating a receiver function is based on analyzing the product of waves that pass from the mantle through the Moho boundary. The large compositional differences between the crust and the mantle cause large differences in seismic waves as they pass through the discontinuity. [5] Receiver functions use Snell's law refraction of P waves and converted S waves to estimate the depth of the Moho. They are only generated clearly if the central angle (angle subtended at the center of the Earth) between the seismic event and the seismograph station is between 30 and 95 degrees (between about 3300 and 8900 km separation). [4] The method is also most effective when the seismic event causing the waves occurs significantly below the Moho, which is important to avoid surface interference. [4]

Over time, several seismic events can occur in the same geographic area, with each event causing a corresponding vertical waveform and horizontal waveform. Once several observations have been collected, the waveforms can be summed together for the vertical component and for the horizontal component. The summed seismogram reduces random noise and makes it easier to see a pattern in the data. By visual inspection, or more commonly by deconvolution of the two waveforms, it is possible to identify each of the relevant phases of the P wave conversions. With the timing of the phases, it is then possible to model the seismic velocities within the crust and the Moho depth beneath the seismic station. [6]

Stacking

Where there are many adjacent seismograph stations, it is possible to "stack" receiver function data across seismograph stations to build a 2D or even 3D model of the depth of the Moho. [7] This is possible because each station can determine the depth of the Moho at its own location (essentially a 1D measurement). Data from multiple individual data points from adjacent stations can be grouped together and plotted side by side to create a unified graph of the Moho depth over a given area. [7]

For deeper interfaces in the mantle, the moveout of Ps converted phases can be corrected for by delaying the time window of the horizontal-component seismograms by a predicted delay time. [8] This delay will align the phases of any pulses on the horizontal components that suffer the predicted time delays.

Applications

Receiver functions contain detailed information on the average seismic velocities within the crust and on the depth of the Moho at a specific location. This data alone can be useful in obtaining information about a specific location. [6] But when receiver function data from one seismic station is combined with data from many other stations, it is possible to build a detailed map of the Moho depth and of seismic velocity across a large geographic area. The dipping top surfaces of subducting lithosphere is often sufficiently sharp to generate P-to-S converted phases that can be detected at depths up to 100 km or more.

This data can be used for a variety of purposes. It can be used to note variations in the depth of the crust. Receiver functions have been used, for example, to discover depressions in the Moho below mountains in southwest Japan. [4] This data can also be used to better understand earthquakes that cause natural disasters. [4] Additionally, maps of seismic velocities and crustal thickness are useful as baseline data for additional seismological studies. [5]

Data from receiver functions can also be used in conjunction with data, such as data from controlled source seismology, to provide higher resolution 3D maps of the Earth's crust. [5]

Related Research Articles

<span class="mw-page-title-main">Seismology</span> Scientific study of earthquakes and propagation of elastic waves through a planet

Seismology is the scientific study of earthquakes and the propagation of elastic waves through the Earth or other planetary bodies. It also includes studies of earthquake environmental effects such as tsunamis as well as diverse seismic sources such as volcanic, tectonic, glacial, fluvial, oceanic, atmospheric, and artificial processes such as explosions. A related field that uses geology to infer information regarding past earthquakes is paleoseismology. A recording of Earth motion as a function of time is called a seismogram. A seismologist is a scientist who does research in seismology.

<span class="mw-page-title-main">Seismic wave</span> Seismic, volcanic, or explosive energy that travels through Earths layers

A seismic wave is a wave of acoustic energy that travels through the Earth or another planetary body. It can result from an earthquake, volcanic eruption, magma movement, a large landslide, and a large man-made explosion that produces low-frequency acoustic energy. Seismic waves are studied by seismologists, who record the waves using seismometers, hydrophones, or accelerometers. Seismic waves are distinguished from seismic noise, which is persistent low-amplitude vibration arising from a variety of natural and anthropogenic sources.

<span class="mw-page-title-main">Seismometer</span> Instrument that records seismic waves by measuring ground motions

A seismometer is an instrument that responds to ground noises and shaking such as caused by quakes, volcanic eruptions, and explosions. They are usually combined with a timing device and a recording device to form a seismograph. The output of such a device—formerly recorded on paper or film, now recorded and processed digitally—is a seismogram. Such data is used to locate and characterize earthquakes, and to study the Earth's internal structure.

<span class="mw-page-title-main">Seismogram</span> Graph output by a seismograph

A seismogram is a graph output by a seismograph. It is a record of the ground motion at a measuring station as a function of time. Seismograms typically record motions in three cartesian axes, with the z axis perpendicular to the Earth's surface and the x- and y- axes parallel to the surface. The energy measured in a seismogram may result from an earthquake or from some other source, such as an explosion. Seismograms can record many things, and record many little waves, called microseisms. These tiny microseisms can be caused by heavy traffic near the seismograph, waves hitting a beach, the wind, and any number of other ordinary things that cause some shaking of the seismograph.

<span class="mw-page-title-main">Mohorovičić discontinuity</span> Boundary between the Earths crust and the mantle

The Mohorovičić discontinuity, usually referred to as the Moho discontinuity, Moho boundary, or just Moho, is the boundary between the Earth's crust and the mantle. It is defined by the distinct change in velocity of seismic waves as they pass through changing densities of rock.

Seismic tomography is a technique for imaging the subsurface of the Earth with seismic waves produced by earthquakes or explosions. P-, S-, and surface waves can be used for tomographic models of different resolutions based on seismic wavelength, wave source distance, and the seismograph array coverage. The data received at seismometers are used to solve an inverse problem, wherein the locations of reflection and refraction of the wave paths are determined. This solution can be used to create 3D images of velocity anomalies which may be interpreted as structural, thermal, or compositional variations. Geoscientists use these images to better understand core, mantle, and plate tectonic processes.

<span class="mw-page-title-main">Andrija Mohorovičić</span> Croatian seismologist and geophysicist

Andrija Mohorovičić was a Croatian geophysicist. He is best known for the eponymous Mohorovičić discontinuity and is considered one of the founders of modern seismology.

<span class="mw-page-title-main">P wave</span> Type of seismic wave

A P wave is one of the two main types of elastic body waves, called seismic waves in seismology. P waves travel faster than other seismic waves and hence are the first signal from an earthquake to arrive at any affected location or at a seismograph. P waves may be transmitted through gases, liquids, or solids.

<span class="mw-page-title-main">Earthscope</span> Earth science program exploring the structure of the North American continent

The EarthScope project was an National Science Foundation (NSF) funded earth science program that, from 2003-2018, used geological and geophysical techniques to explore the structure and evolution of the North American continent and to understand the processes controlling earthquakes and volcanoes. The project had three components: USArray, the Plate Boundary Observatory, and the San Andreas Fault Observatory at Depth. Organizations associated with the project included UNAVCO, the Incorporated Research Institutions for Seismology (IRIS), Stanford University, the United States Geological Survey (USGS) and National Aeronautics and Space Administration (NASA). Several international organizations also contributed to the initiative. EarthScope data are publicly accessible.

Seismic anisotropy is the directional dependence of the velocity of seismic waves in a medium (rock) within the Earth.

<span class="mw-page-title-main">Shadow zone</span> Area not reached by seismic waves from an earthquake

A seismic shadow zone is an area of the Earth's surface where seismographs cannot detect direct P waves and/or S waves from an earthquake. This is due to liquid layers or structures within the Earth's surface. The most recognized shadow zone is due to the core-mantle boundary where P waves are refracted and S waves are stopped at the liquid outer core; however, any liquid boundary or body can create a shadow zone. For example, magma reservoirs with a high enough percent melt can create seismic shadow zones.

Seismic magnitude scales are used to describe the overall strength or "size" of an earthquake. These are distinguished from seismic intensity scales that categorize the intensity or severity of ground shaking (quaking) caused by an earthquake at a given location. Magnitudes are usually determined from measurements of an earthquake's seismic waves as recorded on a seismogram. Magnitude scales vary on what aspect of the seismic waves are measured and how they are measured. Different magnitude scales are necessary because of differences in earthquakes, the information available, and the purposes for which the magnitudes are used.

A synthetic seismogram is the result of forward modelling the seismic response of an input earth model, which is defined in terms of 1D, 2D or 3D variations in physical properties. In hydrocarbon exploration this is used to provide a 'tie' between changes in rock properties in a borehole and seismic reflection data at the same location. It can also be used either to test possible interpretation models for 2D and 3D seismic data or to model the response of the predicted geology as an aid to planning a seismic reflection survey. In the processing of wide-angle reflection and refraction (WARR) data, synthetic seismograms are used to further constrain the results of seismic tomography. In earthquake seismology, synthetic seismograms are used either to match the predicted effects of a particular earthquake source fault model with observed seismometer records or to help constrain the Earth's velocity structure. Synthetic seismograms are generated using specialized geophysical software.

<span class="mw-page-title-main">Shear wave splitting</span> Physical phenomenon

Shear wave splitting, also called seismic birefringence, is the phenomenon that occurs when a polarized shear wave enters an anisotropic medium. The incident shear wave splits into two polarized shear waves. Shear wave splitting is typically used as a tool for testing the anisotropy of an area of interest. These measurements reflect the degree of anisotropy and lead to a better understanding of the area's crack density and orientation or crystal alignment. We can think of the anisotropy of a particular area as a black box and the shear wave splitting measurements as a way of looking at what is in the box.

Recent advances are improving the speed and accuracy of loss estimates immediately after earthquakes so that injured people may be rescued more efficiently. "Casualties" are defined as fatalities and injured people, which are due to damage to occupied buildings. After major and large earthquakes, rescue agencies and civil defense managers rapidly need quantitative estimates of the extent of the potential disaster, at a time when information from the affected area may not yet have reached the outside world. For the injured below the rubble every minute counts. To rapidly provide estimates of the extent of an earthquake disaster is much less of a problem in industrialized than in developing countries. This article focuses on how one can estimate earthquake losses in developing countries in real time.

An ocean-bottom seismometer (OBS) is a seismometer that is designed to record the earth motion under oceans and lakes from man-made sources and natural sources.

A seismic array is a system of linked seismometers arranged in a regular geometric pattern to increase sensitivity to earthquake and explosion detection. A seismic array differs from a local network of seismic stations mainly by the techniques used for data analysis. The data from a seismic array is obtained using special digital signal processing techniques such as beamforming, which suppress noises and thus enhance the signal-to-noise ratio (SNR).

The upper mantle of Earth is a very thick layer of rock inside the planet, which begins just beneath the crust and ends at the top of the lower mantle at 670 km (420 mi). Temperatures range from approximately 500 K at the upper boundary with the crust to approximately 1,200 K at the boundary with the lower mantle. Upper mantle material that has come up onto the surface comprises about 55% olivine, 35% pyroxene, and 5 to 10% of calcium oxide and aluminum oxide minerals such as plagioclase, spinel, or garnet, depending upon depth.

<span class="mw-page-title-main">Travel-time curve</span>

Travel Time means time for the seismic waves to travel from the focus of an earthquake through the crust to a certain seismograph station. Travel-time curve is a graph showing the relationship between the distance from the epicenter to the observation point and the travel time. Travel-time curve is drawn when the vertical axis of the graph is the travel time and the horizontal axis is the epicenter distance of each observation point.

<span class="mw-page-title-main">MERMAID</span>

MERMAID is a marine scientific instrument platform, short for Mobile Earthquake Recorder for Marine Areas by Independent Divers.

References

  1. Langston, Charles A.; Burdick, L. J. (1977-06-01). "Modeling crustal structure through the use of converted phases in teleseismic body-wave forms". Bulletin of the Seismological Society of America. 67 (3): 677–691. ISSN   0037-1106.
  2. Vinnik, L.P. (September 1977). "Detection of waves converted from P to SV in the mantle". Physics of the Earth and Planetary Interiors. 15 (1): 39–45. doi:10.1016/0031-9201(77)90008-5.
  3. Frederiksen, A. W.; Bostock, M. G. (2000-05-01). "Modelling teleseismic waves in dipping anisotropic structures". Geophysical Journal International. 141 (2): 401–412. doi: 10.1046/j.1365-246x.2000.00090.x . ISSN   0956-540X.
  4. 1 2 3 4 5 Yamauchi, Makiko; Hirahara, Kazuro; Shibutani, Takuo (2003-01-01). "High resolution receiver function imaging of the seismic velocity discontinuities in the crust and the uppermost mantle beneath southwest Japan". Earth, Planets and Space. 55: 59–64. doi: 10.1186/BF03352463 . ISSN   1880-5981.
  5. 1 2 3 Wiemer, S.; Agostinetti, N. Piana; Kissling, E.; Bianchi, I.; Spada, M. (2013-08-01). "Combining controlled-source seismology and receiver function information to derive 3-D Moho topography for Italy". Geophysical Journal International. 194 (2): 1050–1068. doi: 10.1093/gji/ggt148 . ISSN   0956-540X.
  6. 1 2 Langston, Charles A.; Burdick, L. J. (1977-06-01). "Modeling crustal structure through the use of converted phases in teleseismic body-wave forms". Bulletin of the Seismological Society of America. 67 (3): 677–691. ISSN   0037-1106.
  7. 1 2 Rondenay, Stéphane (2009-10-01). "Upper Mantle Imaging with Array Recordings of Converted and Scattered Teleseismic Waves". Surveys in Geophysics. 30 (4): 377–405. doi:10.1007/s10712-009-9071-5. ISSN   1573-0956.
  8. Helffrich, George (2006-02-01). "Extended-time multitaper frequency domain crosscorrelation receiver-function estimation". Bull. Seismol. Soc. Am. 98: 344–347. doi:10.1785/0120050098.