Geology topics

Paul A. Bedrosian

Publications and source records attributed to Paul A. Bedrosian.

At least 55 records · Page 3Linked to original sources

USGS2020

High-resolution mapping of the freshwater-brine interface using deterministic and Bayesian inversion of airborne electromagnetic data at Paradox Valley, USA

Salt loads in the Colorado River Basin are a primary water quality concern. Natural groundwater brine discharge to the Dolores River where it passes through the collapsed salt anticline of the Paradox Valley in western Colorado is a significant source of salt to the Colorado River. An airborne electromagnetic survey of Paradox Valley has provided insights into the 3D distribution of brine in the surficial aquifer. A combination of stochastic and deterministic resistivity inversions were used to interpret the top of the freshwater-brine interface and to qualitatively describe the vertical salinity gradients across the interface. Low-resistivity regions indicative of brine occur near the land surface where brine discharges to the Dolores River and increase in depth several kilometers up-gradient along the axis of the valley. The most conductive parts of the brine plume are found in the areas below and adjacent to the river, suggesting that the brine becomes shallower and more concentrated as it reaches its natural discharge location. A significant freshwater lens overlying the brine west of the Dolores River is spatially correlated to the intermittent West Paradox Creek and agricultural irrigation. Below this lens, the transition from freshwater to brine appears to occur abruptly over a few meters and correlates to available well information. However, away from these regions and particularly with distance from the river, the freshwater-brine interface appears to be more diffuse.

Colorado↗

USGS2020

A 100-year geoelectric hazard analysis for the U.S. high-voltage power grid

A once-per-century geoelectric hazard map is created for the United States high-voltage power grid. A statistical extrapolation from 31 years of magnetic field measurements is made by identifying 84 geomagnetic storms with the Kp and Dst indices. Data from 24 geomagnetic observatories, 1079 magnetotelluric survey sites, and 17,258 transmission lines are utilized to perform a geoelectric hazard analysis with the most comprehensive data publicly available. With this data we estimate once-per-century geoelectric fields at the magnetotelluric survey sites and calculate the theoretical voltages within transmission lines in the United States power grid. Once-per-century geoelectric field strengths span more than three orders of magnitude from a minimum of 0.02 V/km at a site in Idaho to a maximum of 26.8 V/km at a site in Maine, with nearly 30% of the surveyed land area ex- ceeding 1 V/km. We show the influence that geoelectric field polarization has on geoelec- tric hazards when viewed on a power transmission network. The calculated transmission line voltages can exceed 1000 V in some transmission lines. Four regions in the United States with particularly notable geoelectric hazards are identified and discussed: the East coast, Pacific Northwest, upper Midwest, and the Denver metropolitan area.

Space Weather↗

USGS2019

Crustal magmatism and anisotropy beneath the Arabian Shield - A cautionary tale

Volcanism in Saudi Arabia includes a historic eruption close to the holy city of Al Madinah. As part of a volcanic hazard assessment of this area, magnetotelluric (MT) data were collected to investigate the structural setting, the distribution of melt within the crust, and the mantle source of volcanism. Interpretation of a new 3‐D resistivity model includes a shallow graben beneath thin lava fields (Harrats), a melt‐free upper crust, and decompression melting in the asthenosphere below thin lithosphere. Within the lower crust the model images elongate conductivity anomalies, one of which was attributed in a previous MT study to melt. The regional MT data, combined with perspective from geology and geophysical modeling, suggest the lower crust is anisotropic with no interconnected melt zones. These divergent interpretations have distinct hazard implications and highlight the importance of large survey aperture and anisotropic modeling to MT studies of volcanic regions. Lower‐crustal anisotropy extends beyond the Harrat, with the most conductive direction oriented N10°E and a factor of 3–5, determined from 2‐D anisotropic inversion, between the most and least conductive directions. The enhanced conductivity is likely due to interconnected grain boundary graphite, while the anisotropy direction reflects either frozen‐in fabric from Neoproterozoic stabilization of the Arabian Shield or modern ductile deformation driven by channelized asthenospheric flow coupled through a thin rigid mantle lid. Asthenospheric melt is interpreted to transect the crust primarily through diking, with limited melt storage and short residence times in the crust.

Arabian Shield↗

USGS2019

Empirical estimation of natural geoelectric hazards

