USGS Science⌕ Search

SEARCH · USGS Science

Results for “Article”

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.

979 records · Page 48Linked to original sources
USGS2026

Bottom-up characterization of geologic methane emissions in the San Juan Basin in the southwestern USA

Methane is a potent greenhouse gas that plays an important role in atmospheric chemistry and global warming. The current global methane budget has large uncertainties, and a better understanding of the budget would help to guide strategies for reducing anthropogenic emissions to fight climate change. Natural geologic methane emissions are a particularly poorly constrained source, with top-down estimates from 14 C in ice cores suggesting much lower geologic emissions than bottom-up scaling of direct flux measurements. Our study aims to contribute to resolving this discrepancy through improved bottom-up characterization of geologic methane seepage in the San Juan Basin in southwestern Colorado and northwestern New Mexico, USA. We performed 983 new flux chamber measurements in this basin during summer 2022 and winter 2023 field campaigns. Our results, in combination with prior measurements, suggest that natural seepage in the San Juan Basin only occurs on or near the Fruitland coal outcrop. Specifically, our new measurements confirm previous measurements of seepage along the northwestern exposure of the Fruitland outcrop in Colorado (a known hydrodynamic overpressure region) and for the first time, identified seepage locations along the southernmost Fruitland outcrop exposure in New Mexico, in association with a coal cleat and a fault. Overall, seepage along the Fruitland coal outcrop is heterogeneously distributed, with both positive and negative (interpreted as microbial soil sink) methane fluxes. Features that are hypothesized to be predictive of seepage (e.g., faults) were not associated with positive methane fluxes in areas outside of the Fruitland outcrop. Our best estimate for total geologic methane seepage in the San Juan Basin from spatial interpolation and statistical upscaling is approximately 0.14 Tg CH 4 /yr, with a range from 0.029 to 0.48 Tg CH 4 /yr. This best-estimate value is lower than a previous bottom-up estimate from a gridded seepage inventory, but higher than a previous top-down estimate.

Colorado, New Mexico↗
USGS2024

3-D geological modeling for numerical flow simulation studies of gas hydrate reservoirs at the Kuparuk State 7-11-12 Pad in the Prudhoe Bay Unit on the Alaska North Slope

Accurate reservoir evaluation requires reliable three-dimensional (3-D) geological models. This study conducted 3-D geological modeling for numerical flow simulation of the B1 sand gas hydrate reservoir at the Kuparuk State 7-11-12 pad, Prudhoe Bay Unit, Alaska North Slope. The model integrates well logs, core, and seismic data to address spatial heterogeneity in geological structures and reservoir properties. Two modeling types were performed: structural framework modeling and petrophysical property modeling. For structural framework modeling, seismic data and well log markers were used to reproduce subsurface structures characterized by a normal fault system. A volume-based modeling algorithm and stair-step gridding were applied. The resulting 3-D model comprised 2,640,000 grid cells across 264 layers, including seven fault grids. For petrophysical property modeling, total porosity was initially modeled using sequential Gaussian simulation with collocated cokriging. To reproduce the upward coarsening of the B1 sand, upscaled log-derived total porosity and a 3-D trend depicting total porosity variation were used as primary and secondary data, respectively. Gas hydrate saturation distribution was modeled similarly, with secondary data from estimated porosity distribution and seismic-derived acoustic impedance map enhancing accuracy. Results indicate higher gas hydrate saturation in the upper part of the B1 sand and areas with higher acoustic impedance. Intrinsic permeability was modeled from the total porosity and clay-bound water volume, and effective permeability was derived from the gas hydrate saturation and intrinsic permeability distributions based on the “Tokyo model”. Effective permeability distributions were influenced by the total porosity, gas hydrate saturation, and intrinsic permeability. Within the same layer, higher gas hydrate saturation leads to decreased effective permeability. In total, 100 sets of multiple scenarios were prepared, providing input data for dynamic flow simulations to evaluate the effects of lateral heterogeneity in reservoir properties and the hydraulic characteristics of faults on production behavior for preassessment before the long-term production test.

Alaska↗
USGS2025

Reservoir thermal energy storage pre-assessment for the United States

