Search USGS⌕ Search

SEARCH · Search USGS

Search Search USGS

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.

At least 1,135 records · Page 63Linked to original sources
USGS2026

Environmental DNA pilot monitoring program for invasive species and biodiversity assessments on Santa Cruz Island: Interim report, September 2025

The U.S. Geological Survey (USGS) and Southern California Coastal Water Research Project supported Channel Islands National Park, The Nature Conservancy’s (TNC) Santa Cruz Island Preserve, and University of California San Diego (UCSD) researchers in using environmental DNA sampling to monitor for invasive Argentine ant ( Linepithema humile ) and to describe spatial biodiversity patterns in a freshwater stream on Santa Cruz Island that is colloquially known as Cañada del Puerto Creek (hydrologic unit code 12 number 180600140201). An eDNA autosampler was deployed May 27–June 02, 2025, which filtered and preserved four 2-L water samples every 12 hours (n = 48 samples). These samples were extracted by USGS and shipped to UCSD for targeted eDNA analyses of Argentine ant. Environmental DNA sampling kits were used to filter and preserve water samples at four sites along the creek spanning ~ 5 km of stream distance. Six 1-L samples were collected at each site May 27–28 and samples were analyzed for eDNA metabarcoding using a comprehensive panel, which encompasses species across the tree of life. In these samples, eDNA from 2,134 unique taxa spanning 30 taxonomic groups (e.g., amoebas, insects, plants) were found. These sampling efforts demonstrated how eDNA autonomous and manual sampling approaches can be applied by National Park Service, TNC, and partners at scale to provide otherwise difficult and expensive to acquire information about biological threats and biodiversity.

California↗
USGS2026

Origin of the Pd/Pt ratio of the J-M Reef, Stillwater Complex, Montana, USA

The J-M reef of the Stillwater Complex is characterized by a high Pd/Pt ratio (mean ~3.8 with a standard error of 0.03) with a homogeneous geospatial distribution at the deposit scale. In this contribution, we demonstrate that the Pd/Pt ratio of the reef is the product of equilibration of an immiscible sulfide liquid with a silicate melt rich in Pd relative to Pt. Despite the high tenors of the J-M reef sulfides (avg 2,700 ppm Pt and 770 ppm Pt), numerical modeling shows that the parental melts did not have extraordinary Pd and Pt concentrations. Instead, the initial composition of a plausible parental silicate melt can have Pd and Pt contents well within the expected range of a normal, mantle-derived partial melt (i.e., ~10–20 ppb for both Pd and Pt with Pd/Pt of ~1). The relative differences in the partitioning behavior of Pt and Pd between sulfide liquid and silicate melt are unlikely to produce a consistent Pd/Pt ratio across a wide range of silicate melt to sulfide liquid mass ratios (i.e., R factors). Instead, the pre-emplacement fractionation of Pt alloy from S-undersaturated silicate magma accounts for the homogeneous and high Pd/Pt ratio of the J-M reef. We show that batch equilibration of sulfide liquid with silicate melt can produce the high Pd/Pt ratios of the reef if the partition coefficients between sulfide liquid and silicate melt for Pd and Pt are extremely high (>10 6 ). In an alternative model, Pd enrichment could be achieved by sulfide upgrading in resident footwall mush even if the partition coefficients between sulfide liquid and silicate melt are relatively small (between 10 4 and 10 6 ) because the instantaneous mass ratio of silicate melt to sulfide liquid is small (R ≈ 100–700), so the partitioning behavior of Pt and Pd has little impact on the composition of sulfide liquid.

Montana↗
USGS2026

Frameworks for assessing tsunami hazard and risk

