USGS ScienceSearch

SEARCH · USGS Science

Results for “JAMA Network Open”

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.

347 records · Page 18Linked to original sources

Scientific opportunities in the National Landscape Conservation System

The National Landscape Conservation System consists of unique and beautiful places across America’s landscapes where identified resources and values are protected and science is highlighted. The mission of the National Landscape Conservation System (NLCS), which is managed by the Bureau of Land Management and is often referred to as the agency’s National Conservation Lands, is to conserve, protect, and restore nationally significant landscapes for their cultural, ecological, and scientific values. This clear inclusion of science in the NLCS mission sets the stage for individual units to serve as places of learning, teaching, discovery, and innovation. Science is an integral part of managing the National Conservation Lands, and science conducted within and across the more than 900 units that make up the NLCS can inform and influence conservation and public land management well beyond its boundaries. Here, we highlight seven core aspects of National Conservation Lands that present valuable science opportunities: (1) the scientific values for which individual units are designated; (2) the many other resources, objects, and values within units; (3) the value of units as “control” sites for understanding the effects of activities such as mineral extraction that commonly occur elsewhere on multiple-use public lands but are often prohibited within National Conservation Lands; (4) the value of units for studying the effects of activities such as recreation that regularly occur and may be intensified on National Conservation Lands; (5) the high visibility of units, which draws strong interest and engagement from scientists, partners, and the public; (6) the functioning of the units as a network managed for a common purpose, which provides an opportunity to explore cross-cutting science questions across widely varying contexts and geographies; and (7) the opportunities units provide to promote and apply Indigenous Knowledge to scientific research to manage natural and cultural resources. Because of all of these characteristics, National Conservation Lands can serve as hubs for basic and applied science that can inform management of all public lands and resources into the future. We highlight these science opportunities through examples from existing units and suggest two actions that could help further science activities and impact on National Conservation Lands.

Parks Stewardship Forum

Influence of inherited structure on flexural extension in foreland basin systems: Evidence from the northern Arkoma basin and southern Ozark dome, USA

Extensional faults are key components of foreland basin systems. They form within the upper crust in response to flexure of the lithosphere and accommodate subsidence within the foredeep and forebulge depozones. Such faults are excellent proxies for orogenic system evolution and control the distribution of natural resources and hazards. However, the spatiotemporal evolution of flexural extension has not been documented previously at a regional scale, thereby limiting our understanding of underlying geodynamic controls. Here, we resolve late Paleozoic flexural extension in the northern Arkoma basin and southern Ozark dome, USA. We synthesize a large database of previous mapping, existing research, subsurface data, and geophysical data into 3D geologic and 2D kinematic models. Mesh surfaces representing several key horizons from the Carboniferous Period (ca. 335-306 Ma) were constructed. These surfaces were built from oil and gas well tops (n = ∼10,000) and surface geologic map contacts using an advanced kriging method. The mesh surfaces are offset by a complex 3D fault network, allowing detailed analysis of along-strike and down-dip variations in fault displacement. Analysis of the 3D model reveals a regular and repeated fault segmentation pattern wherein E -W striking, left- and foreland-stepping en échelon normal faults are segmented by inherited NE striking basement faults. Maximum vertical separation along the E -W normal faults is generally focused between the inherited NE-trending faults. This suggests that the inherited basement faults delocalized extensional strain during late Paleozoic normal faulting. Maximum vertical separation and fault localization may correlate to areas with high-amplitude positive magnetic anomalies interpreted as Mesoproterozoic granitic rocks. Speculative covariance of magnetic anomalies and fault displacements implies that the relatively strong basement granite concentrated stress, leading to localized faulting within the relatively thin sedimentary cover. Lastly, we show that flexural extension migrated southeast to northwest from the Chesterian-Morrowan (ca. 335-319 Ma) to the Desmoinesian (ca. 306 Ma). The migratory flexural extension may be explained by diachronous loading during Pangean assembly, or by synchronous loading but variable load compensation due to inherent factors.

Arkansas, Oklahoma

Groundwater structures fish growth and production across a riverscape