Storing thermal energy underground for later use in electricity production or direct-use heating/cooling is a promising, viable, and economical green energy option. Reservoir thermal energy storage (RTES) is one such option, which stores energy in underutilized permeable strata with low ambient groundwater flow rates and more geochemically evolved (e.g. brackish/saline) waters relative to overlying principal aquifer systems. The U.S. Geological Survey has begun assessing RTES potential nationally by focusing on five generalized geologic regions (Basin and Range, Coastal Plain, Illinois Basin, Michigan Basin, Pacific Northwest) across the United States. Hydrogeologic reservoir models are developed for the following eight metropolitan area cities within those regions to evaluate RTES performance across different climates and subsurface conditions: Albuquerque, New Mexico; Charleston, South Carolina; Chicago and Decatur, Illinois; Lansing, Michigan; Memphis, Tennessee; Phoenix, Arizona; and Portland, Oregon. Evaluated metrics include estimated required well spacing, thermal storage capacity, and thermal recovery efficiency through time. Also considered for each reservoir are potential complicating factors, including reservoir depth, thermally driven free convection, and groundwater salinity. This work focuses on direct-use cooling because the need for cooling modern office buildings greatly exceeds that for heating in most parts of the country (Falta and others, 2016); however, the evaluated metrics are also relevant to heating and electricity applications. Results indicate that favorable RTES conditions exist in each region, with the Coastal Plain and Basin and Range being especially favorable for thermal storage capacity, while the Pacific Northwest and Michigan Basin excel at energy recovery for the evaluated cooling application. The results underscore the utility of developing maps of thermal storage capacity, subsurface temperature models, and volumetric estimates of thermal storage capacity to serve as key RTES resource classification standards. Overall, this pre-assessment provides a basic understanding of RTES potential in several cities and geologic regions throughout the country and will aid ongoing thermal energy storage assessment efforts.

Arizona, Illinois, Michigan, New Mexico, Oregon, S↗
USGS2026

Evidence for organic and igneous carbon sources in the Mesoproterozoic Adirondack Highlands graphite district, New York

This study examined the origin of graphite mineralization in the Mesoproterozoic Adirondack Highlands. Flake graphite was mined from the historic mining district around Ticonderoga, New York, USA, from the 1850s to the 1920s. Most of the production and reserves from the district are disseminated graphite from feldspathic quartz schists (e.g., Dixon-American Graphite Co. mine, Faxon property, Flake Graphite Company). Protolith sediments of the schists were deposited in the Trans-Adirondack backarc basin at ca. 1.25 Ga and were metamorphosed to high grade at 1.18–1.15 Ga during the accretionary Shawinigan orogeny and at 1.09–1.02 Ga during the Ottawan phase of the Grenvillian orogeny. Major elements and carbon isotopes of graphite from schists (δ 13 C = −28.1‰ to −5.4‰ relative to Vienna Peedee belemnite) are consistent with derivation from organic carbon in protolith sandstones that devolatilized during metamorphism. Some small mines and prospects in the district mined graphitic calcite marbles and clinopyroxene ± scapolite skarns associated with a pegmatite suite dated by U-Pb in zircon to 1.04–1.01 Ga (e.g., Lead Hill and Crown Point deposits), and rare graphite veins are observed at several deposits. Hydrothermal flake graphite from pegmatites, skarns, and veins are distinct from graphite in schist-hosted deposits and have δ 13 C values that cluster around −7‰ to −5‰, consistent with deposition from igneous CO 2 . Covariation of carbon and oxygen stable isotopes and U-Pb zircon geochronology are not compatible with precipitation from fluids derived from country rock metapelites or marbles. Some graphite districts elsewhere show compelling evidence for remobilization and redeposition of sedimentary-derived carbon as hydrothermal graphite, sometimes as thick veins producing world-class, high-grade deposits. In the Adirondacks, there is evidence for two generations of graphite mineralization with distinct carbon sources: syngenetic graphite with an organic origin in Shawinigan quartz schists and epigenetic graphite with an igneous origin in Ottawan skarns and pegmatites.

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

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

Near-surface material and topography generate anomalous high-frequency ground motion amplification in Chugiak, Alaska

