abstract: This research investigates the fine scale structure in Earth's mantle, especially for the lowermost mantle, where strong heterogeneity exists. Recent seismic tomography models have resolved large-scale features in the lower mantle, such as the large low shear velocity provinces (LLSVPs). However, differences are present between different models, especially at shorter length scales. Fine scale structures both within and outside LLSVPs are still poorly constrained. The drastic growth of global seismic networks presents densely sampled seismic data in unprecedented quality and quantity. In this work, the Empirical Wavelet construction method has been developed to document seismic travel time and waveform information for a global shear wave seismic dataset. A dataset of 250K high-quality seismic records with comprehensive measurements is documented and made publicly available. To more accurately classify high quality seismic signal from the noise, 1.4 million manually labeled seismic records have been used to train a supervised classification model. The constructed model performed better than the empirical model deployed in the Empirical Wavelet method, with 87% in precision and 83% in recall. To utilize lower amplitude phases such as higher multiples of S and ScS waves, we have developed a geographic bin stacking method to improve signal-to-noise ratio. It is then applied to Sn waves up to n=6 and ScSn wave up to n=5 for both minor and major arc phases. The virtual stations constructed provide unique path sampling and coverage, vastly improving sampling in the Southern Hemisphere. With the high-quality dataset we have gathered, ray-based layer stripping iterative forward tomography is implemented to update a starting tomography model by mapping the travel time residuals along the ray from the surface down to the core mantle boundary. Final updated models with different starting tomography models show consistent updates, suggesting a convergent solution. The final updated models show higher resolution results than the starting tomography models, especially on intermediate-scale structures. The combined analyses and results in this work provide new tools and new datasets to image the fine-scale heterogeneous structures in the lower mantle, which advances our understanding of the dynamics and evolution of the Earth's mantle. / Dissertation/Thesis / Doctoral Dissertation Geological Sciences 2019

Mining induced seismicity can lead to bumps which cause problems at many mines within the United States and around the world. This seismicity, often referred to as bumps or bursts, can result in injuries, fatalities, and expensive capital damage and production interruptions. There are many factors that contribute to mining induced seismicity but there is still no concrete method to forecast future seismic activity around a mine. One of the main precursors to large seismic events is an increase in situ stress. One way to find areas within geological strata that are highly stressed is to measure p-wave propagation velocities. High p-wave propagation velocities are associated with high in-situ stress levels. By using tomography programs a three-dimensional velocity model can be constructed. When seismic activity is present the event arrival times at each geophone, the locations of each geophone, and the three dimensional velocity model are used in conjunction with one another to locate the seismic events. This research compares the locations of seismic events from a deep coal mine longwall panel in the western United States with the associated p-wave propagation velocities from the previous 24 hours. The aim of this comparison is to provide a link between high velocities and seismic activity that could potentially be used to forecast future seismic activity. The comparison is completed both qualitatively through the use of a visual analysis, and quantitatively using various numerical and correlation comparisons on the seismic and velocity data. The qualitative comparison is completed using the event locations from the tomography program SIMULPS. The quantitative comparison is completed twice using two different tomography programs, SIMULPS and TomoDD, which use different methods for locating the seismic events. Before these comparisons were completed the stresses around the longwall panel were first modeled using the boundary element modeling program LAMODEL to study the effects of three backfilled cross panel entries which were located ahead of the mining face. The modeling showed similar vertical stress distributions as a panel without cross panel entries but higher stress magnitudes. The qualitative analysis involved comparing tomograms created with SIMULPS with seismicity plots from the following day. One noticeable feature of these tomograms is the presence of a stressed area directly ahead of the face. This stressed area represents the forward abutment. The results of this qualitative analysis illustrate a correlation between high p-wave velocities and seismic activity 24 hours later for several of the days studied. The other days showed little to no correlation. Additionally, not all high p-wave velocity regions resulted in seismic activity. Due to these inconsistencies visually analyzing velocity plots obtained from the program SIMULPS is not a reliable way to forecast the locations of seismic activity 24 hours later. The result of the quantitative comparisons completed with the programs SIMULPS and TomoDD further highlighted inconsistencies in the correlation between high p-wave velocities and associated seismic activity 24 hours later. TomoDD provided better correlation values than SIMULPS and generally showed that as the level of seismicity increased the p-wave propagation velocities 24 hours prior also increased. Although TomoDD provided good correlations for some of the data pairs studied, the overall inconsistencies prompt the need for further study in this area using TomoDD to find the optimal forecasting time period. / Master of Science

