Academic literature on the topic 'Crustal scattering'

Eastern Indonesia is tectonically complex, formed by different plates and microplates interactions from different origins. This complexity gives geoscientists a challenge to solve the ’jigsaw’ of the complex interactions. The understanding of tectonic processes can lead to a breakthrough in both resource exploration and disaster risk reduction. We utilize teleseismic P wave coda for random coda from scattering and deterministic coda originated from the crust-mantle boundary (Moho) to derive the crustal properties, including thickness, Vp/Vs, and qualitative scattering characteristics. For the scattering properties, we apply Iterative Cross-Correlation and Stacking (ICCS) to align the waveform. At the same time, for the crust characteristic, we employ the Receiver Functions (RF) method alongside H-k stacking. The crustal thickness recovered from the RF and H-k stacking has a good correlation with the crustal origin, where the thickness in older and stable crust originated from Sundaland and Gondwana is thicker than a younger plate of the crust arc and subduction origin. The Vp/Vs is high in a region that is interpreted to be dominated by mafic lower crust originated from oceanic-oceanic subduction during Eocene, anisotropy, or by a magmatic anomaly. The P coda also correlated well with the subsurface magmatic anomaly by providing a unique pattern.

SUMMARY Information on fault zone structure is essential for our understanding of earthquake mechanics, continental deformation and seismic hazard. We use the scattered seismic wavefield to study the subsurface structure of the North Anatolian Fault Zone (NAFZ) in the region of the 1999 İzmit and Düzce ruptures using data from an 18-month dense deployment of seismometers with a nominal station spacing of 7 km. Using the forward- and back-scattered energy that follows the direct P-wave arrival from teleseismic earthquakes, we apply a scattered wave inversion approach and are able to resolve changes in lithospheric structure on a scale of 10 km or less in an area of about 130 km by 100 km across the NAFZ. We find several crustal interfaces that are laterally incoherent beneath the surface strands of the NAFZ and evidence for contrasting crustal structures either side of the NAFZ, consistent with the presence of juxtaposed crustal blocks and ancient suture zones. Although the two strands of the NAFZ in the study region strike roughly east–west, we detect strong variations in structure both north–south, across boundaries of the major blocks, and east–west, parallel to the strike of the NAFZ. The surface expression of the two strands of the NAFZ is coincident with changes on main interfaces and interface terminations throughout the crust and into the upper mantle in the tomographic sections. We show that a dense passive network of seismometers is able to capture information from the scattered seismic wavefield and, using a tomographic approach, to resolve the fine scale structure of crust and lithospheric mantle even in geologically complex regions. Our results show that major shear zones exist beneath the NAFZ throughout the crust and into the lithospheric mantle, suggesting a strong coupling of strain at these depths.

SUMMARY Seismic attenuation is measured from a swarm of 50 earthquakes in Kīlauea volcano in 2018, associated with caldera collapse. The traverse extends at nearly constant azimuth to the saddle between Mauna Loa and Mauna Kea, continuing to Maui beneath the distal flanks of three dormant volcanos. From Maui the traverse then extends seaward to the Aloha Cabled Observatory (ACO) on the seafloor north of O‘ahu. The effective attenuation is measured with respect to an ${\omega ^{ - 2}}$ earthquake source model. Frequency dependent ${Q_P}$ and ${Q_S}$ are derived. The initial path is shallow and uphill, the path to Maui propagates at mid-crustal depths, and the path to ACO extends through oceanic crust. The observations of ${Q_P} \le {Q_S}$ over each traverse are modelled as bulk attenuation ${Q_K}$. Several attenuation processes are observed, including ${Q_\mu }$, ${Q_K}$, $Q\sim f$, constant Q and scattering. The observation of bulk attenuation is ascribed to contrasting physical properties between basalt and water saturated vesicles. The ratio of Q values between shallow and mid-crustal propagation is used to derive an activation energy E* for the undetermined shear attenuation mechanism. A Debye relaxation peak is fit to the ${Q_S}( f )$ and ${Q_K}( f )$ observed for the mid-crustal pathway. A prior high-frequency attenuation study near Wake Island compares well with this Hawaiian Q data set, which in general shows lower values of Q than observed for Wake.

