USGS Science⌕ Search

SEARCH · USGS Science

Results for “Scientific Investigations Map”

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.

862 records · Page 48Linked to original sources
USGS2025

Confocal laser-scanning microscopy (CLSM)-based thermal maturity of Tasmanites and progress in standardization of fluorescence microspectrometry

Evaluation of thermal maturity in vitrinite-free or vitrinite-deficient sediments via fluorescence microspectrometry can provide relevant information related to petroleum exploration and thermal history assessment. However, variation in spectral fluorescence properties of alginite macerals with increasing thermal maturity is largely underexplored. Here, authors of this study have applied confocal laser-scanning microscopy (CLSM) in conjunction with fluorescence microspectrometry to a maturity series of marine Upper Devonian Tasmanites algae from the Ohio Shale (Huron Member) and a single sample from the Marcellus Formation of the Appalachian Basin. Spectral fluorescence properties of Tasmanites were evaluated in relation to orientation, measurement location, and the number of measurements per sample, and were compared to published literature. Emission spectra of Tasmanites from continuous wave laser excitation (405 nm) were acquired from sections perpendicular and parallel to bedding and at different positions within individual Tasmanites bodies. The results showed a progressive red-shift in emission maxima (λ max ) in a large sample sized maturity series ( N = 19), e.g., 493 to 578 nm for the perpendicular section at middle position. Further, blue-shifted apex and mineral-adjacent positions within sections perpendicular to bedding were observed, with the latter being reported here for the first time. While blue-shift at apex positions can be attributed to mechanical deformation-induced reorientation of photoselected fluorophores, the blue-shifted mineral-adjacent positions could result from strain loading and development of a plastic deformation region at the mineral contact zone with Tasmanites . A decrease in standard deviation with increasing number of measured emission maxima is well-observed, and 15 to 20 individual measurements per sample appears sufficient for low standard deviation and coefficient of variance. CLSM-derived thermal maturity parameters indicated that a moderate positive correlation of red/green quotient (Q 650/500 ; R 2 = 0.67) with solid bitumen reflectance (BR o in %) exists. For reproducible results, the determination of λ max and Q 650/500 should be conducted exclusively in the middle position at perpendicular and parallel sections of the polished whole-rock pellets, where the lowest standard deviation in measurement was observed. These results strengthen the suitability and relevance of the CLSM technique in thermal maturity studies of dispersed organic matter (DOM) and contribute to the standardization of fluorescence microspectrometry methods in organic petrology investigation.

Kentucky, Ohio, Virginia, West Virginia↗
USGS2026

Age and colony variation in Adélie penguin metapopulation vital rates: Insights from a 25-year mark–recapture study

Understanding how vital rates vary with age, life-history stage, and among populations is fundamental for predicting the demographic consequences of environmental change, especially in longer-lived species with complex life histories. These species often exhibit delayed maturity and iteroparity, making nuanced demographic insights critical for assessing their long-term viability. This study investigated age- and colony-related variation in Adélie penguin vital rates including survival, recruitment, and breeding propensity. We used mark–recapture data collected over 25 years (1996–2020) from three Adélie penguin breeding colonies that differed in population sizes and trends but comprised a metapopulation located at capes Royds, Bird, and Crozier on Ross Island, Antarctica. We used multi-state models to estimate survival and detection rates relative to reproductive state and breeding colony and estimated transition probabilities reflecting movements between reproductive states and colonies. Apparent survival varied by reproductive state, colony, and age and averaged 0.80 (SD = 0.02) at Bird, and 0.72 (SD = 0.03) and 0.73 (SD = 0.03) at Crozier and Royds, respectively, for pre-breeders age 2–7 years with strong declines in pre-breeder survival after age 8. We observed less age-related variation in survival of breeders and non-breeders, but we observed differences between colonies with lower survival for breeders (0.72 to 0.80) compared to non-breeders (0.75 to 0.82). The average probability of surviving the first 2 years after fledging ranged from 0.43 (SD = 0.14) at Royds and Crozier (0.43, SD = 0.10) to 0.55 (SD = 0.15) at Bird. Movement between colonies was highest for pre-breeders (0.00%–12.00% depending on age and colony) and lowest for breeders (<0.20%). We observed the lowest age-related recruitment rates at Royds, with recruitment at Crozier almost twice as high, and intermediate at Bird. Breeding propensity was highest at Crozier and lowest at Bird. Colony-specific variation in vital rates likely contributed differently to population trajectories, suggesting that care must be taken to extrapolate vital rate estimates across colonies even within a metapopulation. These findings also highlight the importance of considering age, life-history stage, and geographic variation when assessing population vital rates in species with complex life histories.

Frontiers in Ecology and Environment↗
USGS2024

Temporal, spatial, and chemical evolution of Quaternary high-silica rhyolites in the Mineral Mountains, Utah

The Mineral Mountains in southwestern Utah are a structurally controlled core complex at the confluence of the Colorado Plateau and the Basin and Range physiographic provinces. Aside from hosting Utah’s largest batholith, the Mineral Mountains host some of the State’s youngest high-silica rhyolites, which have been linked to a magma source that is presently being utilized as an enhanced geothermal system. The high-silica rhyolites take the form of effusive lavas and domes, and explosive products are rare. Previous K-Ar dating of these Pleistocene rhyolites placed eruptions between about 790 and 500 kilo-annum (ka) with contemporaneous basalts erupting in the valley to the east of the Mineral Mountains. Large uncertainties on these ages obscured the tempo of eruptions and thus hindered attempts to constrain the timescales of the petrogenetic processes that produced the rhyolites. In this study, we build on previous studies conducted in the 1970s and 1980s by using new geochronologic and geochemical data to investigate the temporal and spatial evolution of the youngest phase of volcanism in the Mineral Mountains. We identify two major eruptive periods, from approximately 850 to 750 ka and from approximately 590 to 480 ka. The older phase is characterized by the eruption of several basaltic lavas, two obsidian flows, and a series of coalescing porphyritic rhyolite domes. The younger phase included the eruption of six evolved high-silica rhyolite domes and one pyroclastic deposit, followed by the eruption of trachyandesite in the adjacent valley to the east. Whole-rock geochemical data indicate that the rhyolites can be divided into three chemical groups, with more evolved compositions erupting through time. The youngest rhyolites along the range crest have the lowest total iron and TiO 2 concentrations and the highest incompatible element concentrations, indicative of increasing differentiation with time and elevation. Improved precision on the eruption ages indicates a recurrence interval of approximately 20 thousand years. The eruptive flux for both periods of rhyolitic volcanism is about 0.01 cubic kilometers per thousand years, which is less than the magma resurgence flux rates for syn-caldera and post-caldera eruptions of the Valles Caldera and Yellowstone Caldera volcanic systems. Collectively, these geochemical, geochronological, and volumetric data may facilitate a better understanding of heat flux and the longevity of magmatic sources related to geothermal resources in similar small-volume, silicic systems.

