Search USGS⌕ Search

SEARCH · Search USGS

Results for “Earth Surface Processes and Landforms”

Search indexed USGS publications on groundwater, aquifers, geologic maps, mineral resources and earthquakes. Explore source records by subject and place.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

1,231 records · Page 2Linked to original sources
USGS2025

From hydrated silica to quartz: Potential hydrothermal precipitates found in Jezero crater, Mars

On Earth, silica-rich phases from opal to quartz are important indicators and tracers of geological processes. Hydrated silica, such as opal, is a particularly good matrix for the preservation of molecular and macroscopic biosignatures. Cherts, a type of silica-dominated rocks, provide a unique archive of ancient terrestrial life while quartz is the emblematic mineral of the Earth's continental crust. On Mars, hydrated silica has been detected in several locations based on remote sensing and rover-based studies. In the present article we report on the detection of cobbles made of hydrated silica (opal or chalcedony), as well as well-crystallized quartz. These detections were made with the SuperCam instrument onboard Perseverance (Mars 2020 mission), using a combination of LIBS, infrared and Raman spectroscopy. Quartz-dominated stones are detected unambiguously for the first time on the Martian surface, and based on grain size and crystallinity are proposed to be of hydrothermal origin. Although these rocks were all found as float, we propose that these detections are part of a common hydrothermal system, and represent different depths / temperatures of precipitation. This attests that hydrothermal processes were active in and around Jezero crater, possibly triggered by the Jezero crater-forming impact. These silica-rich rocks, in particular opaline silica, are very promising targets for sampling and return to Earth given their high biosignature preservation potential.

Earth and Planetary Science Letters↗
USGS2025

Rainfall thresholds for postfire debris-flow initiation vary with short-duration rainfall climatology

The size, frequency, and geographic scope of severe wildfires are expanding across the globe, including in the Western United States. Recently burned steeplands have an increased likelihood of debris flows, which pose hazards to downstream communities. The conditions for postfire debris-flow initiation are commonly expressed as rainfall intensity-duration thresholds, which can be estimated given sufficient observational history. However, the spread of wildfire across diverse climates poses a challenge for accurate threshold prediction in areas with limited observations. Studies of mass-movement processes in unburned areas indicate that thresholds vary with local climate, such that higher rainfall rates are required for initiation in climates characterized by frequent intense rainfall. Here, we use three independent methods to test whether initiation of postfire runoff-generated debris flows across the Western United States varies similarly with climate. Through the compilation of observed thresholds at various fires, analysis of the spatial density of observed debris flows, and quantification of feature importance at different spatial scales, we show that postfire debris-flow initiation thresholds vary systematically with short-duration rainfall-intensity climatology. The predictive power of climatological data sets that are readily available before a fire occurs offers a much-needed tool for hazard management in regions that are facing increased wildfire activity, have sparse observational history, and/or have limited resources for field-based hazard assessment. Furthermore, if the observed variation in thresholds reflects long-term adjustment of the landscape to local climate, rapid shifts in rainfall intensity related to climate change will likely induce spatially variable shifts in postfire debris-flow likelihood.

Arizona, California, Colorado, Nevada, New Mexico,↗
USGS2025

Memory and jamming in fault zone sediments

Many subsurface processes involve transitions in granular material states, from arrested to creeping to flowing. Experiments and frameworks for idealized systems reveal that granular fabrics develop during shearing, co-evolve with applied stress, and govern such transitions. We use microtomography to test whether fabrics at two San Andreas fault sites reflect slip history and whether idealized frameworks extend to nature. Near-surface sediments within the fault zone transition between deformation patterns over the seismic cycle, including bulk/localized grain re-arrangements, individual grain fracturing, and localized zones of fracturing. Aseismic and co-seismic shearing produce distinct preferred grain orientations. Co-seismic fabrics can be preserved after centuries of aseismic strain, aseismic fabrics may be overprinted, and grain size and coordination number influence the fabrics. Idealized frameworks, namely anisotropic critical state theory, frictional jamming, and material memory, can explain our observations, and fault zone sediments likely undergo cycles of memory creation and erasure that influence rigidity spatiotemporally.

California↗
USGS2024