SUMMARY Seismic attenuation across the US is estimated using station ML magnitude data from the USArray. Station magnitudes are recalibrated back to amplitude and back projected in a 2-D tomography. Data represent the amplitudes of the horizontal components of the Lg phase. The western US shows regions of very high attenuation and contrasts with the lesser attenuation of the eastern US. Individual attenuation anomalies can be clearly tied to regional geology. Station gains show broad regional variations that match geographic regions. Most of the high-attenuation areas are regions of high geothermal activity suggesting that intrinsic attenuation dominates over scattering attenuation. An exception is the central San Andreas Fault zone because it lacks any localized heat-flow anomaly. The US east of the Rocky Mountains is bland and contains none of the high-attenuation regions of the western US. Instead, the central US has low-attenuation patches that do not obviously correspond to geologic province. Sediments of the Gulf Coast Plain, Willison Basin and Michigan Basin do show up as intermediate attenuation while the Illinois Basin, Appalachian Basin and other basins are not apparent. In Alaska, attenuation is generally less than the western US, but still much greater than the eastern US. In southeast Alaska, the Wrangell Volcanic Field causes a sizeable high-attenuation zone. The volcanic Aleutian Mountains also have high attenuation. However, moderate to high attenuation also correlates with the tertiary sedimentary basins in Alaska. The North Slope Basin does not seem to attenuate. Thicker crust and mountain roots tend to show less attenuation, if anything, but this correspondence is most likely due to differences in temperature and seismic velocity. Heat, scattering and young sedimentary basins create seismic attenuation in the continental crust.

