Seismic tomography or seismotomography is a technique for imaging the subsurface of the Earth using seismic waves. [1] 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.
Tomography is solved as an inverse problem. Seismic data are compared to an initial Earth model and the model is modified until the best possible fit between the model predictions and observed data is found. Seismic waves would travel in straight lines if Earth was of uniform composition, but structural, chemical, and thermal variations affect the properties of seismic waves, most importantly their velocity, leading to the reflection and refraction of these waves. The location and magnitude of variations in the subsurface can be calculated by the inversion process, although solutions to tomographic inversions are non-unique. Most commonly, only the travel time of the seismic waves is considered in the inversion. However, advances in modeling techniques and computing power have allowed different parts, or the entirety, of the measured seismic waveform to be fit during the inversion. [2] [3] [4]
Seismic tomography is similar to medical x-ray computed tomography (CT scan) in that a computer processes receiver data to produce a 3D image, although CT scans use attenuation instead of travel-time difference. Seismic tomography has to deal with the analysis of curved ray paths which are reflected and refracted within the Earth, and potential uncertainty in the location of the earthquake hypocenter. CT scans use linear x-rays and a known source. [5]
In the early 20th century, seismologists first used travel time variations in seismic waves from earthquakes to make discoveries such as the existence of the Moho [6] and the depth to the outer core. [7] While these findings shared some underlying principles with seismic tomography, modern tomography itself was not developed until the 1970s with the expansion of global seismic networks. Networks like the World-Wide Standardized Seismograph Network were initially motivated by underground nuclear tests, [8] but quickly showed the benefits of their accessible, standardized datasets for geoscience. These developments occurred concurrently with advancements in modeling techniques and computing power that were required to solve large inverse problems [9] [10] and generate theoretical seismograms, [11] which are required to test the accuracy of a model. [12] As early as 1972, [13] researchers successfully used some of the underlying principles of modern seismic tomography to search for fast and slow areas in the subsurface. [14]
The first widely cited publication that largely resembles modern seismic tomography was published in 1976 and used local earthquakes to determine the 3D velocity structure beneath Southern California. [15] [14] The following year, P wave delay times were used to create 2D velocity maps of the whole Earth at several depth ranges, [16] representing an early 3D model. The first model using iterative techniques, which improve upon an initial model in small steps and are required when there are a large number of unknowns, was done in 1984. [17] The model was made possible by iterating upon the first radially anisotropic Earth model, created in 1981. [18] A radially anisotropic Earth model describes changes in material properties, specifically seismic velocity, along a radial path through the Earth, and assumes this profile is valid for every path from the core to the surface. This 1984 study was also the first to apply the term "tomography" to seismology, as the term had originated in the medical field with X-ray tomography. [12]
Seismic tomography has continued to improve in the past several decades since its initial conception. The development of adjoint inversions, which are able to combine several different types of seismic data into a single inversion, help negate some of the trade-offs associated with any individual data type. [12] Historically, seismic waves have been modeled as 1D rays, a method referred to as "ray theory" that is relatively simple to model and can usually fit travel-time data well. [19] However, recorded seismic waveforms contain much more information than just travel-time and are affected by a much wider path than is assumed by ray theory. Methods like the finite-frequency method attempt to account for this within the framework of ray theory. [20] More recently, the development of "full waveform" or "waveform" tomography has abandoned ray theory entirely. This method models seismic wave propagation in its full complexity and can yield more accurate images of the subsurface. Originally these inversions were developed in exploration seismology [21] in the 1980s and 1990s and were too computationally complex for global and regional scale studies, [12] but development of numerical modeling methods to simulate seismic waves [22] has allowed waveform tomography to become more common.
Seismic tomography uses seismic records to create 2D and 3D models of the subsurface through an inverse problem that minimizes the difference between the created model and the observed seismic data. [14] Various methods are used to resolve anomalies in the crust, lithosphere, mantle, and core based on the availability of data and types of seismic waves that pass through the region. Longer wavelengths penetrate deeper into the Earth, but seismic waves are not sensitive to features significantly smaller than their wavelength and therefore provide a lower resolution. Different methods also make different assumptions, which can have a large effect on the image created. For example, commonly used tomographic methods work by iteratively improving an initial input model, and thus can produce unrealistic results if the initial model is unreasonable. [14]
P wave data are used in most local models and global models in areas with sufficient earthquake and seismograph density. S and surface wave data are used in global models when this coverage is not sufficient, such as in ocean basins and away from subduction zones. First-arrival times are the most widely used, but models utilizing reflected and refracted phases are used in more complex models, such as those imaging the core. Differential traveltimes between wave phases or types are also used.
Local tomographic models are often based on a temporary seismic array targeting specific areas, unless in a seismically active region with extensive permanent network coverage. These allow for the imaging of the crust and upper mantle.
Regional to global scale tomographic models are generally based on long wavelengths. Various models have better agreement with each other than local models due to the large feature size they image, such as subducted slabs and superplumes. The trade off from whole mantle to whole Earth coverage is the coarse resolution (hundreds of kilometers) and difficulty imaging small features (e.g. narrow plumes). Although often used to image different parts of the subsurface, P and S wave derived models broadly agree where there is image overlap. These models use data from both permanent seismic stations and supplementary temporary arrays.
Seismic tomography can resolve anisotropy, anelasticity, density, and bulk sound velocity. [30] Variations in these parameters may be a result of thermal or chemical differences, which are attributed to processes such as mantle plumes, subducting slabs, and mineral phase changes. Larger scale features that can be imaged with tomography include the high velocities beneath continental shields and low velocities under ocean spreading centers. [31]
The mantle plume hypothesis proposes that areas of volcanism not readily explained by plate tectonics, called hotspots, are a result of thermal upwelling within the mantle. Some researchers have proposed an upper mantle source above the 660km discontinuity for these plumes, [32] while others propose a much deeper source, possibly at the core-mantle boundary. [33]
While the source of mantle plumes has been highly debated since they were first proposed in the 1970s, [34] most modern studies argue in favor of mantle plumes originating at or near the core-mantle boundary. [35] This is in large part due to tomographic images that reveal both the plumes themselves [36] [37] as well as large low-velocity zones in the deep mantle that likely contribute to the formation of mantle plumes. These large low-shear velocity provinces as well as smaller ultra low velocity zones have been consistently observed across many tomographic models of the deep Earth [38]
Subducting plates are colder than the mantle into which they are moving. This creates a fast anomaly that is visible in tomographic images. Tomographic images have been made of most subduction zones around the world and have provided insight into the geometries of the crust and upper mantle in these areas. These images have revealed that subducting plates vary widely in how steeply they move into the mantle. [39] [40] Tomographic images have also seen features such as deeper portions of the subducting plate tearing off from the upper portion. [41]
Tomography can be used to image faults to better understand their seismic hazard. This can be through imaging the fault itself by seeing differences in seismic velocity across the fault boundary [42] or by determining near-surface velocity structure, [43] which can have a large impact on the magnitude on the amplitude of ground-shaking during an earthquake due to site amplification effects. [44] Near-surface velocity structure from tomographic images can also be useful for other hazards, such as monitoring of landslides for changes in near-surface moisture content which has an effect on both seismic velocity and potential for future landslides. [45] [46]
Tomographic images of volcanoes have yielded new insights into properties of the underlying magmatic system. These images have most commonly been used to estimate the depth and volume of magma stored in the crust, [47] [48] but have also been used to constrain properties such as the geometry, temperature, or chemistry of the magma. [49] [50] [51] It is important to note that both lab experiments and tomographic imaging studies have shown that recovering these properties from seismic velocity alone can be difficult due to the complexity of seismic wave propagation through focused zones of hot, potentially melted rocks. [52] [53]
While comparatively primitive to tomography on Earth, seismic tomography has been proposed on other bodies in the solar system and successfully used on the Moon. Data collected from four seismometers placed by the Apollo missions have been used many times to create 1-D velocity profiles for the moon, [54] [55] [56] and less commonly 3-D tomographic models. [57] Tomography relies on having multiple seismometers, but tomography-adjacent methods for constraining Earth structure have been used on other planets. While on Earth these methods are often used in combination with seismic tomography models to better constrain the locations of subsurface features, [58] [59] they can still provide useful information about the interiors of other planetary bodies when only a single seismometer is available. For example, data gathered by the SEIS (Seismic Experiment for Interior Structure) instrument on InSight [60] on Mars has been able to detect the Martian core. [61]
Global seismic networks have expanded steadily since the 1960s, but are still concentrated on continents and in seismically active regions. Oceans, particularly in the southern hemisphere, are under-covered. [62] Temporary seismic networks have helped improve tomographic models in regions of particular interest, but typically only collect data for months to a few years. The uneven distribution of earthquakes biases tomographic models towards seismically active regions. Methods that do not rely on earthquakes such as active source surveys or ambient noise tomography have helped image areas with little to no seismicity, though these both have their own limitations as compared to earthquake-based tomography.
The type of seismic wave used in a model limits the resolution it can achieve. Longer wavelengths are able to penetrate deeper into the Earth, but can only be used to resolve large features. Finer resolution can be achieved with surface waves, with the trade off that they cannot be used in models deeper than the crust and upper mantle. The disparity between wavelength and feature scale causes anomalies to appear of reduced magnitude and size in images. P and S wave models respond differently to the types of anomalies. Models based solely on the wave that arrives first naturally prefer faster pathways, causing models based on these data to have lower resolution of slow (often hot) features. [14] This can prove to be a significant issue in areas such as volcanoes where rocks are much hotter than their surroundings and oftentimes partially melted. [63] Shallow models must also consider the significant lateral velocity variations in continental crust.
Because seismometers have only been deployed in large numbers since the late-20th century, tomography is only capable of viewing changes in velocity structure over decades. For example, tectonic plates only move at millimeters per year, so the total amount of change in geologic structure due to plate tectonics since the development of seismic tomography is several orders of magnitude lower than the finest resolution possible with modern seismic networks. [64] However, seismic tomography has still been used to view near-surface velocity structure changes at time scales of years to months. [65] [45]
Tomographic solutions are non-unique. Although statistical methods can be used to analyze the validity of a model, unresolvable uncertainty remains. [14] This contributes to difficulty comparing the validity of different model results.
Computing power limits the amount of seismic data, number of unknowns, mesh size, and iterations in tomographic models. This is of particular importance in ocean basins, which due to limited network coverage and earthquake density require more complex processing of distant data. Shallow oceanic models also require smaller model mesh size due to the thinner crust. [12]
Tomographic images are typically presented with a color ramp representing the strength of the anomalies. This has the consequence of making equal changes appear of differing magnitude based on visual perceptions of color, such as the change from orange to red being more subtle than blue to yellow. The degree of color saturation can also visually skew interpretations. These factors should be considered when analyzing images. [5]
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.
Geophysics is a subject of natural science concerned with the physical processes and physical properties of the Earth and its surrounding space environment, and the use of quantitative methods for their analysis. Geophysicists, who usually study geophysics, physics, or one of the Earth sciences at the graduate level, complete investigations across a wide range of scientific disciplines. The term geophysics classically refers to solid earth applications: Earth's shape; its gravitational, magnetic fields, and electromagnetic fields ; its internal structure and composition; its dynamics and their surface expression in plate tectonics, the generation of magmas, volcanism and rock formation. However, modern geophysics organizations and pure scientists use a broader definition that includes the water cycle including snow and ice; fluid dynamics of the oceans and the atmosphere; electricity and magnetism in the ionosphere and magnetosphere and solar-terrestrial physics; and analogous problems associated with the Moon and other planets.
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.
A mantle plume is a proposed mechanism of convection within the Earth's mantle, hypothesized to explain anomalous volcanism. Because the plume head partially melts on reaching shallow depths, a plume is often invoked as the cause of volcanic hotspots, such as Hawaii or Iceland, and large igneous provinces such as the Deccan and Siberian Traps. Some such volcanic regions lie far from tectonic plate boundaries, while others represent unusually large-volume volcanism near plate boundaries.
The Iceland hotspot is a hotspot which is partly responsible for the high volcanic activity which has formed the Iceland Plateau and the island of Iceland. It contributes to understanding the geological deformation of Iceland.
The East African Rift (EAR) or East African Rift System (EARS) is an active continental rift zone in East Africa. The EAR began developing around the onset of the Miocene, 22–25 million years ago. It was formerly considered to be part of a larger Great Rift Valley that extended north to Asia Minor.
The EarthScope project (2003-2018) was an National Science Foundation (NSF) funded Earth science program using 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.
A slow earthquake is a discontinuous, earthquake-like event that releases energy over a period of hours to months, rather than the seconds to minutes characteristic of a typical earthquake. First detected using long term strain measurements, most slow earthquakes now appear to be accompanied by fluid flow and related tremor, which can be detected and approximately located using seismometer data filtered appropriately. That is, they are quiet compared to a regular earthquake, but not "silent" as described in the past.
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.
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.
The Society hotspot is a volcanic hotspot in the south Pacific Ocean which is responsible for the formation of the Society Islands, an archipelago of fourteen volcanic islands and atolls spanning around 720 kilometres (450 mi) of the ocean which formed between 4.5 and <1 Ma.
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.
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.
Crustal recycling is a tectonic process by which surface material from the lithosphere is recycled into the mantle by subduction erosion or delamination. The subducting slabs carry volatile compounds and water into the mantle, as well as crustal material with an isotopic signature different from that of primitive mantle. Identification of this crustal signature in mantle-derived rocks is proof of crustal recycling.
Large low-shear-velocity provinces (LLSVPs), also called large low-velocity provinces (LLVPs) or superplumes, are characteristic structures of parts of the lowermost mantle, the region surrounding the outer core deep inside the Earth. These provinces are characterized by slow shear wave velocities and were discovered by seismic tomography of deep Earth. There are two main provinces: the African LLSVP and the Pacific LLSVP, both extending laterally for thousands of kilometers and possibly up to 1,000 kilometres vertically from the core–mantle boundary. These have been named Tuzo and Jason respectively, after Tuzo Wilson and W. Jason Morgan, two geologists acclaimed in the field of plate tectonics. The Pacific LLSVP is 3,000 kilometers across and underlies four hotspots on Earth's crust that suggest multiple mantle plumes underneath. These zones represent around 8% of the volume of the mantle, or 6% of the entire Earth.
The Hellenic subduction zone (HSZ) is the convergent boundary between the African plate and the Aegean Sea plate, where oceanic crust of the African is being subducted north–northeastwards beneath the Aegean. The southernmost and shallowest part of the zone is obscured beneath the deformed thick sedimentary sequence that forms the Mediterranean Ridge accretionary complex. It has a well-defined Wadati–Benioff zone of seismicity, which demonstrates the relatively shallow dip of its southern part, which increases markedly to the north of the non-volcanic part of the Hellenic arc. The descending slab has been imaged using seismic tomography down to the top of the mantle transition zone at 410 km depth.
Intraplate volcanism is volcanism that takes place away from the margins of tectonic plates. Most volcanic activity takes place on plate margins, and there is broad consensus among geologists that this activity is explained well by the theory of plate tectonics. However, the origins of volcanic activity within plates remains controversial.
MERMAID is a marine scientific instrument platform, short for Mobile Earthquake Recorder for Marine Areas by Independent Divers.
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. 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. We can utilize the ambient noise data collected by seismometers 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 to find the velocities of seismic wavefields. 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, monitoring landslides, and mapping the crustal and upper mantle structure.
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.
{{cite web}}
: CS1 maint: bot: original URL status unknown (link)