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 883 records · Page 49Linked to original sources
USGS2024

Hyperspectral (VNIR-SWIR) analysis of roll front uranium host rocks and industrial minerals from Karnes and Live Oak Counties, Texas Coastal Plain

VNIR-SWIR (400–2500 nm) reflectance measurements were made on the surfaces of various cores, cuttings and sample splits of sedimentary rocks from the Tertiary Jackson Group, and Catahoula, Oakville and Goliad Formations. These rocks vary in composition and texture from mudstone and claystone to sandstone and are known host rocks for roll front uranium occurrences in Karnes and Live Oak Counties, Texas. Spectral reflectance profiles, 569 in total, were reduced to 125 representative spectral signatures, which were analyzed using the U.S. Geological Survey's (USGS) Material Identification and Characterization Algorithm (MICA). MICA uses an automated continuum-removal procedure together with a least-squares linear regression to determine the fit of observed sample spectral absorption features to those of reference mineral standards in a spectral library. The reference minerals include various clay, mica, carbonate, ferric and ferrous iron minerals and their mixtures. In addition, absorption feature band-depth analysis was done to identify rock surfaces exhibiting absorption features related to uranium and zeolite minerals, which were not included in the command files used to execute MICA. Rocks from each of the four geologic units produced broadly similar spectral signatures as a result of comparable mineral compositions, but there were some notable differences. For example, Ca- and Na-montmorillonite was matched most frequently to the spectral absorption features in 2-μm (∼2000–2500 nm) wavelengths, while goethite occurred often at 1-μm (∼400–1000 nm) wavelengths. The latter is related to limonitic iron-staining in and around oxidized zones of the uranium roll front as described in previous papers. Rocks of the Jackson Group differed from those of the Catahoula, Oakville and Goliad units in that the former exhibited spectral features we interpret as being due to the presence of lignite-bearing mudstone layers. Goliad rocks exhibit spectral features related to dolomite, gypsum, anhydrite, and an unidentified green clay mineral that is possibly glauconite . Jackson Group rocks also exhibit weak but well-resolved absorption features at 964 and 1157 nm related to either or both zeolite minerals clinoptilolite and heulandite. These zeolite minerals and a few spectra exhibiting hydrous silica absorption features are indicative of alteration of volcanic glass in tuffaceous mudstone and claystone layers. A few sample spectra exhibited strong absorption features at around 1135 nm related to the uranium mineral coffinite. Both the 1135 nm coffinite and 1157 nm zeolite absorption features overlap somewhat, potentially making them difficult to distinguish without additional hyperspectral field, laboratory or remote sensing data. The results of this study were compared to mixtures of minerals described for ore, gangue and alteration minerals in deposit models for sandstone-hosted uranium, sedimentary bentonite and sedimentary zeolite. Use of these spectra can help facilitate mapping of both waste materials from the legacy mining of the above commodities, as well as future exploration and resource assessment activities.

Texas↗
USGS2024

The addition of 144Nd atomic mass to routine ICP-MS analysis as a Quick Screening Tool for Approximating Rare Earth Elements (Q-STAR) in natural waters

Rare earth elements (REEs) are a class of critical minerals, all of which can have supply chain vulnerability that impacts economic security. These elements are widely measured in environmental matrices via inductively coupled plasma mass spectrometry (ICP-MS); however, successful quantification can require time-consuming, sample-specific optimization. While a sample-by-sample approach is appropriate for targeted quantification studies, this approach is not suitable for mineral exploration efforts where rapidly screening thousands of samples for the presence of REEs is desired. Here, we demonstrated the use of a Quick Screening Tool for Approximating REEs (Q-STAR) to detect REEs in surface water and groundwater matrices, collected as part of existing environmental studies. A mass-to-charge ratio of 144 ( m / z = 144) was added to an ICP-MS method to screen for REEs in filtered water samples submitted for metals analyses to the U.S. Geological Survey (USGS) National Water Quality Laboratory. We detected the presence of REEs above a reference threshold of 1200 counts per second in 18 % of pre-selected 6626 samples. Using this screened dataset, we mapped estimated dissolved REE concentrations across the United States in relation to ecoregions and underlying geology. Data are constrained to where sample collection took place but nevertheless show estimated aqueous dissolved REE concentrations on a geographic scale that has not yet been studied. To validate Q-STAR, REEs were measured in a USGS standard reference sample, a subset of 88 archived filtered water samples, and in fresh filtered surface water samples. Our targeted analyses demonstrated a strong linear relationship between Q-STAR predicted and measured values in all archived samples for Nd (r 2 = 0.94), and light REEs (LREEs) such as lanthanum (La) (r 2 = 0.93), praseodymium (Pr) (r 2 = 0.94) and samarium (Sm) (r 2 = 0.94). Using Q-STAR screen values, nine field sites were identified and surface water samples recollected to confirm the continued presence of Nd and LREEs. Q-STAR can be used to screen an unlimited number of water samples for the presence of REEs prior to time-intensive and costly quantitative analyses and to generate large REE datasets for further investigation.

