Geology Reports⌕ Search

Geology topics

Paul A. Bedrosian

Publications and source records attributed to Paul A. Bedrosian.

At least 19 recordsLinked to original sources
USGS2026

The United States Magnetotelluric Array and the National Impedance Map

The United States Magnetotelluric Array (USMTArray) data set, collected in the years 2006–2024, consists of more than 1,700 long-period magnetotelluric stations covering the entirety of the contiguous United States on a quasi-regular 70 km grid. Funding across multiple federal agencies was critical to sustaining this effort to its completion. Important components of the project included active guidance and participation from the MT community, the open and timely availability of all data, and the application of consistent instrumentation and robust data processing. Together with parallel advancement in the development of publicly available three-dimensional (3D) inversion codes, the USMTArray has revitalized the US magnetotelluric community and increased the visibility of magnetotellurics within the Earth-science community. Taken as a whole, these data are visualized as the National Impedance Map, which, together with a 3D synthesis conductivity model of the nation, reveals the electrical architecture of the contiguous US. USMTArray data are used by researchers worldwide for fundamental and applied studies, including investigations of continental architecture and evolution, estimation of hazards to critical infrastructure due to geomagnetic storms, and assessment of the nation's undiscovered geothermal and mineral resources. We here review the history and development of the project, discuss the challenges and successes in its execution, present the National Impedance Map and synthesis conductivity model, and highlight the breadth of research stemming from this rich data set.

conterminous United States↗
USGS2026

Magnetic storms and geoelectric hazards

Magnetic storms induce geoelectric fields at Earth's surface that can interfere with grounded long-line systems. The September 1859 storm disrupted global telegraph operations, the March 1989 storm caused a blackout in Canada and interfered with electric-power-transmission systems in the United States, and other storms have had related impacts. The geographic and temporal dependence of geoelectric fields are functions of both geomagnetic variation and local surface impedance, which differ considerably across different geological regions. These dependencies can be mapped across the contiguous United States by combining magnetotelluric impedance tensors with ground magnetometer time series. This review illustrates such mapping for the 1989 storm and shows that power-system interference was experienced where surface impedance is high, and when and where geoelectric fields were intense. Statistical analyses indicate that storms comparable to that of March 1989 occur roughly once every four solar cycles. Ongoing developments in numerical modeling and real-time monitoring are anticipated to enable prediction of geoelectric hazards. ▪ Magnetic storms can induced electric fields in the solid Earth that interfere with electric-power-transmission systems. ▪ Geoelectric hazards depend on the storm-time geomagnetic disturbance and the electrical conductivity structure of Earth. ▪ Historically, impacts on telecommunication and power-transmission systems in the United States have been concentrated in the East and Midwest. ▪ The future occurrence of a magnetic superstorm could cause widespread disruption of electric-power-transmission systems.

Annual Review of Earth and Planetary Sciences↗
USGS2025

Mapping a Carrington storm

A map is presented of median 1-min-resolution peak geoelectric-field strength across the United States as would be induced by magnetic storms as intense as the 2 September 1859 Carrington storm. The map is constructed from two data sets: Magnetometer time series from 22 ground-based observatories recording 40 magnetic storms, and surface impedance tensors derived from magnetotelluric measurements acquired at 1616 survey sites across the contiguous United States. Carrington-class storm geoelectric fields are likely to be very strong in the United States East and Midwest; > 5.00 V/km at many places. In Virginia, strengths would likely range from 30.30 V/km, with a 68% confidence interval of [19.44,47.20] V/km, to as low as 0.05 [0.03,0.07] V/km. Comparison of model geopotentials with those measured on 30 long lines, indicates errors of about 18%. A Carrington-class storm would likely induce geoelectric fields with strengths 55% greater than for the 13–14 March 1989 storm.

contiguous United States↗
USGS2025

The impact of the May 1921 superstorm on American telecommunication systems

