Geophysics and Seismology Commons^™

Using Shallow Electromagnetic And Magnetic Techniques To Map Soil Texture And Characterize Salinity And Water Dynamics Below Pecan Orchards, El Paso County, Texas, Kristina Sasser

Open Access Theses & Dissertations

Non-invasive, near surface geophysical tools provide a time efficient and cost-effective approach to study the shallow subsurface. Electromagnetic induction (EMI) instruments are a category of these near surface methods that provide spatial and temporal information (2-D to 4-D) about shallow (<6 m) subsurface properties, from which soil salinity, clay content, and water content can be estimated. However, deconstructing soil apparent electrical conductivity (ECa) from EMI methods into its component parts remains a challenge. This research uses EMI and magnetic geophysical techniques to: (1) compare and contrast the distribution and behavior of ECa, both spatially and temporally, at two flood irrigated agricultural sites (pecan orchards) with different soil layers, properties, and controls on electrical conductivities that lie within the same floodplain in far west Texas; and (2) develop a one-dimensional inversion model using ECa and soil texture data at specified layers from established sites of known high and low ECa to determine soil layer resistivities at various dates during the irrigation season. Data are compared to soil sensor bulk EC and volumetric water content data at corresponding depths to validate results. Soils at both sites exhibit large ECa contributions from textural controls, and irrigation frequency. The combination of these two controls dictate where and how quickly secondary minerals precipitate, clog soil pores, reduce infiltration, and salinize the groundwater. Insight from this research aids in ongoing efforts to characterize vadose zone hydrology in dryland critical zones with high infiltration and serves as a guide for locations where future soil sampling will occur.

Approaches To Model Non-Uniqueness And Site Complexity For Non-Invasive Shear-Wave Depth Profiling, Christopher Boucher

Electronic Thesis and Dissertation Repository

Shear-wave velocity (VS) depth profiling and associated seismic site classifications were performed at 15 sites across Metro Vancouver, British Columbia using passive seismic and surface wave methods. Inversion model parameters are constrained at each site using nearby geodata in combination with developing regression models of shear-wave velocity with depth for three primary stratigraphic units. Statistical methods such as a Bayesian Information Criterion are applied post-inversion to evaluate models between and within varying parameterizations. Data evaluation metrics, including the use of microtremor horizontal-to-vertical spectral ratios (MHVSRs), are applied to identify two common deviations from the simple case of normally dispersive laterally …

Changes In Crater Geometry Revealed By Inversion Of Harmonic Infrasound Observations: 24 December 2018 Eruption Of Mount Etna, Italy, Leighton M. Watson, Jeffrey B. Johnson, Mariangela Sciotto, Andrea Cannata

Geosciences Faculty Publications and Presentations

In December 2018, Mount Etna (Italy) experienced a period of increased eruptive activity that culminated in a fissure eruption on the southeast flank. After the onset of the flank eruption, the peak frequency of the summit infrasound signals decreased while resonance increased. We invert infrasound observations for crater geometry and show that crater depth and radius increased during the eruption, which suggests that the flank eruption drained magma from the summit and that eruptive activity led to erosion of the crater wall. By inverting the entire infrasound amplitude spectra rather than just the peak frequency, we are able to place …

Simulation Of A Salt Dome Using 2d Linear And Nonlinear Inverse Modeling Of Residual Gravity Field Data, Soheyl Pourreza, Farnush Hajizadeh

Bulletin of the Mineral Research and Exploration

In gravity field inversion we usually dealing with underdetermined problems and for obtaining realistic solutions can introduce a depth-weighting function to the inversion algorithm. We employ a linear inversion method for determining the underground density distribution of the gravity causative mass. The validation and accuracy of method is tested on two synthetic gravity anomaly from different models, while the data are noise- free and corrupted with noise. In this paper, We also invert the 2D gravity anomaly produced by a salt dome from the northwest of Iran. The salt domes in the region under investigation are a rich source of …

Crustal Composition Beneath Southern Idaho: Insights From Teleseismic Receiver Functions, Thomas Branson Harper

Boise State University Theses and Dissertations

Receiver functions derived from teleseismic earthquakes contain seismic amplitude and velocity information that relate to compositional changes within the Earth’s crust and upper mantle. The receiver function waveform is a combination of P-S converted waves that have reverberated within the lithosphere. Although the largest seismic velocity boundary is found at the base of the crust, I explore the use of lower amplitude receiver function arrivals that represent smaller velocity contrasts within the crust. In my thesis, I calculate and model receiver functions via a Metropolis algorithm approach to extract seismic velocity distributions in the lithosphere. I use the results to …

Crustal Seismic Anisotropy Of The Ruby Mountains Core Complex And Surrounding Northern Basin And Range, Justin T. Wilgus