Journal of Geochemical Exploration↗
USGS2024

ShakeAlert® and schools: Incorporating earthquake early warning in school districts in Alaska, California, Oregon, and Washington

The U.S. Geological Survey-managed ShakeAlert® earthquake early warning system is the first public alerting system in the United States to provide rapid mass notification when an earthquake is detected. Although public alert delivery via mobile phones began in California in 2019 followed by Oregon and Washington in 2021, little is known about what might drive widespread implementation in at-risk institutional settings such as schools. For example, there has been limited research on how to best integrate earthquake early warning into existing emergency plans, alert systems, and drills to keep school children and personnel safe in an earthquake. To address this gap, in the spring of 2022, every school district superintendent in Alaska, California, Oregon, and Washington was sent a 15-min online survey. The survey assessed superintendent knowledge of ShakeAlert, preferences for alert messaging, and perceived opportunities and barriers to incorporating the system in schools. The results showed that superintendents had low awareness of ShakeAlert but held positive perceptions of the system's potential to enable life-saving protective actions. A major barrier to adoption included the perceived financial cost of implementing and maintaining the system. There were some statistically significant differences in state responses, and future research could investigate the specific needs of each state based on school district size and composition, hazard exposure, and earthquake experience. Together these findings can help inform targeted strategies to increase ShakeAlert adoption in schools and ultimately improve the safety of school children and staff during earthquakes.

Alaska, California, Oregon, Washington↗
USGS2024

Will there be water? Climate change, housing needs, and future water demand in California

Climate change in California is expected to alter future water availability, impacting water supplies needed to support future housing growth and agriculture demand. In groundwater-dependent regions like California's Central Coast, new land-use related water demand and decreasing recharge is already stressing depleted groundwater basins. We developed a spatially explicit state-and-transition simulation model that integrates climate, land-use change, water demand, and groundwater gain-loss to examine the impact of future climate and land use change on groundwater balance and water demand in five counties along the Central Coast from 2010 to 2060. The model incorporated downscaled groundwater recharge projections based on a Warm/Wet and a Hot/Dry climate future from a spatially explicit hydrological process-based model. Two urbanization projections from a parcel-based, regional urban growth model representing 1) recent historical and 2) state-mandated housing growth projections were used as alternative spatial targets for future urban growth. Agricultural projections were based on recent historical trends from remote sensing data. Annual projected changes in groundwater balance were calculated as the difference between land-use related water demand, based on historical estimates, and climate-driven recharge plus agriculture return flows. Results indicate that future changes in climate-driven groundwater recharge, coupled with cumulative increases in agricultural water demand, result in overall declines in future groundwater balance, with a Hot/Dry future resulting in cumulative groundwater decline in all but Santa Cruz County. Cumulative declines by 2060 are especially prominent in San Luis Obispo (−2.9 to −5.1 Bm 3 ) and Monterey counties (−6.5 to −8.7 Bm 3 ), despite limited changes in agricultural water demand over the model period. These two counties show declining groundwater reserves in a Warm/Wet future as well, while San Benito and Santa Barbara County barely reach equilibrium. These results suggest future groundwater supplies may not be able to keep pace with regional demand and declining climate-driven recharge, resulting in a potential reduction in water security in the region. However, our county-scale projections showed new housing and associated water demand does not conflict with California's groundwater sustainability goals. Rather, future climate coupled with increasing agricultural groundwater demand may reduce water security in some counties, potentially limiting available groundwater supplies for new housing.

