Geology Reports⌕ Search

SEARCH · Geology Reports

Results for “Journal of Wildlife Management”

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

Non-destructive sampling of larval amphibians shows tail tissues reflect whole-body methylmercury concentrations and trophic-related differences in bioaccumulation

Mercury (Hg), and specifically methylmercury (MeHg), is a contaminant of global concern to humans and wildlife, but there is still limited understanding of its effects on many taxa. Amphibians are closely associated with water, where Hg is converted to MeHg, and high concentrations of MeHg can reduce survival of amphibians in the wild. However, limited non-lethal proxies exist for estimating MeHg bioaccumulation in amphibian larvae, which could be used to reduce destructive sampling and allow sampling of imperiled species. We evaluated whether tail clips were indicative of whole-body MeHg concentrations of larval amphibians, whether dragonflies were an effective bioindicator of Hg concentrations of larval amphibians, and compared MeHg bioaccumulation in predatory caudate (salamanders, newts) and non-predatory anuran (frogs, toads) larvae. Tail-clip and whole-body MeHg concentrations were strongly correlated, especially for caudates. There was high variability among sites and species, but caudates had higher MeHg than anurans, and MeHg concentrations of caudates and anurans were strongly correlated. Dragonflies were a weak indicator of MeHg for anurans but strong indicators for larval caudates, likely because caudates and dragonflies are carnivorous and occupy similar trophic positions. Our study revealed that sampling tails can be an effective non-destructive index of whole-body MeHg for larval amphibians and demonstrates differences in MeHg bioaccumulation related to trophic position in larval amphibians and dragonflies. All post-embryo life stages for a broad suite of amphibian species can now be sampled non-lethally for MeHg, allowing for a broader examination of population-level effects and focus on imperiled species that cannot be sampled destructively.

conterminous United States↗
USGS2026

Groundwater drought in the United States: Spatial and temporal variability

Many communities and ecosystems in the United States that are dependent on groundwater are potentially adversely affected by groundwater drought. We computed yearly groundwater-drought metrics and mean groundwater levels at well locations across the conterminous United States (CONUS), using data from wells and remotely sensed and modeled Gravity Recovery and Climate Experiment Drought Monitor Data Assimilation (GRACE-DADM). We also modeled the probability of low or high human impact at each well location. The spatial distribution of groundwater-drought duration and severity from 2001 to 2020 for 1,510 wells shows longer maximum duration and higher maximum severity events in drier regions like the Southwest than in wetter regions like the Northeast. Based on 613 wells in CONUS from 1981 to 2020, there are many significant decreases in drought duration and severity in the Northeast and many significant increases in annual-mean groundwater levels. In contrast, there are many significant increases in drought metrics and decreases in mean water levels in parts of the Southeast. There are major differences in trends from 2001 to 2020 between well-based and GRACE-DADM-based groundwater metrics in some CONUS regions and a very low correlation between trends at individual locations across CONUS. A potential reason for this disparity is the low GRACE-DADM resolution (∼12 km) and the potential for a large amount of groundwater variation at the local scale. Also, GRACE-DADM represents shallow, unconfined aquifers which may not match the screened interval of the monitoring wells we evaluated. Large spatial gaps in long-term, high frequency, and quality-assured groundwater-well monitoring data present a challenge for understanding groundwater-drought variability across CONUS. Remote sensing tools such as GRACE can help but cannot fully replace well monitoring, as highlighted by our study results. Substantially more long-term monitoring wells would more accurately represent groundwater-drought trends and spatial variability across CONUS, particularly in western regions.

conterminous United States↗
USGS2025

Time-varying rates of organic and inorganic mass accumulation in southeast Louisiana marshes: Relationships to sea-level anomalies and tropical storms

