USGS Publications Warehouse

Hydrogeology, groundwater salinity distributions, and assessment of the effect of oil-production activities on groundwater in the Midway Valley area, western Kern County, San Joaquin Valley, California

Fri, 24 Apr 2026 14:10:27

This study seeks to determine the effects of oil field produced water disposal operations and well mechanical integrity issues on groundwater quality in oil fields in the southwest San Joaquin Valley, California. Whereas previous studies used groundwater wells to study shallow aquifers outside the oil fields, this study demonstrates that future approaches may use oil well geophysical logs to map groundwater head gradients, create salinity profiles and document changes in salinity over time in oil field areas with sparse groundwater well data and at depths greater than 330 m. We also incorporate an analysis of well histories to determine potential effects of compromised wellbore seals on changes in aquifer quality that cannot be explained by water disposal practices. Water quality in the aquifers is naturally brackish across most of the area, with better quality groundwater occurring in the eastern part. Geophysical logs are used to determine salinity variations within aquifers including the depth at which TDS exceeds 10,000 mg/L. This depth ranges from 366 m in the northwest to approximately 1,500 m in the southeast. Oil well porosity logs are used to determine water table elevations. These logs indicate the water table slopes south-southeast, showing the predominant groundwater flow direction is from oil field disposal areas toward better quality groundwater east of the oil fields. Geophysical logs show formation resistivity near some disposal facilities has decreased over time, indicating the salinity of the aquifer has increased due to disposal of saline produced water in injection wells and ponds. Oil well history analysis suggests that increased salinity over time in water-saturated sand intervals >1.5 km from disposal facilities may be caused by mechanical failures and/or incomplete borehole seals in poorly constructed or abandoned wellbores prevalent throughout the study area—particularly wells drilled prior to 1930.

Critical minerals in zinc ore—An update on Earth Mapping Resources Initiative Research in the Boulder Batholith region, Montana

Fri, 17 Apr 2026 14:29:21

Computation of regional groundwater budgets for the Virginia Coastal Plain aquifer system

Pope, Jason; Gordon, Alison; Frederiks, Ryan — Mon, 20 Apr 2026 17:44:26

Computation of detailed groundwater flow budgets for subdivisions of the Virginia Coastal Plain aquifer system has enabled quantification and more thorough understanding of groundwater flow within this important water resource. A zone budget analysis based on previously published groundwater models of the Virginia Coastal Plain and Virginia Eastern Shore indicates that groundwater conditions vary substantially throughout the Coastal Plain aquifer system because of local variations in hydrogeology and historical and ongoing variations in groundwater use and management. Decades of substantial groundwater withdrawal from the Coastal Plain aquifer system have altered groundwater flow from predevelopment conditions. Rates of sustainable withdrawal are limited because the downward groundwater flow rate into confined aquifers is a relatively small part of the total groundwater budget for the aquifer system compared to the rate of recharge at the land surface.

Analyses of groundwater budgets from the Virginia Coastal Plain model indicate that groundwater flow is generally outward from the surficial aquifer to rivers and coastal waterbodies and downward through a series of underlying aquifers and confining units to the Potomac aquifer, which is the deepest aquifer and the source of most groundwater withdrawals. Downward flow into the Potomac aquifer is estimated to be only 7 percent of total net precipitation-derived net recharge at the land surface but makes up about 66 percent of inflow to the aquifer in Virginia, with much of the remaining inflow occurring laterally from outside of defined groundwater budget regions in Virginia. For several decades prior to 2010, high rates of withdrawal from the Potomac aquifer resulted in substantial decline in groundwater storage in the aquifer and in most overlying aquifers and confining units. From 2010 to 2023, rates of withdrawal substantially lower than the historical maximum resulted in small net increases in groundwater storage in the confined aquifer system for most regions of the Virginia Coastal Plain. Nevertheless, for the same period, groundwater storage for the entire model domain continues to incrementally decline, indicating that storage recovery in Virginia is offset by a continued decrease in storage in areas beneath the Chesapeake Bay or adjacent areas of Maryland and North Carolina. Withdrawals from the Potomac aquifer have induced substantial downward flow which is a large part of groundwater budgets for confined aquifers such as the Potomac. For the most recent simulated conditions (2023) downward groundwater flow continues, but because vertical flow rates are a function of the difference between water pressure in the upper surficial systems and lower confined units, rates of downward flow are lower than those in earlier decades as the confined water levels partially recover from larger groundwater withdrawals in the past. Geographically, groundwater flow is generally inward from perimeter regions of the Virginia Coastal Plain toward central regions with the largest withdrawal rates. Groundwater inflow from coastal regions could be contributing to saltwater intrusion, even though that was not measured in this study.

Analyses of groundwater budgets from the Virginia Eastern Shore peninsula, a geographic region of the Virginia Coastal Plain, indicate that groundwater flow for that isolated aquifer system is generally outward from the surficial aquifer to coastal water bodies and downward into the confined Yorktown-Eastover aquifer system, which is the source of most withdrawals. Downward groundwater flow into the confined Yorktown-Eastover aquifer system is estimated to be less than 2 percent of total recharge and less than 9 percent of net recharge at the water table but makes up more than 93 percent of all inflow to the confined aquifer system. Decades of substantial but relatively consistent groundwater withdrawals have induced greater downward flow rates into the confined aquifer system but also have resulted in loss of groundwater from storage. For the most recent simulated period (2023), estimated storage loss accounts for slightly under 7 percent of withdrawals from the confined aquifer system. The reported withdrawal rate for this period from the confined Yorktown-Eastover system is near the highest reported rate for the Virginia Eastern Shore, which means that the storage depletion is expected to continue, even though groundwater levels appear to be relatively stable. Estimated groundwater flow rates upward from the confining unit underlying the Yorktown-Eastover system and low rates of inflow from coastal water bodies underscore ongoing concerns about up-coning and lateral intrusion of salty groundwater.

Analyses of meteorological and hydrological records support Tribal members’ accounts of changing climate on the Fort Apache Reservation, east–central Arizona

Mason, Jon — Mon, 20 Apr 2026 17:40:45