Getting to the surface of the problem: A dynamic rupture benchmark for shallowly dipping faults near Earth’s surface

The Statewide California Earthquake Center (SCEC) Dynamic Rupture Group workshop “Getting to the Surface of the Problem: A Dynamic Rupture Benchmark for Shallowly-Dipping Faults Near Earth’s Surface” was convened virtually, in an online meeting on November 4, 2024. 52 people participated, including two who sent pre-recorded lightning talks earlier, but were unable to join us that day due to the international time-zone difference. One-half of our participants were students or postdocs, and the group included scientists from 11 countries (USA, Canada, China, Czech Republic, France, Germany, Iceland, Japan, New Zealand, Saudi Arabia, Switzerland).

Final Report↗
USGS2026

Continental-scale prediction of hydrologic signatures and processes

Understanding how dominant hydrologic processes and their drivers vary across diverse continental-scale landscapes is critical for hydrologic modeling and water management applications. Our research addresses this question by synthesizing large-sample watershed datasets, Caravan and GAGES-II, and developing random forest models to identify patterns in hydrologic function. We assessed dominant processes by examining hydrologic signatures – summary indicators of watershed function derived from hydroclimatic time series and random forest models across 14 146 gauged United States watersheds. The results reveal clear continental-scale gradients in hydrologic processes, including baseflow, overland flow, storage, and water balance losses. Our map of dominant processes highlights, for example, the transition from baseflow to fast responses and back to baseflow along the elevation gradient from the Appalachian spine, through the Piedmont, to the Eastern Coastal Plain; a distinct outer ring around the Great Lakes region; and sharp contrasts between coastal and inland processes in the West. Variable importance analysis from random forest models show that processes in the western U.S. are primarily controlled by climate, whereas in the eastern U.S., soil, geology, and topography play larger roles, with distinct human influences apparent in urban areas. Our approach of estimating dominant processes and their drivers facilitates extending process knowledge from research watersheds to the continental scale, assessing current hydrological understanding, and evaluating hydrological model structures.

conterminous United States↗
USGS2025

State shifts in the deep Critical Zone drive landscape evolution in volcanic terrains

Understanding the near-surface environment where atmospheric and solid earth processes interact, often termed the “Critical Zone,” is important for assessing resources and building resilient societies. Here, we examine a volcanic landscape in the Oregon Cascade Range, an understudied Critical Zone setting that is host to major regional water resources, pervasive silicate weathering, and significant geohazards. We leverage a bedrock age chronosequence to show that the volcanic Critical Zone undergoes a structural shift, from depth extents of >1 km to meters, over timescales of ~1 My. We map an active groundwater volume comparable to major continental lakes, stored at the Cascade Range crest. This state shift makes volcanic landscape evolution a unique probe of deep coupling between Earth systems.

PNAS↗
USGS2026

Forecasting volcanic activity in Germany—A multi-criteria approach

Igneous activity, including shallow intrusions and volcanism, has the potential to disrupt underground critical infrastructure. Notably, future underground infrastructure projects like high-level radioactive waste repositories must be sited in areas of extremely low disruption probability by igneous activity. In Germany, according to the Repository Site Selection Act of 2017 (Standortauswahlgesetz, or StandAG), areas in which Quaternary volcanism is either present or future volcanic activity is expected within the next 1 million years (m.y.) must be excluded from the site selection process. Although the locations of regions with Quaternary volcanism are reasonably well known in Germany, forecasting potential igneous activity at intraplate volcanic fields is challenging, as many processes and their interactions control the spatial distribution of volcanic centers. Here, a semi-quantitative, multi-criteria approach is proposed for a regional evaluation of the relative potential of future igneous activity in Germany. A variety of geoscientific indicators are used, including seismic anomalies in Earth’s mantle, gravity data, tectonic activity, sutures, ground motion, earthquakes, mantle degassing centers, and geochronological data of volcanic rocks. The indicators describe the sequence of processes from potential melt generation in Earth’s mantle, through ascent and accumulation of melt within the lithosphere, to eruption at Earth’s surface. In total, 15 out of 30 proposed geoscientific indicators are selected and quantified using 20 total assigned parameters. Defined threshold values are used to spatially delimit relevant parameter properties to focus on areas with higher potential of future magmatic activity. To consider uncertainties of parameters and their underlying processes, which are usually more spatially extensive below ground, buffer zones are defined in which values of relevance decrease with increasing distance from the initial lateral shape of a parameter. Normalized parameters are combined into an index, whose spatial value distribution is used to differentiate the relative potential of future igneous activity (within the next 1 m.y.). The sensitivity of the results is shown by varying the weighting factors for the relevant parameters in country-wide index maps. Thereby, profiles illustrate the distribution of the resulting index values and respective index fractions of various parameters. Different index maps for the relative potential of future igneous activity are presented and can be used for hazard assessments.

