USGS Science⌕ Search

SEARCH · USGS Science

Results for “Environmental Science and Technology - Water”

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

Ground-motion aleatory-variability models for Puerto Rico and the U.S. Virgin Islands

I develop independent logic trees for aleatory variability for crustal and subduction-zone (interface and intraslab) earthquakes for seismic hazards analyses in Puerto Rico and the U.S. Virgin Islands (PRVI) from existing suites of ground-motion models (GMMs) and from ground-motion datasets, including a regional PRVI dataset. The aleatory variability models are parameterized using a partially nonergodic partitioning of standard deviation that consists of independently developed between-event ( ), site-to-site ( ), and event-corrected single-station ( ) standard deviation components. The effects of nonlinear site response on aleatory variability are incorporated through additional terms that modify the standard deviation components. Because one goal of this work is to develop independent logic trees for aleatory variability that synthesize the aleatory variability models from GMMs, I make use of the functional forms of the input GMMs. The PRVI dataset contains a limited number of stations with high-quality site metadata and does not contain records from earthquakes with magnitudes greater than 6.1, so I choose not to develop the aleatory variability models from the regional dataset alone. Instead, the standard deviation components from regional ground-motion data are evaluated against the components derived from GMMs and from available global datasets, and regionalized standard deviation components are incorporated where there is evidence that regional effects exhibit substantial differences. The resulting logic trees for aleatory variability consist of models of and that are consistent with semiempirical GMMs for active crustal and subduction-zone regimes, and two alternative models of , including one model that exhibits site-to-site variability informed by PRVI data, with values that exceed global models. The aleatory variability models may be considered in future hazards assessments in PRVI to simplify the hazard calculations, to incorporate regional ground-motion variability effects, and to enable direct logic-tree weighs of aleatory variability.

Puerto Rico, U.S. Virgin Islands↗
USGS2024

Federal lands greenhouse gas emissions and sequestration in the United States: Estimates for 2005–22

In 2016, the Secretary of the U.S. Department of the Interior requested that the U.S. Geological Survey (USGS) produce a publicly available and annually updated database of estimated greenhouse gas emissions associated with the extraction and use of fossil fuels from Federal lands. The first report in this series included emissions estimates from 2005 to 2014 and were reported for 29 States and two offshore areas. Native American and Tribal lands were not included in that analysis. This report recalculates those previous years (2005–14) with updated data and methods and extends the estimates to 2022. Nationwide emissions from fossil fuels produced on Federal lands in 2022 were 1,081.2 million metric tons of carbon dioxide equivalent (MMT CO 2 Eq.) for CO 2 , 33.4 MMT CO 2 Eq. for methane (CH 4 ), and 4.3 MMT CO 2 Eq. for nitrous oxide (N 2 O). Compared to 2005, the 2022 totals represent decreases in emissions for all three greenhouse gases (by 17 percent for CO 2 , 37 percent for CH 4 , and 30 percent for N 2 O). Emissions from fossil fuels produced on Federal lands represent, on average, 21.8 percent of U.S. emissions for CO 2 , 6.1 percent for CH 4 , and 1.3 percent for N 2 O over the 18 years included in this estimate. The trends and relative magnitudes of the greenhouse gas emissions estimated are roughly parallel to the Federal lands fossil fuel production volumes. In 2021, Federal lands of the conterminous United States stored 70,424 MMT CO 2 Eq. in terrestrial ecosystems. Soils stored most of the terrestrial ecosystem carbon (66 percent), followed by live vegetation (25 percent), deadwood (5 percent), and litter (4 percent). Carbon sequestration on Federal lands was highly variable over time, owing primarily to interannual variability in climate and weather, and variability in land use and land cover (LULC) change and disturbances, among these are wildfires and logging. Between 2005 and 2021, Federal lands sequestered an average of 83 MMT CO 2 Eq./yr. By subtracting the cumulative effects of LULC and disturbance-related CO 2 losses to the atmosphere from the total, we estimate that ecosystems at the national level sequestered CO 2 at an annual mean rate of 17 MMT CO 2 Eq./yr in a term called the net ecosystem exchange (NEE). This annual NEE sequestration value represents about 1.4 percent of average fossil fuel emissions over the same period. The USGS estimates presented in this report represent an accounting for the emissions resulting from fossil fuel extraction on Federal lands and the end-use combustion of those fuels, as well as for the sequestration of carbon in terrestrial ecosystems on Federal lands. A combined net CO 2 emissions estimate, which is the difference between the emitted and sequestered CO 2 from both the fossil fuel and ecosystems estimates, provides context for evaluating the greenhouse gas contributions of activities on these lands. The estimates included in this report can provide context for future energy decisions, as well as a basis to track change in the future.