California↗
USGS2024

Indications of preferential groundwater seepage feeding northern peatland pools

Groundwater seepage from underlying permeable glacial sedimentary structures, such as eskers, has been hypothesized to directly feed pools in northern peat bogs. These hypotheses directly contradict classical peat bog models for ombrogenous systems, wherein meteoric water is the sole water input to these systems. Variations in the underlying mineral sediment in contact with the peat imply that unrecognized hydrogeologic connectivity may exist with pools in northern peat bogs, particularly where high permeability materials are in contact with the peat. Seepage dynamics originating from these structural variations were investigated using a suite of thermal and hydrogeophysical methods deployed around pools in a peat bog of northeastern Maine, USA. Thermal characterization methods mapped anomalies that were confirmed as matrix seepage or preferential flow pathways (PFPs). Geochemical methods were employed at identified thermal anomalies to confirm upwelling of solute-rich groundwater. Conduits around pools were associated with surficial terminations of suspected peat pipes, based on the inference of pathways extending down into the peat, that focus flow through PFPs in the peat matrix. Discharge also occurred through the peat matrix adjacent to suspected pipe structures and matrix seepage rates were quantified using analysis of diurnal temperature signals recorded at multiple depths. Seepage rates, with a maximum of nearly 0.4 m/d, were measured at localized points around pools. Periods of synchronized temperatures paired with highly muted diurnal temperature signals, recorded in diurnal temperature with depth data, were interpreted qualitatively as activation of strong upward discharge rates through suspected peat pipes. These time periods correlated strongly with local precipitation events around the peatland. Ground-penetrating radar surveys revealed discontinuities in the low permeability glacio-marine clay at the mineral sediment-peat interface, interpreted to be regional glacial esker deposits, which were located beneath and around pools. Heat tracing, specific conductance contrasts, seepage rates, and trace metal concentrations all imply groundwater seepage originating from underlying permeable glacial esker deposits and directly sourcing pools. Preferential groundwater inputs into northern peat bogs may play a key role in developing and maintaining pool systems, with enhanced solute transport impacting peatland ecology, water resources, and carbon cycling.

Maine↗
USGS2024

Diel temperature signals track seasonal shifts in localized groundwater contributions to headwater streamflow generation at network scale

Groundwater contributions to streamflow sustain aquatic ecosystem resilience; streams without significant groundwater inputs often have well-coupled air and water temperatures that degrade cold-water habitat during warm low flow periods. Widespread uncertainty in stream-groundwater connectivity across space and time has created disparate predictions of energy and nutrient fluxes across headwater networks, hindering predictions of cold-water habitat resilience under climate change scenarios. Recently, annual paired air and water temperature signals have been harnessed to indicate stream water thermal sensitivity and the dominance of deep versus shallow groundwater influence, although the utility of diel air–water temperature signal metrics for hydrologic inference has remained unexplored. Here we analyzed two consecutive years of locally paired, air–water temperature data from 47 headwater stream sites in the Catskill Mountains, New York, USA, and discovered characteristic seasonal patterns in diel temperature signal sinusoid metrics (amplitude ratio, phase lag, and mean ratio) driven by shifts in streamflow generation mechanisms and stream network position. Hydrologic interpretations of observed patterns were supported by stream heat budget model scenarios and additional analysis of paired air–water temperature data from two streams in Shenandoah National Park, Virginia, USA, with well characterized stream-groundwater connectivity. We found that within smaller tributaries, streamflow generation transitions from runoff to groundwater dominance were driven by hillslope drying during seasonal periods of lower precipitation. This was evidenced by significant correlations (p < 0.01) between daily water:air temperature signal amplitudes (non-linear decreases of ∼ 50 %) and derived base-flow index at 22 of the 28 sites, indicating enhanced local groundwater influence on streamflow promotes decoupling of diel air–water temperature signals. Additionally, ratios between daily water:air temperature signal means were lower in tributaries (∼0.68) when compared to main-stem (∼0.8) sites, increasing linearly throughout the observational period. In conceptual stream heat budget models, groundwater inflow had minimal effects on daily phase lags (∼0.2 hr), but increases in fractional groundwater discharge (0–50 %) depressed daily amplitude (∼20 % to 50 %) and mean ratios (∼15 %), supporting the sensitivity of daily metrics to interpreted changes in seasonal groundwater contributions to streamflow. During observational periods (i.e., April through October 2021 and 2022), significant differences (p < 0.01) between tributary and main-stem air–water metrics occurred when base-flow contributions were highest (∼0.93 vs. ∼ 0.68), as sites lower in the network had daily temperature metrics dominated by stream channel thermal inertia, rather than local groundwater connectivity, showing enhanced air–water diel signal coupling during warmer, drier periods. Divergent air temperature coupling across the network was interpreted as being driven by distance from local groundwater source zones, additional lateral groundwater inflows do not contribute a meaningful fraction to channel discharge lower in the network. Given the growing footprint of stream temperature observations, diel air–water temperature signals can provide distributed metrics sensitive to upstream groundwater discharge. Consequently, these metrics can support ongoing efforts by resource managers and researchers seeking to forecast the resilience of cold-water habitat to climate warming and changing precipitation regimes in mountain headwater streams.

