Geology Reports⌕ Search

SEARCH · Geology Reports

Results for “International Journal of Coal Geology”

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

Discrete element investigation of the influence of shallow soil density on the manifestations of strike-slip surface fault rupture

This study investigates the influence of soil relative density on strike-slip surface fault rupture manifestations using three-dimensional numerical simulations performed with the discrete element method (DEM). The simulations capture the formation of distinctive fault strands within complex flower structures using tens of millions of grains. The tendency for dense soils to localize shear manifests as multiple localized shear bands within wide zones of deformation in strike-slip fault rupture, whereas diffuse shear deformation develops within narrow zones in loose soils. The spatial extents of soil deformation are consistent between simulations having similar relative densities but with different quantities of grains. However, individual shears are more distinguishable in assemblages of finer grains than in assemblages of coarser grains. The simulations show the progressive development of new shears within the bounds of previously developed shears. Shear activity transitions inward as fault activity diminishes along the outermost shears and continues along newly developed shears until a vertically dipping throughgoing shear structure develops that accommodates most of the fault displacement thereon. The throughgoing fault develops at smaller fault displacements in looser soils because the first shear rupture propagates closer to the vertical direction and does not undergo as much inward translation of shear activity, as is observed in denser soils. In all simulations, ground surface uplift develops between nonintersecting active shears, and ground surface subsidence tends to develop where new shears intersect previous shears. The surface traces in these simulations are shown to be consistent with analog models and case histories of surface fault rupture occurring in different shallow subsurface materials. Although computationally costly, these modeling results are valuable for providing a strong numerical supplement to traditional analog models used to represent the mechanics of strike-slip zones in soil, and they provide quantifiable stresses and large-strain deformations throughout the model domain.

Journal of Geotechnical and Geoenvironmental Engin↗
USGS2024

Deformation by pressure solution and grain boundary sliding in a retrograde shear zone in southern New England, USA

Alleghanian phyllonites in a shear zone in southern New England were formed by the retrogression and hydration of a high-grade Acadian pelitic schist. The retrogression was locally incomplete, resulting in both heterogeneous mineralogy and mineral compositions, and in many arrested reaction textures. These features, documented by backscattered electron and element map imaging, and by 40 Ar/ 39 Ar dating, provide a record of the processes active as the schist was converted to a phyllonite. On the margins of the shear zone relic biotite and garnet grains survive but are partially replaced by chlorite. Relic high-grade muscovite flakes containing up to 35 % paragonite (pg) also persist but are entrained in swarms of pg-poor muscovite flakes producing folia with an anastomosing network of muscovite grains of variable composition. Single crystals of pg-rich muscovite are truncated by muscovite in successively younger cross-cutting folia that contain decreasing pg contents, falling to <5% pg in the youngest folia. In the core of the shear zone all high-grade minerals are destroyed, but recrystallization continues with earlier greenschist facies folia replaced by lower grade and pg-poorer muscovite in younger folia. Muscovite grains in truncated folia are commonly kinked and selective recrystallization of muscovite in kink bands to pg-poor-compositions demonstrates that strain energy helped drive recrystallization. The boundaries of these new pg-poor muscovite folia are sharp, and truncations indicate that grain boundary sliding (GBS) involved dissolution. GBS also enabled folding by flexural slip along muscovite grain boundaries. Fold amplitudes increased as chlorite dissolved from fold limbs precipitated in fold hinges producing crescent-shaped saddle-reef-like structures. Together these observations of truncation and replacement by dissolution and precipitation demonstrate that pressure solution and GBS facilitated both the formation of the phyllonitic shear zone and slip along it.

Connecticut↗
USGS2024

Following the tug of the audience from complex to simplified hazards maps at Cascade Range volcanoes

Volcano-hazard maps are broadly recognized as important tools for forecasting and managing volcanic crises and for disseminating spatial information to authorities and people at risk. As scientists, we might presume that hazards maps can be developed at the time and with the methods of our discretion, yet the co-production of maps with stakeholder groups, who have programmatic needs of their own, can sway the timing, usability, and acceptance of map products. We examine two volcano hazard map-making efforts by staff at the U.S. Geological Survey. During the 1990s and early 2000s scientists developed a series of hazard assessments and maps with detailed zonations for volcanoes in Washington and Oregon. In 2009, the National Park Service expressed the need for simplified versions of the existing hazard maps for a high-profile visitor center exhibit. This request created an opportunity for scientists to rethink the objectives, scope, content, and map representations of hazards. The primary focus of this article is a discussion of processes used by scientists to distill the most critical information within the official parent maps into a series of simplified maps using criteria specified. We contextualize this project with information about development of the parent maps, public response to the simplified hazard maps, the value of user engagement in mapmaking, and with reference to the abundance of guidance available to the next generation of hazard-mapmakers. We argue that simplified versions of maps should be developed in tandem with any hazard maps that contain technical complexities, not as a replacement, but as a mechanism to broaden awareness of hazards. We found that when scientists endeavor to design vivid and easy-to-understand maps, people in many professions find uses for them within their organization’s information products, resulting in extensive distribution.

