USGS Science⌕ Search

SEARCH · USGS Science

Results for “Data Science in 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.

1,642 records · Page 31Linked to original sources
USGS1995

The distribution of bromide in water in the Floridan aquifer system, Duval County, northeastern Florida

Although Duval County, Florida, has ample ground-water resources for public supply, the potential exists for a problem with excessive disinfectant by-products. These disinfectant by-products result from the treatment of raw water containing low concentrations of bromide and naturally occurring organic compounds. Because of this potential problem, the relation of bromide concentrations to aquifer tapped, well location and depth, and chemical characteristics of water in the Floridan aquifer system underlying Duval County were studied to determine if these relations could be applied to delineate water with low-bromide concentrations for future supplies. In 1992, water samples from 106 wells that tap the Floridan aquifer system were analyzed for bromide and major dissolved constituents. A comparison of bromide concentrations from the 1992 sampling with data from earlier studies (1979-80) indicates that higher bromide concentrations were detected during the earlier studies. The difference between the old and new data is probably because of a change in analytical methodology in the analysis of samples. Bromide concentrations exceeded the detection limit (0.10 milligrams per liter) in water from 28 of the 106 wells (26 percent) sampled in 1992. The maximum concentration was 0.56 milligrams per liter. There were no relations between bromide and major dissolved constituents, well depth, or aquifer tapped that would be useful for determining bromide concentrations. Areal patterns of bromide occurrence are not clearly defined, but areas with relatively high bromide concentrations tend to be located in a triangular area near the community of Sunbeam, Florida, and along the St. Johns River throughout Duval County.

Florida↗
USGS2024

Surficial geology and Quaternary fault map of the Las Vegas Valley, Clark County, Nevada

The Surficial Geology and Quaternary Fault Map of the Las Vegas Valley, Clark County, Nevada is a 1:50,000-scale compilation of published 1:24,000-scale geologic maps integrated with new field and desktop mapping. This geologic map compilation and GIS database are part of a broader study on the Quaternary faults in the Las Vegas Valley, the results of which are summarized in this report. The map compilation utilized a variety of digital base maps to evaluate the characteristics of surficial deposits and Quaternary fault scarps including lidar data in undeveloped parts of the valley and orthorectified historical aerial photos in the urbanized center of the valley. The map distinguishes twenty surficial map units including alluvial-fan and ephemeral-wash deposits, groundwater discharge deposits of the Las Vegas Formation, and alluvium deposits of Las Vegas Wash. Quaternary faults include the Quaternary Las Vegas Valley fault system in the center of the valley and the Frenchman Mountain fault system along the eastern margin. All Quaternary fault traces in the valley were evaluated, and the previously published fault mapping was modified as needed. The new mapping and uniform correlation of surficial deposits displaced by Quaternary faults yields a consistent valley-wide characterization of the recency of fault activity. The map and GIS database include a Quaternary geochronology compilation of 119 previously published surficial deposit ages of various types, 35 new luminescence ages, and 3 new radiocarbon ages. The new age data further bracket the ages of the Las Vegas basin Quaternary stratigraphy and provide new constraints on the timing of Quaternary fault activity.

Nevada↗
USGS2025

Estimated hydrogeologic, spatial, and temporal distribution of self-supplied domestic groundwater withdrawals for aquifers of the Virginia Coastal Plain