Earth and Planetary Sciences ETDs

Metamorphic core complexes (MCC) are distinctive uplifts that expose deeply exhumed and deformed crustal rocks due to localized extensional deformation. Consequently, their detailed structure provide a window into deep crustal mechanics. The North American Cordillera contains numerous MCC, one of which is the Ruby Mountains core complex (RMCC) located in the highly extended northern Basin and Range. To constrain the extent to which anisotropy below the RMCC deviates from the regional Basin and Range average and test the depth dependence of crustal anisotropy we conduct a radial anisotropy investigation below the RMCC and surrounding northern Basin and Range. Data from …

3d Seismic Attribute Analysis And Machine Learning For Reservoir Characterization In Taranaki Basin, New Zealand, Aamer Alhakeem

Doctoral Dissertations

"The Kapuni group within the Taranaki Basin in New Zealand is a potential petroleum reservoir. The objective of the study includes building a sequential approach to identify different geological features and facies sequences within the strata, through visualizing the targeted formations by interpreting and correlating the regional geological data, 3D seismic, and well data by following a sequential workflow. First, seismic interpretation is performed targeting the Kapuni group formations, mainly, the Mangahewa C-sand and Kaimiro D-sand. Synthetic seismograms and well ties are conducted for structural maps, horizon slices, isopach, and velocity maps. Well log and morphological analyses are performed for …

Imaging Wetland Hydrogeophysics: Applications Of Critical Zone Hydrogeophysics To Better Understand Hydrogeologic Conditions In Coastal And Inland Wetlands And Waters, Christine Marie Downs

USF Tampa Graduate Theses and Dissertations

This dissertation consists of three projects utilizing electric and electromagnetic (EM) methods to better understand critical-zone hydrogeologic conditions in select Florida wetlands and waters.

First, a time-lapse electrical resistivity (ER) survey was conducted in section of mangrove forest on a barrier island in southeast Florida to image changes in pore-water salinity in the root zone. ER data show the most variability in the root zone over a 24-hour period, and, generally, the ground is more resistive during the day than overnight.

Second, a suite of three-dimensional forward models, based on varying lateral boundaries and conductivities typical of a coastal wetland, …

Detecting Overpressure Zones By Using Model Based Inversion In Kupe Field, New Zealand, Timucin Cakir

Dissertations, Master's Theses and Master's Reports

Overpressure is a significant parameter during a drilling process in oil exploration. This phenomenon happens when pore pressure exceeds the limit of hydrostatic pressure in a formation and causes drastic hazards such as blow outs while drilling. This thesis focuses on the determination of possible overpressure areas and studies at three wells in Kupe Field, which is considered to be a natural-gas field in Taranaki Basin. In an effort to identify potential overpressure zones at the Kupe Field, acoustic impedance responses of 3D seismic data for three wells, along with their well log responses were evaluated.

Seismic inversion plays an …

Improved 2d And 3d Resistivity Surveys Using Buried Electrodes And Optimized Arrays: The Multi-Electrode Resistivity Implant Technique (Merit), Henok Gidey Kiflu

USF Tampa Graduate Theses and Dissertations

This thesis presents a novel resistivity method called Multi-Electrode resistivity technique (MERIT) that is used for high resolution imaging of complex geologic features at depth and near the edges of survey lines. The MERIT electrodes are especially shaped and designed to be self-driven using a robust-direct push technique. Measurements are taken using optimized arrays that are generated using a modified version of the “Compare-R” optimization algorithm. This work focused on both two-dimensional (MERIT2D) and three-dimensional (MERIT3D) applications of the buried array and show the relevance of the additional information gained by the addition of deep electrodes especially in sites with …

Reservoir Estimation In The Penobscot 3d Seismic Volume Using Constrained Sparse Spike Inversion, Offshore Nova Scotia, Canada, Mengchu Xiao

Dissertations, Master's Theses and Master's Reports

The Penobscot study area is located offshore Nova Scotia, Canada. There are two wells, which penetrate the highest potentially commercial bodies in the Abenaki Formation. In order to investigate the potential for locating additional hydrocarbon reservoirs, well log data was used and the Penobscot 3D seismic dataset was analyzed using Constrained Sparse Spike Inversion. From the well log data, low GR and SP values are an indication of a permeable sand layer, which provides the target zone in this study. Impedance – porosity crossplots gave the relationship between impedance and porosity, where a low impedance sand layer is correlated with …

Targeted Reflection-Waveform Inversion Of Experimental Ground-Penetrating Radar Data For Quantification Of Oil Spills Under Sea Ice, John H. Bradford, Esther L. Babcock, Hans-Peter Marshall, David F. Dickins

Geosciences Faculty Publications and Presentations