New York, Virginia↗
USGS2024

High resolution identification and quantification of diffuse deep groundwater discharge in mountain rivers using continuous boat-mounted helium measurements

Discharge of deeply sourced groundwater to streams is difficult to locate and quantify, particularly where both discrete and diffuse discharge points exist, but diffuse discharge is one of the primary controls on solute budgets in mountainous watersheds. The noble gas helium is a unique identifier of deep groundwater discharge because groundwater with long residence times is commonly enriched in helium. In this study, a portable mass spectrometer was used to measure longitudinal variation in dissolved helium concentrations in two mountainous rivers at high spatial resolution not feasible with traditional sampling techniques. Helium profiles were then simulated using a mass-balance model to quantify longitudinal variation in groundwater discharge to the receiving rivers. Results indicate helium concentrations were enriched by multiple orders of magnitude above atmospheric equilibrium in both rivers and that this persisted for up to 18 km below observed pulse inputs in the Colorado River. Helium mass-balance models match observed longitudinal patterns with the exception of sharp initial increases in helium observed in the rivers. Increased longitudinal groundwater discharge rates correspond to mapped geologic structures in both watersheds that likely transport deep geothermal water. Models show variable sensitivity to spatial assignment of input variables representing the groundwater source, illustrating the importance of collecting data from discrete groundwater discharges where possible. The methodology shows promise for field experiments designed to assess air–water exchange rates and to quantify total groundwater discharge from a combination of discrete and diffuse sources.

Colorado, Utah↗
USGS2024

Event classification, seismicity, and eruption forecasting at Great Sitkin Volcano, Alaska: 1999–2023

The frequency content of volcanogenic seismicity is often used to classify events and their spatial and temporal progression is then used to map subsurface volcanic processes. The progression of volcano-seismic events and associated source processes also plays a critical role in eruption forecasting. Here we develop and evaluate a computerized methodology for characterizing volcano-seismic event types using Frequency Index and Average Peak Frequency. We apply and test this technique at Great Sitkin Volcano, Alaska, classifying over 9000 hypocenters between 1999 and 2023. This 24-year time span covers periods of seismic quiescence, earthquake activity on nearby tectonic (bookshelf) faults, precursory unrest from 2016 to 2021, and the explosive onset in May 2021 of the ongoing effusive eruption. We use the spatial and temporal evolution of classified event types to map the active volcanic and tectonic processes, develop a conceptual model of the subsurface magmatic system, and perform a retrospective analysis of eruption forecasts at Great Sitkin Volcano between 2016 and the present. The classification and progression of hypocenters suggests the subsurface Great Sitkin Volcano magmatic system consists of a mid- to lower- crustal source zone between 10 and 40 km depth and an upper crustal magma storage area between −1 and 10 km depth (hypocenter depth is referenced to sea level and negative depths reflect height above sea level). The earliest precursors occurred in July 2016 and consisted of deep long-period and volcano-tectonic earthquakes at mid-crustal depths suggesting the subsequent unrest and eruption were triggered by a deeper intrusion of magma. This mid-crustal seismic activity was immediately followed by the onset upper-crustal long-period events and volcano-tectonic earthquakes VTs suggesting a strong linkage between the shallow and deeper portions of the magmatic system. The upper crustal area was likely capped by the 1974 lava dome until the magmatic explosion on May 26, 2021.

