Geophysics and Seismology Commons^™

Hydraulic Tomography: 3d Hydraulic Conductivity, Fracture Network, And Connectivity In Mudstone, Claire R. Tiedeman, Warren Barrash

CGISS Publications and Presentations

We present the first demonstration of hydraulic tomography (HT) to estimate the three‐dimensional (3D) hydraulic conductivity (K) distribution of a fractured aquifer at high‐resolution field scale (HRFS), including the fracture network and connectivity through it. We invert drawdown data collected from packer‐isolated borehole intervals during 42 pumping tests in a wellfield at the former Naval Air Warfare Center, West Trenton, New Jersey, in the Newark Basin. Five additional tests were reserved for a quality check of HT results. We used an equivalent porous medium forward model and geostatistical inversion to estimate 3D K at high resolution (K …

Ionoseis: A Package To Model Coseismic Ionospheric Disturbances, Thomas Dylan Mikesell, Lucie Rolland, Rebekah F. Lee, Florian Zedek, Pierdavide Coïsson, Jean-Xavier Dessa

CGISS Publications and Presentations

We present the framework of the modeling package IonoSeis. This software models Global Navigation Satellite System (GNSS) derived slant total electron content (sTEC) perturbations in the ionosphere due to the interaction of the neutral atmosphere and charged particles in the ionosphere. We use a simplified model to couple the neutral particle momentum into the ionosphere and reconstruct time series of sTEC perturbations that match observed data in both arrival time and perturbation shape. We propagate neutral atmosphere disturbances to ionospheric heights using a three-dimensional ray-tracing code in spherical coordinates called Windy Atmospheric Sonic Propagation (WASP3D), which works for a …

Tsunamigenic Splay Faults Imply A Long‐Term Asperity In Southern Prince William Sound, Alaska, L. M. Liberty, D. S. Brothers, P. J. Haeussler

CGISS Publications and Presentations

Coseismic slip partitioning and uplift over multiple earthquake cycles is critical to understanding upper‐plate fault development. Bathymetric and seismic reflection data from the 1964 Mw9.2 Great Alaska earthquake rupture area reveal sea floor scarps along the tsunamigenic Patton Bay/Cape Cleare/Middleton Island fault system. The faults splay from a megathrust where duplexing and underplating produced rapid exhumation. Trenchward of the duplex region, the faults produce a complex deformation pattern from oblique, south‐directed shortening at the Yakutat‐Pacific plate boundary. Spatial and temporal fault patterns suggest that Holocene megathrust earthquakes had similar relative motions and thus similar tsunami sources as in …

Interannual Snow Accumulation Variability On Glaciers Derived From Repeat, Spatially Extensive Ground-Penetrating Radar Surveys, Daniel Mcgrath, Louis Sass, Shad O'Neel, Chris Mcneil, Salvatore G. Candela, Emily H. Baker, Hans-Peter Marshall

CGISS Publications and Presentations

There is significant uncertainty regarding the spatiotemporal distribution of seasonal snow on glaciers, despite being a fundamental component of glacier mass balance. To address this knowledge gap, we collected repeat, spatially extensive high-frequency ground-penetrating radar (GPR) observations on two glaciers in Alaska during the spring of 5 consecutive years. GPR measurements showed steep snow water equivalent (SWE) elevation gradients at both sites; continental Gulkana Glacier's SWE gradient averaged 115 mm 100 m⁻¹ and maritime Wolverine Glacier's gradient averaged 440 mm 100 m⁻¹ (over > 1000 m). We extrapolated GPR point observations across the glacier surface using terrain parameters derived …

Direct Insertion Of Nasa Airborne Snow Observatory-Derived Snow Depth Time Series Into The Isnobal Energy Balance Snow Model, Andrew R. Hedrick, Hans-Peter Marshall

CGISS Publications and Presentations

Accurately simulating the spatiotemporal distribution of mountain snow water equivalent improves estimates of available meltwater and benefits the water resource management community. In this paper we present the first integration of lidar-derived distributed snow depth data into a physics-based snow model using direct insertion. Over four winter seasons (2013–2016) the National Aeronautics and Space Administration/Jet Propulsion Laboratory (NASA/JPL) Airborne Snow Observatory (ASO) performed near-weekly lidar surveys throughout the snowmelt season to measure snow depth at high resolution over the Tuolumne River Basin above Hetch Hetchy Reservoir in the Sierra Nevada Mountains of California. The modeling component of the ASO program …

Revelation Of Early Detection Of Coseismic Ionospheric Perturbations In Gps-Tec From Realistic Modelling Approach: Case Study, Dhanya Thomas, Mala S. Bagiya, Poikayil Sukumaran Sunil, Lucie Rolland, Anakuzhikkal Sudarsanan Sunil, T. Dylan Mikesell, Srinivas Nayak, Subrahmanyam Mangalampalli, Durbha Sai Ramesh

CGISS Publications and Presentations

GPS-derived Total Electron Content (TEC) is an integrated quantity; hence it is difficult to relate the detection of ionospheric perturbations in TEC to a precise altitude. As TEC is weighted by the maximum ionospheric density, the corresponding altitude (hmF2) is, generally, assumed as the perturbation detection altitude. To investigate the validity of this assumption in detail, we conduct an accurate analysis of the GPS-TEC measured early ionospheric signatures related to the vertical surface displacement of the Mw 7.4 Sanriku-Oki earthquake (Sanriku-Oki Tohoku foreshock). Using 3D acoustic ray tracing model to describe the evolution of the propagating seismo-acoustic wave in space …

Seismic Imaging Of The Main Frontal Thrust In Nepal Reveals A Shallow Décollement And Blind Thrusting, Rafael V. Almeida, Judith Hubbard, Lee Liberty, Anna Foster, Soma Nath Sapkota

CGISS Publications and Presentations

Because great earthquakes in the Himalaya have an average recurrence interval exceeding 500 yr, most of what we know about past earthquakes comes from paleoseismology and tectonic geomorphology studies of the youngest fault system there, the Main Frontal Thrust (MFT). However, these data are sparse relative to fault segmentation and length, and interpretations are often hard to validate in the absence of information about fault geometry. Here, we image the upper two km of strata in the vicinity of the fault tip of the MFT in central Nepal (around the town of Bardibas) applying a pre-stack migration approach to two …

Local, Regional, And Remote Seismo‐Acoustic Observations Of The April 2015 Vei 4 Eruption Of Calbuco Volcano, Chile, Robin S. Matoza, David Fee, David N. Green, Alexis Le Pichon, Julien Vergoz, Matthew M. Haney, T. Dylan Mikesell, Luis Franco, O. Alberto Valderrama, Megan R. Kelley, Kathleen Mckee, Lars Ceranna

CGISS Publications and Presentations

The two major explosive phases of the 22–23 April 2015 eruption of Calbuco volcano, Chile, produced powerful seismicity and infrasound. The eruption was recorded on seismo-acoustic stations out to 1,540 km and on five stations (IS02, IS08, IS09, IS27, and IS49) of the International Monitoring System (IMS) infrasound network at distances from 1,525 to 5,122 km. The remote IMS infrasound stations provide an accurate explosion chronology consistent with the regional and local seismo-acoustic data and with previous studies of lightning and plume observations. We use the IMS network to detect and locate the eruption signals using a brute-force, grid-search, cross-bearings …

Diverse Eruptive Activity Revealed By Acoustic And Electromagnetic Observations Of The 14 July 2013 Intense Vulcanian Eruption Of Tungurahua Volcano, Ecuador, J. F. Anderson, J. B. Johnson, A. L. Steele, M. C. Ruiz, B. D. Brand

CGISS Publications and Presentations

During the powerful July 2013 eruption of Tungurahua volcano, Ecuador, we recorded exceptionally high amplitude, long‐period infrasound (1,600‐Pa peak‐to‐peak amplitude, 5.5‐s period) on sensors within 2km of the vent alongside electromagnetic signals from volcanic lightning serendipitously captured as interference. This explosion was one of Tungurahua's most powerful vulcanian eruptions since recent activity began in 1999, and its acoustic wave is among the most powerful volcanic infrasound ever recorded anywhere. We use these data to quantify erupted volume from the main explosion and to classify postexplosive degassing into distinct emission styles. Additionally, we demonstrate a highly effective method of recording lightning‐related …

Ice Core Records Of West Greenland Melt And Climate Forcing, H.P. Marshall, T. Meehan

CGISS Publications and Presentations

Remote sensing observations and climate models indicate that the Greenland Ice Sheet (GrIS) has been losing mass since the late 1990s, mostly due to enhanced surface melting from rising summer temperatures. However, in situ observational records of GrIS melt rates over recent decades are rare. Here we develop a record of frozen meltwater in the west GrIS percolation zone preserved in seven firn cores. Quantifying ice layer distribution as a melt feature percentage (MFP), we find significant increases in MFP in the southernmost five cores over the past 50 years to unprecedented modern levels (since 1550 CE). Annual to decadal …

Forecasting The Eruption Of An Open‐Vent Volcano Using Resonant Infrasound Tones, Jeffrey B. Johnson, Leighton M. Watson, Jose L. Palma, Eric M. Dunham, Jacob F. Anderson

CGISS Publications and Presentations

Open‐vent volcanic systems with active degassing are particularly effective at producing infrasound that exhibits resonant tones controlled by the geometry of the volcano's crater. Changes in the infrasound character can thus provide constraints on a crater's lava level, which may vary dynamically in the lead‐up to an eruption. Here we show that the increasing frequency content and damping characteristics of the resonant infrasound at Volcán Villarrica (Chile) relate to lava lake position in its crater/conduit preceding its 2015 eruption. We model the acoustic response of Villarrica's crater to determine that the lake began to rise on 27 February and reached …

Geothermal Play Fairway Analysis, Phase 3: A Provisional Conceptual Model Of The Camas Prairie, Snake River Plain, Idaho, Lee M. Liberty

CGISS Publications and Presentations

The Snake River Plain (SRP) Geothermal Play Fairway Analysis team identified two regions of interest during Phase 2 studies: the western SRP near Mountain Home, Idaho and Camas Prairie, Idaho. New geological, geochemical, and geophysical (gravity, magnetic, MT, seismic) studies of both areas led to a focus on Camas Prairie for validation during Phase 3. Camas Prairie is an EW-trending half-graben bounded on the north by the Idaho Batholith and on the south by the Mount Bennett Hills. Camas Prairie is bisected by a major NW-trending fault system (The Pothole fault) that separates NW-trending faults to east from ENE-trending faults …

A Geophysical Characterization Of The Structural Framework Of The Camas Prairie Geothermal System, Southcentral Idaho, Lee Liberty

CGISS Publications and Presentations

Play Fairway Analysis methods, utilizing existing geologic, thermal, geochemical, and geophysical data were employed in an initial assessment of geothermal resources in the Snake River Plain. These efforts identified the Camas Prairie in southcentral Idaho as a region with elevated resource potential. Subsequent efforts included structural and geophysical data collection to identify the most favorable structural settings for exploiting resources in the valley. The present work involved high-resolution gravity, magnetic, magnetotellurics (MT), field mapping, and seismic surveys to further characterize the system and target sites for exploration drilling around Barron’s Hot Springs (BHS) in the southwest part of the valley. …

Regional Greenland Accumulation Variability From Operation Icebridge Airborne Accumulation Radar, Gabriel Lewis, Erich Osterberg, Robert Hawley, Brian Whitmore, Hans Peter Marshall, Jason Box

CGISS Publications and Presentations

The mass balance of the Greenland Ice Sheet (GrIS) in a warming climate is of critical interest to scientists and the general public in the context of future sea-level rise. An improved understanding of temporal and spatial variability of snow accumulation will reduce uncertainties in GrIS mass balance models and improve projections of Greenland's contribution to sea-level rise, currently estimated at 0.089 ± 0.03 m by 2100. Here we analyze 25 NASA Operation IceBridge accumulation radar flights totaling > 17 700 km from 2013 to 2014 to determine snow accumulation in the GrIS dry snow and percolation zones over the past …

Probabilistic Inversion With Graph Cuts: Application To The Boise Hydrogeophysical Research Site, Guillaume Pirot, Niklas Linde, Grégoire Mariethoz, John H. Bradford

CGISS Publications and Presentations

Inversion methods that build on multiple-point statistics tools offer the possibility to obtain model realizations that are not only in agreement with field data, but also with conceptual geological models that are represented by training images. A recent inversion approach based on patch-based geostatistical resimulation using graph cuts outperforms state-of-the-art multiple-point statistics methods when applied to synthetic inversion examples featuring continuous and discontinuous property fields. Applications of multiple-point statistics tools to field data are challenging due to inevitable discrepancies between actual subsurface structure and the assumptions made in deriving the training image. We introduce several amendments to the original graph …

Geothermal Play Fairway Analysis Of The Snake River Plain: Phase 2, Lee M. Liberty

CGISS Publications and Presentations

Play Fairway Analysis (PFA) is a methodology adapted from the petroleum industry that integrates data at the regional or basin scale to define favorable plays for exploration in a systematic fashion. Phase 2 of our Play Fairway Analysis of the Western Snake River Plain (WSRP) province in southern Idaho had three primary goals: first, to fill data gaps in critical areas in order to better define potential prospects, second, to integrate these data into new thermal and structural models, and finally, to infer the location of potential resources and drilling targets that could be validated during Phase 3. Prospects in …

Impact Of Spatial Averaging On Radar Reflectivity At Internal Snowpack Layer Boundaries, N. Rutter, H.P. Marshall, K. Tape, R. Essery, J. King

CGISS Publications and Presentations

Microwave radar amplitude within a snowpack can be strongly influenced by spatial variability of internal layer boundaries. We quantify the impact of spatial averaging of snow stratigraphy and physical snowpack properties on surface scattering from near-nadir frequency-modulated continuous- wave radar at 12–18 GHz. Relative permittivity, density, grain size and stratigraphic boundaries were measured in-situ at high resolution along the length of a 9 m snow trench. An optimal range of horizontal averaging (4–6 m) was identified to attribute variations in surface scattering at layer boundaries to dielectric contrasts estimated from centimetre-scale measurements of snowpack stratigraphy and bulk layer properties. Single …

Monitoring Southwest Greenland’S Ice Sheet Melt With Ambient Seismic Noise, Aurélien Mordret, T. Dylan Mikesell, Christopher Harig, Bradley P. Lipovsky, Germán A. Prieto

CGISS Publications and Presentations

The Greenland ice sheet presently accounts for ~70% of global ice sheet mass loss. Because this mass loss is associated with sea-level rise at a rate of 0.7 mm/year, the development of improved monitoring techniques to observe ongoing changes in ice sheet mass balance is of paramount concern. Spaceborne mass balance techniques are commonly used; however, they are inadequate for many purposes because of their low spatial and/or temporal resolution. We demonstrate that small variations in seismic wave speed in Earth’s crust, as measured with the correlation of seismic noise, may be used to infer seasonal ice sheet mass balance. …

Advancements In The Measurement Of The Cryosphere Using Geophysics — Introduction, Andrew D. Parsekian, John Bradford, Georgios Tsoflias, Steven Arcone, Bernd Kulessa

CGISS Publications and Presentations

Frozen regions of the earth are known as the cryosphere. The arctic, Antarctica, permafrost, ice sheets, and glaciers are some of the most challenging places to measure subsurface parameters, but they can also be some of the most important places to science and engineering research due to their susceptibility to environmental change. Ground-based, airborne, and space-borne geophysical methods are deployed to observe targets below the ground or in ice that may be difficult or impossible to measure using conventional direct observations and measurements. The papers in this special section address recent advances in instrumentation development and deployment and computational capabilities …

A Synthetic Study To Assess The Applicability Of Full-Waveform Inversion To Infer Snow Stratigraphy From Upward-Looking Ground-Penetrating Radar Data, Lino Schmid, Jürg Schweizer, John Bradford, Hansruedi Maurer

CGISS Publications and Presentations

Snow stratigraphy and liquid water content are key contributing factors to avalanche formation. Upward-looking ground penetrating radar (upGPR) systems allow nondestructive monitoring of the snowpack, but deriving density and liquid water content profiles is not yet possible based on the direct analysis of the reflection response. We have investigated the feasibility of deducing these quantities using full-waveform inversion (FWI) techniques applied to upGPR data. For that purpose, we have developed a frequency-domain FWI algorithm in which we additionally took advantage of time-domain features such as the arrival times of reflected waves. Our results indicated that FWI applied to upGPR data …

Infrasound From Volcanic Rockfalls, Jeffrey B. Johnson, Timothy J. Ronan

CGISS Publications and Presentations

Proximal infrasound arrays can robustly track rapidly moving gravity-driven mass wasting, which occurs commonly at erupting volcanoes. This study reports on detection, localization, and quantification of frequent small rockfalls and infrequent pyroclastic density currents descending the southeast flanks of Santiaguito’s active Caliente Dome in January of 2014. Such activities are identified as moving sources, which descend several hundred meters at bulk flow speeds of up to ~10 m/s, which is considerably slower than the descent velocity of individual blocks. Infrasound rockfall signal character is readily distinguishable from explosion infrasound, which is manifested by a relatively fixed location source with lower …

Seasonal And Diurnal Cycles Of Liquid Water In Snow—Measurements And Modeling, A. Heilig, C. Mitterer, L. Schmid, N. Wever, J. Schweizer, H.-P. Marshall, O. Eisen

CGISS Publications and Presentations

Evaluating and improving snow models and outflow predictions for hydrological applications is hindered by the lack of continuous data on bulk volumetric liquid water content (��_w) and storage capacity of the melting snowpack. The combination of upward looking ground-penetrating radar and conventional snow height sensors enable continuous, nondestructive determinations of ��_w in natural snow covers from first surficial wetting until shortly before melt out. We analyze diurnal and seasonal cycles of ��_w for 4 years in a flat study site and for three melt seasons on slopes and evaluate model simulations for two different water …

Source Mechanism Of Small Long-Period Events At Mount St. Helens In July 2005 Using Template Matching, Phase-Weighted Stacking, And Full-Waveform Inversion, Robin S. Matoza, Bernard A. Chouet, Phillip B. Dawson, Peter M. Shearer, Matthew M. Haney, Gregory P. Waite, Seth C. Moran, T. Dylan Mikesell

CGISS Publications and Presentations

Long-period (LP, 0.5-5 Hz) seismicity, observed at volcanoes worldwide, is a recognized signature of unrest and eruption. Cyclic LP “drumbeating” was the characteristic seismicity accompanying the sustained dome-building phase of the 2004–2008 eruption of Mount St. Helens (MSH), WA. However, together with the LP drumbeating was a near-continuous, randomly occurring series of tiny LP seismic events (LP “subevents”), which may hold important additional information on the mechanism of seismogenesis at restless volcanoes. We employ template matching, phase-weighted stacking, and full-waveform inversion to image the source mechanism of one multiplet of these LP subevents at MSH in July 2005. The signal-to-noise …

Field Measurements For Remote Sensing Of The Cryosphere, Hans-Peter Marshall, Robert L. Hawley, Marco Tedesco

CGISS Publications and Presentations

Remote sensing observations of the cryosphere, like any other target of interest, require ground-based measurements for both calibration and validation, as inversion algorithms are usually underdetermined and uncertainties in the retrieval are needed for application. Field-based observations are performed in selected representative locations, and typically involve both direct in situ measurements of the physical properties of interest, as well as ground-based remote sensing techniques.

New state-of-the-art modern techniques for measuring physical properties rapidly and at high spatial resolution have recently given us a new view of spatiotemporal variability. These are important, as large variability at scales below the typical footprint …

Data Processing For Oscillatory Pumping Tests, Tania Bakhos, Michael Cardiff, Warren Barrash, Peter K. Kitanidis

CGISS Publications and Presentations

Characterizing the subsurface is important for many hydrogeologic projects such as site remediation and groundwater resource exploration. Methods based on the analysis of conventional pumping tests have the notable disadvantage that at a certain distance, the signal is small relative to the noise due to the effects of recharge, pumping in neighboring wells, change in the level or adjacent streams, and other common disturbances. This work focuses on oscillatory pumping tests in which fluid is extracted for half a period, then reinjected. We discuss a major advantage of oscillatory pumping tests: small amplitude signals can be recovered from noisy data …

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 …

Megathrust Splay Faults At The Focus Of The Prince William Sound Asperity, Alaska, Lee M. Liberty, Shaun P. Finn, Peter J. Haeussler, Thomas L. Pratt, Andrew Peterson

CGISS Publications and Presentations

[1] High-resolution sparker and crustal-scale air gun seismic reflection data, coupled with repeat bathymetric surveys, document a region of repeated coseismic uplift on the portion of the Alaska subduction zone that ruptured in 1964. This area defines the western limit of Prince William Sound. Differencing of vintage and modern bathymetric surveys shows that the region of greatest uplift related to the 1964 Great Alaska earthquake was focused along a series of subparallel faults beneath Prince William Sound and the adjacent Gulf of Alaska shelf. Bathymetric differencing indicates that 12 m of coseismic uplift occurred along two faults that reached the …

Using Svd For Improved Interferometric Green's Function Retrieval, Gabriela Melo, Alison Malcolm, Dylan Mikesell, Kasper Van Wijk

CGISS Publications and Presentations

Seismic interferometry (SI) is a technique used to estimate the Green’s function (GF) between two receiver locations, as if there were a source at one of the receiver locations. However, in many applications, the requirements to recover the exact GF are not satisfied and SI yields a poor estimate of the GF. For these non-ideal cases, we improve the interferometric GFs, by applying singular value decomposition (SVD) to the cross-correlations before stacking. The SVD approach preserves energy that is stationary in the cross-correlations, which is the energy that contributes most to the GF recovery, and attenuates non-stationary energy, which leads …

Aquifer Heterogeneity Characterization With Oscillatory Pumping: Sensitivity Analysis And Imaging Potential, M. Cardiff, T. Bakhos, P. K. Kitanidis, W. Barrash

CGISS Publications and Presentations

[1] Periodic pumping tests, in which a fluid is extracted during half a period, then reinjected, have been used historically to estimate effective aquifer properties. In this work, we suggest a modified approach to periodic pumping test analysis in which one uses several periodic pumping signals of different frequencies as stimulation, and responses are analyzed through inverse modeling using a “steady-periodic” model formulation. We refer to this strategy as multifrequency oscillatory hydraulic imaging. Oscillating pumping tests have several advantages that have been noted, including no net water extraction during testing and robust signal measurement through signal processing. Through numerical experiments, …

Changes In Elastic Wave Velocity And Rock Microstructure Due To Basalt-Co2-Water Reactions, Ludmila Adam, Kasper Van Wijk, Thomas Otheim, Michael Batzle

CGISS Publications and Presentations

The chemical interaction between carbon dioxide, water, and basalt is a common process in the earth, which results in the dissolution of primary minerals that later precipitate as alteration minerals. This occurs naturally in volcanic settings, but more recently basalts have been suggested as reservoirs for sequestration of anthropogenic CO₂. In both the natural and man-made cases, rock-fluid reactions lead to the precipitation of carbonates. Here, we quantify changes in ultrasonic wave speeds, associated with changes in the frame of whole-rock basalts, as CO₂ and basalt react. After 30weeks of reactions and carbonate precipitation, the ultrasonic wave …