Subsurface mapping by ambient noise tomography

Last updated

Subsurface mapping by ambient noise tomography is the mapping underground geological structures under the assistance of seismic signals. Ambient noise, which is not associated with the earthquake, is the background seismic signals. [1] Given that the ambient noises have low frequencies in general, the further classification of ambient noise include secondary microseisms, primary microseisms, and seismic hum, based on different range of frequencies. [2] We can utilize the ambient noise data collected by seismometers (or geophones) to create images for the subsurface under the following processes. Since the ambient noise is considered as diffuse wavefield, we can correlate the filtered ambient noise data from a pair of seismic stations (or seismometers) to find the velocities of seismic wavefields. [3] A 2-dimensional or 3-dimensional velocity map, showing the spatial velocity difference of the subsurface, can thus be created for observing the geological structures. Subsurface mapping by ambient noise tomography can be applied in different fields, such as detecting the underground void space, [4] monitoring landslides, [5] and mapping the crustal and upper mantle structure. [6]

Contents

Characteristics of ambient noise

Characteristic of ambient noise refers to several quantities that can distinguish different ambient noise, including origin, frequency, property, and temporal variation.

Nature of ambient noise

Ambient noise, as the rising star of the seismic source for seismic research other than earthquake, accounts for the naturally and anthropogenically produced seismic vibration of the background. [7] This is different from the active seismic source created solely for seismic research or large seismic source from earthquake. Ocean is the most dominant natural origin of the ambient noise field. [8]

Any seismic source is transmitted as either body waves or surface waves, where ambient noise is no exception. Summary of their properties are shown below.

Seismic Wave
Wave TypeBody/SurfaceDescription
P-wave Body1. Can pass through solid and liquid.

2. Particle movement parallel to wave movement.

S-wave Body1. Can pass through solid only.

2. Particle movement perpendicular to wave movement.

Rayleigh wave Surface1. Include both longitudinal & transverse motions.

2. Amplitude decreases exponentially with increasing distance from surface.

Love wave Surface1. Wave speed lower than P-wave & S-wave, but higher than Rayleigh wave.

2. Horizonal particle movement perpendicular to wave propagation.

Overview Seismic Waves Overview Seismic Waves.jpg
Overview Seismic Waves

The dominance of seismic wave transmission of ambient noise depends on several factors, while the research technique would determine the major type of seismic wave collected for ambient noise. For example, seismologists would often use spatial auto-correlation (SPAC) method which involve the collection and analysis of surface wave. [9]

Frequency of ambient noise

Ambient noise is often known as microseism, where ‘micro’ means very small, [10] and ‘seism’ is an alternative name for earthquake. [11] It can be further classified based on their frequency ranges, namely hum, primary microseism and secondary microseism. [12] The table below shows the comparison of frequency range between the microseisms, arranged from increasing order. Figure a also shows graph of the frequency range of microseisms.

Frequency Range of Microseisms
MicroseismsFrequency Range (Hz)
Hum9×10−3 - 8×10−2
Primary microseism9×10−3 - 8×10−2
Secondary microseism3×10−2 - 1

Origin of ambient noise

Ambient noise can be further classified into two major categories based on the origins of the noise.

Anthropogenic

Anthropogenic ambient noise originates from human activities. Considering the ocean ambient noise source as an example, there are noises that are created unintentionally by human activities, such as shipping and offshore engineering work. [13] During the shipping activity, mechanical waves can be driven up along the water surface and propagate through the ocean. Offshore engineering work can also produce surface waves. Engineering works include but are not limited to borehole drilling, foundation construction and geophysical surveys. [13] Shoreline reclamation has been actively carried out by many countries to create more land for urban development. Those engineering works can thus also be carried out offshore. The processes of offshore drilling and exploration create continuous mechanical waves that can also propagate through the ocean.

In continental urban areas, there are more examples of human activities creating background noise. Other than engineering works, urban traffic is the major component of urban ambient noise. [13] Although the mechanical waves of the continent are not as visible than those from the ocean, they can still be transmitted via the soil and rock layers. Cars travelling on the road can produce repeatable vibration on the road which can then be transmitted through the soil layers.

