USGS Science⌕ Search

SEARCH · USGS Science

Results for “Hydrology and Earth Systems Science”

Search indexed USGS publications on groundwater, aquifers, geologic maps, mineral resources and earthquakes. Explore source records by subject and place.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

1,661 records · Page 63Linked to original sources
USGS2025

Unveiling coseismic deformation from differenced legacy aerial photography and modern lidar topography: The 1983 M6.9 Borah Peak earthquake, Idaho, USA

The 1983 M6.9 Borah Peak, Idaho, earthquake is one of the largest historical normal fault earthquakes in the western United States. We quantified meter-scale vertical change along the 35 km-long rupture using topographic differencing of 1966 aerial imagery and 2019 lidar-derived data. The initial differencing results are largely obscured by horizontal and vertical georeferencing errors and flight-line stripes. Our error corrections are designed to be insensitive to the coseismic deformation and reduced error by 50%. We calculated vertical separation and resolved a maximum of 2.02 ± 0.46 m at Doublespring Pass. Our vertical separation measurements are generally consistent with those from prior studies using field data and post-earthquake topographic data. However, the differencing measurements are a few decimeters lower than these prior measurements, indicating that differencing can isolate historical from prehistoric earthquake deformation. Our study demonstrates that revisiting historical earthquakes can provide new insights into the magnitude and patterns of coseismic deformation.

Idaho↗
USGS2025

Hydrogeology of unconsolidated and bedrock aquifers along the Salmon River, including Malone, Franklin County, New York

The U.S. Geological Survey, in cooperation with the New York State Department of Environmental Conservation, investigated the hydrogeology of the unconsolidated and bedrock aquifers along the Salmon River corridor in northern Franklin County, New York. The study area covers roughly 147 square miles and includes the village of Malone and parts of the Towns of Malone, Bellmont, Burke, Constable, Westville, Bangor, Duane, and Franklin, New York. Groundwater is the primary source for water supply within the study area. Eighty-three percent of all residents use public water supplied from two production wells that draw water from a thick, highly productive sand-and-gravel aquifer likely receiving induced infiltration from the Salmon River. Twenty-four percent of the 187 verified domestic wells in the study area outside of the production well service area boundary are screened in typically discontinuous deposits of stratified sand and sand and gravel. Seventy-six percent of the wells are completed in bedrock aquifers including the Potsdam Sandstone, metamorphic rocks, Theresa Formation, and unknown bedrock. Characterizing and understanding potential groundwater resources is critical for protecting the quality of the groundwater. The information in this report may be used to guide delineation of groundwater contributing areas, assess potential threats to aquifers from both point and nonpoint sources, respond to contamination from spills or leaks from underground storage facilities or other sources, and to support assessments for future development of municipal water supplies.

New York↗
USGS1976

Technique for estimating magnitude and frequency of floods in Kentucky

This report presents flood magnitude and frequency relations applicable to unregulated streams in Kentucky. The relations are based on flood data at 117 gaging stations in Kentucky and 14 in adjacent states having 10 or more years of record not significantly affected by man-made changes. Equations that relate flood magnitude and frequency to contributing drainage area in 16 geographic areas may be used to estimate magnitude of future floods with recurrence intervals of as much as 100 years on gaged and ungaged streams having drainage areas of 10 to 4,300 square miles (25 .9 to 11,100 square kilometers) . Estimating equations are also presented in graphical form for the convenience of the user. Additional graphs are presented to estimate flood magnitude for selected recurrence intervals along the Cumberland, Kentucky, and Ohio Rivers.

Kentucky↗
USGS2025

The δ13C signature of dissolved organic and inorganic carbon reveals complex carbon transformations within a salt marsh