Geoelectric field time series can be estimated by convolving estimates of Earth‐surface impedance, such as those obtained from magnetotelluric survey measurements, with historical records of geomagnetic variation obtained at magnetic observatories. This straightforward procedure permits the mapping of geoelectric field variation during magnetic storms. Statistical analysis of the time series allows extrapolation to extreme‐value amplitudes, such as might be realized during an intense magnetic storm in the future. The development of these products is illustrated for the Mid‐Atlantic United States, using impedances obtained from EarthScope survey data and geomagnetic variation records obtained at the Fredericksburg observatory operated by the U.S. Geological Survey. For this region, 100‐year geoelectric exceedance amplitudes have a range of almost three orders of magnitude (from 0.04 V/km at a site in southern Pennsylvania to 24.29 V/km at a site in central Virginia), and they have significant geographic granularity, which is due to site‐to‐site differences in surface impedance (and subsurface electrical conductivity structure). Maps of 100‐year exceedance amplitudes resemble those of geoelectric amplitudes for the March 1989 magnetic storm, and, in that sense, the March 1989 storm resembles what might be loosely called a “100‐year” event.

Book chapter↗

USGS2019

The first 3D conductivity model of the contiguous US: Reflections on geologic structure and application to induction hazards

Estimation of ground level geoelectric fields has been identified by the National Space Weather Action Plan as a key component of assessment and mitigation of space weather impacts on critical infrastructure. Estimates of spatially and temporally variable electric fields are used to generate statistically based hazard maps and show promise toward monitoring and responding to geomagnetic disturbances in near real‐time. One approach to geoelectric field estimation is to employ three‐dimensional (3D) Earth conductivity models. These data‐constrained conductivity models are the results of regional magnetotelluric inversions based primarily on NSF’s Earthscope USArray impedances, which to date cover ~60% of the contiguous United States. Here, we present the first‐ever composite conductivity model of the contiguous United States and describe its compilation from 3D regional conductivity models, a global mantle conductivity model, offshore bathymetry, and sediment thickness data. We discuss structures within the conductivity model and how they relate to the complex geologic tapestry of the continent. Finally, we discuss the utility of this synthesis model for estimation and mitigation of geomagnetically induced currents.

Book chapter↗

USGS2019

Lithospheric signature of late Cenozoic extension in electrical resistivity structure of the Rio Grande rift, New Mexico, USA

We present electrical resistivity models of the crust and upper mantle from two‐dimensional (2‐D) inversion of magnetotelluric (MT) data collected in the Rio Grande rift, New Mexico, USA. Previous geophysical studies of the lithosphere beneath the rift identified a low‐velocity zone several hundred kilometers wide, suggesting that the upper mantle is characterized by a very broad zone of modified lithosphere. In contrast, the surface expression of the rift (e.g., high‐angle normal faults and synrift sedimentary units) is confined to a narrow region a few tens of kilometers wide about the rift axis. MT data are uniquely suited to probing the depths of the lithosphere that fill the gap between surface geology and body wave seismic tomography, namely the middle to lower crust and uppermost mantle. We model the electrical resistivity structure of the lithosphere along two east‐west trending profiles straddling the rift axis at the latitudes of 36.2 and 32.0°N. We present results from both isotropic and anisotropic 2‐D inversions of MT data along these profiles, with a strong preference for the latter in our interpretation. A key feature of the anisotropic resistivity modeling is a broad (~200‐km wide) zone of enhanced conductivity (<20 Ωm) in the middle to lower crust imaged beneath both profiles. We attribute this lower crustal conductor to the accumulation of free saline fluids and partial melt, a direct result of magmatic activity along the rift. High‐conductivity anomalies in the midcrust and upper mantle are interpreted as fault zone alteration and partial melt, respectively.

New Mexico↗

USGS2018

Crustal inheritance and a top-down control on arc magmatism at Mount St Helens

In a subduction zone, the volcanic arc marks the location where magma, generated via flux melting in the mantle wedge, migrates through the crust and erupts. While the location of deep magma broadly defines the arc position, here we argue that crustal structures, identified in geophysical data from the Washington Cascades magmatic arc, are equally important in controlling magma ascent and defining the spatial distribution and compositional variability of erupted material. As imaged by a three-dimensional resistivity model, a broad lower-crustal mush zone containing 3–10% interconnected melt underlies this segment of the arc, interpreted to episodically feed upper-crustal magmatic systems and drive eruptions. Mount St Helens is fed by melt channelled around a mid-Tertiary batholith also imaged in the resistivity model and supported by potential–field data. Regionally, volcanism and seismicity are almost exclusive of the batholith, while at Mount St Helens, along its margin, the ascent of viscous felsic melt is enabled by deep-seated metasedimentary rocks. Both the anomalous forearc location and composition of St Helens magmas are products of this zone of localized extension along the batholith margin. This work is a compelling example of inherited structural control on local stress state and magmatism.

Washington↗

USGS2018