A compilation is presented of impacts (interference and damage) realized on long-line telegraph and telephone systems across North America during the 13-16 May 1921 magnetic storm. Impacts occurred primarily during local nighttime, after the third of four sudden commencements, and during the storm’s most-prominent main phase. Impacts are attributed to rapid and high-amplitude geomagnetic field variation generated by substorms. This induced potential di erences and between the grounding points of communication networks that were su cient to cause system interference and damage. In the United States, impacts were concentrated in the Midwest and in the East, regions characterized by high electromagnetic surface impedance. Given technological changes, modern telecommunication systems are less exposed to storms like that of May 1921, while power-grid systems are now more exposed to them.

Space Weather↗
USGS2025

Multi-scale geophysical imaging of a hydrothermal system in Yellowstone National Park, USA

Little is known about the local plumbing systems that fuel Yellowstone’s famous hot springs, geysers and mud pots. A multi-method, multi-scale geophysical investigation was carried out in the Obsidian Pool Thermal Area (OPTA) to: (i) delineate the lateral extent of the hydrothermal area and associated surface features; (ii) estimate the dimensions of the upflow zone and identify its main controlling structures; (iii) assess fluids circulation pathways from depth to surface. Ground and airborne geophysical data were acquired to connect local and regional scales, from shallow to large depths. Maps of surface electrical resistivity show a strong correlation with hydrothermal features. At in-termediate depths, electrical resistivity permits delineating the upper limit of the upflow zone, while Poisson’s ratio highlights differences in subsurface fluid content. Combining these results with surface observations and topographic information, we speculate that differential mixing of hydrothermal and fresh water could explain the wide diversity of features observed at OPTA. Low electrical resistivity observed at large depths also suggest that a vast upflow zone, controlled by rhyolite flows and conjugate faults, underlies the OPTA. We speculate that hydrothermal fluids rise along fractures and reach the surface in topographic lows to form hydrothermal features. Our results show that synoptic, multi-scale geophysical measurements provide a roadmap for understanding where and how geologic heterogeneity, topography, fluid-gas separation, and the mixing of thermal and meteoric waters conspire to produce the wide variety of Yellowstone’s renowned hydrothermal features.

Wyoming↗
USGS2025

The progression of basaltic–rhyolitic melt storage at Yellowstone Caldera

Yellowstone Caldera is one of the largest volcanic systems on Earth, hosting three major caldera-forming eruptions in the past two million years, interspersed with periods of less explosive, smaller-volume eruptions 1 . Caldera-forming eruptions at Yellowstone are sourced by rhyolitic melts stored within the mid- to upper crust. Seismic tomography studies have suggested that a broad region of rhyolitic melt extends beneath Yellowstone Caldera, with an estimated melt volume that is one to four times greater than the eruptive volume of the largest past caldera-forming eruption, and an estimated melt fraction of 6–28 per cent 2 , 3 , 4 , 5 . Seismic velocity is strongly influenced by temperature, pressure and melt; however, magnetotelluric data are primarily sensitive to the presence of melt, making these data ideal for constraining volcanic systems. Here we utilize magnetotelluric data to model the resistivity structure of Yellowstone Caldera’s crustal magma reservoir and constrain the region’s potential for producing major volcanic eruptions. We find that rhyolitic melts are stored in segregated regions beneath the caldera with low melt fractions, indicating that the reservoirs are not eruptible. Typically, these regions have melt volumes equivalent to small-volume post-caldera Yellowstone eruptions. The largest region of rhyolitic melt storage, concentrated beneath northeast Yellowstone Caldera, has a storage volume similar to the eruptive volume of Yellowstone’s smallest caldera-forming eruption. We identify regions of basalt migrating from the lower crust, merging with and supplying heat to the northeast region of rhyolitic melt storage. On the basis of our analysis, we suggest that the locus of future rhyolitic volcanism has shifted to northeast Yellowstone Caldera.