Alaska↗
USGS2024

Formation of vertical columnar seismic structures and seafloor depressions by groundwater discharge in the drowned Miami Terrace platform and overlying deep-water carbonates, southeastern Florida

The presence of vertical cross-formational fluid migration passageways within sedimentary basins can profoundly impact aquifer and reservoir fluid-flow and their identification is fundamental to informing management of subsurface fluid resources (groundwater, oil, gas). In an onshore and offshore southeastern part of Florida, 2D/3D seismic-reflection and bathymetry data document ∼153 vertical columnar structures composed of reflection disruptions up to 790 m in the height and averaging 360 m in diameter, and ∼219 subcircular to circular seafloor depressions up to 1334 m wide. Our study focuses on these features found within the offshore shallow-marine carbonate Miami Terrace platform, which drowned approximately at the end of the middle Miocene, and within overlying Plio-Quaternary deep-water carbonate slope and drift deposits. Most columnar structures are rooted in stratiform aquifers of the Miami Terrace platform and associated with faults or fault intersections produced by Eocene and circa late Miocene tectonics. The columns commonly terminate within the platform or as subcircular depressions along an amalgamated karstic and drowning unconformity at the platform top. The columns typically stretch upwards from a zone of deep karst cavity collapse through the Miami Terrace platform with upward decreasing sag on internal reflections. Following drowning and Plio-Quaternary partial burial of the Miami Terrace platform by deep-water deposits, the subcircular depressions and faults along the platform top were points of origin for a second phase of column growth upward into the deep-water deposits. The continuation of deep platform cavity collapse and column evolution produced pockmarks along paleo-seafloors within the deep-water deposits and at the present-day sea floor. The Plio-Quaternary pockmarks formed at water depths too deep to suggest an origin related to meteoric karst above or near sea level, but rather their formation is suggested to be related to cyclic sea level falls that drove increased groundwater head and density gradients, and seafloor discharge of offshore freshened groundwater sourced from the underlying platform. Plausibly, mixing of freshened groundwater and seawater at the seafloor discharge sites drove dissolution of the host deep-water deposits, which together with erosion by groundwater venting and current scouring formed the pockmarks. Seaward of the Plio-Quaternary seafloor pockmarks, at the late-middle Miocene upper slope of the Miami Terrace platform and along the regional karst/drowning unconformity is a slope-parallel band of ∼189 densely distributed subcircular seafloor depressions with diameters up to 1334 m at water depths up to ∼660 m. It is plausible that along the upper slope, faults and fractures produced by gravity-driven slope instability and possibly tectonics formed a dense network of fluid passageways that promoted upward artesian freshened groundwater flow to sites of discharge where mixing with seawater generated limestone dissolution and the depressions. But tectonic uplift may have forced emersion and initial meteoric sinkhole formation circa late Miocene with later enhancement by freshened groundwater discharge and bottom current erosion.

Florida↗
USGS2024

Mapping the resistivity structure of Walker Ridge 313 in the Gulf of Mexico using the marine CSEM method