Louisiana's coastal wetlands are complex systems that require a continuous input of organic and inorganic material to keep pace with relative sea-level rise. Coastal restoration projects such as sediment diversions are being implemented to mitigate land loss and increase availability of inorganic sediment to coastal wetlands, and marshes specifically rely on organic material to build soil volume and maintain surface elevation. Interannual-to-decadal sea-level anomalies such as hurricanes can affect marsh accretion, mineral deposition, and plant productivity. In this light, complex ecogeomorphic feedback controls whether a marsh surface is sustainable or eroded/drowns. This study performs some of the first differential vertical accretion rates (VARs) and organic and inorganic mass accumulation rates (MARs) over time in SE Louisiana marshes determined from the 210 Pb Constant Rate of Supply model, coupled with standard 137 Cs VARs. These accumulation rates over the past ∼100 years were measured from a total of six brackish and salt marsh locations in Barataria Basin near the proposed Mid-Barataria Sediment Diversion. They were then related to interannual sea-surface elevations at Grand Isle, Louisiana, over the last ∼60 years and recorded hurricane activity in the delta. Results show VARs range from 0.63 cm/y to 1.69 cm/y and total MARs range from 0.11 to 0.43 g/cm 2 /y. Temporally, VARs and MARs (total, inorganic, and organic) are characterized by gradual increases in rates with decreasing age along with episodic peaks in VARs and MARs. The findings of this study indicate that no relationship occurs between sea-level anomalies and VARs or organic and inorganic MARs; however, a strong relationship appears to occur between major hurricanes to VAR and MAR contributions. Furthermore, high water content (81 ± 8%) and organic-rich soils in the sediment cores highlight the significance of belowground biomass and associated pore volume in maintaining marsh elevation in the study area.

Louisiana↗
USGS2026

Groundwater quality near an oil field in a stream-dominated recharge setting, California, USA

Alluvial valley aquifers are important sources of water supply in many areas but effects of co-located oil and gas development on these resources have not been widely reported, especially in settings where recharge is dominated by stream infiltration. Interpreting the presence of geochemical indicators in the context of hydrology, geology, and other factors provides a more complete understanding of the relations between groundwater and sources of oil-field fluids and aids in identifying risks associated with oil and gas development. Groundwater and Salinas River water samples were collected in an alluvial valley near the San Ardo Oil Field in Monterey County, California and analyzed for a wide range of dissolved chemical, gas, and isotopic constituents to determine if oil-field fluids (water and gas from oil-producing and non-producing zones) have mixed with fresh groundwater used for supply. Hydraulic gradients, age-dating tracers, and other geochemical indicators show that recharge from the Salinas River has the potential to dilute oil-field fluids that might migrate or seep into the aquifer. Groundwater and Salinas River water collected downgradient of the San Ardo Oil Field showed little or no evidence of mixing with oil-field fluids. Some samples within the oil field contained trace amounts of hydrocarbons or elevated temperatures, indicating that any potential effects from oil-field activities are minor or have been diluted by recharge from the Salinas River. The two samples with the most geochemical evidence of potential mixing with oil-field fluids (SP-18 and GW-17) were collected west of or along the Los Lobos fault, where naturally occurring hydrocarbons are near the land surface. Those samples are also near active or inactive oil-field wells, and so anthropogenic activities and pathways cannot be ruled out as a cause of trace detections of hydrocarbons and elevated temperatures in the aquifer.

California↗
USGS2026

Hydrogeology, groundwater salinity distributions, and assessment of the effect of oil-production activities on groundwater in the Midway Valley area, western Kern County, San Joaquin Valley, California

This study seeks to determine the effects of oil field produced water disposal operations and well mechanical integrity issues on groundwater quality in oil fields in the southwest San Joaquin Valley, California. Whereas previous studies used groundwater wells to study shallow aquifers outside the oil fields, this study demonstrates that future approaches may use oil well geophysical logs to map groundwater head gradients, create salinity profiles and document changes in salinity over time in oil field areas with sparse groundwater well data and at depths greater than 330 m. We also incorporate an analysis of well histories to determine potential effects of compromised wellbore seals on changes in aquifer quality that cannot be explained by water disposal practices. Water quality in the aquifers is naturally brackish across most of the area, with better quality groundwater occurring in the eastern part. Geophysical logs are used to determine salinity variations within aquifers including the depth at which TDS exceeds 10,000 mg/L. This depth ranges from 366 m in the northwest to approximately 1,500 m in the southeast. Oil well porosity logs are used to determine water table elevations. These logs indicate the water table slopes south-southeast, showing the predominant groundwater flow direction is from oil field disposal areas toward better quality groundwater east of the oil fields. Geophysical logs show formation resistivity near some disposal facilities has decreased over time, indicating the salinity of the aquifer has increased due to disposal of saline produced water in injection wells and ponds. Oil well history analysis suggests that increased salinity over time in water-saturated sand intervals >1.5 km from disposal facilities may be caused by mechanical failures and/or incomplete borehole seals in poorly constructed or abandoned wellbores prevalent throughout the study area—particularly wells drilled prior to 1930.

California↗
USGS2024

Near-real-time earthquake-induced fatality estimation using crowdsourced data and few-shot large-language models