Professional Paper↗
USGS2026

Understanding the resource potential of natural hydrogen on Earth: Scientific gaps, uncertainties and recommendations

A comprehensive scientific research roadmap is essential to bridge knowledge gaps and deepen the understanding of key geological, geochemical, and geophysical aspects of natural hydrogen (H 2 ) as a potential new energy resource. This paper reviews major scientific uncertainties on natural H 2 , suggesting research priorities, as a guide for defining exploration strategies, techniques, and data interpretation. The uncertainties concern all phases of the natural H 2 cycle, from generation (source rocks) through migration (advection and diffusion) and accumulation (reservoir and cap rocks) to the application and interpretation of subsurface and surface geochemical and geophysical exploration techniques. Understanding H 2 sources and generation rates (the amount of H 2 generated by a given volume of rock over time) is crucial for determining whether a geological H 2 system operates as a short-term dynamic system with rapid H 2 production and release, or as a conventional gas system with long-term accumulations, analogous to petroleum reservoirs. Preliminary estimates for serpentinisation, radiolysis, and organic matter degradation suggest that H 2 generation is not inherently fast, especially for non-hydrothermal continental systems (crystalline basement of shields, ophiolites, peridotite massifs, sedimentary basins), and long-term accumulations, like those of fossil natural gas systems, represent the most likely scenario. The mechanisms of H 2 migration through geological formations require application of fundamental principles of fluid-flow physics, distinguishing advection and diffusion, as well as their forms (from gas-phase, bubble flows to aqueous solutions). Additional studies of H 2 accumulation and retention in subsurface reservoirs could improve understanding of mechanisms of H 2 migration by focusing on the rock fluid-bearing properties and the factors affecting H 2 preservation, such as the presence of cap rocks impermeable to H 2 , pressure conditions, residence times, and microbial or abiotic consumption. Advanced techniques, including reservoir modelling, flow simulations, 3D imaging (micro-CT) of H 2 -bearing rocks, and extraction and analysis of gas occluded in rocks, can provide insights into the stability and potential recoverability of H 2 accumulations. The interpretation of surface exploration techniques, including gas geochemistry, geophysics, and remote sensing, long employed in mineral and energy resource exploration, is now being adapted for natural H 2 studies, but challenges remain in the data interpretation. Distinguishing H 2 seepage due to geological degassing from H 2 produced near the surface by modern microbial processes or artificial sources, such as hammering or drilling for soil-gas sampling, drilling into aquifers, and corrosion in boreholes, is an essential step in exploration. The simple detection of H 2 in soils, even in morphological structures like sub-circular depressions or “fairy circles”, cannot be cursorily interpreted as a signal of natural H 2 seepage from a deep source. A holistic geochemical approach, including isotopic analyses of gases associated with H 2 , is recommended to distinguish among the variety of possible H 2 origins. Observations of H 2 in wells should be interrogated to rule out possible artifacts such as corrosion and drill bit metamorphism. The integration of multiple geophysical methods, including seismic, gravimetric, magnetic, and electro-magnetic surveys, is recommended to mitigate interpretation ambiguities regarding the structure of a subsurface H 2 system (source and reservoir rocks, including fluid and gas storage), due to the non-uniqueness of rock-specific physical properties.

Earth-Science Reviews↗
USGS2025

Mapping potential sensitivity to hydrogeomorphic change in the UMRS riverscape