Washington↗
USGS2026

The petrogenesis of Þingmúli volcano, East Fjords, Iceland

In this work we revisit Þingmúli volcano (Þ = Th), a classic locality known as an example of a complete tholeiitic differentiation. Þingmúli is a ~ 9.5 Ma extinct central volcano located in the East Fjords of Iceland, in which the whole compositional spectrum from basalt to rhyolites have erupted. These volcanic products have been previously considered as petrogenetically related by an ideal fractionation trend, regardless any temporal relationship or volumetric considerations. Here we report new whole-rock geochemistry, mineral chemistry, isotope analyses, estimation of residence times of the different eruptive deposits, and an update of the original petrogenetic model. Our results highlight that an enriched source, likely spinel lherzolites, generated transitional-alkaline basaltic melts after 15–20% of partial melting at depths of 40–45 km. Many of these basaltic melts erupted at various stages of the volcano's history, while others remained longer in the volcanic plumbing system. These evolved by fractional crystallisation into basaltic andesite magmas with a residence time of ~5 years based on the crystal size distribution of the plagioclase population. Isotopic differences between the basalts/basaltic andesites ( 87 Sr/ 86 Sr ~ 0.7034; 143 Nd/ 144 Nd ~ 0.51315) and the erupted rhyolites ( 87 Sr/ 86 Sr ~ 0.7037; 143 Nd/ 144 Nd ~ 0.51304) indicate that the latter are not petrogenetically related to the former. Therefore, instead of a fractional crystallisation mechanism to generate the rhyolites, we propose the partial melting of ignimbrite layers located beneath the volcano. The broad range of trace element concentrations in andesites and dacites and their different isotopic values compared to the basalts strongly suggest that these magmas have been generated by magma mixing between basaltic and rhyolitic melts, similar to modern day Icelandic volcanoes such as Hekla. These results highlight the need to revisit previously studied Icelandic classic localities and reassess their traditionally proposed petrogenetic models.

Þingmúli volcano, East Fjords↗
USGS2025

Geologic models underpinning the 2024 U.S. Geological Survey assessment of undiscovered oil and gas resources in the Hosston and Travis Peak Formations of the onshore Gulf Coast region, U.S.A.

The Early Cretaceous (Berriasian–Hauterivian) Hosston Formation in Louisiana and eastward is time correlative to the Travis Peak Formation of Texas and southern Arkansas. The formation is a first-order clastic sequence with a regional carbonate transgressive surface as an upper contact. The Hosston and Travis Peak formations contain conventional natural gas and oil accumulations that have been produced for nearly a century. These mature reservoirs contain terrigenous fluvial-deltaic, shore-zone, and paralic deposits across the productive trend; organic-lean mudstone and siltstone lithologies are found outboard of the Lower Cretaceous shelf margin. Producing reservoirs exhibit normal pressure gradients from 0.43 to 0.55 psi/ft (9.7 to 12.4 kpa/m), depths from 4000 to over 20,000 ft (1220 to 6100 m), and temperatures from 150 to 385°F (65 to 196°C). Wells are primarily vertical completions. The number of new field wildcats has been declining since the late 1990s. This paper presents comprehensive geologic models, which include lithofacies maps, structure and isopach maps, burial history models, regional seismic interpretations, and events charts that underpin the recently completed U.S. Geological Survey assessment of undiscovered, technically recoverable hydrocarbons within the Hosston and Travis Peak formations. This study also provides geographic and stratigraphic distributions of Hosston–Travis Peak reservoir properties, including geopressure, reservoir temperature, porosity, permeability, API gravity, and gas-oil ratios. Results indicate estimated undiscovered, technically recoverable mean resources of 28 million barrels of oil and 35.8 trillion cubic ft of gas in conventional and continuous accumulations within the Lower Cretaceous Hosston and Travis Peak formations of the onshore U.S. Gulf Coast region. Quantitative assessment results are detailed in U.S. Geological Survey Fact Sheet 2025–3021 and associated Data Release.

GCAGS Journal↗
USGS2024

Event classification, seismicity, and eruption forecasting at Great Sitkin Volcano, Alaska: 1999–2023

