USGS Science⌕ Search

SEARCH · USGS Science

Search USGS Science

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.

At least 325 records · Page 18Linked to original sources
USGS2026

Groundwater spatial variability within an atoll island: Assessing shallow aquifer heterogeneity with geophysical and physicochemical measurements

This study examines the spatial variability of shallow groundwater on Dhigelaabadhoo Island using electromagnetic induction surveys, groundwater monitoring, and sediment analyses. The research reveals how variations in island morphology—such as differences in elevation, reef flat width, and sediment composition—affect the spatial distribution of groundwater lenses and the overall aquifer dynamics. Saltwater intrusion is especially pronounced in low elevated areas, with narrow reef flat plate and areas where higher hydraulic conductivity—driven by the presence of coarser sediments—is observed, whereas regions characterized by finer sediments, higher elevation, and wider reef flat plates tend to support more symmetric and less saline groundwater lenses. The geophysical investigations reveal that tidal oscillations alter groundwater movement by markedly changing water levels and conductivity, thereby underscoring the critical need to account for temporal dynamics in atoll coastal aquifer systems and the importance of integrating tidal dynamics into the aquifer zone. The findings highlight the significant role of intrinsic morphological and external hydrodynamic factors in shaping groundwater distribution on atoll islands, offering critical insights for sustainable freshwater resource management.

Dhigelaabadhoo Island↗
USGS2026

Computational electromagnetic geophysics for groundwater system studies: A review on established practices and recent advances

Identifying effective solutions for locating groundwater resources and ensuring the quality of drinking water is increasingly urgent, given the challenges posed by climate change and population growth. This review investigates electromagnetic geophysical imaging techniques, in both time- and frequency-domain, that can provide valuable insights for groundwater assessment. We explore computational electromagnetic methods used to evaluate electromagnetic data and several recent hydrogeophysical case studies. As open-source frameworks for modeling electromagnetic geophysical problems become available, a broader range of researchers can interpret their data with computationally advanced software. We provide an overview of documented open-source codes for evaluating electromagnetic data and analyze various hydrological targets in relation to their electromagnetic surveying technique and the computational method applied. Furthermore, we evaluate the potential of advanced computational techniques, including three-dimensional modeling, non-deterministic inversion and machine learning, to couple geophysical with numerical groundwater modeling and apply it in groundwater system studies. Despite obstacles such as complexity and resource demands, our findings indicate that the quantification and integration of predictive uncertainties from both electromagnetic and hydrological data and simulations would significantly improve the reliability of hydrogeophysical models. This can lead to a deeper understanding of groundwater systems and improved management practices.

Journal of Hydrology↗
USGS2026

Seismic characteristics of the April 2024 eruption of Ruang Volcano, North Sulawesi, Indonesia

The two most recent eruptions of Ruang volcano, in 2002 and 2024, have been large (VEI 4) with very short precursory periods, 3 and 6 days, respectively. Background seismicity at Ruang volcano is generally low with 0–2 volcano-tectonic (VT) earthquakes per day. The 2024 eruptions were preceded by a notable increase in VT earthquakes from 11 to 15 April with a rapid escalation from 16 to 17 April in the form of increased VT swarms, and the appearance of low frequency earthquakes and short bursts of VT drumbeats shortly before the eruption. Similar VT seismic unrest was recorded in March 2015, May 2016, and April 2022, although none of these episodes culminated in an eruption. A comparison of seismic energy release between April 2022 and April 2024 shows distinct differences. In 2024, cumulative seismic energy reached its peak within 5 days, just before the onset of eruption on 17 April. In contrast, during the 2022-episode, energy release followed a more gradual increase over 8 days peaking on 16 April, then declining gradually without leading to eruption. Notably, low frequency (LF) and drumbeat earthquakes were absent in non-eruptive unrest episodes. The transition to proximal LFs in 2024 likely signals pressure fluctuations, intensified magma degassing, and conduit formation, acting as an important precursor to eruption. Additionally, the occurrence of short bursts of drumbeat seismicity likely reflects the onset of the shallow magma ascent and was an important short-term precursor to the eruption. Finally, we examine the potential influence of regional tectonic earthquakes on the escalation of precursory unrest in 2024.

