SEARCH · USGS Science

Results for “Limnology and Oceanography”

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

SourceResearch connectionFrom yearThrough year

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

1,656 records · Page 43Linked to original sources

USGS2021

Geology, mineralogy, and cassiterite geochronology of the Ayawilca Zn-Pb-Ag-In-Sn-Cu deposit, Pasco, Peru

The Ayawilca deposit in Pasco, Peru, represents the most significant recent base-metal discovery in the central Andes and one of the largest undeveloped In resources globally. As of 2018, it hosts an 11.7 Mt indicated resource grading 6.9% Zn, 0.16% Pb, 15 g/t Ag, and 84 g/t In, an additional 45.0 Mt inferred resource grading 5.6% Zn, 0.23% Pb, 17 g/t Ag, and 67 g/t In, and a separate Sn-Cu-Ag inferred resource of 14.5 Mt grading 0.63% Sn, 0.21% Cu, and 18 g/t Ag. Newly obtained U–Pb dates for cassiterite by LA-ICP-MS (22.77 ± 0.41 and 23.05 ± 2.06 Ma) assign the Ayawilca deposit to the Miocene polymetallic belt of central Peru. The polymetallic mineralization occurs as up to 70-m-thick mantos hosted by carbonate rocks of the Late Triassic to Early Jurassic Pucará Group, and subordinately, as steeply dipping veins hosted by rocks of the Pucará Group and overlying Cretaceous sandstones-siltstones of the Goyllarisquizga Group. Relicts of a distal retrograde magnesian skarn and cassiterite (stage pre-A) were identified in the deepest mantos. The volumetrically most important mineralization at Ayawilca comprises a low-sulfidation assemblage (stage A) with quartz, pyrrhotite, arsenopyrite, chalcopyrite, Fe-rich sphalerite, and traces of stannite and herzenbergite. Stage A sphalerite records progressive Fe depletion, from 33 to 10 mol% FeS, which is compatible with the observed transition from low- to a subsequent intermediate-sulfidation stage (B) marked by the crystallization of abundant pyrite and marcasite. Finally, during a later intermediate-sulfidation stage (C) sphalerite (up to 11 mol% FeS), galena, native bismuth, Cu-Pb-Ag sulfosalts, siderite, Mn-Fe carbonates, kaolinite, dickite, and sericite were deposited. This paragenetic evolution shows striking similarities with that at the Cerro de Pasco Cordilleran-type polymetallic deposit, even if at Ayawilca stage C did not reach high-sulfidation conditions. The occurrence of an early retrograde skarn assemblage suggests that the manto bodies at Ayawilca formed at the transition between distal skarn and skarn-free (Cordilleran-type) carbonate-replacement mineralization. Mineral assemblages define a T- f S 2 evolutionary path close to the pyrrhotite-pyrite boundary. Buffering of hydrothermal fluids by underlying Devonian carbonaceous phyllites of the Excelsior Group imposed highly reduced conditions during stage A mineralization (log f O 2 < − 30 atm). The low f O 2 favored efficient Sn mobility during stages pre-A and A, in contrast to other known ore deposits in the polymetallic belt of central Peru, in which the occurrence of Sn minerals is minor. Subsequent cooling, progressive sealing of vein walls, and decreasing buffering potential of the host rocks promoted the shift from low- (stage A) to intermediate-sulfidation (stages B and C) states. LA-ICP-MS analyses reveal significant In contents in Fe-rich sphalerite (up to 1.7 wt%), stannite (up to 1908 ppm), and chalcopyrite (up to 1185 ppm). The highest In content was found in stage A sphalerite that precipitated along with chalcopyrite and stannite, thus pointing to the early, low-sulfidation assemblage as prospective for this high-tech metal in similar mineral systems. Indium was likely incorporated into the sphalerite crystal lattice via Cu + + In 3+ ↔ 2 Zn 2+ and (Sn, Ge) 4+ + (Ga, In) 3+ + (Cu + Ag) + ↔ 4 Zn 2+ coupled substitutions. Indium incorporation mechanisms into the stannite and chalcopyrite crystal lattices remain unclear.

Pasco↗

USGS2025

Origin and evolution of mafic volcanism associated with 3 m.y. of andesite production at the Goat Rocks volcanic cluster, southern Washington Cascade Range