Idaho, Montana, Wyoming↗
USGS2025

Synergy between geology and geophysics in graphite mineral resource assessment

Graphite is designated as a critical mineral by the U.S. Government due to its essential role in modern technology and its vulnerability to supply chain disruption. To evaluate domestic graphite resources, the U.S. Geological Survey (USGS) conducted two airborne electromagnetic (AEM) surveys as part of the Earth Mapping Resources Initiative (Earth MRI) over prospective areas in the southern Appalachians of Alabama and the Seward Peninsula of Alaska. These surveys aimed to refine areas of interest by revealing subsurface electrical resistivity patterns indicative of graphite-bearing rocks. The results of the AEM surveys were integrated with existing geological knowledge and supplemented by field observations to assess the likelihood of undiscovered resources in these regions. This multidisciplinary approach combined geophysical techniques with traditional geological methods, enhancing the understanding of graphite potential in areas that may have been previously overlooked. The findings from this initiative are expected to contribute significantly to the ongoing efforts to secure a reliable domestic supply of graphite, thereby reducing import reliance and supporting national interests in mineral resource management.

Conference Paper↗
USGS2024

Geoelectric evidence for a wide spatial footprint of active extension in central Colorado

Three-dimensional magnetotelluric (MT) imaging in central Colorado reveals a set of north-striking high-conductivity tracks at lower-crustal (50–20 km) depths, with conductive finger-like structures rising off these tracks into the middle crust (20–5 km depth). We interpret these features to represent saline aqueous fluids and partial melt that are products of active extensional tectonomagmatism. These conductors are distributed over a wider region than the narrow corridor along which Rio Grande rift structures are traditionally mapped at the surface, and they consequently demarcate regions of the lower crust where accommodation of bulk extensional strain has concentrated conductive phases. Our observations reveal limitations in existing models of Rio Grande rift activity and may reflect unrecognized spatiotemporal variations in rift system evolution globally.

Colorado↗
USGS2023

Magnetotelluric investigation of northern Harrat Rahat, Kingdom of Saudi Arabia

Volcanism within the harrats (Arabic for “volcanic field”) of the Kingdom of Saudi Arabia includes at least one historical eruption occurring close to the holy city of Al Madīnah in 1256 C.E. As part of a volcanic- and seismic-hazard assessment of northern Harrat Rahat, magnetotelluric (MT) data were collected to investigate the structural setting of the area, the presence or absence of melt within the crust, and the mantle-derived magmatic source. Collected MT data were modeled in both two dimensions, where anisotropy can be estimated, and three dimensions. Interpretation of the preferred resistivity model includes a shallow sediment-filled graben beneath northern Harrat Rahat lavas, a melt-free upper crust, and a region of decompression melting in the asthenosphere below 60–70 kilometers depth. Models in two dimensions image the lower crust as anisotropic, demonstrating that a series of elongate conductivity anomalies with a strike of N. 10° E. within the lower crust of the three-dimensional model are artifacts of inverting anisotropic data with an isotropic modeling algorithm. Careful examination of the resistivity models, in combination with regional geological and geophysical data, suggests an anisotropic lower crust that is free of large zones of melt. Azimuthal anisotropy in the lower crust extends well beyond the limits of Harrat Rahat volcanic rocks, with a conductive direction oriented N. 10° E. and an anisotropy factor of 2–5 between the most and least conductive directions. Enhanced conductivity is likely caused by interconnected grain-boundary graphite, where the direction of anisotropy reflects either frozen-in fabric from the Neoproterozoic stabilization of the Arabian Shield or ductile deformation driven by channelized asthenospheric flow coupled with a thin rigid mantle lid. Asthenospheric melt is interpreted to transect the crust largely through diking, with limited melt storage and short residence times within the crustal column.

northern Harrat Rahat↗
USGS2023