In this study Monte Carlo solutions to the radiative transfer equations are used to model translational and rotational motion seismogram envelopes in random elastic media with deterministic background structure assuming multiple anisotropic scattering. The results of the Monte Carlo radiative transfer theory simulations are verified by comparisons with 3D full wave field finite difference simulations. The observation and modeling of the three additional components of rotational ground motions can provide independent information about seismic wave propagation in the Earth’s structure. Rotational motions around the vertical axis observed in the P-wave coda are of particular interest as they can only be excited by horizontally polarized shear waves and therefore indicate the conversion from P- to SH-energy by multiple scattering at 3D-heterogeneities. Scattering and attenuation parameters in south-east Germany beneath the Gräfenberg array and in the Vogtland region are estimated by comparisons of synthesized multi-component seismogram envelopes to seismic data from local and regional swarm earthquakes and to teleseismic events. In a first step, frequency dependent scattering and attenuation parameters from a local data set are estimated for the Vogtland region using nearby swarm earthquakes. The results from the elastic simulations are compared to outcomes from acoustic radiative transfer simulations. Both methods yield similar results and suggest that intrinsic attenuation dominates scattering attenuation. From the elastic simulations it is observable, that forward scattering is required to explain the data. However, the amount of forward scattering strength remains unresolvable. In a second step scattering and attenuation parameters beneath the Gräfenberg array are estimated using a nonlinear genetic inversion of seismogram envelopes from regional events at high frequencies (4–8 Hz). The preferred model of crustal heterogeneity consists of a random medium described by an exponential auto correlation function with a transport mean free path length of ∼ 420 km. The quality factor for elastic S-wave attenuation Q iS is around 700. In a final step simulations of teleseismic P-wave arrivals, using this estimated set of scattering and attenuation parameters, are compared to observed seismogram envelopes from deep events. Simulations of teleseismic events with the parameters found from the regional inversion show good agreement with the measured seismogram envelopes. This includes ringlaser observations of vertical rotations in the teleseismic P-wave coda that naturally result from the proposed model of wave scattering. The model also predicts, that the elastic energy recorded in the teleseismic P-coda is not equipartitioned, unlike the coda of regional events, but contains an excess of shear energy. The combined results from the three different data sets suggest that scattering generating the seismic coda mainly occurs in the crustal part of the lithosphere beneath the receivers. The observations do not require scattering of high frequency waves in the mantle, but weak scattering in the lithospheric mantle cannot be ruled out
In dieser Studie werden Monte Carlo Lösungen für die Energietransfergleichungen genutzt, um Seismogrammeinhüllende von Translations- und Rotationsbewegungen zu modellieren. Die Ergebnisse der Monte Carlo Simulationen werden durch einen Vergleich mit 3D finiten Differenzen Simulationen verifiziert. Diese Modellierung findet in einem elastischen Zufallsmedium mit deterministischer Hintergrundstruktur unter Annahme multipler anisotroper Streuung statt. Die Beobachtung und Modellierung der drei zusätzlichen Komponenten der Rotationsbodenbewegungen kann unabhängige Informationen über die Ausbreitung seismischer Wellen im Erdkörper liefern. Rotationsbewegungen um die vertikale Achse in der P-Wellen Koda sind in diesem Zusammenhang von speziellem Interesse, da sie nur von horizontal polarisierten Scheerwellen angeregt werden können. Die gemessenen Rotationsbewegungen deuten folglich auf Konversionen von P- zu SH-Energie durch multiple Streuung an 3D-Heterogenitäten hin. Für die Bestimmung von Streu- und Dämpfungsparametern im Südosten Deutschlands (Gräfenberg Array, Vogtland) werden synthetisch erzeugte, mehrspurige Seismogrammeinhüllende mit Daten lokaler und regionaler Schwarmbeben und teleseismicher Ereignisse verglichen. In einem ersten Schritt werden frequenzabhängige Krustenparameter für die Vogtlandregion mittels eines lokalen Datensatzes von nahen Schwarmbeben bestimmt. Die Resultate mittels elastischer Energietransfertheorie werden mit Ergebnissen aus Simulationen mittels akustischer Energietransfertheorie verglichen. Beide Methoden liefern ähnliche Parameter und sagen einen größeren Einfluss der intrinsichen Dämpfung im Vergleich zur Streudämpfung voraus. Aus den elastischen Simulationen geht hervor, dass für die Beschreibung der Daten Vorwärtsstreung angenommen werden muss, die Stärke dieser lässt sich jedoch nicht auflösen. In einem zweiten Schritt werden die Streu- und Dämpfungseigenschaften der Erdkruste im Untergrund des Gräfenberg Arrays untersucht. Hierzu wird eine nicht-lineare genetische Inversion von Seismogrammeinhüllenden regionaler Ereignisse bei hohen Frequenzen (4–8 Hz) verwendet. Das bevorzugte Modell der Krustenheterogenität wird durch ein exponentielles Zufallsmedium, einer mittleren freien Transportweglänge von ca. 420 km und einem Qualitätsfaktor für S-Wellen Q iS von ca. 700 beschrieben. Ein letzter Schritt vergleicht Simulationen von teleseismischen P-Welleneinsätzen mit beobachteten Seismogrammeinüllenden von tiefen Erdbeben unter der Nutzung der Parameter aus der regionalen Inversion. Die Simulationen der teleseismischen Ereignisse mit den Parametern der regionalen Inversion zeigen eine gute Übereinstimmung mit den gemessenen Seismogrammeinhüllenden. Dieser Vergleich beinhaltet Ringlaserbeobachtungen der Rotationsbewegungen um die vertikale Achse, welche aus dem angenommenen Streumodell resultieren. Das Modell sagt voraus, dass die elastische Energie in der teleseismischen P-Wellen Koda im Gegensatz zur Koda lokaler oder regionaler Ereignisse nicht gleichverteilt ist, sondern einen Überschuss an Scheerenergie beinhaltet. Die Resultate aus den Untersuchungen der lokalen, regionalen und teleseismischen Datensätze zeigen, dass die Streuereignisse, welche die seismische Koda erklären, hauptsächlich in der Kruste unterhalb der seismischen Empfänger stattfinden. Streuung des Wellenfeldes im Mantel wird für die Erklärung der Daten nicht benötigt, schwache Streuung im lithosphärischen Mantel kann jedoch nicht ausgeschlossen werden

