Geology Reports⌕ Search

SEARCH · Geology Reports

Results for “Journal of Agricultural, Biological and Environmental Statistics”

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,689 records · Page 53Linked to original sources
USGS2025

Global patterns of coseismic landslide runout mobility differ from aseismic landslide trends

Coseismic landslides significantly contribute to human and economic losses during and immediately following earthquakes, yet very little data on the runout of such landslides exist. While well-established behavior of aseismic (e.g., hydrologically triggered) landslide runout mobility suggests strong correlation between landslide size and mobility, limited studies of coseismic landslide runout find conflicting mobility trends. We present a global dataset of runout lengths produced from a new automated method for estimating landslide runout, developed and validated using 1726 manually mapped landslides from five unique earthquakes. We then apply the automated runout tool to 23 global earthquake-induced landslide inventories, producing a compiled database of 73,665 measured and estimated runout lengths of coseismic landslides to assess mobility trends. We find a significant divergence between well-established aseismic mobility trends and that of coseismic landslides, with far greater scatter and more complex mobility patterns in earthquake-triggered landslides. As a function of landslide size, we observe global coseismic landslide mobility patterns are bilinear, becoming increasingly less mobile with increasing size above some threshold. This discordance between aseismic and coseismic landslide mobility may be a function of landslide type, kinematics, hydrology, and or setting that systematically differ between triggering mechanisms and should be explored in more depth to develop predictive models of these unique runout patterns. These results suggest hazard and risk models for coseismic landslides may significantly under-predict or over-predict impacts, depending on the size of triggered landslides.

Engineering Geology↗
USGS2025

PEST++IES how many iterations & realizations, finding the point of diminishing returns

PEST++IES (White 2018; White et al. 2020) is widely used in the groundwater modeling community for its ability to perform computationally efficient history matching and uncertainty analysis in a highly parameterized context. One primary advantage of using an iterative ensemble smoother is that the number of model runs required per iteration depends on the number of realizations in an ensemble, not the number of parameters in each realization. However, this raises the question: what is the optimal number of realizations and iterations to use for any one model before the point of diminishing returns? Using a modified version of the Freyberg model (Freyberg 1988; Hunt et al. 2020), different parameter and observation scenarios were evaluated for four iterations and ensembles of 10, 25, 50, 100, 250, 500, 1000, and 2000 realizations. To match observations, PEST++IES altered hydraulic conductivity ( k ), both globally across the model and locally at three different pilot point densities, as well as global recharge (via a single multiplier), global river conductance, and individual well flow rates. Risk-based well capture zone results (Fienen et al. 2022a) and estimated hydraulic conductivity fields from each scenario were quantitatively and qualitatively compared against the “truth” model and its outputs. Across the cases examined, ensemble sizes of 100 to 250 realizations and two PEST++IES iterations were generally sufficient to achieve good results.

Groundwater↗
USGS2024

Influence of inherited structure on flexural extension in foreland basin systems: Evidence from the northern Arkoma basin and southern Ozark dome, USA

Extensional faults are key components of foreland basin systems. They form within the upper crust in response to flexure of the lithosphere and accommodate subsidence within the foredeep and forebulge depozones. Such faults are excellent proxies for orogenic system evolution and control the distribution of natural resources and hazards. However, the spatiotemporal evolution of flexural extension has not been documented previously at a regional scale, thereby limiting our understanding of underlying geodynamic controls. Here, we resolve late Paleozoic flexural extension in the northern Arkoma basin and southern Ozark dome, USA. We synthesize a large database of previous mapping, existing research, subsurface data, and geophysical data into 3D geologic and 2D kinematic models. Mesh surfaces representing several key horizons from the Carboniferous Period (ca. 335-306 Ma) were constructed. These surfaces were built from oil and gas well tops (n = ∼10,000) and surface geologic map contacts using an advanced kriging method. The mesh surfaces are offset by a complex 3D fault network, allowing detailed analysis of along-strike and down-dip variations in fault displacement. Analysis of the 3D model reveals a regular and repeated fault segmentation pattern wherein E -W striking, left- and foreland-stepping en échelon normal faults are segmented by inherited NE striking basement faults. Maximum vertical separation along the E -W normal faults is generally focused between the inherited NE-trending faults. This suggests that the inherited basement faults delocalized extensional strain during late Paleozoic normal faulting. Maximum vertical separation and fault localization may correlate to areas with high-amplitude positive magnetic anomalies interpreted as Mesoproterozoic granitic rocks. Speculative covariance of magnetic anomalies and fault displacements implies that the relatively strong basement granite concentrated stress, leading to localized faulting within the relatively thin sedimentary cover. Lastly, we show that flexural extension migrated southeast to northwest from the Chesterian-Morrowan (ca. 335-319 Ma) to the Desmoinesian (ca. 306 Ma). The migratory flexural extension may be explained by diachronous loading during Pangean assembly, or by synchronous loading but variable load compensation due to inherent factors.