Scientific Investigations Report↗
USGS2025

High-precision earthquake catalog for Minto Flats fault zone, central Alaska, reveals complex and conjugate faulting

The Minto Flats fault zone (MFFZ) in central Alaska is a left‐lateral strike‐slip fault system situated between the continental‐scale right‐lateral Denali and Kaltag‐Tintina faults. The MFFZ has the potential to generate magnitude 7 earthquakes, and it hosted a magnitude 6 earthquake in 1995. It has also produced exotic events, such as very‐low‐frequency earthquakes and nucleation signals. We use network‐matched filtering and relative earthquake relocation techniques to derive a detailed catalog of earthquake locations for the MFFZ. The catalog spans from August 2014 to December 2019, a time period including 13 temporary seismic stations in the region. Our results provide the most complete catalog for the MFFZ and include deeper events, clusters of shallow seismicity, and a complex and segmented fault structure not observed in the original regional catalog. We document right‐lateral strike‐slip faulting, conjugate to the main northeast‐striking left‐lateral faults of the MFFZ. Below Nenana basin, the relocated seismicity reveals northwest‐dipping left‐lateral faults, supporting the inference that deep crustal active faulting is associated with recent basin deformation.

Alaska↗
USGS2026

Geothermal district energy systems coupled with seasonal underground thermal energy storage: A U.S. techno-economic screening by climate and geology

In the United States, cooling-dominated commercial building loads can cause geothermal heat pump-based district energy systems to accumulate a long-term subsurface thermal imbalance, motivating the incorporation of seasonal underground thermal energy storage. We developed a transferable workflow to evaluate geothermal district systems that pair ground heat exchangers with seasonal underground thermal energy storage. Using standardized hourly loads for seven commercial buildings and a uniform cost framework, we simulated ten U.S. cities with a physics-based ground heat exchanger model, subsurface storage simulations, and economic assessment to isolate the roles of climate and hydrogeology. In cooling-dominated cities, underground thermal energy storage supplied the majority of annual cooling, cutting electricity use and summer peaks substantially while achieving levelized costs comparable to or below conventional chiller-boiler plants. In cooler climates, the storage share shrunk, required borefield size and costs rose, and levelized cost of energy increased nearly linearly with declining underground thermal energy storage fraction, indicating storage fraction as the primary economic lever. Sensitivity analysis showed capital risk dominated by borefield drilling and surface heating, ventilation, and air-conditioning and piping, with underground thermal energy storage costs secondary. This workflow provides a transparent foundation for site-specific design and screening of next-generation geothermal district energy systems.

Arizona, Illinois, Maryland, Michigan, Mississippi↗
USGS2026

Hydrogeologic conceptualization of the Arroyos subarea and Northeast Air Force Research Laboratory subarea of the Air Force Research Laboratory, Edwards Air Force Base, Antelope Valley, California, 2017–21

