Geology topics

Paul A. Bedrosian

Publications and source records attributed to Paul A. Bedrosian.

At least 37 records · Page 2Linked to original sources

USGS2022

Trans-crustal structural control of CO2-rich extensional magmatic systems revealed at Mount Erebus Antarctica

Erebus volcano, Antarctica, with its persistent phonolite lava lake, is a classic example of an evolved, CO 2 -rich rift volcano. Seismic studies provide limited images of the magmatic system. Here we show using magnetotelluric data that a steep, melt-related conduit of low electrical resistivity originating in the upper mantle undergoes pronounced lateral re-orientation in the deep crust before reaching shallower magmatic storage and the summit lava lake. The lateral turn represents a structural fault-valve controlling episodic flow of magma and CO 2 vapour, which replenish and heat the high level phonolite differentiation zone. This magmatic valve lies within an inferred, east-west structural trend forming part of an accommodation zone across the southern termination of the Terror Rift, providing a dilatant magma pathway. Unlike H 2 O-rich subduction arc volcanoes, CO 2 -dominated Erebus geophysically shows continuous magmatic structure to shallow crustal depths of < 1 km, as the melt does not experience decompression-related volatile supersaturation and viscous stalling.

Nature Communications↗

USGS2022

Mapping a magnetic superstorm: March 1989 geoelectric hazards and impacts on United States power systems

A study is made of the relationships between geomagnetic and geoelectric field variation, Earth-surface impedance, and operational interference (anomalies) experienced on electric-power systems across the contiguous United States during the March 13-14, 1989 magnetic storm. For this, a 1-minute-resolution sequence of geomagnetic field maps is constructed from magnetometer time series acquired at ground-based observatories. Induced geoelectric field maps are calculated by convolving the geomagnetic maps with magnetotelluric impedance tensors. During the storm, anomalies were concentrated where the lithosphere is electrically resistive, and when and where geoelectric field amplitudes were high. This was particularly true in the Mid-Atlantic, Northeast, and the upper Midwest. Few anomalies were experienced in other parts of the Midwest and across much of the West, where the lithosphere is more conductive, and when and where geoelectric field amplitudes were low. Peak 1-minute-resolution geoelectric field amplitude ranged from 21.66 V/km in Maine and 19.02 V/km in Virginia to < 0.02 V/km in Idaho. Latitude-dependent organization of geoelectric hazards by auroral-zone electrojet currents is detectable, but it is much weaker than geographic organization due to surface impedance. Hazardous geoelectric fields were induced during different storm phases, at different local times, and, by inference, by a variety of ionospheric currents. Compared to geoelectric field amplitudes realized across the United States during March 1989, hazard maps used by utility companies to estimate systems exposure have much less geographic detail and a much smaller maximum-to-minimum range in geoelectric field amplitude. Future research will benet from denser geomagnetic monitoring, additional magnetotelluric surveying, and access to power-system impact data.

Space Weather↗

USGS2022

Lithospheric conductors reveal source regions of convergent margin mineral systems

The clean energy transition will require a vast increase in metal supply, yet new mineral deposit discoveries are declining, due in part to challenges associated with exploring under sedimentary and volcanic cover. Recently, several case studies have demonstrated links between lithospheric electrical conductors imaged using magnetotelluric (MT) data and mineral deposits, notably Iron Oxide Copper Gold (IOCG). Adoption of MT methods for exploration is therefore growing but the general applicability and relationship with many other deposit types remains untested. Here, we compile a global inventory of MT resistivity models from Australia, North and South America, and China and undertake the first quantitative assessment of the spatial association between conductors and three mineral deposit types commonly formed in convergent margin settings. We find that deposits formed early in an orogenic cycle such as volcanic hosted massive sulfide (VHMS) and copper porphyry deposits show weak to moderate correlations with conductors in the upper mantle. In contrast, deposits formed later in an orogenic cycle, such as orogenic gold, show strong correlations with mid-crustal conductors. These variations in resistivity response likely reflect mineralogical differences in the metal source regions of these mineral systems and suggest a metamorphic-fluid source for orogenic gold is significant. Our results indicate the resistivity structure of mineralized convergent margins strongly reflects late-stage processes and can be preserved for hundreds of millions of years. Discerning use of MT is therefore a powerful tool for mineral exploration.