Rapid spill detection and mapping are needed with increasing levels of oil exploration and production in the Arctic. Previous work has found that ground-penetrating radar (GPR) is effective for qualitative identification of oil spills under, and encapsulated within, sea ice. Quantifying the spill distribution will aid effective spill response. To this end, we have developed a targeted GPR reflection-waveform inversion algorithm to quantify the geometry of oil spills under and within sea ice. With known electric properties of the ice and oil, we have inverted for oil thickness and variations in ice thickness. We have tested the algorithm with data …

Reservoir Characterization Through Pre-Stack Seismic Analysis And Inversion In The Thrace Basin, Northwest Turkey, Emre Doguturk

Dissertations, Master's Theses and Master's Reports

This study analyzes hydrocarbon reservoirs by using Amplitude-Versus-Offset (AVO) analysis and different inversion methods to investigate hydrocarbon reservoir in the Thrace Basin which is located in northwest Turkey.

A 3D seismic survey containing prestack data was provided by Turkish Petroleum Corporation as true amplitude NMO corrected 3D gathers.

The quality of this land data was poor for AVO applications, and steps were taken to make it useful. Radon filtering proved to be extremely useful for this purpose. Although application of this filter affected the AVO characteristics differently depending on the parameters selected for the filter, relative AVO characteristics remained useful …

Reflection Waveform Inversion Of Ground-Penetrating Radar Data For Characterizing Thin And Ultrathin Layers Of Nonaqueous Phase Liquid Contaminants In Stratified Media, Esther Babcock, John H. Bradford

CGISS Publications and Presentations

Accurately quantifying thin-layer parameters by applying a targeted reflection waveform inversion methodology to ground-penetrating radar (GPR) reflection data may provide a useful tool for near-surface investigation and especially for contaminated site investigation where nonaqueous phase liquid (NAPL) contaminants are present. We implemented a targeted reflection waveform inversion algorithm to quantify thin-layer permittivity, thickness, and conductivity for NAPL thin (≤ 1/2 dominant wavelength λ) and ultrathin (≤ 1/8λ) layers using GPR reflection data. The inversion used a nonlinear grid search with a Monte Carlo scheme to initialize starting values to find the global minimum. By taking a targeted approach using a …

Geophysical Analysis And Modelling Of The S‹Mav Basin, Western Anatolia, Ceyhan Ertan Toker

Bulletin of the Mineral Research and Exploration

The various data processing techniques, to illuminate the parameters of the geological structure which are applied in gravity and magnetic potential field methods. Also edge detection procedures are in data processing techniques. In this study, 2D, 3D, inversion and asta technique is one of the new edge detection procedures were applied to clarify correlation between the Simav half graben’s deep position and geometry of the tectonic lineaments. The Asta is obtained using the tilt angle.

Geophysical Model Of The Cu-Mo Porphyry Ore Deposit At Copper Flat Mine, Hillsboro, Sierra County, New Mexico, Adrian Emmanuel Gutierrez

Open Access Theses & Dissertations

A 3D gravity model of the Copper Flat Mine was performed as part of the exploration of new resources in at the mine. The project is located in the Las Animas Mining District in Sierra County, New Mexico. The mine has been producing ore since 1877 and is currently owned by the New Mexico Copper Corporation, which plans o bringing the closed copper mine back into production with innovation and a sustainable approach to mining development. The Project is located on the Eastern side of the Arizona-Sonora-New Mexico porphyry copper Belt of Cretaceous age. Copper Flat is predominantly a Cretaceous …

Hydraulic Conductivity Imaging From 3-D Transient Hydraulic Tomography At Several Pumping/Observation Densities, Michael Cardiff, Warren Barrash, Peter K. Kitanidis

CGISS Publications and Presentations

3-D Hydraulic tomography (3-D HT) is a method for aquifer characterization whereby the 3-D spatial distribution of aquifer flow parameters (primarily hydraulic conductivity, K) is estimated by joint inversion of head change data from multiple partially penetrating pumping tests. While performance of 3-D HT has been studied extensively in numerical experiments, few field studies have demonstrated the real-world performance of 3-D HT. Here we report on a 3-D transient hydraulic tomography (3-D THT) field experiment at the Boise Hydrogeophysical Research Site which is different from prior approaches in that it represents a “baseline” analysis of 3-D THT performance using …

A Modified Delay-Time Method For Statics Estimation With The Virtual Refraction, T. Dylan Mikesell, Kasper Van Wijk, Elmer Ruigrok, Andrew Lamb, Thomas E. Blum

CGISS Publications and Presentations