Research and testing within the Air Force Research Laboratory at Edwards Air Force Base in the southwestern Mojave Desert has been conducted since the 1950s. This testing requires large quantities of groundwater for cooling and operations involving hazardous constituents that have entered the groundwater. To better understand the potential movement of contaminants through the groundwater system, the U.S. Geological Survey entered into a cooperative agreement with the U.S. Air Force Civil Engineer Center with the purpose of (1) gaining a more detailed understanding of the hydraulic connection between the unconsolidated basin-fill sediment and the underlying weathered and crystalline bedrock and (2) investigating faults and their potential effect on groundwater flow. This study was specific to the Arroyos subarea and Northeast Air Force Research Laboratory subarea of the Air Force Research Laboratory. The refined hydrogeologic characterization was derived by evaluating surface geologic maps, a regional gravity model, new subsurface geophysical surveys, lithology data from boreholes, and groundwater-level data. Refinement of the hydrogeologic characterization of the study area involved developing a map of basin-fill sediment thickness, estimating the geometry and depth to bedrock, determining the presence of weathered bedrock material, and identifying the location of faults and subsurface structures. The thickness of basin-fill sediment was derived from a regional gravity model. Basin-fill sediment ranges from less than 20 feet (ft) thick to about 150 ft thick and increases substantially toward a deep structural basin in the northwestern part of the study area. The top of bedrock estimated from the new surface geophysical surveys is generally deeper than the top of bedrock from the regional gravity model. These differences in interpreted bedrock tops indicate that there may be a zone of weathered bedrock in the subsurface, overlying crystalline bedrock. Additional subsurface data evaluating the geometry and thickness of weathered bedrock are necessary for a comprehensive evaluation of weathered bedrock in the study area. Previously unmapped extensions of the Leuhman Fault, Spring Fault Zone, and Boron Gate Fault were identified in the new surface geophysical surveys. Evaluation of groundwater flow between bedrock and basin-fill sediment, and the effects of faults on groundwater flow, was done using groundwater-level data from 2020. The groundwater table was in crystalline bedrock throughout most of the Arroyos and Northeast Air Force Research Laboratory subareas, with a few exceptions where it was in basin-fill sediment and weathered bedrock. Basin-fill sediment, and likely weathered bedrock, was partially saturated in parts of the study area where basin-fill sediment was as thick as 150 ft. Weathered bedrock material may be present near the Spring Fault Zone, and if so, was likely saturated in 2020. The groundwater-level data near and within Spring Fault Zone indicate that the Spring Fault Zone acts as a vertical-along-strike conduit for, and as an across-strike impediment to, groundwater flow. In addition, steep groundwater-level gradients in the southeastern part of the study area indicate that a subsurface structure affecting groundwater flow could be related to the Leuhman Fault. Sufficient groundwater-level data were not available to conclusively determine if the Arroyos, Rich, and Boron Gate Faults, and their inferred extensions from the geophysical surveys, affect the movement or direction of groundwater flow.

California↗
USGS2025

Using gridded seismicity to forecast the long-term spatial distribution of earthquakes for the 2025 Puerto Rico and U.S. Virgin Islands National Seismic Hazard Model

Gridded (or background) seismicity models are a critical component of probabilistic seismic hazard assessments, accounting for off‐fault and smaller‐magnitude earthquakes. They are typically developed by declustering and spatially smoothing an earthquake catalog to estimate a long‐term seismicity rate that can be used to forecast future earthquakes. Here, we present new gridded seismicity models for use in the 2025 National Seismic Hazard Model (NSHM) for Puerto Rico and the U.S. Virgin Islands (PRVI). The previous PRVI NSHM was released in 2003, and our new models incorporate updates to both data and methodology. We utilize an updated earthquake catalog based on improved Puerto Rico Seismic Network data with newly characterized completeness epochs. The catalog is divided into crustal, subduction interface, and intraslab seismicity using new methods and Slab2 subduction zone geometries. To forecast the long‐term spatial distribution of earthquakes, we use an updated methodology developed for the 2023 U.S. 50‐state NSHM, considering three declustering methods and two spatial smoothing methods based on 2D Gaussian kernels. To adapt it for the complex seismotectonics of the region, we also adopt probabilistic methods to account for events with unknown depths and uncertainties in tectonic classification, and develop a new method for spatial scaling to counteract the effects of spatial variability in network coverage while maintaining the use of smaller events. Finally, we test the performance of these spatial models in forecasting the location of M w ≥ 5earthquakes in the region. Our updated methodology improves the representation of epistemic uncertainty relative to the 2003 model, and our results demonstrate the effectiveness of the new measures we have introduced to address heterogeneities in network detection and systematically evaluate forecast performance.

Puerto Rico, U.S. Virgin Islands↗
USGS2024

Quantitative risk of earthquake disruption to global copper and rhenium supply