Water use from private-domestic wells accounts for nearly 40 percent of total groundwater withdrawals in the Virginia Coastal Plain Physiographic Province (henceforth called the Virginia Coastal Plain). However, because self-supplied domestic water use generally falls below the Virginia Department of Environmental Quality (VDEQ) reporting and management threshold of 300,000 gallons per month, quantifying these withdrawals is challenging. This report builds upon the foundation of previous U.S. Geological Survey investigations by providing revised techniques to improve estimates of the aquifer source, spatial distribution, and monthly magnitude of these groundwater withdrawals. The aquifer sources of private-domestic wells in the Virginia Coastal Plain were estimated by cross-referencing 8,264 well records from the VDEQ and the Virginia Department of Health to a digital model of the Virginia Coastal Plain hydrogeologic framework. This analysis highlights the regional importance of the Yorktown-Eastover, Potomac, and surficial aquifers. Collectively, these three aquifers account for 80 percent of self-supplied domestic groundwater withdrawals. The population using self-supplied domestic water was estimated using census blocks, well-use ratios, building footprints, and land-use and land-cover data to produce a high-resolution, disaggregated, raster-based dataset. This approach improves upon previous models at the census-block or road-network scale by reducing the low-density spread of the self-supplied domestic population across undeveloped areas and concentrating the population and its corresponding water use in the areas where it is most likely to occur. Results show that an estimated 475,332 people comprise the 2020 self-supplied domestic population of the Virginia Coastal Plain, an increase of 5.7 percent since 2010, and the greatest concentrations of self-supplied domestic population surround large cities. Estimates could be further refined with the addition of current and complete spatial data on public water-system service areas. The quantity of water used by the self-supplied domestic population was estimated by modifying published state per-capita water-use coefficients with the corresponding monthly variability assessed from Virginia Coastal Plain public water-system withdrawal data. This analysis estimates an average increase of 12 percent from June through August and an average decrease of 8 percent from December through March from the baseline annual average of 80 gallons per day per capita, which generally matches similar studies in the eastern United States. The application of these revised methodologies for the estimation of private-domestic wells and the self-supplied domestic population improves understanding of domestic groundwater use in the Virginia Coastal Plain across hydrogeologic, spatial, and temporal scales. These revisions help better inform water-resource managers and decision makers and support higher resolution groundwater modeling. Furthermore, these methods are transferrable to other areas where self-supplied domestic water withdrawals are important to the overall water budget.

Virginia↗
USGS2026

Key minerals in energy transmission

The electrical grid relies on minerals with unique properties to transform, transmit, and convert electricity. Transmission needs are growing—to support the expansion of data centers and link diverse energy resources into the grid. This expansion will require greater quantities of minerals. Although some minerals are produced in the United States, the United States relies on imports for others. This infographic shows the key minerals that support energy transmission, and how reliant the United States is on imports for these minerals.

General Information Product↗
USGS2025

Using the D-Claw software package to model lahars in the Middle Fork Nooksack River drainage and beyond, Mount Baker, Washington

Lahars, or volcanic mudflows, are the most hazardous eruption-related phenomena that will affect communities living along rivers that originate on Mount Baker. In the past 15,000 years, the largest lahars from Mount Baker have affected the Middle Fork Nooksack River drainage and beyond. Here we use the physics-based D-Claw software package to model nine lahar scenarios that are initiated as water-saturated landslides between Sherman Crater and the Roman Wall on the Mount Baker edifice and flow down the Middle Fork Nooksack River. The scenarios range in volume from 1 to 260 million cubic meters and have an initial hydraulic permeability from 10 −12 to 10 −10 meters squared. Model output includes data such as flow depth, velocity, runout distance, area inundated, arrival time, and sediment concentration as well as information that allows scientists to calculate other important hydrologic characteristics such as lahar discharge. These data are important to officials who have the responsibility to plan for, or take mitigation measures against, future Mount Baker lahars. To check the validity of the D-Claw results, we compare the scenarios to known geologic information. We also compare D-Claw results with empirical models that have been used in the past to determine potential inundation areas, runout distances, and arrival times. These comparisons highlight similarities and differences between empirical and physics-based models. We also present D-Claw scenario-based animations to help scientists, officials, and lay people alike to visualize how future lahars could affect communities.

Washington↗
USGS2026

Microclimate mediates the strength and direction of avian biotic interactions

