USGS Science⌕ Search

SEARCH · USGS Science

Results for “Climate of the Past”

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,687 records · Page 79Linked to original sources
USGS2025

VIPER site analysis

We needed to evaluate available orbital data of NASA’s Volatiles Investigating Polar Exploration Rover (VIPER) mission area in order to derive a variety of maps to help the science team identify scientifically interesting places for the rover to visit and to provide scientific context for our mission. Some of these maps also fulfilled engineering and mission design needs to enable safe and efficient landing and roving. We incorporated data from the Lunar Reconnaissance Orbiter Camera, the Lunar Orbital Laser Altimeter, the Mini-RF instrument, the Chandrayaan-2 Orbital High Resolution Camera, the Korean Pathfinder Lunar Orbiter’s Shadowcam, the Kaguya Spectral Profiler and Multiband Imager, and the Chandrayaan-1 Moon Mineralogy Mapper. We used a variety of techniques to build these maps, including stereogrammetry, shape-from-shading, ice stability depth and surface temperature calculations, and the horizon method for solar illumination and direct-to-Earth communications maps. Altogether, these maps allowed us to survey for boulders, evaluate features in permanently shadowed regions that VIPER might explore, provide mineralogic context for what VIPER’s instruments may learn, estimate the ages and radar properties of craters in the VIPER mission area, and evaluate the potential for gravity traverses with the rover. These data and techniques provided a rich set of information from which both the VIPER science team and engineering teams were able to draw in order to plan a safe landing and to plan a VIPER surface mission that will be both scientifically valuable and robust from an operational perspective.

The Planetary Science Journal↗
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↗
USGS2025

Laboratory assessment for recovery of porcine circovirus 2 and porcine reproductive and respiratory syndrome virus using two types of commercially available hollow-fiber ultrafilters

Groundwater near swine farms is an uninvestigated reservoir for porcine reproductive and respiratory syndrome virus (PRRSV) and porcine circoviruses (PCVs). Enteric microorganisms are often collected from groundwater via dead-end ultrafiltration, but recovery of PRRSV and PCV with this method has not been assessed. We recovered PRRSV2 and PCV2 by dead-end ultrafiltration followed by polyethylene glycol (PEG) precipitation, nucleic acid extraction, and reverse-transcription quantitative real-time PCR. We also compared 2 commercial hemodialysis ultrafilters (Asahi Kasei Rexeed-25A, Nipro Elisio-25H) and compared PRRSV2 recovery in these filters to other waterborne microorganisms. On average, 8 ± 1% of PRRSV2 was recovered by dead-end ultrafiltration and PEG precipitation, compared to 25 ± 6% for adenovirus 41. Full-process recovery of bacteria in the same filters was 5–15%; Cryptosporidium parvum recovery was 42 ± 12%. PCV2 was detected in 4 of 12 replicate filters, but low stock concentrations precluded quantitative recovery estimates. Elisio-25H ultrafilters performed similarly to Rexeed-25A filters for all organisms tested and is an effective replacement for the Rexeed-25A, which is no longer available in the United States. Our recovery of PRRSV2 and PCV2 by dead-end ultrafiltration in the laboratory suggests that PRRSV2 detection limits are as low as 3–50 genomic copies/L in sample volumes of 100–1,500 L. Based on quantitative microbial risk assessment, these concentrations are relevant to PRRSV2 infection rates in the U.S. swine herd.

Journal of Veterinary Diagnostic Investigation↗
USGS2024

Uncertainty and spatial correlation in station measurements for mb magnitude estimation

The body‐wave magnitude (⁠⁠) is a long‐standing network‐averaged, amplitude‐based magnitude used to estimate the magnitude of seismic sources from teleseismic observations. The U.S. Geological Survey National Earthquake Information Center (NEIC) relies on in its global real‐time earthquake monitoring mission. Although waveform modeling‐based moment magnitudes are the modern standard to characterize earthquake size, is important because (1) in many cases, waveform modeling is not possible (e.g., low signal‐to‐noise events), (2) is applicable over a broad range of magnitudes, ∼M 4–7, and (3) there is a many decades‐long history of estimating magnitudes. We use the NEIC Preliminary Determination of Epicenters earthquake catalog to investigate the uncertainty in NEIC station measurements. We show that measurements are spatially correlated, which can bias event ⁠, and we describe an empirical relation between this spatial correlation and station‐to‐station distance. We further describe an approach to mitigate bias from the spatial correlation. Accounting for the spatial covariance of observations can change the event from −0.15 to 0.07 units (10th to 90th percentile) for smaller events (⁠⁠). These smaller events have the largest standard deviations ranging from 0.05 to 0.15 units (10th to 90th percentile).