More than 3 m.y. of mafic volcanism near the Goat Rocks volcanic cluster in the southern Washington Cascade Range, USA, lends insight into the evolution of basalts and the subarc mantle at a long-lived, major arc volcanic locus. We contribute field observations, 40 Ar/ 39 Ar dates, paleomagnetic directions, and bulk rock and mineral compositions to characterize nine mafic units that erupted in association with the Goat Rocks volcanic cluster. The time frame of mafic volcanism, ca. 3.6 Ma to 60 ka, encompasses the lifespan of the central volcanic cluster (3.1 Ma to 115 ka), with a lull from ca. 2.7 Ma to 1.4 Ma. A climactic period of voluminous mafic activity and far-traveled lava flows, including construction of the Hogback Mountain shield volcano, coincided with voluminous andesite eruptions from the central volcanic cluster. The basaltic rocks in the Goat Rocks area are calc-alkaline to barely tholeiitic and have high field strength element depletion relative to large-ion lithophile elements characteristic of calc-alkaline basalts (CAB) of the Cascade volcanic arc. Unlike at neighboring andesitic volcanic centers (Mounts Adams, St. Helens, and Rainier), no other mafic end members such as high-aluminum olivine tholeiite (HAOT) or intraplate-type basalt (IPB) are present at or near the Goat Rocks volcanic cluster, although some of the calc-alkaline basalts in this study have IPB-like affinities. The Goat Rocks mafic units exhibit two main temporal trends in composition: (1) the most primitive basalts erupted earlier, compared to less primitive and more evolved compositions later, and (2) high field strength element concentrations are higher in the younger basalt units relative to the oldest two. In contrast to these temporal trends, the mafic units define two compositional groups that recur through time, a low-Sr and a high-Sr group, each with distinct trace element and Sr and Nd isotope ratios. Although radiogenic isotope ratios are generally aligned with High Cascades CAB and HAOT, some extend toward IPB of Mount Adams and Simcoe Mountains volcanic field. Olivine-dominated crystal fractionation at shallow pressure from a small range of parent magma compositions accounts for much of the variation among the basalts and basaltic andesites. A high-pressure fractionation model is plausible for only one of the youngest basalt units (basalt of Walupt Lake volcano). Mafic recharge and crustal assimilation accounts for the incompatible-element enriched composition of basaltic andesites erupted during construction of the largest andesitic centers, further supporting sustained basalt mass flux and thermal energy driving andesite genesis. We model the most primitive members of the Goat Rocks mafic units as partial melts of successively less depleted mantle in time. Variable degrees of fluxing with fluids and melts from subduction explain the distinction between high-Sr and low-Sr groups. We propose that mantle metasomatism by ancestral subduction and fluid-flux melting is heterogeneously distributed through the local subarc mantle and played a greater role in the genesis of the high-Sr basalt group. The limited range of primitive basalt types around the Goat Rocks volcanic cluster contrasts with the much greater diversity of basalts throughout the southern Washington to northern Oregon Cascade arc. On the other hand, the central volcanic cluster encompasses nearly the entire diversity observed at neighboring composite volcanoes. In the case of the Goat Rocks area at least, and perhaps attributable to the entire region, this means that the genesis of diverse intermediate magmas is independent from and does not require vastly different parental basalt compositions.

California, Oregon, Washington↗

USGS2026

The Great Lakes Geologic Mapping Coalition—Working collaboratively to understand the geology of the Great Lakes Region

