Physical Sciences and Mathematics Commons^™

Groundwater System Characterisation: Fortescue Alluvial Fan, Michael J. Donn, Olga V. Barron, Axel Suckow, Chris Turnadge, John A. Simons, Robert J. Paul, Christopher Schelfhout Dr

Natural resources commissioned reports

This report focuses on groundwater system characterisation in the region north of Newman, based on analysis of pre-existing data and data newly acquired during project activities. Groundwater system characterisation was an important research component supporting the assessment of managed aquifer recharge opportunities, using mine dewatering surplus generated (due to mining below the watertable) at large BHP Billiton Iron Ore operations in the eastern Pilbara mining zone, and aiming to support irrigated agriculture. The assessment area is located north of Ethel Gorge and covers the Upper Fortescue River floodplain and surroundings. The project added much knowledge to this largely ‘data-poor’ region, …

Stratigraphic Architecture Of Pozuelo Mounds As Revealed By Earth Resistivity Tomography, Caeli Connolly

Honors College

This study is a geoarchaeological analysis using earth resistance tomography (ERT) surveys of two of four mounds at Pozuelo (Formative Period, cal yr 3000 BP) in the Chincha Valley of coastal, southern Peru. Layers identified in the subsurface were to determine the presence or absence of regional continuity between the mounds. This effort is part of a larger investigation examining the paleoenvironmental setting of the site, and its influence on site location and use. Ten earth resistance tomography profiles were collected using an ABEM Terrameter LS2 and 81 pin array. These profiles were then topographically corrected using topographic survey data …

Silicon And Oxygen In Earth’S Core: Applications Of Machine Learning To Metal-Silicate Equilibria And Core Formation, Ruben Keane

WWU Honors College Senior Projects

Within Earth’s core, light elements (Si, O, C, S, N, H) are known to make up a small fraction of the total mass of the core with respect to heavy elements. The degree to which these elements exist in the cores of terrestrial planets have geophysical and geochemical implications, most notably the presence of core convection and a geodynamo, thermal conductivity within the core, and core temperature. Comparison of the composition of chondrites to Earth’s mantle composition and the Preliminary Reference Earth Model have given an estimation of about 10 % light elements in Earth’s core. The concentrations of each …

Using Geologic Principles And Data To Create A Symphony For Wind Ensemble, Thomas J. Davis

2022 Symposium

Geology and geologic processes have featured in many pieces of classical music such as Alan Hovhaness’ symphonic works after Mt. St. Helens and Glacier Peak but these pieces are impressionistic and only take general inspiration from the geology of their subjects. This paper describes the creation of Tectonics (2021) a piece which incorporates and describes many geologic principles and processes, such as mantle convection, through the musical form of cannon, the formation of fold and thrust belts whose cross sections are directly visually recreated in the score, and the physical dynamics of earthquake waves demonstrated through variations of pitch and …

Processing Seismic, Gravity And Magnetic Data Over Diebold Knoll On Juan De Fuca Plate, Md Ariful Islam, Irina Filina

UNL Student Research Days Posters, Graduate

The Juan de Fuca plate is subducting beneath the North American plate along the Cascadia Subduction Zone. The CSZ is associated with multiple earthquakes, although compared to other subduction zones, it is less active. Many studies suggest that subduction will cause major mega thrust earthquake in near future, while others oppose this view. This project is focused on the Diebold Knoll on the JdF plate Various geophysical methods will be integrated to determine its origin and tectonic history.

Toward The Understanding Of The 2018 Arnold, Ne Earthquake Cluster: Relocation Of Hypocenters And Establishment Of New Gravity Base Stations, Kris Guthrie

Department of Earth and Atmospheric Sciences: Dissertations, Theses, and Student Research

In 2018, 27 unexpected earthquakes occurred in central Nebraska near Arnold. These earthquakes have no readily apparent cause and no evidence of the responsible fault system is seen at the surface. A joint analysis of vintage, regional gravity and magnetic datasets before the first earthquake revealed a lineament that could represent a pre-existing fault system. New high resolution potential fields data are needed to confidently map the subsurface features around the cluster. To design a survey and collect the data, the following tasks need to be addressed: [1] The earthquakes’ hypocenters need to be improved (relocated) to increase the efficiency …