Natural noise

Natural ambient noise refers to the background noise produced from the natural events. The natural environment is not stationary but constantly changing because nature itself is continuously modified by weather, tectonic movements and biogenic activities. [13] They can also produce low frequency background noise. Some of the most significant events are listed below.

Wind can induce weak ocean waves propagating through the ocean. The varying atmospheric pressure was hypothesized as the origin before but is inadequate to support the existence of all types of microseisms. [2] Instead, ocean waves are proposed as the alternative origin of natural ambient noise. For example, the ocean swells interact with the sea coast to induce hum and primary microseisms, and the interaction of sea waves with opposite direction can produce secondary microseisms. [12]

Variation of ambient noise

To evaluate whether the collected ambient noise source can be further analysed, consider if there are any regular variations or patterns of certain ambient noise source. [13] Referring to the urban noise source, it may experience a daily variation, where the human activities are conducted mostly in daytime and reduced in nighttime. The ambient noise should thus increase in the daytime while reducing at night. Apart from the temporal variation, the spatial variation can also matter. For example, the commercial shipping is usually concentrated on certain routes. The corresponding amplitude of ambient noise should also decrease when moving away from the shipping routes. [13] Nevertheless, it is still difficult to distinguish the ambient noise sources.

Seismic velocity structure modelling

Seismic velocity structure modelling is the modelling technique showing the velocity difference of seismic waves across areas. The modelling process involves some steps, including cross-correlation, Green's function, and inversion. The usage of ambient noise as source of seismic velocity structure modelling rises from 2001 when seismologists tried to correlate the diffused ambient noise wave fields' [14] Velocity structure modelling are complex and require multiple mathematical calculations.

Pre-processing of ambient noise data

The pre-processing of ambient noise data refers to the filtering of the raw data before proceeding to further analysis (cross-correlation, inversion). Raw seismic data can be collected by either geophones, seismometers, or from authorized official bodies. There are also public ambient noise source maps available in recent years. [8] Since the seismometers collect all ambient seismic signals from all directions, the produced seismic waveforms may not reflect the actual background seismic vibrations. Instead, they often contain some occasional seismic signals from earthquakes and other instruments, which is unnecessary in general and thus required to be removed. [15]

Cross-correlation of ambient noise

Figure 1: Simplest situation of ambient noise cross correlation Simplest case of cross correlation.jpg
Figure 1: Simplest situation of ambient noise cross correlation
Figure 2: Resulting Green's functions in different distribution of ambient noise sources Distribution of ambient noise source amended.jpg
Figure 2: Resulting Green's functions in different distribution of ambient noise sources

Ambient noise cross correlation is the process of finding the receiving time lag of ambient noise sources between two nearby stations. Figure 1 illustrates the simplest case of ambient noise cross-correlation. For a pair of receivers (or seismometers or seismic stations) at different locations, the ambient noise signals would be received at a different time, assumed that they travel at the same velocity at the subsurface. The products of cross-correlation of those signals are new seismic waveforms, namely Green's function. Regarding the case with multiple ambient noise sources, the shape of the cross correlation function depends on whether the ambient noise sources are evenly distributed across a certain area. For the most ideal situation where the ambient noise signals are distributed evenly across all directions, the Green’s function would be highly symmetrical (see Figure 2).

Inversion

Inversion is one of the techniques used in ambient noise tomography. Inversion of the Green’s function is used to retrieve the subsurface properties of the Earth, where seismic velocity is one of the important quantities. It is a crucial step in ambient noise tomography. Inversion in seismic analysis can be treated as finding the original factors of the subsurface that induce the current transmission of the ambient noise signals. Inversion of the Green’s function is conducted linearly in the early ambient noise tomography studies with the assumption that the velocity variation is small. [16] Forward model, as the essential process of inversion, is used to estimate the closest quantities of the earth subsurface properties. The cross-correlated seismic waves can be inverted either linearly or non-linearly.

Linkage of geological structure with velocity zone