Coastal wetlands have high rates of atmospheric CO 2 uptake, which is subsequently respired back to the atmosphere, stored as organic matter within flooded, anoxic soils, or exported to the coastal ocean. Transformation of fixed carbon occurs through a variety of subsurface aerobic and anaerobic microbial processes, and results in a large inventory of dissolved carbon. Carbon source and the roles of aerobic respiration, sulfate reduction, and methane cycling were evaluated within salt marsh peat and the underlying sandy subterranean estuary. There is a large increase in dissolved inorganic carbon (DIC, 7,350 ± 3,900 μmol L −1 ), dissolved organic carbon (DOC, 1,040 ± 1,480 μmol L −1 ) and CH 4 (14.5 ± 33.3 μmol L −1 ) within the marsh porewaters compared to creek waters. Alkalinity production (5,730 ± 2,170 μeq L −1 ) and sulfate removal (1,810 ± 1,970 μmol L −1 ) indicate anaerobic respiration, however, relative contributions from the various decomposition pathways cannot be identified due to overlapping geochemical signatures. The δ 13 C of the DOC (−29.0 ± 3.7‰) and DIC (−11.2 ± 1.1‰) produced within the marsh differed from the bulk soil organic matter δ 13 C (−14.5 ± 0.2‰). We explore a variety of mechanisms that could result in co-occurring depleted δ 13 C-DOC and enriched δ 13 C-DIC compared to the bulk soil organic carbon pool and salt marsh vegetation, including selective mineralization, production of δ 13 C-depleted bacterial biomass, and methane-derived DOC. While important questions remain about carbon cycling pathways, we found evidence of a cryptic methane cycle. Alteration of the δ 13 C of carbon species complicates source attribution in solid and dissolved phases and careful consideration should be used when carbon is partitioned between in situ salt marsh production and external marine and terrestrial sources.

Massachusetts↗
USGS2000

Suspended-sediment characteristics in the Housatonic River Basin, western Massachusetts and parts of eastern New York and northwestern Connecticut, 1994-96

Suspended-sediment concentrations, discharges, loads, and yields were determined for eight subbasins in the Housatonic River Basin in western Massachusetts, eastern New York, and northwestern Connecticut from April 1994 through March 1996. Suspended-sediment samples were collected at three continuous-record sediment stations and at four partial-record sediment stations. Suspended-sediment concentrations in samples collected during the period of study ranged from less than 0.5 to 3,400 milligrams per liter, and concurrent streamflows ranged from 0.03 to 126 cubic feet per second per square mile at the seven stations. Median suspended-sediment concentrations in samples collected at each station ranged from 7 to 61 milligrams per liter. Median streamflows during suspended-sediment sampling ranged from 1.86 to 5.88 cubic feet per second per square mile. Instantaneous suspended-sediment yields ranged from less than 0.005 to 185 tons per day per square mile, and medians ranged from 0.03 to 1.12 tons per day per square mile at the seven stations. Total suspended-sediment loads (mass) from April 1994 through March 1996 at the continuous-record sediment stations were 11,603 tons at Housatonic River near Great Barrington, 7,929 tons at Green River, and 54,347 tons at Housatonic River near Ashley Falls. Suspended-sediment load during January 1996 at the Green River station accounted for about 54 percent of the total suspended-sediment load for the Green River during the 2-year study. Suspended-sediment load on January 19 and 20, 1996, at the Green River station accounted for about 50 percent of the January 1996 suspended-sediment load, or about 27 percent of the total suspended-sediment load during the 2-year study. This large suspended-sediment transport was the result of rainfall and snowmelt on January 19 and 20, 1996--the equivalent of a 5- to 6-inch rain storm in the Green River subbasin. Total suspended-sediment loads during the 2-year study at the partial-record sediment stations were 3,052 tons at Williams River, 1,758 tons at Ironworks Brook, and 17,927 tons at Konkapot River. Suspended-sediment yields from April 1994 through March 1996 at the continuous-record sediment stations were 21 (tons/yr)/mi 2 at Housatonic River near Great Barrington, 78 (tons/yr)/mi 2 at Green River, and 58 (tons/yr)/mi 2 at Housatonic River near Ashley Falls. Suspended-sediment yields during the 2- year study at the partial-record sediment stations were 35 (tons/yr)/mi 2 at Williams River, 78 (tons/yr)/mi 2 at Ironworks Brook, and 147 (tons/yr)/mi 2 at Konkapot River. Suspended-sediment yields were estimated for two subbasins in the Housatonic River Basin--Schenob Brook at Sheffield, a partial-record sediment station, and the area adjacent to the Housatonic River between Great Barrington and Ashley Falls. The estimate of suspended-sediment yield at Schenob Brook at Sheffield of 82 (tons/yr)/mi 2 is comparable to the yield determined for the Green River and Ironworks Brook. The estimate of suspended-sediment yield for the area adjacent to the Housatonic River between Great Barrington and Ashley Falls was 395 (tons/yr)/mi 2 . This estimated suspended-sediment yield was 2.7 to 18.8 times greater than that estimated for any of the other subbasins. Several basin and land-use characteristics thought to affect suspended-sediment transport in the subbasins were compared to the suspended-sediment yields. The characteristics that seemed to affect suspended-sediment discharge were dams, which contributes to decreased yields; Hadley, Limerick, Linlithgo, Saco, and Winooski silt loam soils (high erodibility soils), which contributes to increased yields; stratified-drift deposits, which contributes to increased yields; and agricultural and open land, which contributes to increased yields. The effect of stratified-drift deposits on suspended-sediment discharge is thought to be greater when those deposits are of glaciolacustrine (generally clay, silt, and fine sand), rather than glaciofluvial (clay, silt, sand, gravel, and cobbles) origin. The silt loam soils, glaciolacustrine deposits, and agricultural and open land are interrelated, inasmuch as the silt loam soils generally are associated with glaciolacustrine deposits and agricultural activities.