When a damaging earthquake occurs, immediate information about casualties (e.g., fatalities and injuries) is critical for time-sensitive decision-making by emergency response and aid agencies in the first hours and days. Systems such as the Prompt Assessment of Global Earthquakes for Response (PAGER) by the U.S. Geological Survey (USGS) were developed to provide a forecast of such impacts within about 30 min of any significant earthquake globally. However, existing disaster-induced human loss estimation systems often rely on early casualty reports manually retrieved from global traditional media, which are labor-intensive, time-consuming, and have significant time latencies. Recent approaches use keyword matching and topic modeling to identify human casualty-relevant information from social media but tend to be error-prone when dealing with complex semantics in multi-lingual text data and parsing dynamically changing and conflicting human death and injury numbers shared by various unvetted sources in social media platforms. In this work, we introduce an end-to-end framework to significantly improve the timeliness and accuracy of global earthquake-induced human loss forecasting using multi-lingual, crowdsourced social media. Our framework integrates (i) a hierarchical casualty extraction model built upon large language models, prompt design, and few-shot learning to retrieve quantitative human loss claims from social media, (ii) a physical constraint-aware, dynamic-truth discovery model that discovers the truthful human loss from massive noisy and potentially conflicting human loss claims, and (iii) a Bayesian updating loss projection model that dynamically updates the final loss estimation using discovered truths. We test the framework in real-time on a series of global earthquake events in 2021 and 2022 and show that our framework effectively automates the retrieval of casualty information faster but with comparable accuracy to those now retrieved manually by the USGS. The code associated with this work is made available at: https://github.com/SusuXu-s-Lab/Hierarchical-Earthquake-Casualty-Information-Retrieval

International Journal of Disaster Risk Reduction↗
USGS2026

Site-specific amplifications in Northwestern Turkiye: A generic approach

In this study, we derive generic site amplification functions applicable to Northwestern Türkiye and follow guidance on site classifications established by the United States National Earthquake Hazards Reduction Program (NEHRP). We employ the one-dimensional (1-D) equivalent linear ground response analysis method and focus on recordings from a select suite of free-field strong motion (SM) stations. These particular SM stations are associated with a comprehensive geotechnical site conditions dataset comprising 76 shear-wave velocity profiles. Despite their varying resolutions and depths, this diversity allows for a more representative analysis across different site conditions, enhancing the robustness of our findings. For analyses, we utilize a robust set of global ground motion records, selected from both within and outside of Türkiye, to serve as strong and weak input bedrock motions. Site amplification is determined by propagating these bedrock records through 76 site-specific Afet ve Acil Durum Yönetimi Başkanlığı (AFAD) station profiles and calculating the spectral acceleration ratio between the ground surface and the input motion. Our analyses reveal substantial disparities between the amplification functions derived from 1-D site response analyses using strong and weak input motions. These disparities are central to our analysis and inform the subsequent comparison with site amplification functions reported in the literature. We compare the generic site amplification functions derived in this study with those reported in the literature for regions with similar geological and seismic conditions to Northwestern Türkiye. For example, site amplification functions for NEHRP Site Classes C and D have been documented for Greece and the United States utilizing frequency-dependent methods such as the quarter-wavelength approximation. Our findings reveal that the generic site amplification functions derived in this study exhibit differences in peak amplitudes and predominant frequencies compared to those in the literature, reflecting regional variations in geotechnical and seismic properties. This study represents an initial effort to formulate generic site amplification functions for Türkiye, with a particular emphasis on the Northwestern region. A thorough literature review indicated no prior studies have proposed such functions specifically for this area.

Northwestern Turkiye↗
USGS2025

Mapping eelgrass (Zostera marina) cover and biomass at Izembek Lagoon, Alaska, using in-situ field data and Sentinel-2 satellite imagery