Tsunamis are multiscale phenomena resulting from a water column displacement that may be induced by multiple sources, and range from local scale inundation processes to ocean-wide scale wave propagation. Different strategies may be required to model tsunami evolution at different scales and to characterize various intensity measures. Research in tsunami hazard and risk has focused mostly on the tsunami effects such as the wave heights or flow depths. This chapter reviews the evolution of tsunami hazard and risk assessment methodologies, with particular emphasis on the development of probabilistic approaches. Building on advances in numerical modeling and uncertainty analysis, two main frameworks for Probabilistic Tsunami Hazard and Risk Analysis (PTHA/PTRA) are described. Framework 1 (FW1) focuses on quantitative methods, including fully simulation-based assessments (FW1A), integration of hazard with vulnerability and loss models (FW1B), consideration of dynamic processes such as tidal and sea-level variations (FW1C), and approaches using limited scenario sets (FW1D). Framework 2 (FW2) complements this by incorporating indicator-based vulnerability assessments, both physical (FW2A) and social, multi-dimensional (FW2B).

Book chapter↗
USGS2026

Prolonged cooling of andesitic-dacitic lava flows produces optimal groundmass material for 40Ar/39Ar dating

Organizing the eruptive histories and petrogenetic models that underpin hazard assessments at active volcanoes typically requires high-precision ages for late Pleistocene to Holocene lava flows. In this study, features of an existing 40Ar/39Ar dataset for andesitic-dacitic rocks from Ruapehu volcano, Aotearoa New Zealand, were interrogated with new microanalytical data to identify the optimal groundmass material for age determination. To understand the behavior of K and Ar within the cooling melt of lavas during their emplacement, we examined and compared the petrological features of samples from the ice-chilled margin and interior of a lava flow. The groundmass of the rapidly cooled sample comprises non-vesicular rhyolitic glass and microlites of plagioclase, orthopyroxene, and magnetite. The sample from the flow interior has minimal amounts of groundmass glass, and interstices are instead occupied by sanidine and tridymite. Modelling supports the interpretation that sanidine crystallized as a stable groundmass phase during slow cooling of degassed melt in the interior zones of lavas that were insulated behind their glassy margins. Our evaluation of the age data shows that sanidine contributed high yields of radiogenic argon during step-heating experiments for samples collected from lava interiors, which enabled high-precision 40Ar/39Ar ages to be determined. Using Ruapehu as an example of an arc volcano affected by lava-ice interaction, we provide a basic field guide for identifying lava flows that contain groundmass sanidine at the expense of glass. Target outcrops also offer opportunities for reconstructing geomagnetic dynamics and volcano–climate relationships during the Pleistocene and Holocene.

Ruapehu volcano↗
USGS2026

Drinking water arsenic, urinary arsenic biomarkers, and cognitive impairment in the REGARDS study

Background There are several pathways by which inorganic arsenic (iAs) exposure can affect cognition among adults. Few epidemiologic studies evaluate iAs in water and inter-individual differences in urinary arsenic toxicokinetics. We aimed to estimate the association between individual-level urinary arsenic biomarkers, county-level iAs in drinking water, and cognitive impairment in a cohort of Black and White Americans. Methods We evaluated the association between county-level iAs in water and urinary iAs metabolites with incident cognitive impairment in REasons for Geographic and Racial Differences in Stroke (REGARDS). Participants were linked to county-level concentrations of arsenic in public water (n = 15,516) and county-level probabilities of private well arsenic exceeding 10 μg/L (n = 20,448). In addition, urinary concentrations of iAs, monomethylarsonic acid (MMA) and dimethylarsinic acid (DMA) were measured in a subset of participants (n = 1,013). Cognitive impairment was determined by the Six Item Screener (SIS) and the Enhanced Cognitive Battery (ECB). We used multivariable logistic regression to estimate odds ratios. Results Approximately 9 % of participants in REGARDS developed incident cognitive impairment on the SIS over 10 years of follow-up. County-level public water iAs concentrations and private well iAs probabilities were not associated with increased odds of cognitive impairment. Higher concentrations of MMA (OR: 1.74 [95 % CI: 1.22, 2.49]) and iAs (OR: 1.58 [95 % CI: 1.12, 2.22]) in urine were associated with greater odds of cognitive impairment on the SIS. Results for the ECB were similar. Conclusions iAs exposure may be associated with cognitive impairment. Our findings highlight the importance of considering multiple measures of iAs exposure.