Theory predicts that that the strength and direction of species interactions can shift from being competitive in benign environments toward being facilitative in stressful environments. However, the environmental context dependency of species interactions has rarely been tested in animal communities. We capitalized on a 15-year, landscape-scale dataset, collected annually in a relatively stable old-growth forest environment to test the long-held hypothesis that the strength and direction of species interactions might be mediated by climatic conditions. It is generally accepted that competitive and facilitative interactions drive the distributions of many species. Using multi-species dynamic occupancy models applied to long-term data, we tested whether annual settlement by bird species could affect either the persistence or settlement by other phylogenetically related species, and whether these interactions are mediated by microclimate. We found that species interactions were influenced by microclimate for some, but not all avian species pairs. Related species pairs more often showed settlement dynamics that were indicative of attraction rather than repulsion. In some cases, competitive interactions at warmer microclimates flipped to become facilitative in colder ones. However, the reverse was also true: facilitative interactions amplified for some species under warm conditions. Furthermore, species pairs that were closely related were more likely to exhibit competitive relationships along at least part of the microclimatic gradient. Our results highlight the importance of using long-term data to incorporate competitive and facilitative interactions into species distribution models and support the notion that the strength and direction of species interactions can be dependent on microclimatic environmental conditions.

Oregon↗
USGS1998

Water-quality assessment of the Ozark Plateaus study unit, Arkansas, Kansas, Missouri, and Oklahoma — Nutrients, bacteria, organic carbon, and suspended sediment in surface water, 1993-95

Nutrient, bacteria, organic carbon, and suspended- sediment samples were collected from 1993-95 at 43 surface-water-quality sampling sites within the Ozark Plateaus National Water- Quality Assessment Program study unit. Most surface-water-quality sites have small or medium drainage basins, near-homogenous land uses (primarily agricultural or forest), and are located predominantly in the Springfield and Salem Plateaus. The water-quality data were analyzed using selected descriptive and statistical methods to determine factors affecting occurrence in streams in the study unit. Nitrogen and phosphorus fertilizer use increased in the Ozark Plateaus study unit for the period 1965-85, but the application rates are well below the national median. Fertilizer use differed substantially among the major river basins and physiographic areas in the study unit. Livestock and poultry waste is a major source of nutrient loading in parts of the study unit. The quantity of nitrogen and phosphorus from livestock and poultry wastes differed substantially among the river basins of the study unit's sampling network. Eighty six municipal sewage-treatment plants in the study unit have effluents of 0.5 million gallons per day or more (for the years 1985-91). Statistically significant differences existed in surface-water quality that can be attributed to land use, physiography, and drainage basin size. Dissolved nitrite plus nitrate, total phosphorus, fecal coliform bacteria, and dissolved organic carbon concentrations generally were larger at sites associated with agricultural basins than at sites associated with forested basins. A large difference in dissolved nitrite plus nitrate concentrations occurred between streams draining basins with agricultural land use in the Springfield and Salem Plateaus. Streams draining both small and medium agricultural basins in the Springfield Plateau had much larger concentrations than their counterparts in the Salem Plateau. Drainage basin size was not a significant factor in affecting total phosphorus, fecal coliform bacteria, or dissolved organic carbon concentrations. Suspended-sediment concentrations generally were small and indicative of the clear water in streams in the Ozark Plateaus. A comparison of the dissolved nitrite plus nitrate, total phosphorus, and fecal coliform data collected at the fixed and synoptic sites indicates that generally the data for streams draining basins of similar physiography, land-use setting, and drainage basin size group together. Many of the variations are most likely the result of differences in percent agricultural land use between the sites being compared or are discharge related. The relation of dissolved nitrite plus nitrate, total phosphorus, and fecal coliform concentration to percent agricultural land use has a strong positive 2 Water-Quality Assessment-Nutrients, Bacteria, Organic Carbon, and Suspended Sediment in Surface Water, 1993-95 correlation, with percent agricultural land use accounting for between 42 and 60 percent of the variation in the observed concentrations.

Arkansas, Kansas, Missouri, Oklahoma↗
USGS2023

Eruptive history of northern Harrat Rahat—Volume, timing, and composition of volcanism over the past 1.2 million years