The U.S. Geological Survey and the U.S. Fish and Wildlife Service have developed a three-tiered strategy for monitoring eelgrass ( Zostera marina ) beds at Izembek Lagoon, Alaska, that targets different spatial and temporal scales. The broadest-scale monitoring (tier-1) uses satellite imagery about every 5 years to delineate the spatial extent of eelgrass beds throughout the lagoon. This report describes the most recent (mid-2020s) tier-1 eelgrass monitoring at Izembek Lagoon. The monitoring effort began by canvasing all satellite imagery collected during summer, under clear daytime skies and at low-tide, since the last tier-1 effort in 2006. Two eelgrass maps of Izembek Lagoon were generated by first creating maps of spectrally unique classes from two Sentinel-2 satellite images collected on July 1, 2016, and August 14, 2020, then attributing those spectral classes with information about eelgrass conditions based on field data. Specifically, maps depicting various eelgrass metrics, such as percentage of cover and modeled biomass, were generated using summaries of the ground data that spatially intersected each spectral class. Comparisons of the 2016 and 2020 Sentinel-2 maps showing eelgrass distributional extent, as well as a 2006 Landsat map, indicated that areas where eelgrass presence may have declined during 2006–20 were most prevalent in the central part of Izembek Lagoon. More recently, during 2016-20, areas of possible biomass decline were more prevalent in the southern part of the lagoon. Monitoring eelgrass conditions at Izembek Lagoon with satellite imagery and concurrent ground data allows conditions to be compared over time, but the influences of tide levels, growing season phenology, and spatiotemporal co-registration accuracy should be considered when designing and interpreting change detection analyses.

Alaska↗
USGS2026

Shallow differentiation of primitive arc magmas at the Jurassic Emigrant Gap mafic complex, Sierra Nevada, California

The Emigrant Gap composite pluton exposes ultramafic to silicic intrusive rocks that preserve the chemical evolution of primitive mafic arc magmas and their open-system interactions in the upper crust during mid-Jurassic growth of the Sierra Nevada batholith (California). We present field and petrographic observations and mineral and whole-rock chemistry of the ~35-km 2 ultramafic to dioritic Emigrant Gap mafic complex and an adjacent penecontemporaneous ~90-km 2 granodiorite that together make up the composite pluton. In the Emigrant Gap mafic complex, four roughly central masses of dunite, wehrlite, and olivine clinopyroxenite are surrounded by weakly layered gabbronorite and non-layered diorite. The ultramafic rocks are cumulates formed from near-liquidus minerals of primitive arc magmas that accumulated in steep feeder zones with substantial modification by melt–mush reaction as primitive liquids repeatedly transited the mush-filled conduits. The dominant gabbronoritic rocks are the variably accumulative products of more advanced crystallization–differentiation of arc tholeiitic basalts and basaltic andesites. The adjacent granodiorite intrusion originated separately and preserves field and geochemical evidence for assimilation of metasedimentary rocks. Open-system hybridization between the gabbronoritic mushes and the granodioritic magma produced an intervening body of two-pyroxene diorite. We infer that the ultramafic rocks and gabbronorite of the Emigrant Gap mafic complex crystallized from near-primitive arc basaltic to basaltic andesitic magmas at ~0.15–0.3 GPa, with estimated f O 2 of ≥FMQ +1 and dissolved H 2 O concentrations of only ~0.5–2 wt %. Notably, the Emigrant Gap composite pluton is distinct from other Mesozoic plutons in the Sierra Nevada batholith because of (1) its abundance of mafic and ultramafic rocks that crystallized from relatively primitive mafic melts and (2) the low inferred H 2 O concentrations of its parental magmas, indicated by a near absence of igneous amphibole and by the intermediate rather than calcic compositions of plagioclase. A Jurassic regional extension event probably accounts for the formation of relatively dry primitive arc magmas, as well as for their ascending to the upper crust.

California↗
USGS2025

Relating systematic molecular and textural properties of graptolite pyrolyzed via gold tube hydrous pyrolysis: Implications for thermal proxies in lower Paleozoic marine shales

A series of gold tube pyrolysis experiments (72 h, 300–550 °C, 50 MPa) conducted on a graptolite-rich lower Paleozoic marine shale generated pyrolysis residues for a comprehensive evaluation of the molecular and structural variability of three types of graptolite periderm. Organic petrology, Raman spectroscopy, and field emission scanning electron microscopy (FE-SEM) with energy dispersive spectroscopy (EDS) were combined to evaluate the thermal evolution process. The three types of graptolite periderm, namely granular, non-granular, and nodular graptolite, were analyzed by Raman spectroscopy wherein point measurements were obtained after the maceral was identified and the location verified by organic petrology. Distinct thermal evolution pathways among non-granular, granular, and nodular graptolite periderms were recorded. The evolution patterns of the Raman parameters, particularly D1 and G bands, highlight the differences in geochemical composition of the graptolite periderm types and the alteration of molecular structure with increasing thermal maturity. Raman parameters D1 (position of the D1 peak), G-FWHM (full width at half maximum of the G peak), and ratios D1-FWHM/G-FWHM (full width at half maximum of the D1 peak ratioed to G-FWHM) and A D1 /A G (ratio of D1 and G peak intensities) showed effectiveness in assessing thermal maturity. Bireflectance with increasing gold tube pyrolysis temperature followed a hierarchy: non-granular > granular > nodular, reflecting different molecular alignment intensities. Qualitative FE-SEM evaluation showed that fine-grained mineral inclusions (primarily Fe-sulfide as determined via EDS) were associated with the graptolite populations, with granular graptolite containing greater amounts of coarser-grained (e.g., ∼300–1400 nm) mineral inclusions relative to non-granular and nodular graptolite, which contain finer-grained (e.g., ∼100–200 nm) inclusions difficult to resolve with optical microscopy. These findings are investigated to highlight the mechanisms that drive organic matter evolution within graptolite during thermal maturation, as well as to explore some of the limitations of using spectroscopic parameters as thermal maturity proxies.