Arkansas, Oklahoma↗
USGS2024

Earthquake relocations delineate discrete a fault network and deformation corridor throughout Southeast Alaska and Southwest Yukon

Deformation in southeastern Alaska and southwest Yukon is governed by the subduction and translation of the Pacific-Yakutat plates relative to the North American plate in the St. Elias region. Despite notable historical seismicity and major regional faults, studies of the region between the Fairweather and Denali faults are complicated by glacial coverage and the remote setting. In the last decade, significant improvements have been made to the density of regional broadband seismometer networks. We relocate more than 5,000 earthquakes between 2010 and 2021 in the region of southeastern Alaska and southwestern Yukon utilizing these improved seismic networks. With reductions in catalog uncertainty, particularly in depth, we quantify the thickness of the seismogenic layer in the crust throughout the region and locate seismicity on a shallow network of upper-crustal faults. Relocated earthquakes, combined with an updated focal-mechanism catalog, permit estimating and classifying motion of active faults. This includes mapping the Totschunda-Fairweather “Connector” fault, which plays an important role in explaining regional deformation, and identifying new faults like the Kathleen Lake fault. We draw similarities between our seismic observations and simplified conceptual models of regional tectonics, which describe a dominant transpressional regime and localized slip partitioning. Our results support a hypothesis where current deformation is taking place on a well-defined and evolved network of shallow faults in the corridor between the Totschunda-Fairweather “Connector” and Denali faults.

Alaska, Yukon↗
USGS2025

U.S. Geological Survey—Department of the Interior, Region 11, Alaska—2023–24 biennial science report

Introduction U.S. Geological Survey (USGS) Mission—The USGS national mission is to monitor, analyze, and predict the current and evolving dynamics of complex human and natural Earth-system interactions and to deliver actionable information at scales and timeframes relevant to decision makers. Consistent with the national mission, the USGS in Alaska provides timely and objective scientific information to help address issues and inform management decisions across five interconnected focus areas: Energy and Minerals; Geospatial Mapping; Natural Hazards; Water Quality, Streamflow, and Ice Dynamics; and Ecosystems. The USGS in Alaska consists of approximately 350 scientists and support staff working in 3 Alaska-based science centers. USGS science activities are also initiated by the Cooperative Research Unit and USGS centers outside Alaska. In the last 5 years, USGS research in Alaska has produced many scientific benefits resulting from more than 900 publications. Publications relevant to Alaska can be conveniently searched by keyword through the USGS Publications Warehouse at https://pubs.usgs.gov/.

Alaska↗
USGS2024

Deep lithospheric controls on surface deformation and seismicity around the East Anatolian Fault Zone and A3 Triple Junction

The East Anatolian Fault Zone (EAFZ) is a plate-bounding strike-slip fault capable of hosting large earthquakes, as demonstrated by the extremely damaging February 2023 M w 7.8 and M w 7.7 mainshocks of the Kahramanmaraş earthquake sequence. Deformation related to this boundary, part of the Anatolia-Arabia-Africa (A 3 ) Triple Junction, is diffuse, as was shown by part of this earthquake sequence occurring on a northern splay of the EAFZ (the Sürgü-Çardak Fault Zone; SCFZ). Controls on surface deformation are commonly linked to stress in the brittle upper crust, but the complex deformation and seismicity patterns in this region may also reflect deeper processes, such as variations in the location and extent of the strong Arabian Plate lithospheric mantle. Seismic tomography indicates that the Arabian Plate underthrusts Anatolia as far north as the SCFZ and extends as far west as the central Adana Basin, coincident with a zone of relatively deep (>30 km) strike-slip seismogenesis that has produced M w >6 earthquakes. By investigating the relationship between deformation since the inception of the EAFZ (ca. 5 Ma), seismic structure, and seismicity, we infer that the SCFZ will become the future SE boundary of the Anatolian Plate as part of the evolving A 3 Triple Junction.