Utah↗
USGS2025

Estimation of dynamic geologic CO2 storage resources in the Illinois Basin, including effects of brine extraction, anisotropy, and hydrogeologic heterogeneity

Since the vast majority of carbon dioxide (CO 2 ) storage resources in the United States are in deep saline aquifers, optimizing the use of these saline storage resources could be crucial for efficient development of geologic CO 2 storage (GCS) resources and basin- or larger-scale deployment of GCS in the country. Maximum CO 2 injection rates can be enhanced by extracting brine from the CO 2 storage unit. However, disposal of the extracted brine is both a technological and economic challenge. The lowest-cost option would likely be reinjection of the extracted brine into another formation above or below the CO 2 storage unit. Therefore, it is important to estimate brine injectivity as it will constrain the potential to increase CO 2 injectivity at an injection site that has access to multiple geologic storage units where either CO 2 or brine can be injected. Using a simulation-optimization framework, coupled with a non-isothermal, multiphase CO 2 -water-salt equation-of-state module, we developed a computationally efficient method for evaluating optimization of simultaneous CO 2 injection, brine extraction, and brine (re)injection at hypothetical injection sites deployed across a geologic basin. The Illinois basin is ideal for testing our methodology because it contains multiple geologic storage units with seals in between them to isolate injection of CO 2 in one unit from interfering with the injection of either brine or CO 2 in another unit above or below it. In addition, we investigated the relative effects of variation in key geologic parameters as well as two reservoir structures (hydrogeologic heterogeneity/anisotropy and homogeneity/isotropy) on CO 2 injectivities and enhancement of CO 2 injectivity through extracting brine. Results suggest that permeability, depth, and especially thickness of the storage unit could be the most influential parameters determining CO 2 injectivity. They also suggest that only injecting CO 2 into the storage unit with the greatest injectivity, enhancing that unit’s injectivity by extracting brine, and disposing of the produced brine in other suitable units could maximize total CO 2 injectivity in limited regions of the basin. At the majority of simulated injection sites, however, we found that injecting CO 2 into all of the accessible and suitable storage units was more likely to maximize the CO 2 storage resource.

Illinois, Kentucky, Indiana↗
USGS2025

Reaction kinetics and accelerant effects of sulfides in early mature hydrocarbon generation using hydrous pyrolysis

Hydrocarbon generation in organic-rich sediments is influenced by the molecular organic composition and relative abundance of associated minerals. Certain mineral-derived elements act as catalysts and reaction intermediaries, facilitating early-stage hydrocarbon formation in potential source rocks. This study investigated the role of sulfur contributed from pyrite as an accelerant in thermal reaction, focusing on its effects on early maturation and consequent hydrocarbon generation from gilsonite (low-sulfur solid petroleum). Hydrous pyrolysis (HP) experiments were conducted on mixtures of gilsonite and pyrite in varying ratios (1:0.1, 1:0.5, 1:1, 1:2, and 1:10 w/w gilsonite:pyrite) at 320, 350, and 370 °C for 72 h. Untreated and thermally altered residues were analyzed using solid bitumen reflectance (BR o , %), total organic carbon (TOC) content, programmed temperature pyrolysis, scanning electron microscopy with energy-dispersive spectroscopy (SEM-EDS), and X-ray diffraction (XRD) to evaluate the potential accelerant role of pyritic sulfur in hydrocarbon formation. The results show HP residues at 320 and 350 °C with greater pyrite concentrations had higher BR o , while reflectance values were similar in the 370 °C residues, regardless of pyrite concentration, suggesting enhanced reaction at lower thermal conditions. Increasing pyrite content systematically decreased hydrogen index (HI) values while increasing the transformation ratio (TR) and production index (PI), indicating enhanced conversion of organic matter to hydrocarbons with increasing pyrite concentrations. Gas yields increased with pyrite addition, particularly at 350 °C, confirming secondary cracking effects. However, gas production stabilized or declined at higher pyrite loadings (1:10), suggesting alternative reaction pathways such as coke formation. Our data indicate the presence of pyrite lowers the activation energy for thermal cracking, shifting peak experimental hydrocarbon generation temperatures downward by 20–30 °C, with the most pronounced accelerant effects observed at moderate pyrite concentrations (1:0.5 and 1:1). The thermodynamic framework reveals that pyrite stability is influenced by experimental conditions, with pyrrhotite formation favored in the presence of gilsonite due to reduced oxygen fugacity. Pyrite transformation to pyrrhotite, as observed through XRD, SEM-EDS, and predicted by thermodynamic data, further supports the accelerant role of S, as pyrrhotite exhibits a higher hydrogen transfer potential, promoting early oil generation. These findings highlight the importance of pyrite in modulating hydrocarbon generation pathways in organic-rich systems.

Journal of Analytical and Applied Pyrolysis↗
USGS1973

Geohydrology of the Cross-Florida Barge Canal area, with special reference to the Ocala vicinity

