USGS Science⌕ Search

SEARCH · USGS Science

Results for “Journal of African Earth Sciences”

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

Flood-inundation maps for 14.8 miles of Little and Big Papillion Creeks in Omaha, Nebraska, 2023

Digital flood-inundation map libraries for two reaches that constitute 14.8 miles of Little and Big Papillion Creeks in Omaha, Nebraska, were created by the U.S. Geological Survey (USGS) in cooperation with the Papio-Missouri River Natural Resource District. The flood-inundation maps, which can be accessed through the USGS Flood Inundation Mapping Program website at https://www.usgs.gov/mission-areas/water-resources/science/flood-inundation-mapping-fim-program , depict estimates of the areal extent and depth of flooding corresponding to selected water levels (stages) at Little Papillion Creek at Irvington, Nebr. (USGS station 06610750), Little Papillion Creek at Ak-Sar-Ben at Omaha, Nebr. (USGS station 06610765), and Big Papillion Creek at Q Street at Omaha, Nebr. (USGS station 06610770) streamgages. Near-real-time stages at these streamgages may be obtained from the USGS National Water Information System database at https://doi.org/10.5066/F7P55KJN or from the National Weather Service Advanced Hydrologic Prediction Service at https://water.weather.gov/ahps/ . Flood profiles were computed for two different reaches that constitute 14.8 miles of stream length in the study area by using hydraulic models. The models were calibrated by adjusting roughness coefficients to best represent the current (2022) stage-streamflow relation at the streamgages within the study reach. The hydraulic models were then used to compute water-surface profiles at 1-foot stage intervals for selected stage ranges to represent various flooding scenarios at the streamgages in each reach. The simulated water-surface profiles then were combined with a digital elevation model using a geographic information system, which had a 10-foot grid spacing to delineate the flooding extents and water depths for each stage. The availability of these flood-inundation maps, along with information regarding current stage from the USGS streamgages, can provide emergency management personnel and residents with information that is critical for flood response activities and post flood recovery efforts.

Nebraska↗
USGS2026

Regreening, restoring, and reconnecting a southwestern wetland ecosystem – the Zeedyk wetland

Alluvial wetland ecosystems are vital as biodiversity hotspots but are increasingly threatened by anthropogenic stressors and drought. These pressures are especially acute in arid and semi-arid regions, where eco-hydrologic connectivity is fragile and recovery is slow. This study quantifies the efficacy of nature-based solutions, particularly the ‘Zeedyk approach,’ which employs low-tech Natural Infrastructure in Dryland Streams (NIDS)—including rock detention structures—to slow surface water, raise groundwater tables, and restore wetland function at a spring-fed wetland in Cebolla Canyon, New Mexico, U.S.A. Our results depict a Restoration Feedback Loop that captures stages of change from a healthy wetland in 1935, altered by 20th-century agriculture and grazing, to the re-establishment of the historical flow regime by 2024 documented through an 89-year archive of aerial imagery (1935–2024). By the end of our study period, the Spring-Fed Wetland had expanded by roughly 229% of the original 1935 area, to 4.13 ha. Using 40 years of satellite data, we assess changes in vegetation and hydrology with remote sensing indices. Spatial and temporal analyses reveal significant increases in vegetation greenness and wetness, particularly in an Expanded Wetland subregion, which exhibited ∼3.5x higher wetness and ∼1.5x higher greenness trends compared to adjacent areas. Monthly metrics highlight seasonal variability, with increases in greenness linked to monsoonal rainfall and lateral water redistribution, indicating that restoration impacts extend beyond the primary wetland. This study demonstrates the utility of cloud-based platforms like Google Earth Engine and USGS EarthExplorer for long-term monitoring of wetland restoration, while quantifying the efficacy of the ‘Zeedyk approach’ and demonstrating its potential as a scalable method to restore and conserve wetland meadows in other arid and semi-arid landscapes.

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↗
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↗
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

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

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↗
USGS2024

Plan to coordinate post-earthquake investigations supported by the National Earthquake Hazards Reduction Program (NEHRP)