GSA Today↗
USGS2025

The influence of pre-existing structures on geothermal springs: Inferences from potential field mapping in Surprise Valley, CA and other sites In the northwestern Great Basin

Surprise Valley, located in the northwestern Great Basin, is an asymmetric extensional basin that marks a major tectonic transition between the relatively un-extended volcanic Modoc Plateau to the west, and the Basin and Range to the east that has undergone 10-15% extension. In addition, it sits just north of the Walker Lane which accommodates up to 20% of dextral slip associated with Pacific-North American plate interactions. Thermal springs issue from eight areas within Surprise Valley. Most of these occur within the basin and are not situated on the main basin forming range-front faults. As a result, efforts to resolve the structural setting of the valley’s hydrothermal system have relied on geophysics to characterize basin structure and geology. Extensive efforts to map the basin with ground and airborne magnetics have revealed a >35 km-long linear, intra-basin magnetic high, interpreted as a buried dike swarm. Geothermal springs on the eastern side of the valley, including Seifert hot springs, Leonards hot springs, and Surprise Valley hot springs (SVHS), are all situated along the magnetic high and occur at local breaks and bends in the anomaly, suggesting that fracture permeability is enhanced along the feature and particularly at these discontinuities. Recent studies, including drilling over the anomaly near SVHS that likely intersected dike material, as well as subsequent mapping and sampling of dikes outcropping along the anomaly on the playa surface south of SVHS, confirm (as previously inferred) that mafic intrusives are the principal source of the anomaly. Similar interpretations made in two other valleys (in southern Oregon and northwestern Nevada), where inferred intra-basin dikes appear to be spatially correlated with hot springs or prospective geothermal resource areas, suggest that the impact of pre-existing basement structure on hydrothermal activity may pertain more generally to other hydrothermal settings throughout the Great Basin. If so, efforts to map basement may enhance understanding structural controls on some geothermal systems. Furthermore, similarities across these disparate sites suggest that magmatism may play a much larger role in accommodating extension and influencing basin evolution across the western Great Basin than previously recognized.

Conference Paper↗
USGS2025

Scientific opportunities in the National Landscape Conservation System

The National Landscape Conservation System consists of unique and beautiful places across America’s landscapes where identified resources and values are protected and science is highlighted. The mission of the National Landscape Conservation System (NLCS), which is managed by the Bureau of Land Management and is often referred to as the agency’s National Conservation Lands, is to conserve, protect, and restore nationally significant landscapes for their cultural, ecological, and scientific values. This clear inclusion of science in the NLCS mission sets the stage for individual units to serve as places of learning, teaching, discovery, and innovation. Science is an integral part of managing the National Conservation Lands, and science conducted within and across the more than 900 units that make up the NLCS can inform and influence conservation and public land management well beyond its boundaries. Here, we highlight seven core aspects of National Conservation Lands that present valuable science opportunities: (1) the scientific values for which individual units are designated; (2) the many other resources, objects, and values within units; (3) the value of units as “control” sites for understanding the effects of activities such as mineral extraction that commonly occur elsewhere on multiple-use public lands but are often prohibited within National Conservation Lands; (4) the value of units for studying the effects of activities such as recreation that regularly occur and may be intensified on National Conservation Lands; (5) the high visibility of units, which draws strong interest and engagement from scientists, partners, and the public; (6) the functioning of the units as a network managed for a common purpose, which provides an opportunity to explore cross-cutting science questions across widely varying contexts and geographies; and (7) the opportunities units provide to promote and apply Indigenous Knowledge to scientific research to manage natural and cultural resources. Because of all of these characteristics, National Conservation Lands can serve as hubs for basic and applied science that can inform management of all public lands and resources into the future. We highlight these science opportunities through examples from existing units and suggest two actions that could help further science activities and impact on National Conservation Lands.

Parks Stewardship Forum↗
USGS2025