Environmental Research↗
USGS2026

Hakalau’s moving castle: How climate change and restoration are shifting an island fortress for forest birds

Hakalau Forest Unit of the Big Island National Wildlife Refuge Complex (hereafter, Hakalau) protects the largest area with the highest endemic forest bird diversity in Hawaiʻi, including four federally listed species. Hakalau’s higher elevation montane forest provides refuge from avian malaria ( Plasmodium relictum ), a primary driver of Hawaiian honeycreeper extinctions. However, recent declines in Hakalau’s birds at lower elevations could indicate that conditions have become suitable for disease vector Culex quinquefasciatus . We evaluated the statuses of Hakalau’s bird populations in the context of recent climatic changes using new survey data from point-transect distance sampling, producing abundance estimates from 1999 to 2024. We stratified our analysis across four elevation ranges (<1500 m, 1500–1700 m, 1700–1900 m, and >1900 m) and assessed trends for each species using state-space models (SSMs). We constrained population trajectories to be biologically realistic by incorporating population dynamic models within the SSMs. We observed highly species-specific abundance trends below 1500 m, predominantly stable to upward trends within 1500–1700 m, stable trends within 1700–1900 m, and upward trends above 1900 m. Declines in Hawaiʻi ʻamakihi ( Chlorodrepanis v. virens ) and endangered ʻakiapōlāʻau ( Hemignathus wilsoni ) abundance coincided with lengthening warm seasonal temperatures indicative of shrinking disease-free habitat below 1700 m. Above 1900 m, however, increases in nearly all species indicate that reforestation has likely restored disease-free habitat since 1999. While most species were stable to increasing overall, surveillance for mosquitoes and disease at lower elevations, documenting changes in habitat, and continuing bird population monitoring can help to gauge their long-term persistence at Hakalau.

Hawaii↗
USGS2026

A regional simulation modeling framework for evaluating invasive annual grass management across the sagebrush biome

Invasive annual grasses (IAG) continue to spread within the sagebrush biome of the western United States, degrading plant communities and wildlife habitat, decreasing forage for ranching livelihoods, and heightening wildfire risk. Effective management of IAGs requires action and long-term strategic planning across the sagebrush biome, but the cumulative effects of IAG treatments over time and space are not well understood, especially over broad extents defined for strategies like the Sagebrush Conservation Design. We developed a simulation model and sampling framework that allow local-scale actions to be ‘scaled up’ to evaluate large-scale regional and biome-wide management strategy outcomes. We worked with natural resource managers and experts to co-develop a spatially explicit state-and-transition simulation model of IAG dynamics in sagebrush landscapes that can be used to evaluate alternative management strategies. We evaluated our framework by contrasting two baseline scenarios in terms of their long-term effects on the sagebrush biome. We show that focusing management efforts on moderate to high IAG cover was effective at reducing full conversion to IAGs but failed to prevent widespread establishment of IAGs in core sagebrush areas, exposing them to increased risk of wildfire and wildlife habitat degradation. The results of our model help quantify the extent of the problem that IAGs pose to sagebrush ecosystems given current knowledge and management efforts. Our framework provides a platform to explore alternative management strategy outcomes and can help managers develop informed conservation plans with realistic expectations for return on investment of resources committed to sagebrush landscapes.

Arizona, California, Colorado, Idaho, Montana, Nev↗
USGS2026

Carbon sequestration along a gradient of tidal marsh degradation in response to sea level rise