Testing A Drone-Based Magnetic Field Surveying System, Erik Jacobson, Irina Filina

UCARE Research Products

Aeromagnetic surveys are conducted by geoscientists to study subsurface geologic structures, such as faults. This type of survey uses a magnetometer mounted upon an airborne vehicle to collect magnetic field data. Magnetic anomalies are caused by variations in subsurface geology, namely in magnetic properties of subsurface rocks. Jacobson and Filina (2019) reported on the development of a new low cost drone, based magnetic field surveying system by the UNL Geophysics Team. This drone-based magnetic system is capable of collecting high resolution data at low speeds and low altitudes. The current study focuses on testing this system by conducting two flights …

Improving Aquifer Characterization Through Integration Of Airborne Electromagnetics (Aem) And Well Hydrographs, Jacqueline Polashek

School of Natural Resources: Dissertations, Theses, and Student Research

The objective of this study is to evaluate methods of hydrostratigraphic modeling using geophysics and well hydrographs at the eastern edge of the High Plains aquifer (HPA) in Platte and Colfax counties within Nebraska, USA. The HPA is very heterogeneous in the study area, being hosted by architecturally complex glacial sediments and having many irregular hydraulic boundaries. Further, the HPA exhibits local variations between unconfined and confined conditions. Pumping in such bounded aquifers can be unsustainable because of cost increases and lost agricultural productivity. Moreover, the large drawdowns typical of confined aquifers can contribute to well interference during heavy pumping. …

Geophysical Analysis Of The Midcontinent Rift’S Subsurface Structure In Southeastern Nebraska, Patrick Szopinski

Honors Theses

The Midcontinent Rift System (MCRS) is a 1.1 billion-year-old failed rift system that spans much of the North American continental interior. The MCRS is exposed at Lake Superior and is buried in the subsurface along its southwest-extending arm through southeastern Nebraska. Due to the presence of buried volcanic rocks, the MCRS has characteristic highly-pronounced potential field anomalies (gravity and magnetic). Despite these large anomalies, not much is known about the subsurface faulting associated with the rift zone in the Midwest. The goal of this project is to attempt to use integrated analysis of collected geophysical data from multiple methods to …

Integrated Geophysical Analysis In The Northeastern Gulf Of Mexico, Mei Liu

Department of Earth and Atmospheric Sciences: Dissertations, Theses, and Student Research

The formation of the Gulf of Mexico (GoM) relates to the breakup of Pangea and opening of the central Atlantic Ocean. The tectonic history of the basin is still being debated due to lack of geological constraints. This project addresses the crustal architecture in the northeastern GoM from integrative analysis of multiple geophysical datasets to provide constraints for the tectonic reconstruction.

The objectives of this study are: 1) to delineate various tectonic zones (continental and oceanic domains) and map the boundary between them, 2) to derive physical properties of the subsurface rocks, 3) to map the major tectonic structures in …

Combining Hydraulic Head Analysis With Airborne Electromagnetics To Detect And Map Impermeable Aquifer Boundaries, Jesse T. Korus

Conservation and Survey Division

Impermeable aquifer boundaries affect the flow of groundwater, transport of contaminants, and the drawdown of water levels in response to pumping. Hydraulic methods can detect the presence of such boundaries, but these methods are not suited for mapping complex, 3D geological bodies. Airborne electromagnetic (AEM) methods produce 3D geophysical images of the subsurface at depths relevant to most groundwater investigations. Interpreting a geophysical model requires supporting information, and hydraulic heads offer the most direct means of assessing the hydrostratigraphic function of interpreted geological units. This paper presents three examples of combined hydraulic and AEM analysis of impermeable boundaries in glacial …

Aeromagnetic, Gravity, And Differential Interferometric Synthetic Aperture Radar Analyses Reveal The Causative Fault Of The 3 April 2017 MW 6.5 Moiyabana, Botswana, Earthquake, Folarin Kolawole, Estella A. Atekwana, S. Malloy, Dorothy Sarah Stamps, Raphael Grandin, Mohamed G. Abdel Salam, Khumo Leseane, Elisha M. Shemang