Introduction The Great Lakes Geologic Mapping Coalition (GLGMC), commonly referred to as the “Coalition,” is a partnership between the U.S. Geological Survey (USGS), the U.S. States of Illinois, Indiana, Michigan, Minnesota, New York, Ohio, Pennsylvania, and Wisconsin and the Canadian province of Ontario. The member States receive funding for geologic mapping work from the USGS National Cooperative Geologic Mapping Program (NCGMP), whereas Ontario participates as a nonfunded partner. The mission of the GLGMC is to produce three-dimensional (3D) geologic maps that depict unconsolidated sediments and near-surface bedrock in the Great Lakes region of North America. Geologic maps are the basis of most earth science investigations and help support resource exploration (energy, minerals, groundwater), natural hazard mitigation, infrastructure development, and land-use planning, all of which can be used to advance economic development and strengthen national security in the Great Lakes region. During the last few million years, the Great Lakes region has experienced repeated glacial advances and retreats, leaving behind extensive sediments, abundant natural resources, and widespread effects on the underlying bedrock geology (Swezey and others, 2022). Linked by shared histories of past glaciations, industrial agriculture, and legacy automotive, coal, steel, and manufacturing industries, the GLGMC member States collaborate to improve the understanding of the 3D distribution of the sediments overlying the region’s bedrock (fig. 1). Developing a comprehensive subsurface 3D framework of this glaciated terrain can provide earth science data to policymakers at all levels. These insights facilitate informed decisions on the exploration, use, and protection of vital resources, such as critical minerals, industrial materials, and aquifers, thereby supporting economic prosperity and the well-being of the citizens of this region. Since its inception in 1998, the Coalition has completed more than 100 geologic mapping projects across the Great Lakes region. Each project aims to deliver geologic maps, 3D datasets, and other information that improves understanding of the geology of the Great Lakes region, with an emphasis on economic and water resources. Key deliverables include 3D geologic maps and models typically portraying sediment thickness, often derived from top-of-bedrock and borehole data. These products are developed through a combination of fieldwork, subsurface modeling, and the collection and analysis of rock and sediment cores. To support Coalition goals, member States collaborate with scientists working on related STATEMAP, EDMAP, and FEDMAP projects. Coalition scientists also engage with Tribal Nations in the Great Lakes region to ensure that Tribal interests pertaining to Coalition work are addressed. Through this collaboration, the Coalition unites the efforts of State, Federal, and Tribal Nation stakeholders to advance geologic data production and enhance understanding of the geologic resources of the Great Lakes region.

Illinois, Indiana, Michigan, Minnesota, New York, ↗

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↗

USGS1996

Pesticides in surface water from three agricultural basins in south-central Georgia, 1993-97

Twenty-two of 43 pesticides analyzed were detected in 128 water samples collected from the Tucsawhatchee Creek, the Little River, and the Withlacoochee River. These streams drain agricultural basins in south-central Georgia and were sampled from March 1993 through June 1995. Herbicides were detected more frequently than insecticides. The most frequently detected herbicides were atrazine and metolachlor and the most frequently detected insecticide was carbaryl. Pesticide concentrations in the three streams were low and did not exceed U.S. Environmental Protection Agency drinking water standards. The maximum pesticide concentration was 2.6 mg/L (micrograms per liter) for propargite, a miticide detected in only one sample. The maximum concentrations of the remaining 21 pesticides were less than 0.25 mg/L. The median concentrations were equal to the method detection limit for all pesticides except atrazine (0.008 mg/L) and metolachlor (0.012 mg/L). The ratio of herbicide detections to nondetections was largest in the planting season, smaller in the harvest season and smallest in the fallow season for the three basins. The same pattern existed for the insecticide ratios in the Little River and the Withlacoochee River. Pairwise correlations between concentrations of atrazine and metolachlor, and four parameters (discharge, and concentrations of dissolved organic carbon, suspended organic carbon, and suspended sediment) were evaluated for each stream. The strongest correlations existed between metolachlor and mean daily discharge, and metolachlor and sediment in the Withlacoochee River. The only significant correlation for the Little River was between atrazine and suspended sediment.

Georgia↗

USGS2024

Complex landslide patterns explained by local intra-unit variability of stratigraphy and structure: Case study in the Tyee Formation, Oregon, USA