Landscapes are composed of habitat patches and conditions that vary across space and time. While habitat variability and complexity can support important ecological processes and ecosystem services, the dynamic nature of habitats can also constrain organismal growth and production as optimal conditions are fleeting. In riverine ecosystems, groundwater discharge to streams stabilises water temperature and flow regimes, thus mediating how habitat complexity is expressed. Yet, how stable habitats structure growth and production within the broader landscape matrix is not well understood. In this study, we explored the effects of groundwater on spatiotemporal variation in growth and production for juvenile Yellowstone cutthroat trout ( Oncorhynchus virginalis bouvieri ) across the upper Snake River catchment, Wyoming, USA. We combined machine learning techniques and remotely sensed landscape data to estimate groundwater availability across the river network, which we linked to stream temperature regimes and conspecific density. We then used Bayesian hierarchical models to quantify the effects of temperature, density and groundwater on spatiotemporal variation in fish growth and production in 52 focal reaches. Finally, we predicted body size trajectories and trends in total production continuously over both space and time to understand the effect of groundwater at the riverscape scale. Groundwater discharged to streams where topography changes abruptly in valley-bottom areas underlain by coarse glacial deposits. Groundwater stabilised temperature regimes and was associated with high trout densities. Temperature and density, in turn, interacted to influence growth rates: growth increased strongly with temperature, but this effect was reduced when density was high. Accordingly, variation in groundwater availability among stream reaches diversified growth and production regimes. In reaches with low groundwater availability, growth and production declined over time from summer maxima. In contrast, in reaches with high groundwater availability, temporal trends in growth and production were hump-shaped—peaking in autumn—and mean production was greater. At the riverscape scale, temporal asynchrony in growth rates generated convergent spatial variation in growth capacity, but—when combined with density—led to the formation of distinct hotspots of production. Our results demonstrate how groundwater, an important driver of aquatic ecosystem heterogeneity, structures trout growth and production across space and time. Importantly, rare, but stable habitats may disproportionately affect ecological processes and serve as key sources of population diversity at larger spatial scales.

Wyoming

Geologic map of MTM −15032 and −20032 quadrangles, western Ladon basin, Mars

Our mapping efforts focused on the geomorphic terrains and geologic units contained within Mars Transverse Mercator (MTM) −15032 and −20032 quadrangles. These two quadrangles are located along the west side of Ladon basin and span lat −12.5° N. to −22.5° N. and long 325° E. to 330° E. The western part of Ladon basin and its bounding basin ring structures to the west preserved features that help to understand the long history of drainage across the Margaritifer Terra region of Mars. Our geologic map provides new insight into the extent of aqueous activity across western Ladon basin and the adjoining highlands, which includes the northern part of Ladon Valles, Arda Valles, numerous small valleys in the western highlands, and light-toned layered deposits associated with these fluvial features. The Mars Odyssey Thermal Emission Imaging System (THEMIS) infrared (IR) daytime mosaic (100 meters per pixel) was used as the primary base map. We constructed the geologic map of western Ladon basin at 1:1,000,000 scale. We identified 20 geologic units in the map area, which we divided into the following groups: crater units; volcanic units; chaotic units; basin fill units; crater, valley, and channel units; and plateau and highlands units. We observed the following sequence of events: (1) formation of Ladon and Holden impact basins in the middle Noachian, producing a mountainous unit; (2) landscape degradation and infilling during the Late Noachian to Late Hesperian forming terra units; (3) formation of Ladon Valles in the Late Noachian to Early Hesperian by catastrophic flooding, producing channel units; (4) accumulation of sediments in the Late Noachian to Early Amazonian, forming basin fill units; (5) formation of a smooth crater fill unit during the Early to Late Hesperian; (6) deposition during the Late Hesperian to Early Amazonian, creating light-toned layered units within Ladon Valles, Ladon basin, and other smaller valley networks along the western uplands, as well as formation of a light-toned unit along some crater floors; (7) formation a of chaotic unit in the Early to Middle Amazonian, formation of alluvial fans in the Early Amazonian, and eruption of a volcanic unit in the Middle Amazonian; and (8) formation of craters throughout the geologic history of the map region.

Scientific Investigations Map

Ductile and brittle Rio Grande Rift deformation in Oligocene granite records a two-stage rift history in southern Colorado