A marine controlled source electromagnetic (CSEM) campaign was carried out in the Gulf of Mexico to further develop marine electromagnetic techniques in order to aid the detection and mapping of gas hydrate deposits. Marine CSEM methods are used to obtain an electrical resistivity structure of the subsurface which can indicate the type of substance filling the pore space, such as gas hydrates which are more resistive. Results from the Walker Ridge 313 study (WR 313) are presented in this paper and compared with the Gulf of Mexico Gas Hydrate Joint Industry Project II (JIP2) logging while drilling (LWD) results and available seismic data. The hydrate, known to exist within sheeted sand deposits, is mapped as a resistive region in the two dimensional (2D) CSEM inversion models. This is consistent with the JIP2 LWD resistivity results. CSEM inversions that use seismic horizons provide more realistic results compared to the unconstrained inversions by providing sharp boundaries and architectural control on the location of the resistive and conductive regions in the CSEM model. The seismic horizons include: 1) the base of the gas hydrate stability zone (BGHSZ), 2) the top of salt, and 3) the top and bottom of a fine grained marine mud interval with near vertical hydrate filled fractures, to constrain the CSEM inversion model. The top of salt provides improved location for brines, water saturated salt, and resistive salt. Inversions of the CSEM data map the occurrence of a ‘halo’ of conductive brines above salt. The use of the BGHSZ as a constraint on the inversion helps distinguish between free gas and gas hydrate as well as gas hydrate and water saturated sediments.

Louisiana↗
USGS2024

Factors contributing to pesticide contamination in riverine systems: The role of wastewater and landscape sources

Wastewater treatment plant (WWTP) discharges can be a source of organic contaminants, including pesticides, to rivers. An integrated model was developed for the Potomac River watershed (PRW) to determine the amount of accumulated wastewater percentage of streamflow (ACCWW) and calculate predicted environmental concentrations (PECs) for 14 pesticides in non-tidal National Hydrography Dataset Plus Version 2.1 stream segments. Predicted environmental concentrations were compared to measured environmental concentrations (MECs) from 32 stream sites that represented a range of ACCWW and land use to evaluate model performance and to assess possible non-WWTP loading sources. Statistical agreement between PECs and MECs was strongest for insecticides, followed by fungicides and herbicides. Principal component analysis utilizing optical fluorescence and ancillary water quality data identified wastewater and urban runoff sources. Pesticides that indicated relatively larger sources from WWTPs included dinotefuran, fipronil, carbendazim, thiabendazole, and prometon whereas imidacloprid, azoxystrobin, propiconazole, tebuconazole, and diuron were more related to urban runoff. In addition, PECs generally comprised a low proportion of MECs, which indicates possible dominant loading sources beyond WWTP discharges. Cumulative potential toxicity was higher for sites with greater ACCWW and/or located in developed areas. Imidacloprid, fipronil, and carbendazim accounted for the largest portion of predicted potential toxicity across sites. The chronic aquatic life toxicity benchmarks for freshwater invertebrates were exceeded for 82 % of the imidacloprid detections ( n = 28) and 47 % of the fipronil detections ( n = 19). These results highlight the ecological implications of pesticide contamination from WWTP discharges and also the potential legacy effects from accumulated soil and groundwater sources. Pesticide management strategies that mitigate both current and historical impacts may improve the health of aquatic ecosystems.

Potomac River watershed↗
USGS2024

Produced water geochemistry from hydraulically stimulated Niobrara Formation petroleum wells: Origin of salinity and temporal perspectives on treatment and reuse

Produced water (i.e., a mixture of returned injection fluids and geologic formation brines) represents the largest volumetric waste stream associated with petroleum production in the United States. As such, produced water has been the focus of intense study with emphasis on understanding the geologic origin of the fluids, environmental impacts of unintended or intentional release, disposal concerns, and their commodity (e.g., lithium) potential. However, produced water geochemistry from many active petroleum plays remain poorly constrained leading to knowledge gaps associated with the origin of brine salinity and parameters (e.g., radium levels) that can impact treatment, disposal, and possible reuse. Here we evaluate the major ion geochemistry, radium concentrations, and stable water isotope composition of ~120 produced water samples collected from 17 producing unconventional petroleum wells in Weld County, Colorado from the Late Cretaceous Niobrara Formation. This sample set encompasses eight produced water time series from four new wells across production days 0 to ~365 and from four established wells across production days ~1000 to ~1700. Additionally, produced water from nine other established Niobrara Formation wells were sampled at discrete time points ranging from day 458 to day 2256, as well as hydraulic fracturing input fluids. These results expand the available Niobrara Formation produced water geochemical data, previously limited to few wells sampled within the first year of production, allowing for the heterogeneity of major ions and radium to be evaluated. Furthermore, we explore the geochemical relationships between major ion ratios and stable water isotope composition to understand the origin of salinity in Niobrara Formation brines from the Denver-Julesburg Basin. These findings are discussed with perspective toward potential treatment and reuse of Niobrara produced water prior to disposal.