Geosciences and Geological and Petroleum Engineering Faculty Research & Creative Works

On 3 April 2017, a M_w 6.5 earthquake struck Moiyabana, Botswana, nucleating at >20 km focal depth within the Paleoproterozoic Limpopo-Shashe orogenic belt separating the Archean Zimbabwe and Kaapvaal Cratons. We investigate the lithospheric structures associated with this earthquake using high-resolution aeromagnetic and gravity data integrated with Differential Interferometric Synthetic Aperture Radar (DInSAR) analysis. Here we present the first results that provide insights into the tectonic framework of the earthquake. The ruptured fault trace delineated by DInSAR aligns with a distinct NW striking and NE dipping magnetic lineament within the Precambrian basement. The fault plane solution and numerical modeling …

St. Louis Area Earthquake Hazards Mapping Project: Seismic And Liquefaction Hazard Maps, Chris H. Cramer, Robert A. Bauer, Jaewon Chung, J. David Rogers, Larry Pierce, Vicki Voigt, Brad Mitchell, David Gaunt, Robert A. Williams, David J. Hoffman, Gregory L. Hempen, Phyllis J. Steckel, Oliver Salz Boyd, Connor M. Watkins, Kathleen B. Tucker, Natasha S. Mccallister

Geosciences and Geological and Petroleum Engineering Faculty Research & Creative Works

We present probabilistic and deterministic seismic and liquefaction hazard maps for the densely populated St. Louis metropolitan area that account for the expected effects of surficial geology on earthquake ground shaking. Hazard calculations were based on a map grid of 0.005°, or about every 500 m, and are thus higher in resolution than any earlier studies. To estimate ground motions at the surface of the model (e.g., site amplification), we used a new detailed near-surface shear-wave velocity model in a 1D equivalent- linear response analysis. When compared with the 2014 U.S. Geological Survey (USGS) National Seismic Hazard Model, which uses …

Imaging Bioturbation In Supratidal Carbonates: Non-Invasive Field Techniques Enhance Neoichnological And Zoogeomorphological Research, San Salvador, The Bahamas, Karen Kopcznski, Ilya Buynevich, H. Allen Curran, Jon Caris, Jonathan Nyquist

Geosciences: Faculty Publications

A case study in unconsolidated carbonates on San Salvador Island, The Bahamas, utilized high-frequency (800 MHz) georadar imaging to augment existing methodologies (burrow counts and measurements, casting) in brachyuran bioturbation research (Ocypode quadrata and Gecarcinus lateralis), and as part of a new dataset characterizing blue land crab (Cardisoma guanhumi) burrows. Non-invasive techniques such as ground-penetrating radar (GPR) can complement traditional field surveys aimed at quantifying mesoscale bioturbation in modern settings. These methods can establish diagnostic features for tracemaker identification and refine existing ichnofacies models. Drone-mounted aerial coverage provided the first high-resolution images of the micro-topography and …

An Iterative Approach To Ground Penetrating Radar At The Maya Site Of Pacbitun, Belize, Sheldon Skaggs, Terry G. Powis, Clara R. Rucker, George Micheletti

Publications and Research

Ground penetrating radar (GPR) surveys provide distinct advantages for archaeological prospection in ancient, complex, urban Maya sites, particularly where dense foliage or modern debris may preclude other remote sensing or geophysical techniques. Unidirectional GPR surveys using a 500 MHz shielded antenna were performed at the Middle Preclassic Maya site of Pacbitun, Belize. The survey in 2012 identified numerous linear and circular anomalies between 1 m and 2 m deep. Based on these anomalies, one 1 m × 4 m unit and three smaller units were excavated in 2013. These test units revealed a curved plaster surface not previously found at …

Seismic-Hazard Map Of Southeast Missouri And Likely Magnitude Of The February 1812 New Madrid Earthquake, Jaewon Chung, J. David Rogers

Geosciences and Geological and Petroleum Engineering Faculty Research & Creative Works