The Fort Apache Reservation in east–central Arizona, home to the White Mountain Apache Tribe of the Fort Apache Reservation, Arizona, contains several climate zones because of the large variation in surface elevation within the reservation. This study was carried out in cooperation with the White Mountain Apache Tribe of the Fort Apache Reservation, Arizona, to raise awareness of how the changing climate affects the Fort Apache Reservation. This report documents the evaluation of existing multidecadal meteorological and hydrological datasets for the Fort Apache Reservation, used to evaluate the effects of a changing climate on the reservation. In this evaluation, near-surface air temperature, snow depth, snow water equivalent, precipitation, and streamflow datasets were analyzed for monotonic trends indicative of changing climatic conditions during specified periods of time. The results of these trend analyses were then compared with the Tribal community's memories of the changing climate.

Trend analysis of near-surface air temperatures from a U.S. Historical Climatological Network station on the Fort Apache Reservation at Whiteriver, Arizona, indicated that mean annual air temperatures have increased by an average of 2.48 degrees Fahrenheit from 1980 to 2023. Records from the same station also indicated that average monthly maximum temperatures recorded for March increased by 5.39 degrees Fahrenheit for the same time period.

Annual precipitation at the five precipitation stations used in this study decreased greatly from the 1980s to 2023. The largest total decrease was 10.07 inches, or 34.7 percent. However, only one of the two precipitation stations with longer term data available prior to 1980 had a significant negative trend when data from the entire period of record, from 1901 to 2023, were analyzed.

Trend analyses show a decrease in the annual maximum snow water equivalent and an earlier disappearance of the snowpack at two Natural Resources Conservation Service snow telemetry stations in the mountainous region just east of the Fort Apache Reservation from 1981 to 2023. Based on the trend analyses, the average annual maximum snow water equivalent has decreased by more than 40 percent at both stations, and the average date when the snowpack was fully melted at the stations in the spring has moved earlier in time from late April to early April or late March. However, a statistically significant trend was not determined for the early April snow water equivalent measured at a nearby Natural Resources Conservation Service snow course across its period of record, indicating that the history of mountain snowpack in this area is not fully understood. Analysis of snowfall data from a National Oceanic and Atmospheric Administration Cooperative Observer Program network station on the Fort Apache Reservation at McNary 2N, AZ (station 025412) indicated that, on average, the measured total annual snowfall at the station decreased 42.4 percent from 1935 to 2023.

Streamflow data from six U.S. Geological Survey streamgages on the Fort Apache Reservation were analyzed for trends. For most streamflow gages, statistically significant trends were not determined for tested parameters when the entire streamflow period of record was used for stations with records going back to at least the 1960s. However, when the data from 1980 to 2023 was tested, most of the streamflow parameters had statistically significant negative trends. All six streamgages showed a decrease in average annual runoff of at least 50 percent from 1980 to 2023; one streamgage showed an 81.8 percent decrease.

A similar statistical finding was observed in the analysis of the annual spring snowmelt peak from one of the six streamgages used in the study and located in an area receiving measurable amounts of snowmelt runoff. When data from the entire period of record (1958–2023) was used, no trend in streamflow was determined; however, a significant negative trend was determined from 1980 to 2023, indicating a decrease in average annual springtime runoff of 62.6 percent. Statistical analysis on the timing of the annual spring snowmelt peak at the same streamgage indicated the snowmelt peak is happening on average about 12 days earlier now (2023) than it did in the past. The trend results for the timing of the annual spring snowmelt peak were the same and statistically significant for both periods tested (1958–2023 and 1980–2023). Two of the streamflow records from the Fort Apache Reservation were compared to the Palmer Hydrological Drought Index computed for Arizona Climate Division 4 (East Central) by the National Centers for Environmental Information. The comparison showed that the streamflow records generally tracked the Palmer Hydrological Drought Index.

In interviews, Tribal community members living on the Fort Apache Reservation described the changes in climate that they observed during their lifetimes. Common themes reported were that air temperatures have become warmer, and the weather is less predictable with changes in seasonal patterns. Drier conditions, lower snowfall, shorter winters, and lower river levels were also reported. These community member observations align with the results of this study.

Arsenic and isotope concentrations in the lower Platte River valley of eastern Nebraska, early 1970s to 2023

Moser, Matthew; Cherry, Mikaela; Hall, Brent — Thu, 16 Apr 2026 17:22:01

The City of Lincoln, Nebraska, has been monitoring concentrations of arsenic in their source water and evaluating their options for treatment and removal since at least 2002. In 2022, the City of Lincoln, Nebr., with funding assistance from the Nebraska Water Sustainability Fund, began cooperating with the U.S. Geological Survey to examine arsenic concentrations in surface water and groundwater in the lower Platte River valley and the area around City of Lincoln Water System (LWS) well field. Arsenic data collected from the Platte River since 1974 were examined using the “weighted regression on time, discharge, and season” model, which compared the streamflow (also referred to as “discharge”), time of year, and season to estimate concentrations of arsenic. Annual mean arsenic concentrations modeled for more than 49 years at the Platte River at Louisville, Nebr., U.S. Geological Survey streamgage (station 06805500), indicated a significant increasing trend. Arsenic concentrations in the Platte River were seasonal, with the highest concentrations being observed during mid- to late summer. When seasonal patterns and streamflow were combined with arsenic concentrations in the Platte River during low streamflow conditions, groundwater contributions, which can have higher arsenic concentrations, make up a larger portion of the streamflow. Arsenic samples were collected from upstream rivers in 2022 and 2023 and were paired to analyze the arsenic contributions at the U.S. Geological Survey streamgage on the Platte River near Ashland, Nebr. (station 06801000), near the City of Lincoln well field. The arsenic concentrations from the streamgage on the Platte River near Ashland, Nebr., location, were higher than the U.S. Geological Survey streamgage on the Elkhorn River at Waterloo, Nebr. (station 06800500), and significantly lower than at the U.S. Geological Survey streamgage on the Platte River near Leshara, Nebr.(station 06796500), indicating that the Platte River usually contributes a higher concentration of arsenic than does the Elkhorn River as they join near Ashland, Nebr. During 1991–2023, six groundwater monitoring wells were analyzed to identify trends in arsenic concentrations. Two of the six wells had a positive trend during the 33-year period. One monitoring well did not reveal a long-term trend during this period but showed a trend during 2019–23, correlating to a period when the island in the middle of the Platte River was connected to the east bank of the river when manganese reducing conditions were present and groundwater levels were declining in the well. Across all wells the oxidation and reduction (redox) condition during the time of sampling was assessed. Mixed anoxic and (or) oxic redox condition was the most common redox process and the highest sampled arsenic concentrations in monitoring wells were observed in anoxic conditions driven by manganese reduction. Groundwater arsenic concentrations had seasonal variation around the City of Lincoln well field, with higher arsenic concentrations tending to be further south in comparison to samples collected further north. Isotope samples were collected and analyzed in surface water and groundwater around the LWS well field. The samples indicate that the proportion of surface water present in the LWS production wells can be higher in the spring and lower in the summer. With higher arsenic concentrations observed in the stream water during the summer period, the LWS source water can be affected by these elevated arsenic concentrations even though the proportion of surface water is lower.