Introduction This report presents a plan supported by the National Earthquake Hazards Reduction Program (NEHRP) to coordinate domestic and international post-earthquake investigations (herein called “the Plan”). Post-earthquake scientific and engineering investigations are undertaken to capture critical information to understand the causes and impacts of the event, lessons from which can substantially improve the Nation’s resilience after future earthquakes. NEHRP is the Federal Government’s coordinated nationwide program to reduce risks to life and property from earthquakes. The Plan describes the activation and coordination of the four designated NEHRP Agencies in the Federal Government: Federal Emergency Management Agency, National Institute of Standards and Technology, U.S. National Science Foundation, and U.S. Geological Survey (USGS). The Plan also describes coordination between NEHRP Agencies and other organizations that may participate in pre-event and post-earthquake investigations, including non-NEHRP Federal agencies; State, regional, local, Tribal, and territorial agencies; domestic nongovernmental organizations; academic institutions and affiliated organizations; private companies; foreign governmental agencies and nongovernmental organizations; and international organizations. The Plan delineates the coordination of NEHRP post-earthquake scientific and engineering investigations to document the direct, indirect, and cascading physical and societal impacts from fault rupture and ground shaking hazards and from secondary hazards such as landslides, liquefaction, and tsunamis. In addition, the Plan identifies pre-event activities necessary to ensure that post-earthquake investigations are executed effectively. The USGS is the lead NEHRP Agency for activating and coordinating NEHRP post-earthquake investigations and for implementing this Plan. The USGS also leads coordination of the NEHRP Agencies in completing the pre-event activities identified in the Plan. The Plan has new information and supersedes USGS Circular 1242, “The Plan to Coordinate NEHRP Post-Earthquake Investigations,” which was published in 2003. This second-generation Plan was developed with the assistance of the Applied Technology Council (ATC) of Redwood City, California, under USGS contract 140G0121P0309, ATC-155. A 17-member Project Review Panel provided guidance on plan development, and input was solicited from subject matter experts representing key stakeholder groups and from participants at a public workshop.

Circular↗
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↗
USGS2026

Mountain goat declines in a protected, interior, native population

A shifting climate poses threats to alpine-adapted species including mountain goats. We used long-term (12 years) citizen science monitoring data and Bayesian N-mixture modeling to estimate population trends and drivers of population metrics among mountain goats in Glacier National Park (GNP). Median goats per site ( n = 37 sites) declined by 45% (95% credible interval [CRI] = 32%, 57%) from 77.8 (95% CRI = 64.4, 95.1) in 2008 to 42.3 (95% CRI = 34.3, 52.2) in 2019, with consistent declines from 2008 until 2015, when the number of estimated goats stabilized. The decline exceeds IUCN criteria for classifying a population as vulnerable, >30% declines over only two generations. Across years, relatively few goats occupied northwestern GNP. Goat numbers declined the most at northeastern sites, trended toward decline in most southern sites, and increased at only two west-central sites. The proportion of permanent snow and glaciers, the presence of natural mineral licks, and habituation strongly increased the initial abundance of goats in the area. Weather variables had the greatest influence on population growth rates, particularly precipitation between May 15 and June 15 of the previous summer, the neonatal period. Lower growth occurred with less snow water equivalent and lower mean winter temperature, early summer temperature, and early summer precipitation. Projected reductions of permanent snow, increasing spring and summer temperatures, and insufficient and variable spring precipitation raise concerns for the future of native goats in this region. Our analyses reveal ways to improve detection rates of goats during surveys, which is important for optimizing the precision of estimates and the power to detect future trends. Detection increased with goat habituation, retention of observers with experience, use of binoculars, and conducting surveys at lower temperatures and earlier dates. Improving detection will be particularly important given the lower number of goats currently observed in the park. Research to estimate park-wide population size, evaluate genetic structure and diversity, assess changing habitat, human recreation levels and forage, and forward-project climate effects on persistence will be crucial to understanding the context of these results and conserving this iconic, metapopulation at the southern edge of the distribution of native mountain goats.