The Seismic Record↗
USGS2026

Streamflow, water quality, and constituent loads and yields, Scituate Reservoir drainage area, Rhode Island, water year 2023

As part of a long-term cooperative program to monitor water quality within the Scituate Reservoir drainage area, the U.S. Geological Survey, in cooperation with Providence Water (formerly the Providence Water Supply Board), collected streamflow and water-quality data in tributaries to the Scituate Reservoir, Rhode Island. Streamflow and concentrations of chloride and sodium estimated from records of specific conductance for 16 tributaries were used to calculate loads of chloride and sodium during water year 2023 (October 1, 2022, through September 30, 2023). Water-quality samples were collected by Providence Water at 38 sampling stations on tributaries to the Scituate Reservoir during water year 2023. These water-quality data are summarized by using values of central tendency and are used, in combination with measured (or estimated) streamflows, to calculate loads and yields of selected water-quality constituents for water year 2023. Annual mean streamflows for monitoring stations in this study ranged from 0.47 to 34.4 cubic feet per second during water year 2023. At the 16 continuous-record streamgages, tributaries transported about 3,000 metric tons of chloride and 1,900 metric tons of sodium to the Scituate Reservoir; annual chloride yields for the tributaries ranged from 19 to 200 metric tons per square mile, and annual sodium yields ranged from 13 to 120 metric tons per square mile. At the stations where water-quality samples were collected by Providence Water, the medians of the median daily loads were 65,000 million colony forming units per day for coliform bacteria, 870 million colony forming units per day for Escherichia coli , 200 kilograms per day for chloride, 8.2 grams per day as nitrogen for nitrite, 520 grams per day as nitrogen for nitrate, and 300 grams per day as orthophosphate for phosphate. The medians of the median yields were 25,000 million colony forming units per day per square mile for coliform bacteria, 600 million colony forming units per day per square mile for Escherichia coli , 65 kilograms per day per square mile for chloride, 4.6 grams per day per square mile as nitrogen for nitrite, 180 grams per day per square mile as nitrogen for nitrate, and 180 grams per day per square mile as orthophosphate for phosphate.

Rhode Island↗
USGS2026

Morphodynamics of snow barchans at Concordia Station, Antarctica

Barchans are crescent-shaped dunes that form and migrate under unimodal wind regimes in areas of low sediment availability. While sand barchan morphodynamics has been extensively studied, snow barchans remain poorly documented. Here, using hourly time-series images taken near Concordia Station, Antarctica, we measure the sizes (lengths 0.3 m to 13 m, widths 0.4 m to 7.5 m) and migration rates (0.13 m h −1 to 3.5 m h −1 ) of 44 populations of snow barchans during 2017–2022. We compare this dataset with established height–velocity and length-width relationships for sand barchans in both air and water. Our results show that snow barchans follow the same scaling as sand barchans, allowing us to estimate a characteristic sintering (bonding of snow particles) time of about one day. This sintering process limits the size of snow barchans when active, reduces their aspect ratios, and, once wind falls below the transport threshold, explains their preservation over a few years. During wind events, migrating barchans are too small to maintain their shape and orientation, making them highly responsive to fluctuations in wind direction. Thus, studying snow barchans improves understanding of an elementary bedform encountered in multiple environments and also of the surface properties of snow on different timescales, offering insights into snow accumulation and wind-blown snow in ice cores.

Geomorphica↗
USGS2026

An entropic explanation for Gutenberg-Richter scaling