Incorporating data sets with multiple sources of uncertainty in integrated species distribution models

Lunt, Fiona; Scher, C.; Mummah, Riley; Miller, David — Tue, 14 Apr 2026 16:29:45

Data integration methods aim to improve species distribution estimates by incorporating multiple sources of uncertainty across datasets. Two major sources of uncertainty are: (1) variation in sampling effort across space and within datasets, and (2) variation in reliability associated with data collection protocols or timing among datasets. Our goal was to evaluate how different approaches to address these uncertainties influence predictive performance of integrated models. We modeled distributions of four bird species using three datasets that differed in sampling design. We examined three strategies to reduce uncertainty: (1) filtering data, (2) incorporating functions that account for uncertainty in observation models, and (3) varying how datasets are integrated into a single estimate. We first examine methods to account for variable effort in observations, focusing on both spatial differences in sampling intensity and effort given to a single observation record. We then examine approaches to account for data sets with differing reliability. Sampling effort was best addressed through conservative filtering, including spatial thinning and excluding observations with highly variable effort. Next, we considered how to account for potential false positive detections—due to either misidentification or changes in distributions. We found that treating less reliable data as a covariate, an approach previously suggested for data integration that can greatly speed up model fitting, performed well. Other effective approaches included directly modeling false positive rates and complete exclusion of less reliable data sets. Our results provide insights into best practices in integrated modeling for handling uncertainty in integrated models. We demonstrate the flexible options available when using integrated models to address uncertainty.

Mineral chemistry perspective on remobilization of stored magma at Kamakai'a Hills, Southwest Rift Zone of Kilauea, Island of Hawai'i, USA

Downs, Drew; Sas, May — Mon, 13 Apr 2026 15:00:24

Differentiated magmas stored in the rift zones of Kīlauea have received more attention in recent years following eruption of andesite during the early phase of 2018 lower East Rift Zone activity. Despite this growing interest, some of the most voluminous eruptions of differentiated rift zone magmas remain poorly studied. One such eruption, and the most voluminous exposed differentiated flow field at Kīlauea, is the Kamakaiʻa Hills. This eruption took place in the Southwest Rift Zone of Kīlauea, a region that is hypothesized to contain a long-lived rift zone reservoir. The Kamakaiʻa Hills flow field encompasses >250 × 10⁶ m³ of basaltic andesite and basalt compositions with a mineral assemblage of orthopyroxene + clinopyroxene + plagioclase during its early ʻaʻā phase and clinopyroxene + plagioclase + olivine during its late pāhoehoe phase. To better understand storage conditions and magma accumulation, this study focuses on major, minor, and trace elements from the mineral assemblage present within the early ʻaʻā and late pāhoehoe phases. The diversity of clinopyroxene and plagioclase compositions within the early ʻaʻā and late pāhoehoe phases, as well as diverse compositions of plagioclase and orthopyroxene within the early ʻaʻā phase, suggest multiple magma bodies and limited pre-eruption magma mixing within the broader Kamakaiʻa Hills reservoir. Oscillatory zoning patterns (particularly in clinopyroxene) imply processes such as recharge events, magma mixing or mingling, or convection within a differentially cooling, chemically stratified reservoir over protracted time intervals, whereas only limited resorbed mineral textures indicate incomplete mixing of heat and chemically distinct magmas during the dike intrusion that triggered the eruption. Mineral-mineral and mineral-melt thermobarometry indicate predominantly shallow (≤2.5 km depth) crustal storage conditions of the cooled, differentiated magma (∼1100 °C and cooler for the basaltic andesites) to hotter temperatures for the basalts (all >1100 °C). Despite the known large standard errors estimated for mineral-melt and mineral-mineral barometry (10s to >100 MPa), the calculated pressures and depths broadly correspond with earthquake swarm depths beneath the Kamakaiʻa Hills, and drill core and fluid inclusion barometry storage depths of differentiated magmas within the lower East Rift Zone. The Kamakaiʻa Hills differentiated magmas have H₂O contents (∼0.5 wt%, using plagioclase-melt hygrometry) equivalent to typical Kīlauea basalts. Our data and interpretations demonstrate a complex, long-lived rift zone storage system that consisted of multiple magma bodies and was mobilized into eruption through intrusion of a hotter and more primitive summit-derived (uprift) magma.

Deep groundwater total dissolved solids mapping in the Dakota Group, Williston Basin, USA

Stephens, Michael; Hoogenboom, Bennett; Ball, Lyndsay; Chang, Will — Tue, 7 Apr 2026 14:23:14