North America, South America↗

USGS2022

Characteristics and sources of intense geoelectric fields in the United States: Comparative analysis of multiple geomagnetic storms

Intense geoelectric fields during geomagnetic storms drive geomagnetically induced currents in power grids and other infrastructure, yet there are limited direct measurements of these storm-time geoelectric fields. Moreover, most previous studies examining storm-time geoelectric fields focused on single events or small geographic regions, making it difficult to determine the typical source(s) of intense geoelectric fields. We perform the first comparative analysis of (a) the sources of intense geoelectric fields over multiple geomagnetic storms, (b) using 1-s cadence geoelectric field measurements made at (c) magnetotelluric survey sites distributed widely across the United States. Temporally localized intense perturbations in measured geoelectric fields with prominences (a measure of the relative amplitude of geoelectric field enhancement above the surrounding signal) of at least 500 mV/km were detected during geomagnetic storms with Dst minima ( Dst min ) of less than −100 nT from 2006 to 2019. Most of the intense geoelectric fields were observed in resistive regions with magnetic latitudes greater than 55° even though we have 167 sites located at lower latitudes during geomagnetic storms of −200 nT ≤ Dst min < −100 nT . Our study indicates intense short-lived (<1 min) and geoelectric field perturbations with periods on the order of 1–2 min are common. Most of these perturbations cannot be resolved with 1-min data because they correspond to higher frequency or impulsive phenomena that vary on timescales shorter than that sampling interval. The sources of geomagnetic perturbations inducing these intense geoelectric fields include interplanetary shocks, interplanetary magnetic field turnings, substorms, and ultralow frequency waves.

Space Weather↗

USGS2022

Geophysical imaging of the Yellowstone hydrothermal plumbing system

The nature of Yellowstone National Park’s plumbing system linking deep thermal fluids to its legendary thermal features is virtually unknown. The prevailing concepts of Yellowstone hydrology and chemistry are that fluids reside in reservoirs with unknown geometries, flow laterally from distal sources and emerge at the edges of lava flows . Here we present a high-resolution synoptic view of pathways of the Yellowstone hydrothermal system derived from electrical resistivity and magnetic susceptibility models of airborne geophysical data . Groundwater and thermal fluids containing appreciable total dissolved solids significantly reduce resistivities of porous volcanic rocks and are differentiated by their resistivity signatures . Clay sequences mapped in thermal areas and boreholes typically form at depths of less than 1,000 metres over fault-controlled thermal fluid and/or gas conduits . We show that most thermal features are located above high-flux conduits along buried faults capped with clay that has low resistivity and low susceptibility. Shallow subhorizontal pathways feed groundwater into basins that mixes with thermal fluids from vertical conduits. These mixed fluids emerge at the surface, controlled by surficial permeability, and flow outwards along deeper brecciated layers. These outflows, continuing between the geyser basins, mix with local groundwater and thermal fluids to produce the observed geochemical signatures. Our high-fidelity images inform geochemical and groundwater models for hydrothermal systems worldwide.

Idaho, Wyoming↗

USGS2022

Graphite as an electrically conductive indicator of ancient crustal-scale fluid flow within mineral systems

Magnetotelluric (MT) imaging results from mineral provinces in Australia and in the United States show an apparent spatial relationship between crustal-scale electrical conductivity anomalies and major magmatic-hydrothermal iron oxide-apatite/iron oxide-copper-gold (IOA-IOCG) deposits. Although these observations have driven substantial interest in the use of MT data to image ancient fluid pathways, the exact cause of these anomalies has been unclear. Here, we interpret the conductors to be the result of graphite precipitation from CO 2 -rich magmatic fluids during cooling. These fluids would have exsolved from mafic magmas at mid- to lower-crustal depths; saline magmatic fluids that could drive mineralization were likely derived from related, more evolved intrusions at shallower crustal levels. In our model, the conductivity anomalies then mark zones that once were the deep roots of ancient magmatic-hydrothermal mineral systems.

