USGS Science⌕ Search

SEARCH · USGS Science

Results for “Communications Earth & Environment”

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.

512 records · Page 27Linked to original sources
USGS2025

Widespread occurrence of former anhydrite phenocrysts in Laramide-age magmas related to porphyry-skarn Cu mineralization at Santa Rita and Hanover-Fierro, New Mexico, USA

Reports of magmatic anhydrite are relatively rare, with only ~30 occurrences documented worldwide so far. However, magmatic anhydrite saturation is difficult to recognize because anhydrite decomposes rapidly in near-surface environments. In most cases, only anhydrite inclusions shielded within other phenocryst phases were able to survive. Alternatively, since anhydrite phenocrysts preserved in fresh volcanic rocks are characteristically intergrown with apatite phenocrysts, the former presence of anhydrite phenocrysts can be recognized based on the occurrence of lath-shaped cavities that show a strong spatial association with apatite phenocrysts. These cavities can be either empty or filled with low-temperature, secondary minerals such as zeolites, carbonates, or microcrystalline silica. A systematic search for the occurrence of such cavities, combined with optical and Raman-spectroscopic identification of anhydrite inclusions preserved within apatite, hornblende and quartz phenocrysts, demonstrates that most of the Laramide-age magmas associated with the Santa Rita and Hanover-Fierro porphyry-skarn Cu (Zn, Mo, Au, Pb) deposits were saturated in magmatic anhydrite. The anhydrite typically coexisted with monosulfide solid solution (MSS), suggesting oxygen fugacities of ~2.0±0.5 log units above the fayalite-magnetite-quartz buffer. The magmas range from andesitic to rhyodacitic in composition, and from shortly pre-mineralization (~61 Ma) to shortly post-mineralization (~57 Ma) in age. In three samples with particularly well-recognizable former anhydrite phenocrysts, their modal abundance could be quantified based on high-resolution scans of polished hand specimens. The observed modal anhydrite abundances of 0.63–1.8 vol% translate into minimum magma sulfur contents of 0.20–0.56 wt% S. The highest sulfur content of 0.56 wt% S is difficult to reconcile with available anhydrite solubility models, but it could be reproduced in an anhydrite solubility experiment performed at 950 °C and 1.15 GPa on a natural latite containing 13.1 wt% dissolved H2O. The sample with the second-highest sulfur content of 0.26 wt% S requires ~10 wt% H2O in the silicate melt, and, consequently, a minimum pressure of ~0.5 GPa. Taken together, the results suggest that the magmas of the Central Mining District were extremely hydrous and thus originated from great depth. Indeed, their major element compositions and reconstructed H2O and S contents agree well with experimentally observed and numerically predicted compositions of residual silicate melts after 50–70 wt% crystallization of ordinary arc basalts at high pressure and high oxygen fugacities.

New Mexico↗
USGS2025

The 3D Elevation Program—Supporting the economy of Massachusetts

