Geophysics and Seismology Commons^™

Clustering Revisited: A Spectral Analysis Of Microseismic Events, Deborah Fagan, Kasper Van Wijk, James Rutledge

CGISS Publications and Presentations

Identifying individual subsurface faults in a larger fault system is important to characterize and understand the relationship between microseismicity and subsurface processes. This information can potentially help drive reservoir management and mitigate the risks of natural or induced seismicity. We have evaluated a method of statistically clustering power spectra from microseismic events associated with an enhanced oil recovery operation in southeast Utah. Specifically, we were able to provide a clear distinction within a set of events originally designated in the time domain as a single cluster and to identify evidence of en echelon faulting. Subtle time-domain differences between events were …

Recognizing And Modeling Variable Drawdown Due To Evapotranspiration In A Semiarid Riparian Zone Considering Local Differences In Vegetation And Distance From A River Source, Brady Johnson, Bwalya Malama, Warren Barrash, Alejandro N. Flores

CGISS Publications and Presentations

Riparian zones in semiarid regions often exhibit high rates of evapotranspiration (ET) in spite of low-soil moisture content due to the presence of phreatophytic vegetation that is able to withdraw water from shallow aquifers. This work seeks to better define the relationship between ET, the saturated zone and the river boundary by comparing observed water table drawdown records to analytically modeled drawdown in fully penetrating wells of an unconfined aquifer in response to daily ET flux. ET at the Boise Hydrogeophysical Research Site (BHRS), a riparian zone in a temperate, semiarid environment, is calculated using a radiation-based method to provide …

Compressional And Em Wave Velocity Anisotropy In A Temperate Glacier Due To Basal Crevasses, And Implications For Water Content Estimation, John H. Bradford, Joshua Nichols, Joel T. Harper, Toby Meierbachtol

CGISS Publications and Presentations