The Cross-Florida Barge Canal route commences at Palatka on the St. Johns River, about 75 miles upstream from the Atlantic Ocean, and extends 110 miles southwestward across Peninsular Florida into deep water in the Gulf of Mexico near Yankeetown. The canal will be equipped with five locks, each 600 feet long and 84 feet wide, and the channel will be a minimum of 12 feet deep and 150 feet wide. From near Ocala northeastward, the canal channel will replace much of the natural channel of the Oklawaha River, and will be excavated into beds of the so-called shallow sand aquifer of Miocene age and younger, which overlies limestone of the Floridan aquifer. Westward from Ocala, most of the canal will be excavated below the potentiometric surface into limestone and dolomite of the Floridan aquifer. Water levels of Rodman, Eureka, and Inglis Pools will be controlled by dams and spillways with the limited exchange of water between the pools and the aquifers. The water levels in the Summit Pools will fluctuate with the natural changes in the ground-water level of the Floridan aquifer, although the stage of the pool will be controlled partly by the stage held in the Eureka Pool. A dynamic inflow-outflow relationship will exist between the Summit Pool and the Floridan aquifer. The Floridan aquifer in the canal area is 1,000 to 1,200 feet thick and consists of limestone and dolomite of middle Eocene Miocene age, including from older to younger, the Lake City, Avon Park, and Ocala limestones plus permeable sandy, dolomitic limestone in the lower part of the Hawthorn Formation. It is possible that most of the flow to the two major springs in the area occurs in the upper 100 feet or so of the aquifer in the Ocala Limestone. The aquifer is underlain by the Oldsmar limestone of early Eocene age and is overlain by sand, clayey sand, clay and shell beds of Miocene through Holocene age, in thickness from a few feet to 300 feet. The permeable beds overlying the Floridan aquifer constitute the shallow aquifer, while the poorly permeable ones act as confining beds where the Floridan aquifer is under artesian conditions. A north-south line drawn separating the head of Silver Springs on the west from the Oklawaha River on the east marks the approximate western limit of a continuous blanket of materials of Miocene-Pliocene(?) age covering the rocks of the Floridan aquifer. East of the line, much of the aquifer is under artesian conditions, particularly in the Oklawaha River valley, although in some areas east of the valley, direct recharge through thick permeable Miocene-Pliocene(?) sands occurs. West of the line, only scattered remnants of a once continuous Miocene-Pliocene(?) cover remain. Lack of the cover is a result of erosion on the crest and flank of the Ocala Uplift, a broad northwest-southeast trending anticlinal upwarp, the axis of which is crossed by the canal route in the Dunnellon area. Over most of this area the Floridan aquifer is unconfined and receives direct recharge through a cover of a few tens of feet of sand and clayey sand of Quaternary age. Tensional stresses during the structural evolution of the Ocala Uplift produced an intersecting system of fractures and normal faults in rocks of the Floridan aquifer. The fractures and faults are important controls for orientation of solution channels and, therefore, for development of ground-water circulation patterns. When the system surface streams, which once drained the Barge Canal area, eroded the poorly permeable Miocene-Pliocene(?) cover from the flanks of the Ocala Uplift, surface runoff was reduced and precipitation began to directly infiltrate the underlying limestones. Now only principal streams remain, such as the Oklawaha and Withlacoochee Rivers and a few short tributaries, while one of the most highly developed subsurface drainage systems in the world has evolved in cavernous limestone of the Floridan aquifer. Two of the larger freshwater springs in the world now discharge from the Floridan aquifer in the canal area. Silver Springs near Ocala discharges an average 531 mgd (million gallons per day) down the 4-mile long Silver River, which flows on poorly permeable beds to the Oklawaha River. Rainbow Springs near Dunnellon discharges on average 468 mgd from numerous orifices in the bed of the 5-mile-long Rainbow River, which flows into the Withlacoochee River. The heads of the springs have migrated to their present positions partly because of a tendency of ground-water levels to decline as permeability in the aquifer is increased due to removal of limestone by solution, and because of mechanical erosion of the limestone in the vicinity of the spring heads. Also, points of principal spring discharges have shifted in the past due to changes in ground-water levels in response to changes in sea level. The subsurface drainage system is continuing to evolve today, as evidenced in part by frequent occurrence of new sinkholes and by the presence of significant amounts of calcium bicarbonate in the spring waters. Rodman Pool, at the east end of the canal, is separated from the Floridan aquifer by poorly permeable materials. The pool's operating water level will be only a few feet above the potentiometric surface at the downstream end, and at or slightly below the potentiometric surface at the upstream end. Little exchange of water between the Rodman Pool and the Floridan aquifer is expected. Eureka Pool, just upstream from Rodman Pool, will also be separated from the Floridan aquifer by poorly permeable beds. However, the stage of the pool will be about 15 feet higher than the natural potentiometric surface at the pool's downstream end, and some seepage into the Floridan aquifer is anticipated through faults and leaky parts of the poorly permeable beds, with a consequent rise in ground-water levels in areas adjacent to the lower end of the pool. Possibilities for particulate contamination of the aquifer will tend to be minimized because of the filtering capacity of the materials through which water must pass to reach the aquifer, although the natural filter will not preclude movement into the aquifer of contaminants which might become dissolved in the pool waters. No significant interchange of water between the pool and the aquifer is expected at the upstream end of the Eureka Pool. Present construction plans indicate an operating stage for Eureka Pool which will range between 38 and 40 feet above mean sea level, although it is possible to dredge the pool deep enough to permit a range in stage of 36 to 40 feet. A backwater effect extending up Silver River from the Eureka Pool is expected to regulate the stage at the head of Silver Springs between 39 and 44 feet above mean sea level if the pool ranges between 36 and 40 feet. If Eureka Pool ranges only between 38 and 40 feet, the range of stage at the head of the springs should be about 41 to 44 feet above mean sea level. From Inglis Lock west, the canal will have direct connection with the Gulf, and canal stage will fluctuate with the Gulf tide. Since the canal stage will be slightly lower than the adjacent ground-water levels along much of the reach, there will be some ground-water inflow to the canal. No significant changes in the existing ground-water regime are expected in the vicinity of Inglis Pool, the first step up in the canal east of the Gulf. Existing ground-water and surface-water levels in the area will not change appreciable, and the natural stage and flow of Rainbow Springs, which will flow by way of Rainbow River into Inglis Pool, should not be affected by canal operations. A possible adverse effect of the canal on the Inglis Pool area could result if sea water is locked up from the Gulf through Inglis Lock. However, the high step of 25 feet at the lock, flushing action of continuous flow from Inglis Pool to the lower reaches of the Withlacoochee River, and use of possible preventive locking procedures should minimize the problem. The potential for adverse effects on the ground-water regime is greatest in the area of the Summit Pool. Through most of the length of the pool, the canal channel will be excavated into limestone of the Floridan aquifer to depths of 12 to 27 feet below the potentiometric surface. Changes that will take place in the ground-water flow system in the Silver Springs drainage area, once the canal is completed, were estimated by flow-net analysis. Variation in aquifer transmissivity was determined by calculating transmissivity in 25 different flow cells surrounding Silver Springs. Transmissivity in the 25 cells averages about 15,600,000 gpd/ft [2,090,000 ft2/day (feet squared)], but transmissivity in the six cells through which the Summit Pool passes ranges from 9,000,000 to 44,000,000 gpd/ft (1,210,000 to 5,900,000 ft2/day). Transmissivity was used to compute static stage of the Summit Pool under given ground-water level conditions. Had the canal existed in May 1968 and had the stage of Eureka Pool been held at 36 feet at the time, the static stage in Summit Pool would have been about 42.1 feet above mean sea level. Thus, a conceptual model of the changes in the potentiometric surface wrought by the finished canal was drawn, and zones of ground-water inflow and outflow were delineated. Most outflow from the Summit Pool to the aquifer should be limited to one 4-mile-long zone along the north side of the pool, about 5 miles south of Silver Springs. It is estimated that a water volume equivalent to about 8 percent of the daily flow of Silver Springs will enter the Summit Pool each day from the southern one-third of the Silver Springs drainage area. A like amount will reenter the aquifer at the main zone of outflow and move toward Silver Springs at an estimated average velocity of about 200 feet per day, if something close to the natural static stage of the pool is maintained by return pumpage of the lockage losses. At a velocity of 200 feet per day, water from the Summit Pool would discharge at Silver Springs about 140 days later. However, any estimate of velocity in the highly cavernous limestone aquifer in the area should be used with caution, because difficult to measure changes in porosity and thickness of major zones of flow may cause large variations in velocity. If all lockage losses are returned to the Summit Pool by pumping from Eureka Pool, no net loss from the Silver Springs drainage area, except for some evaporation from the water surface in the canal and possible leakage around locks, will result from canal operations. The zone of outflow from the Summit Pool to the aquifer will be in a natural potentiometric trough, and the zone of inflow will be in a potentiometric ridge area. The equilibrium water level in the pool will tend to be about 1 foot higher than the altitude of the lowest level in the pre-canal potentiometric trough, and about 2 feet lower than the highest level on the pre-canal potentiometric ridge. West of the Silver Springs drainage area just east of Dunnellon Lock, in the area of a local potentiometric high, the water level in the Summit Pool is expected to be about 15 feet below the natural potentiometric surface. In most areas 2 to 3 miles away from the Summit Pool, effects of the canal on the natural potentiometric surface should be slight. The stage of the Summit Pool, judging from the 36-year record for ground-water level changes and the anticipated indirect effect of the controlled stage in Eureka Pool, should have a maximum of about 10.5 feet with a maximum water level of about 51.5 feet above mean sea level and a minimum of about 41.0 feet above mean sea level. Of particular importance in the Summit Pool is an implementation of well planned construction and operational procedures designed to minimize risks of ground-water contamination.