Magnetotelluric insights into the formation and reactivation of trans-crustal shear zones in Precambrian basement of the eastern U.S. Midcontinent

Three-dimensional inversion of regional long-period magnetotelluric (MT) data reveals the presence of two distinct sets of high-conductivity belts in the Precambrian basement of the eastern U.S. Midcontinent. One set, beneath Missouri, Illinois, Indiana, and western Ohio, is defined by northwest–southeast-oriented conductivity structures; the other set, beneath Kentucky, West Virginia, western Virginia, and eastern Ohio, includes structures that are generally oriented northeast–southwest. The northwest-trending belts occur mainly in Paleoproterozoic crust, and we suggest that their high conductivity values are due to graphite precipitated within trans-crustal shear zones from intrusion-related CO 2 -rich fluids. Our MT inversion results indicate that some of these structures dip steeply through the crust and intersect the Moho, which supports an interpretation that the shear zones originated as “leaky” transcurrent faults or transforms during the late Paleoproterozoic or the early Mesoproterozoic. The northeast-trending belts are associated with Grenvillian orogenesis and also potentially with Iapetan rifting, although further work is needed to verify the latter possibility. We interpret the different geographic positions of these two sets of conductivity belts as reflecting differences in origin and/or crustal rheology, with the northwest-trending belts largely confined to older, stable, pre-Grenville cratonic Laurentia, and the northeast-trending belts largely having formed in younger, weaker marginal crust. Notably, these high-conductivity zones spatially correlate with Midcontinent fault-and-fold zones that affect Phanerozoic strata. Stratigraphic evidence indicates that Midcontinent fault-and-fold zones were particularly active during Phanerozoic orogenic events, and some remain seismically active today, so the associated high-conductivity belts likely represent long-lived weaknesses that transect the crust.

Illinois, Indiana, Iowa, Kentucky, Michigan, Misso↗
USGS2023

Correction to “The March 1940 superstorm: Geoelectromagnetic hazards and impacts on American communication and power systems”

In Love et al. (2023), 5 estimated geopotentials on long lines for the March 1989 storm are incorrect. This was the result of incorrect calculation of line integrals. This error affects the two paragraphs of Section 11 and Figure 10 of that paper. The most significant impact of this error is that we overestimated the average ratio of the geopotentials measured the long lines during the March 1940 storm to the estimated geopotentials for the March 1989 storm. Instead of being a factor of 9 higher, the average ratio is only 9% higher. This means the inductive geoeffectiveness of the March 1940 storm was comparable to, or just slightly greater than, the geoeffectiveness of the March 1989 storm. Corrected paragraphs and a corrected figure are provided here.

Space Weather↗
USGS2023

Beyond conductive targets: Characterizing lithium-prospective lacustrine evaporite mineral systems of North America’s Basin and Range Province with regional-scale AEM

The Basin and Range province of North America hosts substantial lacustrine evaporite mineral systems prospective for lithium, a critical mineral currently listed for mineral resource assessment by the U.S. Geological Survey. Airborne electromagnetic (AEM) surveys are being conducted to support these assessments by identifying shallow clays and brines, as well as through improving the shallow subsurface geologic framework of the regional fluid flow system. In 2022-2023, three focus areas with proven lithium resources or considered highly prospective for lithium are being surveyed. Results from this effort can help to improve our understanding of the geologic conditions and geophysical signatures associated with known resource regions and benefit future lithium resource assessments by identifying regions with similar geophysical and geologic characteristics.

Arizona, California, Idaho, Nevada, Oregon, Utah↗
USGS2023

Importance of understanding bottom-up control when characterizing geothermal systems

