USGS Science⌕ Search

SEARCH · USGS Science

Results for “Journal of Coastal Research”

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,660 records · Page 72Linked to original sources
USGS2026

Global pegmatite-hosted lithium, cesium, and rubidium resources: A dataset for grade and tonnage modeling

Quantitative mineral resource assessments of potential undiscovered deposits can inform future mineral supply scenarios, but their accuracy is conditional on building robust grade and tonnage models of known deposits. This study presents an up-to-date global compilation and analysis of recently discovered and original, in-situ pegmatite-hosted Li, Cs, and Rb resources prior to historic production. Our analysis yields a median tonnage of 21.2 million tons (Mt) and grade of 1.12% Li 2 O, respectively, for global Li pegmatite deposits ( n = 73). The grades and tonnages of Li pegmatite resources vary depending on the age of the bedrock host domain, pegmatite crystallization age, and primary ore mineralogy. Lithium pegmatite resources hosted in Archean to transitional Archean-Paleoproterozoic domains have the largest median tonnage (29.8 Mt; n = 38), and those hosted in Paleoproterozoic to Mesoproterozoic domains have smaller median tonnages (6.5 Mt; n = 16). Cesium deposits where pollucite is the primary ore mineral have a bimodal grade distribution, with modes of 2.40 and 0.035 wt% Cs 2 O for high- and low-grade deposits, respectively, while Rb deposits are more unimodal with a median grade of 0.247 wt% Rb 2 O. Pegmatite-hosted Cs and Rb resources have median tonnages of 7.6 and 6.3 Mt, respectively. Covariation between ore mineralogy and the degree of crustal enrichment in pegmatite-hosted deposits is diagnostic of petrogenetic differences, including melt source characteristics, magma evolution, or variable degrees of volatile solubility. The Li pegmatite compilation is suitable for fitting robust numerical models to support quantitative assessments. More well-defined Rb and Cs pegmatite resources are required for quantitative assessments, but these data provide useful information about original in-place resources for framing supply discussions.

Ore Geology Reviews↗
USGS1995

Streamflow, ground-water recharge and discharge, and characteristics of surficial deposits in Buzzards Bay Basin, southeastern Massachusetts

Streamflow measurements at 14 low-flow partial-record stations and two discontinued streamflow-gaging stations were related to concurrent streamflows at six long-term streamflow-gaging stations to estimate streamflows at selected flow durations from 50 to 99 percent for water years 1967-91. At low flows, stream discharge per square mile generally increased with increasing percentage of subbasin underlain by stratified-drift deposits. At the 70-, 90-, and 99-percent flow durations, subbasins underlain primarily by stratified-drift deposits had stream discharges per square mile at least 2, 4, and 8 times greater, respectively, than subbasins underlain primarily by till and bedrock deposits. Streamflow measured three times at the Paskamanset River showed that streamflows downstream of municipal pumped wells increased at a lower rate per unit drainage area than streamflows upstream of the pumped wells. When pumpage from wells near the river was added to measured streamflow at the most downstream station, the sum was similar to the potential streamflow at the station if its streamflow per unit area was the same as that upstream from the pumped wells. Mean ground-water recharge and discharge rates were computed from continuous records of daily mean discharge during water years 1967-91 for six streamflow-gaging stations in southeastern Massachusetts and Rhode Island. Estimates of mean ground-water recharge were 19.7 to 22.6 inches per year for stations with drainage areas primarily underlain by till and bedrock deposits, and 23.8 to 25.2 inches per year for stations with drainage areas primarily underlain by stratified-drift deposits. During drought years, annual ground-water recharge can be less than 50 percent of the mean ground-water recharge rate for water years 1967-91. Estimates of ground-water discharge rates generally were 1 to 4 inches per year less than estimates of ground-water recharge rates. Ground-water discharge-duration curves were calculated from hydrograph separation of streamflow records at the six continuous streamflow-gaging stations for water years 1967-91. When streamflows were less than the median streamflow (50-percent duration), ground-water discharges averaged 85.7 percent of the total streamflow at the six streamflow-gaging stations. Ground-water discharges were estimated at the 14 low-flow partial-record stations and two discontinued streamflow-gaging stations by multiplying streamflows at the selected flow durations from 50 to 99 percent for water years 1967-91 by 85.7 percent. Data on surficial deposits in the basin were obtained by drilling 10 wells, by conducting 12 seismic-refraction surveys, and by collecting well log data for 313 wells from consulting firms and the Massachusetts Department of Environmental Management. These data indicated that a few minor corrections could be made to existing saturated thickness maps for stratified-drift deposits in Buzzards Bay Basin. Few well logs were detailed enough to estimate transmissivities; thus existing transmissivity maps were not updated for stratified-drift deposits in the basin.