Tidal marshes are considered one of the world's most efficient ecosystems for belowground organic carbon sequestration and hence climate mitigation. Marsh systems are however also vulnerable to degradation due to climate-induced sea level rise, whereby marsh vegetation conversion to open water often follows distinct spatial patterns: levees (i.e. marsh zones < 10 m from tidal creeks) show lower vulnerability of vegetation conversion to open water than basins (i.e. interior marsh zones > 30 m from creeks). Here, we use sediment cores to investigate spatial variations in organic carbon accumulation rates (OCAR) in a microtidal system (Blackwater marshes, Maryland, USA): (1) across a gradient of marsh zones with increasing marsh degradation, assessed as increasing ratio of unvegetated versus vegetated marsh area and (2) by comparing levees versus basins. We show that OCAR is up to four times higher on marsh levees than in adjacent basins. The data suggest that this is caused by spatial variation in three processes: sediment accretion rate, vegetation productivity, and sediment compaction, which are all higher on levees. Additionally, OCAR was observed to increase with increasing degree of marsh degradation in response to sea level rise. We hypothesize this may be due to more soil waterlogging in more degraded marsh zones, which may decrease carbon decomposition. Our results highlight that tidal marsh levees, in a microtidal system, are among the fastest soil organic carbon sequestration systems on Earth, and that both levees and basins sustain their carbon accumulation rate along gradients of increasing marsh degradation in response to sea level rise.

Maryland↗
USGS2026

Strategic approach for bird restoration evaluation - FY 2025 annual report

The Deepwater Horizon oil spill (hereafter “ spill ”) resulted in an estimated loss of 56,141–102,399 birds across 93 species along the U.S. Gulf of America coast. Injured species included those that nest locally along the northern Gulf of America coast as well as species that nest hundreds or thousands of kilometers outside of the region. The breadth of injury and the expansive geographic distribution of birds affected by the spill necessitated implementation of numerous and diverse projects designed to restore the injury of all affected bird species. As of September 2025, the Trustees, who are authorized to act on behalf of the public to restore injuries caused by the spill , have approved 39 projects, amounting to >60% of the total funds specifically allocated for bird restoration (i.e., Bird Restoration Type funding). Strategic allocation of remaining Bird Restoration Type funding is necessary to ensure that restoration of all injured bird species is achieved. An essential component of planning for future restoration is understanding the progress made from restoration actions conducted to date. The purpose of this project is to develop a plan that assesses the current restoration progress for each of the 93 injured species based on available data. The restoration assessment plan will be derived from a data gap analysis that identifies species for which data may be insufficient to assess restoration progress. This project will also allow Trustees to review and update the indicators of restoration success developed for injured bird species (e.g., population counts, reproduction metrics) that are used to evaluate restoration progress.

Annual Report↗
USGS2026

Restoration based on cost-benefit optimization: A grasslands pilot study

Ecological restoration is essential to meeting global biodiversity conservation goals. Given limited conservation budgets, deciding where to restore habitat is a key challenge for the coming decade. We developed a spatially explicit framework to optimize ecological restoration site selection by integrating land use history, species distributions, and economic costs. The framework includes the following steps: identify potential restoration area based on relevant environmental measures like land use; identify species of interest; calculate restoration benefits by modeling habitat and climate suitability and estimating reduced extinction risk associated with restoring a particular land parcel based on a modified species–area relationship; aggregate benefits across species; and compare to parcel-level land acquisition costs. We applied linear programming to maximize conservation benefit/restoration cost ratios to identify optimized restoration sites. We illustrate this approach using a case study for highly threatened grassland ecosystems in the Great Plains region of Kansas, USA. We selected five grassland animal species (greater prairie chickens [ Tympanuchus cupido ], lesser prairie chickens [ Tympanuchus pallidicinctus ], swift fox [ Vulpes velox ], pronghorn [ Antilocapra americana ], and regal fritillary [ Speyeria idalia ]) as indicators of restoration benefit across taxa. For the indicator species that we chose, shortgrass and mixed-grass prairies had the highest conservation benefit to cost ratio. Setting a minimum restoration threshold for each habitat type allowed us to identify high-priority tallgrass prairie sites. Despite increasing interest in ecological restoration, optimizing restoration site selection is challenging because one must consider habitat features that do not currently exist. The modeling approach described here is flexible and can be updated for different ecosystems, species, and conservation priorities. We outline potential alterations that can be made in future analyses, depending on desired restoration goals.

Kansas↗
USGS2026