Florida↗
USGS2025

Assimilation of reduced carbon triggers platinum alloy saturation in mafic and ultramafic magmas

It is generally observed that magmatic sulfide ores have higher ratios of Pd/Pt than the mantle-like values of their parental magmas. This discrepancy has defied simple explanation because the partitioning behavior of both elements between sulfide and silicate liquids is very similar. Assimilation of sulfur- and carbon-rich country rocks by mafic and ultramafic magmas is considered a critical, if not essential, step in the formation of magmatic base metal sulfide deposits. Although there is general consensus that the assimilation of external sulfur and carbon promotes sulfide saturation, the effect of carbon assimilation on the solubilities of platinum-group elements in natural S-bearing silicate melt has been overlooked. In this study, we investigate the variations of platinum and palladium solubilities during assimilation of graphite and methane through thermodynamic modeling, in comparison with data from an array of highly distinctive magmatic sulfide ore systems representing ages from Archean to Paleozoic, melt compositions from komatiite to basalt, and magmatic settings including lavas, hypabyssal intrusions, plutonic continental arc roots, and plutonic layered intrusions, namely: Raglan, Norilsk-Talnakh, Lac des Iles, and the J-M Reef of the Stillwater Complex. We model assimilation-fractional crystallization processes to estimate the reduction of oxygen fugacity ( ⁠fO2⁠ ) of the melt due to incorporation of graphite and methane. The simulations show that although Pd remains highly soluble during the progressive assimilation of reduced carbon, Pt solubility decreases significantly as the silicate melt becomes increasingly reduced. With less than 8 % of sediment assimilation, Pt alloy may saturate and then deviate from sulfide-undersaturated silicate melts, concomitantly increasing the Pd/Pt value of the remaining melts of the Raglan and Norilsk-Talnakh systems. For the Lac des Iles and Stillwater systems, a higher extent of assimilation is needed to reach Pt saturation because of the relatively carbon-poor nature of the lower crustal rocks. The assimilation of methane volatiles is shown to be more effective than graphite assimilation, and it provides a pathway to Pt alloy fractionation in the absence of detectable amounts of bulk host-rock assimilation. High Pd/Pt values have been documented in many world-class magmatic sulfide deposits whose parental magmas have demonstrably experienced crustal contamination. Our model suggests that although anomalous Pd/Pt values may be explained by other mechanisms such as incongruent melting of preexisting sulfide or differences in the diffusivities of the metals within achieving equilibration, the assimilation of graphite or methane may play an important role in the global occurrence of magmatic sulfide ores with elevated Pd/Pt values.

Economic Geology↗
USGS2026

Small-volume tephra deposits of the May 1924 explosions from Halemaʻumaʻu, Kīlauea volcano, and their origin

More than 50 explosive eruptions occurred from Halemaʻumaʻu at Kīlauea volcano over 17 days from May 11 to 27, 1924. Ballistics weighing as much as 14,000 kg were ejected and most landed within 2 km of the vent. Fine ash made up a major component of the tephra and was dispersed tens of kilometers downwind. Draining of the Halemaʻumaʻu lava lake occurred in late February 1924, with the crater floor eventually subsiding by a further ∼70 m (to ∼180 m below the crater rim) by the time the first explosions took place during the night of May 10–11. The largest explosions occurred on May 17–18 and smaller explosions continued until May 27, at which point Halemaʻumaʻu had more than doubled in width and depth. The explosions generated plumes reaching up to ∼10 km high with ballistics ejected up to 2 km from the crater. Almost 100 years later, we investigate and characterize the preserved tephra deposits within ∼3 km of the 1924 crater rim. Grain size and shape analyses were performed on 202 samples collected from 34 tephra profiles using dynamic image analysis, with a subset of layers from nine tephra profiles used for componentry (200 grains per layer in the 0.5–1 mm size fraction). Additionally, we characterize the average diameters (using the five largest clasts) at 216 locations and measure the average diameters of 2291 ballistics (largest per ∼100 m 2 area). Physical descriptions from fieldwork and grain size distributions were used to subdivide the tephra layers into five lithofacies: coarse homogeneous, fine homogenous, red ash, accretionary lapilli-bearing, and finely laminated. Grain size versus shape data show a range of values that demonstrate most grains are dense, smooth, and equant, in alignment with lithic clasts dominating the tephra componentry. The fine grained and accretionary lapilli-bearing nature of some of these lithofacies confirms that water influenced the style of the explosions. However, we also note juvenile clasts within many of the tephra layers, indicating that many of the layers were formed during phreatomagmatic explosions (sensu stricto), despite the eruptive mechanism being dominantly phreatic. Juvenile clasts are more abundant higher in the tephra profiles, suggesting that juvenile magma was more involved later in the explosive sequence. Thermal and hydrologic modeling indicate that groundwater inflow into a short-lived, small-diameter volcanic conduit (10-m to 120-m-diameter used for modeling) during the 78–85 days preceding the first explosion provides a physically plausible mechanism for this eruptive sequence.

Hawaii↗
USGS2025

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

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

Baja California Sur↗
USGS2025

Pre-eruptive characteristics of “suspect” silicic magmas in Carlin-type Au-forming systems