Structural setting and geothermal potential of northeastern Reese River Valley, north-central Nevada: Highly prospective detailed study site for the INGENIOUS project

The northeastern part of the Reese River basin situated ~15 km southeast of Battle Mountain, Nevada, scored highly in the Nevada geothermal play fairway analysis (PFA) for hosting potential hidden geothermal systems. This site (also referred to as Argenta Rise) was therefore chosen for detailed study in the INGENIOUS project (INnovative Geothermal Exploration through Novel Investigations Of Undiscovered Systems). The high PFA scores resulted primarily from favorable structural settings (e.g., fault intersections and pull aparts) with relatively high slip rates on Quaternary faults. The INGENIOUS project is utilizing additional parameters and more rigorous analytical techniques to further advance exploration at this site. This includes integration of geological (e.g., Quaternary fault mapping) and new geophysical datasets (e.g., gravity, magnetics, MT data, and five reprocessed seismic reflection profiles) to build a structural model and to identify specific favorable sites for potential geothermal upwellings. Two-meter temperature surveys were also conducted in the area (139 measurements). This part of north-central Nevada is characterized by systems of intersecting northerly and ENE-striking faults within the broader Humboldt structural zone, a poorly understood belt of ENE-striking faults and relatively high heat flow extending across northern Nevada. Kinematic analysis of exposed fault surfaces shows that ENE-striking faults have accommodated sinistral-normal slip, and normal slip characterizes N- to NNE-striking faults. Northeastern Reese River Valley lies within a broad left step between major ENE-striking fault zones on the northern flanks of the Argenta Rim and Shoshone Range and thus corresponds to a broad pull-apart in the ENE-striking sinistral-normal fault system. Notably, the nearby Beowawe geothermal system in Whirlwind Valley (with abundant sinter, hot springs, and a geothermal power plant) occupies a fault intersection in a relatively small left step in a major ENE-striking sinistral-normal fault and may serve as an analogue for a potential hidden system in northeastern Reese River Valley. Existing geological maps, high-resolution lidar, and seismic reflection data demonstrate that northeastern Reese River Valley is structurally complex with multiple intersections between the ENE- and N- to NNE-striking fault systems. Some of these fault intersections correspond to low resistivity anomalies, magnetic lows, and/or very subtle 2-m temperature anomalies, which may indicate hidden geothermal upwellings. Three-dimensional modeling and temperature-gradient drilling are planned to further evaluate these sites for geothermal activity.

Nevada↗
USGS2025

A 700-year rupture sequence of great eastern Aleutian earthquakes from tsunami modeling of stratigraphic records

Great Aleutian underthrusting earthquakes produced destructive tsunamis impacting Hawaiʻi in 1946 and 1957. Prior modeling of the 1957 tsunami deposit and runup records on eastern Aleutian and Hawaiian Islands jointly with tide-gauge observations across the Pacific Ocean constrained a rupture model with shallow slip up to 26 m along 600 km of the plate boundary. Here we implement this modeling approach to older deposits and show alternating deep and shallow megathrust slip up to 26, 32, and 22 m for great earthquakes along the same segment in the 18 th , 15 th , and 14 th centuries. All three modeled prehistoric Aleutian earthquakes produce tsunami inundation in Hawaiʻi with the most severe, 14 th century event having impacts exceeding the 1957 event. The along-dip variability of these four ruptures spanning seven centuries provides insights on earthquake cycles for engineering design and hazard assessment. The 15 th century and 1957 rupture models provide evidence for recurrence of tsunami earthquakes, which can produce disproportionately large tsunamis for a given moment magnitude due to reduced rigidity in the shallow megathrust. The 14 th and 18 th century events likely ruptured deeper regions that did not slip in 1957, suggesting potential for corresponding deeper failure in the next great eastern Aleutian earthquake.

Alaska, Hawaii↗
USGS2025

A comprehensive geologic framework of the National Crustal Model for seismic hazard studies in the conterminous United States