Montana↗
USGS2025

Overview of The SCEC/USGS Community Stress Drop Validation Study using the 2019 Ridgecrest earthquake sequence

We present initial findings from the ongoing Community Stress Drop Validation Study to compare spectral stress‐drop estimates for earthquakes in the 2019 Ridgecrest, California, sequence. This study uses a unified dataset to independently estimate earthquake source parameters through various methods. Stress drop, which denotes the change in average shear stress along a fault during earthquake rupture, is a critical parameter in earthquake science, impacting ground motion, rupture simulation, and source physics. Spectral stress drop is commonly derived by fitting the amplitude‐spectrum shape, but estimates can vary substantially across studies for individual earthquakes. Sponsored jointly by the U.S. Geological Survey and the Statewide (previously, Southern) California Earthquake Center our community study aims to elucidate sources of variability and uncertainty in earthquake spectral stress‐drop estimates through quantitative comparison of submitted results from independent analyses. The dataset includes nearly 13,000 earthquakes ranging from M 1 to 7 during a two‐week period of the 2019 Ridgecrest sequence, recorded within a 1° radius. In this article, we report on 56 unique submissions received from 20 different groups, detailing spectral corner frequencies (or source durations), moment magnitudes, and estimated spectral stress drops. Methods employed encompass spectral ratio analysis, spectral decomposition and inversion, finite‐fault modeling, ground‐motion‐based approaches, and combined methods. Initial analysis reveals significant scatter across submitted spectral stress drops spanning over six orders of magnitude. However, we can identify between‐method trends and offsets within the data to mitigate this variability. Averaging submissions for a prioritized subset of 56 events shows reduced variability of spectral stress drop, indicating overall consistency in recovered spectral stress‐drop values.

California↗
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↗
USGS2025

Geochemistry and radiogenic isotopes constrain the mantle source region of the Mountain Pass Intrusive Suite, California

The Mountain Pass carbonatite stock is the largest rare earth element (REE) deposit and only active REE mine in the United States. The carbonatite intrusion and spatially associated alkaline silicate intrusions constitute the Mountain Pass Intrusive Suite, which is located within the Mojave Province in California. Both the carbonatite and the alkaline silicate rocks are enriched in large ion lithophile elements and light REEs, and less enriched to depleted in high field strength elements, indicating the mantle source region was metasomatically enriched in incompatible trace elements. The cause of this metasomatic mantle enrichment and the genetic relationship between the carbonatite and the alkaline silicate stocks are poorly understood. In this study, major and trace element geochemical data and isotopic (Rb-Sr, Sm-Nd, and Lu-Hf) data are presented to constrain genesis of the Mountain Pass Intrusive Suite, from mantle source region to the intrusion of the stocks. Our geochemical data are consistent with derivation of the alkaline silicate and carbonatite melts through partial melting from a shared mantle source region rather than through liquid immiscibility or fractional crystallization and separation of a carbothermal fluid. Although the Rb-Sr isotopic system in the Mountain Pass Intrusive Suite is disturbed at the whole-rock scale, the isotopic systems for whole-rock Sm-Nd (εNd i = ‐2.2 ± 0.8) and zircon Lu-Hf (εHf i = 0.1 ± 1.1) are robust and support mantle derivation of the magmas. Geochemical modeling using experimentally derived partition coefficients was used to identify possible causes of enrichment in incompatible elements through metasomatism in the mantle source region. Modeling of metasomatism by melts derived by partial melting of deeply subducted carbonated sediments approximates observed Mountain Pass Intrusive Suite trace element chemistry. Scattered εHf i in inherited zircon (2.8 ± 2.6) is consistent with derivation from an arc-related environment with substantial crustal contamination. Paleotectonic studies in the Mojave Province indicate that regional subduction preceded emplacement of the Mountain Pass Intrusive Suite by ∼300 Ma. Melting of the Mountain Pass source region may have been caused by post-collisional thermal relaxation and extension.