Massachusetts↗
USGS1977

Leachate plumes in ground water from Babylon and Islip landfills, Long Island, New York

Landfills operated by the Towns of Babylon and Islip in southwest and central Suffolk County contain urban refuse, incinerated garbage, and scavenger (cesspool) waste; some industrial refuse is deposited at the Babylon site. The Islip landfill was started in 1933, the Babylon landfill in 1947. The landfills are in contact with and discharge leachate into the highly permeable upper glacial aquifer (hydraulic conducitvity 190 to 500 ft/d [feet per day]). The aquifer is 74 feet thick at the Babylon landfill and 170 feet thick at the Islip landfill. The leachate-enriched water occupies the entire thickness of the aquifer beneath both landfills, but hydrologic boundaries retard downward migration of the plumes to deeper aquifers. The Babylon plume is 1,900 feet wide at the landfill and narrows to about 700 feet near its terminus 10,000 feet from the landfill. The Islip plume is 1,400 feet wide at the landfill and narrows to 500 feet near its terminus 5,000 feet from the landfill. Hydrochemical maps and sections show the distribution of the major chemical constituents of the plumes. The most highly leachate-enriched ground water obtained was from the Babylon site; it contained 860 mg/L (milligrams per liter) sodium, 110 mg/L potassium, 565 mg/L calcium, 100 mg/L magnesium, 2,700 mg/L bicarbonate, and 1,300 mg/L chloride. Sulfate was notably absent or in low concentration in most parts of both plumes. Nitrogen in plume water was mostly in the form NH4+, and concentrations as high as 90 mg/L were found; concentrations of nitrogen as N in the plume were less than 10 mg/L. As much as 440 mg/L iron and 190 mg/L manganese were found in the leachate-enriched water. Samples were also tested for arsenic, boron, cadmium, cobalt, chromium, copper, mercury, nickel, lead, selenium, strontium, and zinc. Boron was more or less ubiquitous and was found in concentrations as high as 2 mg/L. Organic carbon was found in concentrations as high as 2,250 mg/L in the most highly leachate-enriched water but attenuated rapidly to less than 20 mg/L. Dissolved-solids concentrations near the landfills were between 400 and 3,000 mg/L at Babylon and between 500 and 1,500 mg/L at Islip. Ground-water temperatures near the landfills exceed those in ambient water by as much as 7°C at Babylon and 16°C at Islip. Heat contributed by the landfills was mostly dissipated within 0.4 mi of the landfill, butat Islip, the warm leachate-enriched water extended 0.5 mi downgradient. The entrance of leachate into the less dense ground water as pulsations after rainfall may explain the presence of high leachate enrichment near the bottom of the aquifer. A comparison of the physical characteristics of leachate-enriched ground water with those of ambient water suggests that the downward movement of leachate results from its greater density. Simulation of the movement and dispersion of the Babylon plume with a mathematical dispersion model indicated the coefficient of longitudinal dispersion to be about 60 ft 2 /d (feet squared per day) and the ground-water velocity to be 1 ft/d. However, the velocity determined from the hydraulic gradient and public-supply wells in the area was 4 ft/d, which would cause a plume four times as long as that predicted by the mathmatical dispersion model. At the Islip site, the plume was one-third the length calculated on the basis of the age of the landfill. The shortness of the plumes has not been explained; it may be a result of the leachate's having been too dilute to form a plume during the early years of the landfills.