Thesis (Ph.D.)--Indiana University, Dept. of Physics, 2008.
Title from PDF t.p. (viewed on Jul 29, 2009). Source: Dissertation Abstracts International, Volume: 69-12, Section: B, page: 7583. Adviser: Charles J. Horowitz.

Array seismology is an useful tool to perform a detailed investigation of the Earth’s interior. Seismic arrays by using the coherence properties of the wavefield are able to extract directivity information and to increase the ratio of the coherent signal amplitude relative to the amplitude of incoherent noise. The Double Beam Method (DBM), developed by Krüger et al. (1993, 1996), is one of the possible applications to perform a refined seismic investigation of the crust and mantle by using seismic arrays. The DBM is based on a combination of source and receiver arrays leading to a further improvement of the signal-to-noise ratio by reducing the error in the location of coherent phases. Previous DBM works have been performed for mantle and core/mantle resolution (Krüger et al., 1993; Scherbaum et al., 1997; Krüger et al., 2001). An implementation of the DBM has been presented at 2D large-scale (Italian data-set for Mw=9.3, Sumatra earthquake) and at 3D crustal-scale as proposed by Rietbrock & Scherbaum (1999), by applying the revised version of Source Scanning Algorithm (SSA; Kao & Shan, 2004). In the 2D application, the rupture front propagation in time has been computed. In 3D application, the study area (20x20x33 km3), the data-set and the source-receiver configurations are related to the KTB-1994 seismic experiment (Jost et al., 1998). We used 60 short-period seismic stations (200-Hz sampling rate, 1-Hz sensors) arranged in 9 small arrays deployed in 2 concentric rings about 1 km (A-arrays) and 5 km (B-array) radius. The coherence values of the scattering points have been computed in the crustal volume, for a finite time-window along all array stations given the hypothesized origin time and source location. The resulting images can be seen as a (relative) joint log-likelihood of any point in the subsurface that have contributed to the full set of observed seismograms.
Università di Bologna
Unpublished
6T. Sismicità indotta e caratterizzazione sismica dei sistemi naturali
open

The Geoscience Laser Ranging System (GLRS), a spaceborne picosecond laser ranging and altimetry system, has been proposed as an instrument on the Earth Observing System (EOS), an unmanned remote sensing platform to be launched by NASA in the 1990s. In the ranging mode, the GLRS instrument measures the times-of-flight between the spacecraft and an array of ground-based retroreflectors. The total data set of individual range measurements can be processed to yield the 3-D relative coordinates of each target with centimeter accuracies. By revisiting the target sites at periodic intervals on subsequent orbits, relative motions can be observed. Applications include the measurement of tectonic plate or ice sheet velocities, crustal deformation, ground subsidence caused by oil or water withdrawal, large scale engineering projects (highways, pipelines, etc.), and survey densification. In the altimetry mode, GLRS measures time-of-flight to an underlying diffuse scattering surface. The altimetry capabilities of GLRS are greatly enhanced by its unique ability to provide centimeter accuracy updates to its own orbital ephemeris while retroranging. Potential applications include measurement of terrain profiles, biomass (tree canopy heights), ice-sheet thickness, marine geoid, atmospheric surface pressure over oceans, and oceanographic features (wave height, surface roughness, etc.). The presentation reviews the design, applications, simulations, hardware status, and instrument heritage of the GLRS system.