International Journal of Coal Geology↗
USGS2026

Sedimentological and geochemical characterization of lacustrine deposits of the Babouri-Figuil basin, northern Cameroon: Implications for source rocks distribution and petroleum exploration

The West and Central African Rift System (WCARS) refers to the series of Cretaceous rift basins where commercial hydrocarbon accumulations have been discovered. Some of the WCARS frontier basins are currently being investigated to increase our understanding of these basins in light of new commercial discoveries. The present study was performed in the Babouri-Figuil Basin (BFB), which is genetically related to the WCARS and constitutes an area of interest in terms of petroleum prospecting, where the distribution of petroleum source rocks and potential targets for petroleum exploration across the entire basin is poorly understood. For the current study, an integrated facies analysis along with organic and inorganic geochemical techniques were applied to the basin's Cretaceous deposits with the aim of reconstructing the paleodepositional environment, assessing factors that triggered the input of organic matter, and providing a spatial overview of the organic matter accumulation in the basin based on outcrop samples. An alluvial fan-lacustrine-braided river system is inferred from the facies analysis of the stratigraphic sequence consisting of conglomerate, sandstone, siltstone, limestone, marlstone, and claystone. Bulk analysis of organic matter reveals that black shale and massive claystone are the main prospective petroleum source rocks in the basin. Inorganic geochemical analyses reveal the influence of anoxic conditions, moderate to high primary productivity, and low terrigenous inputs in organic matter enrichment. The formations rich in organic matter are predominantly concentrated in the western and eastern parts of the basin which may represent areas with depressions, characterized by high accommodation space. In terms of the regional context of the WCARS rift basins, typical hydrocarbon exploration in the BFB may target basal-conglomerate, sandstone beds situated directly above or/and interbedded with the Lower Cretaceous source rocks, and the Upper Cretaceous sandstone beds. Basement rocks (granite, granodiorite, and gneisses) and oil shale deposits may represent potential unconventional hydrocarbon exploration. The current integrated study provides an insight that should guide future hydrocarbon exploration campaigns in the basin.

Babouri-Figuil basin↗
USGS2024

Indications of preferential groundwater seepage feeding northern peatland pools

Groundwater seepage from underlying permeable glacial sedimentary structures, such as eskers, has been hypothesized to directly feed pools in northern peat bogs. These hypotheses directly contradict classical peat bog models for ombrogenous systems, wherein meteoric water is the sole water input to these systems. Variations in the underlying mineral sediment in contact with the peat imply that unrecognized hydrogeologic connectivity may exist with pools in northern peat bogs, particularly where high permeability materials are in contact with the peat. Seepage dynamics originating from these structural variations were investigated using a suite of thermal and hydrogeophysical methods deployed around pools in a peat bog of northeastern Maine, USA. Thermal characterization methods mapped anomalies that were confirmed as matrix seepage or preferential flow pathways (PFPs). Geochemical methods were employed at identified thermal anomalies to confirm upwelling of solute-rich groundwater. Conduits around pools were associated with surficial terminations of suspected peat pipes, based on the inference of pathways extending down into the peat, that focus flow through PFPs in the peat matrix. Discharge also occurred through the peat matrix adjacent to suspected pipe structures and matrix seepage rates were quantified using analysis of diurnal temperature signals recorded at multiple depths. Seepage rates, with a maximum of nearly 0.4 m/d, were measured at localized points around pools. Periods of synchronized temperatures paired with highly muted diurnal temperature signals, recorded in diurnal temperature with depth data, were interpreted qualitatively as activation of strong upward discharge rates through suspected peat pipes. These time periods correlated strongly with local precipitation events around the peatland. Ground-penetrating radar surveys revealed discontinuities in the low permeability glacio-marine clay at the mineral sediment-peat interface, interpreted to be regional glacial esker deposits, which were located beneath and around pools. Heat tracing, specific conductance contrasts, seepage rates, and trace metal concentrations all imply groundwater seepage originating from underlying permeable glacial esker deposits and directly sourcing pools. Preferential groundwater inputs into northern peat bogs may play a key role in developing and maintaining pool systems, with enhanced solute transport impacting peatland ecology, water resources, and carbon cycling.