A three-dimensional (3D) geologic framework has been developed for the conterminous United States (U.S.) as part of the U.S. Geological Survey National Crustal Model to enhance seismic hazard modeling. The geologic framework is created from geologic maps and multiple subsurface geologic unit boundaries including the base of the Miocene, Cenozoic, Phanerozoic, and the Mohorovičić discontinuity. Modifications are made to surficial geologic maps to remove discontinuities across state and country borders. The subsurface distribution of rock type and age is extrapolated from the surface, seeded with subsurface geologic information, and constrained by a map of basement geology. The framework provides the basis for estimates of subsurface seismic velocity and density that is needed to improve estimates of earthquake ground shaking and seismic hazard. The present framework greatly expands and updates a previously published 3D geologic framework of the western part of the U.S. that was itself a first-of-its-kind digital 3D portrayal of the nation.

conterminous United States↗
USGS2026

Estimation of parental abundance using hierarchical Bayesian modeling with data augmentation

Pedigree-based estimation methods leverage the fact that each offspring in a cohort is genotypically “marked” by its parents and represent a recent and promising toolset for estimating population dynamics. This includes pedigree accumulation estimators that model the “accumulation” of inferred unique parents within a given cohort to estimate parental abundance. Unlike close-kin mark-recapture approaches, which rely on intercohort comparisons, pedigree accumulation modeling can be completed solely using intracohort samples. This is particularly advantageous for semelparous species, where intercohort pairs are impossible and adult life stages can be difficult to sample without affecting their likelihood of successfully reproducing. Previous work has evaluated a range of estimators for such datasets, concluding that the non-parametric Chao estimator provides the most accurate and precise estimates for feasible levels of sampling effort. We used simulated data to evaluate an alternative estimator based on hierarchical modeling and data augmentation in a Bayesian framework. Results indicate that estimates from the hierarchical Bayesian estimator had comparable accuracy and better precision than both the previously tested Chao1 estimator and the improved iChao formulation across a range of sample sizes and sex ratios. Furthermore, the Bayesian estimator was far more robust to simulated errors in pedigree reconstruction, especially the presence of false negatives. Hierarchical Bayesian pedigree accumulation models can also provide additional insight into underlying reproductive ecology through their use of an explicit observation process, allowing for the incorporation or estimation of species- and population-specific reproductive dynamics. More broadly, the parametric nature of these models offers opportunities to efficiently pool information among datasets as well as to propagate uncertainty within more complex models.

Ecology and Evolution↗
USGS2025

Site response and wave propagation effects in the eastern United States

Fourier amplitude spectra from regional earthquakes in the eastern United States are used in a parametric inversion for source, path, and site effects. Five earthquakes are selected for analysis during the installation of the United States National Seismic Network (US), Earthscope’s USArray Transportable Array (TA), and other temporary arrays to maximize station coverage. A global search algorithm is used to solve for site response from 0.1 to 15 Hz, corner frequency, geometrical spreading ( r - γ ), and frequency dependent anelastic attenuation in the form Q(f) = Q o f α . Tradeoff between moment and geometric spreading is handled by fixing the moment. The tradeoff between corner frequency and Q(f) is solved by selecting the value of corner frequency that minimizes an objective function defined over all stations. Values of site response and attenuation parameters show a strong spatial correlation with the physiographic provinces of the eastern United States. Site response for the Atlantic Coastal Plain is consistent with previous work using spectral ratios relative to a reference site, defined by strong resonance peaks correlated with the thickness of sediments. Site response for the other physiographic provinces is markedly different from the coastal plain, with a lack of distinct resonance peaks and a broad moderate high at frequences from 0.1 to 0.5 Hz consistent with the hard-rock geology of the regions. Like site response, Q(f) has a strong correlation with physiographic province, showing lower values on the coastal plain and higher values inland. Geometric spreading exponent, γ, decreases with increasing hypocenter distance from just above 1 at a few tens of kilometers to 0.9 at 500 km. The limited range in geometric spreading values is attributed to starting the Fourier transform window at the S ‐wave arrival for all distances and averaging over multiple wave types.

eastern United States↗
USGS2024

Distribution of ancient carbon in groundwater and soil gas from degradation of petroleum near the Red Hill Bulk Fuel Storage Facility, O‘ahu, Hawai‘i