World-class Carlin-type Au deposits hosted in sedimentary rock were formed when profuse Eocene silicic magmatism swept across northern Nevada in response to arc migration. Carlin-type Au deposits formed along with porphyry/skarn Cu-Mo-W-Au deposits, epithermal Ag-Au deposits, and distal disseminated Ag-Au deposits. But unlike these other Au-bearing deposits that have clear associations with igneous intrusions, Carlin-type ore deposits appear to have formed distant from concealed plutons, and their origin remains controversial. Despite decades of abundant geophysical, geochronological, and geochemical studies suggesting the involvement of magmas, concrete evidence for magmatic involvement is still lacking. Consequently, the involvement of contemporaneous igneous systems remains inferred based on age, proximity, and variable isotopic, geochemical, and geophysical clues. A recent synthesis of deposit models postulates that Carlin-type Au deposits are intrusion-related, but that the causative magmas reside deeper (∼6–12 km) than in typical porphyry and peripheral systems (∼3–5 km), meaning that Carlin-type deposits are perhaps more distal expressions of igneous intrusions. We investigate a collection of “suspect” magmatic systems over a ∼7 m.y. timespan (∼41–34 Ma) that are contemporaneous with and near known Carlin-type ore deposits. We report results of a multifaceted array of in situ geochemical analyses (FTIR, EMP, SHRIMP-RG, LA-ICP-MS) of quartz-hosted melt inclusions, biotite, and quartz to better characterize the pre-eruptive characteristics of these magmas. We also report results of thermobarometry and thermodynamic phase equilibria modeling to help place constraints on magmatic reservoir depths and processes. Rather than a single “flavor” of silicic magma, we observe a surprisingly broad compositional spectrum of rhyolites, with one end of the spectrum exhibiting more arc-like (I-type) characteristics and the other end displaying more post-subduction, thick-crust extensional (A-type) characteristics. This broad compositional spectrum suggests a more complex picture of silicic crustal magmatism operating over a narrow span of time during slab rollback. Despite this spectrum, magmatic systems in this study are consistently ferroan and generally peraluminous, which we interpret as an expression of the relatively elevated geotherm at the time and incorporation of variable amounts of highly peraluminous metasedimentary crustal components. The silicic magma spectrum encompasses a range of mineralization associations, including subduction-related Cu-Mo-W-Au-Ag and post-subduction, thick-crust extensional rare-metal Mo-Sn-W-F-Be-Ag-Au, consistent with the prolific and diverse array of ore deposits that formed during this time. Carlin-type Au deposition appears to be associated with nearly the entire magmatic spectrum. This apparent indifference to silicic magma “flavor” would seem to imply that if magmas are involved in Carlin-type Au deposit genesis, they perhaps do not need to be compositionally specialized and/or possibly are only relevant as heat sources driving circulation to remobilize and redistribute metals.

Nevada↗
USGS2025

Chemostratigraphy of the Cretaceous Hue Shale in Arctic Alaska: Exploring paleoceanographic controls on trace element enrichment, organic matter accumulation, and source-rock evolution

We document chemostratigraphy in an outcrop of late Albian to early Campanian (∼103–82 Ma) marine source rocks to better understand paleoenvironmental controls on trace element (TE) enrichment and organic matter accumulation in the distal Colville foreland basin of Arctic Alaska and how those drivers are linked to arc volcanism and successions of Cretaceous oceanographic and climatic biogeochemical events. This unique, 113-m-thick section of Cretaceous Hue Shale deposited during a series of previously undocumented Arctic Cretaceous oceanic anoxic events (Lease et al., 2024) is the only known exposure of thermally immature (0.48–0.52% R o , random vitrinite reflectance) Hue Shale in Arctic Alaska. Strata comprise mainly clay-rich mudstone with elevated total organic carbon (TOC) and hydrogen index values reaching 26.3 wt% (mean = 7.5 wt%) and 689 mg hydrocarbon (HC)/g TOC (mean = 385 mg HC/g TOC), respectively. Maceral composition consists predominantly of fluorescent amorphous organic matter, with abundant brightly fluorescent alginite, including Tasmanites , acritarchs, and Leiosphaeridia . Discrete layers of volcanic ash (preserved as bentonite) are present throughout the section and provide quantitative age control based on U–Pb dates. Chemostratigraphic trends are investigated to advance our understanding of local oceanographic conditions and controls on stratigraphic and temporal heterogeneity of Brookian source rocks. Concurrent sedimentary enrichment in Mo, U, V, Pb, and Cu across the Albian–Cenomanian boundary of the exposed basal gamma-ray zone, may reflect anoxic to euxinic benthic redox conditions favoring organic matter accumulation and preservation. Fluctuating degrees of anoxia-euxinia are inferred throughout the overlying Hue Shale succession, reflected by varying patterns of TE enrichment and TE–TOC covariation. Elevated C org /P molar ratios (>250) across most of the section, with several values exceeding 690, signify that enhanced biological productivity is sustained throughout deposition. Enhanced productivity, recorded by both C org /P and excess Ba, also parallels increases in source rock richness (elevated TOC and S 2 values) during the late Albian–early Cenomanian and late Cenomanian – Turonian. Enhanced productivity and variations in oceanic circulation/stratification likely both drove changes in benthic redox conditions that favored organic carbon accumulation and preservation. Increased continental arc volcanism (e.g., Okhotsk–Chukotka volcanic belt) and High Arctic Large Igneous Province magmatic eruptions throughout the Cretaceous, inferred to have influenced nutrient cycling and local aqueous nutrient availability, also have been invoked as potential drivers of organic carbon burial and source-rock development across the sedimentary sequence. Results presented here document the organic-rich and oil-prone source-rock quality of the Hue Shale in the distal part of the Colville foreland basin and bolster the potential for a Cretaceous petroleum system beneath the eastern North Slope.

Alaska↗
USGS2026

Determining Volcanic Risk in Auckland (DEVORA) Research Programme—A transdisciplinary approach to address the challenge of distributed volcanism in an urban environment