Effects of carbamazepine to visual function in early life stage fish

The frequent detection of pharmaceuticals and personal care products (PPCPs) in the environment raises concern for aquatic systems. Carbamazepine (CBZ), an antiepileptic drug, is among the most detected PPCP globally, with concentrations in surface water exceeding those that induce toxicity to aquatic organisms. Non-targeted transcriptomic profiling was conducted in zebrafish ( Danio rerio ) larvae exposed to 0, 1, 5, 10, or 50 μg/L CBZ from 2 h post fertilization (hpf) through hatching, and then sampled at 48, 72, or 144 hpf. Transcriptomic profiles were annotated and characterized with in silico bioinformatic software to assess top enriched pathways and identify targets of environmentally relevant concentrations of CBZ and anchor molecular effects to higher levels of biological organization. Based on this analysis, CBZ was predicted to impair visual perception and sensory system development. The number of eye saccades, determined with a visually mediated behavioral assay, optokinetic response, was significantly reduced in 144 hpf larvae exposed to concentrations as low as 1 μg/L CBZ. These results indicate that environmentally relevant concentrations of CBZ may target and impact processes involved in visual function in fish.

Environmental Research↗
USGS2026

Net widening of Southern California beaches

Human impacts from dams reduce river sediment fluxes and are primary causes of coastal erosion worldwide. Here we provide new satellite-derived shoreline observation techniques to examine beach area trends across the diverse coastal settings of California. Contrary to global trends, these data reveal that the most heavily urbanized and dammed region of southern California experienced net beach growth of over 2 million m 2 during 1984-2024. While several beaches experienced severe erosion, overall widening is explained by sufficient sediment supply and concentrated widening from longshore transport captured at coastal structures and in littoral convergence zones. These results indicate that adequate sediment sources exist in this human-modified landscape to mitigate coastal erosion, but that this sediment is not effectively distributed to vulnerable beaches. This highlights the critical role that longshore sediment transport plays in long-term beach trends and illuminates management opportunities for coastal sustainability at the regional scale.

California↗
USGS2026

Magmatic source of the opening phase of Kīlauea’s 2018 Lower East Rift Zone eruption

The 2018 eruption of Kīlauea volcano in its Lower East Rift Zone began with the discharge of evolved high-Ti basalt as weak lava fountains and short, slow-moving lava flows. The lavas were quickly geochemically recognized as being derived from magmas stored within the rift zone and remobilized by a new intrusion, a sequence that is common at Kīlauea. This initial phase of the 2018 eruption, referred to as phase 1a, lasted for 6 days and was followed by extrusion of mixed magma after a 3-day pause. Even though remobilization of older rift zone magmas is common within Kīlauea’s rift zones, it is difficult to determine which past intrusion(s) may have initially emplaced those stored magmas. This difficulty stems from the tendency for Kīlauea magmas to follow very similar differentiation paths without significant variations in major, minor, or even trace element chemistry. We investigate possible magma sources for the lavas erupted during phase 1a of the 2018 eruption using whole-rock, mineral, and glass major and trace element compositions from historical East Rift Zone eruptions with adjacent fissures. We consider two primary hypotheses for the phase 1a source: magmas associated with the 1955 Lower East Rift Zone eruption or the nine eruptions in the Middle and Upper East Rift Zone during the 1960s. Our results suggest that magma associated with the earliest phases of Kīlauea’s 1955 eruption was the most likely source of the 2018 phase 1a remobilized magma. We determine volatile saturation pressures from melt inclusion chemistry and find similar storage depths for the 2018 phase 1a and early 1955 magmas. The phase 1a and early 1955 lavas are nearly indistinguishable in all of the compositional criteria considered, implying that the leftover 1955 magma body barely cooled and differentiated in the 63 years between eruptions (cooling rates of ~0.1 °C/year). This study sheds light on the potential for protracted storage of eruptible magmas in rift zones at Kīlauea, and highlights some of the challenges and solutions to identifying genetic relationships between magmas at Kīlauea.