Harrat Rahat, one of several large, basalt-dominated volcanic fields in the western part of the Kingdom of Saudi Arabia, is a prime example of continental, intraplate volcanism. Excellent exposure makes this an outstanding site to investigate changing volcanic flux and composition through time. We present 93 40 Ar/ 39 Ar ages and 6 36 Cl surface-exposure ages for volcanic deposits throughout northern Harrat Rahat that, integrated with a new geologic map, define 12 eruptive stages. Exposed volcanic deposits in the study area erupted less than 1.2 million years ago (Ma), and 214 of 234 identified eruptions occurred less than 570 thousand years ago (ka). Two eruptions were in the Holocene, including a historically described basaltic eruption in 1256 C.E. and a trachyte eruption newly recognized as Holocene (4.2±5.2 ka). An estimated approximately 82 cubic kilometers (km 3 ; dense rock equivalent) of volcanic products can be documented as having erupted since 1.2 Ma, though this is a lower limit because of concealment of deposits older than 570 ka. Over the last 570 thousand years (k.y.), the average eruption rate was 0.14 cubic kilometers per thousand years (km 3 /k.y.), but volcanism was episodic with periods alternating between low (0.04–0.06 km 3 /k.y.) and high (0.1–0.3 km 3 /k.y.) effusion rates. Before 180 ka, eruptions vented from the volcanic field’s dominant eastern vent axis and from a subsidiary, diffuse, western vent axis. After 180 ka, volcanism focused along the eastern vent axis, and the composition of volcanism varied systematically along its length from basalt dominated in the north to trachyte dominated in the south. We hypothesize that these compositional variations younger than 180 k.y. reflect the growth of a mafic intrusive complex beneath the southern part of the vent axis, which led to the development of evolved magmas. Lastly, these new age data allow for a reassessment of the volcanic recurrence interval at northern Harrat Rahat. Based on available data, volcanism in northern Harrat Rahat over the last 180 k.y. is poorly described using a Poisson distribution with a single recurrence interval. Instead, data for northern Harrat Rahat are better described using a mixed exponential distribution that is applicable for volcanic systems characterized by two different eruptive states, where one state with a longer recurrence interval corresponding to periods of low eruption frequency and one state with a shorter recurrence interval corresponding to periods of high eruption frequency. The preferred model for northern Harrat Rahat over the last 180 k.y. uses a long recurrence interval of 4.0 k.y. and a short recurrence interval of 0.22 k.y.

Professional Paper↗
USGS2026

Magnitude conversion relations create substantial differences in seismic hazard models

Earthquake catalogs are essential data inputs for seismic hazard modeling. Because earthquake magnitudes are reported in a variety of types (e.g., local magnitudes and moment magnitudes), magnitude conversion relationships must be used to convert the different magnitude types present in a catalog to a uniform magnitude type to avoid biases in the hazard computation. However, these conversion relationships are often uncertain and have been shown to sometimes perform poorly. Here, we investigate the sensitivity of the gridded seismicity component of the National Seismic Hazard Model (NSHM) to the catalog conversion equations in the Eastern United States. In the 2023 NSHM, magnitudes of various types were converted to moment magnitudes using equations developed by the Central and Eastern United States Seismic Source Characterization for Nuclear Facilities (CEUS‐SSCn), based on least‐squares (LS) regressions made using data from a catalog containing events up through 2008. We recompute these equations using events in the Advanced National Seismic System Comprehensive Earthquake Catalog with multiple magnitudes from 2000 to 2023. Although we prefer the use of orthogonal regressions for our datasets, LS regressions produce broadly similar results, with both approaches exhibiting large deviations from the CEUS‐SSCn conversions, especially at smaller magnitudes. We compare the spatial distribution of annual rates using three different models: (1) the 2023 NSHM conversions, (2) our updated conversions, and (3) no conversions. We find that the choice of conversions leads to substantial differences in the rate forecasts, which can greatly impact the seismic hazard model, particularly in regions with low‐seismicity rates such as the Eastern United States, where the hazard is dominated by gridded seismicity rather than a fault model.

Seismological Research Letters↗
USGS2026

Channel morphology and large wood control postfire debris-flow erosion and deposition