Earthquakes have the potential to substantially affect mining operations, potentially leading to supply chain disruptions and adversely affecting the global economy. This study explores the quantification of earthquake risk to copper and rhenium commodity supply by examining the spatial concentration of high earthquake hazard areas and the commodity-specific mining, smelting, and refining operations across the globe. Because many of the largest facilities are concentrated geographically near the highly seismic regions of South America, East Asia, and the Pacific, there is a potential for cascading effects on the entire supply chain. The analysis indicates that the expected annual disruption of global production is 0.3–1.1 percent for copper mines, 1.8–4.0 percent for smelters, and 1.5–3.3 percent for refineries. Expected annual disruption of global rhenium production capacity is 0.32–1.32 percent. The research highlights that the potential lost revenue from earthquake disruptions is from $315 million to $1.29 billion for copper mining, $1.92 billion to $4.33 billion for copper smelting, $2.06 billion to $4.52 billion for copper refining, and $337,000 to $1.40 million for rhenium production capacity.

Open-File Report↗
USGS2025

Geologic map of MTM −15032 and −20032 quadrangles, western Ladon basin, Mars

Our mapping efforts focused on the geomorphic terrains and geologic units contained within Mars Transverse Mercator (MTM) −15032 and −20032 quadrangles. These two quadrangles are located along the west side of Ladon basin and span lat −12.5° N. to −22.5° N. and long 325° E. to 330° E. The western part of Ladon basin and its bounding basin ring structures to the west preserved features that help to understand the long history of drainage across the Margaritifer Terra region of Mars. Our geologic map provides new insight into the extent of aqueous activity across western Ladon basin and the adjoining highlands, which includes the northern part of Ladon Valles, Arda Valles, numerous small valleys in the western highlands, and light-toned layered deposits associated with these fluvial features. The Mars Odyssey Thermal Emission Imaging System (THEMIS) infrared (IR) daytime mosaic (100 meters per pixel) was used as the primary base map. We constructed the geologic map of western Ladon basin at 1:1,000,000 scale. We identified 20 geologic units in the map area, which we divided into the following groups: crater units; volcanic units; chaotic units; basin fill units; crater, valley, and channel units; and plateau and highlands units. We observed the following sequence of events: (1) formation of Ladon and Holden impact basins in the middle Noachian, producing a mountainous unit; (2) landscape degradation and infilling during the Late Noachian to Late Hesperian forming terra units; (3) formation of Ladon Valles in the Late Noachian to Early Hesperian by catastrophic flooding, producing channel units; (4) accumulation of sediments in the Late Noachian to Early Amazonian, forming basin fill units; (5) formation of a smooth crater fill unit during the Early to Late Hesperian; (6) deposition during the Late Hesperian to Early Amazonian, creating light-toned layered units within Ladon Valles, Ladon basin, and other smaller valley networks along the western uplands, as well as formation of a light-toned unit along some crater floors; (7) formation a of chaotic unit in the Early to Middle Amazonian, formation of alluvial fans in the Early Amazonian, and eruption of a volcanic unit in the Middle Amazonian; and (8) formation of craters throughout the geologic history of the map region.

Scientific Investigations Map↗
USGS2024

Satellite interferometry landslide detection and preliminary tsunamigenic plausibility assessment in Prince William Sound, southcentral Alaska

Regional mapping of actively deforming landslides, including measurements of landslide velocity, is integral for hazard assessments in paraglacial environments. These inventories are also critical for describing the potential impacts that the warming effects of climate change have on slope instability in mountainous and cryospheric terrain. The objective of this study is to identify slow-moving landslides in the Prince William Sound region, southcentral Alaska, United States, which has had rapid deglaciation since the mid-1800s, and assess their tsunamigenic plausibility. We use an automated time series persistent scatterer interferometric synthetic aperture radar processing method with 7 years of Sentinel-1 data (2016–22) to identify 43 slow-moving slopes with average velocities ranging from approximately 0.2 to 21 millimeters per year. Landslide presence is confirmed using aerial imagery and previous landslide inventory records. We assess the tsunamigenic plausibility of the landslides using empirically derived estimates of landslide mobility based on modeled landslide volumes. Of the identified landslides, our preliminary analysis suggests that 11 have tsunamigenic potential if they were to fail rapidly and catastrophically. Although our estimate of tsunamigenic plausibility is preliminary and can be refined with additional observations and analyses, it can be used to prioritize ongoing and future hazard assessment, surveillance, and research efforts.