North Sulawesi, Ruang volcano↗
USGS2026

The petrogenesis of Þingmúli volcano, East Fjords, Iceland

In this work we revisit Þingmúli volcano (Þ = Th), a classic locality known as an example of a complete tholeiitic differentiation. Þingmúli is a ~ 9.5 Ma extinct central volcano located in the East Fjords of Iceland, in which the whole compositional spectrum from basalt to rhyolites have erupted. These volcanic products have been previously considered as petrogenetically related by an ideal fractionation trend, regardless any temporal relationship or volumetric considerations. Here we report new whole-rock geochemistry, mineral chemistry, isotope analyses, estimation of residence times of the different eruptive deposits, and an update of the original petrogenetic model. Our results highlight that an enriched source, likely spinel lherzolites, generated transitional-alkaline basaltic melts after 15–20% of partial melting at depths of 40–45 km. Many of these basaltic melts erupted at various stages of the volcano's history, while others remained longer in the volcanic plumbing system. These evolved by fractional crystallisation into basaltic andesite magmas with a residence time of ~5 years based on the crystal size distribution of the plagioclase population. Isotopic differences between the basalts/basaltic andesites ( 87 Sr/ 86 Sr ~ 0.7034; 143 Nd/ 144 Nd ~ 0.51315) and the erupted rhyolites ( 87 Sr/ 86 Sr ~ 0.7037; 143 Nd/ 144 Nd ~ 0.51304) indicate that the latter are not petrogenetically related to the former. Therefore, instead of a fractional crystallisation mechanism to generate the rhyolites, we propose the partial melting of ignimbrite layers located beneath the volcano. The broad range of trace element concentrations in andesites and dacites and their different isotopic values compared to the basalts strongly suggest that these magmas have been generated by magma mixing between basaltic and rhyolitic melts, similar to modern day Icelandic volcanoes such as Hekla. These results highlight the need to revisit previously studied Icelandic classic localities and reassess their traditionally proposed petrogenetic models.

Þingmúli volcano, East Fjords↗
USGS2026

Mineral chemistry perspective on remobilization of stored magma at Kamakai'a Hills, Southwest Rift Zone of Kilauea, Island of Hawai'i, USA

Differentiated magmas stored in the rift zones of Kīlauea have received more attention in recent years following eruption of andesite during the early phase of 2018 lower East Rift Zone activity. Despite this growing interest, some of the most voluminous eruptions of differentiated rift zone magmas remain poorly studied. One such eruption, and the most voluminous exposed differentiated flow field at Kīlauea, is the Kamakaiʻa Hills. This eruption took place in the Southwest Rift Zone of Kīlauea, a region that is hypothesized to contain a long-lived rift zone reservoir. The Kamakaiʻa Hills flow field encompasses >250 × 10 6 m 3 of basaltic andesite and basalt compositions with a mineral assemblage of orthopyroxene + clinopyroxene + plagioclase during its early ʻaʻā phase and clinopyroxene + plagioclase + olivine during its late pāhoehoe phase. To better understand storage conditions and magma accumulation, this study focuses on major, minor, and trace elements from the mineral assemblage present within the early ʻaʻā and late pāhoehoe phases. The diversity of clinopyroxene and plagioclase compositions within the early ʻaʻā and late pāhoehoe phases, as well as diverse compositions of plagioclase and orthopyroxene within the early ʻaʻā phase, suggest multiple magma bodies and limited pre-eruption magma mixing within the broader Kamakaiʻa Hills reservoir. Oscillatory zoning patterns (particularly in clinopyroxene) imply processes such as recharge events, magma mixing or mingling, or convection within a differentially cooling, chemically stratified reservoir over protracted time intervals, whereas only limited resorbed mineral textures indicate incomplete mixing of heat and chemically distinct magmas during the dike intrusion that triggered the eruption. Mineral-mineral and mineral-melt thermobarometry indicate predominantly shallow (≤2.5 km depth) crustal storage conditions of the cooled, differentiated magma (∼1100 °C and cooler for the basaltic andesites) to hotter temperatures for the basalts (all >1100 °C). Despite the known large standard errors estimated for mineral-melt and mineral-mineral barometry (10s to >100 MPa), the calculated pressures and depths broadly correspond with earthquake swarm depths beneath the Kamakaiʻa Hills, and drill core and fluid inclusion barometry storage depths of differentiated magmas within the lower East Rift Zone. The Kamakaiʻa Hills differentiated magmas have H 2 O contents (∼0.5 wt%, using plagioclase-melt hygrometry) equivalent to typical Kīlauea basalts. Our data and interpretations demonstrate a complex, long-lived rift zone storage system that consisted of multiple magma bodies and was mobilized into eruption through intrusion of a hotter and more primitive summit-derived (uprift) magma.