Growing concern about the quantity of available freshwater around the world has led to interest in surveying groundwater total dissolved solids (TDS) below water well depths. Deep TDS has not been systematically mapped, and there is much to learn about the distribution and controls on deeper groundwater. In sedimentary basins across the United States, groundwater resources often overlie hydrocarbon resources, providing an opportunity to use borehole geophysical data collected for hydrocarbons to characterize groundwater and pore space resources. This study adapts a recently developed subsurface geostatistical and geophysical modeling approach to continuously map groundwater TDS, porosity, and temperature in the Dakota Group of the Williston Basin—an undercharacterized regional aquifer system overlying deeper hydrocarbon reservoirs. Groundwater TDS in the Dakota Group ranges from approximately 4800 to 26,900 mg/L. TDS patterns are stratified with higher TDS in the lower and upper Dakota Group, and relatively lower TDS in the middle Dakota Group. The lower TDS in the middle zone may represent a preferential regional flow path for lower-TDS meteoric recharge from the west. The alternating pattern of TDS may also be evidence of higher-TDS inflows into the Dakota Group from underlying and potentially from overlying aquifers. Porosity is lower near the center of the Williston Basin and tends to be higher to the east, which may be related to grain size distributions. The new regional TDS and porosity modeling serves as a quantitative reference for water users and provides supporting evidence for hypotheses on Dakota Group recharge.

Seabed maps showing topography, ruggedness, backscatter intensity, sediment mobility, and the distribution of geologic substrates in quadrangle 3 of the Stellwagen Bank National Marine Sanctuary region offshore of Boston, Massachusetts

Valentine, Page; Cross, VeeAnn — Fri, 3 Apr 2026 17:31:25

The U.S. Geological Survey, in cooperation with the National Marine Sanctuary Program of the National Oceanic and Atmospheric Administration, has conducted seabed mapping and related research in the Stellwagen Bank National Marine Sanctuary (SBNMS) region since 1993. The area being mapped using geophysical and geological data includes the SBNMS and the surrounding region, which totals approximately 3,700 square kilometers (km²) and is subdivided into 18 quadrangles. The seabed is a glaciated terrain that is topographically and texturally diverse. Quadrangle 3, the subject of this scientific investigations map, has a mapped area of 185 km² and has water depths that range from about 30 meters (m) on the Stellwagen Bank crest to about 135 m in a basin east of South Ninety Bank, which lies off the eastern margin of Stellwagen Bank. Seven map types, each at a scale of 1:25,000, depict seabed topography, ruggedness, backscatter intensity, distribution of geologic substrates, sediment mobility, distribution of fine- and coarse-grained sand, and substrate mud content. These maps show the distribution of geologic substrates on the southeastern part of Stellwagen Bank, on adjacent banks and basins in deeper water to the east, in the eastern part of Race Point Channel to the south of the bank, and on the northern slope of Cape Cod. Interpretations of multibeam sonar bathymetric and seabed backscatter imagery, photographs, video imagery, and grain-size analyses were used to create the geology-based maps. Data from 309 stations were analyzed, including 279 sediment samples. The geologic substrate maps of quadrangle 3 show the distribution of 21 geologic substrates that represent a wide range of textures, such as rippled sand, immobile sand, immobile muddy sand, sand that partially veneers gravel, and boulder ridges. Mapped substrates are characterized by sediment grain-size composition, surface morphology, substrate layering, the mobility or immobility of substrate surfaces, and water depth range. This scientific investigations map portrays the major geological elements (substrates, topographic features, and processes) of environments in quadrangle 3. It is intended to provide a foundation for research into present and past sediment transport processes in a complex terrain, provide insights into the ecological requirements of invertebrate and vertebrate species that use the various substrates, and support seabed management in the region.

Geologic map of the Emmons Lake volcanic center, Alaska

Thu, 2 Apr 2026 14:28:05

Introduction

The Emmons Lake volcanic center is a spatially clustered group of stratovolcanoes and calderas in the southwestern part of the Alaska Peninsula, Alaska. The volcanic center is characterized by several ice- and snow-clad stratovolcanoes located within and along the margins of a nested-caldera complex that includes Emmons Lake. A shieldlike ancestral edifice (ancestral Mount Emmons) is truncated by the caldera complex and forms a broad volcanic platform around the center. The main stratovolcanoes of the Emmons Lake volcanic center are Pavlof Sister, Pavlof Volcano, Little Pavlof, Double Crater, Mount Hague, and Mount Emmons. Several small unnamed cinder cones and vents also are located within Emmons Lake volcanic center and on the east flank of Pavlof Volcano. Many of these cones and vents have been the source of the young lava flows that mantle the floor of the caldera. Pavlof Volcano, in the northeastern part of the Emmons Lake volcanic center, is one of the most historically (that is, the past about 300 years) active volcanoes in Alaska, and eruptions from Pavlof Volcano pose the greatest hazards to the region.

Volcanic rocks of the Emmons Lake volcanic center overlie continental and marine sedimentary rocks of chiefly Late Jurassic to early Tertiary age. The oldest rocks in the area are those of the Naknek Formation, consisting of volcaniclastic sandstone, siltstone, and conglomerate of Late Jurassic age. The southern part of the area includes rocks of the Belkofski Formation, a thick sequence of volcaniclastic sandstone, siltstone, and conglomerate of middle Tertiary age. Lava flows, volcanic breccia, and fluvial volcaniclastic rocks of late Miocene age, which unconformably overlie the Belkofski Formation south of the Emmons Lake volcanic center, are primarily exposed on the islands just south of the Alaska Peninsula.

The Emmons Lake volcanic center was affected multiple times by glaciation associated with the glacier expansion that characterized the Quaternary. Glaciation has played a key role in shaping the present-day landscape, and much of the eruptive history of the Emmons Lake volcanic center has involved interactions with glacier ice. Thus, a brief review of the Quaternary glacial history of the area is provided to establish the physical context for Emmons Lake volcanic center eruptive activity.

Low streamflows in Massachusetts: Variability over space and time and relations with climatic and basin variables

Chamberlin, Catherine; Hodgkins, Glenn — Mon, 6 Apr 2026 14:19:40