The New Madrid seismic zone lies beneath the upper Mississippi Embayment, straddling the border between southeastern Missouri and northwestern Tennessee. In late 1811 and early 1812, it produced five earthquakes of magnitudes >6.5, violently shaking the central and eastern United States (CEUS). Its magnitude and recurrence are of concern to today's central United States regions. By considering the effects of local geology, deterministic scenario maps (M_w 7.3 and 7.7) were produced for ground motions intended to simulate the 7 February 1812 event (NM3), which was the largest felt. These maps include spatial estimates of peak ground acceleration and of …

Seismic Slip Deficit In The Kashmir Himalaya From Gps Observations, Celia Schiffman, Bikram Singh Bali, Walter Szeliga, Roger Bilham

Faculty Scholarship for the Cascadia Hazards Institute

GPS measurements in Kashmir Himalaya reveal rangenormal convergence of 11±1 mm/yr with dextral shear of 5±1 mm/yr. The transition from a fully locked 170 km wide décollement to the unrestrained descending Indian plate occurs at ~25 km depth over an ~23 km wide transition zone. The convergence rate is consistent with the lower bounds of geological estimates for the Main Frontal Thrust, Riasi, and Balapora fault systems, on which no surface slip has been reported in the past millennium. Of the 14 damaging Kashmir earthquakes since 1123, none may have exceeded M_w = 7.6. Therefore, either a seismic moment …

Analyzing Tidal Fluctuations In The Big Pine Key Freshwater Lens With Time-Lapse Resistivity, Nicole M. Tucker

FIU Electronic Theses and Dissertations

The tidal influence on the Big Pine Key saltwater/freshwater interface was analyzed using time-lapse electrical resistivity imaging and shallow well measurements. The transition zone at the saltwater/freshwater interface was measured over part of a tidal cycle along three profiles. The resistivity was converted to salinity by deriving a formation factor for the Miami Oolite. A SEAWAT model was created to attempt to recreate the field measurements and test previously established hydrogeologic parameters. The results imply that the tide only affects the groundwater within 20 to 30 m of the coast. The effect is small and caused by flooding from the …

Evolution Of Dike Opening During The March 2011 Kamoamoa Fissure Eruption, Kīlauea Volcano, Hawai'i, Paul Lundgren, Michael Poland, Asta Miklius, Tim Orr, Sang-Ho Yun, Eric Fielding, Zhen Liu, Scott Hensley, Susan Owen, Akiko Tanaka, Walter Szeliga

Faculty Scholarship for the Cascadia Hazards Institute