In 2020 the U.S. Geological Survey (USGS), as part of the Upper Mississippi River Restoration (UMRR) Program, began a new project to characterize potential hydrogeomorphic change associated with hydrogeomorphic units (HGUs) and their catenae (units linked by their association with sediment sources and flow origins). The goal of the project was to develop a geographic information system (GIS) database of HGUs for the Upper Mississippi River System (UMRS) available to both scientists and river managers working on UMRR studies and HREP planning and design studies. The characterization was based on a hydrogeomorphic change hierarchical classification developed previously for the UMRS. The products were generated with automated techniques in a GIS using systemic datasets. Landforms were mapped from the 2015 UMRS topobathymetric dataset with geomorphon (shape-based) tools tailored for the large riverscape, valley bottom environments in the UMRS. A clustering analysis was applied to the resulting landforms to identify HGUs associated with zones of perennial low flows, bankfull flows, and overbank floodplains. Catenae were assembled based on the proximity of the units to the main channel, tributary mouths, and side channels from previously published aquatic areas (USACE, 2018) coupled with least-cost flowpath linkages between potential sediment origins and planform change units developed by Rogala, Fitzpatrick, and Henderson (2020). These GIS-based analyses were successful at identifying a range of HGUs using an automated technique with available data across the entire riverscape, with emphasis on those that have the potential for hydrogeomorphic change. Most of the resulting features are depositional, as expected in a large river system. However, this is the first attempt of linking tributary inputs, side channel erosion and levee breaches with their depositional counterparts. The approach was successfully piloted in Pools 8 and 10 in the Upper Impounded Reach and Pool 14 in the Lower Impounded Reach, with next steps for application in reaches of the unimpounded section and Illinois River. This report emphasizes results from Pool 10, which was the focus of most of our attention during the pilot phase.

Illinois, Indiana, Iowa, Minnesota, Missouri, Sout↗
USGS2025

Earth Mapping Resources Initiative protocols—Sampling hard-rock mine waste and perpetual mine water sources

Supporting the overarching goal to evaluate critical minerals nationwide, the mine waste characterization effort in the U.S. Geological Survey (USGS) Earth Mapping Resources Initiative has created a series of protocols to standardize sampling carried out under this effort by the participating State geological surveys and their cooperators. The protocols are based on published, reviewed methods that can be deployed in the field. The protocols include (1) collecting and processing composite samples of mine and mill waste, including tailings, waste rock, gangue, heap leach piles, ore stockpiles, slag, or other mineralized and processed materials and (2) collecting and preserving water samples from perpetual or long-term mine water sources. The protocols also specify information to document on field sheets and detail the collection of geospatial data. The analytical methods used by the USGS and USGS contract laboratories are described in this report, including the data delivery pathway for USGS-derived data.

Scientific Investigations Report↗
USGS2025

Tracing metal sources and groundwater flow paths in the Upper Animas River watershed using rare earth elements and stable isotopes

Groundwater flow paths and processes that govern metal mobility and transport are difficult to characterize in mountainous bedrock watersheds. Despite the difficulty in holistic characterization, conceptual understanding of subsurface hydrologic and geochemical processes is key to developing remediation plans for locations affected by acid mine drainage, such as the Upper Animas River watershed in southwestern Colorado, USA. Stable isotopes of water and rare earth elements were utilized to evaluate groundwater flow and metal sources within this complex catchment. Stable isotope samples collected from draining mine adits and springs display systematic spatial variation wherein sample sites at higher elevations have greater seasonal variability than sites at lower elevations. The Upper Cement Creek watershed, where multiple draining mines are present, displays the lowest seasonal variation in stable isotopic signatures, potentially indicating the presence of a large, well-mixed volume of groundwater storage or interbasin groundwater flow. Rare earth elements display statistically significant variation between different alteration styles in the catchment. Overprinting of regional propylitic alteration is evident based on enrichment of middle rare earth elements in acidic springs and mines that are not spatially associated with surficial exposures of acid generating alteration styles. Europium anomaly and middle rare earth enrichment signatures from two flooded mine tunnels on opposite sides of a watershed divide indicate connections to the same subsurface flooded mine workings.

Colorado↗
USGS2025

Special issue "Effects of surface geology on seismic motion (ESG): General state-of-research"