Geoelectric hazard maps for the Pacific Northwest

Maps of extreme value, horizontal component geoelectric field amplitude are constructed for the Pacific Northwest United States (and parts of neighboring Canada). Multidecade long geoelectric field time series are calculated by convolving Earth surface impedance tensors from 71 discrete magnetotelluric survey sites across the region with historical 1‐min (2‐min Nyquist) geomagnetic variation time series obtained from two nearby observatories. After fitting statistical models to 1‐min geoelectric amplitudes realized during magnetic storms, extrapolations are made to estimate threshold amplitudes that are only exceeded, on average, once per century. One hundred‐year geoelectric exceedance amplitudes range from 0.06 V/km at a survey site in western Washington State to 9.47 V/km at a site in southeast British Columbia; 100‐year geoelectric exceedance amplitudes equal 7.10 V/km at a site north of Seattle and 2.28 V/km at a site north of Portland. Systematic and random errors are estimated to be less than 20%, much less than site‐to‐site differences in geoelectric amplitude that arise from site‐to‐site differences in surface impedance. Maps of 100‐year exceedance amplitudes are compared with the peak geoelectric amplitudes realized during the March 1989 magnetic superstorm; it is noted that some storms of relatively modest intensity can generate localized geoelectric fields of relatively high amplitude. The geography of geoelectric hazard across the Pacific Northwest is closely related to known geologic and tectonic structures.

Pacific Northwest↗

USGS2018

On the feasibility of real-time mapping of the geoelectric field across North America

A review is given of the present feasibility for accurately mapping geoelectric fields across North America in near-realtime by modeling geomagnetic monitoring and magnetotelluric survey data. Should this capability be successfully developed, it could inform utility companies of magnetic-storm interference on electric-power-grid systems. That real-time mapping of geoelectric fields is a challenge is reflective of (1) the spatiotemporal complexity of geomagnetic variation, especially during magnetic storms, (2) the sparse distribution of ground-based geomagnetic monitoring stations that report data in realtime, (3) the spatial complexity of three-dimensional solid-Earth impedance, and (4) the geographically incomplete state of continental-scale magnetotelluric surveys.

Open-File Report↗

USGS2018

Combining multiphase groundwater flow and slope stability models to assess stratovolcano flank collapse in the Cascade Range

Hydrothermal alteration can create low‐permeability zones, potentially resulting in elevated pore‐fluid pressures, within a volcanic edifice. Strength reduction by rock alteration and high pore‐fluid pressures have been suggested as a mechanism for edifice flank instability. Here we combine numerical models of multiphase heat transport and groundwater flow with a slope‐stability code that incorporates three‐dimensional distributions of strength and pore‐water pressure to address the following questions: (1) What permeability distributions and contrasts produce elevated pore‐fluid pressures in a stratovolcano? (2) What are the effects of these elevated pressures on flank stability? (3) Finally, what are the effects of magma intrusion on potential flank failure in an edifice? Simulation results show that under a range of plausible parameters, water tables in a stratovolcano can be elevated or perched. These elevated water tables result in universally lower stability (lower factor of safety) compared with equivalent dry edifices, indicating a higher likelihood of flank collapse. Low‐permeability (<1 × 10 −17 m 2 ) layers such as altered pyroclastic deposits or breccias can result in locally saturated regions (perched water) and lower factors of safety near the ground surface but may actually reduce liquid water saturation and pore pressures in the core of the edifice and thus may favor small, shallow collapses over larger, deeper collapses. Magma intrusion into the base of the edifice increases pore‐fluid pressures and decreases the factor of safety. However, the shear strength of edifice rocks also exerts a significant control on stability, so both mechanical properties and pore‐fluid pressures are important for stability assessments.

Washington↗

USGS2018

Geoelectric hazard maps for the Mid-Atlantic United States: 100 year extreme values and the 1989 magnetic storm

Maps of extreme value geoelectric field amplitude are constructed for the Mid‐Atlantic United States, a region with high population density and critically important power grid infrastructure. Geoelectric field time series for the years 1983–2014 are estimated by convolving Earth surface impedances obtained from 61 magnetotelluric survey sites across the Mid‐Atlantic with historical 1 min (2 min Nyquist) measurements of geomagnetic variation obtained from a nearby observatory. Statistical models are fitted to the maximum geoelectric amplitudes occurring during magnetic storms, and extrapolations made to estimate threshold amplitudes only exceeded, on average, once per century. For the Mid‐Atlantic region, 100 year geoelectric exceedance amplitudes have a range of almost 3 orders of magnitude (from 0.04 V/km at a site in southern Pennsylvania to 24.29 V/km at a site in central Virginia), and they have significant geographic granularity, all of which is due to site‐to‐site differences in magnetotelluric impedance. Maps of these 100 year exceedance amplitudes resemble those of the estimated geoelectric amplitudes attained during the March 1989 magnetic storm, and, in that sense, the March 1989 storm resembles what might be loosely called a “100 year” event. The geoelectric hazard maps reported here stand in stark contrast with the 100 year geoelectric benchmarks developed for the North American Electric Reliability Corporation.