Connecticut, Massachusetts, New York↗
USGS2026

Tectonic implications of transitional melting regimes from petrological, geochronological, and compositional characterization of the ophiolitic Seventymile terrane, Alaska, USA

New geochemical, U-Pb geochronology, and Sr-Nd-Hf isotope data provide evidence for the tectonic evolution of the Seventymile terrane in interior Alaska, USA. Ultramafic and mafic rocks of the Seventymile terrane are thought to represent components of a dismembered ophiolite and provide unique constraints on regional terrane evolution and accretion. The Seventymile ophiolite represents fragments of the Devonian to Permian Slide Mountain Ocean (SMO) that separated allochthonous and parautochthonous continental fragments of western North America. It now occurs as multiple thrust sheets containing Permian mafic and ultramafic rocks overlying and/or possibly imbricated with amphibolite-facies supracrustal rocks of allochthonous Yukon-Tanana terrane and parautochthonous North America. Seventymile klippen contain variably serpentinized peridotite, primarily harzburgite, low-grade meta-mafic rocks, and minor oceanic sedimentary rocks (argillite, chert, limestone, and metasandstone). Mafic rocks include gabbro to diabase, typically as dikes, veinlets, or rare massive stocks intruding peridotite. Mafic rocks also include greenstones of the Seventymile assemblage in klippen structurally underlying, and in shear zone contact with, ultramafic klippen. New trace element and radiogenic isotope data from mafic magmatic rocks associated with the Seventymile ultramafic bodies show evidence for a weakly subduction-modified mantle source, like the mantle source of normal mid-ocean-ridge basalt (N-MORB) or back-arc basin basalt (BABB). Seventymile assemblage greenstones are more heterogeneous. They range from N-MORB to enriched mid-ocean-ridge basalt (E-MORB) and ocean-island basalt (OIB), with a subset of samples indicative of continental arc affinity. These geochemistry results indicate that distinct tectonic environments are represented by at least two, and possibly three, lithological and structural units comprising the Seventymile terrane. Hf-Nd isotope systematics are consistent with a depleted MORB mantle (DMM)−like component that overlaps with Pacific MORB. Primary zircon is rare, but new in situ U-Pb data for gabbro and greenstone indicate ca. 274−272 Ma peak zircon and titanite crystallization. Scattered younger zircons define a ca. 255 Ma zircon peak and correspond to secondary crystallization associated with baddeleyite reaction of high-Si fluids during low-grade metamorphism. If Seventymile suites are contemporaneous, obduction associated with the closure of the SMO resulted in the stacking of ophiolitic packages representing distinct tectonomagmatic settings across the transition from pericontinental, to epicontinental, to distal ocean back-arc. Intrusions hosted in klippe of ultramafic rocks, plus the least subduction-modified greenstones underlying them, geologically and compositionally resemble Slide Mountain rocks of the Campbell Range formation in eastern Yukon and may provide a new piercing point across the Tintina fault.