We have conducted a series of experiments designed to investigate elastic and electromagnetic (EM) velocity anisotropy associated with a preferentially aligned fracture system on a temperate valley glacier in south-central Alaska, USA. Measurements include a three-dimensional compressional wave (P-wave) seismic reflection survey conducted over a 300m × 300m survey patch, with uniform source grid and static checkerboard receiver pattern. Additionally, we acquired a multi-azimuth, multi-offset, polarimetric ground-penetrating radar (GPR) reflection experiment in a wagon-wheel geometry with 94° of azimuthal coverage. Results show azimuthal variation in the P-wave normal-moveout velocity of >3% (3765 and 3630ms^-1 in the fast and slow …

Megathrust Splay Fault Geometry In Prince William Sound, Alaska, Shaun Patrick Finn

Boise State University Theses and Dissertations

I present marine seismic reflection results from Prince William Sound, Alaska that document the location of active faults related to the subduction zone processes.

Subduction zones along convergent margins experience large earthquakes, magnitude >8, with recurrence intervals on the order of centuries. Smaller magnitude earthquakes with shorter recurrence intervals are probable along the same subduction zone fault zones. Convergent margin earthquakes also are associated with high uplift rates and tsunami generation, yet the location and uplift history of most tsunamigenic faults are unknown. In this thesis, I present the processed results of high resolution marine seismic reflection data of Prince …

Variability Of Pennsylvanian-Permian Carbonate Associations And Implications For Nw Pangea Palaeogeography, East-Central British Columbia, Canada, K. D. Zubin-Stathopoulos, B. Beauchamp, V. I. Davydov, C. M. Henderson

Geosciences Faculty Publications and Presentations

Different stages of Pennsylvanian-Permian carbonate sedimentation in east-central British Columbia record a complex history of changing environments influenced by evolving palaeogeography and climate. Newly recognized tectonically controlled features affected the distribution and variability of carbonate associations, providing new interpretations for this portion of the west coast of Pangea. Both a heterozoan (cool water) and photozoan (warm-water) association were identified on either side of a palaeogeographic high here informally termed “Tipinahokan Peninsula”. Cool water carbonates were located outboard, or to the west of this high, an area influenced by upwelling waters. Inboard of this high, a warm, protected sea developed, here …

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 …

Extension Of The Spatial Autocorrelation (Spac) Method To Mixed-Component Correlations Of Surface Waves, Matthew M. Haney, T. Dylan Mikesell, Kasper Van Wijk, Hisashi Nakahara

CGISS Publications and Presentations

Using ambient seismic noise for imaging subsurface structure dates back to the development of the spatial autocorrelation (SPAC) method in the 1950s. We present a theoretical analysis of the SPAC method for multicomponent recordings of surface waves to determine the complete 3 × 3 matrix of correlations between all pairs of three-component motions, called the correlation matrix. In the case of isotropic incidence, when either Rayleigh or Love waves arrive from all directions with equal power, the only non-zero off-diagonal terms in the matrix are the vertical–radial (ZR) and radial–vertical (RZ) correlations in the presence of Rayleigh waves. Such combinations …

Northern Thailand Geothermal Resources And Development: A Review And 2012 Update, Fongsaward Suvagondha Singharajwarapan, Spencer H. Wood, Natthaporn Prommakorn, Lara Owens

Geosciences Faculty Publications and Presentations

Northern Thailand has 16 hot spring systems with surface temperatures near or greater than 80°C with potential for binary plant power generation. Presently only Fang system generates power from wells flowing a total of 8.3 1/s of 116°C water to a 300 kW single module Ormat binary plant. Current production is 150-250 KW, which potentially can be increased by constructing new wells and increasing flow by pumping. Of the other 15 systems, 4 are in national parks and not considered for development. Several of the hot springs systems have silica geothermometry >130°C suggesting significant undeveloped resources exist in northern Thailand. …

Probing Local Wind And Temperature Structure Using Infrasound From Volcan Villarrica (Chile), Jeffrey B. Johnson, J. Anderson, O. Marcillo, S. Arrowsmith

Geosciences Faculty Publications and Presentations

We use the continuous and intense (∼10⁷ W) infrasound produced by Volcan Villarrica (Chile) to invert for the local dynamic wind and temperature structure of the atmosphere. Infrasound arrays deployed in March 2011 at the summit (2826 m) and on the NNW flank (∼8 km distant at 825 m) were used to track infrasound propagation times and signal power. We model an atmosphere with vertically varying temperature and horizontal winds and use propagation times (ranging from 23 to 24 s) to invert for horizontal slowness (2.75–2.94 s/km) and average effective sound speeds (328–346 m/s) for NNW propagating infrasound. The …

Hydrologic Data Assimilation With A Hillslope-Scale-Resolving Model And L Band Radar Observations: Synthetic Experiments With The Ensemble Kalman Filter, Alejandro N. Flores, Rafael L. Bras, Dara Entekhabi

Geosciences Faculty Publications and Presentations

Soil moisture information is critical for applications like landslide susceptibility analysis and military trafficability assessment. Existing technologies cannot observe soil moisture at spatial scales of hillslopes (e.g., 10⁰ to 10² m) and over large areas (e.g., 10² to 10⁵ km²) with sufficiently high temporal coverage (e.g., days). Physics-based hydrologic models can simulate soil moisture at the necessary spatial and temporal scales, albeit with error. We develop and test a data assimilation framework based on the ensemble Kalman filter for constraining uncertain simulated high-resolution soil moisture fields to anticipated remote sensing products, specifically NASA's Soil …

Statistical Clustering Of Microseismic Event Spectra To Identify Subsurface Structure, Deborah Kay Fagan

Boise State University Theses and Dissertations

Understanding subsurface structure by studying microseismicity influences a wide range of activities, including energy extraction, aquifer storage, carbon sequestration, and seismic hazard assessment. Identifying individual fractures in a larger fault system is key to characterizing, understanding, and potentially mitigating risks of natural or induced seismicity.

A year-long study associated with a carbon dioxide (CO₂) sequestration project was conducted at the Aneth oil field in southeast Utah to record microseismicity at a single downhole geophone array. A previous analysis located events by first identifying event multiplets consisting of highly correlated time-domain waveforms on receivers shallower than the depth of …

Monitoring Co2 Sequestration In Basalt With Elastic Waves, Larry Thomas Otheim

Boise State University Theses and Dissertations

The continued burning of fossil fuels as a source of energy is contributing to greater concentrations of carbon dioxide (CO₂) in the atmosphere. Increased levels of CO₂ in the atmosphere have been linked to an increase of global mean surface temperature. To mitigate the continued release of CO₂, projects to capture this gas at large point sources and sequester it in geologic formations are in place. Carbon dioxide sequestration in basalts promises permanent trapping of the fluid as these rocks react with carbonic acid and precipitate carbonate minerals. It is important to monitor the injection …

A Feasibility Study Of Time-Lapse Seismic Monitoring Of Co2 Sequestration In A Layered Basalt Reservoir, Murari Khatiwada, Ludmila Adam, Michael Morrison, Kasper Van Wijk

CGISS Publications and Presentations

We investigate the potential of scattered seismic waves to remotely sense geological sequestration of CO₂ in basalt. Numerical studies in horizontally layered models suggest that strong scattering quickly complicates the wave fields, but also provides a sensitive tool to monitor physical changes in and around the reservoir. These results go hand-in-hand with recent laboratory work and rock-physics modeling that has shown significant changes in the seismic properties of a reservoir undergoing CO₂ sequestration, due to fluid substitution and mineral precipitation.

Scanning For Velocity Anomalies In The Crust And Mantle With Diffractions From The Core-Mantle Boundary, Elmer Ruigrok, T. Dylan Mikesell, Kasper Van Wijk

CGISS Publications and Presentations

A novel method, based on differential arrival times of diffractions from the core-mantle boundary, swiftly scans for seismic velocity anomalies in the crust and mantle below an array of seismometers. The method is applied to data from the USArray and the large-scale structural features in the western United States are resolved. High lateral resolution is achieved, but structure is averaged over depth. As such, this method is complementary to surface-wave and tomographic body-wave methods, where averaging takes place in the lateral sense. Processing and data-volume requirements involved are minimal. Therefore, this method can be applied during the early stages of …

The Emergence And Future Of Near-Surface Geophysics, William E. Doll, Richard D. Miller, John Bradford

CGISS Publications and Presentations

Over the past 30 years, geophysical methods have assumed a much more prominent and integral role in many investigations where subsurface features have environmental and engineering importance. In fact, the field once referred to as "environmental and engineering geophysics" has broadened to include other applications (e.g., archeology, forensics), and is now commonly referred to more generally as "near-surface geophysics." It is difficult to precisely define near-surface geophysics, and the definition will likely depend on whom you ask. However, we define it as the use of geophysical methods to investigate the zone between the surface and hundreds of meters into the …

Three-Dimensional Stochastic Estimation Of Porosity Distribution: Benefits Of Using Ground-Penetrating Radar Velocity Tomograms In Simulated-Annealing-Based Or Bayesian Sequential Simulation Approaches, Baptiste Dafflon, Warren Barrash

CGISS Publications and Presentations

Estimation of the three-dimensional (3-D) distribution of hydrologic properties and related uncertainty is a key for improved predictions of hydrologic processes in the subsurface. However it is difficult to gain high-quality and high-density hydrologic information from the subsurface. In this regard a promising strategy is to use high-resolution geophysical data (that are relatively sensitive to variations of a hydrologic parameter of interest) to supplement direct hydrologic information from measurements in wells (e.g., logs, vertical profiles) and then generate stochastic simulations of the distribution of the hydrologic property conditioned on the hydrologic and geophysical data. In this study we develop and …

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 …

A Field Proof-Of-Concept Of Aquifer Imaging Using 3-D Transient Hydraulic Tomography With Modular, Temporarily-Emplaced Equipment, Michael Cardiff, Warren Barrash, Peter K. Kitanidis

CGISS Publications and Presentations

Hydraulic tomography is a field scale aquifer characterization method capable of estimating 3-D heterogeneous parameter distributions, and is directly sensitive to hydraulic conductivity (K), thus providing a useful data source for improving flow and transport models. We present results from a proof-of-concept field and modeling study in which we apply 3-D transient hydraulic tomography (3DTHT) to the relatively high-K and moderately heterogeneous unconfined aquifer at the Boise Hydrogeophysical Research Site. Short-duration (20 min) partially penetrating pumping tests, for which observed responses do not reach steady state, are used as the aquifer stimulation. To collect field data, we utilize a system …

Subsurface Characterization Using Head-Wave Artifacts In Seismic Interferometry, Thomas Dylan Mikesell

Boise State University Theses and Dissertations

Seismologists continually work to improve images of the Earth's interior. One new approach is seismic interferometry, which involves cross-correlating the seismic wave field recorded at two receivers to generate data as if one of the receivers was a source. Over the past decade, seismic interferometry has become an established technique to estimate the surface-wave part of the impulse response between two receivers; however, practical limitations in the source-energy distribution have made body-wave recovery difficult and causes spurious energy in the estimated impulse response. Rather than suppress such spurious energy, it can be useful to analyze coherent spurious events to help …

Non-Destructively Mapping The In-Situ Hydrologic Properties Of Snow, Firn, And Glacial Ice With Georadar, Joel Matthew Brown

Boise State University Theses and Dissertations

Ground penetrating radar (GPR) is a useful tool for studying the in-situ properties of glacial ice, firn, and snowpacks. The main focus of this dissertation is improving and expanding methods employed when collecting, processing, and understanding GPR data collected in the Cryosphere, or the snow and ice covered areas of the earth. The data used herein were collected on the Greenland Ice Sheet (GrIS) and on seasonal snowpacks of Western Montana, USA. This document has three sub-topics.

The first sub-topic is comparing the spatial variability of GPR data to the spatial variability of core data collected in two locations within …

Hillslope Asymmetry Maps Reveal Widespread, Multi-Scale Organization, Michael J. Poulos, Jennifer L. Pierce, Alejandro N. Flores, Shawn G. Benner

Geosciences Faculty Publications and Presentations

Hillslope asymmetry is the condition in which oppositely-facing hillslopes within an area have differing average slope angles, and indicates aspect-related variability in hillslope evolution. As such, the presence, orientation and magnitude of asymmetry may be a useful diagnostic for understanding process dominance. We present a new method for quantifying and mapping the spatial distribution of hillslope asymmetry across large areas. Resulting maps for the American Cordillera of the Western Hemisphere and the western United States reveal that hillslope asymmetry is widespread, with distinct trends at continental to drainage scales. Spatial patterns of asymmetry correlate with latitude along the American Cordillera, …

Georadar-Derived Estimates Of Firn Density In The Percolation Zone, Western Greenland Ice Sheet, Joel Brown, John Bradford, Joel Harper, W. Tad Pfeffer, Neil Humphrey, Ellen Mosley-Thompson

CGISS Publications and Presentations

Greater understanding of variations in firn densification is needed to distinguish between dynamic and melt-driven elevation changes on the Greenland ice sheet. This is especially true in Greenland’s percolation zone, where firn density profiles are poorly documented because few ice cores are extracted in regions with surface melt. We used georadar to investigate firn density variations with depth along a ~70 km transect through a portion of the accumulation area in western Greenland that partially melts. We estimated electromagnetic wave velocity by inverting reflection traveltimes picked from common midpoint gathers. We followed a procedure designed to find the simplest velocity …

Sequence Stratigraphy And Onlap History Of The Donets Basin, Ukraine: Insight Into Carboniferous Icehouse Dynamics, J. M. Eros, I. P. Montañez, D. A. Osleger, Vladimir Davydov, T. I. Nemyrovska, V. I. Poletaev, M. V. Zhykalyak

Geosciences Faculty Publications and Presentations

The degree to which Permo-Carboniferous cyclothemic successions archive evidence for long-term variations in ice volume during the Late Paleozoic Ice Age is insufficiently resolved. Here we develop the sequence stratigraphy and onlap-offlap history for a 33-my interval of the Carboniferous using the U-Pb calibrated succession of the Donets Basin, Ukraine, in order to assess the relationship between sea-level, high-latitude changes in glacial extent, and climate. Integrated subsurface and outcrop data permit meter-scale correlation of 242 biostratigraphically constrained limestones and coals, and in turn individual cyclothems, across ~250 km of the Donets Basin. Rapid uniform subsidence and basinwide continuity of marker …

Partitioning Of Rare Earth And High Field Strength Elements Between Titanite And Phonolitic Liquid, Paul H. Olin, John A. Wolff

Geosciences Faculty Publications and Presentations

We present the results of a LA–ICPMS study of titanites and associated glasses from the mixed-magma phonolitic Fasnia Member of the Diego Hernández Formation, Tenerife, Canary Islands. We employ a method of identifying equilibrium mineral–melt pairs from natural samples using REE contents and a linear form of the lattice strain model equation (Blundy and Wood, 1994), where the Young's modulus (E_M) for the 7-fold coordinated site is an output variable. For felsic magmas that contain crystals potentially derived from a variety of environments within the system, this approach is more rigorous than the use of solely textural criteria …

Near-Surface Imaging Of A Hydrogeothermal System At Mount Princeton, Colorado Using 3d Seismic, Self-Potential, And Dc Resistivity Data, Andrew P. Lamb, Lee M. Liberty, Kasper Van Wijk, André Revil, Chuck Diggins

CGISS Publications and Presentations

The Upper Arkansas Valley in the Rocky Mountains of central Colorado is the northernmost extensional basin of the Rio Grande Rift(Figure 1). The valley is a half graben bordered to the east and west by the Mosquito and Sawatch ranges,respectively. The Sawatch Range is home to the Collegiate Peaks,which include some of the highest summits in the Rocky Mountains. Some Collegiate Peaks over 4250 m (14,000 ft) from north to south include Mount Harvard, Mount Yale, Mount Princeton,and Mount Antero. The Sawatch range-front normal fault strikes north-northwest along the eastern margin of the Collegiate Peaks and is characterized by a …

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 …

Hydrological Parameter Estimations From A Conservative Tracer Test With Variable-Density Effects At The Boise Hydrogeophysical Research Site, B. Dafflon, Warren Barrash, Michael Cardiff, T. C. Johnson

CGISS Publications and Presentations

Reliable predictions of groundwater flow and solute transport require an estimation of the detailed distribution of the parameters (e.g., hydraulic conductivity, effective porosity) controlling these processes. However, such parameters are difficult to estimate because of the inaccessibility and complexity of the subsurface. In this regard, developments in parameter estimation techniques and investigations of field experiments are still challenging and necessary to improve our understanding and the prediction of hydrological processes. Here we analyze a conservative tracer test conducted at the Boise Hydrogeophysical Research Site in 2001 in a heterogeneous unconfined fluvial aquifer. Some relevant characteristics of this test include: variable-density …

Co2 Sequestration In Basalt: Carbonate Mineralization And Fluid Substitution, Thomas L. Otheim, Ludmila Adam, Kasper Van Wijk, Michael L. Batzle, Travis Mcling, Robert Podgorney

CGISS Publications and Presentations

Geologicalsequestration of carbon dioxide in deep reservoirs may provide alarge-scale option for reducing the emissions of this gas intothe atmosphere. The effectiveness of sequestration depends on the storagecapacity and stability of the reservoir and risk of leakageinto the overburden. Reservoir rocks can react with a CO₂-watermixture, potentially resulting in the precipitation of minerals in theavailable matrix pore space and within pre-existing fractures. This inducedmineralization may form internal seals that could help mitigate theleakage of CO₂ into the overburden. For basaltic host rocks,carbonic acid partially dissolves minerals in the host rock, suchas the calcium plagioclase mineral, freeing various …

A Simplified Approach For Estimating Soil Carbon And Nitrogen Stocks In Semi-Arid Complex Terrain, Melvin L. Kunkel, Alejandro N. Flores, Toni J. Smith, James P. Mcnamara, Shawn G. Benner

Geosciences Faculty Publications and Presentations

We investigated soil carbon (C) and nitrogen (N) distribution and developed a model, using readily available geospatial data, to predict that distribution across a mountainous, semi-arid, watershed in southwestern Idaho (USA). Soil core samples were collected and analyzed from 133 locations at 6 depths (n=798), revealing that aspect dramatically influences the distribution of C and N, with north-facing slopes exhibiting up to 5 times more C and N than adjacent southfacing aspects. These differences are superimposed upon an elevation (precipitation) gradient, with soil C and N contents increasing by nearly a factor of 10 from the bottom (1100 m elevation) …

Theory And Laboratory Experiments Of Elastic Wave Scattering By Dry Planar Fractures, Thomas E. Blum, Roel Snieder, Kasper Van Wijk, Mark E. Willis

CGISS Publications and Presentations

Remote sensing of fractures with elastic waves is important in fields ranging from seismology to nondestructive testing. In many geophysical applications, fractures control the flow of fluids such as water, hydrocarbons or magma. While previous analytic descriptions of scattering mostly deal with very large or very small fractures (compared to the dominant wavelength), we present an analytic solution for the scattering of elastic waves from a fracture of arbitrary size. Based on the linear slip model for a dry fracture, we derive the scattering amplitude in the frequency domain under the Born approximation for all combinations of incident and scattered …