Streamflows in Massachusetts have set record lows in recent years despite generally wetter conditions than during the drought of the 1960s, and the reasons for this are not known. To analyse potential drivers of low streamflows in Massachusetts, six low-flow metrics were computed at 107 streamgages. These metrics represent low-flow magnitude, magnitude normalized to median flows, and duration. Multiple linear regressions were used to analyse the variability of low flows over space and time. Potential explanatory variables were computed using climatic, land use, water use, and basin data. For all low-flow metrics, the ratio of precipitation to potential evapotranspiration (P/PET) in July–August explained the most variability, with decreasing P/PET largely explained by lower precipitation. Water/wetland area was a significant explanatory variable in all the normalized-magnitude and duration models, with greater area associated with lower normalized magnitudes and with shorter durations of low flows. Human influence (characterized by development, population, water use, and artificial water storage) had mixed effects. Trends from 1983 to 2022 in summer P/PET and human influence have been strongest in the eastern part of the state where the strongest decreases in flows are observed. Low flows in Massachusetts seem to be driven by a combination of low summer precipitation and human effects, though the specific mechanisms of human influence on flow likely vary between basins.

Early Miocene volcanic rocks and associated tectonics, Lava Hills and southern Bristol Mountains, California

Miller, David M.; Harvey, Janet; Buesch, David; Gans, Phillip B. — Thu, 23 Apr 2026 14:17:09

Volcanic rocks of latest Oligocene to early Miocene age form an east-west belt across part of the central eastern Mojave Desert from the Whipple Mountains on the east to the Rosamond Hills on the west. We term this the central belt because it is separated from northern and southern belts by swaths with no volcanic rocks. Limited geochronologic data indicate that much of the belt is latest Oligocene and early Miocene in age, about 24 to 19 Ma, a finding that is consistent with these rocks being overlain by the 18.8 Ma Peach Spring Tuff in many places.

We describe Miocene geology in a central area of the belt, in the Lava Hills, southern Bristol Mountains, and southern Old Dad Mountains. Sedimentary basins formed coeval with early andesite to rhyolite volcanism, progressing from fluvial and lacustrine tuffaceous sandstone to volcanic lavas, tuffs, and breccias, indicating that early basins formed proximal to volcanic edifices. Higher strata are fluvial and lacustrine with lavas punctuating the sequence. Although basins may partly have been formed within topographic lows bounded by volcanic domes, plateaus, and stratovolcanoes, consistent stratigraphic sections over wide areas indicate that tectonic basin evolution affected broad areas. The volcanic section is capped by local basalt flows and the regional Peach Spring Tuff. Limited data on normal faults support interpretations of early extensional basin development caused by northeast-southwest oriented stretching. Later extension caused stratal rotations, tilting early deposits down to the southwest. This tilted and subsequently beveled basin architecture was overlain by the youngest volcanic deposits, primarily rhyolite and basalt. The Peach Spring Tuff, 18.8 Ma, lies within this upper unit. Similar stratigraphic and structural relations are exposed in the nearby Marble Mountains and Van Winkle Mountain sections, reinforcing that a broad area underwent similar volcanism and tectonism. In our study area the upper unit is only very gently tilted except near dextral strike-slip faults of the eastern California shear zone. These late Miocene to Recent faults are represented as four main fault zones spaced about 5 km apart, representing distributed shear north of the Bristol Lake basin.

Small cumulative survival costs of enzootic disease could suppress long-term population size

Mon, 20 Apr 2026 14:05:28

Fungal pathogens can cause epizootics that result in widespread mortality and rapid population declines in some species. However, even in the absence of high disease-induced mortality, enzootic mycoses could have large-scale impacts on host population dynamics. Here, we examined the effects of ophidiomycosis, an enzootic fungal disease, on a Louisiana snake community over a 3-year period using a multi-state Jolly–Seber model with disease-state misclassification. We did not detect a difference between the average weekly apparent survival probability of uninfected and infected hosts for either Nerodia species or Thamnophis proximus. We also found a strong positive association between snout-to-vent length and weekly apparent survival probability across all species. We found that recruitment of infected hosts was slightly higher than recruitment of uninfected hosts for two of the three species. Population projections suggested divergent trajectories between disease-present and disease-absent scenarios, where disease-absent populations had higher abundance than disease-present populations. Our results highlight that small differences in survival can accumulate over time, as well as the challenges of quantifying population-level impacts of enzootic diseases when survival differences are not readily detected, underscoring the importance of continued long-term monitoring to assess whether ophidiomycosis affects snake population dynamics.

Sources and streambed storage of soft sediment and sediment-bound phosphorus in an agricultural Great Lakes tributary

Thu, 23 Apr 2026 14:47:18

The East River, an agricultural tributary to the Lower Fox River and Lake Michigan in Wisconsin, USA, has excessive phosphorus (P) and suspended-sediment loads that contribute to downstream eutrophication and habitat-related impairments. Spatial variations and connectivity in the sources and streambed storage of soft, fine-grained (silt and clay) sediment and related sediment-bound P (sed-P) were examined, from first-order ephemeral channels to a downstream water-monitoring station. Analysis included field inventories, a channel corridor sediment and sed-P budget applied to an extended channel network, and geochemical fingerprinting. Corridor inventories included mass wasting along valley sides, eroding streambanks, gullying along perennial and ephemeral channels, and streambed storage volumes in perennial reaches; each converted to masses. Erosion results estimate 7400 Mg/yr of fine-grained sediment, similar to the mean annual suspended sediment load of 5400 Mg/yr. Corridor erosion contributed 7200 kg/yr of sed-P, less than the mean annual particulate-P load of 10,000 kg/yr P. Soft sediment storage was 1400 Mg, with 1500 kg sed-P. Apportionment of soft sediment as streambank sourced was spatially variable, contributing ≥95 % in high order reaches with high storage and as little as 20 % in upstream reaches, where gully, crop, and forest provided the remainder. Two nearby tributaries showed similarity in the predominance of streambank-sourced material in stored soft sediment but differences in geomorphic setting affected its spatial distribution. The results of this study show the importance of including corridor erosion as a source of sediment and sed-P in agricultural basins, which can be helpful in decision-making regarding conservation practices.

Moving toward a more human-oriented analysis of urban heat: Examining differences of heat exposure intensity at busy commuting locations

Tue, 24 Mar 2026 16:52:30