Runoff-generated debris flows are a known response to wildfire, and accurately predicting the volume of these debris flows is important for estimating the magnitude of downstream hazards. Prior data collection efforts have focused on debris-flow volume measurements at catchment outlets, but few studies have considered how erosion and deposition modulate the volume of sediment arriving at catchment outlets. This study takes advantage of a high-resolution dataset to document the factors that control the total debris-flow volume reaching the catchment outlet during a fatal postfire debris flow. Using pre- and post-event airborne lidar, satellite imagery and field mapping, we found that a postfire debris flow in the Black Hollow catchment in northern Colorado eroded 136,000 ± 30,000 m 3 and redeposited 27,000 ± 7,500 m 3 in the main channel. Most of the in-channel deposition (52% by volume) occurred where a confined channel reach transitioned to an unconfined channel reach downstream, allowing the flow to widen and deposit material. Wood jams played multiple roles in the debris-flow dynamics, both nucleating deposition (25% of the deposit volume was stored behind wood jams) and exacerbating erosion (50% of the total erosion occurred downstream from a wood dam break). The remaining deposition occurred due to spatial changes in channel slope as well as deposition observed at newly formed channel bars. Using these data in this study, we identified topographic and vegetation metrics that can be used (pre-event) to anticipate where deposition may occur in channels prior to a debris flow.

Colorado↗
USGS2025

Hydraulic conductivity and transmissivity estimates from slug tests in wells within the Mississippi Alluvial Plain, Arkansas and Mississippi, 2020

During the spring and summer of 2020, the U.S. Geological Survey conducted single-well slug tests on selected observation wells within the Mississippi Alluvial Plain in Arkansas and Mississippi to estimate hydraulic conductivity and transmissivity values for the Mississippi River Valley alluvial and middle Claiborne aquifers. Well and aquifer data were collected from field measurements, well-construction reports, and published aquifer-thickness information. A total of 324 slug-in and slug-out tests were conducted on 48 wells by using mechanical slugs to displace the water column and submersible pressure transducers to record changes in water levels in the wells. Hydraulic conductivity of the aquifers in which the wells are screened was estimated by curve fitting the water-level-change data using aquifer test analysis software. Estimates of aquifer transmissivity were made by multiplying the estimated hydraulic conductivity value by the aquifer thickness at well locations. Mean hydraulic conductivity estimates for 44 observation wells screened in the Mississippi River Valley alluvial aquifer range from 3 to 401 feet per day, and mean transmissivity estimates range from 285 to 80,559 feet squared per day. Mean hydraulic conductivity estimates for four observation wells screened in units of the middle Claiborne aquifer range from 0.14 to 183 feet per day, and mean transmissivity estimates range from 55 to 67,913 feet squared per day. The results from these tests can be used to improve the understanding of water availability and groundwater migration, to refine groundwater models, and to ultimately provide stakeholders and decisionmakers better information for management of the groundwater resources within the Mississippi Alluvial Plain.

Arkansas, Mississippi↗
USGS2026

Top-down targeted network analysis of critical mineral commodities applied to international geochemistry database

The global demand for critical mineral commodities is rapidly increasing, making domestic production an important factor in supporting the economy and national security. Large scale, publicly available geochemical databases allow for the application of data informatics methods to interrogate critical mineral commodities data for correlations in deposit formation and distribution, particularly for identifying enrichment of multiple critical mineral commodities at the same deposit. In this study, we applied network analysis to the Critical Minerals Mapping Initiative (CMMI) ore geochemistry (Critical Minerals in Ores, CMiO) database to identify the high concentration (defined as 10× bulk crustal abundance) co-occurrence of different critical mineral commodities across a mineral system hierarchy from deposit environments to individual deposits. Identifying patterns or unique outliers in enrichment in network communities will allow for the location of secondary critical mineral commodity resources from under-utilized deposits. We find trends in the enrichment of critical mineral commodities in network-communities between the elements praseodymium (Pr), neodymium (Nd), terbium (Tb), and dysprosium (Dy) across multiple CMiO database deposit environments and groups down to specific deposit types and sites. A separate trend in network community deposition is observed as well between iridium (Ir) and platinum (Pt) in deposit environments, groups, types, and sites. Network analysis focused on critical minerals in magmatic-hydrothermal deposits identified multiple deposit sites from different deposit types within the CMiO database with concentrations of Dy, Nd, Tb, Pr, Ir, and Pt that are at least ten times greater than the crustal average. This approach can be applied to any target element(s) or deposit(s) of interest, allowing broad investigation of co-enriched critical mineral commodities.

Journal of Geochemical Exploration↗
USGS2024