Alaska↗
USGS2025

Integrated stratigraphic and geochemical analysis of organic-rich intervals of the Lewis Shale in the eastern Washakie Basin, Wyoming

Geological studies in the Cretaceous Western Interior Seaway (KWIS) in North America reveal highly variable sedimentological conditions on its western shore caused by rapidly changing sea level and detrital input during the seaway’s closure. Here we examine a 601-foot (183 meters) continuous core through the lower part of the Maastrichtian Lewis Shale in the eastern Washakie Basin, Wyoming, through integrating stratigraphic and geochemical analyses to better understand organic matter deposition and preservation during the final marine transgression within the seaway. The core penetrates eight organic-rich flooding intervals (F1-F8) and a regional condensed section, informally known as the Asquith marker. The lower portions of the core record sediment input from a southern source, likely the Sierra Madre/Park Ranges, while the upper part records sediment input from a northern source, likely the Granite Mountains. This provenance transition is supported by regional stratigraphic analysis and changes in bulk mineral and trace metal composition. The Asquith marker and early flooding surfaces are enriched in oil-prone, marine organic matter deposited under dysoxic to euxinic conditions, whereas younger flooding surfaces show increased terrigenous input and poorer preservation conditions. The Asquith marker is identified as a prime oil-prone source rock. Although younger flooding surfaces also exhibit favorable source-rock properties, their generative potential is reduced due to increased clastic and terrigenous organic matter deposition caused by regression. These results emphasize the importance of sediment source variability, organic matter preservation, and changing redox conditions to provide insight into the sediment provenance and petroleum potential of the Lewis Shale in the context of the final transgression within the KWIS.

Colorado, Wyoming↗
USGS2026

CRESCENT earthquake dynamic rupture, earthquake cycle, and tsunami code verification platform

Physics-based simulations are critical for understanding natural hazards. The increasing complexity of numerical codes requires benchmark exercises to verify that different computational methods yield consistent results when solving the same governing equations. Here, we present an open-access web platform designed for the verification of earthquake dynamic rupture, seismic cycle, and tsunami simulations. The platform architecture utilizes a modular, serverless backend on Amazon Web Services (AWS) to provide scalable file processing and visualization. A lightweight static web application provides a secure interface for uploading and managing results, while the browser-based data visualization enables interactive analysis of time series and surface grid data. By using structured JavaScript Object Notation (JSON) text files to define benchmark structures, the system remains fully extensible, allowing the addition of new scenarios without modifying the underlying software logic. The platform hosts the "The Tsunami Problem Versions" (TTPV) 1 & 2, two benchmarks for 3D fully coupled earthquake dynamic rupture and tsunami generation, and provides a framework for earthquake cycle models. This community resource aims to build trust in numerical simulations and facilitate long-term collaborative code verification as modeling software continues to evolve.

Seismica↗
USGS2026

The Sand AbrasioN Device for Aeolian Research (SANDAR): A new experimental device for investigating how wind transport affects sand on Earth and Mars

On Earth and Mars, aeolian transport causes sand grains to become abraded, resulting in mineralogic and textural changes. Understanding how sands evolve, or mature, with transport via experimental studies is important for understanding the origins, geologic history, and cycling of sediments, as well as dust production. Previous experimental works have used a variety of methods to simulate aeolian transport in the laboratory, but practical limitations and similitude concerns have limited such research. Here, we present and validate the Sand AbrasioN Device for Aeolian Research (SANDAR), a modified air mill that uses pressurized air to circulate sand around a small abrasion chamber, simulating the effects of aeolian transport. This device is re-circulating to simulate long-distance transport, and it allows for repeated analyses of well-constrained sediment samples, revealing their evolution over time. It is compatible with the grain sizes (74–500 μm) and grain impact velocities (∼0.6–3.7 m/s) typically expected for natural aeolian environments, and is also adaptable for diverse applications simulating different wind conditions. We show that the SANDAR achieves similitude of kinetic energy with respect to saltating sand on both Earth and Mars. SEM and optical microscope imaging reveal that the SANDAR produces microtextures on the surfaces of sand grains similar to those found with natural aeolian transport, demonstrating that it effectively simulates the mechanical effects of aeolian processes. Thus, the SANDAR is a valid tool for use in experimental research to improve our understanding of sedimentary processes across the Solar System.