Alaska, British Columbia↗
USGS2008

Late Jurassic igneous rocks in south-central Arizona and north-central Sonora: Magmatic accompaniment of crustal extension

In Middle Jurassic time, the region of south-central Arizona and north-central Sonora was part of a continental magmatic arc. In early Late Jurassic time, this magmatic arc gradually gave way to a extensional regime. Late Jurassic rifting was accompanied by magmatism, typically small in volume relative to the preceding arc phase, and with mildly alkaline compositions. These igneous rocks, the 158- to 146-Ma Ko Vaya Suite, make up much or most of the Comobabi, Artesa, Quijotoa, and Brownell Mountains and Sierra del Cobre. The bimodal Ko Vaya Suite constitutes volcanic to shallow plutonic complexes comprising chiefly trachyandesite and volcanic wacke; monzodiorite; rhyolite porphyry; distinctive, compositionally and texturally heterogeneous, quartzpoor granite (the Ko Vaya Granite); A-type perthite granite and quartz syenite; aphyric intrusive rhyolite; and locally abundant hematite veins. Alteration, particularly potassic alteration, is widespread. Minimally altered samples of granite and quartz syenite have the composition of trachyte or trachydacite and alkali rhyolite. The Ko Vaya Suite is nonconformably overlain by the latest Jurassic Sand Wells Formation, approximately correlative with the Glance Conglomerate, the basal unit of the Bisbee Group. The Ko Vaya Suite is in some ways analogous to the coeval Independence dike swarm of eastern California, but less diffuse. We infer that the Ko Vaya Suite represents especially intense or focused Late Jurassic extensional magmatism.

Arizona, Sonora↗
USGS2025

Estimated ultimate recovery (EUR) Prediction for Eagle Ford Shale using integrated datasets and artificial neural networks

The estimated ultimate recovery (EUR) is an important parameter for forecasting oil and gas production and informing decisions regarding field development strategies. In this study, we combined site-specific geologic, completion, and operational parameters with the predictive capabilities of machine learning (ML) models to predict EURs of the wells for the Eagle Ford Marl Continuous Oil Assessment Unit. We developed an extensive dataset of wells that have produced from the lower and upper Eagle Ford Shale intervals and reduced the model complexity using principal component analysis. We tested the ML models and estimated the sensitivities of ML-predicted EURs to changes in the values of different input variables. The results of applying the optimized ML model to the Eagle Ford suggest that the approach developed in this study could be promising. The ML estimates of the EURs fit the DCA-based values with an R 2 ~ 0.9 and a mean absolute error of ~36 × 10 3 bbl. In the lower Eagle Ford Shale, the EUR estimates were found to be most sensitive to changes in porosity, net thickness of the interval, clay volume, and the API gravity of the oil; and that in the upper Eagle Ford Shale they were most sensitive to changes in the total organic carbon and water saturation, which suggests that it could be important to consider these parameters in assessing these intervals or close analogs.

Louisiana, Mississippi, Texas↗
USGS2024

The geometry of fault reactivation and uplift along the central part of the Maacama fault zone, northern California Coast Ranges (USA)

Fault reactivation of bedrock structures in active fault zones influences stress state and earthquake rupture phenomena through the introduction of weak slip surfaces that impact fault zone geometry and width. Yet, geometric relationships between modern faults and older reactivated faults are difficult to quantify in rocks that have experienced multiple deformation episodes. We used new geologic mapping, geomorphic tools, and structural modeling to quantify rock uplift and subsurface fault geometry of the central part of the Maacama Fault Zone near Ukiah, California, USA, and the surrounding area. Results suggest that the northern Mayacamas Mountains are in a tectonically driven disequilibrium, with differential rock uplift focused on the western side of the range. Steeply east-dipping fault surfaces and splays characterize the geometry of the Maacama Fault Zone. We mapped two newly identified faults to the east of the main Maacama Fault, the Cow Mountain–Mill Creek Fault, and Willow Creek Fault, which align with a moderately east-dipping cluster of microseismicity between 4–10 km depth beneath the Mayacamas Mountains. Static stress modeling on the Maacama Fault Zone and newly identified faults to the east quantify slip tendency values of 0.5–0.4, which suggests that the faults are moderately to poorly suited for slip in the modern stress field and may be weak. We infer that modern uplift is driven by oblique reverse, up-to-the-east, dip-slip motion on the reactivated Cenozoic Cow Mountain–Mill Creek and Willow Creek Faults as material is advected through a restraining bend on the Maacama Fault. This study shows that reactivated bedrock faults increase the fault zone width and introduce fault surfaces that contribute a component of vertical deformation and uplift in major strike-slip fault zones. Deformation is accommodated on an interconnected network of new and reactivated faults that delineate a complex seismic hazard.