An ∼3 km long nodal array oriented approximately east–west was deployed in Chugiak, Alaska, by the U.S. Geological Survey during 2021. The array intersects with the permanent NetQuakes station NP.ARTY, where peak ground acceleration (PGA) value of 1.98 g was recorded during the 2018 M w 7.1 Anchorage, Alaska, earthquake, in sharp contrast to the PGA of ∼0.3 g at a site just 4 km to the west. Seismic data for M w 1.8–4.3 aftershocks from the M w 7.1 event recorded by the nodal array confirm the anomalously large ground motions obtained at NP.ARTY as well as similar amplifications at nodes within ∼1 km to the east. Here, we performed 0–10 Hz 3D finite‐difference simulations, including high‐resolution surface topography, to explore the cause of the unexpectedly large amplification. As expected, the simulations computed with a regional 3D tomography velocity model severely underpredict the 0–10 Hz acceleration records at almost all sites. Adding a near‐surface low‐velocity taper to 300 m depth amplifies the accelerations by up to a factor of 5 and enables a reasonable match between the nodal data and simulations at sites to the west of NP.ARTY. However, this model still underpredicts the spectral energy in the area covered by glacial sediments by up to an order of magnitude. The addition of a till layer using a depth‐dependent shear‐wave velocity ( ⁠⁠ V s ) profile along with a homogeneous, 8 m thick low‐velocity layer with V s = 250 m/s representing the kame terraces improves the fit to data to within a factor of 2 at nodes located on top of the glacial sediments. Our study shows that the anomalously large high‐frequency amplification recorded at and near NP.ARTY can be explained by a combination of topographic effects and near‐surface low‐velocity material with amplification effects on the high‐frequency ground motion by up to about 40% and an order of magnitude, respectively.

Alaska↗
USGS2025

Decomposing the Tea Bag Index and finding slower organic matter loss rates at higher elevations and deeper soil horizons in a minerogenic salt marsh

Environmental gradients can affect organic matter decay within and across wetlands and contribute to spatial heterogeneity in soil carbon stocks. We tested the sensitivity of decay rates to tidal flooding and soil depth in a minerogenic salt marsh using the Tea Bag Index (TBI). Tea bags were buried at 10 and 50 cm depths across an elevation gradient in a subtropical Spartina alterniflora marsh in Georgia (USA). Plant and animal communities and soil properties were characterized once, while replicate tea bags and porewaters were collected several times over 1 year. TBI decay rates were faster than prior litterbag studies in the same marsh, largely due to rapid green tea loss. Rooibos tea decay rates were more comparable to natural marsh litter, potentially suggesting that is more useful as a standardized organic matter proxy than green tea. Decay was slowest at higher marsh elevations and not consistently related to other biotic (e.g., plants, crab burrows) or abiotic factors (e.g., porewater chemistry), indicating that local hydrology strongly affected organic matter loss rates. TBI rates were 32 %–118 % faster in the 10 cm horizon than at 50 cm. Rates were fastest in the first 3 months and slowed 54 %–60 % at both depths between 3 and 6 months. Rates slowed further between 6 and 12 months, but this was more muted at 10 cm (17 %) compared to 50 cm (50 %). Slower rates at depth and with time were unlikely due to the TBI stabilization factor, which was similar across depths and decreased from 6 to 12 months. Slower decay at 50 cm demonstrates that rates were constrained by environmental conditions in the deeper horizon rather than the composition of this highly standardized litter. Overall, these patterns suggest that hydrological setting, which affects oxidant introduction and reactant removal and is often overlooked in marsh decomposition studies, may be a particularly important control on organic matter loss in the short term (3–12 months).

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

Uncertainty quantification of geologic energy storage in depleted gas reservoirs using material balance equations embedded in a hierarchical errors-in-variables model