Earth and Planetary Science Letters↗

USGS2022

Geoelectric constraints on the Precambrian assembly and architecture of southern Laurentia

Using images from an updated and expanded three-dimensional electrical conductivity synthesis model for the contiguous United States (CONUS), we highlight the key continent-scale geoelectric structures that are associated with the Precambrian assembly of southern Laurentia. Conductivity anomalies are associated with the Trans-Hudson orogen, the Penokean suture, the ca. 1.8–1.7 Ga Cheyenne belt and Spirit Lake tectonic zone, and the Grenville suture zone; the geophysical characteristics of these structures indicate that the associated accretionary events involved the closure of ancient ocean basins along discrete, large-scale structures. In contrast, we observe no large-scale conductivity anomalies through the portion of southern Laurentia that is generally viewed as composed of late Paleoproterozoic–early Mesoproterozoic accretionary crust. The lack of through-going conductors places constraints on the structure, petrology, and geodynamic history of crustal growth in southern Laurentia during that time period. Overall, our model highlights the enigmatic nature of the concealed Precambrian basement of much of southern Laurentia, as it in some places supports and in other places challenges prevailing models of Laurentian assembly. The revised CONUS electrical conductivity model thus provides important constraints for testing new models of Precambrian tectonism in this region.

continental United States↗

USGS2022

The nature of active magma reservoirs and storage underneath Cascade volcanoes

Volcanic activity is dictated by crustal pathways and reservoirs through which magma ascends and collects. The Cascades are a natural laboratory to study the variable influence of the crust on erupted magmas and subterranean reservoirs. The interaction of ascending magmas with the variable subduction geometry and tectonics of the overriding North American plate have given rise to a diversity of magmatic storage conditions. The confluence of geochemical and geophysical investigations emphasizes that most magmatic systems in the Cascades have been built at multiple levels in the crust as determined by tectonics, pre-existing structure, and magmatic flux from the mantle.

British Columbia, California, Oregon, Washington↗

USGS2021

Down to Earth with nuclear electromagnetic pulse: Realistic surface impedance aﬀects mapping of the E3 geoelectric hazard

An analysis is made of Earth-surface geoelectric fields and voltages on electricity transmission power-grids induced by a late-phase E3 nuclear electromagnetic pulse (EMP). A hypothetical scenario is considered of an explosion of several hundred kilotons set several hundred kilometers above the eastern-midcontinental United States. Ground-level E3 geoelectric fields are estimated by convolving a standard parameterization of E3 geomagnetic field variation with magnetotelluric Earth-surface impedance tensors derived from wideband measurements acquired across the study region during a recent survey. These impedance tensors are a function of subsurface three-dimensional electrical conductivity structure. Results, presented as a movie-map, demonstrate that localized differences in surface impedance strongly distort the amplitude, polarization, and variational phase of induced E3 geoelectric fields. Locations with a high degree of E3 geoelectric polarization tend to have high geoelectric amplitude. Uniform half-space models and one-dimensional, depth-dependent models of Earth-surface impedance, such as those widely used in government and industry reports informing power-grid vulnerability assessment projects, do not provide accurate estimates of the E3 geoelectric hazard in complex geological settings. In particular, for the Eastern-Midcontinent, half-space models can lead to (order-one) overestimates/underestimates of EMP-induced geovoltages on parts of the power grid by as much as 1,000 volts (a range of 2,000 volts)—comparable to the amplitudes of the geovoltages themselves.

Earth and Space Sciences↗

USGS2021

Magnetotelluric sampling and geoelectric hazard estimation: Are national-scale surveys sufficient?