The timing and nature of early deformation in the Rio Grande Rift remains poorly constrained. We present evidence for the earliest structural signature of rift extension in the Sangre de Cristo Range, southern Colorado, based on new geologic mapping, structural analysis, rock magnetic data, and thermochronology. These analyses focus on the ~30.0 Ma granite of Chokecherry Canyon, which hosts discrete low-angle mylonitic shear zones and a distributed, gently SW-dipping protomylonitic fabric. Incremental stretching axes, stretching lineations, and Kmax magnetic lineations plunge gently WSW. Quartz microstructures and crystallographic orientations indicate dominantly coaxial strain in the protomylonite and general shear in the discrete shear zones. Quartz c-axis opening-angle thermometry suggests deformation at ~420–540°C. Thermal modeling of ⁴⁰Ar/³⁹Ar K-feldspar data indicates rapid post magmatic cooling below the brittle–plastic transition, supporting shear-zone formation immediately after emplacement. Slow cooling from ~20–13 Ma was followed by renewed rapid cooling at ~13 Ma, interpreted as the onset of extensional exhumation along the Sangre de Cristo Fault System. These results show that extension in the northern Rio Grande Rift was active by ~30 Ma, earlier than previously recognized. We propose a two-stage model for northern Rio Grande Rift evolution: Stage I (30–23 Ma) records ENE–WSW extension localized in low-angle mylonitic shear zones associated with mid-crustal intrusions; Stage II (≤18 Ma) reflects brittle high-angle normal faulting, focused exhumation, and rift narrowing. Stage I magmatism and deformation along the western range front likely established crustal weaknesses that guided later fault development.

Colorado

Comparative assessment of STIC sensors, streamflow and rain gauges for quantifying river connectivity in intermittent systems

In intermittent stream systems, including those occurring in Texas, USA, the severity of low-flow conditions, duration of seasonal disconnection, and frequency of no-flow events have been amplified by drought. Documentation of these no-flow events is necessary to evaluate ecosystem health. However, many intermittent reaches remain un-gauged given that perennial river sec-tions are often prioritized for gauge placement. Our objectives were to 1) document stream flow using Stream Temperature, Intermittency, and Conductivity (STIC) loggers to determine the frequency and duration of no-flow events in intermittent tributaries of the Colorado River, Texas and 2) compare logger data to publicly available data from streamflow discharge and precipitation gauge networks to understand differences among these data types for drying event characterization. We use these comparisons to summarize benefits and limitations of the application of in-stream data loggers. STIC loggers were deployed at 19 sites, one in each pool and riffle habitat of a stream reach. STIC loggers recorded a measurement of relative conductance every six hours from June 2022 to March 2024, which was used to determine the presence or absence of flow connectivity in a reach. No-flow duration among intermittent reaches varied between 37 and 270 days across tributaries during an ongoing drought in the study area. Overall, logger data was more precise than discharge data for characterizing no-flow events or precipitation data when documenting presence of water in the stream channel due to runoff. Lack of discharge gauges in intermittent tributaries left large sections of stream reaches undocumented and resulted in mischaracterization of flow patterns. Drought severity across the tributaries did not follow longitudinal patterns that would be expected by the climatic precipitation gradient of the study area. More research is needed to determine if factors such as population size affect severity. Likewise, precipitation data did not correlate well with logger water presence data, lacking consideration for groundwater recharge, soil hydrophobicity, and surface compaction. This study shows that to monitor no-flow events, detailed spatial datasets are necessary and that STIC loggers are useful tools that provide data to fill spatial information gaps and facilitate more accurate flow characterization and water presence data in intermittent systems.

Texas

Earthquake-triggered ground-failure inventory associated with the M7.1 2018 Southcentral Alaska earthquake