New York↗
USGS2025

Igneous differentiation in arcs

This chapter assesses crystallization-differentiation of convergent-margin magmas both from the standpoints of natural rock suites and from laboratory phase-equilibrium experiments. Crystallization-differentiation of dry subalkaline and alkaline basalts are summarized at the beginning for reference and context. Experiments crystallizing arc basalts through andesites at middle to upper crustal pressures (100–400 MPa), saturated either with H 2 O or mixed H 2 O-CO 2 , produce liquid lines of descent similar to arc igneous suites, with higher f O 2 s giving results closer to calc-alkaline magma series. Early enrichment of melts in SiO 2 with little or no enrichment in FeO* (total Fe as FeO) comes about by the combined effects of H 2 O suppressing crystallization of plagioclase relative to olivine, shifting plagioclase to more calcic – therefore lower SiO 2 – compositions, and suppressing silicate minerals to greater extents than Fe-rich spinel. Also contributing to SiO 2 enrichment is crystallization near the liquidus of higher proportions of olivine relative to calcic clinopyroxene, compared with MORB, but it is unclear if this is brought about by H 2 O or elevated pressure. Exhumed arc terranes, however, generally fail to expose evidence for primitive liquids having reached and differentiated at middle to upper crustal depths, at least in any abundance. Instead, most differentiation in arcs appears to take place in the lower crust or near the crust-mantle transition, and although hydrous experiments at such pressures produce liquids broadly similar to arc suites in their SiO 2 enrichment, the experimental liquids differ in becoming peraluminous at early stages of differentiation while melts have basaltic, basaltic andesitic, or andesitic SiO 2 concentrations. Crystallization of abundant calcic clinopyroxene close to the liquidus causes the high-pressure experimental liquids to become peraluminous with only small increases in SiO 2 concentration. In contrast, natural arc volcanic and plutonic suites become gradually less metaluminous with increasing SiO 2 , crossing to peraluminous at dacitic to rhyolitic SiO 2 concentrations. Hybridism, encompassing magma-mixing, assimilation of old crust, and incorporation of near-solidus melts from only shortly older deep, antecedent arc intrusions may accompany natural crystallization-differentiation, preventing melts from attaining peraluminous compositions at low-SiO 2 concentrations. If so, then nearly all arc suites are to substantial extents composites, rather than preserving simple, progressive liquid lines of descent. Some comments and assessments are also made on the participation of garnet in arc differentiation and on the production of arc batholiths.

Book chapter↗
USGS2025

Critical minerals in mine waste

Introduction Critical minerals are commodities with vulnerable supply chains that play a vital role in supporting the United States’ economy, national defense and security, emerging technologies, and energy independence. The prosperity of our Nation depends on generating a resilient supply of domestic critical minerals; mine waste may be an untapped source of these commodities. Mine waste from centuries of legacy mining persist on the landscape and may contain critical minerals and other valuable commodities previously deemed uneconomic to recover. At modern mines, the financial viability of recovering byproduct critical minerals, which are not the primary target, may be marginal and can ultimately destine them to mine waste. Further, mine waste can be a liability for the mining company or, at legacy mines, the taxpayer because of its effect on the landscape. The U.S. Geological Survey (USGS) has several initiatives to evaluate critical mineral resources in various types of waste. This factsheet highlights studies of mine waste carried out by USGS scientists at the Geology, Energy & Minerals Science Center in collaboration with other science centers funded through the USGS Mineral Resources Program. Recovery of critical minerals from mine waste can aid in remediation efforts and increase domestic supply of vital mineral resources.

Fact Sheet↗
USGS2026

Ground-motion characterization for the 2025 U.S. National Seismic Hazard Model for Puerto Rico and the U.S. Virgin Islands