At present, the most reliable information for inferring storm-time ground electric fields along electrical transmission lines comes from coarsely sampled, national-scale magnetotelluric (MT) data sets, such as that provided by the EarthScope USArray program. An underlying assumption in the use of such data is that they adequately sample the spatial heterogeneity of the surface relationship between geomagnetic and geoelectric fields. Here, we assess the degree to which the density of MT data sampling affects geoelectric hazard assessments. For electrical transmission networks in each of four focus regions across the contiguous United States, we perform two parallel band-limited (10 1 –10 3 s) hazard analyses: one using only USArray-style (∼70-km station spacing) MT data, and one incorporating denser (≪70-km station spacing) MT data. We find that the use of USArray-style MT sampling alone provides a useful first-order estimate of integrated geoelectric fields along electrical transmission lines. However, we also find that the use of higher density MT data can in some areas lead to order-of-magnitude differences in line-averaged electric field estimates at the level of individual transmission lines and can also yield significant differences in subregional hazard patterns. As we demonstrate using variogram plots, these differences reflect short-spatial-scale variability in Earth conductivity, which in turn reflects regional lithotectonic structure and history. We also provide a cautionary example in the use of electrical conductivity models to predict dense MT data; although valuable for hazard applications, models may only be able to reproduce surface geoelectric fields as captured by the MT data from which they were derived.

AGU Space Weather↗

USGS2021

When Wyoming became Superior: Oblique convergence along the southern Trans-Hudson orogen

The Trans-Hudson orogen (THO) is one of the best-preserved Proterozoic orogens on Earth, largely unaffected by subsequent tectonism, yet its southern extent lies concealed beneath the North American Central Plains. A new 3D resistivity model over the southern orogen is developed and interpreted alongside borehole, potential field, and seismic reflection data. We present the first synoptic crustal view of the southern THO and a new tectonic model of the orogen. Our model reveals high-conductivity belts marking paleo-subduction zones while the orogen center consists of deeply exhumed relatively dense and mostly magnetic juvenile crust preserved between the deformed margins of the Wyoming and Superior cratons. Complex structure along the western margin suggests convergence began with oblique subduction and the northward transport of severed fragments of the Wyoming Province. High-conductivity belts are in places offset from upper-crustal geophysical boundaries, consistent with the thrusting of accreted arcs over the Archean margins during terminal closure.

Geophysical Research Letters↗

USGS2021

Airborne geophysical surveys of the lower Mississippi Valley demonstrate system-scale mapping of subsurface architecture

The Mississippi Alluvial Plain hosts one of the most prolific shallow aquifer systems in the United States but is experiencing chronic groundwater decline. The Reelfoot rift and New Madrid seismic zone underlie the region and represent an important and poorly understood seismic hazard. Despite its societal and economic importance, the shallow subsurface architecture has not been mapped with the spatial resolution needed for effective management. Here, we present airborne electromagnetic, magnetic, and radiometric observations, measured over more than 43,000 flight-line-kilometers, which collectively provide a system-scale snapshot of the entire region. We develop detailed maps of aquifer connectivity and shallow geologic structure, infer relationships between structure and groundwater age, and identify previously unseen paleochannels and shallow fault structures. This dataset demonstrates how regional-scale airborne geophysics can close a scale gap in Earth observation by providing observational data at suitable scales and resolutions to improve our understanding of subsurface structures.

Mississippi River Valley↗

USGS2021

Airborne geophysical imaging of weak zones on Iliamna Volcano, Alaska: Implications for slope stability

Water‐saturated, hydrothermally altered rocks reduce the strength of volcanic edifices and increase the potential for sector collapses and far‐traveled mass flows of unconsolidated debris. Iliamna Volcano is an andesitic stratovolcano located on the western side of the Cook Inlet, ∼225 km southwest of Anchorage and is a source of repeated avalanches. The widespread snow and ice cover on Iliamna Volcano make surface alteration difficult to identify. However, intense hydrothermal alteration significantly reduces both the electrical resistivity and magnetization of volcanic rock and can therefore be identified with airborne geophysical measurements. We use airborne electromagnetic and magnetic data to map snow and ice thickness and identify underlying alteration zones at Iliamna Volcano, Alaska. Resistivities were calculated to an average depth of >300 m, and a 3‐D susceptibility model extends from the surface to the base of the volcano, about 3,000 m below the summit. Geophysical models image low resistivity (<30 ohm‐m) and low susceptibilities near the summit of Iliamna and below its older vent complex, with the low susceptibilities indicating alteration up to ∼800 m in thickness. Thin conductors (∼50–100 m thick) on the edifice slopes coincide with recorded locations of repeated debris avalanches over the past ∼60 years and are attributed to saturated zones at high elevation. Three‐dimensional slope stability models based upon the geophysically constrained alteration distribution suggest the edifice of Iliamna is unstable and could lead to collapse scars ∼400 m deep near the current and former vent complexes.