Colorado, Nebraska, Wyoming↗
USGS2024

Geophysical characterization of an alkaline‑carbonatite complex using gravity and magnetic methods at Magnet Cove, Arkansas, USA

The Magnet Cove alkaline‑carbonatite complex (MCC), located in the Ouachita Mountains of south-central Arkansas in the United States, hosts an extensive variety of rare rock types and critical mineral resources with physical properties (density and magnetization) that contrast significantly with the sedimentary rocks into which they have intruded. Newly acquired ground-based gravity and magnetic data were used to develop two-dimensional and three-dimensional geophysical models of the Cretaceous-aged Magnet Cove intrusive complex. The models reveal that the MCC: (1) widens out at middle crustal depths to as much 22 km across, and may reach a depth of 20 km; (2) has a total volume (exposed and subsurface) that may be over 800 km 3 ; (3) is likely connected at depth to other intrusions in the Arkansas alkaline province; and (4) has a geometry that is aligned with pre-existing structures such as the Reelfoot rift and the Ouachita orogenic belt, some of which were likely structurally controlled by the Precambrian crystalline basement and the continent-ocean transition zone buried beneath the Ouachita orogen. For the first time, the magnetic models of the MCC account for the presence of strong remanent magnetization. This results in a geophysical workflow necessary to accurately interpret magnetic anomalies over the much larger Arkansas alkaline province, its geologic and structural framework, and critical mineral potential.

Arkansas↗
USGS2024

Subducting plate structure and megathrust morphology from deep seismic imaging linked to earthquake rupture segmentation at Cascadia

The origin of rupture segmentation along subduction zone megathrusts and linkages to the structural evolution of the subduction zone are poorly understood. Here, regional-scale seismic imaging of the Cascadia margin is used to characterize the megathrust spanning ~900 km from Vancouver Island to the California border, across the seismogenic zone to a few tens of kilometers from the coast. Discrete domains in lower plate geometry and sediment underthrusting are identified, not evident in prior regional plate models, which align with changes in lithology and structure of the upper plate and interpreted paleo-rupture patches. Strike-slip faults in the lower plate associated with oblique subduction mark boundaries between regions of distinct lower plate geometry. Their formation may be linked to changes in upper plate structure across long-lived upper plate faults. The Juan de Fuca plate is fragmenting within the seismogenic zone at Cascadia as the young plate bends beneath the heterogeneous upper plate resulting in structural domains that coincide with paleo-rupture segmentation.

Oregon, Washington↗
USGS2024

Evaluation of the lithium resource in the Smackover Formation brines of southern Arkansas using machine learning

Global demand for lithium, the primary component of lithium-ion batteries, greatly exceeds known supplies, and this imbalance is expected to increase as the world transitions away from fossil fuel energy sources. High concentrations of lithium in brines have been observed in the Smackover Formation in southern Arkansas (>400 milligrams per liter). We used published and newly collected brine lithium concentration data to train a random forest machine-learning model using geologic, geochemical, and temperature explanatory variables and create a map of predicted lithium concentrations in Smackover Formation brines across southern Arkansas. Using these predicted lithium maps with reservoir parameters and geologic information, we calculated that there are 5.1 to 19 million tons of lithium in Smackover Formation brines in southern Arkansas, which represents 35 to 136% of the current US lithium resource estimate. Based on these calculations, in 2022, 5000 tons of dissolved lithium were brought to the surface within brines as waste streams of the oil, gas, and bromine industries.

Arkansas↗
USGS2024

Predictions of groundwater PFAS occurrence at drinking water supply depths in the United States

Per- and polyfluoroalkyl substances (PFAS), known colloquially as “forever chemicals”, have been associated with adverse human health effects and have contaminated drinking water supplies across the United States owing to their long-term and widespread use. People in the United States may unknowingly be drinking water that contains PFAS because of a lack of systematic analysis, particularly in domestic water supplies. We present an extreme gradient boosting model for predicting the occurrence of PFAS in groundwater at the depths of drinking water supply for the conterminous United States. Our model results indicate that 71 to 95 million people in the conterminous United States potentially rely on groundwater with detectable concentrations of PFAS for their drinking-water supplies prior to any treatment.