The 30 November 2018, magnitude (Mw) 7.1 earthquake in Southcentral Alaska triggered substantial landslides, liquefaction, and ground cracking throughout the region, resulting in widespread geotechnical damage to buildings and infrastructure. Despite a challenging reconnaissance and remote-sensing environment, we constructed a detailed digital inventory of ground failure associated with the event from several sources. Sources included information derived from remotely sensed data, and data compiled from literature, social media postings, and earthquake damage information compiled by local, state, and federal agencies. Each instance of ground failure within the inventory contains information on the location and type of observed ground failure, and the methods and data used to document the occurrence. Where high-quality data, such as LIDAR or satellite imagery, were available and showed the ground-failure instance clearly, the extent is mapped as a polygon or polyline. All other locations are mapped as points. There are a total of 886 ground-failure instances documented within the inventory (400 landslides, 286 liquefaction features, and 200 features unattributed to specific processes). A semi-quantitative confidence scheme is used to describe mapping certainty associated with each ground-failure feature. This inventory represents a relatively moderate ground-failure-triggering event that occurred in a subarctic environment. This data paper describes the content within the inventory, the inventory data collection procedures, and limitations of the data. Events of this type are not often documented in detail; thus, adding the inventory data to the US Geological Survey Open Repository of Earthquake-Triggered Ground-Failure Inventories further diversifies the datasets available to the scientific community to be used to better understand and model earthquake-triggered ground failure.

Alaska

High-resolution transboundary vegetation community maps of the Sonoran and Mojave Desert ecoregion to support critical landscape conservation planning and habitat management needs