Aeolian Research↗
USGS2025

Geologic map of the Greater Antilles and Virgin Islands

Introduction The geologic map of the Greater Antilles and Virgin Islands is a compilation of information from the literature, integrated to provide a seamless geologic map of the region. This map was prepared to serve as a base map for a mineral resource assessment of the region. Several small-scale regional geologic maps of the region have been prepared in the past. This report supersedes an earlier version of the geologic map of the Greater Antilles and the Virgin Islands that was released as U.S. Geological Survey (USGS) Open-File Report 2019–1036. For this report, the regional geologic overview shown on sheet 1 covers Cuba; the island of Hispaniola, which includes Haiti and Dominican Republic; Jamaica; the Cayman Islands; Puerto Rico; and the U.S. and British Virgin Islands. Sheet 2 shows the geology of Cuba and Cayman Islands, sheet 3 the geology of Jamaica, sheet 4 the geology of Hispaniola, and sheet 5 the geology of Puerto Rico and Virgin Islands. Accompanying the maps are three tables: table 1 lists the map units (in alphabetical order by map-unit label), their ages, and the countries or territories in which they can be found. These are listed in alphabetical order for ease in locating unit names from the symbols on the maps. Table 2 is a list of sources for the geologic map compilation of Puerto Rico. Table 3 shows formations in the Devil’s Race Course Group of Jamaica.

Greater Antilles, Virgin Islands↗
USGS2024

New occurrences of the rare, REE minerals daqingshanite, törnebohmite, biraite, sahamalite, and ferriperbøeite from the Sheep Creek area, Montana, USA

Over 30 small, discontinuous, tabular carbonatite bodies are located in the Sheep Creek area, Ravalli County, southwest Montana. The age and origin of these REE-Nb-rich deposits are currently being investigated. The purpose of this paper is to document the occurrence of several rare minerals, including daqingshanite, törnebohmite, biraite, sahamalite, and ferriperbøeite, in two of the carbonatite bodies. These minerals are found in association with monazite, hydroxylbastnäsite, ferriallanite, calcite, dolomite, baryte, quartz, actinolite, apatite, celsian, and Sr-rich aragonite. Automated SEM-EDS was used to target the areas of interest in polished specimens for more detailed spot SEM-EDS and electron probe microanalysis. Raman spectra were also acquired for each of the rare minerals. The complex mineralogy of the Sheep Creek carbonatites is most likely due to several overlapping thermal events, including primary magmatic, overprinting hydrothermal, and supergene weathering stages. The rare minerals described in this study are believed to be hydrothermal and/or carbothermal in origin, although no estimates of temperature are available at this time.

Montana↗
USGS2025

Oligocene–Miocene development and evolution of the south Dome Rock Mountains basin, lower Colorado River corridor, Arizona, USA

Sedimentary basins in the Colorado River extensional corridor record large-magnitude Basin and Range extension and younger dextral shear deformation in the evolving Pacific−North America plate boundary. The south Dome Rock Mountains basin is located in west-central Arizona (USA), where the history of basin evolution, style of deformation, and timing of the transition between extension and dextral shear are not well constrained. We integrate new geologic mapping of the south Dome Rock Mountains basin with zircon U-Pb geochronology of six marker beds to characterize the timing of basin evolution and the slip history of the south Dome Rock Mountains normal fault. Structural analysis defines structures consistent with regional extension and younger dextral shear deformation. We use sedimentological and petrological analysis to interpret the depositional environments of three basin sequences. A lower basin sequence consists of fluvial strata deposited in an internally drained intermontane basin from ca. 35 Ma to 24.4 Ma. A sequence of volcanic rocks was emplaced in the basin between 24.4 Ma and 23.3 Ma, before the onset of local extensional faulting. An upper basin sequence of coarse conglomerate and sedimentary breccia was deposited on the hanging wall of a half graben structure after 23.3 Ma and through ca. 12.8 Ma, synchronous with 3−7 km of dip slip on the west-dipping, listric south Dome Rock Mountains fault. The basin was subsequently deformed by dextral shear after 12.8 Ma and before 4.8 Ma. These results document how distributed deformation related to the evolving Pacific−North America plate boundary occurred >100 km from the primary plate boundary.