To understand and predict the behaviors of strong ground motions incurred during devastating earthquakes, studies focused on the “effects of surface geology (ESG) on seismic motions” (following tradition, herein contracted as: “ESG”) have progressed steadily in the last three decades. These ESG studies typically involve in situ measurements of the wavefield and commonly apply both analytical and computational approaches. Concurrently, improvements in ESG related research can be readily attributed to the proliferation of openly accessible strong motion data, as well as advances in computational power. Nevertheless, there remain significant shortfalls in our understanding of the epistemic and aleatory uncertainties associated with the ESG as demonstrated by phenomena from recent deadly earthquake-related site effects. Thus, investigations toward addressing these deficiencies should be underscored in future earthquake related disaster mitigation efforts. This special issue of Earth, Planets, and Space (EPS) is dedicated to the ongoing efforts by the International Association of Seismology and Physics of the Earth’s Interior (IASPEI) and the International Association of Earthquake Engineering (IAEE) that champion activities promoting studies related to the ESG.

Earth, Planets and Space↗
USGS2025

Rare earth element-mineralized carbonatite in the Bear Lodge Alkaline Complex, USA—Ore genesis implications from fluid inclusion characterization

Rare earth element (REE) resources of the Bear Lodge Alkaline Complex, Wyoming, are hosted in variably leached carbonatite spatially related to diatreme breccia pipes. We investigated the genesis of REE and lesser-known gold resources through fluid inclusion analysis of carbonatite, fluorite breccia, and smoky quartz vein samples. Physicochemical characteristics of inclusion-trapped fluids were evaluated using petrography, microthermometry, Raman spectroscopy, decrepitate mound analysis, energy-dispersive spectroscopy, laser ablation inductively coupled plasma mass spectrometry, and noble gas isotope analysis. Microthermometry results reveal three fluid types that affected carbonatite dikes within deeper zones that escaped near-surface, ore-grade REE enrichment: (1) high-temperature (330–432°C) magmatic fluid captured in fine-grained calcite; (2) REE-enriched alkali bicarbonate-sulfate brine; and (3) low-temperature (117–182°C) diluted magmatic or meteoric water. Multiphase brine-like inclusions contain burbankite, nahcolite, strontianite, celestine and alkali sulfate daughter crystals, linking them to early burbankite mineralization. Peripheral smoky quartz and fluorite occurrences at Smith Ridge, 1.5 km from the central carbonatite dike swarm, contain primary inclusions that are Cl-poor and rich in Na-HCO 3 -SO 4 , similar to secondary and pseudosecondary inclusions in carbonatites. Helium isotopes reveal a MORB-like source for carbonatite samples and an older crust signature at Smith Ridge, consistent with the proximal ridge-top exposures of Archean granite. Results from this fluid inclusion study coupled with previous studies of carbonatite mineral paragenesis, show that light REEs (LREEs) were not mobilized great distances. Instead, burbankite crystallized within carbonatite from alkali bicarbonate fluids. With sodium retained in early burbankite, outward-emanating fluids enriched in potassium relative to sodium (higher K:Na) resulted in potassium–ferric iron metasomatism of silicate host rocks. This alkali fractionation was accompanied by fractionation of LREEs and heavy REEs (HREEs), with LREEs dominating the central carbonatite resources. In contrast, areas of peripheral REE mineralization at Bear Lodge are commonly characterized by higher HREE:LREE ratios. The K:Na ratio of associated fenites or alteration assemblages could be indicative of early crystallized burbankite in carbonatites and REE fractionation processes potentially leading to areas of concentrated HREEs with greater supply vulnerabilities.

South Dakota, Wyoming↗
USGS2026

Compact seismicity bursts have different characteristics from regional seismicity

Earthquakes tend to cluster, developing into sequences driven by stress perturbations and transient fault-zone processes. Depending on the driving process, earthquake sequences show differing behaviors. This variability challenges our ability to observe or distinguish these driving processes in high resolution. Here we systematically identify seismicity bursts throughout southern California using new statistical methods and examine their causes with other independent geophysical observations. Seismicity bursts are defined as statistically significant seismicity rate anomalies. Our approach has the advantage of being data driven, depending on few earthquake occurrence assumptions. We find abundant seismicity bursts across southern California, most frequently occurring along the San Jacinto Fault and in the Salton Sea and Coso geothermal fields. These seismicity bursts are highly compact in space and time, often encompassed by a 5 km radius and have durations less than 10 hr. Many of the seismicity bursts have their largest earthquake near the beginning of the sequence, but the largest earthquake is usually not an obvious mainshock. We utilize a variety of independent geophysical data sets to analyze the characteristics of the seismicity bursts, finding that the seismicity bursts have low b-values, low spectral stress drops, and varied stress ratios compared to regional seismicity. These differences suggest that seismicity bursts are driven by transient processes acting frequently across fault networks.