California↗
USGS1999

Estimation of magnitude and frequency of floods for streams in Puerto Rico: New empirical models

Flood-peak discharges and frequencies are presented for 57 gaged sites in Puerto Rico for recurrence intervals ranging from 2 to 500 years. The log-Pearson Type III distribution, the methodology recommended by the United States Interagency Committee on Water Data, was used to determine the magnitude and frequency of floods at the gaged sites having 10 to 43 years of record. A technique is presented for estimating flood-peak discharges at recurrence intervals ranging from 2 to 500 years for unregulated streams in Puerto Rico with contributing drainage areas ranging from 0.83 to 208 square miles. Loglinear multiple regression analyses, using climatic and basin characteristics and peak-discharge data from the 57 gaged sites, were used to construct regression equations to transfer the magnitude and frequency information from gaged to ungaged sites. The equations have contributing drainage area, depth-to-rock, and mean annual rainfall as the basin and climatic characteristics in estimating flood peak discharges. Examples are given to show a step-by-step procedure in calculating a 100-year flood at a gaged site, an ungaged site, a site near a gaged location, and a site between two gaged sites.

Puerto Rico↗
USGS2026

High-resolution magnetic survey using an unoccupied aerial vehicle to constrain buried lava flow geometry, volume, and eruptive history of Little Cones, Crater Flat, Nevada

Magnetic surveys are an important tool used to augment geologic mapping in distributed volcanic fields. Using magnetic anomalies, it is possible to model the geometry of shallowly buried volcanic features, such as conduits, sills, and lava flows. This subsurface mapping is important for understanding eruption dynamics and emplacement of lava flows, and it sometimes reveals buried volcanoes no longer visible at the surface. These data are critical to better interpret the numbers, styles, and magnitudes of eruptions in distributed volcanic fields and their associated volcanic hazards. New advances in unoccupied aerial vehicles (UAVs) offer an attractive middle range of resolution and aerial coverage between ground-based magnetic surveys and aeromagnetic surveys. Here, we present the results of a UAV fluxgate magnetic survey of the Little Cones, Nevada, scoria cones, which have been the target of previous ground and aeromagnetic surveys. The magnetic anomalies at Little Cones are of interest because the surrounding alluvium conceals lava flows that erupted from Little Cones, making it very difficult to understand the volume and morphology of lava flows from geologic mapping alone. Nonlinear inversion of UAV-collected magnetic data were used to model the thickness and morphology of buried Little Cones’ lava flows with higher precision than achieved previously. The sequence of events and calculated flow characteristics are then interpreted. The total volume of Little Cones, including concealed lava flows, is approximately 0.016 cubic kilometer, and the initial sheet flow erupted in less than 24 hours. The findings presented herein demonstrate that UAV-based magnetic surveys are a reliable method of data collection and an efficient alternative to other survey methods, facilitating development of a three-dimensional perspective of distributed volcanic fields.

Nevada↗
USGS2026

Suitability of ground-motion models for seismic hazard assessment in Puerto Rico and the U.S. Virgin Islands