Thesis (M.S.)--West Virginia University, 2005. / Title from document title page. Document formatted into pages; contains xi, 145 p. : ill. (some col.), maps (some col.). Includes abstract. Includes bibliographical references (p. 141-145).

Seismic tomography is a technique for illuminating the physical dynamics of the Earth by seismic waves generated by earthquakes or explosions. In both industry and academia, the seismic exploration does not yet have the capability of imaging seismic tomography in real-time and with high resolution. There are two reasons. First, at present raw seismic data are typically recorded on sensor nodes locally then are manually collected to central observatories for post processing, and this process may take months to complete. Second, high resolution tomography requires a large and dense sensor network, the real-time data retrieval from a network of large-amount wireless seismic nodes to a central server is virtually impossible due to the sheer data amount and resource limitations. This limits our ability to understand earthquake zone or volcano dynamics. To obtain the seismic tomography in real-time and high resolution, a new design of sensor network system for raw seismic data processing and distributed tomography computation is demanded. Based on these requirements, three research aspects are addressed in this work. First, a distributed multi-resolution evolving tomography computation algorithm is proposed to compute tomography in the network, while avoiding costly data collections and centralized computations. Second, InsightTomo, an end-to-end sensor network emulation platform, is designed to emulate the entire process from data recording to tomography image result delivery. Third, a sensor network testbed is presented to verify the related methods and design in real world. The design of the platform consists of hardware, sensing and data processing components.

The largest earthquakes and the majority of the seismic energy are released on megathrust faults in subduction zones. The goal of this dissertation is to characterize the seismic behavior, structural and physical properties of the megathrust fault, so that we can better understand the controls on slip behavior and large earthquakes. To address this goal, I analyzed seismicity data collected by a local seismic network deployed in southern Alaska and multi-channel seismic (MCS) data from an active-source survey offshore of the Alaska Peninsula. This dissertation work revealed seismicity patterns associated with a large asperity, downdip transitions in megathrust fault structure, and along-strike variations in the properties of subducting sediment on the shallow part of the subduction zone. All of these observations have important implications for seismic behavior of the megathrust.

[Truncated abstract] Seismic tomography has been extensively used in geophysics for different purposes such as geological mapping and prospecting for oil and gas. In geophysics, ultrasound or electromagnetic waves are normally used to provide the tomographic information. In the geotechnical area, seismic tomography is emerging as a promising technique that can be used to determine the spatial variability of shear wave velocities and hence the small strain stiffness of geomaterials. Although some studies have been undertaken to incorporate seismic measurement into centrifuge modelling, there has been to date no attempt to build a complete seismic tomography facility with high resolution for use in a geotechnical centrifuge. Such a powerful facility can help in better understanding of soil behaviour by providing a complete picture of the spatial variation of the soil property of concern. The main aim of this study was to develop a high-resolution seismic tomography (ST) system for the beam centrifuge at the University of Western Australia (UWA) by which the shear wave velocity and hence maximum shear modulus could be determined anywhere in the centrifuge model. ... This limitation was the requirement to use an a priori model. The exact solutions in the different examples presented in this chapter were known, and they were used as a priori models into the inversion process. However, in practice the exact solution is unknown, and the aim of the tomographic inversion is to obtain a solution that best describes the measured data. Carrying out inversion without using an a priori model can yield an output model that hints at the nature of the model. This output can then be used as the starting point in an iterative process, in which the output from one step is used as an a priori model for reinverting the original data in a subsequent step. In this case, this process slightly improved the output tomogram and decreased the value of root mean squares of travel time residuals (Rrms). An alternative inversion strategy was proposed based on the results obtained in this study. It involves using a searching algorithm. A searching process can be carried out based on the output from the first iteration (without using an a priori model). The search can involve varying the parameters that describe buried anomalies, such as the size of the anomaly, the velocity value in the anomaly, and the location of the anomaly. The aim is to search for the combination of anomaly parameters that minimises the resulting error parameters (mainly Rrmx, but also the average error and the standard deviation of the error). For more subtle cases, such as the velocity model under a footing, where inversion without using an a priori model did not recover the input model, a searching algorithm involving applying perturbations to the exact Boussinesq model can be performed. Not only can the searching procedure involve adding perturbation to the velocity values in the Boussinesq model, but it can also add perturbation to the shape of the velocity distribution below the footing. The searching process can continue until a model that fits the data with a minimum error is found, i.e., minimising Rrms.