Introduction Massachusetts extends from the mountains of the Appalachian system in the west of the State to the sandy beaches and rocky shorelines of the Atlantic coast in the east. Inland topographic data support a wide range of important activities, including geologic mapping, transportation planning, forest and wildlife management, quantifying ecological services, water supply protection, commonwealth-wide infrastructure planning, local site planning, and flood-plain management. Nearshore bathymetry can be used to support coastal portions of the Commonwealth by addressing the combined threats of ocean warming, strong storm surge, and rising sea levels. The maintenance and (or) expansion of Massachusetts ports (for instance, Boston, New Bedford) and Cape Cod sediment management depends upon the accurate mapping of bathymetry and the frequent influx of sediment and redeposition. Critical applications that address the broad range of requirements depend on light detection and ranging (lidar) data that provide a highly detailed three-dimensional (3D) model of the Earth’s surface and aboveground features. The 3D Elevation Program (3DEP) is managed by the U.S. Geological Survey (USGS) in partnership with Federal, State, Tribal, U.S. territorial, and local agencies to acquire consistent lidar coverage at quality level 2 or better to meet the many needs of the Nation and Massachusetts. The status of available and in-progress 3DEP baseline lidar data in Massachusetts is shown in figure 1. 3DEP baseline lidar data include quality level 2 or better, 1-meter or better digital elevation models, and lidar point clouds, and must meet the Lidar Base Specification version 1.2 ( https://www.usgs.gov/3dep/lidarspec ) or newer requirements. The National Enhanced Elevation Assessment identified user requirements and conservatively estimated that availability of lidar data would result in at least $1.23 million in new benefits annually to Massachusetts. The top 10 Massachusetts business uses for 3D elevation data, which are based on the estimated annual conservative benefits of 3DEP, are shown in table 2.

Massachusetts↗
USGS2025

SURF: An automated method for building nonplanar 3D fault models from earthquake hypocenters

Accurately characterizing 3D fault geometry is vital for improving our understanding of earthquake behavior and informing the development of seismic hazard models. Despite their importance, subsurface fault structures tend to be poorly constrained because of limitations in observational data. Improvements to the seismic networks and earthquake detection algorithms have increased the precision and volume of earthquake catalogs, which help illuminate detailed subsurface fault structure and provide the most direct information available about fault geometries at depth. We present a Python package to automate generating 3D fault geometries directly from hypocentral seismicity patterns. This method begins with clustering events based on their spatial density, identifying coherent patterns. Nearby clusters are then merged based on the similarity of their orientations. We fit nonplanar surfaces using support vector regression to balance surface accuracy with minimal deviations from planarity. The fault models are output as quadrilateral meshes at user‐defined resolution. In the process of generating the 3D fault surfaces, we compute the spatial density of seismicity around the surface and the planarity as quantitative metrics of the model outputs. As a proof of concept, we apply this approach to the San Andreas–Calaveras fault junction region and the 2019 Ridgecrest earthquake sequence, both in California, which contain complex subparallel faults well defined at the Earth’s surface and abundant microseismicity. These case studies demonstrate the method’s ability to model complex fault structures, including long continuous fault surfaces, crossing faults, variably dipping segments, and subparallel faults. We test the method on both standard network catalogs and double‐difference relocated catalogs. We find that our seismicity‐based fault model results align with published 3D models that incorporate additional constraints and interpretations ( Plesch et al. , 2020 ; Aagaard and Hirakawa, 2021 ). This workflow provides a low‐user‐input solution for estimating fault geometries at depth from earthquake catalogs.

California↗
USGS2026

Simulation-based scenario ShakeMaps for large magnitude (MW6.5+) crustal earthquakes on the Seattle, Tacoma, and southern Whidbey Island faults, Washington, USA

Scenario ground‐motion maps based on empirical ground‐motion models (GMMs) provide a rapid and generally reliable means of estimating the amplitude and distribution of earthquake shaking. However, because GMMs are designed for broad applicability, they often rely on simplified representations of Earth structure, which can limit their accuracy in regions with complex source, path, and site effects. This can substantially impact the accuracy of predicted shaking in areas like western Washington State, where deep, interconnected basin structure exerts a strong influence on seismic‐wave propagation. In this study, we present a new suite of simulation‐based scenario ShakeMaps that characterize ground shaking from large‐magnitude ( ⁠ M W 6.5–7.5) earthquakes on the Seattle, Tacoma, and southern Whidbey Island faults. These maps are developed using results from recent 3D wave propagation simulations ( Stone et al. , 2022 , 2023 , 2025 ) that incorporate realistic rupture geometries, variable slip distributions, and a regional 3D seismic velocity model with shallow soils. Broadband ground motions are estimated by combining the low‐frequency (<1 Hz) deterministic seismograms from these studies with high‐frequency (1–10 Hz) stochastic seismograms. Simulated ground motions are corrected to account for the enforced minimum shear‐wave velocity and nonlinear site response. The resulting ShakeMaps represent median ground‐shaking estimates derived from multiple rupture scenarios with varying slip distributions and hypocenter locations for each fault. To extend ShakeMap coverage beyond the simulation domain (i.e., into eastern Washington, northern Oregon, and southwestern British Columbia), we scale GMM‐based ground‐motion estimates using amplification patterns observed in the simulations. These new ShakeMaps reveal the substantial influence of deep basin structure on shaking intensity, underscoring the importance of considering crustal structure complexity in regional hazard assessments for the Pacific Northwest.

Washington↗
USGS2026

Recent scientific contributions by the U.S. Geological Survey in the San Francisco Bay and Sacramento–San Joaquin Delta Estuary

Introduction The San Francisco Bay and Sacramento–San Joaquin Delta Estuary (hereafter, Bay-Delta) is the largest estuary on the West Coast of the United States. The Bay-Delta covers more than 1,600 square miles and drains a watershed of more than 75,000 square miles, which is greater than 40 percent of California. The region surrounding the Bay-Delta is home to about 10 million people, and its habitats (fig. 1) support more than 800 plant and animal species. The waterways of the Bay-Delta are the central hub of California’s extensive freshwater delivery system, supplying water to more than 27 million Californians and 4 million acres of farmland in the Central Valley. The U.S. Geological Survey (USGS) is the primary science agency of the U.S. Department of the Interior, providing a broad range of Earth, water, biological, and mapping data and expertise to inform natural resource management across the country. This fact sheet focuses on research conducted by the USGS in the Bay-Delta region, mostly within the past 5 years. The fact sheet is organized across five major themes: water flow and water quality, fish and wildlife, wetland restoration, invasive species, and hazards. In each of these areas, the USGS works closely with Federal, State, and local agencies; academia; and non-governmental organizations to provide objective, evidence-based science. The data and knowledge gained from USGS research and monitoring are publicly available and are routinely used by agencies, including the U.S. Bureau of Reclamation, U.S. Army Corps of Engineers, U.S. Fish and Wildlife Service, California Department of Water Resources, California Department of Fish and Wildlife, California State Water Resources Control Board, Delta Stewardship Council, and many additional organizations and stakeholders to inform management of the Bay-Delta region and its resources.

California↗
USGS1999

Water quality of Rob Roy Reservoir and Lake Owen, Albany County, and Granite Springs and Crystal Lake Reservoirs, Laramie County, Wyoming, 1997-98

The water quality of four reservoirs was assessed during 1997 and 1998 as a cooperative project between the Cheyenne Board of Public Utilities and the U. S. Geological Survey. The four reservoirs, Rob Roy, Lake Owen, Granite Springs, and Crystal Lake, provide approximately 75 percent of the public water supply for Cheyenne, Wyoming. Samples of water and bottom sediment were collected and analyzed for selected physical, chemical, and biological characteristics to provide data about the reservoirs. Water flows between the reservoirs through a series of pipelines and stream channels. The reservoirs differ in physical characteristics such as elevation, volume, and depth. Profiles of temperature, dissolved oxygen, specific conductance, and pH were examined. Three of the four reservoirs exhibited stratification during the summer. The profiles indicate that stratification develops in all reservoirs except Lake Owen. Stratification developed in Rob Roy, Granite Springs, and Crystal Lake Reservoirs by mid-July in 1998 and continued until September, with the thickness of the epilimnion increasing during that time. Secchi disk readings indicated Rob Roy Reservoir had the clearest water of the four reservoirs studied. The composition of the phytoplankton community was different in the upper two reservoirs from that in the lower two reservoirs. Many of the species found in Rob Roy Reservoir and Lake Owen are associated with oligotrophic, nutrient-poor conditions. In contrast, many of the species found in Granite Springs and Crystal Lake Reservoirs are associated with mesotrophic or eutrophic conditions. The total number of taxa identified also increased downstream. The chemical water type in the reservoirs was similar, but dissolved-solids concentrations were greater in the downstream reservoirs. Water in all four reservoirs was a calcium-bicarbonate type. In the fall of 1997, Rob Roy Reservoir had the lowest dissolved-solids concentration (19 milligrams per liter), whereas Crystal Lake Reservoir had the highest concentration (63 milligrams per liter). Relatively little differences in the concentrations of major-ion species were noted between samples collected near the surface and near the bottom of the same reservoir. In contrast, iron and manganese concentrations generally were higher in samples collected near the bottom of a reservoir than in near-surface samples collected from the same reservoir. Composite bottom-sediment samples from all four reservoirs contained similar concentrations of bulk constituents such as aluminum, iron, phosphorus and titanium, but varied in concentrations of trace elements. Trace-element concentrations in Rob Roy Reservoir and Lake Owen were similar to the crustal average, whereas in Granite Springs and Crystal Lake Reservoirs the concentrations were similar to granitic rocks.

Wyoming↗
USGS2026

Archean to Mesozoic–Cenozoic seismic crustal structure: Implications for geological and biological evolution

We use >4500 measurements of crustal structure to investigate the seismic structure of continental crust, Archean to Mesozoic–Cenozoic. The mean crustal thickness of continents, including their margins, is 36.5 km. We find that Archean, Paleoproterozoic, and Mesoproterozoic crust have similar mean crustal thickness (41 km), seismic velocities, Vp/Vs ratio, and density. Crusts of these ages span 3 Ga, from 4.0 to 1.0 Ga, and their similarity of physical properties (thickness, density, and seismic velocities) suggests that the process of crustal formation may have been similar during this time period, which covers 66% of Earth's history. Neoproterozoic crustal properties differ significantly from Mesoproterozoic and older crusts. The mean thickness of Neoproterozoic crust is 32 km, some 9 km thinner than the mean thickness of Archean, Paleoproterozoic, and Mesoproterozoic crust. A lithospheric root with a thickness of 150–200 km underlies Archean, Paleoproterozoic, and Mesoproterozoic crusts, and this root resists lithospheric rifting and crustal extension and thinning. Based on these observations, particularly the lithospheric thickness, we conclude that Archean, Paleoproterozoic, and Mesoproterozoic lithosphere are unique and together form the stable cratonic nuclei, defined as the thick (150–250 km), long-lived cores of continents. Higher mantle temperatures during the Archean, Paleoproterozoic, and Mesoproterozoic may have played a key role in the formation of the thick lithospheric roots. The second key finding is that the similarities of mean seismic properties indicate that the process of crustal formation operating in the Archean continued during the Paleoproterozoic and Mesoproterozoic. The thin (mean value < 135 km) lithospheric root beneath Neoproterozoic and younger crust may be related to the steady decrease in mantle temperature through time. Neoproterozoic and Paleozoic crust have similar physical properties, and these eras are characterized by pronounced biodiversification, including the renowned Garden of Ediacara, the Cambrian Explosion, and the Great Ordovician Biodiversification. Mesozoic–Cenozoic crust is the most diverse and reflects the current tectonic and magmatic processes of crustal formation.

Book chapter↗
USGS2026

Tectonic implications of transitional melting regimes from petrological, geochronological, and compositional characterization of the ophiolitic Seventymile terrane, Alaska, USA

New geochemical, U-Pb geochronology, and Sr-Nd-Hf isotope data provide evidence for the tectonic evolution of the Seventymile terrane in interior Alaska, USA. Ultramafic and mafic rocks of the Seventymile terrane are thought to represent components of a dismembered ophiolite and provide unique constraints on regional terrane evolution and accretion. The Seventymile ophiolite represents fragments of the Devonian to Permian Slide Mountain Ocean (SMO) that separated allochthonous and parautochthonous continental fragments of western North America. It now occurs as multiple thrust sheets containing Permian mafic and ultramafic rocks overlying and/or possibly imbricated with amphibolite-facies supracrustal rocks of allochthonous Yukon-Tanana terrane and parautochthonous North America. Seventymile klippen contain variably serpentinized peridotite, primarily harzburgite, low-grade meta-mafic rocks, and minor oceanic sedimentary rocks (argillite, chert, limestone, and metasandstone). Mafic rocks include gabbro to diabase, typically as dikes, veinlets, or rare massive stocks intruding peridotite. Mafic rocks also include greenstones of the Seventymile assemblage in klippen structurally underlying, and in shear zone contact with, ultramafic klippen. New trace element and radiogenic isotope data from mafic magmatic rocks associated with the Seventymile ultramafic bodies show evidence for a weakly subduction-modified mantle source, like the mantle source of normal mid-ocean-ridge basalt (N-MORB) or back-arc basin basalt (BABB). Seventymile assemblage greenstones are more heterogeneous. They range from N-MORB to enriched mid-ocean-ridge basalt (E-MORB) and ocean-island basalt (OIB), with a subset of samples indicative of continental arc affinity. These geochemistry results indicate that distinct tectonic environments are represented by at least two, and possibly three, lithological and structural units comprising the Seventymile terrane. Hf-Nd isotope systematics are consistent with a depleted MORB mantle (DMM)−like component that overlaps with Pacific MORB. Primary zircon is rare, but new in situ U-Pb data for gabbro and greenstone indicate ca. 274−272 Ma peak zircon and titanite crystallization. Scattered younger zircons define a ca. 255 Ma zircon peak and correspond to secondary crystallization associated with baddeleyite reaction of high-Si fluids during low-grade metamorphism. If Seventymile suites are contemporaneous, obduction associated with the closure of the SMO resulted in the stacking of ophiolitic packages representing distinct tectonomagmatic settings across the transition from pericontinental, to epicontinental, to distal ocean back-arc. Intrusions hosted in klippe of ultramafic rocks, plus the least subduction-modified greenstones underlying them, geologically and compositionally resemble Slide Mountain rocks of the Campbell Range formation in eastern Yukon and may provide a new piercing point across the Tintina fault.

Alaska, British Columbia↗
USGS2026

Ground-motion simulations for the 2024 Mw 4.8 Tewksbury, New Jersey, earthquake

Ground-motion simulations of notable earthquakes in the central and eastern United States are limited and typically assume one-dimensional (1D) Earth structure. In this study, we use a three-dimensional (3D) seismic velocity model to better constrain the depth and focal mechanism of the April 5th, 2024, moment magnitude 4.8 Tewksbury earthquake and investigate the spatial variability of earthquake ground motions and the effects of nearby sedimentary basins. We perform earthquake ground-motion simulations up to 0.5 Hz using the 3D spectral-element wave-propagation solver SPECFEM3D over a region 280-km wide by 260-km long by 77-km deep. Topography and subsurface geophysical structure are assigned using the U.S. Geological Survey National Crustal Model with a minimum shear-wave velocity of 200 m/s. We use earthquake time series from 13 broadband seismic stations in the region that have a uniform azimuthal distribution and epicentral distances ranging from 76 to 131 km to compare with synthetics and explore the effects of 1D versus 3D seismic structure on focal mechanism and depth solutions. Ground-motion intensity metrics are also presented relative to the NGA-East ground-motion models (GMMs) currently used in seismic hazard assessments for the region. We find that the 3D model, which reveals a wide spatial variability of period-dependent ground motions, yields better predictions of earthquake ground motions relative to the 1D model and the NGA-East ergodic ground-motion model, with 76 percent reduction of residual variance in observed ground motions averaged over 3-, 5-, 7-, and 10-second periods. Use of the 3D model to solve for a focal mechanism yields a shallower focal depth at 4 km and a shallower east-dipping focal plane relative to the U.S. Geological Survey regional moment tensor and Global Centroid Moment Tensor. Our study demonstrates that use of 3D seismic velocity models can improve estimates of earthquake focal mechanisms, ground motions, and seismic hazard.

New Jersey↗
USGS2025

Assessment of coastal and fluvial morphodynamic changes using Structure-for-Motion: A case study of the Sfȃntu Gheorghe Mouth (Danube Delta, Romania)

The ability to accurately map erosion, flooding, and habitat loss in coastal environments is crucial for formulating national strategies aimed at preventing and mitigating the impacts of natural disasters. A fundamental component of this process is the implementation of coastal morphodynamics monitoring through Structure-from-Motion (SfM) techniques, utilizing high-resolution 2D/3D data obtained from aerial photogrammetry. To assess morphodynamic changes over a three-year period (2022 – 2024), several SfM-based photogrammetric studies were conducted, each year, in the Romanian sector of the Danube-Black Sea coastal zone, specifically at the mouth of one of the Danube River distributaries (Sf Gheorghe branch) into the Black Sea, and along the left bank, near Sf Gheorghe locality, located within the Danube Delta Biosphere Reserve (DDBR). The essential equipment for aerial photogrammetry comprises Unmanned Aerial Vehicles (UAVs) and Global Navigation Satellite Systems (GNSS). In this study, the UAV used was a DJI Mavic 3T (Enterprise/Thermal) drone, complemented by two Trimble R12i and R4 GNSS systems, as well as approximately 10 Ground Control Points (GCPs). Data acquisition and processing were carried out using specialized photogrammetric software (Agisoft Metashape) along with various GIS tools (e.g., Blue Marble Geographics Global Mapper and ESRI ArcMap). The photogrammetric products generated for the study, as detailed in this paper, include Digital Elevation Models (DEMs), Digital Terrain Models (DTMs), orthomosaics (orthophotos), and others. At Sfântu Gheorghe beach, a comparison between 2023 and 2024 photogrammetric surveys revealed that the left bank of the Sf. Gheorghe Arm, at the river mouth into the Black Sea, suffered from a twist (erosion) of up to 64 metres. Additionally, on the selected perimetre (total area of 31,910 square meters ) from the beach and dune zone of Sf. Gheorghe, an area of up to 16,202 square meters was eroded between 2023 and 2024. This contrasts with the period between 2022 and 2023, during which deposition predominated. Erosion at the Danube mouths and the adjacent Black Sea coastline is driven by a complex interaction of natural and anthropogenic factors. Natural processes, including subsidence, sea-level rise, and episodic extreme storm events, contribute significantly to coastal dynamics. Meanwhile, human-induced factors, such as upstream hydrotechnical works that limits sediment transport, cutting of navigation canals, as well as the exacerbating effects of climate change, further accelerate erosion. The recent Structure-from-Motion (SfM) surveys provide essential quantitative data, enabling a detailed analysis of both short-term and long-term morphodynamic changes influenced by seasonal variations and extreme hydrometeorological events in this highly dynamic coastal system.

Sfȃntu Gheorghe Mouth (Danube Delta)↗
USGS2023

Sediment thickness map of United States Atlantic and Gulf Coastal Plain Strata, and their influence on earthquake ground motions

With the recent successful accounting of basin depth ground-motion adjustments in seismic hazard analyses for select areas of the western United States, we move toward implementing similar adjustments in the Atlantic and Gulf Coastal Plains by constructing a sediment thickness model and evaluating multiple relevant site amplification models for central and eastern United States seismic hazard analyses. We digitize and combine existing sediment thickness data sets into a composite surface that delineates the base of Cretaceous sediments under the Atlantic Coastal Plain and the base of Mesozoic sediments under the Gulf Coastal Plain. Amplification models dependent on sediment thickness, site natural period, and source-to-site path length are compared with data sets of observed ground motions to evaluate the ability of the new models to improve ground motion estimates. We find that the amplification models can account for observed trends in sediment-thickness and period-dependent residuals, but some tuning is required. For example, the model of Chapman and Guo requires a reference V S 30 , the time-averaged shear-wave velocity within 30 m of the Earth’s surface, for non-Coastal Plain sites, which we estimate to be between about 1 and 2 km/s. Along with our sediment thickness model, we estimate a velocity profile for application to the Harmon et al. site-natural-period-based model in order to best match the Chapman and Guo period dependence for a broad range of sediment thicknesses. The Next Generation of Attenuation models for the eastern United States Gulf Coast path-based adjustment models can also account for seismic attenuation in the Coastal Plain sediments and reduce the standard deviation of total residuals. If enacted in the U.S. Geological Survey National Seismic Hazard Model, these amplification models will reduce predicted short-period (<1 s) and increase predicted long-period (>1 s) ground motions in the Coastal Plains appreciably.

Earthquake Spectra↗
USGS2025

Regional characterization of coal resources in the U.S. Gulf Coast

There is increasing interest in extracting critical minerals (CM), including rare earth elements (REE), from coals in the United States to address the overreliance on imported REE. The U.S. Gulf Coast and the Williston basins are the two major lignite-bearing basins within the country. Recent REE and CM studies of the lignite in these basins have indicated that the coals may be a viable source material for REE and CM extraction. To evaluate in-place coal as a potential source of REE and CM, the coal resources need to be quantified. This study presents the results of a regional analysis of the U.S. Gulf Coast lignite and bituminous coal resources that might be available as potential sources of REE and CM. The resource analysis used kriging methods to develop isopleth maps of cumulative coal thickness throughout the region using data from 31,181 drill holes and other data points. The estimated total coal resource in the Gulf Coast is about 83 billion metric tons in the upper 90 m (~ 300 ft) of the subsurface. Texas accounted for 40 percent (32 billion metric tons) of the total resource, followed by Mississippi (24 %, 20 billion metric tons), Louisiana (14 %, 12 billion metric tons), Tennessee (10 %, 8.5 billion metric tons), and Arkansas (6 %, 5.1 billion metric tons). The remaining states each accounted for less than 5 percent of the total resource. Georgia had the smallest resource estimated at 7 million metric tons. Here we report the first known state-wide lignite resource estimates for Georgia, Kentucky (820 million metric tons), and Missouri (1,800 million metric tons). A comparison of the results of this study with those of previous Gulf Coast and Williston Basin resource studies is difficult because each study used different data sources, assessment methodologies, overburden depths, and qualifying coal thicknesses. Coal-power electric generation has sharply decreased in past decades and mining of these coals for CM and REE could provide additional co-products such as activated carbon and other uses such as fertilizer (soil enhancer).

EarthArXiv↗
USGS2025

Preserving and increasing water resources—Natural infrastructure in dryland streams in Baja California Sur, Mexico

The Los Planes watershed of Baja California Sur, Mexico, and its underlying aquifer are experiencing groundwater decline owing to low average annual rainfall (28.1 centimeters per year) and rising water demand from population growth and agricultural activities. This decline in water availability can lead to desertification—a process that changes arable land to desert by degrading soil and vegetation—and can pose serious challenges to livelihoods that depend on the land. To address these issues, a ranch in the Los Planes watershed has installed many natural infrastructures in dryland streams (NIDS) in channels for soil and water conservation. In 2022, the U.S. Geological Survey (USGS) began working with regional researchers and land managers to investigate the effects of NIDS on natural biological, geochemical, and physical processes and determine the efficacy of NIDS for water augmentation in the Los Planes watershed. The USGS also worked with local academic institutions and nonprofit organizations to create public educational opportunities focused on the area’s hydrogeology. These and other collaborative efforts with the U.S. Water Partnership and Innovaciones Alumbra aim at enhancing water resources in the Baja California Sur region and promoting water security and safeguarding community well-being. Resume La cuenca de Los Planes, ubicada en Baja California Sur, México, y su acuífero subyacente, están sufriendo una disminución de las aguas subterráneas debido a la baja precipitación media anual (28.1 centímetros por año) y la alta demanda de agua por parte de una población creciente y la actividad agrícola. Esta disminución de la disponibilidad de agua puede conducir a la desertificación—un proceso que por medio de la degradación del suelo y la vegetación convierte a la tierra cultivable en desierto—representando un serio desafío para los medios de vida de las personas. Para abordar estos problemas, un rancho en la cuenca de Los Planes ha instalado numerosas obras de Infraestructura Natural en Arroyos de Tierras Áridas (INATS) para conservación del suelo y del agua. En 2022, el Servicio Geológico de los Estados Unidos (USGS, por sus siglas en inglés) comenzó a trabajar con investigadores regionales y gestores de tierras para estudiar los efectos de INATS en los procesos biológicos, geoquímicos y físicos, y determinar su eficacia en el aumento de los recursos hídricos en la cuenca de Los Planes. El USGS se ha asociado con instituciones académicas y organizaciones locales sin fines de lucro para crear oportunidades educativas públicas centradas en la hidrogeología de la zona. Estos y otros esfuerzos colaborativos con la Asociación del Agua de Estados Unidos (U.S. Water Partnership) e Innovaciones Alumbra, tienen como objetivo mejorar el uso de los recursos hídricos en la región de Baja California Sur, promover la seguridad hídrica y proteger el bienestar de la comunidad.

Baja California Sur↗
USGS2026

Insights into Mountain Pass carbonatite formation from in-situ sulfur isotopes and geochemistry of sulfate and sulfide minerals

The Mountain Pass carbonatite stock hosts a world-class rare earth element deposit and may be classified as a carbonate-sulfate igneous rock, as it contains on average > 50 volume percent carbonate minerals and 20 to 30 volume percent sulfate minerals. The sulfates range in composition from barite to celestine and locally occur with sparse sulfide minerals. We investigate the origin of sulfur enrichment and the occurrence of sulfur-bearing minerals in the Mountain Pass carbonatite with in-situ sulfur isotope and mineral chemistry. Barite cores with δ 34 S of 1 to 3‰ do not coexist with sulfides, whereas celestine rims with δ 34 S of > 3‰ are associated with sulfides with δ 34 S < -10‰. We propose a model in which sulfur-bearing sediments were subducted during episodes of plate convergence in the Mojave Province that preceded Mountain Pass magmatism. Metasomatism of the overlying mantle by melts derived from the subducted sediments generated an unusually carbon- and sulfur-rich source to yield carbonatite magmas. Sulfur from primary carbonatite magmas and ~ 1 to 7% sulfur from subducted sediment melts yielded a slightly enriched δ 34 S composition (relative to depleted mantle δ 34 S of -1‰) for early crystallizing barite. Celestine rims on magmatic barite cores formed at low, hydrothermal temperatures (< 350 °C) based on S isotope thermometry for equilibrium celestine-galena and celestine-pyrite pairs. The sparse sulfides in the carbonatite stock are not in equilibrium with the primary barite cores and therefore do not permit S isotope thermometry estimates of magmatic temperatures. The S/Se ratios of sulfide minerals (> 3,400) typically exceed primitive mantle values (S/Se of 3,340), also consistent with their derivation from hydrothermal fluids. Trace occurrences of sulfide and sulfate minerals in alkaline silicate stocks related to the carbonatite stock have similar δ 34 S compositions and yield similarly low formation temperatures, suggesting regionally extensive and chemically similar sulfur-bearing hydrothermal fluids that imparted lithologically diverse rocks with a consistent sulfur isotope fingerprint.

California, Nevada↗
USGS2025

Pre-Acadian tectonics of the eastern Orange-Milford Belt, south-central Connecticut

This excursion presents a reinterpretation of mapping and new analytical data from the eastern Orange-Milford belt (OMB) in south-central Connecticut. The OMB is a fault-bound terrane of argillites and mafic rocks of anomalously low metamorphic grade—and of poorly constrained ages and tectonic affinity—wedged between kyanite/sillimanite-grade peri- Laurentian rocks to the west and anatectic peri-Gondwanan rocks to the east (Fig. 1A). Our data demonstrate that Ordovician(?) igneous and sedimentary rocks of the OMB were variably metamorphosed in the Ordovician and Silurian but escaped regionally pervasive, high-grade Devonian and later metamorphism. Previous interpretations (Fritts 1963a, 1965a, 1965b; Burger, 1967; Burger and others, 1968; Rodgers, 1985) described these rocks as a conformable sequence of low-grade, Ordovician to Devonian metasediments and metavolcanics. Our results reveal that the “metavolcanics” are not extrusive rocks but rather slivers of lower oceanic crust with complicated high- and low-grade metamorphic fabrics, intruded by a swarm of Silurian sheeted basalt dikes, and in fault contact with the surrounding metasediments. The purpose of this trip is to show evidence of early Paleozoic (Taconic) deformation and metamorphism preserved in rocks of the eastern OMB. These rocks remained shallow, cool, and sufficiently dry during the regionally dominant Acadian and Alleghanian orogenies to have avoided significant overprinting. As such, these rocks serve as windows into a geologic history otherwise unavailable between anatectic rocks of the peri-Gondwanan Bronson Hill, Avalon, and Gander terranes east of the Hartford basin and sillimanite-grade rocks of the peri-Laurentian Hartland and gneiss dome belts west of the OMB. We present major and trace element geochemistry including rare-earth element patterns for all mafic units in the eastern OMB as well as 40 Ar/ 39 Ar age spectra of amphibole, muscovite, and K-feldspar from rocks of the Maltby Lakes complex (of Deasy and others, 2017), Savin Schist, and Wepawaug Schist. Our evidence demonstrates that the units of the OMB have been assembled by faulting or intrusion, and that no stratigraphic relationships exist between the argillaceous schists and the metaigneous rocks.

Connecticut↗
USGS2023

Oligocene–Miocene northward growth of the Tibetan Plateau: Insights from intermontane basins in the West Qinling Belt, NW China

Growth of the Tibetan Plateau, Earth’s broadest and highest elevation collisional system, shapes orographic barriers, reorganizes drainage networks, and influences surface erosion and sediment delivery, whose changes in space and provenance feed back to intracontinental tectonic processes. Studies of interior basins within the northern Tibetan Plateau provide new sediment accumulation, provenance, paleodrainage, and deformation timing data that enable a reconstruction of the far-field tectono-geomorphic evolution of the rising Tibetan Plateau. Along the northern plateau margin, topographic growth in the West Qinling Belt is inferred to have initiated in the Eocene, nearly coeval with the India-Asia collision, as well as in the late Miocene. However, geological knowledge about the intervening period remains at present enigmatic, and the kinematics and dynamics are uncertain. This study presents a multidisciplinary data set from the intermontane Anhua-Huicheng Basin (AHB; Gansu Province, China) to fill this gap. Magnetostratigraphic dating, regional mapping, and sedimentological analysis imply that contractional deformation and thrust-top basin systems formed within the West Qinling Belt in the Oligocene (not later than ca. 24 Ma). A combination of observations including paleocurrent changes, detrital zircon U-Pb age variations, and appearance of growth strata along the Anhua-Huicheng Basin reveal the rapid uplift of the West Qinling Belt at ca. 15 Ma. Sedimentation in the intermontane basins ended after the late Miocene (ca. 8 Ma), when the region experienced intrabasinal deformation, uplift, and erosion with the establishment of an external drainage system. Since the late Miocene, the growth of the West Qinling Belt reached a climax with the lack of substantial contractional deformation in Cenozoic sequences heralding the onset of the modern kinematic regime and attainment of high elevation. Observed transitions in the tectonostratigraphy and paleodrainage define different phases of deformation and plateau-wide shifts in stress reorganization, which led to the northward growth and later lateral expansion of the Tibetan Plateau.

GSA Bulletin↗
USGS2026

Ten new insights in climate science 2025

Interdisciplinary understanding is vital for delivering sound climate policy advice. However, navigating the ever-growing and increasingly diverse scholarly literature on climate change is challenging for any individual researcher. This annual synthesis highlights and explains recent advances across a variety of fields of climate change research. This year, the 10 insights focus on: (1) the record-warmth of 2023/2024 and the elevated Earth energy imbalance; (2) acceleration of ocean warming and intensifying marine heatwaves; (3) northern land carbon sinks under strain; (4) reinforcing feedback between biodiversity loss and climate change; (5) accelerated depletion of groundwater; (6) global dengue incidence; (7) global income losses and labour productivity declines; (8) strategic scaling of CDR; (9) integrity challenges in carbon credit markets and emerging responses; and (10) effective policy mixes for emissions reductions. The insights have been written to be accessible to researchers from different fields, serving as entry-points to specific topics, as well as providing an overview of the evolving landscape of climate change research. In the final section, the insights are used to develop overarching policy-relevant messages. This paper provides the basis for a science-policy report that was shared with all Party delegations ahead of COP30 in Belém, Brazil.

Global Sustainability↗
USGS1999

Simulation of the effects of development of the ground-water flow system of Long Island, New York

Extensive development on Long Island since the late 19th century and projections of increased urbanization and ground-water use makes effective water-resource management essential for preservation of the island's hydrologic environment and maintenance of a reliable source of water supply. This report presents results of a ground-water flow simulation analysis of the effects of development on the Long Island ground-water system. It describes ground-water levels, stream-flow, and the ground-water budget for the predevelopment period (pre-1900), the 1960's drought, and a more recent (1968-83) period with significant hydrologic stress. The report also presents estimated effects of a proposed water-supply strategy for the year 2020. Long Island has three major aquifers-the upper glacial (water-table), the Magothy, and the Lloyd aquifers-that are separated to varying degrees by confining units. Before development, recharge from precipitation entered the ground-water system at a rate of more than 1.1 billion gallons per day. An equal amount discharged to streams (41 percent), the shore (52 percent), and subsea boundaries (7 percent) . Urbanization and withdrawal of more than 400 Mgal/d (million gallons per day) from wells have resulted in local effects that include declines in ground-water levels, drying up and burial of streams and wetlands, reduction of ground-water recharge by increased overland flow to the ocean, a general decrease in ground-water discharge, and salt water intrusion. In some areas, the reduction in recharge is mitigated by leakage from water-supply and wastewater disposal lines, and infiltration of storm water through recharge basins. During 1968-83, a net loss of 240 Mgal/d from the ground-water system caused a decrease in ground-water discharge to streams (135 Mgal/d), to the shore (82 Mgal/d), and to subsea boundaries (23Mgal/d).The greatest adverse effects have been in western Long Island, where the most severe development has occurred. This analysis shows stream base flow to be highly sensitive to water-table fluctuations, and long streams to be more sensitive than short ones. A water-supply scenario for the year 2020 was simulated that employs redistribution of pumping centers to mitigate extreme local effects . Although the net stress on the ground-water system was projected to increase 57 Mgal/d (24 percent) above that of 1968-83, redistribution of ground-water withdrawals across the island would allow recovery of cones of depression in western Long Island, thereby reducing the threat of salt water intrusion and increasing base flow of some streams . The increased stress would cause a net decrease in base flow island wide of 44 Mgal/d; total base flow would be 281 Mgal/d - 39 percent below predevelopment levels or 14 percent below 1968-83 levels. The most severe effects would be in Nassau and western Suffolk Counties.

New York↗