We produced a 30-m resolution binational land cover map of Bird Conservation Region 33 (BCR 33) for the U.S. North American Bird Conservation Initiative. The region covers large portions of the Sonoran and Mojave Deserts. The map can support the U.S. Fish and Wildlife Service (FWS) Migratory Bird Program’s recovery planning efforts and constitutes the first known binational land cover dataset spanning sections of the United States–Mexico border and using a consistent classification system for both countries. The mapped region includes 152 distinct land cover classes, covering a total area of 38,421,453 ha (148,345 mi 2 ), of which 13,148,345 ha (52,706 mi 2 ) are located in Mexico and 24,770,640 ha (95,639 mi 2 ) in the United States. We primarily used Landsat 8 (OLI) imagery, supplemented by limited ground surveys from two field campaigns, drone-based aerial data, and existing vegetation classification frameworks from both countries. The classification applied a data-fusion approach integrating 30-m Landsat 8 imagery, decadal phenology metrics from vegetation indices, and a random forest model trained mainly with datasets from a comprehensive national mapping project from the U.S. Geological Survey (USGS) GAP Analysis Project (GAP) and federal wildland fire agencies’ Landscape Fire and Resource Management Planning Tools (LANDFIRE) (GAP/LANDFIRE) [United States side] and the National Institute of Statistics and Geography (INEGI) [Mexico side] as well as land cover maps and opportunistic open-access and field observations. Mapping of the full BCR 33 region was carried out in two phases: 1) Phase I, the prototype map, covered a smaller portion of the transboundary area and identified 31 land cover classes, and 2) Phase II, the full BCR 33 map (refer to Figure 1), which resulted in 152 land cover classes. Using a Random Forest classifier, we achieved an overall prediction accuracy of 92% for the Phase I map and 87% for the Phase II full region map. This slight decrease can be attributed to working on a larger, more complex area with a greater number of land cover classes. No formal validation was conducted, aside from using a subset of the collected field observations and training data to assess model performance during and after training. The training sites were further verified using Google Earth (Google, 2026) imagery. Two undergraduate students who worked for over a year visually inspected imagery and open access public images to confirm each training site during model training using in-house developed, online, visual tools. A portion of this field training data was reserved for model validation, and the corresponding results are to be presented in later sections. The project developed an end-to-end, medium- and fine-resolution remote sensing–based data fusion mapping approach. This effort produced a map (Nagler et al., 2025) and the online tools to support a dynamic, live, online map for visualizing the transboundary vegetation communities in BCR 33. The toolset is currently hosted by the University of Arizona (UofA) Vegetation Index and Phenology (VIP) Lab to support FWS partners (https://vip.arizona.edu/viplab_data_explorer?LCM_BCR33). The online map is designed to allow rapid updates using new training, validation, or correction data, making it dynamic and maintainable. The approach we took established a framework for rapid updating and correction of land cover maps, as the model can be quickly retrained with new field observations, updated training data, or other sources. This enables dynamic mapping and change detection of the region’s vegetation. This framework is an advance in data fusion and crowdsourced mapping of complex, vulnerable regions, providing support to regional stakeholders and the wider user community. This transboundary map can inform the protection, conservation, and restoration of vegetation, habitat, and ecosystems, particularly for threatened and endangered species across the two nations using consistent and harmonized binational mapping systems. Beyond supporting land management decisions and stakeholders in the transboundary desert ecoregions, this BCR 33 mapping effort establishes a foundation for future rapid, low-cost, cross-border land cover mapping that can benefit and advance ecosystem management.

Arizona, Baja California, California, Nevada, Sina

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

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

Connecticut, Massachusetts, New York

The impact of burial diagenesis on soil-formed minerals in paleosols using stable isotopes of phyllosilicates and carbonate clumped isotopes

To understand the effects of burial diagenesis on the stable isotope geochemistry of soil-formed clay and carbonate minerals in paleosols, samples were collected from seven cores, spanning middle- to upper-Pennsylvanian strata of the Illinois Basin, with varied maximum burial depths of 1–3 km. Mixed-layer illite-smectite and kaolinite mixtures give δ 2 H and δ 18 O values of −83 ‰ to −36 ‰ and 11.9 ‰ to 21.1 ‰ (VSMOW), respectively. After carbonates were screened petrographically for diagenetic textures using transmitted light and cathodoluminescence, measured clumped isotope Δ 47 values range from 0.504 to 0.563 ‰ (I-CDES). Resulting mineral formation temperatures for phyllosilicate mineral mixtures are 28 to 66 °C (mean = 47 °C), whereas T(Δ 47 ) estimates for calcites are 36 to 61 °C (mean = 45 °C). Calculated δ 18 O water values from which phyllosilicate minerals and calcites precipitated under isotopic equilibrium ranges from −7.1 to −1.2 ‰ and − 1.4 to +4.9 ‰, respectively. Closed and open-system phyllosilicate-fluid exchange modeling indicates that phyllosilicate alteration occurred in the presence of a low temperature brine or meteoric water and is interpreted to occur in a layer-by-layer illitization transformation. Due to the lack of diagenetic textures and positively correlated T(Δ 47 ) and δ 18 O water , calcites are interpreted to have undergone solid-state bond reordering. Despite low to moderate temperatures (<125 °C) and varying depths of shallow burial (1–3 km), solid-state transformation of phyllosilicates and calcites indicates paleosols had prolonged exposure to burial conditions which has implications for the use of paleosol minerals for paleoenvironmental reconstructions.

Illinois, Indiana, Kentucky

Biome-scale spatial patterns of avian abundance reveal proactive conservation opportunities in North American grasslands

North American grassland birds have experienced steeper population declines than any other avian guild, yet conservation efforts remain largely reactive and fragmented. We used nearly four decades of North American Breeding Bird Survey data to identify biome-scale spatial patterns (clustering) of grassland bird abundance for the Great Plains. Our results reveal an ecological core in the north-central Plains where community-level abundance is either increasing by >100% or remains high and stable, providing a strategic roadmap for a “Defend the Core” conservation approach. This approach flips the script from reactive triage centered on isolated population fragments to a proactive strategy of maintaining large-scale ecosystem integrity. Conversely, we found that population losses are more spatially clustered than wins, reflecting the relentless, one-way movement of woody encroachment and agricultural conversion. This asymmetry supports prioritizing intact landscapes, as current restoration rates are often outpaced by the scale of habitat loss. Notably, we found that community-level spatial clustering is a more robust indicator of biome condition than trends of individual flagship species, suggesting that managing for ecosystem integrity provides a more effective multi-species umbrella. Given our results, there is an opportunity for operationalizing a Great Plains Conservation Design that is ecosystem-centric and rooted in the sustainability of the private-land cattle production that maintains these open spaces. By leveraging avian abundance as a biological sensor, managers and producers can deploy a shared vision that matches the spatial scale of the threats, moving from reactive triage to proactive defense of core working grasslands in North America.

Great Plains biome

Reassessing water availability in the Nile River Basin using satellite-based irrigation water use accounting

Existing estimates of irrigation areas and water use in the Nile River Basin vary widely and are based on inconsistent data and methodologies. This study leverages advances in remote sensing to provide an up-to-date and consistent basin-wide reassessment of water use in the Nile River Basin. Land cover and actual evapotranspiration data from 2013 to 2022 were used to quantify irrigated areas, irrigation water use, and the average naturalized water yield. Across the basin, 7.1 million hectares of irrigated land consume an average of 72.5 ± 2.6 billion cubic meters (BCM) of blue water annually. Egypt and Sudan together account for 92% of the irrigated area and 97% of the irrigation water use. The basin’s average naturalized yield, evaluated at Dongola station, is determined to be between 108 and 114.4 ± 3.2 BCM. Although conservative, this updated analysis indicates the basin’s average naturalized yield is 14%–20% higher than the previous naturalized yield estimate (95 BCM) and 29%–36% higher than the commonly cited observed flow of 84 BCM, which is frequently treated as the total basin-wide available water. Importantly, the study underscores that the 84 BCM represents flow at the Aswan gauging station after upstream consumptive uses, not the total water available in the basin under naturalized conditions. This distinction is critical for accurate water accounting, planning, and governance in this water-scarce basin. This study demonstrates the value of open-source remote sensing resources in data-scarce regions while emphasizing the need for region specific validation and bias correction for improved accuracy.

Nile River Basin

Identifying precursors and tracking pulses of magma ascent in multidisciplinary data during the 2018–2023 phreatomagmatic eruption at Semisopochnoi Island, Alaska

The 2018–2023 phreatomagmatic eruptions at Semisopochnoi Island, Alaska produced abundant long-period (LP) seismicity, harmonic and broadband tremor, and explosion signals over several well-monitored periods of eruption and quiescence. The corresponding dataset provides an excellent opportunity to investigate precursory and syn -eruptive geophysical signals of long-lived phreatomagmatic eruptions using multiparameter observations. We generated explosion and LP event catalogs through novel implementations of the REDPy ( Hotovec-Ellis, 2024 ) repeating event detector in mid-2021 following a network upgrade and the onset of a new phase of the eruption. The hundreds of detected explosions show a high degree of infrasound waveform similarity over more than a year, indicating a repeating source mechanism likely associated with explosive magma-water interaction. The seismic LP catalog shows that events began over a month prior to renewed explosive activity at the beginning of August 2021, and that lower frequency index (FI) LPs were generated in the week prior to the onset of explosions. We applied a recently developed machine learning tool (VOISS-Net, Tan et al., 2024 ) to catalog abundant broadband and harmonic seismic tremor recorded before and during the renewed explosive activity, along with LPs and explosions. The tremor catalogs complement the LP and explosion catalogs by filling out the seismic sequence with the dominant signal types. Together, these catalogs reveal a seismic sequence of renewed unrest that started with several weeks of LP events, followed by LPs with lower FI values and harmonic tremor in the days prior to explosive activity, and finally the onset of discrete explosions and broadband eruption tremor. We interpret this sequence as the ascent of a new pulse of magma that first interacted with the hydrothermal/groundwater system to produce LPs, followed by harmonic tremor, and that ultimately drove explosive magma-water interactions and periods of continuous ash emissions. The 2021 seismic sequence, in combination with long-term records of satellite SO 2 emissions, deformation from interferometric synthetic aperture radar (InSAR) analysis, ash sample analysis, infrasound, and volcano tectonic seismicity, allows us to interpret the entire 9-year period of unrest and eruption that began with an intrusion and earthquake swarm in 2014.

Alaska

Nutrients and pesticides in ground water of the Ozark Plateaus in Arkansas, Kansas, Missouri, and Oklahoma

A total of 229 ground-water samples were collected from 215 sites as part of the Ozark Plateaus study unit of the National Water-Quality Assessment Program. These samples were collected from 1993 through 1995 using a network of springs and wells with three scale-dependent components. The first component, the study-unit survey, consisted of 99 randomly selected springs and domestic wells in the Springfield Plateau and Ozark aquifers. The second component, two land-use studies, consisted of 42 springs and domestic wells in a poultry-dominated agricultural area and 40 springs and domestic wells in a cattle-dominated agricultural area overlying the Springfield Plateau aquifer. The third component, the small-watershed study, consisted of 4 springs, 18 domestic wells, and 11 monitoring wells in a small basin within the poultry land-use study area. Samples were analyzed for major ions, nutrients, dissolved organic carbon, methylene blue active substances, tritium, and 88 pesticides and metabolites.The water-quality data from these samples were analyzed with descriptive and statistical methods. Nitrite plus nitrate, which was detected more often and in greater concentrations than any of the other nutrients, ranged from less than 0.05 to 25 milligrams per liter as nitrogen. Nitrite plus nitrate concentrations positively correlated to percent agricultural land use around each site. Median nitrite plus nitrate concentrations generally were greater in samples from springs than in samples from wells. Concentrations of nitrite, ammonia, and ammonia plus organic nitrogen were also affected by land use and also by concentrations of dissolved oxygen in the ground water. Concentrations of phosphorus and orthophosphate probably were affected by land use and also by phosphorus solubility. Pesticides were detected in 80 of 229 samples from 73 of 215 sites. A total of 20 pesticides were detected with a maximum of 5 pesticides detected in any 1 sample. The most commonly detected pesticides were tebuthiuron, atrazine, prometon, desethylatrazine, and simazine. Maximum concentrations ranged from 0.003 to 1.0 microgram per liter. The occurrence and distribution of pesticides were related to land use. Percent agricultural land use was greater for samples with pesticides detected than for samples with no pesticides detected. Pesticides were detected more often in samples from springs than in samples from wells. The occurrence of pesticides also was related to seasonality and chemical characteristics, such as solubility and persistence, of the compounds.

Arkansas, Kansas, Missouri, Oklahoma

Evolution of permeability and strength recovery of shear fracture under hydrothermal conditions

Geothermal energy is a clean and renewable resource that depends on the ability to move water through hot rock. In many locations, the ability to move water through rock requires the presence of extensive natural or human-made systems of fractures. However, these fracture systems are influenced by a variety of complex processes that occur at the temperature and pressure conditions found in geothermal reservoirs. To make geothermal systems more efficient, it is important to understand how these fractures evolve over time in response to these processes. This study explored how fractures in rock change under high-temperature and pressure conditions, like those found in geothermal reservoirs, through laboratory experiments and numerical simulations. In general fractures, tend to close due to pressure from the surrounding rock preventing the flow of water. One set of laboratory experiments revealed that fractures in granitic rock can weaken over time at elevated temperatures. These weak fractures can then slip slowly potentially keeping the fracture open for a longer period. This study also showed that at the highest temperature examined (250 °C) large differences between the chemical composition of the rock and the fluid moving through it could improve water flow temporarily likely due to the water dissolving parts of the rock. Flow rate predicted by models shows overall good agreement with the flow rate measured in our experiments but with some long time-scale variations that may be due to mechanical wear of the fracture surface. A second set of experiments found that, over the temperature range examined, the closure of fractures was caused by brittle failure driven by high stresses at points where fracture surfaces were in contact. This result suggests that fracture closure in geothermal reservoirs could be counteracted by maintaining high pore pressures. Additionally, we also tested a new method using electrical resistance to measure fracture closure in controlled systems with high precision.

Final Technical Report

Seismic response of rock towers at the Trona Pinnacles (U.S.A.) to the 2019 Ridgecrest earthquake sequence: Theory, observations, and models

We analyze the seismic response of a class of fragile geologic features (FGFs), referred to as rock towers (RTs) at the Trona Pinnacles, a group of RTs in southern California that suffered strong shaking during the 2019 Ridgecrest earthquake sequence. FGFs, including RTs, may provide maximum constraints on past earthquake shaking intensity, and thereby support probabilistic seismic hazard assessments (PSHAs). In a rare case study, we explore the hypothesis that RT structural integrity is time dependent, as damage accumulates progressively. We develop finite‐element method (FEM) models of the RTs using photogrammetric shape models. We validate the models by comparing numerical simulations of their response to broadband ground shaking with low‐intensity seismic recordings obtained at the Pinnacles. Results of our simulations are in good agreement with the seismic recordings of actual earthquake aftershocks. We next use the results of the FEM models to analyze the response and evolution of RTs. Our analyses elucidate the influence of geometry over their seismic response, providing a rationale that may explain the rarity of slender RTs at Trona: high‐aspect‐ratio structures that respond in bending develop detrimental tensile stresses that crack the rock, whereas low‐aspect‐ratio ones’ response also includes shearing, which does not compromise material integrity as much as tension. Field measurements with a rebound hammer support this finding, suggesting that the material around the base of slender rocks has been weakened relative to other parts of the RT. We also study how to define simplified mechanical models (“archetypes”) to predict the natural frequencies of RTs. Results from our work illuminate the fundamental mechanisms of seismic response and progressive failure of RTs, and open new avenues of research to potentially incorporate these geologic features as long‐return period constraints on PSHA, in ways analogous to those of the widely used precariously balanced rocks.

California

Using a temporary emigration model to estimate abundance of stream fishes from hybrid removal surveys with and without block nets

Monitoring programs are often faced with a decision to allocate resources into either robust spatiotemporal coverage to estimate a population index (e.g., not true abundance) or confirming closed sampling conditions (e.g., with block nets) for an unbiased population estimate at the cost of spatiotemporal coverage. However, making accurate and precise abundance estimates at robust spatiotemporal scales is possible when combining open and closed sampling designs with integrated modeling techniques. We used simulations and a case study of backpack electrofishing surveys in the Santa Ana River, California to test the efficacy of an integrated abundance model (temporary emigration model, TE) to estimate abundance of fishes using removal sampling methods with a hybrid sampling design (sampling with and without block nets during removal sampling). We found that the TE model performed well under most modeling scenarios (sample size, amount of closure violation, number of samples collected during closure), although at least a few samples with block nets were necessary for all parameters to be estimable. When applied to fish surveys in the Santa Ana River, we found that catch of the fishes fit to the TE model (Santa Ana Sucker, Arroyo Chub, Channel Catfish, Largemouth Bass, Yellow Bullhead) showed little evidence that the closure assumption was violated when block nets were not used. Additionally, we found that the abundance of non-native fishes negatively affected the abundance of the native Santa Ana Sucker, which was also found to adversely affect the native fish’s access to critical habitat consisting of gravel and cobble substrate. Our results indicate that the TE model presents a viable solution to common sampling problems that impact many monitoring programs, where precise and accurate population estimates can be made at large spatiotemporal scales even when most samples violate the closure assumption.

California

Pediment formation and subsequent erosion in Gale crater: Clues to the climate history of Mars

Evidence of paleo-rivers, fans, deltas, lakes, and channel networks across Mars has prompted much debate about what climate conditions would permit the formation of these surface water derived features. Pediments, gently sloping erosional surfaces of low relief developed in bedrock, have also been identified on Mars. On Earth, these erosional landforms, often thought to be created by overland flow and shallow channelized runoff, are typically capped by fluvial sediments, and thus in exceptionally arid regions, pediments are interpreted to record past wet periods. Here we document the Greenheugh pediment in Gale crater, exploiting the observational capability of the Curiosity rover. The absence of a fluvial cap suggests that the pediment was likely cut by wind erosion, not fluvial processes. The pediment was then buried by an aeolian deposit (Stimson sandstone) that mantled the lower footslopes of Aeolis Mons (informally known as Mt. Sharp). This burial terminated active wind erosion, preserving the pediment surface (as an angular unconformity). Groundwater was present prior-to, during, and shortly after Stimson deposition, perhaps contributing to lithification and certainly to early diagenesis. Post lithification, wind erosion first cut canyons in the northern most footslopes (north of Vera Rubin ridge). Unlithified gravels were deposited in these canyons, likely due to runoff from Mt. Sharp. Boulder-rich fluvial and debris flow deposits built a > 70 m thick sequence (Gediz Vallis ridge) on the southern Greenheugh pediment. Continued wind erosion left elevated patches of gravel on the northern footslopes, and exposure age dating shows that erosion essentially ceased before 1 Ga (but possibly much earlier). Erosion to the south led to emergence of Vera Rubin ridge, retreat of the Greenheugh pediment, and the formation of Glen Torridon valley. Hence, this footslope environment of Mt. Sharp records climate-driven periods of wind erosion, aeolian deposition (and groundwater activity), surface runoff and sediment deposition, followed by further significant wind erosion that declined to present very slow rates. This likely occurred during the late Hesperian and possibly into the Amazonian.

Icarus