Hawaii↗
USGS2026

PFAS remediation in a bioelectrochemical system inoculated with the west branch consortium (WBC-2)

Groundwater contamination by per - and polyfluoroalkyl substances (PFAS) poses a persistent environmental and public health concern. This study evaluates a two-chambered bioelectrochemical system (BES) inoculated with the West Branch Consortium (WBC-2) for PFAS remediation. Under an applied cathodic potential of −450 mV (versus Ag/AgCl), the BES with active WBC-2 achieved >99.0% perfluorooctanesulfonic acid (PFOS) removal within 21 days in deionized water with culture medium and > 98.9% removal of PFOS, perfluorooctanoic acid (PFOA), perfluorohexanoic acid (PFHxA), and perfluorohexanesulfonic acid (PFHxS) in contaminated groundwater after 102 days. Intermediate formation (e.g., PFOA, 6:2 fluorotelomer sulfonate (6:2 FTS), perfluoropropionic acid (PFPrA), perfluorobutanoic acid (PFBA)) and background-corrected fluoride release were consistent with PFOS transformation under anaerobic reducing conditions potentially involving defluorination. Following repeated PFOS spikes (100 μg/L on Days 0, 50, and 399), PFOA, PFPrA, and PFBA accumulated over 664 days. Despite being the dominant accumulated compound, PFOA accounted for <1.8% of the total spiked PFOS mass. Minimal PFOS transformation occurred in controls without active WBC-2, highlighting the importance of microbial metabolism. Biofilm analysis revealed dense colonization of rod-shaped bacteria on carbon fiber brushes. Enrichment of Bacillus , Agrobacterium , and other low-abundance taxa suggests selective adaptation to BES and PFAS conditions. These findings highlight BES driven by electrochemically stimulated microbial activity as a promising strategy for PFAS remediation.

Journal of Water Process Engineering↗
USGS2026

Generating geochemical and mineralogy distributions of soil in the conterminous United States using Bayesian hierarchical spatial models

Characterizing geochemical and mineralogical soil distributions across large spatial extents is essential for understanding mineral resources, ecosystem processes, and environmental risks. Rasters of soil geochemical distributions for the conterminous United States, however, are limited. We present a Bayesian modeling workflow and tool for generating predictive geochemical and mineralogy distribution maps for the conterminous United States using integrated nested Laplace approximation (INLA) with the stochastic partial differential equation approach. By modeling soil geostatistical data with environmental covariates (soil properties, topography, climate, and land cover), we generate predictive distributions of soil geochemistry that can be mapped or extracted for further analyses. As an example, we model the spatial distribution of trace elements in soil relevant to vertebrate health (cobalt, copper, iron, manganese, selenium, and zinc) and provide a workflow that can be used to generate and visualize predictive distributions of 39 other major and trace elements and 21 minerals of the soil survey, supporting a variety of ecological, environmental, and agricultural applications.

MethodsX↗
USGS2026

Spatial variation in effectiveness and costs of sagebrush restoration treatments across the western United States

