Geology topics

Anna Kelbert

Publications and source records attributed to Anna Kelbert.

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 intense geoelectric and geomagnetic field perturbations observed after an interplanetary magnetic field turning

Intense geoelectric fields during geomagnetic storms generate geomagnetically induced currents in power grids and other infrastructure, necessitating an understanding of their causes, for example, through coordinated space and ground observations. This study investigates localized intense geoelectric ( E ) and geomagnetic ( B ) field perturbations following an Interplanetary Magnetic Field (IMF) turning during a geomagnetic storm on 25 October 2011. Observations from EarthScope magnetotelluric sites in the upper Midwest United States revealed shorter period ( ∼ 1 min) ultra-low-frequency (ULF) waves superimposed on longer period ( ∼ 10 min) perturbations in both E and B fields. These sites, located at ∼ 19 hr magnetic local time and 56 −57 ° magnetic latitude, recorded large amplitude E and B perturbations. Ground-based all-sky imagers showed auroral brightening with sunward and poleward propagation, while upstream spacecraft linked the perturbations to an IMF turning and solar wind dynamic pressure impulse. The longer-period E and B field perturbations likely stem from localized ionospheric currents tied to substorm auroral activity post-IMF turning. The combination of ionospheric currents, ULF waves, and the Earth's varying conductivity produces intense geoelectric fields of ≥ 2 V/km in the upper Midwest. A comparison using input data and software compatible with the NOAA/USGS geoelectric field nowcast model revealed its limitations in capturing such events due to the temporal and spatial resolution of the underlying data. Using 1-s geomagnetic field data can improve geoelectric field models by capturing short-period and large spatial scale waves, although localized magnetic perturbations remain underestimated due to insufficient ground magnetometer density.

Space Weather↗

USGS2024

The solar cycle, geology, and geoelectric hazards for power grids

When sunspots are large and numerous, intense magnetic storms are likely to occur on the Earth. Magnetic storms can generate electric fields in the Earth, and these fields can, in turn, interfere with electric power transmission grids that are grounded at the Earth’s surface. Across the contiguous United States, geoelectric hazards are highest in the Upper Midwest and in the East. These regions correspond to geological structures that are electrically resistive, and they have, historically, experienced the most interference to electric power systems.

Fact Sheet↗

USGS2024

Hybrid CPU-GPU solution to regularized divergence-free curl-curl equations for electromagnetic inversion problems

The Curl-Curl equation is the foundation of time-harmonic electromagnetic (EM) problems in geophysics. The efficiency of its solution is key to EM simulations, accounting for over 95% of the computation cost in geophysical inversions for magnetotelluric or controlled-source EM problems. However, most published EM inversion codes are still central processing unit (CPU)-based and cannot utilize recent computational developments on the graphic processing units (GPUs). Based on a previously proposed divergence-free algorithm developed on CPUs, this study demonstrates the current limits of the CPU-based inversion procedure. To exploit the high throughput capability of GPUs, we propose a hybrid CPU-GPU framework to solve forward and adjoint problems required for EM inversions. The large sparse linear systems arising from the staggered-grid finite difference approximation of the Curl-Curl equation are solved with a mixed-precision Krylov subspace solver implemented on a GPU. The algorithm is then tested in EM forward and adjoint calculations, with real-world three-dimensional numerical examples. Test results show promising 30× kernel-level speed-ups over the conventional CPU algorithm. This approach may further take the complex frequency domain EM inversions onto the next, practical stage on small affordable GPU platforms.

Computers and Geosciences↗

USGS2024

Geoelectric field model validation in the southern California Edison system: Case study

Geomagnetic storms are a natural phenomenon that cause magnetic field variations at the surface of the Earth. These variations induce electrical current in natural and artificial conductors at and below the surface, resulting in geomagnetically induced currents (GIC) in power systems. The key to modeling GIC is to estimate the geoelectric field in the region of the power grid. The estimation of GIC follows by integrating the geoelectric field along the power transmission lines. In this study, we test the accuracy of the NOAA/USGS geoelectric field model during two geomagnetic storms in a regional power grid of the Southern California Edison system. Our comparison of modeled GIC with measurements shows reasonable agreement and gives confidence that the modeling is physically realistic. We also demonstrate the improved accuracy of modeled GIC relative to a simpler geoelectric field model based on regional 1D conductivity profiles.

Conference Paper↗

USGS2024

Aurora: An open-source Python implementation of the EMTF package for magnetotelluric data processing using MTH5 and mt-metadata

The Aurora software package robustly estimates single station and remote reference electromagnetic transfer functions (TFs) from magnetotelluric (MT) time series. Aurora is part of an open-source processing workflow that leverages the self-describing data container MTH5, which in turn leverages the general mt-metadata framework to manage metadata. These pre-existing packages simplify the processing by providing managed data structures, allowing for transfer functions to be generated with only a few lines of code. The processing depends on two inputs -- a table defining the data to use for TF estimation, and a JSON file specifying the processing parameters, both of which are generated automatically, and can be modified if desired. Output TFs are returned as mt_metadata objects, and can be exported to a variety of common formats for plotting, modeling and inversion.

Journal of Open Source Software↗

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

Estimating geomagnetically induced currents in southern Brazil using 3-D Earth resistivity model

Geomagnetically induced currents (GICs) result from the interaction of the time variation of ground magnetic field during a geomagnetic disturbance with the Earth's deep electrical resistivity structure. In this study, we simulate induced GICs in a hypothetical representation of a low-latitude power transmission network located mainly over the large Paleozoic Paraná basin (PB) in southern Brazil. Two intense geomagnetic storms in June and December 2015 are chosen and geoelectric fields are calculated by convolving a three-dimensional (3-D) Earth resistivity model with recorded geomagnetic variations. The dB / dt proxy often used to characterize GIC activity fails during the June storm mainly due to the relationship of the instantaneous geoelectric field to previous magnetic field values. Precise resistances of network components are unknown, so assumptions are made for calculating GIC flows from the derived geoelectric field. The largest GICs are modeled in regions of low conductance in the 3-D resistivity model, concentrated in an isolated substation at the northern edge of the network and in a cluster of substations in its central part where the east-west (E-W) oriented transmission lines coincide with the orientation of the instantaneous geoelectric field. The maximum magnitude of the modeled GIC was obtained during the main phase of the June storm, modeled at a northern substation, while the lowest magnitudes were found over prominent crustal anomalies along the PB axis and bordering the continental margin. The simulation results will be used to prospect the optimal substations for installation of GIC monitoring equipment.

Space Weather↗

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

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

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

MTH5: An archive and exchangeable data format for magnetotelluric time series data

Magnetotellurics (MT) is a passive electromagnetic geophysical method that measures variations in subsurface electrical resistivity. MT data are collected in the time domain and processed in the frequency domain to produce estimates of a transfer function representing the Earth’s electrical structure. Unfortunately, the MT community lacks metadata and data standards for time series data. As the community grows and findability, accessibility, interoperability, and reuse of digital assets (FAIR) data principles are enforced by government and funding agencies, a standard is needed for time series data. Presented here is a hierarchical data format (MTH5) that is logically formatted to how MT data are collected. Open-source Python packages are also described to read, write, and manipulate MTH5 files. These include a package to deal with metadata ( mt_metadata ) based on standards developed by the Working Group for Magnetotelluric Data Handling and Software assembled by the Incorporated Research Institutions for Seismology (IRIS), and mth5 : a package to interact with MTH5 files that uses mt_metadata . Example code and workflows are presented.

Computers & Geosciences↗

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↗

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↗