We perform linear mixed-effects analyses with a ground-motion dataset to evaluate how well ground-motion models (GMMs) fit active crustal, subduction interface, and subduction intraslab earthquakes in Puerto Rico and the US Virgin Islands. Most of the GMMs reproduce the spatial variation in peak ground motions with earthquake magnitude and rupture distance but predict ground motions 0.3–1.0 natural log units (35%–270%) greater than observations. Two GMMs developed for Puerto Rico that are based on ground-motion records from mostly magnitude 4–5 earthquakes do not perform as well. We attribute the period-dependent overprediction in observed ground motions to differences between the observed site response and the linear site response in the GMMs. Consequently, we developed region-specific GMMs by adjusting the period-dependent linear site response coefficients and period-dependent constant coefficients to remove most of the bias between predicted ground motions and observations. For the analysis, we compile ground-motion records and process waveforms to build a dataset with 10,127 records at 72 stations from 849 magnitude 4.0–6.4 earthquakes between 1 January 2006 and 31 March 2024. The earthquakes include active crustal, subduction interface, and subduction intraslab events. We evaluate the GMMs using the time-averaged shear wave speed in the top 30 m ( ), which we compile from site surveys and proxy values computed from horizontal to vertical spectral ratios. Site terms exhibit strong consistency across GMMs and crustal and subduction earthquakes, indicating that the linear mixed-effects analysis successfully isolates the effects of local site response. The event terms show little spatial correlation and more substantial variability than in other regions, which we attribute to uncertainties in the earthquake magnitudes. This analysis guides the selection of the GMMs for the 2025 update of the National Seismic Hazard Model for Puerto Rico and the US Virgin Islands.

Puerto Rico, U.S. Virgin Islands↗
USGS2024

Assessment and characterization of ephemeral stream channel stability and mechanisms affecting erosion in Grand Valley, western Colorado, 2018–21

The Grand Valley in western Colorado is in the semiarid Southwest United States. The north side of the Grand Valley has many ungaged ephemeral streams, which are of particular interest because (1) the underlying bedrock geology, Late Cretaceous Mancos Shale, is a sedimentary rock deposit identified as a major salinity contributor to the Colorado River and (2) despite infrequent streamflows of short duration, monsoon-derived floods in these ephemeral streams can carry substantial amounts of sediment downstream, affecting upstream and downstream banks and channel cross sections. The study area is of interest, because salinity, or the total dissolved solids concentration, in the Colorado River causes an estimated $300 million to $400 million per year in economic damages in the United States, and it is estimated 62 percent of the Upper Colorado River Basin’s total dissolved solid loads originate from geologic sources. In an effort to minimize salt contributions to the Colorado River from public lands administered by the Bureau of Land Management, a comprehensive salinity control approach is typically used to reduce nonpoint sources of salinity through land management techniques and practices. In 2018, the U.S. Geological Survey, in cooperation with the Bureau of Land Management, began an assessment of ephemeral streams located on the north side of the Grand Valley, western Colorado, to characterize stream channel stability and identify mechanisms affecting erosion. The U.S. Geological Survey developed a method for automatically extracting channel cross-section geometry from existing remotely sensed terrain models. Based on estimated flood stage and surrogate streamflows, hydraulic characteristics were calculated. Furthermore, the channel geometries and hydraulic characteristics were used to estimate channel stability using a statistical model. Cross-section stabilities were determined from a stream channel stability assessment for a subset of 1,406 visited (field observed) locations out of 13,415 cross sections, which were delineated from remotely sensed terrain models. The application of Manning’s resistance equation in combination with multiple logistic regression models demonstrated channel stability can be estimated with a 0.845 goodness of fit for a validation dataset when using a combination of drainage area, width-to-depth ratio, sinuosity, and shear stress as the explanatory variables. Stream channel stability was extrapolated for 13,415 unvisited (not field observed) cross sections using the multiple logistic regression model and defined explanatory variables. Mapping of the ephemeral streams and their associated stabilities may be used by the Bureau of Land Management to prioritize areas for remediation or changes in management strategies to reduce sediment and salinity loading to the Colorado River. The study found channel stability within the ephemeral streams to be spatially variable, longitudinally discontinuous, and dictated by changes in channel bed slope. The stable ephemeral streams were relatively wide and shallow and often had smaller drainage areas with less potential for producing shear stresses capable of overcoming channel adhesion. A change in channel bed slope can provide the means necessary to generate shear stresses appropriate to initiate erosion and a subsequent stability transition to incising channels. Channel widening happens when either or both banks of an incising channel reach a critical height for mass wasting, or when channel curvature causes higher sidewall stress. Regardless, widening channels can promote increases in sinuosity and subsequently reduce steep channel bed slopes. Consequently, stable and widening channels can have comparable bed slopes, making channel bed slope a poor explanatory variable to predict channel stability overall, despite its function to initiate channel instability. The results were based on a surrogate 0.10 annual exceedance probability (AEP; return period equal to the 10-year flood) interval streamflow, although it was recognized fluctuations in streamflow would also affect channel stability. Past and current changes within the study area affect streamflow; therefore, mechanisms affecting erosion include land use disturbances, soil compaction, loss of vegetation cover, drought, less frequent and more extreme precipitation, and fires—which all intensify the potential runoff and erosion within the study area.