Natural resource managers strive to improve restoration efficiency across the vast network of public lands of the United States, seeking the best return on investment through innovation, geographic prioritization, and adaptive management. This challenge is exemplified by ecosystems dominated by sagebrush ( Artemisia spp.), which occur across the western United States and have been degraded by increasing wildfire frequency driven by invasive plant spread, among other factors. However, both costs and effectiveness of sagebrush restoration can vary spatially due to management practices in addition to biotic and abiotic factors, and characterizing this variation could inform broad-scale planning. We leveraged published models of sagebrush recovery and treatment costs to predict and evaluate cost-effectiveness of aerial and drill seeding Artemisia spp. across 429 718 km 2 that overlapped 12 US states. Compared to natural recovery, effectiveness for both methods was generally low, yet effectiveness was often greater, and costs were always greater, for drill than aerial seeding, resulting in slightly higher cost-effectiveness on average for drill seeding. Cost-effectiveness for both treatments increased near major roads and in areas with more repeated burns. Cost-effectiveness also increased with predisturbance cover of sagebrush but decreased with herbaceous cover. However, we also identified areas where aerial seeding was more cost-effective than drill seeding, particularly with greater slope. Our results provide spatially explicit estimates and potential mechanisms of cost-effectiveness for two common seeding methods for sagebrush, which can help prioritize limited resources, guide land use, and improve restoration efficiency and effectiveness across public rangelands of the western United States.

western United States↗
USGS2026

Regreening, restoring, and reconnecting a southwestern wetland ecosystem – the Zeedyk wetland

Alluvial wetland ecosystems are vital as biodiversity hotspots but are increasingly threatened by anthropogenic stressors and drought. These pressures are especially acute in arid and semi-arid regions, where eco-hydrologic connectivity is fragile and recovery is slow. This study quantifies the efficacy of nature-based solutions, particularly the ‘Zeedyk approach,’ which employs low-tech Natural Infrastructure in Dryland Streams (NIDS)—including rock detention structures—to slow surface water, raise groundwater tables, and restore wetland function at a spring-fed wetland in Cebolla Canyon, New Mexico, U.S.A. Our results depict a Restoration Feedback Loop that captures stages of change from a healthy wetland in 1935, altered by 20th-century agriculture and grazing, to the re-establishment of the historical flow regime by 2024 documented through an 89-year archive of aerial imagery (1935–2024). By the end of our study period, the Spring-Fed Wetland had expanded by roughly 229% of the original 1935 area, to 4.13 ha. Using 40 years of satellite data, we assess changes in vegetation and hydrology with remote sensing indices. Spatial and temporal analyses reveal significant increases in vegetation greenness and wetness, particularly in an Expanded Wetland subregion, which exhibited ∼3.5x higher wetness and ∼1.5x higher greenness trends compared to adjacent areas. Monthly metrics highlight seasonal variability, with increases in greenness linked to monsoonal rainfall and lateral water redistribution, indicating that restoration impacts extend beyond the primary wetland. This study demonstrates the utility of cloud-based platforms like Google Earth Engine and USGS EarthExplorer for long-term monitoring of wetland restoration, while quantifying the efficacy of the ‘Zeedyk approach’ and demonstrating its potential as a scalable method to restore and conserve wetland meadows in other arid and semi-arid landscapes.

New Mexico↗
USGS2026

Comparative crop yield forecasting using satellite-derived biophysical and agro-climatic predictors in Sub-Saharan Africa

Timely and accurate crop yield forecasting is central to food security early warning systems, particularly in climate-vulnerable regions. While operational forecasting frameworks commonly rely on precipitation and vegetation indices such as NDVI, their ability to provide actionable lead time remains limited. Here, we evaluate the added value of satellite-derived biophysical Essential Climate Variables (ECVs): Leaf Area Index (LAI) and Fraction of Photosynthetically Active Radiation (FAPAR), for forecasting millet yield in Burkina Faso (BF) and maize yield in South Africa (ZA) and Malawi (MW). Using Random Forest models, we quantify forecast skill across the growing season at both national and sub-national scales. Results show that LAI and FAPAR provide effective forecast lead times of approximately 4 months in BF, 2 months in ZA, and up to 6 months in MW relative to harvest. At peak performance, Mean Absolute Percentage Error (MAPE) reaches 19.8% (LAI) and 23.8% (FAPAR) in BF, 12.0% and 9.8% in ZA, and 21.8% and 20.8% in MW, respectively. Across countries, biophysical parameters often outperform NDVI and precipitation, particularly in arid and semi-arid regions. At the sub-national level, LAI and FAPAR enable classification of administrative units into high and moderate-skill forecast units, revealing strong spatial heterogeneity linked to crop dominance. However, forecast skill declines where the target crop is not the dominant type, highlighting an important limitation for operational deployment. Overall, the findings suggest that satellite-derived biophysical parameters can provide earlier and more spatially resolved yield signals than commonly used predictors, with potential to improve the timeliness and effectiveness of food security early warning systems.