Maine↗
USGS2026

Magmatic source of the opening phase of Kīlauea’s 2018 Lower East Rift Zone eruption

The 2018 eruption of Kīlauea volcano in its Lower East Rift Zone began with the discharge of evolved high-Ti basalt as weak lava fountains and short, slow-moving lava flows. The lavas were quickly geochemically recognized as being derived from magmas stored within the rift zone and remobilized by a new intrusion, a sequence that is common at Kīlauea. This initial phase of the 2018 eruption, referred to as phase 1a, lasted for 6 days and was followed by extrusion of mixed magma after a 3-day pause. Even though remobilization of older rift zone magmas is common within Kīlauea’s rift zones, it is difficult to determine which past intrusion(s) may have initially emplaced those stored magmas. This difficulty stems from the tendency for Kīlauea magmas to follow very similar differentiation paths without significant variations in major, minor, or even trace element chemistry. We investigate possible magma sources for the lavas erupted during phase 1a of the 2018 eruption using whole-rock, mineral, and glass major and trace element compositions from historical East Rift Zone eruptions with adjacent fissures. We consider two primary hypotheses for the phase 1a source: magmas associated with the 1955 Lower East Rift Zone eruption or the nine eruptions in the Middle and Upper East Rift Zone during the 1960s. Our results suggest that magma associated with the earliest phases of Kīlauea’s 1955 eruption was the most likely source of the 2018 phase 1a remobilized magma. We determine volatile saturation pressures from melt inclusion chemistry and find similar storage depths for the 2018 phase 1a and early 1955 magmas. The phase 1a and early 1955 lavas are nearly indistinguishable in all of the compositional criteria considered, implying that the leftover 1955 magma body barely cooled and differentiated in the 63 years between eruptions (cooling rates of ~0.1 °C/year). This study sheds light on the potential for protracted storage of eruptible magmas in rift zones at Kīlauea, and highlights some of the challenges and solutions to identifying genetic relationships between magmas at Kīlauea.

Hawaii↗
USGS2026

Overstorey mortality promotes juvenile piñon pine growth during favourable weather at cooler, wetter sites

Hotter droughts have resulted in widespread tree die-off events globally, frequently leading to regeneration failure. Dry forest recovery often depends on the growth and survival of extant juvenile trees. However, it is unclear how microenvironmental changes following overstorey tree die-off affect juvenile trees, particularly in dryland systems where tree recruitment is typically limited by water availability and heat stress. We simulated an overstorey tree die-off event by girdling trees in piñon-juniper woodlands across the south-western United States. We sampled juvenile piñon pine growth from live and dead overstorey treatments across six study sites spanning a regional latitudinal gradient and local elevational gradients. We examined how juvenile branch and needle growth differed between live and dead overstorey treatments, and whether responses varied with weather conditions and juvenile tree size following overstorey mortality. We found greater juvenile branch and needle growth under dead compared with live overstorey trees for 2 years following overstorey mortality at mid- and high-elevation sites which are typically cooler and wetter than the other sites. These observed growth releases were contingent on favourable post-mortality weather conditions. Higher growth under dead overstorey occurred at sites experiencing near-average climatic water deficits compared with sites experiencing above-average climatic water deficits. Growth at all sites increased from the first to second year after overstorey mortality. Across sites, growth was unrelated to juvenile tree size. Synthesis . Our results underscore differentiation in juvenile responses to overstorey tree die-off driven by local site conditions and weather across the range of Pinus edulis . Overstorey mortality resulted in consistently higher juvenile growth only at climatically favourable sites and during favourable weather, while unmeasured microsite differences could help account for variation observed at the hottest and driest site. Results from less climatically favourable sites suggest that overstorey trees neither directly limit nor facilitate juvenile growth, though further study over longer timeframes is needed to resolve the pace and magnitude of potential recovery or decline. Overall, juvenile vigour may be promoted following overstorey mortality only in a narrow spatial (site) and temporal (weather) environmental context, suggesting additional vulnerabilities for piñon populations under more arid conditions.