The storage potential for gas in geologic settings, such as depleted hydrocarbon reservoirs and solution-mined salt caverns, is becoming salient to future energy infrastructure planning. Technologies such as carbon capture, utilization, and storage, carbon dioxide-enhanced oil recovery, and natural gas and hydrogen storage help meet growing energy demands, reduce carbon emissions to meet climate goals, and provide energy security amid geopolitical uncertainties. Therefore, estimates of underground gas storage capacity could be useful for efficiently navigating the energy transitions. Material balance is a fundamental method in reservoir engineering for estimating original gas in place and potential storage capacity at the scale necessary for national assessments of subsurface pore space resources. However, the deterministic method cannot accommodate multiple data sources or quantify uncertainty in predictions. In this study, a method that embeds material balance equations within a hierarchical errors-in-variables model is proposed which allows the estimation of the distributions of reservoir properties needed for assessments. Uncertainties associated with these reservoir properties have traditionally been expert-elicited, whereas the uncertainty estimates from the proposed models are data-driven. Capacity and uncertainty estimates can be used in a probabilistic resource assessment, supplementing information traditionally used by assessors or even replacing this expert elicitation step when data are unavailable. Various regression models are compared in a case study of the Michigan Basin, a large contributor to the United States’ current natural gas storage capacity. In particular, errors-in-variables models help ameliorate regression dilution and can quantify uncertainty in predictions of pressure in addition to storage capacity. Overfitting is addressed by quantifying generalization error and model averaging in simple and stratified cross-validation against reported working gas capacity, representing the varying quality and quantity of available data. Incorporating a statistical framework into existing numerical methods in reservoir engineering can improve the quality of estimation, and in particular, this method brings rigor to uncertainty quantification as part of a larger effort by the U.S. Geological Survey to assess domestic energy gas storage resources in depleted hydrocarbon reservoirs.

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

Geologic input databases for the 2025 Puerto Rico – U.S. Virgin Islands National Seismic Hazard Model update: Crustal faults component

The last National Seismic Hazard Model (NSHM) for Puerto Rico and the U.S. Virgin Islands (PRVI) was published in 2003. In advance of the 2025 PRVI NSHM update, we created three geologic input databases to summarize new onshore and offshore fault source information in the northern Caribbean region between 62°–70° W and 16°–21° N. These databases, of fault sections, fault‐zone polygons, and geologic estimates of fault activity (fault‐slip rate and earthquake recurrence intervals) at specific sites, document updates to fault parameters used in prior seismic hazard models in PRVI. Fault sources were reviewed from published studies since 2003, which document substantial changes to the understanding of fault location, geometry, or activity. New fault section sources were added for features that meet the criteria of (1) length ≥7 km, (2) unequivocal evidence of recurrent tectonic Quaternary activity, and (3) documentation that is publicly available in a peer‐reviewed source. In addition, we revised several broad areal sources, such as the Mona and Anegada extensional zones. The 2003 model included three fault sections and two fault‐zone polygons (areal sources). These databases include 35 fault sections, 6 fault‐zone polygons, and 51 earthquake geology sites. To characterize fault activity rates, slip‐rate bins were assigned based on landscape expression and paleoseismic trench observations for faults without published slip‐rate sites. Additional fault sources were evaluated but not included in these databases due to a lack of published information about fault location, geometry, or recurrent Quaternary activity. The PRVI NSHM 2025 geologic input databases describe crustal faulting; the geometries and coupling of Puerto Rico subduction zone and Muertos Trough models are considered in a separate database. Updates to the fault sections, fault‐zone polygons, and earthquake geology databases can help inform the location and recurrence rate of damaging earthquakes in the PRVI NSHM implementation.

Puerto Rico, U.S. Virgin Islands↗
USGS2024

A comparative analysis of OpenET for evaluating evapotranspiration in California almond orchards

The almond industry in California faces water management challenges that are being exacerbated by droughts, climate change, and groundwater sustainability legislation. The Tree-crop Remote sensing of Evapotranspiration eXperiment (T-REX) aims to explore opportunities to improve precision irrigation management for woody perennial cropping systems. Almond orchards in the California Central Valley were equipped with eddy covariance flux measurements to evaluate satellite remote sensing-based evapotranspiration (RSET) models. OpenET provides high-resolution (30-m spatial and daily temporal) RSET data, synthesizing decades of research for practical water management. This study provides an evaluation of OpenET performance at six almond sites covering a large range in soils, age, and variety. It also compares OpenET ensemble evapotranspiration (ET) data with applied irrigation and precipitation records over an additional 148 almond orchards located in the Central Valley of California. Results show OpenET models, including the ensemble ET value, produced reasonable and actionable ET values, with overall coefficient of determination (R 2 ) and mean absolute error values of 0.73- and 0.95-mm d −1 at the daily time step, respectively. However, given the temporal sampling of Landsat (8-day revisit) and the interpolation methods used, the assessed ET models had difficulty in capturing short-term variability in almond ET; for example, the rapid decline in measured ET observed as a response to lack of irrigation preceding and during almond harvest. The study also drew attention to the spatial complexity in scenarios where irrigated orchards are surrounded by hot/dry areas, causing discrepancies between measured and modeled ET values. In comparison with irrigation records, OpenET ensemble ET was capable of quantifying water input (applied irrigation + precipitation) in almond orchards to within 13 % when evaluating monthly data. Initial results presented here reinforce the idea that RSET models, such as in OpenET, are powerful tools, yet their application requires nuanced understanding and careful consideration of local conditions.

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