Geophysical Research Letters↗

USGS2018

Extreme-event geoelectric hazard maps: Chapter 9

Maps of geoelectric amplitude covering about half the continental United States are presented that will be exceeded, on average, once per century in response to an extreme-intensity geomagnetic disturbance. These maps are constructed using an empirical parameterization of induction: convolving latitude-dependent statistical maps of extreme-value geomagnetic disturbances, obtained from decades of 1-minute magnetic observatory data, with local estimates of Earth-surface impedance obtained at discrete geographic sites from magnetotelluric surveys. Geoelectric amplitudes are estimated for geomagnetic waveforms having a 240-s (and 1200-s) sinusoidal period and amplitudes over 10 min (1 h) that exceed a once-per-century threshold. As a result of the combination of geographic differences in geomagnetic variation and Earth-surface impedance, once-per-century geoelectric amplitudes span more than two orders of magnitude and are a highly granular function of location. Specifically for north-south 240-s induction, once-per-century geoelectric amplitudes across large parts of the United States have a median value of 0.34 V/km; for east-west variation, they have a median value of 0.23 V/km. In Northern Minnesota, amplitudes exceed 14.00 V/km for north-south geomagnetic variation (23.34 V/km for east-west variation), while just over 100 km away, amplitudes are only 0.08 V/km (0.02 V/km). At some sites in the northern-central United States, once-per-century geoelectric amplitudes exceed the 2 V/km realized in Québec during the March 1989 storm.

Book chapter↗

USGS2018

Extreme‐value geoelectric amplitude and polarization across the northeast United States

Maps are presented of extreme‐value geoelectric field amplitude and horizontal polarization for the Northeast United States. These maps are derived from geoelectric time series calculated for sites across the Northeast by frequency‐domain multiplication (time‐domain convolution) of 172 magnetotelluric impedance tensors, acquired during a survey, with decades‐long, 1‐min resolution time series of geomagnetic variation, acquired at three magnetic observatories. The maps show that, during intense magnetic storms, high geoelectric amplitude hazards are realized across electrically resistive, igneous and metamorphic rock of the Appalachian Mountains and the New England Highlands, while low geoelectric hazards are realized across electrically conductive, sedimentary rock of the Appalachian Plateau and the Mid‐Atlantic Coastal Plain. From statistical extrapolation, once‐per‐century (100‐year) geoelectric amplitudes are highest at a site in Virginia at 25.44 V/km (followed by a site in Maine at 21.75 V/km and a site in Connecticut at 19.39 V/km); 100‐year geoelectric amplitude exceeds 10 V/km at 12 sites across the northeast; geoelectric amplitude is lowest at a site in Virginia at 0.05 V/km. Average errors for these values are estimated to be about 38%, or much less than the more than 2 orders of magnitude range seen in geoelectric amplitudes from one survey site to another across the northeast. It is noteworthy that geoelectric fields tend to be most (least) polarized at locations with high (low) geoelectric hazard. Furthermore, geoelectric fields over the Appalachians tend to be polarized southeast‐to‐northwest, or generally in a direction orthogonal to the southwest‐to‐northeast geological strike. Results reported here inform utility companies in projects for evaluating and managing the response of power grid systems to the deleterious effects of geomagnetic disturbance.

Connecticut, Delaware, Maine, Maryland, Massachuse↗

USGS2017

Magnetotelluric imaging of lower crustal melt and lithospheric hydration in the Rocky Mountain Front transition zone, Colorado, USA

We present an electrical resistivity model of the crust and upper mantle from two‐dimensional (2‐D) anisotropic inversion of magnetotelluric data collected along a 450 km transect of the Rio Grande rift, southern Rocky Mountains, and High Plains in Colorado, USA. Our model provides a window into the modern‐day lithosphere beneath the Rocky Mountain Front to depths in excess of 150 km. Two key features of the 2‐D resistivity model are (1) a broad zone (~200 km wide) of enhanced electrical conductivity (<20 Ωm) in the midcrust to lower crust that is centered beneath the highest elevations of the southern Rocky Mountains and (2) hydrated lithospheric mantle beneath the Great Plains with water content in excess of 100 ppm. We interpret the high conductivity region of the lower crust as a zone of partially molten basalt and associated deep‐crustal fluids that is the result of recent (less than 10 Ma) tectonic activity in the region. The recent supply of volatiles and/or heat to the base of the crust in the late Cenozoic implies that modern‐day tectonic activity in the western United States extends to at least the western margin of the Great Plains. The transition from conductive to resistive upper mantle is caused by a gradient in lithospheric modification, likely including hydration of nominally anhydrous minerals, with maximum hydration occurring beneath the Rocky Mountain Front. This lithospheric “hydration front” has implications for the tectonic evolution of the continental interior and the mechanisms by which water infiltrates the lithosphere.