Remote Sensing Applications: Society and Environme↗
USGS2026

Hydrologic variability drives environmental and geospatial relationships in Smallmouth Bass (Micropterus dolomieu) distribution

Hydrologic variation is a primary driver of stream ecosystems. Changing hydrology can lead to assemblage shifts and alterations in suitable habitat for freshwater species. As climate change is predicted to alter flow patterns in addition to increasing water temperatures, insight into relationships between species occupancy, hydrology, and temperature is critical for understanding current and future distributions. We examined how hydrologic variability, temperature, and other environmental variables interact to influence Micropterus dolomieu (Smallmouth Bass) occurrence. We used Spatial Stream Network models, allowing for the incorporation of spatial autocorrelation along streams' unique dendritic network, to examine Smallmouth Bass occupancy across a range of hydrologic variation in the Ozark-Ouachita Interior Highlands, USA. Hydrologic variation was the main driver of Smallmouth Bass occurrence, with occurrence more likely in groundwater streams with low hydrologic variation and high flow permanence. For groundwater streams, occurrence was positively associated with summer stream temperature and negatively associated with annual stream temperature. As variation increased, more variables showed significant relationships with occurrence. Distance metrics were important for all models, however as hydrologic disturbance increased, flow connected distance played a lesser role and stream distance played a greater role. Hydrologic variability was the overarching determinant of Smallmouth Bass occurrence and strongly influenced the predictive importance of environmental variables and geospatial relationships. Greater hydrologic variability resulted in stronger statistical relationships between occurrence and environmental variables and an increased importance of system connectivity. As climate change alters hydrologic processes and streams become more variable, understanding and accounting for these shifting relationships is essential.

Arkansas, Kansas, Missouri, Oklahoma↗
USGS2026

Legacy of the fumigant 1,2-dibromo-3-chloropropane (DBCP) in California groundwater

The fumigant pesticide 1,2-dibromo-3-chloropropane (DBCP) was widely used in California agriculture during the 1960s and 1970s before being banned in 1979. Despite this ban, DBCP continues to contaminate groundwater due to its persistence and mobility. This study evaluates the distribution, historical trends, and projected persistence of DBCP in California using data from over 13,000 public supply wells and additional domestic, irrigation, and observation wells (1980-2022). Since 2010, DBCP has been detected in 9% of public supply wells statewide, with higher frequencies in the San Joaquin Valley (21%) and upper Santa Ana River watershed (13%), where DBCP use was most prevalent. Approximately 70% of wells had decreasing concentration trends, whereas increases were more common in deeper wells, indicating downward vertical migration of the DBCP front. Groundwater age estimates show that recharge timing aligns with the 1960s–1970s loading period, enabling reconstruction of peak inputs and providing a basis for age based modeling. To estimate future persistence, we applied a one dimensional advection–dispersion model that simulates long term declines in peak concentrations based on groundwater age, historical loading, and a 38 year degradation half life. Model projections suggest that concentrations above the maximum contaminant level may persist in a declining number of wells until approximately 2080 (range: 2048–2109), with longer persistence in the San Joaquin Valley. The simplified modeling framework, based on age distributions typical of wells capturing peak concentrations, can provide practical regional scale assessment of non-point source contaminants where long-term monitoring exists. This study highlights how the legacy of DBCP contamination will likely affect California's groundwater resources throughout the 21st century.

California↗
USGS2026

GST-1: A high-resolution global sediment thickness model