The frequency content of volcanogenic seismicity is often used to classify events and their spatial and temporal progression is then used to map subsurface volcanic processes. The progression of volcano-seismic events and associated source processes also plays a critical role in eruption forecasting. Here we develop and evaluate a computerized methodology for characterizing volcano-seismic event types using Frequency Index and Average Peak Frequency. We apply and test this technique at Great Sitkin Volcano, Alaska, classifying over 9000 hypocenters between 1999 and 2023. This 24-year time span covers periods of seismic quiescence, earthquake activity on nearby tectonic (bookshelf) faults, precursory unrest from 2016 to 2021, and the explosive onset in May 2021 of the ongoing effusive eruption. We use the spatial and temporal evolution of classified event types to map the active volcanic and tectonic processes, develop a conceptual model of the subsurface magmatic system, and perform a retrospective analysis of eruption forecasts at Great Sitkin Volcano between 2016 and the present. The classification and progression of hypocenters suggests the subsurface Great Sitkin Volcano magmatic system consists of a mid- to lower- crustal source zone between 10 and 40 km depth and an upper crustal magma storage area between −1 and 10 km depth (hypocenter depth is referenced to sea level and negative depths reflect height above sea level). The earliest precursors occurred in July 2016 and consisted of deep long-period and volcano-tectonic earthquakes at mid-crustal depths suggesting the subsequent unrest and eruption were triggered by a deeper intrusion of magma. This mid-crustal seismic activity was immediately followed by the onset upper-crustal long-period events and volcano-tectonic earthquakes VTs suggesting a strong linkage between the shallow and deeper portions of the magmatic system. The upper crustal area was likely capped by the 1974 lava dome until the magmatic explosion on May 26, 2021.

Alaska↗
USGS2024

Framework for mapping liquefaction hazard–Targeted design ground motions

Liquefaction-induced ground failure poses substantial challenges to geotechnical earthquake engineering design. Current approaches for designing against liquefaction hazards, as specified in most seismic provisions, focus on estimating a liquefaction factor of safety ( 𝐹⁢𝑆𝐿 ) and typically characterize earthquake loading using design parameters based on probabilistic or deterministic ground motion levels. Because 𝐹⁢𝑆𝐿 is estimated deterministically, this basis of design neglects considerable uncertainties for estimating liquefaction triggering and its consequences and results in a lack of liquefaction-specific design criteria, particularly as structural design has advanced toward risk-targeted performance objectives. This study presents a framework for developing liquefaction-targeted design criteria based on a minimum acceptable return period of liquefaction, informed by probabilistic liquefaction hazard analysis (PLHA). PLHA quantifies annualized rates of liquefaction by considering contributions from (1) the full ground-motion probability space, and (2) uncertainties in liquefaction triggering using probabilistic models. PLHA is used in this study to characterize the current, effective return periods of 𝐹⁢𝑆𝐿 ( 𝑇𝑅,𝐹⁢𝑆 ) obtained from conventional liquefaction hazard analysis (CLHA) using uniform-hazard ground motions. 𝑇𝑅,𝐹⁢𝑆 is evaluated in a parametric study of nearly 100 sites throughout the conterminous United States. The results indicate large geographic variations in acceptable liquefaction hazard levels, with implied 𝑇𝑅,𝐹⁢𝑆 ranging between approximately 1,000 to 3,000 years. To address these inconsistencies without the computational demands of full PLHA, a framework is proposed for developing a liquefaction-targeted design peak ground acceleration, 𝑃⁢𝐺⁢𝐴𝐿 , for use in liquefaction models that result in consistent liquefaction design levels across all geographic locations. The mapped 𝑃⁢𝐺⁢𝐴𝐿 is shown to be somewhat sensitive to site-specific properties, and adjustment factors are developed and presented. The proposed 𝑃⁢𝐺⁢𝐴𝐿 mapping procedure produces 𝐹⁢𝑆𝐿 estimates that are consistent with those obtained from full PLHA at a target 𝑇𝑅,𝐹⁢𝑆 , providing a promising roadmap to incorporating PLHA concepts into current liquefaction design methods.

Journal of Geotechnical and Geoenvironmental Engin↗
USGS2026

Nonstationary demographic state-space models using unreplicated counts for species undergoing environmental stressors