"July 2004" / Bibliography: leaves 215-230. / viii, 230 leaves : ill., maps ; 30 cm. / Title page, contents and abstract only. The complete thesis in print form is available from the University Library. / Thesis (Ph.D.)--University of Adelaide, School of Chemistry and Physics, Discipline of Physics, 2004

Many scientific studies have been conducted on the Cretaceous-Paleogene boundary (KTB) in the Gulf coast region and, in particular, the Brazos River section in Falls County, Texas. Despite this, there remains much to be learned about the KTB and its depositional environment. Study of the KTB has been multidisciplinary, primarily in the fields of sedimentology and paleontology. Some researchers in these disciplines have questioned the consensus view of the placement of the KTB and subsequent interpretation of the timing of depositional events and mass extinction events. Geophysical methods have potential to provide additional understanding of the physical properties of the KTB. To date, study of the KTB has relied on point data and borehole information to create cross sections of the complex. Seismic refraction surveys can provide spatially continuous information on susburface horizons located adjacent to the KTB. In this study, seismic first-arrival traveltimes are processed with a tomographic modeling program to map the top of the hummocky cross-bedded sandstone (HCS), which is a key indicator of the deposition environment at the time of KTB boundary complex placement. The survey area is located at Cottonmouth Creek, a tributary of the Brazos River. Three seismic lines were surveyed, one across Cottonmouth Creek, and two parallel to the creek on either side. The data from the two parallel lines were processed using the 2-D seismic refraction tomography algorithm of Zelt and Smith. The reconstructed depth to the HCS in the survey area is approximately 6 m, with layer seismic velocities of 364, 1800, and 2200 m/s, respectively. Seismic tomography successfully mapped the HCS layer and reveals approximately 1 m amplitude undulations vertically and undulations on the order of several m horizontally. These variations are consistent with exposed surfaces of the HCS in the creek bed. Seismic refraction has been utilized successfully herein to map a key buried indicator, namely the top of the HCS layer, associated with the KTB complex. A detailed 3-D seismic refraction survey at this site is recommended to generate a high-resolution 2-D terrain map of the top of the HCS layer.

No description available.

This dissertation explores subduction zone-related deformation both on short time scales in the form of subduction zone earthquakes and over larger time and geographical scales in the form of subduction rollback or detachment of the subducting lithosphere. The study presented here is composed of two parts. First, we analyzed the source-rupture processes of the April 1, 2007 Solomon Islands Earthquake (Mw=8.1) using a body-wave inversion technique. Our analysis indicated that the earthquake ruptured approximately 240 km of the southeast Pacific subduction zone in two sub-events.In the second part of this study, we used shear-wave splitting analysis to investigate the effects of the subducting African lithosphere on the upper-mantle flow field beneath the Anatolian Plate in the Eastern Mediterranean region. Our shear-wave splitting results are consistent with relatively uniform southwest-directed flow towards the actively southwestward-retreating Aegean slab. Based on spatial variations in observed delay times we identified varying flow speeds beneath Anatolia and we attribute this variation to the differential retreat rates of the Aegean and the Cyprean trenches.Finally, we used teleseismic P-wave travel-time tomography to image the geometry of the subducting African lithosphere beneath the Anatolia region. Our tomograms show that the subducting African lithosphere is partitioned into at least two segments along the Cyprean and the Aegean trenches. We observed a gap between the two segments through which hot asthenosphere ascends beneath the volcanic fields of western Anatolia. Our results show that the Cyprean slab is steeper than the Aegean slab. We inferred that this steep geometry, in part, controls the flow regime of asthenosphere beneath Anatolia causing variations in flow speeds inferred from shear-wave splitting analysis.