Alaska↗

USGS2021

3-D Modeling of the Duluth Complex from geophysical data

The Mesoproterozoic Duluth Complex in northeastern Minnesota is one of the major plutonic components of the Midcontinent Rift System and hosts a variety of copper-nickel sulfide and platinum-group element deposits. The Duluth Complex is composed of a series of individual mafic and felsic intrusions emplaced 1110-1098 Ma within Paleoproterozoic sedimentary rocks of the Animikie basin and volcanic flows of the Midcontinent Rift. Prior work has included 2-D modeling and qualitative geologic interpretations of gravity and magnetic data (e.g., Chandler, 1990; Chandler and Ferderer, 1989), much of which is still preliminary (V. Chandler, written commun., 2020). Three-dimensional modeling has been limited, with only one 3-D model created using Bouguer gravity data constrained by seismic-reflection interpretations as part of a PhD thesis (Allen, 1994). Given the complex geology of the area, 3-D modeling is useful for providing a complete picture of the variable densities, susceptibilities, and electrical resistivities throughout the Duluth Complex and associated volcanic rocks as well as their depth extent beneath sedimentary cover. Models of these geophysical properties at depth enable more accurate geologic mapping in the subsurface which can lead to an improved understanding of the formation history of the Duluth Complex. In this study, we use aeromagnetic data acquired between 1979-1991 (Chandler, 2007), Bouguer gravity data collected since 1950 (Chandler and Lively, 2019), and magnetotelluric data collected in 2019 to create new 2-D and 3-D geophysical models of the Duluth Complex constrained by seismic reflection, geologic, and rock property data. An inversion of the Bouguer gravity data for thickness of the Duluth Complex using constant densities of 3110 kg/m3 and 2670 kg/m3 for the Duluth Complex and surrounding crustal rocks, respectively, results in thicknesses ranging from ~3-28 km for the Duluth Complex and related intrusions and volcanic rocks (Figure 1A). A 3-D model of the magnetotelluric data reveals low resistivity anomalies at ~5-10 km depth below the northern margin of the Duluth Complex and below the Greenwood Lake intrusion (Figure 1B). We expect to encounter low resistivities at depth associated with the Paleoproterozoic Animikie basin, which makes up the floor of the Duluth Complex, and therefore interpret these anomalies as either the base of the complex or as fragments of Animikie sediments interfingered with igneous intrusive rocks. Finally, 3-D voxel models of density and susceptibility illuminate the subsurface distribution of rock properties below the Duluth Complex which, in combination with resistivity and thickness models, can be used to create a 3-D geologic map of this area.

Minnesota↗

USGS2021

Geophysical mapping of the eastern arm of the Midcontinent Rift in Upper Michigan