Topography and near-surface heterogeneities lead to traveltime perturbations in surface land-seismic experiments. Usually, these perturbations are estimated and removed prior to further processing of the data. A common technique to estimate these perturbations is the delay-time method. We have developed the “modified delay-time method,” wherein we isolate the arrival times of the virtual refraction and estimate receiver-side delay times. The virtual refraction is a spurious arrival found in wavefields estimated by seismic interferometry. The new method removes the source term from the delay-time equation, is more robust in the presence of noise, and extends the lateral aperture compared to the …

Monitoring Glacier Surface Seismicity In Time And Space Using Rayleigh Waves, T. D. Mikesell, Kasper Van Wijk, Matthew M. Haney, John H. Bradford, Hans-Peter Marshall, Joel T. Harper

CGISS Publications and Presentations

Sliding glaciers and brittle ice failure generate seismic body and surface wave energy characteristic to the source mechanism. Here we analyze continuous seismic recordings from an array of nine short-period passive seismometers located on Bench Glacier, Alaska (USA) (61.033°N, 145.687°W). We focus on the arrival-time and amplitude information of the dominant Rayleigh wave phase. Over a 46-hour period we detect thousands of events using a cross-correlation based event identification method. Travel-time inversion of a subset of events (7% of the total) defines an active crevasse, propagating more than 200 meters in three hours. From the Rayleigh wave amplitudes, we estimate …

Infrasound As Upper Atmospheric Monitor, Jelle Daniel Assink

Electronic Theses and Dissertations

Understanding and specification of the higher altitudes of the atmosphere with global coverage over all local times is hampered by the challenges of obtaining direct measurements in the upper atmosphere. Methods to measure the properties of the atmosphere above the stratopause is an active area of scientific research. In this thesis, we revisit the use of infrasound as a passive remote sensing technique for the upper atmosphere. Signals from the Tungurahua volcano in Ecuador are used to investigate the behavior of the upper atmosphere. Depending on the atmospheric conditions, stratospheric, mesospheric and thermospheric arrivals are observed during intervals of explosive …

3-D Transient Hydraulic Tomography In Unconfined Aquifers With Fast Drainage Response, Michael Cardiff, Warren Barrash

CGISS Publications and Presentations

We investigate, through numerical experiments, the viability of three-dimensional transient hydraulic tomography (3DTHT) for identifying the spatial distribution of groundwater flow parameters (primarily, hydraulic conductivity K) in permeable, unconfined aquifers. To invert the large amount of transient data collected from 3DTHT surveys, we utilize an iterative geostatistical inversion strategy in which outer iterations progressively increase the number of data points fitted and inner iterations solve the quasi-linear geostatistical formulas of Kitanidis. In order to base our numerical experiments around realistic scenarios, we utilize pumping rates, geometries, and test lengths similar to those attainable during 3DTHT field campaigns performed at …

Effects Of High-Impedance-Contrast Boundary Upon Multi-Modal Seismic Surface Wave Data, Prajwol Tamrakar

UNLV Theses, Dissertations, Professional Papers, and Capstones

Synthetic studies and analyses of an experimental dataset were conducted to address the use of Rayleigh-type surface waves for estimation of shear wave velocity (VS) profiles of shallow bedrock sites. The shallow bedrock presents a high impedance contrast boundary which causes surface wave energy to be partitioned to higher modes. Idealized studies of hypothetical datasets and root-mean-squared calculations of error surfaces showed that if reliable dispersion data are available over a broad frequency spectrum, the VS profile can be recovered using the fundamental mode alone. However, when dispersion data are limited to a relatively narrow frequency band representing what might …

Software Development For A Three-Dimensional Gravity Inversion And Application To Study Of The Border Ranges Fault System, South-Central Alaska, Rolando Cardenas

Open Access Theses & Dissertations

The Border Ranges Fault System (BRFS) bounds the Cook Inlet and Susitna Basins, an important petroleum province within south-central Alaska. A primary goal in the research is to test several plausible models of structure along the Border Ranges Fault System using a novel three-dimensional inversion utilizing gravity and magnetic data, constrained with other geophysical, borehole and surface geological information. This research involves the development of inversion modeling software using a Borland C++ compiler as part of the Rapid Application Development (RAD) Studio. The novel inversion approach directly models known geology, and "a priori" uncertainties on the geologic model to allow …

Incorporating Geostatistical Constraints In Nonlinear Inverse Problems, T. C. Johnson, Partha S. Routh, Tom Clemo, Warren Barrash, William P. Clement

CGISS Publications and Presentations

In this paper we present a method of incorporating semivariogram constraints into nonlinear inversion problems. That is, we describe a method of sampling the space of inverse solutions that honor a specified semivariogram or set of semivariograms and also explain a set of state data. The approach can be considered a method of conditional simulation where model conditioning is based upon state data (as opposed to parameter data). The difference between this approach and other simulation approaches is that the simulation is posed as an optimization problem with the joint objective of matching the semivariograms and honoring the state data. …