Methods designed to identify favorable areas for geothermal resources have traditionally been focused on near-surface information, namely data that can be compiled into a 2D map. However, these methods fail to account for the third dimension: depth. As a result, they do not incorporate deep crustal and mantle features like heat sources. Geophysical methods with multi-scale capabilities, such as magnetotellurics (MT), provide tools to image deeper structures and bottom-up control on location of near surface hydrothermal systems in 3D. This study briefly demonstrates the advantage of understanding bottom-up control on hydrothermal systems to aid assessment and characterization. A regional 3D electrical resistivity model of the Great Basin is developed from MT data that images the near surface down to the mantle. From the 3D model, electrical conductance (depth integrated electrical conductivity) maps are created at logical depth intervals to identify anomalies. The conductance maps image discrete zones of high-conductance between 15-20 km depth indicative of fluid collection at the brittle-ductile transition; high-conductance zones around Moho depths (30-50 km) suggestive of partial melt; high-conductance zones in the upper mantle indicative of higher temperature and larger melt fraction, and low-conductance zones indicative of lithospheric material descending in the mantle. One anomalous zone of low-conductance in the mantle is under the north central Nevada, suggesting vertical mantle flow transport of heat to the crust related to sinking lithospheric material.

Conference Paper↗
USGS2022

Fluid transport and storage in the Cascadia forearc influenced by overriding plate lithology

Subduction of hydrated oceanic lithosphere can carry water deep into the Earth, with consequences for a range of tectonic and magmatic processes. Most of the fluid is released in the forearc where it plays a critical role in controlling the mechanical properties and seismic behaviour of the subduction megathrust. Here we present results from three-dimensional inversions of data from nearly 400 long-period magnetotelluric sites, including 64 offshore, to provide insights into the distribution of fluids in the forearc of the Cascadia subduction zone. We constrain the geometry of the electrically resistive Siletz terrane, a thickened section of oceanic crust accreted to North America in the Eocene, and the conductive accretionary complex underthrust along the margin. We find that fluids accumulate over timescales exceeding 1 My above the plate in metasedimentary units, while the mafic rocks of Siletzia remain dry. Fluid concentrations tend to peak at slab depths of 17.5 and 30 km, suggesting control by metamorphic processes, but also concentrate around the edges of Siletzia, suggesting that this mafic block is impermeable, with dehydration fluids escaping up-dip along the megathrust. Our results demonstrate that the lithology of the overriding crust can play a critical role in controlling fluid transport in a subduction zone.

British Columbia, California, Oregon, Washington↗
USGS2022

Geophysical extent of the Wyoming Province, western USA: Insights into ancient subduction and craton stability

A new 3-D resistivity model, estimated from inversion of magnetotelluric data, images crustal and upper-mantle structure of the Wyoming Province and adjacent areas. The Archean province is imaged as a coherent resistive domain, in sharp contrast to active tectonic domains of the western U.S. Prominent high-conductivity belts define the northern, eastern, and southern margins of the Wyoming Province and are interpreted as sutures marking the remnants of Paleoproterozoic orogens. The model results suggest the northern boundary of the Wyoming Province is located 150 km south of its traditional placement and adjacent to a composite orogen separating the Wyoming Province and Medicine Hat block. The eastern province boundary is clearly imaged along the Black Hills, whereas the western margin is obscured by Cenozoic extension and magmatism. An internal boundary within the Wyoming Province is interpreted to represent a Neoarchean suture; in stark contrast to Proterozoic sutures, though, it is not marked by a high-conductivity belt. This difference in conductivity is speculated to reflect changes in the subduction process through time. The absence of high-conductivity along Archean sutures appears to be global in nature and related to reduced continental freeboard in the Archean which limited continental weathering and the delivery of carbon-rich sediments to the seafloor. Although the entire Wyoming Province has been proposed to have undergone lithospheric modification that lessened its stability, the resistivity model suggests a thick lithospheric root remains in place except along its western margin. These results suggest that Archean cratons may be more resistant to lithospheric modification by influx of heat and fluids associated with extension and plumes than previously thought, and that metasomatism does not necessarily weaken the lithosphere and set a craton on the path to destruction.

Colorado, Montana, Utah, Wyoming↗