We develop a simple explanation for Gutenberg-Richter (G-R) size scaling of earthquakes on a single fault. We discretize the fault and consider all possible contiguous ruptures at that level of discretization. In this static model, we assume that slip scales with rupture length, and that the rupture rates at each point along the fault are consistent with an a priori long-term slip rate. These simple assumptions define an (under-determined) non-negative least-squares inverse problem. Each solution to this inverse problem is a set of earthquake rates that matches the slip-rate constraint. We use a Markov Chain Monte Carlo (MCMC) algorithm to uniformly sample the solution space assuming constant slip rates along the fault. At finer discretizations, deviations from G-R behavior decrease, which is consistent with an entropic pressure towards G-R solutions. When the fault is discretized into 10 or more segments, random solutions found by the MCMC algorithm have G-R size scaling, even though there are trivial solutions that, for example, have earthquakes of only one size. This is because there are simply far more solutions that have G-R scaling; as the problem size increases, the strong degeneracy of GR solutions results in other solutions becoming improbably rare. Also, the entropically favored G-R distribution has a b -value of approximately 1, which agrees with measured b -values in real earthquake catalogs.

JGR Solid Earth↗
USGS2025

Dynamic feedbacks between river meandering and landsliding in northwestern Washington glacial terraces

Landsliding in river valleys poses unique risks for cascading hazards and can damage infrastructure and cause fatalities. In postglacial valleys, many landslides are posited to occur in relation to lateral river erosion, but the dynamics of fluvial-hillslope interactions are not well understood. Here, we investigate a section of the Nooksack River in western Washington State where the channel is flanked by landslide-prone glacial terraces similar to those that failed in the 2014 State Route 530 “Oso” landslide. We map 216 landslides through time across 17 aerial imagery data sets (1933–2022) and analyze them in relation to river meandering and curvature. We observe dynamic feedbacks between lateral river meandering and valley-adjacent landsliding. Terrace lateral retreat rates of up to 25 m/year owing to combined fluvial erosion and slope failure occur on pinned, outer meander bends immediately downstream from peaks in river curvature (>0.0075 1/m); these locations are predisposed to both shallow and deep-seated landslides. Deep-seated landslides extending 17%–32% of the active valley width into the floodplain can displace the river away from the floodplain margin and change the channel planform. River-displacing landslides relocate meanders up- or downstream, thereby conditioning the location of subsequent landslides. This conceptual model of coupled landslide-driven meander displacement and valley-adjacent landsliding is exemplified across western Washington river systems. The distance between up- and downstream valley-adjacent landsliding scales with valley width, meander wavelength, and terrace height. Our results can advance our understanding of the river-hillslope interface in landscape evolution and can be used to inform hazard management in river corridors.

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

Completion summary for monitor wells NRF-17 and NRF-18 at the Naval Reactors Facility, Idaho National Laboratory, Idaho

The U.S. Geological Survey (USGS)—in cooperation with the U.S. Department of Energy (DOE) for the Naval Reactors Laboratory Field Office that supports operations for the Naval Reactors Facility (NRF) located at the Idaho National Laboratory (INL)—drilled and constructed well NRF-17 (formerly borehole USGS 151) and well NRF-18 (formerly borehole USGS 152) for stratigraphic framework analyses and water-quality monitoring at the Idaho National Laboratory (INL) near the NRF, in southeastern Idaho. Borehole USGS 151 was continuously cored from about 48 to 1,070 feet (ft) below land surface (BLS); rotary drilled from approximately 1,070 to 1,720 ft BLS; and re-drilled to complete construction as a monitor well NRF-17, completed to 461 ft BLS. Borehole USGS 152 was continuously cored from approximately 19 to 1,259 ft BLS; rotary drilled from approximately 1,259 to 1,630 ft BLS; and re-drilled to complete construction as a monitor well NRF-18, completed to 450 ft BLS. Geophysical data were examined with photographed core material to record lithologic descriptions and to suggest zones where groundwater flow was anticipated. Basalt flows varied from highly fractured to dense, with high-to-low vesiculation. Well NRF-17 generally was constructed in mostly dense basalt (greater than 75 percent), and well NRF-18 was constructed in primarily fractured and (or) vesicular basalt. In well NRF-17, the well capacity is directly affected by the limited amount of fractured basalt, which serves as the primary pathway for groundwater. This effect was observed during the pumping test conducted after the well's final construction. Single-well aquifer tests were done at wells NRF-17 and NRF-18 to provide estimates of transmissivity and hydraulic conductivity after final well construction and initial well development. Estimated values of transmissivity and hydraulic conductivity for well NRF-17 were 8.81 feet squared per day (ft 2 /d) and 1.04×10 -2 feet per day (ft/d), respectively. Estimated values of transmissivity and hydraulic conductivity for well NRF-18 were 4.77×10 3 ft 2 /d and 5.61 ft/d, respectively. The NRF-17 pump test resulted in 19.41 ft of measured drawdown at a sustained average pumping rate of 3.3 gallons per minute (gal/min). The NRF-18 pump test resulted in 0.55 ft of measured drawdown at a sustained average pumping rate of 31.0 gal/min. Water-quality samples collected from the two wells were analyzed for cations, anions, metals, nutrients, volatile organic compounds, stable isotopes, and radionuclides. Water samples for select inorganic constituents showed concentrations consistent with signatures from tributary valley groundwater with influences from ephemeral surface-water recharge from the Big Lost River. Water-quality samples analyzed for stable isotopes of oxygen and hydrogen are consistent with signatures from tributary valley groundwater and surface-water recharge inputs to the aquifer. No measured water-quality results were greater than their respective maximum contaminant levels for public drinking-water supplies. Inorganic and nutrient water-quality results for well NRF-17 and well NRF-18 suggest the groundwater in this area is potentially affected by industrial wastewater disposal.