Peri-Gondwanan sediment in the Arkoma Basin derived from the north: The detrital zircon record of a uniquely concentrated non-Laurentian source signal in the late Paleozoic

During the assembly of Pangea, peri-Gondwanan terranes collided with the eastern and southern margins of Laurentia and brought with them unique detrital zircon U-Pb signatures. Discriminating between individual peri-Gondwanan terranes in the detrital record is difficult due to their similar geologic histories. However, characterization of this provenance is critical for understanding late Paleozoic sediment routing during development of Pangea. Along southeastern Laurentia, in the Arkoma Basin (present-day Arkansas and eastern Oklahoma, southeastern United States), we identified Middle Pennsylvanian (Desmoinesian) strata that exhibit a concentrated peri-Gondwanan detrital zircon signature (e.g., ca. 800–550 Ma). Although several southern peri-Gondwanan terranes (e.g., Maya, Suwannee) are closer to the Arkoma Basin, geologic data, such as predominantly north-to-south paleocurrents and proximal-to-distal facies relationships in these Desmoinesian strata, support a northern source (e.g., Ganderia, Avalonia, Meguma). Further evidence of a northern source comes from detrital zircon source mapping, which reveals the persistence of this peri-Gondwanan signal in depocenters to the north of the basin after the signal had diminished in the Arkoma Basin. To this end, bottom-up detrital zircon source modeling, source mapping, regional stratigraphy, paleocurrent data, and sandstone petrography allow us to reconstruct the evolution of this Middle Pennsylvanian (Desmoinesian) sediment pathway in the context of intraplate and plate-margin tectonic activity. This reconstruction documents processes affecting Earth’s surface (e.g., tectonics, climate) during the assembly of Pangea and describes in detail part of a dynamic continental-scale drainage system.

Geosphere↗
USGS2025

Quantifying groundwater response and uncertainty in beaver-influenced mountainous floodplains using machine learning-based model calibration

Beavers ( Castor canadensis ) alter river corridor hydrology by creating ponds and inundating floodplains, and thereby improving surface water storage. However, the impact of inundation on groundwater, particularly in mountainous alluvial floodplains with permeable gravel/cobble layers overlain by a soil layer, remains uncertain. Numerical modeling across various floodplain structures considers topographic and sediment complexity and multidirectional flow, linking inundation to groundwater response. This study develops a model-data integration workflow to address uncertainty in groundwater response to beaver-induced inundations in a mountainous alluvial floodplain in the Upper Colorado River Basin. Uncertain factors include seasonal hydrologic dynamics, hydraulic conductivities, floodplain structures, and meteorological forcings. We employed an ensemble of groundwater models, based on geophysical and hydrologic data, with machine learning-based calibration using a neural density estimator. This allowed us to quantify the vertical flux from the soil layer to the permeable gravel bed, the down-valley underflow within the gravel bed, and their ratios. Results show a significant increase in the vertical flux relative to down-valley underflow, from 2% during dry pond periods to 20% during wet periods, serving as an analogy for conditions without and with beaver ponds. The study highlights the influence of floodplain structure on groundwater storage, water balance, and water quality impacted by beaver ponds. A thick gravel bed layer, with a large down-valley underflow, minimizes the effect of beaver-induced inundation on water quality. We emphasize the need for field-scale measurements of floodplain structure and improved characterization of evapotranspiration changes to reduce uncertainty in groundwater response.

Colorado↗
USGS2026

Integrating detrital magnetite geochemistry and (U-Th)/He chronometry as a sediment provenance tool in geologic and metallogenic terranes