The 5–9 March 2011 Kamoamoa fissure eruption along the east rift zone of Kīlauea Volcano, Hawai`i, followed months of pronounced inflation at Kīlauea summit. We examine dike opening during and after the eruption using a comprehensive interferometric synthetic aperture radar (InSAR) data set in combination with continuous GPS data. We solve for distributed dike displacements using a whole Kīlauea model with dilating rift zones and possibly a deep décollement. Modeled surface dike opening increased from nearly 1.5 m to over 2.8 m from the first day to the end of the eruption, in agreement with field observations of surface fracturing. …

Gps Constraints On Indo-Asian Convergence In The Bhutan Himalaya: Segmentation And Potential For A 8.2-8.8 Mw Earthquake, Dowchu Drukpa, Phuntsho Pelgay, Anjan Bhattacharya, Phillipe Vernant, Walter Szeliga, Roser Bilham

Faculty Scholarship for the Cascadia Hazards Institute

The seismogenic setting of Bhutan is unusual due to its lower-than-average 20th century seismic moment release (Drukpa et al. 2006), its absence of a reliable historical record, and its unusual location near the Shillong plateau where a great earthquake in 1897 resulted in ≈10 m of N/S shortening of the Indian plate to its south (Gahalaut et al. 2011). Despite these indicators that lower than normal convergence velocities should currently prevail, the GPS velocity between Shillong and Lhasa suggests that convergence in Sikkim and Bhutan occurs at velocities exceeding 20 mm/yr. GPS points between the Greater Himalaya and the Shillong …

Interseismic Strain Accumulation Along The Western Boundary Of The Indian Subcontinent, Walter Szeliga, Roger Bilham, Din Mohammad Kakar, Sarosh H. Lodi

Faculty Scholarship for the Cascadia Hazards Institute

Despite an overall sinistral slip rate of approximately 3 cm/yr, few major earthquakes have occurred in the past 200 years along the Chaman fault system, the western boundary of the India Plate with the Eurasia Plate. GPS and InSAR data reported here indicate sinistral shear velocities of 8-17 mm/yr across the westernmost branches of the fault system, suggesting that a significant fraction of the plate boundary slip is distributed in the fold and fault belt to the east. At its southernmost on-land segment, near the triple junction between the Arabia, Eurasia, and India Plates, we find the velocity across the …

Magnetic Susceptibility As A Proxy For Investigating Microbially Mediated Iron Reduction, Farag M. Mewafy, Estella A. Atekwana, D. Dale Werkema, Lee D. Slater, Dimitrios Ntarlagiannis, Andre Revil, Magnus E. Skold, Geoffrey N. Delin

Geosciences and Geological and Petroleum Engineering Faculty Research & Creative Works

We investigated magnetic susceptibility (MS) variations in hydrocarbon contaminated sediments. Our objective was to determine if MS can be used as an intrinsic bioremediation indicator due to the activity of iron-reducing bacteria. A contaminated and an uncontaminated core were retrieved from a site contaminated with crude oil near Bemidji, Minnesota and subsampled for MS measurements. The contaminated core revealed enriched MS zones within the hydrocarbon smear zone, which is related to iron-reduction coupled to oxidation of hydrocarbon compounds and the vadose zone, which is coincident with a zone of methane depletion suggesting aerobic or anaerobic oxidation of methane is coupled …

Future Cascadia Megathrust Rupture Delineated By Episodic Tremor And Slip, James S. Chapman, Timothy I. Melbourne

Faculty Scholarship for the Cascadia Hazards Institute

A suite of 15 episodic tremor and slip events imaged between 1997 and 2008 along the northern Cascadia subduction zone suggests future coseismic rupture will extend to 25 km depth, or approximately 60 km inland of the Pacific coast, rather than stopping offshore at 15 km depth. An ETS-derived coupling profile accurately predicts GPS measured interseismic deformation of the overlying North American plate, as measured by approximately 50 continuous GPS stations across western Washington State. When extrapolated over the 550-year average recurrence interval of Cascadia megathrust events, the coupling model also replicates the pattern and amplitude of coseismic coastal subsidence …

Seismic And Geodetic Constraints On Cascadia Slow Slip, Timothy I. Melbourne, Aaron G. Wech, Kenneth C. Creager

Faculty Scholarship for the Cascadia Hazards Institute

Automatically detected and located tremor epicenters from episodic tremor and slip (ETS) episodes in northern Cascadia provide a high-resolution map of Washington’s slow slip region. Thousands of epicenters from the past four ETS events from 2004 to 2008 provide detailed map-view constraints that correlate with geodetic estimates of the simultaneous slow slip. Each of these ETS events exhibits remarkable similarity in the timing and geographic distribution of tremor density and geodetically inferred slip. Analysis of the latest 15-month inter-ETS period also reveals ageodetic tremor activity similar both in duration and extent to ETS tremor. Epicenters from both ETS and inter- …

Fold And Thrust Partitioning In A Contracting Fold Belt: Insights From The 1931 Mach Earthquake In Baluchistan, Walter Szeliga, Roger Bilham, Daniel Schelling, Din Mohamed Kakar, Sarosh Lodi

Faculty Scholarship for the Cascadia Hazards Institute

Surface deformation associated with the 27 August 1931 earthquake near Mach in Baluchistan is quantified from spirit-leveling data and from detailed structural sections of the region interpreted from seismic reflection data constrained by numerous well logs. Mean slip on the west dipping Dezghat/Bannh fault system amounted to 1.2 m on a 42 km x 72 km thrust plane with slip locally attaining 3.2 m up dip of an inferred locking line at approximately 9 km depth. Slip also occurred at depths below the interseismic locking line. In contrast, negligible slip occurred in the 4 km near the interseismic locking line. …

Temporal Geophysical Signatures From Contaminant-Mass Remediation, Vukenkeng Che-Alota, Estella A. Atekwana, Eliot A. Atekwana, William August Sauck, D. Dale Werkema

Geosciences and Geological and Petroleum Engineering Faculty Research & Creative Works

We have previously documented changes in bulk electrical conductivity, self-potential (SP), and ground-penetrating-radar (GPR) reflections in a field setting caused by biogeochemical transformations of hydrocarbon-contaminated media. These transformations are associated with hydrocarbon biodegradation. The results of surface geophysical surveys acquired in 1996, 2003, and 2007 document changes in geophysical signatures associated with removing hydrocarbon mass in the contaminated zone. Initial investigations in 1996 showed that relative to background, the contaminated area was characterized by higher bulk electrical conductivity, positive SP anomaly, and attenuated GPR reflections. Repeated surveys in 2003 and 2007 over the contaminated area showed that in 2007, the …

Understanding The Depth And Nature Of Flow Systems In The Nashoba Terrane, Eastern Massachusetts, U.S.A., John P. Diggins

Masters Theses 1911 - February 2014

Igneous and metamorphic rock units have long been considered marginal aquifers yet they are a significant source for potable drinking water in many areas worldwide. Additionally, use of these systems is on the rise due to many factors including, contamination and overuse of surficial systems, as well as expanding population and drought. The Nashoba Terrane is a fault-bounded block of high-grade, steeply dipping metavolcanic and metasedimentary rock located in eastern Massachusetts, U.S.A. The Nashoba is northeast trending, extending from Oxford, MA to the Gulf of Maine south of Newburyport, MA. Seventeen previously drilled wells throughout the Nashoba were selected for …

Gps Constraints On 34 Slow Slip Events Within The Cascadia Subduction Zone, 1997–2005, Walter Szeliga, Timothy I. Melbourne, Marcelo Santillan, Meghan Miller

Faculty Scholarship for the Cascadia Hazards Institute

Refinements to GPS analyses in which we factor geodetic time series to better estimate both reference frames and transient deformation resolve 34 slow slip events located throughout the Cascadia subduction zone from 1997 through 2005. Timing of transient onset is determined with wavelet transformation of geodetic time series. Thirty continuous stations are included in this study, ranging from northern California to southwestern British Columbia. Our improvements in analysis better resolve the largest creep events and also identify many smaller events. At 48.5 degrees N latitude, a 14-month average recurrence interval has been observed over eight events since 1997. Farther north …

The Microbial Community Structure In Petroleum-Contaminated Sediments Corresponds To Geophysical Signatures, Jonathan P. Allen, Estella A. Atekwana, Eliot A. Atekwana, Joseph W. Duris, D. Dale Werkema, Silvia Rossbach

Geosciences and Geological and Petroleum Engineering Faculty Research & Creative Works

The interdependence between geoelectrical signatures at underground petroleum plumes and the structures of subsurface microbial communities was investigated. For sediments contaminated with light non-aqueous-phase liquids, anomalous high conductivity values have been observed. Vertical changes in the geoelectrical properties of the sediments were concomitant with significant changes in the microbial community structures as determined by the construction and evaluation of 16S rRNA gene libraries. DNA sequencing of clones from four 16S rRNA gene libraries from different depths of a contaminated field site and two libraries from an uncontaminated background site revealed spatial heterogeneity in the microbial community structures. Correspondence analysis showed …

Active-Passive Array Surface Wave Inversion And Comparison To Borehole Logs In Southeast Missouri, Alexei A. Malovichko, Neil Lennart Anderson, Dmitriy A. Malovichko, Denis Yu Shylakov, Pavel G. Butirin

Geosciences and Geological and Petroleum Engineering Faculty Research & Creative Works

In May 2002, both active- and passive-source surface wave data were acquired using 4-channel arrays at six selected bridge sites in southeast Missouri. Processing of acquired data (increase of signal-to-noise ratio, estimation of phase velocities) was carried out and dispersion curves of Rayleigh wave phase velocities were constructed. Each fundamental mode dispersion curve was then inverted by linearised optimization to a layered shear-wave velocity profile to depths of up to 60 m. The estimated shear-wave velocity profiles were compared to other geotechnical data that had been previously acquired at each test site for the Missouri Department of Transportation (MoDOT) including …