Global Sediment Thickness 1 (GST-1) is a high-resolution sedimentary thickness model calculated on a 0.125° x 0.125° grid. It modifies the sediment thickness of the 1° x 1° Earth Crustal Model 1 (ECM1) by means of 3D inversions of free air gravity anomalies. GST-1 is calculated by performing structural inversions on high-density contrasts across two crustal boundaries: the sediment – basement interface and the crystalline crust – upper mantle interface. The inversions are calculated in each of ten overlapping 3D models that span the globe. These ten models are merged to obtain the GST-1 global model, providing an eight-fold increase in lateral spatial resolution in comparison with ECM1 and CRUST 1.0. Our sediment thickness model exploits the nearly continuous sampling of gravity data when compared to the irregular, sparse sampling of seismic refraction data. Sediment thickness values in GST-1 are in excellent agreement with independently derived cross sections from well-studied sedimentary basins, and within expected resolution limits of seismic refraction data. GST-1 offers a robust, high resolution global model of sedimentary thickness to support studies of sedimentary basins.

Tectonophysics↗
USGS2026

Comparative assessment of STIC sensors, streamflow and rain gauges for quantifying river connectivity in intermittent systems

In intermittent stream systems, including those occurring in Texas, USA, the severity of low-flow conditions, duration of seasonal disconnection, and frequency of no-flow events have been amplified by drought. Documentation of these no-flow events is necessary to evaluate ecosystem health. However, many intermittent reaches remain un-gauged given that perennial river sec-tions are often prioritized for gauge placement. Our objectives were to 1) document stream flow using Stream Temperature, Intermittency, and Conductivity (STIC) loggers to determine the frequency and duration of no-flow events in intermittent tributaries of the Colorado River, Texas and 2) compare logger data to publicly available data from streamflow discharge and precipitation gauge networks to understand differences among these data types for drying event characterization. We use these comparisons to summarize benefits and limitations of the application of in-stream data loggers. STIC loggers were deployed at 19 sites, one in each pool and riffle habitat of a stream reach. STIC loggers recorded a measurement of relative conductance every six hours from June 2022 to March 2024, which was used to determine the presence or absence of flow connectivity in a reach. No-flow duration among intermittent reaches varied between 37 and 270 days across tributaries during an ongoing drought in the study area. Overall, logger data was more precise than discharge data for characterizing no-flow events or precipitation data when documenting presence of water in the stream channel due to runoff. Lack of discharge gauges in intermittent tributaries left large sections of stream reaches undocumented and resulted in mischaracterization of flow patterns. Drought severity across the tributaries did not follow longitudinal patterns that would be expected by the climatic precipitation gradient of the study area. More research is needed to determine if factors such as population size affect severity. Likewise, precipitation data did not correlate well with logger water presence data, lacking consideration for groundwater recharge, soil hydrophobicity, and surface compaction. This study shows that to monitor no-flow events, detailed spatial datasets are necessary and that STIC loggers are useful tools that provide data to fill spatial information gaps and facilitate more accurate flow characterization and water presence data in intermittent systems.

Texas↗
USGS2026

Upper Triassic igneous rocks of the southern Kenai Peninsula, Alaska – Prelude to Early Jurassic subduction along the western Wrangellia composite terrane margin

New U–Pb zircon geochronology identifies a latest Triassic (ca 214–201 Ma) igneous suite of tuff, hypabyssal dikes, and a pluton on the southern Kenai Peninsula, Alaska. The igneous suite was emplaced within Upper Triassic sedimentary rocks along the southern margin of Western Wrangellia, the western-most fragment of the Wrangellia composite terrane. The igneous rocks range from mafic (50.6% SiO 2 ) to felsic (78.3% SiO 2 ), characteristically have less than 1.55% K 2 O, and generally have low trace element abundances. The tonalitic and trondhjemitic magmas were largely sourced in mafic-rich lower crust and incompletely assimilated quartz and other mineral xenocrysts are common. Fractionation involving plagioclase and amphibole is indicated for some magmas and composite intrusions and igneous xenoliths indicate magma mixing was possible. Paleozoic and Precambrian inherited zircons and initial 87 Sr/ 86 Sr (0.704103–0.705609) and 143 Nd/ 144 Nd (0.512396–0.512777) ratios indicate that the Western Wrangellia crustal sources are heterogeneous and contain sialic components. The latest Triassic magmatism reflects processes that preceded Early Jurassic subduction along the Wrangellia composite terrane and Pacific Ocean plate boundary. These processes involved heating and melting of mantle lithosphere and lower crust as mantle instabilities accompanied the breaking of the plate boundary linkages. The Late Triassic transition to subduction along the Wrangellia composite terrane margin coincided with the transition to subduction cessation in the Late Triassic arcs of the western Intermontane terranes of Canada. The shift to subduction along the outboard Wrangellia composite terrane margin marks the beginning of the Pacific Ocean–Cordillera plate interactions that came to dominate the tectonic evolution of the northern Cordillera from the Early Jurassic to today.