Magnetite is ubiquitous in porphyry Cu systems and in sediment sourcing both barren and mineralized regions, with potential as an indicator mineral in concealed and coarsely-mapped terranes. We develop and test a workflow for integrated geochemistry and (U–Th)/He (He) dating for inferring detrital magnetite (DMt) provenance in these settings. The ca. 70 Ma Taurus porphyry Cu–Mo(–Au) district in eastern interior Alaska serves as a test case. DMt from streams draining porphyry-related mineralization was characterized by geochemistry and mineral inclusion and microstructure observations, complemented by similar data for potential porphyry and host rock sources. Principal component analysis and clustering of DMt geochemical data resolve multiple populations in our samples geochemically and texturally compatible with derivation from metamorphic, porphyry-related hydrothermal, and igneous sources. Hydrothermal magnetite comprises ∼16–50% of DMt nearest porphyry mineralization but diminishes to ∼4% ∼15 km downstream. Subsampled grains within populations yield ∼160–110 Ma, ∼70 Ma, ∼55 Ma, and ∼20 Ma magnetite He date modes. Combined with provenance, He dates capture Early Cretaceous regional exhumation of metamorphic host rock and Late Cretaceous porphyry Cu mineralization. DMt grains showing partial hematite replacement yield ca. 55–20 Ma dates regardless of source, overlapping regional warm/wet climatic intervals. We interpret Cenozoic dates to reflect exhumation to near-surface oxidizing conditions and(or) supergene weathering. Magnetite is thus a promising target phase for (1) tracking the spatiotemporal distribution of porphyry systems, and (2) linking the formation and exhumation of these systems to a regional geologic history, both in Alaska and globally.

Alaska, Yukon↗
USGS2026

Ductile and brittle Rio Grande Rift deformation in Oligocene granite records a two-stage rift history in southern Colorado

The timing and nature of early deformation in the Rio Grande Rift remains poorly constrained. We present evidence for the earliest structural signature of rift extension in the Sangre de Cristo Range, southern Colorado, based on new geologic mapping, structural analysis, rock magnetic data, and thermochronology. These analyses focus on the ~30.0 Ma granite of Chokecherry Canyon, which hosts discrete low-angle mylonitic shear zones and a distributed, gently SW-dipping protomylonitic fabric. Incremental stretching axes, stretching lineations, and Kmax magnetic lineations plunge gently WSW. Quartz microstructures and crystallographic orientations indicate dominantly coaxial strain in the protomylonite and general shear in the discrete shear zones. Quartz c-axis opening-angle thermometry suggests deformation at ~420–540°C. Thermal modeling of ⁴⁰Ar/³⁹Ar K-feldspar data indicates rapid post magmatic cooling below the brittle–plastic transition, supporting shear-zone formation immediately after emplacement. Slow cooling from ~20–13 Ma was followed by renewed rapid cooling at ~13 Ma, interpreted as the onset of extensional exhumation along the Sangre de Cristo Fault System. These results show that extension in the northern Rio Grande Rift was active by ~30 Ma, earlier than previously recognized. We propose a two-stage model for northern Rio Grande Rift evolution: Stage I (30–23 Ma) records ENE–WSW extension localized in low-angle mylonitic shear zones associated with mid-crustal intrusions; Stage II (≤18 Ma) reflects brittle high-angle normal faulting, focused exhumation, and rift narrowing. Stage I magmatism and deformation along the western range front likely established crustal weaknesses that guided later fault development.

Colorado↗
USGS2025

An early Holocene wet period in the southwestern United States

Multiple generations of spring-fed streams traversed ∼800 km 2 of the Las Vegas Valley in southern Nevada between ca. 10.9 ka and 8.5 ka, depositing an extensive tufa network. The scale of this network and diversity of tufa morphologies is novel in North America and offers an opportunity to obtain quantitative paleoclimate data for the region during the early Holocene. We determined isotopic compositions and estimated past temperatures using clumped isotope data from early Holocene tufa on the valley floor (698 m) as well as tufa forming today at higher elevation in the nearby Spring Mountains at Cold Creek Spring (1856 m). Modern and fossil tufa yielded comparably low δ 18 O values, implying that source waters for both were derived from high-elevation winter precipitation. Clumped isotope temperatures of modern tufa average 15.8 ± 2.5 °C, aligning with mean summer temperatures of the emergent spring water, and indicate equilibrium conditions of tufa formation. The early Holocene tufa yielded similar clumped isotope temperatures, averaging 15.2 ± 3.9 °C, meaning it precipitated at temperatures that occur at much higher elevations today. The Las Vegas tufa record, combined with nearby and temporally correlative paleospring and lacustrine records, suggest that cool/wet conditions prevailed throughout the Mojave Desert during the early Holocene. These records also demonstrate that spring ecosystems responded to millennial-scale hydroclimate variations that supersede climate change driven solely by insolation. The previously unrecognized pattern of ecosystem response to hydroclimate documented here may assist in understanding climate drivers for the early Holocene and provide critical information for the fate of groundwater-dependent ecosystems in the southwestern United States.