A fundamental task in ecological statistics is to estimate abundance and growth rate distributions from wildlife monitoring data to inform conservation management. Modeling time series of wildlife populations presents a number of challenges from both statistical and ecological perspectives, including discreteness; lack of replication; nonstationarity; and observation, demographic, and other phenomenological processes. Nonstationary dynamics are often exhibited by populations undergoing environmental stressors. Models must account for these characteristics to produce reliable estimates of abundance and trends, yet estimation can be challenging with unreplicated data. We propose nonstationary demographic state-space models using unreplicated counts for populations undergoing environmental stressors. A reduced growth rate model matches the complexity of the unreplicated count data, and a fecundity bound on growth rate distributions allows the separation of processes affecting growth rates like environmental stressors from those affecting abundance external to growth rates like migration. NDSSMs allow for the embedding of nonstationary model components, and we explore the use of changepoints, volatility clustering, and migration processes. We apply the proposed nonstationary models in case studies of herons affected by predator/competitor reestablishment and three bat species affected by a fungal pathogen causing white-nose syndrome. Nonstationary models outperform stationary models and generalized linear mixed effects models according to model scoring and visual inspection of predictions, and provide estimates more consistent with published values. Incorporating migration improves model fit universally, even with approximate one-way immigration, most likely because populations are extirpated, recolonized, and increase multiple-fold over the upper bound set by species fecundity. In addition, estimates of the timing and severity of the environmental stressor differed for models with migration. Including nonstationary and demographic components in a fecundity-bounded growth rate model improves inference and benefits interpretability of hyperparameters. In turn, this adjusts uncertainties in predictions of abundance and growth rates over time, providing the ingredients needed for informed conservation analysis and for directing future monitoring of at-risk species.

Journal of Agricultural, Biological and Environmen↗
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

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

Solution-collapse breccia pipe uranium deposits of the southern Colorado Plateau, northwestern Arizona, USA

Some of the highest-grade uranium deposits in the United States occur in breccia pipes that formed by solution and collapse of sedimentary strata, which occur in the southern portion of the Colorado Plateau in northwestern Arizona. The host breccia pipes are up to 1200 m in vertical extent, average about 90 m in diameter, and can cross-cut strata from their base in the Mississippian Redwall Limestone to as stratigraphically high on some plateaus as the Triassic Chinle Formation. These uranium-base metal deposits are up to 600 m thick and formed within the breccia pipes where they transect the Permian Coconino Sandstone, Hermit Formation, and the Esplanade Sandstone. Of the hundreds of breccia pipes identified across this region, only a small percentage are known to contain mineralization. The main uranium ore mineral is uraninite that is intergrown with at least 20 base-metal sulfide minerals, which contribute Fe, Cu, Co, As, Pb, Zn, Ni, and Ag to the deposits. This study considered regional stratigraphy, sulfur isotope systematics, mineralogy, in situ dating, and compilation and analysis of previous work on the deposits. A comprehensive deposit model has not been published for these deposits. This analysis identified new additions to update the deposit model for these unusual, possibly unique deposits. Proposed modifications to the model include: (1) the source, mechanisms, timing of the base-metal sulfide mineral assemblages, and (2) the source, mechanism, and timing of the uranium mineralization. Sulfide and uranium deposition are shown to be separate mineralization events. The study proposes the possible role of gypsum as a source of sulfur for the sulfide minerals in the deposits. Groundwaters carrying uranium encountered the preexisting sulfides in breccia pipes, reducing the uranyl ions, and precipitating U oxide (as uraninite). Analysis of the regional stratigraphy recognized that numerous beds of gypsum are in the strata that lie only tens of meters above the breccia pipe deposits. In the breccia pipe region, if these stratigraphic units (Toroweap and Kaibab Formations) do not contain gypsum layers then the underlying pipes are not mineralized; where these Permian gypsum layers do occur, breccia pipes can host mineralization. This new understanding should be useful in identifying the prospective region for mineralized pipes.

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

Cold blood in warming waters: Effects of air temperature, precipitation, and groundwater on Gulf Sturgeon thermal habitats in a changing climate