Colorado↗
USGS2024

Telecommunications fiber for sensing earthquake aftershocks: Progress and hurdles

Aftershocks offer valuable clues to earthquake behavior. The challenge: quickly deploying sensors to capture the early details of earthquake ruptures within the zone of aftershocks. Telecommunication fibers might be an answer, providing denser networks in otherwise difficult areas, potentially faster than traditional methods.

Conference Paper↗
USGS2026

Modeled groundwater and surface-water interactions surrounding Mobile Bay, Alabama, 2008–15

The U.S. Geological Survey, in cooperation with the Gulf Coast Ecosystem Restoration Council, has used MODFLOW 6 to develop a groundwater-flow model to simulate groundwater and surface-water interactions in the Mobile Bay, Alabama, area. The model results indicated that, on average, groundwater discharge near the coastline is equal to 2.5 percent of the surface water that flows into the bay. The model was also used to determine how changes in recharge, sea level, and groundwater pumping affect groundwater levels and discharge rates. The results indicate that more groundwater discharge occurred in the winter and spring when recharge was higher, sea level was lower, and groundwater pumping was lower than during the summer and fall. Additionally, the amount of emergent groundwater was closely related to sea level; when sea level was higher, there was more area with emergent groundwater. Furthermore, the depth of nonemergent groundwater was related to trends in recharge and pumping. During periods of increased recharge and reduced pumping, a significant portion of the model area exhibited a depth to the water table of less than 1 meter below the land surface. Conversely, during periods of decreased recharge and heightened pumping, much of the model area showed a depth to the water table ranging from 1 to 5 meters.

Alabama, Mississippi↗
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

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

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

Journal of Agricultural, Biological and Environmen↗
USGS2025

Manipulation of the Symbiodiniaceae microbiome confers multigenerational impacts on symbioses and reproductive ecology of its Exaiptasia diaphana host

Symbiodiniaceae-associated microbiota strongly influence cnidarian symbioses. We systematically reduced the bacterial and fungal communities associated with Symbiodiniaceae to study potential effects on the cnidarian holobiont Exaiptasia diaphana (Aiptasia). Clonal anemones were inoculated with xenic Breviolum minutum (SSB01) and microbiome-manipulated cultures after antibacterial or antifungal treatment. The asexual reproduction of pedal laceration allowed for three generations of clonal aposymbiotic Aiptasia to be utilised in this study, from the initial adult generation (G0), to the first (G1), and second (G2) generation. We inoculated small and large G1 Aiptasia with SSB01 algae and monitored onset of symbiosis, rate of algal proliferation, and holobiont characteristics. Sequencing the 16S and 18S rRNA gene regions identified significant differences in the bacterial and fungal communities of the G0 and G1 generations, alongside differences between the size classes of small and large G1 anemones. The microbiome of larger G1 individuals was distinct to the smaller G1 anemones, suggesting a microbiome maturation process. Control Breviolum minutum cultures exhibited a significantly greater proliferation rate in large G1 anemones when compared to antibacterial or antifungal treated cultures, whereas the opposite trend was documented in the small G1 anemones. Although no differences were observed between algal photochemical parameters, or the growth and polyp activity of G1 juveniles, we observed a significant influence in the production of G2 clones between treatments. Overall, we provide strong ecological implications of manipulating Symbiodiniaceae microbiome, not for the algae themselves, but for the maturation of the host Aiptasia, as well as for the cnidarian holobiont over multiple generations.

The ISME Journal↗