Idaho↗
USGS2025

Groundwater salinity: Applying the specific conductance and water type proxy

Groundwater is increasingly needed for water supplies but may have limited utility in some locations because of its salinity. Salinity, often expressed as total dissolved solid (TDS), is frequently estimated using specific conductance (SC) measurements. However, the commonly used proxy (0.65 multiplied by SC to indicate TDS, common in many handheld meters) can result in inaccurate TDS estimates. First, the TDS–SC relationship is not linear over the entire concentration range of groundwater. Furthermore, the TDS (and salinity)–SC relationships vary substantially depending on the major-ion composition. Here we develop a proxy method utilizing SC and major-ion water type to estimate TDS and salinity specifically for groundwaters. Compared to most surface waters, groundwater tends to have a wider range of salinity (fresh to highly saline) and higher concentrations of bedrock-derived solutes such as carbonate ions, silica, and many other ions. The dataset used to develop the proxies includes water chemistry data from 149,059 discrete groundwater samples. The groundwater proxies, which employ nonlinear log–log relations, utilize five water types (HCO 3 , Cl, Ca-Mg-SO 4 , Na-K-SO 4 , and mixed waters), are accurate (median percent difference between TDS and salinity determined using the proxy compared to discrete measurements was <±0.8%) over a wide range of SC (up to 200 mS/cm), rapid, cost-effective, and can be measured on-site.

Groundwater↗
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↗
USGS1998

Trace elements in streambed sediment and fish liver at selected sites in the Upper Colorado River Basin, Colorado 1995-96

Trace elements were analyzed in streambed-sediment samples collected from 16 sites and in fish-liver samples collected from 14 sites in the Upper Colorado River Basin in Colorado as part of the National Water-Quality Assessment program. Sites sampled represented agricultural, mining, mixed, and urban/recreation land uses and background conditions. The results for 15 trace elements in streambed-sediment and in fish-liver samples are presented in this report. Fourteen of the selected trace elements were detected in streambed-sediment samples collected at all sites. Twelve of the selected trace elements were detected in fish liver at more than 50 percent of the sites. Cadmium, copper, selenium, and zinc were selected for a more detailed analysis. Cadmium, copper, and zinc concentrations in streambed sediment were highest at mining land-use sites in the Southern Rocky Mountains physiographic province. Selenium concentrations in streambed sediment were highest at an agricultural land-use site in the Colorado Plateau physiographic province. The concentration of trace elements in streambed sediment generally increased as particle size decreased. Concentrations of trace elements in fish liver generally did not follow the same relation to land use as concentrations in streambed sediment; however, cadmium concentrations in fish liver were highest at a mining land-use site in the Southern Rocky Mountains physiographic province, and selenium concentrations in fish liver were highest at an agricultural land-use site in the Colorado Plateau physiographic province. Copper and zinc concentrations in fish liver were highest at mixed land-use sites. Comparison of streambed-sediment and fish-liver concentrations to two other similar NAWQA studies in the Rocky Mountain region generally indicated similar patterns in relation to land use for streambed sediment, but not for fish liver. Cadmium, copper, and zinc concentrations in streambed sediment were highest at sites affected by mining in all three study units. Selenium concentrations in streambed sediment did not indicate relations among the three study units when compared to land use. Cadmium in fish liver was highest at sites affected by mining in all three study units. Copper, selenium, and zinc in fish liver did not indicate relations among the three study units when compared to land use.