The Determining Volcanic Risk in Auckland (DEVORA) Research Programme was launched in 2008 to address the challenges associated with monogenetic volcanism in an urban setting and to enhance volcanic risk management in Tāmaki Makaurau Auckland in Aotearoa New Zealand. It is a multi-agency, increasingly transdisciplinary (defined here as research that transcends traditional disciplinary boundaries by integrating diverse types of knowledge, perspectives, and methods from academic and non-academic participants to create novel solutions to complex problems), and collaborative research program jointly led by Waipapa Taumata Rau University of Auckland and Earth Sciences New Zealand (ESNZ; formerly GNS Science), with core funding from Natural Hazards Commission Toka Tū Ake (NHC; formerly the Earthquake Commission, EQC) and Te Kaunihera o Tāmaki Makaurau Auckland Council (AC). The primary research focus of DEVORA is to investigate the geologic history, volcanic hazards, and risk posed by the basaltic intraplate Auckland Volcanic Field. Disruption from ash fall and gas from other Aotearoa New Zealand volcanoes is also considered. DEVORA’s work to explore exposure and vulnerability in Tāmaki Makaurau Auckland is also useful for assessing risks from other non-volcanic natural hazards, such as seismic and tsunami hazards. The greater Tāmaki Makaurau Auckland region has an ethnically and socio-economically diverse population of approximately 1.7 million, representing about one-third of the Aotearoa New Zealand population, and hosts critical infrastructure of national significance. The size and nature of the populace, consequential economic base, and important infrastructure within Tāmaki Makaurau Auckland mean that the effects of a volcanic eruption would be felt nationally, including through the disruption of air travel to Aotearoa New Zealand. The hazards from such an eruption could potentially affect hundreds of thousands of people, businesses, and lifelines (critical infrastructure). A considerable challenge for emergency and risk managers is the monogenetic nature of the volcanic field. It is not known where or when the next eruption will occur, how much warning we may get before an eruption, nor how an eruption and its effects might unfold. In this contribution, we highlight the concept and collaborative intent of the DEVORA Programme and show how it has evolved over the 16 years since its inception. We describe how DEVORA has unified more than 100 researchers (including more than 50 graduate students) and numerous stakeholders to address key issues facing Tāmaki Makaurau Auckland and describe how research findings are being implemented into policy and communicated to stakeholder agencies and the public. We also illustrate the broader influence of the DEVORA Programme and provide some learnings that might benefit others embarking on similar integrated projects, especially those focused on distributed volcanism in and near populated areas.

Auckland↗
USGS2025

Tracking diagenetic alteration of magnetic susceptibility in thrust ridge and slope basin sediments of the Cascadia margin (ODP Sites 1249 and 1252; IODP Site U1325)

We investigated sediment core records from the Cascadia Margin (Ocean Drilling Program Sites 1249 and 1252 at Hydrate Ridge; Integrated Ocean Drilling Program Site U1325 offshore Vancouver Island) using a Zr/Rb heavy mineral proxy from X-ray fluorescence (XRF) core scanning to identify intervals of primary detrital magnetic susceptibility (κ) and predict intervals where diagenesis caused magnetite dissolution by hydrogen sulfide. We also measured total sulfur (TS) content, grain size distributions, total organic carbon (TOC) content, and the magnetic mineral assemblage to further constrain the role of diagenesis on κ. Understanding how κ can be used to better characterize the varied effects of detrital and diagenetic signals in marine settings is important for understanding biogeochemical cycling and records of paleoenvironmental change. The upper 100 m of slope basin Site 1252 contains multiple intervals (> 90 m total) of decreased κ correlated with elevated TS content, consistent with dissolution of magnetite and precipitation of pyrite, iron monosulfides, and/or elemental sulfur. Similarly at the other slope basin site, Site U1325, κ is lower and TS is elevated in the interval between 24 and 51 mbsf, due to sulfide formation. At both slope basin sites, these low κ intervals correspond with high TOC, suggesting the possibility that organoclastic sulfate reduction (OSR) is likely a major driver of diagenetic alteration of κ at these sites. High TS:TOC ratios at Site U1325 suggests anaerobic oxidation of methane (AOM) during sulfate-methane transition zone (SMTZ) migration may have contributed to alteration of κ. In contrast, within the upper 90 m of Site 1249, a methane seep site at the summit of Hydrate Ridge, κ is almost entirely altered by diagenetic processes, with much of the low κ explained by a high degree of iron sulfide formation, while some intervals are affected by precipitation of magnetic iron sulfides that maintain or even increase κ. The presence of abundant methane seepage and gas hydrate as well as chemosynthetic seafloor fauna at this site, suggests that sulfide is released to the water column and AOM, rather than OSR, drives diagenetic alteration of κ at this site. Overall, the slope basin sites show episodic variation of κ that is influenced by TOC content, likely driven by changes in marine primary productivity and sedimentation rate, while the seep site shows consistently altered κ with lower TS content and no correlation with TOC. Methane seep environments likely experience loss of hydrogen sulfide to the water column and oxidation of hydrogen sulfide by seafloor seep fauna, which limits the amount of solid phase sulfur (pyrite, iron monosulfides, elemental sulfur) that can be precipitated within the sediments. In contrast, the migration of a buried SMTZ at slope basin sites results in enhanced sulfur precipitation within the sediments. This integrated magnetic and geochemical approach reveals the diagenetic production pathway and residence time of sulfide with the sediment column ultimately controls the style and degree of diagenetic loss of κ in marine sediments. This approach works best in environments with unaltered reference intervals/sites, consistent magnetic mineralogy, and clay-to-silt grain sizes.

Cascadia margin, Pacific Ocean↗
USGS1976

Helium detection as a guide for uranium exploration

Helium, a byproduct of radioactive decay, may prove to be a valuable indicator of the presence and distribution of uranium deposits. Recent technological advances permit the development of instrumentation not previously adapted for this purpose. Commercially available equipment can provide high sensitivity at low cost and allow reasonable mobility for field use. A truck-mounted mass-spectrometer, tuned for helium-4, permits immediate adjustment or modification of sampling patterns in response to accumulating data. The inlet system of the spectrometer has been designed to allow flexibility in gas analyses from various sample types--soil gas, atmosphere, or gases in water. Sensitivity of the instrument is better than 50 parts of helium per 10 9 parts of gas. Replicate samples and standards can be analyzed in only 3 minutes. Helium in soil gas is being used for the initial evaluation of the technique. A hollow probe, as much as 2 m long, is driven into the ground; and a 10-cc syringe is used to purge the probe and extract the gas. This sample is then injected into the inlet system and introduced into the spectrometer at constant pressure. The output signal is displayed on a chart recorder. Investigations can be performed on a qualitative, relative basis, or on a quantitative basis by comparison to calibrated helium standards. Preliminary field testing includes studies of the responses to variations in wind speed, temperature, barometric pressure, moisture, and sampling depth over extended time periods, as well as studies of geologic controls on the helium content in soil gas. Surveys over known uranium occurrences reveal some anomalous helium distributions. This report describes an updated technological approach to an old idea: that of using helium, a byproduct of uranium radioactive decay, as an exploration tool for uranium. Helium is the sixth most abundant gas found in the earth's atmosphere (Table 1). It is an inert gas and very mobile, mixing rapidly. When the alpha particles from the radioactive decay of uranium and thorium pick up two electrons, atoms of the isotope helium-4 are formed. Several helium-4 atoms are produced from each decay series of U-238, U-235, and Th-232 (fig. 1). The half-lives of the parent isotopes are5also shown in figure 1. One gram of uranium will produce ~10 5 atoms of helium-4 per second, and one gram of thorium produces ~2.5 x 10 4 atoms of helium-4 per second (fig. 2). Calculations of the helium produced from the crust and mantle reveal that 1,125 x 10 30 atoms per year are produced, but only 7 x 10 30 atoms per year escape from the earth (Damon and Kulp, 1958). More is being produced than is being lost; in fact, the total atmospheric content of helium could be produced in only 2 million years (MacDonald, 1963)! However, all the crustal and mantle helium-4 does not degas into the atmosphere as it is formed. It is trapped in crystal lattices and in pore spaces within the earth. What this means is that there is excess helium-4 in the earth--excess, that is, compared to the atmospheric concentration in equilibrium with the helium-4 escape rate into space and the helium-4 flux from the crust and mantle (Nicolet, 1957). Near-surface pockets of high helium-4 concentrations are known, and some natural gas fields are so enriched that they are the source for commercial production of helium. Work by Roberts and others (1975) has shown high helium concentration in the soil gas associated with geothermal areas; and work by Clark and Kugler (1973), Dyck (1975), and Goldak (1974), for example, has noted high helium concentrations in soils and waters in the vicinity of uranium deposits. These latter studies certainly indicate the potential of helium detection for use as an exploration tool for uranium.