Nevada↗
USGS2025

Factors affecting the distribution of water-bearing fractures in the bedrock aquifers of West Virginia

Bedrock aquifers cover 23,601 square miles within the State of West Virginia and comprise 97.4 percent of the surficial area within the State; the remaining 2.6 percent (621 square miles) consists of alluvial sand-and-gravel and glacial outwash aquifers bordering the State’s major rivers. While West Virginia’s alluvial aquifers have been studied extensively, bedrock aquifers have only been characterized for studies completed in a few areas in Jefferson, McDowell, and Monroe Counties. Bedrock aquifers are water supplies for public supply, agriculture, industry, and residential homeowner use. In this study, the U.S. Geological Survey, in cooperation with the West Virginia Department of Environmental Protection Division of Water and Waste Management, provides a statewide assessment of the occurrence and distribution of fractures within bedrock aquifers of the State and the various topographic, physiographic, and lithologic influences controlling the occurrence and distribution of bedrock fractures. The results of this study provide an increased understanding of the distribution of fractures in bedrock aquifers in West Virginia and help to verify trends that have been suspected for many years but were never well documented or verified by data. The results confirmed that the density of fractures and those that were determined to be water bearing decrease significantly with depth. A statistically significant difference in the density of fractures was observed at a depth of 215 feet for wells in the Appalachian Plateaus Physiographic Province’s and in the Valley and Ridge Physiographic Province’s aquifers; a higher density of fractures and water-bearing fractures were above a depth of 215 feet than below that depth. This is an important consideration when drilling wells for residential, commercial, industrial, or agricultural water supply. Abandoned underground coal mines are commonly believed to form large pools of water in the interconnected mine entries in abandoned room and pillar coal mines. Such pools of water can and do exist in abandoned underground coal mines, but many mines lack open entries and are held up by overburden strata and pillars that can collapse and form aquifers comprised of vast interconnected rubble zones (gob), especially in older mines. Data assessed for this study showed that shale-corrected values of effective porosity for limestone aquifers in West Virginia had a median value of 2 percent and an average value of 4 percent and generally are mineralized with low effective porosity. Argillaceous or sandy limestone has a median shale-corrected porosity of 4 percent and an average shale-corrected porosity of 5 percent. The median and average shale-corrected porosity of sandstone aquifers was estimated to be 14 percent, but the median shale-corrected porosity for argillaceous or calcareous sandstone was 5 percent and the average shale-corrected porosity for argillaceous or calcareous sandstone was 6 percent. Even though shale has a relatively high total sonic porosity compared to other lithologies, shale and siltstone had relatively low shale-corrected porosity, ranging from 0 to 2 percent. Well yields were previously documented to be highest in valley settings, lowest on hilltops, and intermediate on hillsides. Transmissivity data provided by this study confirm this general pattern within the Appalachian Plateaus Province; however, the Valley and Ridge Province does not follow this pattern. While still lowest on hilltop settings, the highest well yields were in hillside settings. The trend for the Valley and Ridge Province was likely skewed because of 9 high-yield wells specifically targeting deeper thin limestone units, such as the Tonoloway and Helderberg Limestones, at depths with transmissivity in excess of 2,000 feet squared per day in Mineral County, West Virginia, or targeting karst aquifers in Berkeley, Jefferson, or Greenbrier Counties, West Virginia. Finally, water-bearing fractures have been hypothesized to comprise a small number of all fractures within a typical bedrock well in West Virginia. Data collected for this study support this theory. A total of 3,403 fractures were identified during this study; 3,151 (92.6 percent) of those fractures are low-transmissive, and only 252 (7.4 percent) fractures are water-bearing. Even though a well may contain many fractures, less than 8 percent are considered water-bearing fractures.

West Virginia↗