We develop the ground-motion characterization (GMC) for the 2025 U.S. National Seismic Hazard Model for Puerto Rico and the U.S. Virgin Islands (NSHM-PRVI) for earthquakes in active crustal, subduction interface, and subduction intraslab regimes. Using ground-motion models (GMMs) from the Next-Generation Attenuation (NGA)-West2 and NGA-Subduction projects, the GMC is parameterized by scaled-backbone models for median ground motions and by independent logic trees of aleatory variability. We introduce several novel GMC features into the U.S. NSHM: (1) use of regional ground-motion data for modeling median ground motions; (2) development of scaled-backbone models for median ground motions; and (3) development of independent logic trees of aleatory variability from variance components of GMMs and computed from multiple ground-motion datasets, and incorporating regional ground-motion data effects on variability. We compute probabilistic seismic hazard curves and maps to evaluate the contributions from the GMC components and for comparison with the 2003 GMC. Contributions of the GMC to epistemic uncertainty in seismic hazard are evaluated through spatial variations in epistemic uncertainty in hazard maps, comparisons between mean hazard curves and fractiles, as well as investigations of the contributions of logic-tree branches to hazard maps and curves. Comparisons between seismic hazard from the 2025 and 2003 GMCs allow examination of the changes in hazard curves and mapped values with 2% and 10% probabilities of exceedance in 50 years. The 2025 GMC exhibits modest changes in median predictions, relative to the 2003 GMC; overall, values of aleatory variability are higher, except at long periods ( ), resulting in changes to probabilistic ground motions at low probabilities of exceedance (<10% probability of exceedance in 50 years). Changes in hazard at 2% and 10% probabilities of exceedance in 50 years are also relatively modest (within 20%) at most sites, with the impacts from the 2025 GMC exhibiting minor reductions and increases at 0.2 and 1.0 s periods, respectively, relative to the 2003 GMC.

Puerto Rico, U.S. Virgin Islands↗
USGS2026

High-resolution transboundary vegetation community maps of the Sonoran and Mojave Desert ecoregion to support critical landscape conservation planning and habitat management needs