Lithology and geologic structure are important controls on landslide susceptibility and are incorporated into many regional landslide hazard models. Typically, metrics for mapped geologic units are used as model input variables and a single set of values for material strength are assumed, regardless of spatial heterogeneities that may exist within a map unit. Here we describe how differences in bedding thickness, grain size, inferred uniaxial compressive strength, and bedding dip control the inherent susceptibility of slopes to deep-seated failure within a single mapped geologic unit - the Tyee Formation of Oregon, USA. The Tyee, which covers over 15,000 km2 and underlies much of the Oregon Coast Range, comprises gently folded alternating beds of sandstone and siltstone deposited as turbidites, forming a 2-km thick Eocene submarine fan which has been uplifted and exhumed through the Cenozoic. Deep-seated landslides are widespread in the Tyee, but form a complex spatial pattern such that landslide density ranges from 0 to 24% of the total landscape area. These slides are often extensive and sufficiently deep to reduce local hillslope gradients, resulting in a strong negative correlation between landslide density and mean local slope. Mean annual precipitation and predicted strong ground motions from Cascadia earthquake scenarios also fail to explain the spatial distribution of deep-seated landslides. Consequently, landslide stability models, which are strongly influenced by landscape slope, pore-water pressure, and seismic acceleration, yield landslide susceptibility maps which are broadly anti-correlated with mapped deep-seated landslide density. Through a multivariable linear regression model, we show that much of the variance in deep-seated landslide density can be explained by variability of intra-unit stratigraphic and structural characteristics, which we measure at 128 sites across two study areas totaling ∼3000 km2. Our results suggest bedding dip is only weakly correlated to landslide density, but strongly influences landslide failure style. Subtle increases in bedding dip, even in the gently folded Tyee Formation, result in a substantially higher likelihood of a landslide being cataclinal, or parallel to bedding. Overall, we find a slight majority of landslides fail within these cataclinal slopes, and that these landslides tend to be larger than non-cataclinal landslides. We also show that the lithological and structural properties that influence landslide susceptibility are distinct for these two populations of landslides. Our results demonstrate how localized, intra-unit, geologic variability can exert strong control on landslide susceptibility and failure style. This suggests that in some locations, landslide hazard models could be significantly improved by incorporating detailed, spatially variable, geologic properties rather than relying solely on generalized geologic map units.

Oregon↗

USGS2024

Hypothetical CO2 leakage into, and hydrological plume management within, an underground source of drinking water at a proposed CO2 storage facility, Kemper County, Mississippi, USA

A large Geologic Carbon Sequestration (GCS) hub has been proposed in Kemper County, Mississippi. The target injection interval consists of numerous Cretaceous-aged deep saline aquifers overlain by a competent and extensive regional sealing layer. Above the seal, the deepest Underground Source of Drinking Water (USDW) at the site is the Eutaw aquifer of the Eutaw Group and McShan Formation, undifferentiated. To assess potential risks of leakage from the deep sequestration reservoir, a model of a portion of the Cretaceous Eutaw Group was constructed in this study. Simulations tested various permeabilities, hypothetical leakage rates, and plume mitigation strategies utilizing existing wells. Results suggest that, under the influence of regional groundwater flow fields, leaking CO 2 would effectively bypass the existing wells, and to influence this migration would require very large water extraction rates. Therefore, to ensure plume detection, monitoring for leakage at the injection wells themselves is very important.

Mississippi↗

USGS2024

Analyzing spatial distributions and alignments of pitted cone features in Utopia Planitia on Mars

Martian geomorphology and surface features provide links to understanding past geologic processes such as fluid movement, local and regional tectonics, and feature formation mechanisms. Pitted cones are common features in the northern plains basins of Mars. They have been proposed to have formed from upwelling volatile-rich fluids, such as magma or water-sediment slurries. In this study, we map the spatial distributions of pitted cone features across the Utopia Planitia (UP) basin. Using the Average Nearest Neighbor technique, we find that pitted cone features appear to be clustered across the basin, occurring in a narrow band around the circumferential rim of UP. Additionally, we find that pitted cone features also appear to be aligned in chains that are sub-parallel with the basin margin. Parallel bands of cones generally follow elevation contours of the UP basin, which suggests elevation or a correlated factor plays a major role in pitted cone and cone chain formation. We propose that pitted cones and cone chains may be related to vertical fractures formed around the UP basin rim from subsidence of infilling basin material.

Icarus↗

USGS2025

Constraining large magnitude event source and path effects using ground motion simulations