California↗
USGS2025

Groundwater-storage change in the north Phoenix aquifer, Arizona, 2020–23

The city of Phoenix, Arizona, relies primarily on surface water for municipal water supply. The city also maintains wells to withdraw groundwater, particularly in times of drought and reduced surface-water supply, and to recharge groundwater when excess surface water is available. As of 2023, withdrawals from the aquifer in the northeastern part of the city are a small volume of water, less than 3,000 acre-feet in most years. Each year a similar volume of water is recharged through injection wells. The withdrawal wells are permitted to produce in total more than 20,000 acre-feet per year; increased pumping could lead to future groundwater-storage declines. To better understand groundwater-storage change in the north Phoenix aquifer, a repeat microgravity monitoring network was established in 2020. Measurements of changes in Earth’s gravity provide a direct, non-invasive measurement of subsurface mass change. Groundwater-storage changes were small during the 2020–23 study period, consistent with the relatively small volumes of pumping and recharge. Groundwater levels measured in monitoring wells were stable or increased slightly during this period, although the number of monitoring wells within the area of the gravity monitoring network is sparse. In total, about 15,000 acre-feet of water were pumped and 31,000 acre-feet recharged through injection in the north Phoenix aquifer during the 2020–22 period. Within the monitored area, groundwater storage increased by a small amount, about 1,000 acre-feet, in 2020, and decreased by a larger amount, about 6,000 acre-feet, each year in 2021 and 2022. Storage decreased at 89 of 102 stations from 2020 to 2023. Groundwater-storage decreases greater than the volume of net pumping indicate down-gradient subflow from the aquifer is greater than recharge plus incoming subflow, drying of the unsaturated zone resulting from decreased land-surface recharge, or both. At present (2023), the aquifer appears able to store and supply the relatively small amounts of water needed without excessive drawdown or harmful effects, such as land subsidence and surface inundation. If pumping or recharge increases in the future, the established repeat microgravity network is well suited to capture the spatial extent and magnitude of groundwater-storage changes in the aquifer.

Arizona↗
USGS2000

Quality assurance and analysis of water levels in wells on Pahute Mesa and vicinity, Nevada Test Site, Nye County, Nevada

Periodic and continual water-level data from 1963 to 1998 were compiled and quality assured for 65 observation wells on Pahute Mesa and vicinity, Nye County, Nevada. As part of the quality assurance of all water levels, ancillary data pertinent to computing hydraulic heads in wells were compiled and analyzed. Quality-assured water levels that were not necessarily in error but which did not represent static heads in the regional aquifer system, or required some other qualification, were flagged. Water levels flagged include those recovering from recent pumping or well construction, water levels affected by nuclear tests, and measurements affected by borehole deviations. A cursory examination of about 30 wells with available water-level and down-hole temperature data indicate that water levels in most wells on Pahute Mesa would not be significantly affected by temperature if corrected to 95 degrees Fahrenheit. Wells with large corrections (greater than 10 feet) are those with long water columns (greater than 1,500 feet of water above the assumed point of inflow) in combination with mean water-column temperatures exceeding 105 degrees Fahrenheit. Water-level fluctuations in wells on Pahute Mesa are caused by several factors including infiltration of precipitation, barometric pressure, Earth tides, ground-water pumpage, and seismic events caused by tectonic activity and underground nuclear testing. No observed water-level fluctuations were attributed to a naturally occurring earthquake. The magnitude and duration of changes in water levels caused by nuclear tests are affected by the test size and the distance from a well to the test. Identifying water levels that might be affected by past nuclear tests is difficult because pre-testing water-level data are sparse. Hydrologically significant trends were found in 13 of 25 wells with multiple years of water-level record. The largest change in water levels (1,029 feet in 25 years) occurred in well U-19v PS 1D as a result of the Almendro nuclear test. Likely explanations for trends in most of the wells are either changes in precipitation patterns that affect recharge rates to the ground-water system, pumping effects from water-supply well U-20 WW, or a combination of these two factors.

Nevada↗