Hawaii↗
USGS2026

Groundwater-level elevations in the bedrock aquifers of the Denver Basin aquifer system, Elbert County, Colorado, 2015–23

Water users in Elbert County, Colorado, rely on groundwater from bedrock aquifers in the Denver Basin aquifer system (upper Dawson, lower Dawson, Denver, Arapahoe, and Laramie-Fox Hills aquifers) for approximately half of their water uses. Withdrawals from the bedrocks aquifers have increased to meet the water use needs of expanding regional population growth and development. The U.S. Geological Survey, in cooperation with the Elbert County Board of County Commissioners, began a study in 2015 to monitor groundwater levels within Elbert County. The primary purpose of this report is to present a summary of groundwater levels measured during the study period (2015–23) and present results from statistical analyses of changes in groundwater-level elevations through time. Discrete groundwater levels were measured at 36 wells within Elbert County. Seven of those wells contained equipment to make and record continuous groundwater-level measurements at hourly intervals. All aquifers, except the lower Dawson aquifer, had only declining groundwater-level elevations in discrete measurements for wells with statistically significant trends. Of the eight statistically significant trends in the lower Dawson aquifer, two wells indicated increasing groundwater-level elevation from discrete measurements. The groundwater-level elevation trend medians in the upper Dawson, lower Dawson, Denver, Arapahoe, and Laramie-Fox Hills aquifers were −0.23, −0.66, −0.64, −0.39, and −0.63 feet per year, respectively, for discrete groundwater-level elevation measurements. Trends in continuous groundwater-level elevations were in agreement with statistically significant trends in discrete groundwater-level elevations for all wells. The groundwater-level elevation trend medians in this study, compared to the overall trends in a 2015−2018 study, both indicated declining groundwater-level elevations except in the upper Dawson aquifer, where the trend direction was opposite, a positive trend from 2015 to 2018 and a negative trend (declining groundwater elevations) from 2015 to 2023. The change in trends within the upper Dawson aquifer may be affected by differences in the study period and the trend analysis applied. Trends during the 2015–23 study period were compared to departures from the median 2015 groundwater-level elevation for each site in each aquifer. In general, the departures from the 2015 median supported trends observed at each site and correlated spatially with greater departures near the western border of Elbert County. Additionally, 30-year precipitation data showing wet and dry periods were overlaid with the departure from the 2015 median to assess groundwater-level patterns in wells in the five aquifers. Departures from the 2015 median groundwater-level elevations appeared greatest during the dry period between 2020 and 2023. Potentiometric-surface maps of the upper and lower Dawson aquifers created from static April 2023 groundwater elevations indicated groundwater-flow direction is generally from the south to the north. Results of this study could be used to guide additional groundwater monitoring in Elbert County and could aid in long-term planning of water resources.

Colorado↗
USGS2026

Differentiating persistent and intermittent euxinia from the molecular derivatives of green sulfur bacteria carotenoids

Green sulfur bacteria biomarkers that indicate euxinia within the photic zone sometimes co-occur with evidence of contradictory depositional redox conditions, such as oxygen-requiring fossils or bioturbation. Intermittent euxinia may explain this apparent contradiction, and recent studies of modern environments show that green sulfur bacteria dwell in transiently euxinic settings. As a result, new approaches are needed to distinguish green sulfur bacteria biomarkers indicative of persistent versus intermittent euxina in ancient sedimentary samples, which this study addresses by investigating how the distribution of isorenieratene derivatives relates to depositional redox conditions. To accomplish this objective, this study focuses on two drill cores through the Upper Cretaceous Mowry Shale and the Eagle Ford Group. These drill cores are comparably thermally immature according to vitrinite reflectance data, and these formations in these cores capture a spectrum of depositional redox conditions according to multiple organic and inorganic proxies, including newly acquired iron speciation and kerogen organic sulfur data. The results presented here reveal that higher molecular weight diagenetic products of isorenieratene are preferentially preserved under persistent euxinia compared to intermittently euxinic intervals that contain isorenieratene derivatives that are shifted to lower molecular weights. Further, the total inventory of aromatic carotenoid diagenetic products contains features that can be used to identify green sulfur bacteria biomarkers from reworked petrogenic sources. Accordingly, the diagenetic fate of isorenieratene and the distribution of its diagenetic products distinguish persistent versus intermittent euxinia, which can be used to sharpen our evaluation of euxinia in the geologic record.