Arizona, Colorado↗
USGS2025

Rare earth element-mineralized carbonatite in the Bear Lodge Alkaline Complex, USA—Ore genesis implications from fluid inclusion characterization

Rare earth element (REE) resources of the Bear Lodge Alkaline Complex, Wyoming, are hosted in variably leached carbonatite spatially related to diatreme breccia pipes. We investigated the genesis of REE and lesser-known gold resources through fluid inclusion analysis of carbonatite, fluorite breccia, and smoky quartz vein samples. Physicochemical characteristics of inclusion-trapped fluids were evaluated using petrography, microthermometry, Raman spectroscopy, decrepitate mound analysis, energy-dispersive spectroscopy, laser ablation inductively coupled plasma mass spectrometry, and noble gas isotope analysis. Microthermometry results reveal three fluid types that affected carbonatite dikes within deeper zones that escaped near-surface, ore-grade REE enrichment: (1) high-temperature (330–432°C) magmatic fluid captured in fine-grained calcite; (2) REE-enriched alkali bicarbonate-sulfate brine; and (3) low-temperature (117–182°C) diluted magmatic or meteoric water. Multiphase brine-like inclusions contain burbankite, nahcolite, strontianite, celestine and alkali sulfate daughter crystals, linking them to early burbankite mineralization. Peripheral smoky quartz and fluorite occurrences at Smith Ridge, 1.5 km from the central carbonatite dike swarm, contain primary inclusions that are Cl-poor and rich in Na-HCO 3 -SO 4 , similar to secondary and pseudosecondary inclusions in carbonatites. Helium isotopes reveal a MORB-like source for carbonatite samples and an older crust signature at Smith Ridge, consistent with the proximal ridge-top exposures of Archean granite. Results from this fluid inclusion study coupled with previous studies of carbonatite mineral paragenesis, show that light REEs (LREEs) were not mobilized great distances. Instead, burbankite crystallized within carbonatite from alkali bicarbonate fluids. With sodium retained in early burbankite, outward-emanating fluids enriched in potassium relative to sodium (higher K:Na) resulted in potassium–ferric iron metasomatism of silicate host rocks. This alkali fractionation was accompanied by fractionation of LREEs and heavy REEs (HREEs), with LREEs dominating the central carbonatite resources. In contrast, areas of peripheral REE mineralization at Bear Lodge are commonly characterized by higher HREE:LREE ratios. The K:Na ratio of associated fenites or alteration assemblages could be indicative of early crystallized burbankite in carbonatites and REE fractionation processes potentially leading to areas of concentrated HREEs with greater supply vulnerabilities.

South Dakota, Wyoming↗
USGS2026

Rare earth element potential in coal and coal ash in the U.S. Gulf Coast

United States heavy reliance on imports of critical minerals (CMs), including rare earth elements (REEs), underscores the importance of development of domestic sources. The study objective was to quantify CM and REE concentrations in coal and coal ash in the US Gulf Coast region. CM and REE concentrations were measured for 118 samples from outcrops and 14 mines in the Gulf Coast. Results show that total REE + Yttrium (REY) concentrations (dry coal basis) are comparable to those of the upper continental crust (UCC) with localized hot spots, such as the Texas Gibbons Creek mine (REY ≤ ~ 2860 ppm). When normalized to UCC REY concentration (169 ppm, dry coal basis), REY to UCC ratios for Gulf Coast coal samples range from 0.1 to 17 (median ratio 0.6). REE extractability from lignites is high (median: 63%–93%) using environmentally benign weak acid. In addition to raw coal, coal ash from power plants could also serve as an REE source with a median ratio of REY in ash relative to coal of 4; however, extractability from coal ash is generally much lower (≤ 5% using the same weak acid as in coal). The median basket price for extracted REY as oxides from coal, assuming 70% extractability, is $3.2 per tonne of coal and $186 billion based on 58 billion metric tonnes of dry coal in the Gulf Coast. REEs important for magnets (Pr + Nd + Tb + Dy) account for ~ 80% of the total value. The corresponding median basket price for extracted REY as oxides from coal ash, assuming ~ 30% extractability, is ~$4.4 per tonne of ash and $1.2 billion based on 258 million tonnes of ash. REE production from coal would likely require co-products, such as activated carbon or humic acids, to attain economic viability. Production of REEs from coal ash could offset remediation costs related to potential water contamination. This reconnaissance study shows the potential for REE production from coal and coal ash in the Gulf Coast; however, carbon coproducts and/or societal benefits would likely be required for socioeconomic viability.