Examining urban thermal environments has become a critical area of research spanning epidemiology, urban planning, and ecology. While traditional metrics like air temperature (T_air) and satellite-derived surface temperature dominate urban heat studies, these measures often fail to reflect how people actually experience thermal exposure intensity. More human-oriented metrics, such as mean radiant temperature (MRT), and the wet bulb globe temperature (WBGT), better capture this lived experience, particularly at locations where people are likely to encounter outdoor heat, such as bus stops. Human demographics further complicate heat exposure, as access to cooling resources like trees and greenspaces can vary by neighborhood income. Our study addresses these complications by collecting thermal data across 60 commuting locations in Denver, Colorado in the summer. We evaluate (1) the extent to which more human-oriented metrics capture thermal exposure compared to T_air and LST, and (2) how heat exposure varies by neighborhood income levels. We observed that bus stops in low-income neighborhoods had an MRT increase 2.8 °C compared wealthier neighborhoods, and that income-driven differences in MRT and WBGT strongly depended on bus stop aspect. After accounting for solar orientation, differences in MRT increased to as much as 6.3 °C at north-facing stops. Our results suggest tree canopy shade explains some observed heat exposure patterns, with south facing bus stops seeing a MRT and WBGT decrease of 0.42 °C and 0.11 °C from a percent increase in tree canopy. Interestingly, depending on bus stop aspect, nearby buildings can increase MRT and WBGT (facing east), or decrease MRT and WBGT (facing south) If planners aim to address this issue, consideration of bus stops, and land covers configuration may help.

Spatial and temporal geochemical variations of lava flows and tephra deposits from the December 2020 to September 2024 eruptions of Kīlauea volcano

Downs, Drew; Lynn, Kendra; Winslow, Heather; Lundblad, Steven; Decker, Meghann — Tue, 24 Mar 2026 15:58:48

Kīlauea volcano underwent dramatic morphological changes in 2018. That year recorded the end of the 35-year-long eruption of Puʻuʻōʻō (1983–2018) and 10-year-long (2008–2018) Halemaʻumaʻu lava lake and emplacement of the ~4-month-long lower East Rift Zone lava flows that coincided with ~500 m of summit caldera collapse. Starting on December 20, 2020, eruptions resumed at Kīlauea’s summit. There were five summit eruptions between December 2020 and September 2023, which ranged in duration from more than a year to as short as a week. Following these summit eruptions, seismicity and deformation increased in the upper Southwest Rift Zone in 2024, culminating in a ~8.5-h-long eruption in this region on June 3, 2024. Increased seismicity and deformation then shifted to the upper and middle East Rift Zone and after several months culminated in an eruption just west of, and within, Nāpau Crater in the middle East Rift Zone from September 15 to 20, 2024. Despite vast morphological changes at Kīlauea’s summit, the geochemical compositions (i.e., whole rock and glass) that erupted from December 2020 to September 2023 are all remarkably similar to each other. Whole-rock compositions appear distinct from the preceding 2008–2018 Halemaʻumaʻu lava lake and phase 3 (i.e., summit or uprift-derived mafic lavas) of the 2018 lower East Rift Zone lava flows, although glass compositions appear to have more overlap with 2018 lower East Rift Zone glasses. The June 3, 2024, upper Southwest Rift Zone spatter and lava flows exhibit a dramatic enrichment in whole-rock MgO that is not recorded in glass, which reflects accumulation of olivine (e.g., antecrysts or xenocrysts) during dike emplacement, and is consistent with the abundance of olivine in the lava flows (5–10%). June 2024 Southwest Rift Zone whole-rock and glass compositions overlap with those erupted at the summit from December 2020 to September 2023, whereas some whole-rock trace (i.e., Sc, Sr, and Zr) and major elements (i.e., CaO) are suggestive of mixing with a magmatic component that had fractionated plagioclase and pyroxene and/or a new parental magma influencing the summit reservoir system. The September 15–20, 2024, eruption at Nāpau Crater in the middle East Rift Zone involved the most differentiated magma since eruptive activity resumed in December 2020, with its magma fractionating olivine + plagioclase + pyroxene. The September 15–20, 2024, composition resembles Puʻuʻōʻō lava flows that erupted in, or near, Nāpau Crater in 1983 (episode 1), 1997 (episode 54), and 2011 (episode 59), with episode 59 having a compositional cluster that is most similar to that of the September 2024 lava flows. The data presented and provided herein open new research perspectives for long-term analyses of geochemical variations following caldera collapse at Kīlauea volcano and facilitate comparisons with other basaltic caldera systems worldwide.

Finding the (small) cores: Spatial covariance tracks grassland bird community occupancy in fragmented grasslands

Berry, Lauren; DeGregorio, Brett; Uden, Daniel; Roberts, Caleb — Thu, 2 Apr 2026 18:04:33

Grasslands are an imperiled ecosystem, and grassland bird abundance is declining across North America. One of the strongest drivers for these declines is woody plant encroachment of grasslands. In the Great Plains and Sagebrush biomes of North America, spatial covariance—a remote-sensing metric for tracking boundaries between vegetation types—is emerging as a new method to identify and strategize conservation of grassland cores in the face of woody plant encroachment. However, the relationship between spatial covariance and grassland bird community occupancy is unknown. Here, we used Bayesian multispecies occupancy models to understand how occupancy probability of six declining grassland species responded to spatial covariance at three scales (0.81, 7.29, and 65.61 ha) and tree cover in fragmented grasslands of Arkansas, USA. Model selection revealed that the smallest spatial scale (0.81 ha) best explained grassland bird occupancy. Tree cover alone was a poor predictor of grassland bird occupancy compared to models that included spatial covariance at the 0.81- and 7.29-ha scales. Grassland bird occupancy declined at tree-grass boundaries (negative spatial covariance at the 0.81-ha scale) and increased in grassland cores (near-zero or slightly positive spatial covariance at the 0.81-ha scale). At low tree cover, Dickcissel (Spiza americana), Eastern Kingbird (Tyrannus tyrannus), Loggerhead Shrike (Lanius ludovicianus), Northern Bobwhite (Colinus virginianus), and Scissor-tailed Flycatcher (Tyrannus forficatus) occupancy probability more than doubled in grassland cores (where spatial covariance approached zero). Eastern Meadowlark (Sturnella magna) had the weakest relationship with spatial covariance. Our results suggest that spatial covariance can identify grassland cores and serve as a powerful predictor of grassland bird community occupancy, even in highly fragmented grasslands. Identifying grassland cores empowers defending core grasslands from woody plant encroachment and then growing cores via active restoration.