Before interpreting the velocity zone, it is necessary to understand how the seismic velocity varies. In general, P and S wave travel faster in high density medium. Only P wave can travel through any medium while S wave can only travel through solids. Therefore, a low velocity zone can indicate some vacuum space in the subsurface layer, such as void space and faults. Conversely, a high velocity zone may refer to the lithology with closely packed rocks, such as igneous rock. To correlate the velocity zone with geological structure, it is necessary to consider the size and shape of the velocity zones, and more importantly, the resolution of the subsurface velocity image. The resolution of the image can affect the scale of the subsurface we can interpret. Sometimes, fieldwork is also needed in order to better correlate the velocity map.[ citation needed ]

Example of subsurface structure

Here are some examples of the subsurface structures and features, including but not limited to the following.

See also

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 planetary bodies. It also includes studies of the environmental effects of earthquakes such as tsunamis; other seismic sources such as volcanoes, plate tectonics, glaciers, rivers, oceanic microseisms, and the atmosphere; and artificial processes such as explosions.

Seismic tomography or seismotomography is a technique for imaging the subsurface of the Earth using seismic waves. The properties of seismic waves are modified by the material through which they travel. By comparing the differences in seismic waves recorded at different locations, it is possible to create a model of the subsurface structure. Most commonly, these seismic waves are generated by earthquakes or man-made sources such as explosions. Different types of waves, including P, S, Rayleigh, and Love waves can be used for tomographic images, though each comes with their own benefits and downsides and are used depending on the geologic setting, seismometer coverage, distance from nearby earthquakes, and required resolution. The model created by tomographic imaging is almost always a seismic velocity model, and features within this model may be interpreted as structural, thermal, or compositional variations. Geoscientists apply seismic tomography to a wide variety of settings in which the subsurface structure is of interest, ranging in scale from whole-Earth structure to the upper few meters below the surface.

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

A geophone is a device that converts ground movement (velocity) into voltage, which may be recorded at a recording station. The deviation of this measured voltage from the base line is called the seismic response and is analyzed for structure of the Earth.

Exploration geophysics is an applied branch of geophysics and economic geology, which uses physical methods at the surface of the Earth, such as seismic, gravitational, magnetic, electrical and electromagnetic, to measure the physical properties of the subsurface, along with the anomalies in those properties. It is most often used to detect or infer the presence and position of economically useful geological deposits, such as ore minerals; fossil fuels and other hydrocarbons; geothermal reservoirs; and groundwater reservoirs. It can also be used to detect the presence of unexploded ordnance.

<span class="mw-page-title-main">Magnetotellurics</span> Electromagnetic geophysical technique

Magnetotellurics (MT) is an electromagnetic geophysical method for inferring the earth's subsurface electrical conductivity from measurements of natural geomagnetic and geoelectric field variation at the Earth's surface.

<span class="mw-page-title-main">Seismic source</span> Device that generates controlled seismic energy used for seismic surveys

A seismic source is a device that generates controlled seismic energy used to perform both reflection and refraction seismic surveys. A seismic source can be simple, such as dynamite, or it can use more sophisticated technology, such as a specialized air gun. Seismic sources can provide single pulses or continuous sweeps of energy, generating seismic waves, which travel through a medium such as water or layers of rocks. Some of the waves then reflect and refract and are recorded by receivers, such as geophones or hydrophones.

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

Geophysical imaging is a minimally destructive geophysical technique that investigates the subsurface of a terrestrial planet. Geophysical imaging is a noninvasive imaging technique with a high parametrical and spatio-temporal resolution. It can be used to model a surface or object understudy in 2D or 3D as well as monitor changes.

Seismic migration is the process by which seismic events are geometrically re-located in either space or time to the location the event occurred in the subsurface rather than the location that it was recorded at the surface, thereby creating a more accurate image of the subsurface. This process is necessary to overcome the limitations of geophysical methods imposed by areas of complex geology, such as: faults, salt bodies, folding, etc.