Colorado↗

USGS2017

Helicopter electromagnetic data map ice thickness at Mount Adams and Mount Baker, Washington, USA

Ice-thickness measurements critical for flood and mudflow hazard studies are very sparse on Cascade Range (North America) volcanoes. Helicopter electromagnetic (HEM) data collected to detect hydrothermal alteration are used to determine ice thickness over portions of Mount Baker and Mount Adams volcanoes. A laterally continuous inversion method provides good estimates of ice <100 m thick over water-saturated and altered regions where the resistivity of the basement is <200 Ωm. For areas with ice overlying fresh, resistive rocks with small resistivity contrasts between ice and rock, ice thickness is not well resolved. The ice thicknesses derived from HEM data are consistent with the previous drillhole data from Mount Adams and radar data from both volcanoes, with mean thicknesses of 57 m for Mount Adams and 68 m for Mount Baker. The thickest ice on Mount Baker rests on the gentle lower slopes whereas the thickest ice at Mount Adams lies on the flat summit. Ice volume calculations suggest that Mount Baker contains ∽710 × 10 6 m 3 of ice in the HEM survey area, with a crude estimate of ∽1800 × 10 6 m 3 for the entire volcano. Ice volume on Mount Adams is 65 × 10 6 m 3 in parts of the HEM survey area and ∽200 × 10 6 m 3 overall.

Washington↗

USGS2017

Down to Earth with an electric hazard from space

In reaching across traditional disciplinary boundaries, solid-Earth geophysicists and space physicists are forging new collaborations to map magnetic-storm hazards for electric-power grids. Future progress in evaluation storm time geoelectric hazards will come primarily through monitoring, surveys, and modeling of related data.

Space Weather↗

USGS2017

Mineral potential mapping in an accreted island-arc setting using aeromagnetic data: An example from southwest Alaska

The distribution of volcanogenic massive sulfide (VMS), porphyry-epithermal, Alaska-type ultramafic-mafic complexes, intrusion-related Au, and granitoid Sn-W ore deposits in southwest Alaska supports current metallogenic models linking the formation of these deposit types to the emplacement of different suites of igneous rocks during the evolution of this convergent plate margin. Regional-scale aeromagnetic data provide a continuous set of observations over the deposits and show contrasting patterns over the igneous rock suites hosting the various deposit types. Combined with surface geologic data and regional metallogenic constraints, aeromagnetic data—filtered to enhance the anomalous magnetic field and map magnetic domains—were used to produce a mineral potential map across this accreted island-arc setting. The reduced-to-pole, upward continuation, and total horizontal gradient transform maps show anomalies that could represent porphyry-epithermal deposits within the intraoceanic- and continental-arc terranes. The tilt derivative transform highlights lineaments within the back arc that may represent zones with potential for VMS deposits. The truncations of tilt derivative lineaments outline a major magnetic domain boundary between the back-arc and craton margin, which is prospective for granitoid Sn-W deposits. Annular tilt derivative highs outline granitoids that could be associated with intrusion-related Au deposits within the craton margin. Shallow, magnetite-rich Alaska-type ultramafic-mafic complexes are mapped by their short-wavelength, high-amplitude anomalies. Successful mineral potential mapping across southwestern Alaska as performed in the present study suggests that filtered aeromagnetic data can be effectively used in mineral exploration in convergent continental margin settings.

Alaska↗

USGS2017

Geoelectric monitoring at the Boulder magnetic observatory

Despite its importance to a range of applied and fundamental studies, and obvious parallels to a robust network of magnetic-field observatories, long-term geoelectric field monitoring is rarely performed. The installation of a new geoelectric monitoring system at the Boulder magnetic observatory of the US Geological Survey is summarized. Data from the system are expected, among other things, to be used for testing and validating algorithms for mapping North American geoelectric fields. An example time series of recorded electric and magnetic fields during a modest magnetic storm is presented. Based on our experience, we additionally present operational aspects of a successful geoelectric field monitoring system.

Colorado↗