Accumulation of per- and polyfluoroalkyl substances (PFAS) and their association with immune parameters in nestling ospreys (Pandion haliaetus) from Chesapeake and Delaware Bays, USA

Mon, 23 Mar 2026 12:53:32

Per- and polyfluoroalkyl substances (PFAS) are a class of widespread, environmentally persistent compounds that pose a potential threat to wildlife and human health. Despite recent efforts to reduce the use of long-chain PFAS in industrial practices and commercial/consumer products, the persistence and solubility of PFAS have led to their detection in wildlife on a global scale. Osprey (Pandion haliaetus) have long been used as a sentinel species with an extensive history of serving as an effective bioindicator of contamination. Here we report on a large-scale evaluation of PFAS and potential health effects in osprey from the Chesapeake and Delaware Bays, USA. In 2011 and 2015, we collected plasma samples from osprey nestlings throughout the Chesapeake and Delaware Bay watersheds. We quantified 40 PFAS congeners in osprey plasma via liquid chromatography-mass spectrometry and analyzed plasma for indicators of immune and thyroid function, and plasma biochemistry. In all birds, perfluorooctanesulfonic acid (PFOS) was the most commonly detected PFAS, followed by perfluoroundecanoic acid, (PFUnA) and perfluorodecanoic acid (PFDA). In nestling plasma from Chesapeake Bay, PFOS tended to be a higher average contributor to PFAS profiles compared to samples from Delaware Bay. In contrast, long-chain perfluoroalkyl carboxylic acids (PFCAs) such as PFUnA and PFDA comprised larger percentages of total PFAS in osprey plasma from Delaware Bay relative to Chesapeake Bay. While some PFAS concentrations were associated with plasma health indicators, the proportion of variation explained was low. Overall, our study provides a more thorough understanding of PFAS presence in the Chesapeake and Delaware Bays and is one of the first to examine whether PFAS exposure is associated with adverse health effects in wildlife.

Stream macroinvertebrate responses vary with region, land use and management practice type

Wed, 22 Apr 2026 14:34:11

Intensive land use alters hydrology and water quality, threatening freshwater benthic macroinvertebrates. Over 200,000 management practices (MPs) have been implemented across the Chesapeake Bay watershed since the 1980s, yet biological responses remain inconsistent. We synthesized 29 studies from 4 physiographic provinces covering 8 MP categories and evaluated macroinvertebrate responses along MP gradients using structural (richness), functional (biomass), tolerance, and biotic metrics. We hypothesized that MPs enhancing habitat complexity or restoring flow regimes would benefit taxa sensitive to sediment, hydrologic instability and organic pollution, with outcomes shaped by regional context, land use, and chosen metrics. Four themes emerged. (i) Agricultural Riparian Forest Buffers (RFBs) consistently improved sensitive metrics related to abundance, biomass and richness. (ii) Urban streams with Stream Habitat Improvement and Management (SHIM) showed improved richness and diversity, but biomass and tolerance metrics declined or remained neutral, indicating unresolved hydrologic and pollutant stress. (iii) Structural and functional responses diverged: effect sizes for total and feeding-group biomasses (functional metrics) were negative, whereas genus-level Ephemeroptera-Plecoptera-Trichoptera (EPT) richness (structural metric) was positive, indicating that structural shifts may not track underlying production changes. (iv) Physiographic comparisons showed counterintuitive patterns, as RFBs improved EPT richness in Piedmont streams but had negative effects in the Coastal Plain. Evaluating MP effectiveness requires distinguishing a no-MP pathway (stressors → instream conditions → assemblages → responses) from an MP-mediated pathway (practice regime → modified stressors → instream conditions → assemblages → responses), underscoring the need for region-specific, multi-metric monitoring and improved understanding of MP density thresholds and recovery lags.

Diverse novel and avian-associated viruses in the ileal viromes of northern mockingbird (Mimus polyglottos)

Thu, 2 Apr 2026 13:42:58

Viruses are the most abundant and diverse organisms on Earth, though only a small portion cause disease. Understanding viral diversity is key to understanding and predicting pathogen emergence and zoonotic spillover. Here, we use meta-transcriptomic sequencing to examine the viral communities in the ileum of 25 Northern Mockingbirds (Mimus polyglottos) from various locations across Texas. We assembled high-quality genomes of 43 viral species (40 species identified to 13 families, one to kingdom, and two to realm), 38 of which were novel. They tentatively represent avian- (n = 3), arthropod- (n = 21), plant- (n = 5) and fungi- (n = 4) associated, or other (n = 10) viruses. The arthropod-associated Dicistroviridae family was the most dominant, comprising known and potentially new species. Of potential epidemiological importance were three novel and avian-associated viruses: members of the families Hepeviridae and Picornaviridae, and a new Matryoshka RNA virus. The Matryoshka RNA virus 8 (MaRNAV-8) is sister to other Matryoshka RNA viruses, and its co-occurrence with haemosporida further supports the nested virus-parasite-vector-vertebrate host relationship of this group of viruses, with potential implications for parasite evolution, fitness and load and vector competence. The Picornaviridae virus is a member of an avian hepatovirus clade, found nested within a clade containing both the mammalian pathogens Hepatovirus A – I and the avian Tremovirus pathogens, suggestive of a newly discovered pathogen of Northern Mockingbird. Although the recovered Hepeviridae virus is of unknown pathology, its family members include the Hepatitis E viruses. With the great diversity and novelty described from ileal viromes, discriminating potential pathogens and commensal microbiota from viruses associated with food items remains challenging. A deeper understanding of virus transmission and the risk of potential zoonosis can be enhanced by tracking viruses through the food web and via inter-specific and predator-prey interactions, particular in areas subject to land-use change, where human-wildlife interactions are increased and the risks from emerging pathogens of veterinary and medical importance are more pronounced.

Mercury cycling across a U.S. semi-arid mountain ecosystem elevation gradient

Tue, 24 Mar 2026 14:10:42