The groundwater below the Red Hill Bulk Fuel Storage Facility (the facility) in Oʻahu, Hawaiʻi, contains fuel compounds from past spills. This study used carbon-14 analyses to distinguish fuel-derived carbon from background carbon, along with other biodegradation indicators, to address two goals: (1) determine the extent and migration direction of groundwater affected by residual fuel below the facility and (2) determine if residual fuel locations in the subsurface could be identified by analyzing soil gas at the surface above the facility. Groundwater from 19 wells was sampled between September 2022 and April 2023. Nonvolatile dissolved organic carbon (NVDOC) from a well presumed to be unaffected by past spills contained 38 percent ancient carbon indicating a natural source of ancient carbon in the subsurface. The NVDOC concentrations and ancient carbon percentages indicate fuel biodegradation products are likely present on the north and south of Red Hill with the greatest effects at well RHMW02 near the 2014 spill site. The NVDOC concentrations are almost three times higher than diesel range organic (DRO) concentrations in groundwater from the same sites. Major ion data indicate that iron reduction is an important biodegradation process. Soil probe samples and soil carbon traps were used to determine the carbon-14 content of soil carbon dioxide. Ancient carbon from fuel biodegradation was not detected at any soil probe or carbon trap site in contrast to a 2017 study which reported ancient carbon detections. A reanalysis of the 2017 results using a range of local values for background carbon-14 indicates that ancient carbon from fuel biodegradation was probably only detected in lower tunnel exhaust system samples and not in any soil carbon trap samples. Measurements of carbon dioxide efflux with a dynamic closed chamber were highly variable. The soil gas results indicate that soil gas measurements at land surface were not useful for detecting residual fuel at the facility.

Hawaii↗
USGS2026

Mid-Cretaceous porphyritic magmatism in the Beartooth Mountains of the northern Laramide foreland and its connection to Frontier Formation sediment dispersal in Bighorn basin, Montana-Wyoming (USA), and Cordilleran tectonism

Porphyritic intrusions in the southeast Beartooth Mountains and porphyritic cobbles from the Torchlight conglomerate of the Cenomanian Frontier Formation in Bighorn basin (Montana-Wyoming, USA) bear striking resemblance. We utilize geologic mapping, petrography, geochemistry, and geochronology to characterize rocks from both locations and test whether the Beartooth Mountains area was the provenance for Bighorn basin cobbles. Mapping on Line Creek Plateau in the Beartooth Mountains shows three units with porphyritic phenocrysts consisting dominantly of plagioclase, plagioclase + potassium feldspar, and plagioclase + quartz. Petrography shows embayed quartz, opacitic rims, and zoned plagioclase. Geochemistry characterizes the intrusions as trachydacite to trachyte and trace elements exhibit enrichment in light rare earth elements and slight depletion in heavy rare earth elements relative to primitive mantle. Zircon U-Pb dating shows substantial Archean inheritance and ca. 104−92 Ma dates; some samples display multiple mid-Cretaceous date populations. Bighorn basin porphyritic cobbles exhibit similar modal mineralogy, near-identical geochemistry, and Archean inheritance and mid-Cretaceous geochronology. We suggest the porphyritic cobbles were sourced from basement-hosted porphyritic intrusions now exposed in the Beartooth Mountains, given that geochemistry does not agree with coeval volcanic sources along the Frontier Formation depositional fairway (Idaho batholith) or sedimentary-hosted Beartooth Mountains porphyritic sills intruded into Cambrian shale. Stratigraphic evidence does not support kilometer-scale mid-Cretaceous Beartooth Mountains exhumation to expose basement-hosted porphyritic intrusions, and we therefore present a volcanic hypothesis to bring rocks to the surface. This mid-Cretaceous multimillion-year magmatism is unexplained by tectonic models for the north-central Laramide foreland, and thus we further explore the emplacement processes, spatiotemporal significance, and tectonic implications of these rocks.

Montana, Wyoming↗
USGS2025

Anatectic origin of Mississippian spodumene-bearing pegmatites in western Maine during orogenic plateau collapse