Colorado↗
USGS2024

Classification of lakebed geologic substrate in autonomously collected benthic imagery using machine learning

Mapping benthic habitats with bathymetric, acoustic, and spectral data requires georeferenced ground-truth information about habitat types and characteristics. New technologies like autonomous underwater vehicles (AUVs) collect tens of thousands of images per mission making image-based ground truthing particularly attractive. Two types of machine learning (ML) models, random forest (RF) and deep neural network (DNN), were tested to determine whether ML models could serve as an accurate substitute for manual classification of AUV images for substrate type interpretation. RF models were trained to predict substrate class as a function of texture, edge, and intensity metrics (i.e., features) calculated for each image. Models were tested using a manually classified image dataset with 9-, 6-, and 2-class schemes based on the Coastal and Marine Ecological Classification Standard (CMECS). Results suggest that both RF and DNN models achieve comparable accuracies, with the 9-class models being least accurate (~73–78%) and the 2-class models being the most accurate (~95–96%). However, the DNN models were more efficient to train and apply because they did not require feature estimation before training or classification. Integrating ML models into benthic habitat mapping process can improve our ability to efficiently and accurately ground-truth large areas of benthic habitat using AUV or similar images.

Michigan, Wisconsin↗
USGS2024

A methodology to estimate CO2 and energy gas storage resources in depleted conventional gas reservoirs

Depleted hydrocarbon reservoirs are subsurface geological structures capable of sequestering vast quantities of carbon dioxide (CO 2 ) as well as storing other energy gases for later usage, such as natural gas, and potentially hydrogen (H 2 ). Here we outline a methodology to quantify multi-gas storage resources in depleted conventional gas reservoirs for usage in assessments by the United States Geological Survey (USGS) at the scale of sedimentary basins. The methodology consists first of quantifying accessible pore volume in a depleted reservoir for natural gas storage using up to three equations. Input data are derived from commonly reported or estimated reservoir parameters and natural gas production volumes, and equations may be combined in linear models to improve pore volume estimates. Storage estimates from these equations are tested and validated for 31 reservoirs in the Michigan Basin Province, USA that were previously converted to underground gas storage facilities and have known (federally reported) natural gas storage capacities. Secondly, natural gas storage capacities can be transformed via fluid substitution calculations to estimate the storage resources for non-native fluids, applied here for, CO 2 , H 2 , and methane-H 2 blends, accounting for molecule-specific deviations from ideal gas behavior at reservoir pressures and temperatures as well as differing storage efficiencies. Importantly, the storage of non-native fluids may not be appropriate in all depleted gas reservoir settings due to potential risks like leakage, in particular in the case of H 2 storage, requiring additional knowledge of caprock sealing capacity. Given this caveat, we demonstrate the fluid substitution method for natural gas reservoirs of the Northern Niagaran Reef and Southern Niagaran Reef USGS plays in the Michigan Basin Province, as these trends of Silurian pinnacle reefs are capped with tight-sealing evaporite facies. The deterministic equations outlined from this methodology can be incorporated into future probabilistic USGS gas storage assessments for CO 2 , H 2 , and natural gas resources in the United States.

Michigan↗
USGS2025

A method to obtain remotely sensed grain size distributions from nonplanar granular deposits