Arizona, California↗
USGS2025

Mineral commodity summaries 2025

Introduction Each mineral commodity chapter of the 2025 edition of the U.S. Geological Survey (USGS) Mineral Commodity Summaries (MCS) includes information on events, trends, and issues for each mineral commodity as well as discussions and tabular presentations on domestic industry structure, Government programs, tariffs, 5-year salient statistics, and world production, reserves, and resources. The MCS is the earliest comprehensive source of 2024 mineral production data for the world. More than 90 individual minerals and materials are covered by two-page synopses. Abbreviations and units of measure and definitions of selected terms used in the report are in Appendix A and Appendix B, respectively. Reserves and resources information is in Appendix C, which includes “Part A—Resource and Reserve Classification for Minerals” and “Part B—Sources of Reserves Data.” A directory of USGS minerals information country specialists and their responsibilities is in Appendix D. The USGS continually strives to improve the value of its publications to users. Constructive comments and suggestions by readers of the MCS 2025 are welcomed.

Mineral Commodity Summaries↗
USGS2026

Mineral commodity summaries 2026

Introduction Each mineral commodity chapter of the 2026 edition of the U.S. Geological Survey (USGS) Mineral Commodity Summaries (MCS) includes information on events, trends, and issues for each mineral commodity as well as discussions and tabular presentations on domestic industry structure, Government programs, tariffs, 5-year salient statistics, and world production, reserves, and resources. The MCS is the earliest comprehensive source of 2025 mineral production data for the world. More than 90 individual minerals and materials are covered by two-page synopses. Abbreviations and units of measure and definitions of selected terms used in the report are in Appendix A and Appendix B, respectively. Reserves and resources information is in Appendix C, which includes “Part A—Resource and Reserve Classification for Minerals” and “Part B—Sources of Reserves Data.” A directory of USGS minerals information country specialists and their responsibilities is in Appendix D. The USGS continually strives to improve the value of its publications to users. Constructive comments and suggestions by readers of the 2026 MCS are welcomed.

Mineral Commodity Summaries↗
USGS2024

Exploring landscape and geologic controls on spatial patterning of streambank groundwater discharge in a mixed land use watershed

Preferential groundwater discharge features along stream corridors are ecologically important at local and stream network scales, yet we lack quantification of the multiscale controls on the spatial patterning of groundwater discharge. Here we identify physical attributes that best explain variation in the presence and lateral extent of preferential groundwater discharges along two 5th order streams, the Housatonic and Farmington Rivers, and 32 1st to 4th order reaches across the Farmington River network. We mapped locations of preferential groundwater discharge exposed along streambanks using handheld thermal infrared cameras paired with high-resolution topographic and land use land cover datasets, surficial soil characteristic maps, and depth-to-bedrock geophysical measurements. The unconfined Housatonic River, MA, USA (12 km) had fewer discharge locations and less lateral extent (41 discharge locations with 38 m of active discharge/km of river) compared to the partially confined Farmington River, CT, USA (26 km; 169 discharge locations with 129 m of active discharge/km of river). Using a moving window analysis, we found along both rivers that discharge was more likely to occur where bank slopes were steeper, floodplain extent was narrower, and degree of confinement was higher. Along the Farmington River, groundwater discharge was more likely to occur where saturated hydraulic conductivity was higher and depth-to-bedrock was shallower. Among the 32 stream reaches surveyed (33.2 km of total stream length) within the Farmington River watershed, preferential discharge was observed in all but two stream reaches, varied from 0 to 25% of lateral extent along stream banks (mean = 6%), and was more likely to occur where stream reach slopes were steep, saturated hydraulic conductivity was high, and watershed urbanization was low. Our results show that, though both surface (e.g., topographic, land use land cover) and subsurface (e.g., soil characteristics, bedrock depth) factors control the prevalence of streambank preferential groundwater discharge, the dominant controls vary across valley settings and stream sizes.

Connecticut, Massachusetts↗