In acoustics, microbaroms, also known as the "voice of the sea", are a class of atmospheric infrasonic waves generated in marine storms by a non-linear interaction of ocean surface waves with the atmosphere. They typically have narrow-band, nearly sinusoidal waveforms with amplitudes up to a few microbars, and wave periods near 5 seconds. Due to low atmospheric absorption at these low frequencies, microbaroms can propagate thousands of kilometers in the atmosphere, and can be readily detected by widely separated instruments on the Earth's surface.

In seismology, a microseism is defined as a faint earth tremor caused by natural phenomena. Sometimes referred to as a "hum", it should not be confused with the anomalous acoustic phenomenon of the same name. The term is most commonly used to refer to the dominant background seismic and electromagnetic noise signals on Earth, which are caused by water waves in the oceans and lakes. Characteristics of microseism are discussed by Bhatt. Because the ocean wave oscillations are statistically homogenous over several hours, the microseism signal is a long-continuing oscillation of the ground. The most energetic seismic waves that make up the microseismic field are Rayleigh waves, but Love waves can make up a significant fraction of the wave field, and body waves are also easily detected with arrays. Because the conversion from the ocean waves to the seismic waves is very weak, the amplitude of ground motions associated to microseisms does not generally exceed 10 micrometers.

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.

In geophysics, geology, civil engineering, and related disciplines, seismic noise is a generic name for a relatively persistent vibration of the ground, due to a multitude of causes, that is often a non-interpretable or unwanted component of signals recorded by seismometers.

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

Interferometry examines the general interference phenomena between pairs of signals in order to gain useful information about the subsurface. Seismic interferometry (SI) utilizes the crosscorrelation of signal pairs to reconstruct the impulse response of a given media. Papers by Keiiti Aki (1957), Géza Kunetz, and Jon Claerbout (1968) helped develop the technique for seismic applications and provided the framework upon which modern theory is based.

<span class="mw-page-title-main">Surface wave inversion</span>

Seismic inversion involves the set of methods which seismologists use to infer properties through physical measurements. Surface-wave inversion is the method by which elastic properties, density, and thickness of layers in the subsurface are obtained through analysis of surface-wave dispersion. The entire inversion process requires the gathering of seismic data, the creation of dispersion curves, and finally the inference of subsurface properties.

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

Michel Campillo is a French seismologist and geophysicist who is currently a professor at Grenoble Alpes University.

Fluvial seismology is the application of seismological methods to understand river processes, such as discharge, erosion, and streambed evolution. Flowing water and the movement of sediments along the streambed generate elastic (seismic) waves that propagate into the surrounding Earth materials. Seismometers can record these signals, which can be analyzed to illuminate different fluvial processes such as turbulent water flow and bedload transport. Seismic methods have been used to observe discharge values that range from single-digits up through tens of thousands of cubic feet per second (cfs).

<span class="mw-page-title-main">MERMAID</span> Marine scientific instrument platform

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