We produced a 30-m resolution binational land cover map of Bird Conservation Region 33 (BCR 33) for the U.S. North American Bird Conservation Initiative. The region covers large portions of the Sonoran and Mojave Deserts. The map can support the U.S. Fish and Wildlife Service (FWS) Migratory Bird Program’s recovery planning efforts and constitutes the first known binational land cover dataset spanning sections of the United States–Mexico border and using a consistent classification system for both countries. The mapped region includes 152 distinct land cover classes, covering a total area of 38,421,453 ha (148,345 mi 2 ), of which 13,148,345 ha (52,706 mi 2 ) are located in Mexico and 24,770,640 ha (95,639 mi 2 ) in the United States. We primarily used Landsat 8 (OLI) imagery, supplemented by limited ground surveys from two field campaigns, drone-based aerial data, and existing vegetation classification frameworks from both countries. The classification applied a data-fusion approach integrating 30-m Landsat 8 imagery, decadal phenology metrics from vegetation indices, and a random forest model trained mainly with datasets from a comprehensive national mapping project from the U.S. Geological Survey (USGS) GAP Analysis Project (GAP) and federal wildland fire agencies’ Landscape Fire and Resource Management Planning Tools (LANDFIRE) (GAP/LANDFIRE) [United States side] and the National Institute of Statistics and Geography (INEGI) [Mexico side] as well as land cover maps and opportunistic open-access and field observations. Mapping of the full BCR 33 region was carried out in two phases: 1) Phase I, the prototype map, covered a smaller portion of the transboundary area and identified 31 land cover classes, and 2) Phase II, the full BCR 33 map (refer to Figure 1), which resulted in 152 land cover classes. Using a Random Forest classifier, we achieved an overall prediction accuracy of 92% for the Phase I map and 87% for the Phase II full region map. This slight decrease can be attributed to working on a larger, more complex area with a greater number of land cover classes. No formal validation was conducted, aside from using a subset of the collected field observations and training data to assess model performance during and after training. The training sites were further verified using Google Earth (Google, 2026) imagery. Two undergraduate students who worked for over a year visually inspected imagery and open access public images to confirm each training site during model training using in-house developed, online, visual tools. A portion of this field training data was reserved for model validation, and the corresponding results are to be presented in later sections. The project developed an end-to-end, medium- and fine-resolution remote sensing–based data fusion mapping approach. This effort produced a map (Nagler et al., 2025) and the online tools to support a dynamic, live, online map for visualizing the transboundary vegetation communities in BCR 33. The toolset is currently hosted by the University of Arizona (UofA) Vegetation Index and Phenology (VIP) Lab to support FWS partners (https://vip.arizona.edu/viplab_data_explorer?LCM_BCR33). The online map is designed to allow rapid updates using new training, validation, or correction data, making it dynamic and maintainable. The approach we took established a framework for rapid updating and correction of land cover maps, as the model can be quickly retrained with new field observations, updated training data, or other sources. This enables dynamic mapping and change detection of the region’s vegetation. This framework is an advance in data fusion and crowdsourced mapping of complex, vulnerable regions, providing support to regional stakeholders and the wider user community. This transboundary map can inform the protection, conservation, and restoration of vegetation, habitat, and ecosystems, particularly for threatened and endangered species across the two nations using consistent and harmonized binational mapping systems. Beyond supporting land management decisions and stakeholders in the transboundary desert ecoregions, this BCR 33 mapping effort establishes a foundation for future rapid, low-cost, cross-border land cover mapping that can benefit and advance ecosystem management.

Arizona, Baja California, California, Nevada, Sina↗
USGS2025

Water use across the conterminous United States, water years 2010–20

Withdrawals of water for human use are fundamental to the evaluation of the Nation’s water availability. This chapter provides an analysis of public supply, crop irrigation, and thermoelectric power water use for the conterminous United States (CONUS) during water years 2010–20. These three categories account for about 90 percent of water withdrawals in the Nation. The values presented here are based on modeling approaches that estimate water use at temporal (monthly) and spatial scales (12-digit hydrologic unit code—small watersheds sized 50–100 square kilometers) compatible for integration into a broader national assessment of water availability. Models also provide an understanding of factors that influence water use. An estimated 244,817 million gallons per day (Mgal/d; 28,677 million cubic meters per month [Mm 3 /mo]) were withdrawn on average within the CONUS during water years 2010–20 from fresh water and saline water for crop irrigation, public supply, and thermoelectric power, with shares of 43, 14.5, and 42.5 percent for each of these categories, respectively. In the same period, estimated withdrawals and consumptive use (1) for public supply were 35,400 and 4,219 Mgal/d (4,081 and 486 Mm 3 /mo), respectively; (2) for crop irrigation were 105,497 and 75,698 Mgal/d (12,147 and 8,716 Mm 3 /mo), respectively; and (3) for thermoelectric power from fresh water were 82,656 and 2,904 Mgal/d (9,952 and 345 Mm 3 /mo), respectively. Withdrawals for these categories of water use are highly spatially variable, with western States dominated by crop irrigation and eastern States dominated by thermoelectric-power water use. Public supply accounts for the largest percentage of water use in several heavily populated northeastern States. Reliance on groundwater compared to surface water depends on the availability of water sources and the type of water use. For public supply, withdrawals from groundwater are greater than withdrawals from surface water in the Western aggregated hydrologic regions, whereas the balance shifts to more surface water for the rest of the CONUS. In all aggregated hydrologic regions, the predominant source of water for crop irrigation is groundwater. Most thermoelectric power facilities in the eastern half of the CONUS use surface water from freshwater and saline sources; most thermoelectric power facilities in the western half of the CONUS use groundwater.

conterminous United States↗
USGS2025

How do hydrological variability and human activities control the spatiotemporal changes of riverine nitrogen export in the Upper Mississippi River Basin?

Excessive nitrogen export from agricultural watersheds remains a critical water quality challenge, with the Upper Mississippi River Basin (UMRB) significantly contributing to downstream eutrophication and hypoxia in the Gulf. This study investigates the spatiotemporal dynamics of riverine nitrate plus nitrite (NO 3 – + NO 2 – -N) export across the UMRB at high spatial resolution (12-digit Hydrologic Unit Codes or HUC12 subwatershed scale) during 2001–2020 and quantifies the effects of anthropogenic activities and hydrological variability on riverine NO 3 – + NO 2 – -N export changes in the region between 2001–2005 and 2016–2020. Our results revealed hotspots of substantial increases in NO 3 – + NO 2 – -N yields across the UMRB, with distinct regional patterns in driving factors. Over the entire UMRB, NO 3 – + NO 2 – -N yields increased by 9.7 kg/ha/yr on average from 2001–2005 to 2016–2020, with anthropogenic activities contributing 4.8 kg/ha/yr and hydrological variability contributing 4.9 kg/ha/yr. The northern and western UMRB had combined influences from both anthropogenic activities and hydrological variability, while the east-central regions had predominantly hydrologically driven changes. Agricultural sources, including fertilizer, manure, and biological nitrogen fixation, collectively contributed over 80% of NO 3 – + NO 2 – -N loading throughout the basin. This framework for disentangling human and hydrological impacts provides critical insights for developing effective and targeted watershed management strategies to reduce nutrient losses and improve water quality.

Upper Mississippi River Basin↗
USGS2025

A probabilistic assessment methodology for the evaluation of geologic energy storage capacity—Natural gas storage in depleted hydrocarbon reservoirs

The need for energy storage, particularly underground, where capacity and duration may far exceed battery storage technologies, is especially relevant given the increasing demands for reliable power alongside the development of intermittent renewable electricity sources. Geologic energy storage facilities already exist, and expanded use would enable storing gases such as methane and hydrogen. In 2018, a National Academies of Sciences, Engineering, and Medicine report, “Future Directions for the U.S. Geological Survey's Energy Resources Program,” recommended that the U.S. Geological Survey (USGS) prioritize assessing underground energy storage in geologic formations in the United States. The U.S. Geological Survey has since developed a methodology for assessing natural gas storage capacities in depleted hydrocarbon reservoirs on a national scale. The methodology introduced in this report prescribes three approaches for calculating gas storage capacity. This methodology relies on the availability of input data, including cumulative hydrocarbon production records, reservoir petrophysical properties, and reservoir pressure data. Assessment inputs can be obtained from public, State-level databases and propriety national-scale databases, although the use of analogs could be warranted for estimating input parameters. Probabilistic assessment results are aggregated to play, petroleum province, regional, and national scales. The steps defined in this report are demonstrated on the Michigan Basin Province, which includes the Mississippian Sandstone Gas Play and the Clinton Structural Play. This methodology could be used to systematically and consistently assess hydrocarbon plays and provinces for natural gas storage capacity across the United States.

Michigan↗
USGS2026

New constraints on northeast Seattle basin structure from converted seismic waves

The Seattle basin is a deep sedimentary basin in the Seattle–Bellevue, Washington metropolitan area within the Puget Lowland of Washington State. We determine the structure of a portion of the basin and the underlying basement using analysis of P waves converted from direct S incident from below. A deep local crustal event beneath Monroe, about 35 km northeast of Seattle, was recorded by a 100‐station nodal array deployed in 2019. The event produced a variety of coherent seismic phases, including converted waves from the sediment—basement boundary, internal structure within the basin, and additional crustal discontinuities. Using observed Sp converted waves, we apply an adjoint‐based full waveform inversion (FWI) method to determine the amplitude and extent of seismic discontinuities at depth. We find the strongest source of converted waves for this event lies ∼6 to 7 km depth below northern Lake Washington, interpreted to be the local depth to basement rock. The newly imaged shallow basement structure may be part of a deformation zone associated with the Siletzia eastern boundary. Our results highlight the utility of converted seismic waves recorded by a dense array, combined with an FWI method, to illuminate crustal structure.

Washington↗
USGS2026

Widespread abyssal turbidites record megathrust earthquake-triggered landslides and coseismic deformation in the Cascadia subduction zone

Abyssal marine turbidites provide some of the longest and most spatially extensive records of subduction zone earthquake recurrence globally; however, correlation of these deposits over long distances and interpretation of synchronous emplacement requires both an understanding of the turbidite generating systems and precise dating. Here, we present an integrated suite of high-resolution bathymetry, subbottom profiles, and sediment cores from combined autonomous underwater vehicle, remotely operated vehicle, and ship-based studies at a key paleoseismic site in the southern Cascadia subduction zone. We demonstrate how widespread, earthquake-triggered landslides on the lower slope deposit discrete, proximal mass transport deposits (MTDs) that grade offshore into complex, interfingered abyssal turbidites, which correspond to records of megathrust earthquake history. We propose accretion and oversteepening of thrust folds on the lower slope both preconditions the slope to fail and provides a perpetual source of unstable material to fail during every earthquake cycle. Furthermore, we suggest the periodic and pervasive landsliding indicates coseismic deformation of the outer accretionary wedge during megathrust rupture.

California, Oregon↗
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

Simulation of groundwater flow in Wake County, North Carolina, 2000 through 2070

In 2019, the U.S. Geological Survey and Wake County Environmental Services began a collaborative study to evaluate groundwater resources and long-term groundwater availability in the county’s fractured-rock groundwater system. Wake County, in central North Carolina, is experiencing rapid population growth, associated land development, and changing water use. Hydrogeologic data including groundwater levels, aquifer testing, borehole fracture flow measurements, water-quality samples, and groundwater age-dating tracers were collected, along with findings from previous investigations, to help inform a conceptual model of the flow system used to develop a modular three-dimensional finite-difference groundwater-flow model (MODFLOW) for simulating historical and future groundwater conditions from 2000 to 2070. Hydraulic conductivity and transmissivity ranges were estimated from 17 slug tests and 21 borehole-flow measurements. Groundwater-quality analytical results from 19 sampling sites indicate that oxidation-reduction (redox) conditions varied within the regolith and bedrock and that minimal evaporation occurred before recharge entered the groundwater system. Age dating revealed mixtures of older and younger water, ranging from the 1940s to the 1990s—indicating variable flow pathways of recharge within permeable bedrock fracture zones. To simplify the complex fractured-rock groundwater system, two layers representing the regolith and the fractured bedrock were used in the MODFLOW model. Model calibration included parameter estimation and provided a reasonable fit to observed groundwater levels and estimated stream base flows. The model forecast scenarios incorporated future climate-model data for two emissions scenarios with land cover change projections to simulate potential impacts to future groundwater levels, recharge, and base flows. Recharge and base flow projections were largely within historical ranges, with no apparent long-term trends, but did indicate a slight downward shift in median values—likely, in part, because of differences in spatial resolution of input climate datasets. Seasonal patterns were consistent with historical data, with projections of possible increases in future winter recharge. Model limitations are discussed, and additional monitoring and model refinement needs are highlighted to support decision making for local groundwater management.

North Carolina↗
USGS2025

Evidence for nonlocal sediment transport on hillslopes from fault scarp morphology

Hillslope sediment transport processes such as bioturbation, rainsplash, and granular mechanics occur across the entire planet. Yet, it remains uncertain how these small-scale processes act together to shape landscapes. Longstanding hillslope diffusion theory posits that hillslope processes are spatially limited, whereas new concepts of nonlocal sediment transport argue otherwise. However, each theory produces subtly different, but distinct, predictions for the evolution of fault scarps. We use the topographic change of fault scarps to demonstrate that hillslope processes produce nonlocal sediment transport. Analysis of a global compilation of 340 dated single-earthquake scarp profiles reveals a statistically significant ( p < 0.05 ) relationship between scarp age and scarp asymmetry, here defined as the ratio of imaginary to real components of the Fourier transform of absolute slope. Numerical simulations show that nonlocal models predict this relationship, whereas hillslope diffusion models do not. To further investigate this result, we examined the depositional geometry of a well-exposed colluvial wedge along the Wasatch fault in central Utah, United States. Our quantitative comparison between the exposure and numerical simulations reveals better agreement with the nonlocal model. Nonlocal sediment transport theory appears to best capture the physics of how hillslope processes shape fault scarps, yet hillslope diffusion provides a useful approximation in many cases. As the processes that act on fault scarps are nearly identical to those acting on hillslopes, our results provide evidence supporting nonlocality as a generalized model of hillslope sediment transport.

Utah↗
USGS2024

3-D geological modeling for numerical flow simulation studies of gas hydrate reservoirs at the Kuparuk State 7-11-12 Pad in the Prudhoe Bay Unit on the Alaska North Slope

Accurate reservoir evaluation requires reliable three-dimensional (3-D) geological models. This study conducted 3-D geological modeling for numerical flow simulation of the B1 sand gas hydrate reservoir at the Kuparuk State 7-11-12 pad, Prudhoe Bay Unit, Alaska North Slope. The model integrates well logs, core, and seismic data to address spatial heterogeneity in geological structures and reservoir properties. Two modeling types were performed: structural framework modeling and petrophysical property modeling. For structural framework modeling, seismic data and well log markers were used to reproduce subsurface structures characterized by a normal fault system. A volume-based modeling algorithm and stair-step gridding were applied. The resulting 3-D model comprised 2,640,000 grid cells across 264 layers, including seven fault grids. For petrophysical property modeling, total porosity was initially modeled using sequential Gaussian simulation with collocated cokriging. To reproduce the upward coarsening of the B1 sand, upscaled log-derived total porosity and a 3-D trend depicting total porosity variation were used as primary and secondary data, respectively. Gas hydrate saturation distribution was modeled similarly, with secondary data from estimated porosity distribution and seismic-derived acoustic impedance map enhancing accuracy. Results indicate higher gas hydrate saturation in the upper part of the B1 sand and areas with higher acoustic impedance. Intrinsic permeability was modeled from the total porosity and clay-bound water volume, and effective permeability was derived from the gas hydrate saturation and intrinsic permeability distributions based on the “Tokyo model”. Effective permeability distributions were influenced by the total porosity, gas hydrate saturation, and intrinsic permeability. Within the same layer, higher gas hydrate saturation leads to decreased effective permeability. In total, 100 sets of multiple scenarios were prepared, providing input data for dynamic flow simulations to evaluate the effects of lateral heterogeneity in reservoir properties and the hydraulic characteristics of faults on production behavior for preassessment before the long-term production test.

Alaska↗
USGS2024

Methods for quantifying interactions between groundwater and surface water

Driven by the need for integrated management of groundwater (GW) and surface water (SW), quantification of GW–SW interactions and associated contaminant transport has become increasingly important. This is due to their substantial impact on water quantity and quality. In this review, we provide an overview of the methods developed over the past several decades to investigate GW–SW interactions. These methods include geophysical, hydrometric, and tracer techniques, as well as various modeling approaches. Different methods reveal valuable information on GW–SW interactions at different scales with their respective advantages and limitations. Interpreting data from these techniques can be challenging due to factors like scale effects, heterogeneous hydrogeological conditions, sediment variability, and complex spatiotemporal connections between GW and SW. To facilitate the selection of appropriate methods for specific sites, we discuss the strengths, weaknesses, and challenges of each technique, and we offer perspectives on knowledge gaps in the current science.

Annual Review of Environment and Resources↗