The purpose of this study is to use ground motion simulations to investigate ways in which source and path effects for large magnitude events can be represented in non-ergodic GMMs. While we initially developed computation techniques using CyberShake simulations, the range of magnitudes and source-site combinations is not adequate to replicate what is observed empirically. We therefore designed a new ground motion simulation study, which includes earthquakes with a large range of magnitudes distributed uniformly on a fault plane, and sites covering a large range of rupture distances and azimuths. After running a large suite of simulations (M4-M7), we then develop a non-ergodic GMM with the simulation data. We find that the within-site residuals are dominated by the radiation pattern, rupture directivity, and slip patterns. Next, we modify an existing rupture directivity model to fit and remove the observed radiation pattern and rupture directivity from the residuals. We also minimize the contributions of slip patterns by averaging the within-site residuals among multiple source realizations. Finally, after removing the source effects from the within-site residuals, we compare the path effects computed with different magnitude groups using two approaches. The first approach only considers the small events that have the same shortest path to a site with the large events, while the second approach considers all small events on the fault plane. The results indicate that the path effects of large events cannot be satisfactorily approximated with that of small events using either approach.

Conference Paper↗

USGS2026

Neutron scattering reveals fractionation of natural gas mixtures in unconventional petroleum reservoir pores: Perspectives on energy resource recovery and storage

In unconventional petroleum reservoirs hydrocarbon fluids are hosted by both mineral and organic matter pores. These pores can have diameters that range from microns to less than a single nanometer and, for unconventional reservoirs, there is evidence that small pores ( <20 nm diameter) may constitute a large proportion of the available space. Understanding subsurface volumes and how fluids behave in them can be helpful for predicting hydrocarbon production and storage in the subsurface. One area with knowledge gaps regarding hydrocarbon behavior in small pores is the possibility for mixtures to fractionate (i.e., unmix) based on pore size or pore type. Mixture fractionation as a function of pore size could impact recovery of hydrocarbons, drive compositional shifts during production, and limit fluid storage within candidate reservoirs. To investigate natural gas fractionation in small geologic pores, we applied total neutron scattering to probe methane-ethane mixtures at reservoir pressures (up to ≈30 MPa) and temperature (60°C) within a sample from the Upper Cretaceous Niobrara Formation. Neutron scattering data reveal only minor fractionation occurs between methane and ethane in 20-nm diameter sample mesopores. Increased fractionation is observed for sample micropores, with up to 72% (±1% at 1-sigma) methane found in 2 nm diameter pores following injection of a 50%-50% methane-ethane mixture. These data provide rarely available direct experimental observations of hydrocarbon mixture behavior under nanoconfinement in a sample from an important unconventional petroleum reservoir. Our results are discussed in the context of evaluating hydrocarbon resources in unconventional reservoir meso- and micropores, reconciling observed gas composition changes during production, and more broadly, understanding subsurface pore volumes within an energy storage framework.

Fuel↗

USGS2025

Induced earthquakes are generally not tidally triggered in Oklahoma and Kansas

Human-induced earthquakes occur along critically stressed faults as injected wastewater simultaneously heightens fluid pressure and pushes faults to failure. We investigate the possibility that small stresses imposed by Earth tides could trigger earthquakes in the induced seismicity region of Oklahoma and Kansas from 2011 to 2018. We decluster a catalog consisting of ∼110,000 earthquakes using three methods (Reasenberg, nearest-neighbor distance, and phase-bin). We find no significant tidal earthquake triggering using Schuster's p -value test for the declustered catalogs as a whole. We search for localized triggering using discretized space-time cells and find ∼0–6% of cells have significant tidal triggering which is close to what is randomly expected (5%) and indicates there is an insignificant amount of tidal triggering for the full study region. One area that has significant p -values across multiple time windows, ∼2014–2016 is ∼15 km from a region of large wastewater injection volume. It is possible that localized tidal triggering occurs for this time and area because faults remain critically stressed and are particularly susceptible to slip under the small stress load from the semidiurnal tide. Possible explanations for the lack of tidal triggering in our broader study are that the pre-seismic stressing rate in the earthquake nucleation area is faster than the tidal stressing rate (∼3 kPa/day), faults are not close enough to critically stressed to be affected by tidal forcing, and that nucleation occurs over longer periods than the tides considered in this study (∼1, ∼14 days). Fluid injection could be the source of a higher pre-seismic stress rate.

Kansas, Oklahoma↗

USGS2026

Flood-inundation maps for Río Grande De Loíza in and near Caguas, Puerto Rico, 2026