Spodumene pegmatites are an important lithium source, but the processes and tectonic settings in which they form are poorly understood. The Rumford pegmatite district surrounding Plumbago Mountain, western Maine, is host to numerous spodumene pegmatites, including the Plumbago North pegmatite (a world-class spodumene resource). Competing petrogenetic models for these spodumene pegmatites include (1) highly fractionated melts of the Mooselookmeguntic igneous complex and (2) anatexis. We tested these hypotheses by constraining the geologic, magmatic, metamorphic, and tectonic history of the Plumbago Mountain area with detailed geologic mapping and U-(Th)-Pb geochronology. The Silurian Rangeley Formation records initial isoclinal folding prior to, and contact-related metamorphism synchronous with, the intrusion of the 417 ± 4 Ma Plumbago Mountain pluton. Peak amphibolite facies metamorphism and crustal melting occurred during the ca. 410 to 400 Ma Acadian orogeny. Pulsed emplacement of the Mooselookmeguntic igneous complex occurred between ca. 389 and 356 Ma. Cassiterite U-Pb dates of spodumene pegmatites (333–327 Ma) are ≥23 m.y. younger than nearby granitic plutons, strongly arguing against the fractional crystallization model. Metamorphic monazite and xenotime (346–328 Ma) and 330 to 308 Ma 40 Ar/ 39 Ar hornblende dates indicate metamorphism coeval with spodumene pegmatite emplacement, supporting anatectic models. Reheating, anatexis, and spodumene pegmatite emplacement occurred during collapse of the 380 to 330 Ma Acadian orogenic plateau. Lithium enrichment may be linked to one or more stages of partial melting of metasedimentary and plutonic rocks during the formation, tenure, and collapse of the Acadian altiplano and emphasizes the role of anatexis in producing spodumene pegmatites of economic significance.

Maine↗
USGS2025

Estimating the probability of export restrictions to inform mineral criticality

To assess risks associated with advanced technologies’ supply chain disruptions, governmental agencies and others have developed mineral “criticality” assessments, with criticality described using the economic impact and probability of supply chain disruptions. Previous work developed subjective supply risk indicators to approximate this probability, typically combining several factors such as supply diversity and trading partners’ political stability, where indicator weightings can substantially impact results. This work explicitly quantifies export barrier probability using an ensemble of machine learning classifiers, with probability estimates informed by exogenous variables, including prior barrier implementation and global export dominance. Major differences in high-probability countries and commodities are observed across models, but the ensemble method highlights Indonesia, China, Tanzania, and the United States as particularly high risk. The Supplementary Data File provides export barrier probability estimates for each analyzed country-commodity pair, enabling a direct, quantitative, objective contribution to assessing mineral criticality, enhancing risk identification and prioritization for policymakers.

Resources, Conservation, and Recycling↗
USGS2025

Hydrogeologic investigation, framework, and conceptual flow model of the Antlers aquifer, southeastern Oklahoma, 1980–2022

The 1973 Oklahoma Groundwater Law (Oklahoma Statute §82–1020.5) requires that the Oklahoma Water Resources Board conduct hydrologic investigations of the State’s groundwater basins to support a determination of the maximum annual yield for each groundwater basin. Every 20 years, the Oklahoma Water Resources Board is required to update the hydrologic investigation on which the maximum annual yield determinations were based. The maximum annual yield allocated per acre of land is used to set the equal-proportionate share pumping rate. The maximum annual yield of 5,913,600 acre-feet per year and equal-proportionate-share of 2.1 acre-feet per acre per year currently (2025) in place for the Antlers aquifer were issued by the Oklahoma Water Resources Board on February 14, 1995. Because more than 20 years have elapsed since the 1995 final order for the Antlers aquifer was issued, the U.S. Geological Survey, in cooperation with the Oklahoma Water Resources Board, completed an in-depth hydrologic study that included a hydrogeologic framework and conceptual groundwater-flow model for the 1980–2022 study period. The results of an analysis of land use, long-term climate patterns, streamflow and base-flow patterns, historical groundwater use, as well as groundwater-level fluctuations across the Antlers aquifer are described. In addition, groundwater quality was analyzed for total dissolved solids concentrations and major ions for the Antlers aquifer. An updated hydrogeologic framework was developed that included refining the aquifer boundary in Oklahoma, the creation of new potentiometric surface and saturated thickness of fresh groundwater maps, one multiple-well aquifer test, slug tests, and an analysis of lithologic logs across the aquifer. A conceptual groundwater flow model and water budget were developed by incorporating estimates of recharge from precipitation, saturated-zone evapotranspiration, streambed seepage, lateral groundwater flows, vertical leakage, and withdrawals from groundwater wells.

Oklahoma, Texas↗