Alaska↗
USGS2026

Geophysical architecture and geochronology of the Neoarchean Mentor anorthosite intrusive complex, northwestern Minnesota: Largest anorthosite complex of the Superior Province?

The Neoarchean Mentor anorthosite intrusive complex (MAIC) lies within the Wawa subprovince in northwestern Minnesota, in a region where the Wawa, Quetico, and Wabigoon subprovinces are juxtaposed in close proximity. Archean rocks are entirely concealed, and interpretations are developed from aeromagnetic, gravity, and borehole samples. The MAIC includes both anorthosite (dense, weakly magnetized) and oxide-rich gabbro (strongly magnetized) lithologies. Anorthosite is coarse-grained to megacrystic, intensely altered, and locally brecciated. Pervasive epidote alteration enhanced the density of the anorthosite via introduction of higher-density mineral assemblages, explaining why the MAIC produces a significant gravity high. Oxide-rich gabbro forms a border phase of the MAIC and has potential for vanadium, chromium, and titanium mineralization, and produces a strong aeromagnetic high. The MAIC is interpreted to extend over an area of 640 km 2 , making it the largest known anorthosite complex of the Superior Province, as measured by preserved areal extent. Modeling indicates the MAIC extends more than 6 km into the subsurface. A new Pb–Pb zircon age of 2737.2 ± 4.5 Ma is interpretated as the crystallization age of anorthosite within the MAIC, showing that the MAIC formed well before the ca. 2690 Ma Shebandowanian orogeny, and raising new questions about correlations with other parts of the Wawa subprovince. A low-density batholith, here informally called the Fertile batholith, is interpreted to intrude the southern part of the MAIC. A new Pb–Pb zircon age of 2701.1 ± 6 Ma is interpreted as the magmatic age of the Fertile batholith.

Minnesota↗
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↗
USGS2026

New measurements indicate that natural geologic methane emissions from microseepage in the Michigan Basin are likely negligible

The magnitude of natural geologic methane (CH 4 ) emissions to the atmosphere (including emissions of fossil CH 4 from offshore and onshore gas and oil seeps, diffuse microseepage, mud volcanoes, volcanic vents, and geothermal areas) is highly uncertain. The largest component of geologic emissions is thought to be microseepage, which is the diffuse flux of CH 4 from soils across large areas of productive hydrocarbon basins. The accuracy of existing bottom-up estimates of microseepage is limited by low spatial coverage of published microseepage measurements. We present the first soil—atmosphere CH 4 flux measurements from Michigan Basin, USA. Results from 335 measurements taken during summer and winter seasons across a large portion of the basin suggest microseepage is nonexistent in the sampled region. Even areas with predictive features for microseepage (e.g., underlying mature, organically rich source rocks, proven gas accumulations, faults, and lineaments) yield null or negative fluxes, suggesting that CH 4 emissions from microseepage are negligible throughout our study region. A Monte Carlo method was used to place an upper limit on the regional-mean microseepage flux, in which synthetic patchy microseepage distributions were generated and tested against our observations to assess the impact of possible emission hot spots that were missed by our sampling strategy. Our analysis finds it is very unlikely that regional-mean emissions are as high as assumed in a previous global microseepage study. The observed lack of seepage may be explained by groundwater flow, active methanotrophy, glacial sediments, and bedded salt deposits, which could inhibit vertical gas migration and release to atmosphere.

Michigan↗