Open-File Report↗
USGS2025

Phase 1 technical implementation plan for the expansion of the ShakeAlert earthquake early warning system to Alaska

Executive Summary The conference report accompanying the fiscal year (FY) 2022 Consolidated Appropriations Act (Public Law 117–103) for the U.S. Department of the Interior and related agencies directed the U.S. Geological Survey (USGS) to “work with the State of Alaska to develop an implementation plan to be completed within two years in order to put ShakeAlert/Earthquake Early Warning in Alaska” (p. 29). Congress included $1 million in the FY 2022 appropriation to conduct this effort. The USGS Earthquake Hazards Program, along with partner organizations, has developed the ShakeAlert earthquake early warning (EEW) system for the West Coast, which currently operates in California, Oregon, and Washington. The purpose of the system and its alert delivery partners is to reduce the impact of earthquakes and save lives and property by delivering ShakeAlert-powered alerts that are transmitted to the public via mass notification technologies, and by providing more detailed data streams to institutional users and commercial service providers to trigger automated, user-specific, protective actions. ShakeAlert was designed in such a way that it could be expanded to other U.S. regions with high earthquake risk, after the build-out of seismic and geodetic networks to support ShakeAlert in a specified region is completed and the necessary funding is secured for long-term operation and maintenance. When an earthquake occurs, seismic waves radiate from the rupturing fault like waves on a pond. It is these waves that people feel as earthquake shaking and that can cause damage to structures. Using networks of ground-motion sensors and sophisticated computer algorithms, ShakeAlert can detect an earthquake seconds after it begins, calculate its location and magnitude, and estimate the resulting intensity of shaking. Early warnings of impending shaking are then sent to people and systems that may experience damaging shaking, allowing them to take appropriate protective actions. Depending on the user’s distance from the earthquake, alerts may be delivered before, during, or after the arrival of strong shaking. There will almost always be a region near the earthquake epicenter where alerts arrive after damaging shaking has begun. The ShakeAlert system updates its ground-motion estimates as an earthquake grows larger. In response to the FY 2022 congressional direction, the USGS worked with the State of Alaska to devise this implementation plan for ShakeAlert expansion to Alaska. The USGS engaged with the Alaska Division of Homeland Security and Emergency Management (DHS&EM) and the Alaska Division of Geological and Geophysical Surveys (DGGS). A cooperative agreement was awarded to the Alaska Earthquake Center (AEC) at the University of Alaska Fairbanks (UAF) for their contributions to the plan and their work coordinating with other networks in Alaska. The USGS engaged with the Alaska Seismic Hazards Safety Commission (ASHSC) throughout the process. The USGS also held a series of Alaska stakeholder engagements. The process of developing the implementation plan was facilitated by contracted staff from Corner Alliance, which is a government consulting firm. This implementation plan describes the details and estimates the costs for a Phase 1 expansion of the ShakeAlert system to Alaska. A geographically limited Phase 1 goal was chosen that covers the highest risk and most populated areas of Alaska. The areas proposed encompass the State’s main population centers and 90 percent of the State’s population. This Phase 1 design is considered very challenging and ambitious from the viewpoint of network operators. The lessons learned if this plan is implemented could be used to consider subsequent phases to expand EEW beyond Phase 1 in Alaska in the future. ShakeAlert is built on the foundation of the sensor networks and data processing infrastructure of the USGS-led Advanced National Seismic System (ANSS). This implementation plan calls for a total of 450 high-quality, real-time EEW-capable ANSS seismic stations in Alaska: 270 new stations, 160 upgraded stations, and 20 existing stations. These seismic station numbers are based on a station spacing of 10 kilometers (km) in urban areas, 20 km in seismic source areas that endanger population centers, and 40 km in other areas. The associated costs also include support for some EEW-capable global navigation satellite system (GNSS) stations, with a focus on improving warnings for large subduction zone earthquakes. For effective EEW, ShakeAlert requires low-latency, high-availability, robust telemetry links to deliver continuous, real-time data from field stations to the data centers. The Alaska data processing hardware infrastructure would follow the general design for fail-safe operation that is used for the ShakeAlert system on the West Coast. The ShakeAlert architecture uses two independent layers: the production layer for earthquake processing and the alert layer to make alerting decisions and serve alerts to users. This implementation plan includes two geographically separated data centers in Alaska, each with two fully independent production and alert layers using the same system design developed for the West Coast. As of March 2024, the ShakeAlert system is at version 3.0.1, with more advanced versions in the development and testing pipeline. ShakeAlert originally used two algorithms to determine the location and magnitude of earthquakes using seismic data. A third algorithm that can calculate very large magnitudes of very large earthquakes with geodetic data was added in March 2024. ShakeAlert publishes several data and alert products to meet the needs of different users. All messages include the location of the earthquake, either as a point or a line, and its magnitude. Ground-shaking estimates are published in two forms, as ground-motion contours and a map grid. Providing adequate warning time for strong shaking (the “target threshold”) requires sending alerts at a threshold lower than that strong shaking level (the “alert threshold”). The thresholds for public alerting in Alaska would be a joint USGS and State decision. To have the greatest benefit, ShakeAlert-powered alerts would be delivered to institutional users and individuals by all practical pathways. The USGS alert layer can support thousands of institutional users and alert redistributors, but the USGS does not have the mission nor the infrastructure and expertise to perform mass notifications to the public or implement automatic actions for end users of the alerts. To meet this need, ShakeAlert recruits private sector “technology enablers” that have the necessary expertise to develop end-user implementations using EEW alerts with the goal of stimulating an EEW industry. Earthquake early warning alerts are useless if people do not know how to respond to them. Although the alert messages include instructions about what to do (drop, cover, and hold on), alerts are more effective if people have been trained in advance. Messages about ShakeAlert’s capabilities, limitations, and benefits could be integrated with existing earthquake education programs, including State-run programs. Therefore, ShakeAlert would coordinate with both public and private partners and stakeholders through various partnerships and agreements to accomplish consistent and ongoing public earthquake hazard education. The estimated capital cost of completing the computing infrastructure and sensor networks for the Phase 1 ShakeAlert expansion to Alaska is approximately $66 million in 2024 dollars. The annual operation and maintenance cost of the completed system is estimated to be $12 million per year in 2024 dollars when fully built out.