Application of non-stationary shear-wave velocity randomization approach to predict 1D seismic site response and its variability at two downhole array recordings

Accounting for uncertainties in seismic site response is crucial to improving the performance of one-dimensional (1D) ground response analyses (GRAs) at downhole array recording sites. In addition to site effects, uncertainties in 1D-GRAs can also be contributed from the seismic source and/or path. Though often representing not more than one percent of the distance (path) from the source, site conditions are known to have an enormous influence on ground shaking. In this study, we focus on the site shear-wave velocity ( V S ) structure, which is the main ingredient for estimating the variability of site response. As such, V S can manifest aleatory uncertainties related to the effects of small-scale spatial heterogeneities within the near surface, thus V S can substantially modify ground shaking during earthquakes. We apply a novel V S randomization approach to propagate the small-scale heterogeneities of V S to estimate seismic site response within a non-stationary probabilistic framework. The randomization approach generates samples of V S profiles that are used to perform several 1D-GRAs and obtain an averaged site response and related variability. The proposed method is implemented on data recorded at two downhole array sites with different subsurface soil conditions: a soft soil site on Treasure Island (California, United States of America) and a rock outcrop site in Cadarache (South-East France). We show that synthetic surface-to-borehole transfer functions from 1D-GRAs provide an acceptable fit to the empirical transfer functions from low-motion earthquake records and succeed in reproducing most of the site-specific seismic response variability. The remaining mismatch between transfer functions is likely due to insufficient precision on the seismic bedrock and the impedance contrast. The variability in site response is discussed with emphasis on the role of V S small-scale heterogeneities, attenuation, and input motion incidence angle in ground motion variability for the site and soil conditions at both locations.

California↗
USGS2025

Leveraging wildfire to augment forest management and amplify forest resilience

Successive catastrophic wildfire seasons in western North America have escalated the urgency around reducing fire risk to communities and ecosystems. In historically frequent-fire forests, fuel buildup as a result of fire exclusion is contributing to increased fire severity. The probability of high-severity fire can be reduced by active forest management that reduces fuels, prompting federal and state agencies to commit significant resources to increase the pace and scale of fuel reduction treatments. However, lower severity areas of wildfires also have the potential to act as “treatments,” and even catastrophic fires with large areas of high severity can still have substantial areas of lower severity fire that may be improving forest conditions locally. We quantified active management and wildfire severity across yellow pine and mixed conifer (YPMC) forests in the Sierra Nevada of California over a 22-year period (2001–2022). We did not detect increases in the area treated through time, but the area of beneficial wildfire (low to moderate severity) increased substantially, exceeding active treatment area in 8 of 22 years. Overall, beneficial wildfire treated ~17% more area than all treatments combined, and roughly four times more area than fire-related treatments alone. We then used disturbance history to evaluate resistance to high-severity wildfire and forest loss across the YPMC range. Of the 2.3 million ha YPMC of forests in 2001, 20% lost mature forests due to high-severity fire by 2022, which is nearly half of all YPMC area burned. Most of the landscape (47%) remains at risk of high-severity fire because it had no restorative disturbances, but 33% of the study area has some level of resistance to high-severity wildfire. In these areas, resistance will need to be enhanced and maintained over time via active management or managed wildfire, but these treatment needs will likely outpace capacity even under optimistic implementation scenarios. Given limited resources for implementing active management and the likelihood of a more fiery future, incorporating beneficial wildfire into landscape-level treatment planning has the potential to amplify the impact of active management treatments.

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