The western arm of the Midcontinent Rift system (MRS), extending from western Lake Superior into Kansas, is somewhat understood due to exposures of MRS igneous and sedimentary rocks, scattered drillholes, and seismic imaging. By comparison, the eastern rift arm (ERA), extending from eastern Lake Superior through lower Michigan is poorly understood, with almost no outcropping volcanic rocks south of Michipicoten and Caribou Islands in Lake Superior, few boreholes that intercept MRS rocks, and the majority of the rift concealed beneath the Michigan basin. Seismic reflection sections in the eastern Lake Superior basin (Behrendt et al., 1990; Mariano and Hinze, 1994) and in Lake Michigan (Cannon et al., 1991) provide the primary constraints on the gross structure of the ERA. Unlike the western rift arm, which has undergone considerable post-rift shortening (Cannon et al., 1993), the ERA appears less deformed, with perhaps as much as 30 km of rift clastic and volcanic rocks preserved within a broadly symmetric basin (Behrendt et al., 1988) and about 5 km of anticlinal relief on the top of the volcanics (Mariano and Hinze, 1994). Under Paleozoic cover in the eastern part of Michigan’s Upper Peninsula, the ERA western margin is inferred from linear northwest-trending potential-field anomalies (Fig. 1a, b) and the occurrence of MRS basalt and rhyolite beneath ~2 km of Paleozoic and MRS clastic rocks in the St. Amour drillhole (Ojakangas and Dickas, 2002).

Michigan↗

USGS2020

Simultaneous observations of geoelectric and geomagnetic ﬁelds produced by magnetospheric ULF waves

Geomagnetic perturbations ( B G E O ) related to magnetospheric ultralow frequency (ULF) waves induce electric fields within the conductive Earth—geoelectric fields ( E G E O )—that in turn drive geomagnetically induced currents. Though numerous past studies have examined ULF wave B G E O from a space weather perspective, few studies have linked ULF waves with E G E O . Using recently available magnetotelluric impedance and E G E O measurements in the contiguous United States, we explore the relationship between ULF waves and E G E O . We use satellite, ground‐based radar, B G E O , and E G E O measurements in a case study of a plasmaspheric virtual resonance (PVR), demonstrating that the PVR E G E O has significant spatial variation in contrast to a relatively uniform B G E O , consistent with spatially varying Earth conductivity. We further show ULF wave E G E O measurements during two moderate storms of ∼ 1 V/km. We use both results to highlight the need for more research characterizing ULF wave E G E O .

Geophysical Research Letters↗

USGS2020

Quantifying model structural uncertainty using airborne electromagnetic data

The ability to quantify structural uncertainty in geological models that incorporate geophysical data is affected by two primary sources of uncertainty: geophysical parameter uncertainty and uncertainty in the relationship between geophysical parameters and geological properties of interest. Here, we introduce an open-source, trans-dimensional Bayesian Markov chain Monte Carlo (McMC) algorithm GeoBIPy—Geophysical Bayesian Inference in Python—for robust uncertainty analysis of time-domain or frequency-domain airborne electromagnetic (AEM) data. The McMC algorithm provides a robust assessment of geophysical parameter uncertainty using a trans-dimensional approach that lets the AEM data inform the level of model complexity necessary by allowing the number of model layers itself to be an unknown parameter. Additional components of the Bayesian algorithm allow the user to solve for parameters such as data errors or corrections to the measured instrument height above ground. Probability distributions for a user-specified number of lithologic classes are developed through posterior clustering of McMC-derived resistivity models. Estimates of geological model structural uncertainty are thus obtained through the joint probability of geophysical parameter uncertainty and the uncertainty in the definition of each class. Examples of the implementation of this algorithm are presented for both time-domain and frequency-domain AEM data acquired in Nebraska, USA.

Nebraska↗

USGS2020

Linking magma storage and ascent to eruption volume and composition at an arc caldera

Conceptual models of magma storage and transport under calderas favor a connected system of sills and dikes. These features are individually below the resolution of standard seismic tomography, but radial seismic anisotropy can reveal where they exist in aggregate. We model radial anisotropy at Okmok caldera, Alaska, to demonstrate the presence of a caldera-centered stacked sill complex and surrounding dike system. We show that ascending magma, inferred from seismicity, either intersects the sill complex, resulting in a larger volume eruption of evolved magma, or bypasses the overlying sill complex via dikes, resulting in a low-volume mafic eruption. Our results exemplify how the locations of magma storage and paths of transport impact eruption size and composition. As this type of crustal storage is likely common to many calderas, this analysis offers a potential new framework for volcano observatories to forecast the size of impending eruptions.

Alaska↗