Mountains comprise ∼30% of the Earth's surface, but mercury (Hg) cycling in these regions remains understudied, particularly in the semi-arid western U.S. where strong climatic and ecological gradients in mountainous landscapes influence Hg deposition, retention, and bioaccumulation. In this study, we quantified growing season inputs, storage, and bioaccumulation of Hg along a ∼2,000 m elevation gradient in the Colorado Rocky Mountains, spanning the plains to the alpine. We measured Hg in atmospheric deposition, vegetation, soil, and 12-day-old chickadees. Accounting for percent canopy cover, open precipitation was the largest component of atmospheric deposition at all elevations, followed by throughfall and litterfall fluxes. Atmospheric Hg fluxes peaked at mid-elevations, likely due to cloud-cap dynamics and denser canopy cover. Total gaseous Hg and precipitation fluxes were highest at low elevations, likely reflecting local emissions and meteorological pooling. Surface soil Hg storage was more strongly predicted by organic matter content (R² = 0.49; p < 0.01) and water retention (R² = 0.45; p < 0.01) than by elevation (R² = 0.21; p < 0.05). Alpine soils (66.3 ± 25.3 ng g⁻¹) had significantly higher total Hg concentrations than lower elevations (<41.0 ± 12.7 ng g⁻¹; p < 0.01), likely reflecting slower organic matter turnover. Soils on north-facing slopes also retained significantly higher pools of Hg in surface soils compared with south- and east-facing slopes. Vegetation Hg pools were greatest in the alpine region, likely due to long-lived plant species. Methylmercury (MeHg) concentrations in chickadee feathers peaked at mid-elevations (205 ± 155 ng g⁻¹), corresponding to higher ecosystem Hg inputs via throughfall. Our results show that deposition, canopy cover, and meteorological conditions—not elevation alone—predict Hg retention and bioaccumulation.

Terrestrial ecosystem response to changing temperature and seasonality in the Paleocene-Eocene Thermal Maximum: Shallow marine records from the Salisbury Embayment, USA

Tue, 10 Mar 2026 13:41:31

The Paleocene-Eocene thermal maximum (PETM, ∼56 Ma) is marked by a massive and rapid rise in atmospheric CO₂ and ∼5°C of global warming. It is globally characterized by a negative carbon isotope excursion (CIE), and, at least locally, is preceded by a pre-onset excursion (POE). We present palynological and bioclimatic analyses from stratigraphically expanded marginal marine sediment sections from the eastern United States. Late Paleocene forests were dominated by needle-leaved gymnosperms and broad-leaved angiosperms characteristic of warm climates. The POE is marked by a minor expansion of angiosperms and pteridophytes, warmer winters, and altered seasonal precipitation, followed by a return to pre-POE conditions. Increased terrestrial palynomorph concentrations before the CIE are suggestive of increased fluvial discharge before the PETM. Early PETM assemblages are characterized by dominance of ferns, loss of conifers, and expansion of broad-leaved angiosperm forests. Bioclimatic analyses indicate warmer mean atmospheric temperatures in early PETM time, driven primarily by winter warming of ∼3°C. A shift in seasonality, associated with increased severity of storms and floods that scoured the late Paleocene floodplain, facilitated establishment of riparian fern communities at the CIE onset. These flooding events persisted through the early part of the PETM and were severe enough to transport Westphalian-age (Middle Pennsylvanian) reworked material from the central Appalachian Basin and flush large amounts of terrestrial material and carbon onto the continental shelf, resulting in decreased salinity, increased productivity, and water-column stratification.

Boxed in or branching out? Movement and resource selection of eastern box turtles (Terrapene carolina carolina) in an urban green space

Jones, Max; Ferebee, Kenneth; Ford, W. Mark; Hunter, Elizabeth — Thu, 2 Apr 2026 18:18:03

The eastern box turtle (Terrapene carolina carolina) is a long-lived terrestrial turtle species distributed throughout the eastern United States that has experienced widespread population decline. Many eastern box turtle populations are persisting as remanent populations in small, fragmented urban green spaces. We investigated the movement and resource selection of eastern box turtles within a mid-Atlantic region urban forest in the eastern United States. We used a combination of turtle occurrence data (via visual encounter surveys) and radio telemetry to create resource selection functions. Additionally, we applied a simulation modeling approach and modeled activity areas via dynamic Brownian Bridge Movement Models to quantify interactions between turtles and roads or trails. We also used these models to determine the propensity for turtles to move outside of the managed urban forest boundary and into surrounding development. We observed that turtles selected for deciduous forest patches and avoided roads and trails despite the urban forest having very little available areas where anthropogenic features could be avoided. We also demonstrated observed (and probable) movements outside of the urban forest boundary. Although eastern box turtles are persisting within the urban green space we examined, our work determined that interactions with roads and trails, and movements outside of protected boundaries into developed areas present challenges to individuals navigating the urban forest.

Magnetic storms and geoelectric hazards

Thu, 26 Mar 2026 18:08:52

Magnetic storms induce geoelectric fields at Earth's surface that can interfere with grounded long-line systems. The September 1859 storm disrupted global telegraph operations, the March 1989 storm caused a blackout in Canada and interfered with electric-power-transmission systems in the United States, and other storms have had related impacts. The geographic and temporal dependence of geoelectric fields are functions of both geomagnetic variation and local surface impedance, which differ considerably across different geological regions. These dependencies can be mapped across the contiguous United States by combining magnetotelluric impedance tensors with ground magnetometer time series. This review illustrates such mapping for the 1989 storm and shows that power-system interference was experienced where surface impedance is high, and when and where geoelectric fields were intense. Statistical analyses indicate that storms comparable to that of March 1989 occur roughly once every four solar cycles. Ongoing developments in numerical modeling and real-time monitoring are anticipated to enable prediction of geoelectric hazards.

▪ Magnetic storms can induced electric fields in the solid Earth that interfere with electric-power-transmission systems.
▪ Geoelectric hazards depend on the storm-time geomagnetic disturbance and the electrical conductivity structure of Earth.
▪ Historically, impacts on telecommunication and power-transmission systems in the United States have been concentrated in the East and Midwest.
▪ The future occurrence of a magnetic superstorm could cause widespread disruption of electric-power-transmission systems.