Digital flood-inundation maps for a 2.7-mile reach of Río Grande De Loíza in Caguas, Puerto Rico, were created by the U.S. Geological Survey. Water-surface profiles were computed for the stream reach by using a one-dimensional, steady-state, step-backwater model. The model was calibrated to the current (2025) stage-streamflow relation (rating curve) for the U.S. Geological Survey streamgage 50055000, Río Grande De Loíza, Puerto Rico. The resulting hydraulic model was then used to compute 16 water-surface profiles for water levels (flood stages) ranging from 19.00 to 34.00 feet at the streamgage; these flood stages range from “moderate flood stage” to above “major flood stage” as defined by the National Weather Service. The 34.00-foot stage exceeds the historical maximum peak stage of 33.20 feet, recorded at the streamgage in 1945. The simulated water-surface profiles were used in combination with a digital elevation model derived from light detection and ranging (lidar) data to map the inundated areas associated with each flood profile. The flood-inundation maps and the supporting hydraulic model produced by this study can be used by emergency managers and local officials to assess flood-mitigation strategies and to define flood-hazard areas to help protect life and property, to coordinate flood-response activities such as evacuations and road closures, and to aid post-flood recovery efforts.

Puerto Rico↗

USGS2026

The effect of temperature, flow, density and disease on smallmouth bass (Micropterus dolomieu) early life growth rate

Individual growth rates of fish are influenced by gradually changing environmental conditions and rapid-onset events. However, there is a paucity of information on drivers of spatiotemporal variability in growth rates of fishes across large spatial extents. Using a 36-year dataset of smallmouth bass ( Micropterus dolomieu ) length-at-age data across 10 river reaches in Pennsylvania, USA ( n = 54,068 individuals), we estimated annual early life growth rates (mm year −1 ) and assessed the effects of relative abundance, summer temperature, summer streamflow rate and a disease outbreak on growth rate. We found that temperature had a positive effect on growth rate, varied spatially, and had a stronger effect in diseased reaches. There was also evidence for an impact of disease on early life growth through density-dependent mechanisms, with growth rates increasing from ∼85 to 101 mm·year −1 following disease-related decreases in abundance. This study adds to our understanding of the factors that shape early life growth, including rapid-onset events like disease which can have immediate and lasting effects on growth rate trajectory.

Pennsylvania↗

USGS1976

Transport and dispersion of fluorescent tracer particles for the dune-bed condition, Atrisco Feeder Canal near Bernalillo, New Mexico