Alaska↗
USGS2026

Wind River subbasin restoration: Annual report of U.S. Geological Survey activities January 2024 through December 2024

We sampled juvenile wild Oncorhynchus mykiss (Steelhead Trout) in headwater streams of the Wind River, Washington, to characterize population attributes and investigate life-history metrics, particularly migratory patterns, and early life-stage survival. Additionally, this sampling provides an index of abundance of non-native Brook Trout Salvelinus confluentus . We used passive integrated transponder (PIT) tagging and a series of instream PIT-tag interrogation systems (PTISs) to track juveniles and adults. The Wind River subbasin is considered a wild Steelhead refuge by Washington Department of Fish and Wildlife (WDFW). No hatchery Steelhead Trout have been released in the Wind River subbasin since 1997, and hatchery adults are estimated at less than one percent of spawners in most years. Over 25 years of Steelhead Trout status and trend monitoring and research in the subbasin is contributing to understanding of population response to numerous restoration actions in the subbasin, including removal of Hemlock Dam from Trout Creek in 2009, which had an inadequate adult ladder and contributed to increased water temperatures, to near lethal levels for salmonids. Data from our study, and companion work by WDFW, are contributing to the Columbia Basin Fish and Wildlife Program (2009) Research, Monitoring, and Evaluation (RM&E) Strategy of Fish Population Status Monitoring. Specifically, this work addresses the sub-strategies of 1) Assessing the Status and Trends of Diversity of Natural Origin Fish Populations and Uncertainties Research regarding differing life histories of a wild Steelhead Trout population, 2) Assessing the Status and Trend of Adult Natural Origin Fish Populations, and 3) Monitoring and Evaluating the Effectiveness of Tributary Habitat Actions Relative to Environmental, Physical, or Biological Performance Objectives. During summer and fall 2024, we PIT-tagged 2,075 Steelhead Trout parr (age-0 and age-1), in the Trout Creek and upper Wind River watersheds. Age-0 and age-1 parr were at high abundance due to strong spawner numbers in 2023 and 2024 (Estimate of 932 in 2023; and a moderate number of spawners in 2024 though an estimate is not available at this time; 24-year median = 619; Charlie Cochran, WDFW Fish Biologist, personal commun., 2025). An additional 100 parr (age-0 and age-2 or older) were tagged to provide fish for estimating detection efficiencies at PTISs. Steelhead Trout parr were recaptured and detected through repeat headwater sampling, smolt trapping, instream PTISs and Columbia River mainstem PIT-tag detection. We maintained and upgraded six instream PTISs to detect PIT-tagged Steelhead Trout parr, smolts, and adults, providing data for population assessments, and life-cycle research. Detection data from PIT-tagged adult Steelhead Trout at our four primary PTISs allow assessment of adult escapement to tributary watersheds within the Wind River subbasin. Most adult Steelhead detected were tagged by WDFW at an adult trap in the lower Wind River and some as smolts at screw traps. Adult Steelhead Trout detection efficiency estimates at our primary PTIS in Trout Creek have been greater than 99 percent during eight of the past ten years and have exceeded 97% at our primary PTIS in the Wind River during nine of the past ten years. Adult escapement estimates to tributary watersheds are helping evaluate the efficacy of the 2009 removal of Hemlock Dam from rkm 2.0 of Trout Creek. Detections at the instream PTISs have shown trends of age-0 and age-1 Steelhead Trout parr emigration from natal areas during summer and fall, in addition to the expected movement of parr and smolts in spring. Our data suggest that most fish from a cohort that migrate downstream will do so at age-1 for additional rearing in the subbasin. We have estimated that from 15% to 56% of parr tagged as age-0 fish in headwater areas make downstream migrations at age-1 for additional rearing (median = 29%). We have estimated that from 1% to 27% of Steelhead Trout parr, tagged as age-1 fish, make downstream migrations during fall (median = 3.5%): this is especially pronounced in the upper Wind River portion of the watershed (median 10%). These findings raise questions about parr rearing habitat and whether migrations are density- or habitat-quality driven, and answers to such questions are long-term goals of this study as active and passive habitat restoration occurs. Repeat sampling at sites in the subbasin within and between years has enabled assessment of juvenile Steelhead Trout growth patterns. Growth rates (relative change in weight) of age-0 PIT-tagged parr during summer have been similar across the subbasin, though slightly lower in the Trout Creek watershed (Trout Creek sites, 0.57 - 1.25%/day, Wind River sites = 0.94 – 2.48%/day). Summer growth rates have been lower for age-1 parr in the Trout Creek watershed than the upper Wind River watershed (Trout Creek sites, = 0.03 - 0.14%/day, Wind River sites = 0.05 – 0.19%/day). Yearly relative growth was similar across the subbasin for both age-0 and age-1 tagged parr. During 2024 range of mean daily relative change in weight from the previous year for age-0 tagged fish was: Trout Creek sites = 0.96 – 1.41%/day, Wind River sites = 0.80 – 1.30 %/day; and for age1 tagged fish was: Trout Creek sites = 0.37 - 0.50 %/day, Wind River sites = no age-1 tagged recaptures. Non-native Brook Trout are present in the subbasin, chiefly the Trout Creek watershed. Mean percent-of-catch that is Brook Trout, at four sample sites in Trout Creek, has declined from the period 1998 – 2003 to the period 2011 – 2024 (range of decrease = 19 - 82%). Percent-of-catch (of 64 samples, 52 < 10%; 3 > 20%) and number of Brook Trout at the Trout Creek sites from 2011 through 2024, though variable, have generally declined. Habitat restoration evaluation efforts are critical to ensure efficient use of money and resources but require long-term monitoring of population metrics. Detailed viable salmonid population and life-history data, such as that provided by PIT-tagging and instream PTIS networks inform fisheries policy and management and enable assessment of long-term effects of habitat restoration actions such as the removal of Hemlock Dam on Trout Creek and proposed major instream habitat restoration in the upper Wind River.

Washington↗