Constraining the grain size distribution of granular deposits with complex surfaces is difficult with existing approaches. Field and laboratory techniques are time consuming and limited by the maximum grain size that laboratories can accommodate. In this study, we present a new method to identify the coarse fraction of the grain size distribution at a debris-flow fan deposit surveyed with terrestrial laser scanning (TLS) in Glenwood Canyon, Colorado, USA. This method is a novel grain segmentation algorithm developed for application to point cloud data of deposits with complex surfaces and angular grains ranging in size from centimeters to a meter. This approach combines an existing random forest machine learning method with a novel iterative clustering algorithm. We compared the grain size distribution from our algorithm with a Wolman pebble count conducted in the field, and found a root mean squared error of less than 2 cm from the 5th to 95th percentile of the grain size distribution of grains ranging from cobble to boulder sized (6.3–78 cm in our application). Finally, we compared our new algorithm with an existing open-source grain segregation algorithm, and our method outperformed the selected alternative when applied to the debris-flow deposit point cloud.

Colorado↗
USGS2026

Riverscape heterogeneity shapes population diversity for a migratory fish

Habitat patch dynamics can scale up to influence population demography and diversity with implications for resilience to environmental stochasticity. But how the spatial arrangement and size of habitat patches interact with other components of habitat heterogeneity to shape population diversity at larger spatial scales is not well understood. For riverine fishes, there is increasing evidence that tributary streams provide critical demographic support to main stem rivers. However, the extent to which main stem rivers rely on demographic contributions from tributaries, and the factors underlying this dependence, have not been assessed. Here, we used genetic stock identification to evaluate the effect of tributaries on population diversity of Yellowstone cutthroat trout ( Oncorhynchus virginalis bouvieri ) occupying the main stem Snake River, Wyoming, USA. We found that the main stem relied almost entirely on tributaries for demographic support, but main stem composition varied spatially among river sections. Distance between habitat patches, catchment area, and groundwater availability acted in concert to determine the contribution of specific tributaries to the main stem, but contributions were ultimately modulated by habitat connectivity. We also found evidence for multi-scale spatial structure in tributary contributions, providing insight into untested drivers of main stem river population diversity. Our results demonstrate how spatially discrete and distributed riverscape attributes influence population diversity at broader spatial scales, illustrating how ecosystem resilience emerges from the dynamic, two-way exchange of individuals and energy across habitat networks. Management plans for large rivers that address the ecological contributions of tributaries may be needed to achieve optimal outcomes. Similarly, conservation strategies that exclusively focus on headwater streams may fail to capture the broader habitat requirements necessary to maintain robust cold-water fish populations and associated recreational fisheries, particularly under global environmental change.

Wyoming↗
USGS2021

Critical review of mercury methylation and methylmercury demethylation rate constants in aquatic sediments for biogeochemical modeling

Mercury is a toxin that causes neurological impairments in adults, is particularly harmful for fetuses and children, and is deadly in severe cases, making it a worldwide health concern. Methylmercury (MeHg) is the environmentally relevant form of mercury (Hg) because it biomagnifies along the food chain. Methylmercury is mainly produced in aquatic sediments via methylation of inorganic Hg (Hg(II)) and transformed back via demethylation. Because transformation rates determine MeHg concentrations, quantification of methylation and demethylation rates is needed to inform management of MeHg. Published rate constants for Hg(II) methylation ( 𝑘 𝑚 ) and MeHg demethylation ( 𝑘 𝑑 ) vary greatly, stemming partly from differences in experimental methods. We conducted a comprehensive review of rate laws, evaluated published rate constants, and performed biogeochemical simulations to assess variability in reported 𝑘 𝑚 and 𝑘 𝑑 . Based on selected studies employing the same pseudo-first-order rate law and similar experimental methods, we found that 𝑘 𝑚 = 0.04 ± 0.03 d −1 is a reasonable range for wetland sediments. Over a number of environments, maximum 𝑘 𝑑 was smaller at sites without Hg source ( 𝑘 𝑑 = 0.5 d −1 ) than at sites with identified Hg source ( 𝑘 𝑑 = 1.8 d −1 ). Larger variability and higher uncertainty in 𝑘 𝑑 compared to 𝑘 𝑚 highlight the need for more research on MeHg demethylation rates. This critical review: (a) aids the design of future experimental studies of 𝑘 𝑚 and 𝑘 𝑑 ; (b) provides guidance for comparing rate constants from different studies; (c) presents a biogeochemical reaction model to assess rate constants; and (d) informs selection of 𝑘 𝑚 and 𝑘 𝑑 values from the literature for use in model simulations.

Critical Reviews in Environmental Science and Tech↗