California↗
USGS2026

Recycling radiogenic osmium by crustal foundering in subduction zones: Evidence from pyroxenite xenoliths in the north Andean arc

Removal of mafic-ultramafic lower crust (e.g., via delamination) is fundamental to making andesitic continental crust, yet direct evidence of this process remains elusive. A unique suite of garnet clinopyroxenite and hornblendite (arclogite) xenoliths from the northern volcanic zone of the Andes, erupted to the surface in the mid-Pleistocene, have bulk-rock geochemistry and osmium (Os) isotopic compositions that demonstrate their crustal affinity, yet have equilibration pressures and temperatures below the arc Mohorovičić discontinuity (~53 km) and some as deep as within the sub-arc mantle wedge (~105 km). Garnet websterites from the same xenolith suite, sourced from depths approaching the Wadati-Benioff zone (~140 km), have elemental and isotopic compositions indicative of a mantle origin, likely formed as products of peridotite-liquid reactions above the subducting slab. Variability in bulk-rock 143 Nd/ 144 Nd values and garnet oxygen isotope ratios ( 𝛿 18 O) for these samples is attributed to minor assimilation of subducted components and/or older crustal material, but assimilation alone cannot account for the highly radiogenic Os isotopic compositions. These results provide direct petrologic evidence for modern arclogite formation and foundering in the Andean orogen, the archetypal active continental subduction system. Elevated 187 Os/ 188 Os values of Mercaderes arclogites imply that recycling of gravitationally unstable arc cumulates during continental crust formation introduces highly radiogenic Os into the convective mantle, which has implications for mass fluxes across the crust-mantle boundary as well as the sources and evolution of mantle heterogeneity as seen in orogenic peridotite massifs and in the sources of oceanic basalts

Earth and Planetary Science Letters↗
USGS2025

Probing lower-crustal fault properties with frequency-dependent tidal tremor triggering

The way seismicity responds to periodic stress perturbations offers crucial insights into the processes that can trigger an earthquake. Laboratory and theoretical analyses have shown that the period of imposed forcing and source properties affect the sensitivity to triggering, but frequency-dependent triggering of tectonic faults is poorly understood. The rate of low-frequency earthquakes (LFEs) near Parkfield, California has been found to be strongly correlated with solid earth tides. Tidal forcing acts over multiple frequencies, and the sensitivity to tidal triggering of LFEs therefore provides a unique opportunity to probe the physics of earthquake triggering and underlying fault properties. Here, we find that the response of LFEs to solid earth tides at diurnal and semi-diurnal frequencies is highly variable but spatially coherent along the San Andreas Fault. Using rate-state friction modeling, we find that the variation of the amplitude of tidal modulation is mainly affected by the spatial variation of the background effective normal stress, whereas the spatially varied tidal modulation at diurnal and semi-diurnal frequencies is mainly affected by the fault frictional property. The spatial variations of the LFE response reveal a heterogeneous lower-crustal geologic structure and complex physical faulting processes below the rupture zone of eventual great earthquakes.

California↗
USGS2025

Bayesian mapping of regionally grouped, sparse, univariate earth science data

Some earth science data are naturally grouped by region, and it is often desirable to map these data by region. However, if there are only a few samples within each region, then the map should be smoothed in an appropriate way to mitigate the problems that arise from having only a few samples. A smoothing algorithm based on a Bayesian hierarchical model is developed and presented in this report. This algorithm has several features that make it especially suitable for mapping earth science data: it can account for measurements that are censored, it can process multiple datasets with different measurement errors and different censoring thresholds, and it can calculate the uncertainty in any statistic that is mapped. The algorithm is demonstrated by mapping gold concentrations that are measured in streambed sediments in the Taylor Mountains quadrangle in southwestern Alaska.

Alaska↗
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↗