Alabama, Florida, Louisiana, Mississippi, Texas↗
USGS2024

The addition of 144Nd atomic mass to routine ICP-MS analysis as a Quick Screening Tool for Approximating Rare Earth Elements (Q-STAR) in natural waters

Rare earth elements (REEs) are a class of critical minerals, all of which can have supply chain vulnerability that impacts economic security. These elements are widely measured in environmental matrices via inductively coupled plasma mass spectrometry (ICP-MS); however, successful quantification can require time-consuming, sample-specific optimization. While a sample-by-sample approach is appropriate for targeted quantification studies, this approach is not suitable for mineral exploration efforts where rapidly screening thousands of samples for the presence of REEs is desired. Here, we demonstrated the use of a Quick Screening Tool for Approximating REEs (Q-STAR) to detect REEs in surface water and groundwater matrices, collected as part of existing environmental studies. A mass-to-charge ratio of 144 ( m / z = 144) was added to an ICP-MS method to screen for REEs in filtered water samples submitted for metals analyses to the U.S. Geological Survey (USGS) National Water Quality Laboratory. We detected the presence of REEs above a reference threshold of 1200 counts per second in 18 % of pre-selected 6626 samples. Using this screened dataset, we mapped estimated dissolved REE concentrations across the United States in relation to ecoregions and underlying geology. Data are constrained to where sample collection took place but nevertheless show estimated aqueous dissolved REE concentrations on a geographic scale that has not yet been studied. To validate Q-STAR, REEs were measured in a USGS standard reference sample, a subset of 88 archived filtered water samples, and in fresh filtered surface water samples. Our targeted analyses demonstrated a strong linear relationship between Q-STAR predicted and measured values in all archived samples for Nd (r 2 = 0.94), and light REEs (LREEs) such as lanthanum (La) (r 2 = 0.93), praseodymium (Pr) (r 2 = 0.94) and samarium (Sm) (r 2 = 0.94). Using Q-STAR screen values, nine field sites were identified and surface water samples recollected to confirm the continued presence of Nd and LREEs. Q-STAR can be used to screen an unlimited number of water samples for the presence of REEs prior to time-intensive and costly quantitative analyses and to generate large REE datasets for further investigation.

Journal of Geochemical Exploration↗
USGS2025

Evaluation of solid bitumen created from marine oil shale bituminite under hydrous and anhydrous pyrolysis conditions

To test the influence of environmental conditions on aromaticity of solid bitumen generated during petroleum generation, four organic-rich (26–36 wt% total organic carbon) oil shale samples collected from the Neoproterozoic–Lower Cambrian restricted marine Salt Range Formation in the upper Indus Basin, Pakistan, were pyrolyzed under hydrous and anhydrous conditions. Experiments used closed system batch reactors at subcritical water temperatures between 300 and 370°C for 72 h. Thermal conversion of bituminite in the Salt Range oil shales created a newly formed solid bitumen, similar to previous observations from experiments on the Eocene lacustrine Green River Formation Mahogany zone oil shale. Newly formed solid bitumen in the Salt Range Formation oil shales generally has higher reflectance (R o ) in hydrous residues compared to anhydrous experiments subjected to the same time-temperature conditions, also similar to prior observations. This finding supports the hypothesis that radical disproportionation is favored in the presence of hydrogen contributed by water, promoting aromatization in the solid bitumen residue with concomitant increased expulsion of saturated hydrocarbons. Indigenous solid bitumen (and vitrinite in a comparison sample) also shows higher reflectance values in hydrous versus anhydrous residues, indicating that additional aromatization in the presence of exogenous hydrogen occurs both in newly formed organic matter and in organic matter that is present throughout the experiment. Despite similarities in their bulk rock geochemical screening parameters, R o evolution shows different trajectories amongst the four Salt Range oil shales, suggesting as-yet undetermined differences in kinetic properties which are probably related to differences in a priori chemical composition. These results have implications for the use of solid bitumen reflectance (BR o ) as a thermal proxy, suggesting BR o values and appearance could vary as a function of the concentration of water. Variation in water concentration may be present at the reservoir or formation scale, but may also be present at a much finer scale in tight oil shales where permeability is several orders of magnitude lower than conventional reservoirs. Therefore, local variations in the presence of water potentially could explain substantial variation in BR o values and appearance in closely spaced source rock samples and even within an individual microscope field.

Journal of Analytical and Applied Pyrolysis↗