A tracer technique in which mineral particles were coated with fluorescent dyes was used to study the rates of transport and dispersion of sediment particles`of various diameters and specific gravities for a dune-bed condition in an alluvial channel. The experiment was conducted in the Atrisco Feeder Canal near Bernalillo, N. Mex., between May 1 and July 14, 1967. A continuous point source of tracers, approximated by injections at 5- or 10-minute intervals, was maintained for 7 days so that the steady-dilution procedure could be used to calculate the transport rate of bed material. After termination of the injection process, the spatial-integration procedure was used to follow the movement of the tracers downstream and to calculate the transport rate. Samples of the bed material in transport and the accompanying tracers moving along the surface of the dune bed were obtained periodically throughout the study with the "dustpan" sampler especially designed for fluorescent tracer studies. In addition, the spatial distributions of the tracers in the dune bed were determined three times during the study by core sampling.

New Mexico↗

USGS2025

Reply to, “Comment on ‘The 1886 Charleston, South Carolina, earthquake: Relic railroad offset reveals rupture,’ by Roger Bilham and Susan E. Hough”

We welcome this opportunity to respond to Pratt et al. (2024) (hereinafter P24). Bilham and Hough (2023) proposed a “first-cut” elastic deformation model for the 1886 earthquake, a quantitative source model constrained by identified coseismic constraints. A key observation was the measurement of a lateral offset of a railroad line south of Summerville, leading to a model with predominately dextral slip and minor convergence, from which we concluded that active faulting had raised the Penholoway Marine Terrace >6 m since ∼770 ka. P24 questioned these constraints and proposed an alternative rupture model with predominantly reverse slip. This alternative model is neither consistent with coseismic constraints nor with other geophysical data. In a revised model presented here, we recognize that uplift of the Penholoway Terrace is confined to the eastern edge of the terrace, which we conclude results from active folding and tectonic transpression centered on the dextral fault that offset the railroad in 1886.

South Carolina↗

USGS2026

Altitude of the potentiometric surface and depth to water in the Mississippi River Valley alluvial aquifer, spring 2024

Potentiometric-surface and depth-to-water maps for spring 2024 were created for the Mississippi River Valley alluvial aquifer using groundwater-altitude data from 1,151 wells completed in the Mississippi River Valley alluvial aquifer and from the altitude of the top of the water surface in area rivers from 160 U.S. Geological Survey and U.S. Army Corps of Engineers streamgages. The potentiometric-surface and depth-to-water maps for 2024 were created to support investigations to characterize the Mississippi River Valley alluvial aquifer as part of the U.S. Geological Survey Water Availability and Use Science Program. Sufficient data were available to map the potentiometric surface and depth to water of the Mississippi River Valley alluvial aquifer for spring 2024 for about 83 percent of the aquifer area. The potentiometric contours ranged from 0 to 330 feet (ft) above the North American Vertical Datum of 1988. The regional direction of the groundwater gradient was generally to the south-southwest, except in areas of groundwater-altitude depressions, where the groundwater gradient direction was into the depression, and near rivers, where the groundwater gradient direction was either from the aquifer to the river or from the river into the aquifer. Groundwater depressions in the potentiometric-surface map are in the lower one-half of the Cache region and in most of the Grand Prairie and Delta regions. Depth to water by well in the Mississippi River Valley alluvial aquifer in spring 2024 ranged from 0.27 ft above land surface to 145.65 ft below land surface.

Arkansas, Illinois, Kentucky, Louisiana, Mississip↗

USGS2024

Will there be water? Climate change, housing needs, and future water demand in California

Climate change in California is expected to alter future water availability, impacting water supplies needed to support future housing growth and agriculture demand. In groundwater-dependent regions like California's Central Coast, new land-use related water demand and decreasing recharge is already stressing depleted groundwater basins. We developed a spatially explicit state-and-transition simulation model that integrates climate, land-use change, water demand, and groundwater gain-loss to examine the impact of future climate and land use change on groundwater balance and water demand in five counties along the Central Coast from 2010 to 2060. The model incorporated downscaled groundwater recharge projections based on a Warm/Wet and a Hot/Dry climate future from a spatially explicit hydrological process-based model. Two urbanization projections from a parcel-based, regional urban growth model representing 1) recent historical and 2) state-mandated housing growth projections were used as alternative spatial targets for future urban growth. Agricultural projections were based on recent historical trends from remote sensing data. Annual projected changes in groundwater balance were calculated as the difference between land-use related water demand, based on historical estimates, and climate-driven recharge plus agriculture return flows. Results indicate that future changes in climate-driven groundwater recharge, coupled with cumulative increases in agricultural water demand, result in overall declines in future groundwater balance, with a Hot/Dry future resulting in cumulative groundwater decline in all but Santa Cruz County. Cumulative declines by 2060 are especially prominent in San Luis Obispo (−2.9 to −5.1 Bm 3 ) and Monterey counties (−6.5 to −8.7 Bm 3 ), despite limited changes in agricultural water demand over the model period. These two counties show declining groundwater reserves in a Warm/Wet future as well, while San Benito and Santa Barbara County barely reach equilibrium. These results suggest future groundwater supplies may not be able to keep pace with regional demand and declining climate-driven recharge, resulting in a potential reduction in water security in the region. However, our county-scale projections showed new housing and associated water demand does not conflict with California's groundwater sustainability goals. Rather, future climate coupled with increasing agricultural groundwater demand may reduce water security in some counties, potentially limiting available groundwater supplies for new housing.

California↗

USGS2025

Bathymetric and acoustic-backscatter mapping of Lake Sammamish, Washington, during USGS field activity 2021-656-FA

In 2021, scientists from the U.S. Geological Survey (USGS) acquired bathymetric and acoustic-backscatter data within Lake Sammamish, Washington. Mapping was completed to find evidence of past earthquakes such as underwater landslides and is part of a larger USGS project to understand the overall geologic hazards history of the Cascadia Margin region. The survey was conducted using the USGS research vessel Parke Snavely outfitted with an interferometric sidescan-sonar system for swath mapping. Three estimates of bathymetric uncertainty were calculated from the post-processed bathymetric data. This report provides a summary of the mapping mission.

Washington↗