<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. Sleeman, Reinoud. "Ambient Earth noise and instrumental noise" . Retrieved 2023-11-01.
  2. 1 2 Nishida, Kiwamu (2017-08-02). "Ambient seismic wave field". Proceedings of the Japan Academy, Series B. 93 (7): 423–448. Bibcode:2017PJAB...93..423N. doi:10.2183/pjab.93.026. PMC   5713174 . PMID   28769015.
  3. Benson, G. D.; Ritzwoller, M. H.; Barmin, M. P.; Levshin, A. L.; Lin, F.; Moschetti, M. P.; Shapiro, N. M.; Yang, Y. (2007-06-01). "Processing seismic ambient noise data to obtain reliable broad-band surface wave dispersion measurements". Geophysical Journal International. 169 (3): 1239–1260. arXiv: 2007.03374 . Bibcode:2007GeoJI.169.1239B. doi: 10.1111/j.1365-246X.2007.03374.x . S2CID   229068738.
  4. Wang, Yao; Khorrami, Mohammad; Tran, Khiem T.; Horhota, David (2023). "Application of ambient noise tomography for deep void detection". Journal of Applied Geophysics. 209. Bibcode:2023JAG...20904922W. doi:10.1016/j.jappgeo.2022.104922. S2CID   255338249.
  5. Le Breton, Mathieu; Bontemps, Noelie; Guillemot, Antoine; Baillet, Laurent; Larose, Eric (2021-01-28). "Landslide monitoring using seismic ambient noise correlation: challenges and applications". Earth-Science Reviews. 216: 1–26. Bibcode:2021ESRv..21603518L. doi:10.1016/j.earscirev.2021.103518. S2CID   234037825.
  6. Yao, Huajian; Beghein, Caroline; Hilst, Robert D. Van Der (2008-04-01). "Surface wave array tomography in SE Tibet from ambient seismic noise and two-station analysis - II. Crustal and upper-mantle structure". Geophysical Journal International. 173 (1): 205–219. Bibcode:2008GeoJI.173..205Y. doi: 10.1111/j.1365-246X.2007.03696.x . S2CID   29548841.
  7. Oakley, David O. S.; Forsythe, Brandon; Gu, Xin; Nyblade, Andrew A.; Brantley, Susan L. (2021-01-13). "Seismic Ambient Noise Analyses Reveal Changing Temperature and Water Signals to 10s of Meters Depth in the Critical Zone". Journal of Geophysical Research: Earth Surface. 126 (2). Bibcode:2021JGRF..12605823O. doi: 10.1029/2020JF005823 . S2CID   234198739.
  8. 1 2 Igel, Jonas K. H.; Bowden, Daniel C.; Fichtner, Andreas (2023). "SANS: Publicly Available Daily Multi-Scale Seismic Ambient Noise Source Maps". Journal of Geophysical Research: Solid Earth. 128 (1). Bibcode:2023JGRB..12825114I. doi: 10.1029/2022JB025114 . hdl: 20.500.11850/591152 . S2CID   255123230.
  9. Nthaba, Bokani; Ikeda, Tatsunori; Nimiya, Hiro; Tsuji, Takeshi; Lio, Yoshihisa (2022-06-20). "Ambient noise tomography for a high-resolution 3D S-wave velocity model of the Kinki Region, Southwestern Japan, using dense seismic array data". Earth, Planets and Space. 74 (1): 96. Bibcode:2022EP&S...74...96N. doi: 10.1186/s40623-022-01654-x .
  10. "micro". Cambridge Dictionary. Retrieved 2023-11-08.
  11. "seism". Merriam-Webster. Retrieved 2023-11-05.
  12. 1 2 Tanimoto, Toshiro; Anderson, Aaron (2023-09-11). "Seismic noise between 0.003 Hz and 1.0 Hz and its classification". Progress in Earth and Planetary Science. 10 (1): 56. Bibcode:2023PEPS...10...56T. doi: 10.1186/s40645-023-00587-7 .
  13. 1 2 3 4 5 6 Hildebrand, John A. (2023-12-03). "Anthropogenic and natural sources of ambient noise in the ocean". Marine Ecology Progress Series. 395: 5–20. doi: 10.3354/meps08353 .
  14. Sager, Korbinian; Ermert, Laura; Boehm, Christian; Fichtner, Andreas (12 July 2021). "Towards full waveform ambient noise inversion". Geophysical Journal International. 212: 566–590. doi: 10.1093/gji/ggx429 . hdl: 20.500.11850/224948 .
  15. da Silva, Cicero Costa; Poveda, Esteban; Dantas, Renato Ramos da Silva; Julia, Jordi (2021-04-22). "Ambient Noise Tomography with Short-Period Stations: Case Study in the Borborema Province". Pure and Applied Geophysics. 178 (5): 1709–1730. Bibcode:2021PApGe.178.1709D. doi:10.1007/s00024-021-02718-x. S2CID   233330462.
  16. Perez, Ivan Cabrera; D' Auria, Luca; Soubestre, Jean; Barrancos, Jose; Padilla, German D.; Perez, Nemesio M. (2021). "A nonlinear multiscale inversion approach for ambient noise tomography". Geophysical Journal International. 225 (2): 1158–1173. doi: 10.1093/gji/ggaa574 .