Geochemica et Cosmochimica Acta↗
USGS2026

Millennial-scale climatic and cultural impacts on vegetation and fire at the southern edge of the Rocky Mountains, USA

Mountain forests and woodlands in semiarid regions of the world are threatened by climatic change and other human impacts. In the southwestern USA, climate and culturally driven changes to the structure and fire regimes of dry coniferous forests over recent centuries are well documented by tree-ring archives. However, the roles of climate and people as drivers of millennial-scale changes are less understood. We present a new record of vegetation dynamics and regional fire activity inferred from pollen, microscopic charcoal, and sediment geochemistry from Santa Fe Lake, NM (3532 m a.s.l.), the southernmost natural lake in the Rocky Mountains. Declining elemental intensities (e.g., Ti, Fe, Si) record local deglaciation after 14,000 cal yr BP followed by upslope expansion of trees indicated by increasing Pinus and Picea pollen percentages as the climate warmed during the Late Glacial. Cool, dry growing seasons probably limited the establishment of dense forests until after 12,200 cal yr BP when we hypothesize that Pinus ponderosa (ponderosa pine) expanded regionally (i.e., within 10–100 km) due to continued warming and a possible increase in monsoonal precipitation. Subalpine forests established near Santa Fe Lake by 10,300 cal yr BP as indicated by abundant conifer stomates and increasing Picea and Pinus aristata pollen percentages, which are highest after 5100 cal yr BP. Increasing Juniperus -type, Quercus , and Pinus edulis -type pollen record changes in vegetation belts at lower elevations after 10,300 cal yr BP. Maize pollen, an unambiguous indicator of Puebloan agriculture, first appears at 800 cal yr BP (1150 cal yr CE) and is present through 1600 cal yr CE when Spanish colonization brought cultural upheaval and population loss to local indigenous communities. Regional fire activity, inferred from microscopic charcoal influx, remained relatively constant for most of our record. However, average microscopic charcoal influx more than doubled relative to the rest of the Holocene during the last 2000 years, with the highest influx after 1550 cal yr BP (400 cal yr CE). The increased influx of microscopic charcoal coincides with archaeological evidence for dramatic regional population growth in the northern Rio Grande region and probably reflects extensive cultural burning associated with the development of an agricultural economy. Although climate was an important driver of vegetation change and fire regimes in the mountains of the southwestern USA during the Holocene, our results suggest that human impacts on regional burning during the last 2000 years far exceeded earlier climate-driven variation.

Rocky Mountains↗
USGS2026

Thinking outside the rocks: Subsurface water storage, topography, and land cover are key modulators of large-scale riverine dissolved silicon dynamics

Riverine dissolved silicon (DSi) dynamics reflect integrated geologic, hydrologic, climatic, and ecological controls. We compiled annual DSi data for 337 rivers across four continents and trained interpretable machine-learning models to predict concentrations and yields from 28 watershed variables. Both models reproduced testing data ( R 2 = 0.85 for concentration and 0.96 for yield) and withheld-site validation ( R 2 = 0.91 and 0.93). Lithology, especially volcanic rock fraction, strongly controlled DSi while subsurface storage, topography, and land cover further shaped DSi dynamics. DSi concentrations and yields exhibited nonlinear responses to basin slope, recession-curve slope, proportion of open-water cover, and nutrient availability. Concentrations showed sharper threshold responses to hydrologic and biotic variables, whereas yields varied more gradually with climate and lithology. These results provide a framework for forecasting DSi under land cover and climate change and for embedding realistic, nonlinear processes in mechanistic models.

Geophysical Research Letters↗