Objective In a changing climate, the effects of air temperature, precipitation, and groundwater on water temperature and thermal habitat suitability for Gulf Sturgeon Acipenser desotoi , listed as threatened under the U.S. Endangered Species Act, are not well understood. Hence, we incorporated these factors into thermal habitat models to forecast how Gulf Sturgeon may be affected by wide‐ranging climate change scenarios in 2024–2074. Methods Using data from the Choctawhatchee River, Florida, we developed precipitation‐ and groundwater‐corrected air–water temperature models, compared their accuracy with that of conventional air–water temperature models used in fisheries management, and projected future Gulf Sturgeon thermal habitat suitability for normal physiological functioning and fieldwork (i.e., population sampling and telemetry surgeries) in summer (May–August) under 16 climate change scenarios. Result Precipitation‐ and groundwater‐corrected models were more accurate than conventional air–water temperature models (mean improvement in adjusted R 2 = +0.45; range = +0.09 to +0.75). Water temperature was projected to warm at widely variable rates across climate change scenarios encompassing different air temperature, precipitation, and groundwater regimes. Importantly, Gulf Sturgeon summer aggregation areas were cooler and influenced more by precipitation and groundwater and less by air temperature than were non‐aggregation areas. If precipitation and groundwater—as drivers of cooling—become warm in a changing climate, summer aggregation areas were projected to exhibit thermal habitat degradation equivalent to or greater than that of non‐aggregation areas. Conclusion Our results add hydrological context to the premise that aggregation areas provide cool water and energetic savings for Gulf Sturgeon during summer, underscoring the importance of protecting these habitats through groundwater conservation, water quality monitoring, and riparian/watershed habitat management. Our findings indicate that identifying thermally appropriate times for fieldwork activities will be increasingly important and time‐restricted as climate change intensifies. However, our research provides managers with a portfolio of water temperature models and an accurate, cost‐effective, management‐relevant approach to forecasting thermal habitat conditions for Gulf Sturgeon and other species in a changing climate.

Alabama, Florida↗
USGS2025

Forecasting water levels using the ConvLSTM algorithm in the Everglades, USA

Forecasting water levels in complex ecosystems like wetlands can support effective water resource management, ecological conservation, and understanding surface and groundwater hydrology. Predictive models can be used to simulate the complex interactions among natural processes, hydrometeorological factors, and human activities. The Greater Everglades in the USA is a well-known example of an ecosystem where complexity has motivated adoption of machine learning algorithms in water level prediction studies. This paper aims to contribute to extending existing machine learning algorithms by integrating spatiotemporal data with deep-learning algorithms in the forecasting process. In this study, a deep-learning model is developed to predict water levels on a regional scale, covering a large area of approximately 9,138 square kilometers in the Everglades ecosystem. This model has the architecture of Convolutional Long Short-Term Memory which can deal with spatiotemporal data by capturing both spatial and temporal dependencies in the training data. The forecasting capabilities of this model (referred to as the global model) are assessed by comparing the global model to two Artificial Neural Networks developed at two different gaging stations, referred to here as local models. One local model is developed at a gaging station directly influenced by nearby water control structures, whereas the other is developed at a gaging station located farther away from these structures. By leveraging data from the Everglades Depth Estimation Network spanning from January 2002 to May 2023, the global and local models were trained to forecast water levels with a two-day lead time. Our findings suggest that both the global and local models perform with approximately the same level of accuracy, with Mean Absolute Relative Error values ranging from 0.38% to 1.4% at the selected stations. The developed global model has demonstrated strong potential as a standalone forecasting tool for the entire study area in the Everglades and could eliminate the need for developing multiple local models. This finding also highlights how machine learning can capture complex spatial and temporal relationships to generate accurate water level predictions on a regional scale.

Florida↗
USGS2024

Assessing potential effects of oil and gas development activities on groundwater quality near and overlying the Elk Hills and North Coles Levee Oil Fields, San Joaquin Valley, California

Groundwater resources are utilized near areas of intensive oil and gas development in California’s San Joaquin Valley. In this study, we examined chemical and isotopic data to assess if thermogenic gas or saline water from oil producing formations have mixed with groundwater near the Elk Hills and North Coles Levee Oil Fields in the southwestern San Joaquin Valley. Major ion concentrations and stable isotope compositions were largely consistent with natural processes, including mixing of different recharge sources and water-rock interactions. Trace methane concentrations likely resulted from microbial rather than thermogenic sources. Trace concentrations of benzene and other dissolved hydrocarbons in three wells had uncertain sources that could occur naturally or be derived from oil and gas development activities or other anthropogenic sources. In the mid-1990s, two industrial supply wells had increasing Cl and B concentrations likely explained by mixing with up to 15 percent saline oil-field water injected for disposal in nearby injection disposal wells. Shallow groundwater along the western margin of Buena Vista Lake Bed had elevated Cl, B, and SO 4 concentrations that could be explained by accumulation of salts during natural wetting and drying cycles or, alternatively, legacy surface disposal of saline oil-field water in upgradient ephemeral drainages. This study showed that groundwater had relatively little evidence of thermogenic gas or saline water from oil and gas sources in most parts of the study area. However, the evidence for groundwater mixing with injected disposal water, and possibly legacy surface disposal water, demonstrates produced water management practices as a potential risk factor for groundwater-quality degradation near oil and gas fields. Additional studies in the San Joaquin Valley and elsewhere could improve understanding of such risks by assessing the locations, volumes, and types of produced water disposal practices used during the life of oil fields.

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