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 generation and 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 microseism, atmospheric, and artificial processes such as explosions and human activities. A related field that uses geology to infer information regarding past earthquakes is paleoseismology. A recording of Earth motion as a function of time, created by a seismograph is called a seismogram. A seismologist is a scientist works in basic or applied 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 mechanical 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">Epicenter</span> Point on the Earths surface that is directly above the hypocentre or focus in an earthquake

The epicenter, epicentre, or epicentrum in seismology is the point on the Earth's surface directly above a hypocenter or focus, the point where an earthquake or an underground explosion originates.

<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 called the Moho discontinuity, Moho boundary, or just Moho – is the boundary between the crust and the mantle of Earth. It is defined by the distinct change in velocity of seismic waves as they pass through changing densities of rock.

Seismic tomography or seismotomography 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">Reflection seismology</span> Explore subsurface properties with seismology

Reflection seismology is a method of exploration geophysics that uses the principles of seismology to estimate the properties of the Earth's subsurface from reflected seismic waves. The method requires a controlled seismic source of energy, such as dynamite or Tovex blast, a specialized air gun or a seismic vibrator. Reflection seismology is similar to sonar and echolocation.

<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 (2003-2018) 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.

<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.

The transition zone is part of the Earth's mantle, and is located between the lower mantle and the upper mantle, between a depth of 410 and 660 km. The Earth's mantle, including the transition zone, consists primarily of peridotite, an ultramafic igneous rock.

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 based 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.

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.

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">MERMAID</span>

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

Seismic wide-angle reflection and refraction is a technique used in geophysical investigations of Earth's crust and upper mantle. It allows the development of a detailed model of seismic velocities beneath Earth's surface well beyond the reach of exploration boreholes. The velocities can then be used, often in combination with the interpretation of standard seismic reflection data and gravity data, to interpret the geology of the subsurface.

<span class="mw-page-title-main">Seismic velocity structure</span> Seismic wave velocity variation

Seismic velocity structure is the distribution and variation of seismic wave speeds within Earth's and other planetary bodies' subsurface. It is reflective of subsurface properties such as material composition, density, porosity, and temperature. Geophysicists rely on the analysis and interpretation of the velocity structure to develop refined models of the subsurface geology, which are essential in resource exploration, earthquake seismology, and advancing our understanding of Earth's geological development.

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. Bibcode:1977BuSSA..67..677B. doi:10.1785/BSSA0670030677. 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. Bibcode:1977PEPI...15...39V. 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. Bibcode:2000GeoJI.141..401F. 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 (1): 59–64. Bibcode:2003EP&S...55...59Y. 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 . hdl: 20.500.11850/70318 . 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. Bibcode:1977BuSSA..67..677B. doi:10.1785/BSSA0670030677. 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. Bibcode:2009SGeo...30..377R. 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 (1): 344–347. Bibcode:2006BuSSA..96..344H. doi:10.1785/0120050098.