conterminous United States↗
USGS2024

Mesoproterozoic to Paleozoic tectonics, Pleistocene landforms, and Holocene seismicity in the Blue Ridge: Results from integrated studies of the 9 August 2020, Mw 5.1 earthquake area near Sparta, North Carolina, USA

This field trip examines the results of integrated geologic studies of the 9 August 2020, M w 5.1 earthquake near Sparta, North Carolina, USA. The earthquake generated ~4 km of coseismic surface rupture of the Little River fault and uplifted a surface area of ~11 km 2 . The Little River fault is a thrust fault oriented 110–130°/45–70°SW, and mapped fault segments are en echelon with scarp heights from <5–30 cm. The epicenter is in polydeformed rocks of the Ashe and Alligator Back Metamorphic Suites in the eastern Blue Ridge. Bedrock structure formed during multiple Paleozoic orogenies; the regional foliation strikes NE-SW and dips SE (mean orientation 063°/52°SE). Mapping identified late Paleozoic veins and shear zones, a regional joint set striking 330–340° and 250–240°, and brittle faults that cut the Paleozoic foliation. Brittle faults oriented similar to the Little River fault are mapped up to 4 km along strike from the coseismic rupture along Bledsoe Creek valley, and the combined length of the Little River fault system is ~8 km. Paleoseismic trenches across the Little River fault corroborate the reactivation of an older fault by the 2020 earthquake and reveal two events during late Pleistocene (<50 ka). Surficial mapping identified several terrace deposits, including a deposit along Bledsoe Creek that yielded a 26 Al/ 10 Be isochron burial age of 0.46 ± 0.13 Ma and overlies a brittle fault, thus constraining the timing of movement of the fault at that location. Paleoliquefaction studies document soft-sediment deformation features in alluvium that may represent paleoseismic events. Collectively, these results highlight long-lived paleoseismicity of the Blue Ridge and that the 9 August 2020 earthquake reactivated an older, suitably oriented brittle fault in the bedrock. The Little River fault is an example of a previously unknown but active fault lying outside of known seismic zones with demonstrated recurrence of paleo-ruptures, raising questions about the assumption that damaging earthquakes are limited to areas of ongoing background seismicity, which is counter to seismic hazard assessments in the eastern United States. Bedrock mapping separates eastern Blue Ridge lithostratigraphy of the Lynchburg Group and Ashe and Alligator Back Metamorphic Suites into separate fault-bound packages juxtaposed over various 1.3–1.0 Ga basement rocks of the northern French Broad massif by the Gossan Lead fault.

North Carolina↗
USGS2024

Unscrambling the Proterozoic supercontinent record of northeastern Washington State, USA

The time interval from Supercontinent Nuna assembly in the late Paleoproterozoic to Supercontinent Rodinia breakup in the Neoproterozoic is considered by some geologists to comprise the “Boring Billion,” an interval possibly marked by a slowdown in plate tectonic processes. In northeastern Washington State, USA, similar to much of western Laurentia, early workers generally thought the tectonostratigraphic framework of this interval of geologic time consisted of two major sequences, the (ca. 1480–1380 Ma) Mesoproterozoic Belt Supergroup and unconformably overlying (<720 Ma) Neoproterozoic Windermere Supergroup. However, recent research indicates that strata considered by early workers as Belt Supergroup equivalents are actually younger, and a post-Belt, pre-Windermere record is present within the <1360 Ma Deer Trail Group and <760 Ma Buffalo Hump Formation. Thus, the northeastern Washington region perhaps comprises the most complete stratigraphic record of the “Boring Billion” time interval in the northwestern United States and holds important insights into global Proterozoic supercontinent tectonic processes. In light of these exciting developments, this field guide will address the early historic economic geology and original mapping of these Proterozoic sequences in the northeastern Washington region, and from that foundation explore more recent isotopic provenance data and their regional to global context. Finally, the guide will end with a discussion of remaining questions with a goal of stimulating interest in these relatively understudied, yet important, rocks.

Washington↗