USGS Science⌕ Search

SEARCH · USGS Science

Results for “Hydrology and Earth System Sciences (HESS)”

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

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

1,661 records · Page 78Linked to original sources
USGS2026

Site response models based on geometric parameters for southern California sedimentary basins

Site response in sedimentary basins is influenced by complex three-dimensional (3D) features, including trapping of seismic waves, focusing of seismic energy and basin resonance. Current ground motion models (GMMs) incorporate basin effects using one-dimensional parameters like V S30 and shear wave velocity isosurface depths, which are limited in capturing lateral and 3D effects. To address these limitations, we develop seismic site response models based on novel parameters that represent multi-dimensional properties of the Los Angeles Basin (LAB) geometry and shear wave velocity. We define a basin shape for the LAB using depth to subsurface geologic interfaces associated with the oldest sedimentary deposits (depth to a particular shear wave velocity horizon, i.e., 1.5 km/s - z 1.5 ) and the depth to the crystalline basement ( z cb ) which are determined using geologic cross sections and community seismic velocity model profiles. We explore a suite of geometric descriptors computed for the LAB and southern California, from which three parameters with the greatest predictive potential are selected and evaluated using empirical ground motion residual analyses in combination with the Boore et al. GMM. The results demonstrate that the zonal heterogeneity index ( ), standard deviation of the absolute difference between z 1.5 and z cb ( ) and standard deviation of z cb ( ) each provide a reduction in site-to-site variability ( ϕ S2S ) of empirical GMMs. The reduction in ϕ S2S is period-dependent, with average decreases of 3%, 26% and 6% for , , and , respectively. Although these reductions are modest from an engineering application perspective, they are statistically significant, underscoring the inherent difficulty in fully characterising complex basin effects. Collectively, these findings indicate that the inclusion of basin-specific geometric parameters yields measurable, albeit incremental, improvements in site response prediction and establishes a framework for the progressive refinement of seismic hazard characterisation within sedimentary basins.

California↗
USGS2025

Revisiting an enigma on California's north coast: The Mw6.5 Fickle Hill earthquake of 21 December 1954

Many earthquakes occur along the North Coast of California in the vicinity of the Mendocino Triple Junction (MTJ), where the Pacific, Gorda, and North American (NA) plates meet, and on the adjacent plate boundaries. The MTJ marks the nexus of the Mendocino and San Andreas faults with the Cascadia subduction zone (CSZ). Historically, most large earthquakes around the MTJ have been within the offshore Gorda plate and its subducted portion beneath the NA plate. North of the MTJ, active faults mapped in the NA plate are part of the CSZ fold‐and‐thrust belt. Although some events have been detected in the NA plate, no large historic events have been associated with mapped surface faults. The 21 December 1954 M w 6.5 earthquake in Humboldt County is one possible exception. Using published data from catalogs and articles, unpublished data from Berkeley’s archives, and S ‐ P times interpreted from two U.S. Coast and Geodetic Survey (USCGS) accelerometers, we determine a probability cloud for the earthquake’s hypocenter using NonLinLoc. The highest probability location lies beneath Fickle Hill just east of the city of Arcata, California, at 40.87° N, 124.03° W, and ∼11 km depth. Using P ‐wave polarities from Berkeley stations and the digitized waveforms from the accelerometers, we find that the focal mechanism most consistent with the data indicates thrust movement with strike, dip, and rake of 350°, 10°, and 90°, respectively, at a depth of 14 km. Given the depth uncertainties of both this event and the megathrust, this implies that the earthquake most likely took place on the subduction interface rather than on the mapped faults in the Mad River fault zone that trend 322° and dip to the northeast. The revisited intensity in the epicentral region also supports a location beneath Fickle Hill to the east of the city of Arcata, California.

California, Idaho, Nevada, Oregon, Washington↗
USGS2024

Geomorphic map of the Umatilla River corridor, Oregon

This map portrays the distribution of landforms along the Umatilla River in northeastern Oregon and covers a corridor 127 kilometers long from the confluence of the Umatilla River with the Columbia River upstream to Meacham Creek. The map encompasses the valley bottom and extends about 1 kilometer up the adjoining hillslopes. Map data are intended to support water quality and fisheries enhancement efforts pursuant to the First Foods, a resource-management approach that focuses on traditionally gathered foods including water, fish, big game, roots, and berries and calls attention to the reciprocity between people and the foods upon which humans depend. The Umatilla River drains about 6,300 square kilometers on the northwest slope of the Blue Mountains in northeast Oregon. Most of the drainage basin is underlain by Miocene basalt flows of the Columbia River Basalt Group. Younger, weakly lithified, late Miocene and early Pliocene gravel deposits of local origin (for example, McKay Formation) are mapped in a few places. Upland surfaces are mantled with windborne silt (loess) correlative with deposits elsewhere known as the Palouse Formation. Surfaces below an elevation of about 340 meters were inundated repeatedly by large Pleistocene glacial outburst floods, most emanating from glacial Lake Missoula in western Montana. In backflooded areas such as the lower Umatilla River valley, Missoula floods deposited extensive slack-water silt. Areas mapped as open water, active channel and tie channel, flood basin, valley bottom, and modified land constitute the geomorphic floodplain: the area subject to occasional inundation by the Umatilla River. Deposits and landforms within the floodplain are inset into Missoula flood deposits and hence postdate the 20–15-kilo-annum Missoula floods. Some floodplain deposits are no more than a few centuries old, as indicated by substantial erosion and deposition during the Umatilla River flood of February 2020, the largest since systematic measurements began in October 1903. Deposits and landforms of the floodplain are transient features within the longer-term incision of the Umatilla River into mid-Miocene flood basalts and younger gravel of the McKay Formation.

Oregon↗
USGS2025

Risk implications of Poisson assumptions and declustering inferred from a fully time-dependent earthquake forecast

We use the Third Uniform California Earthquake Rupture Forecast Epidemic Type Aftershock Sequence model, which is fully time-dependent in terms of including spatiotemporal clustering, to evaluate the effects of the Poisson assumption and declustering algorithms on statewide loss exceedance curves. The model is simulation based, meaning it produces synthetic catalogs that exhibit realistic behavior with respect to aftershocks and multi-fault earthquakes. A Poisson version of the model was constructed by randomizing event times, and the influence of two declustering algorithms was examined as well. We demonstrate that the probability of one-or-more loss exceedances (occurrence exceedance probability) is greater for the Poisson model because it has fewer seismically quiet time windows. The discrepancy between dollar loss estimates with a given exceedance probability is up to a factor of 32% but varies depending on the loss threshold (the x-axis value) and the forecast duration (we examined a range between 24 h and 50 years, with the discrepancy for the latter being negligible). We discuss how the one-or-more loss exceedance metric is questionable because it ignores all but the maximum loss experienced in each timeframe. An alternative metric based on total aggregate loss in each time window (aggregate exceedance probability) was therefore also examined, for which the Poisson model again implies higher risk at intermediate losses but lower risk at higher losses (because large, triggered events now contribute to total aggregate losses for the fully time-dependent model). We also argue that declustering is not a scientifically justifiable way to deal with full time dependence, in agreement with a chorus from other recent studies. It is difficult to draw generally applicable conclusions from our study, in part because application specific details will likely be important, but our results highlight how full time dependence can be reckoned with once authoritative forecast models are made available.

California↗
USGS2025

Structural controls on splay fault rupture dynamics during Cascadia megathrust earthquakes

Great subduction earthquakes ( M w ≥ 8.0) can generate devastating tsunamis by rapidly displacing the seafloor and overlying water column. These potentially tsunamigenic seafloor offsets result from coseismic fault slip and deformation beneath or within the accretionary wedge. The mechanics of these shallow rupture phenomena and their dependence on subduction zone properties remain unresolved, partly due to the sparsity of offshore observations of shallow megathrust earthquake deformation. Here, we analyze how offshore structure influences shallow rupture mechanics and slip partitioning using 3D dynamic earthquake simulations of the Cascadia subduction zone (CSZ) megathrust with and without variably dipping seaward- or landward-vergent splay faults in the wedge that sole into the megathrust. Resulting tradeoffs between splay and megathrust slip reveal structural controls on rupture partitioning, with greater splay slip leading to less shallow megathrust slip updip. Gently dipping and seaward-vergent splays host more slip than those with steeper, landward-vergent splays. To isolate the underlying mechanisms, we compare models with Andersonian and plunging principal stresses. Results suggest distinct static and dynamic processes control the dip- and vergence-dependence of splay rupture: static (mis)alignment relative to far-field tectonic loading favors slip on more optimally oriented, shallowly dipping splay faults. In contrast, dynamic stress interactions of an updip-propagating megathrust rupture front with the free surface and potential branch faults favor forward branching onto seaward-vergent splays and inhibit backward branching onto landward-vergent splays. Resulting seafloor displacements suggest splay fault structure may influence coseismic tsunami source processes, highlighting the importance of dynamically viable rupture scenarios in subduction hazard assessments.

Cascadia subduction zone↗
USGS2026

Simulation of groundwater flow to evaluate hydrogeologic controls on a PFAS plume, Coakley Landfill Superfund Site, Rockingham County, New Hampshire

Per- and polyfluoroalkyl substances (PFAS), including perfluorooctanoic acid (PFOA) and perfluorooctanesulfonic acid (PFOS), have been detected at combined concentrations above 2,000 nanograms per liter (ng/L) at groundwater seep locations near the Coakley Landfill Superfund site, in North Hampton, New Hampshire. The landfill was active from 1972 to 1985. An impermeable cap was placed on the landfill in 1998. The adjacent area to the Coakley Landfill has many water supply wells, and transport of PFAS compounds to the wells is a concern. Fracture anisotropy in the underlying bedrock aquifer complicates the understanding of PFAS transport because groundwater preferentially travels along fractures that may not align with the prevailing groundwater flow direction. In 2018, the U.S Environmental Protection Agency and the U.S. Geological Survey began an investigation of the groundwater flow from the Coakley Landfill site. This report describes the modification of a numerical groundwater-flow model for the local area around the Coakley Landfill and summarizes findings of the investigation. In addition, this report includes a brief description of PFOA and PFOS occurrence, a discussion of model construction, evaluation of model performance through calibration, and discussion of simulation results for two periods (before and after capping). Limitations are also discussed. Results show that simulated groundwater flow moves from the Coakley Landfill to the west and north. Advective transport modeling using particle tracking shows that groundwater from the landfill discharges primarily to streams to the west and north, and a small amount is transported to distal wells. Dilution of contaminants through advection and dispersion likely plays a role in whether PFAS compounds from the landfill will be detected above laboratory reporting levels at distal wells.

EarthArXiv↗
USGS2024

Joint pilot fish habitat framework

This story map will take you through the process of exploring and testing methods necessary for a higher resolution, seamless fish habitat assessment across both inland and estuarine waters through the lens of our joint pilot assessment Fish habitat assessments attempt to relate past, current, or future landscape conditions to the state of fish species occurrence, distribution, abundance, or community and habitat condition in streams, rivers, or estuaries. Previous fish habitat assessments, such as the National Fish Habitat Assessment, conducted separate and disconnected assessments for inland waters and estuaries. In this project, National Oceanic and Atmospheric Administration ( NOAA ) and U.S. Geological Survey ( USGS ) researchers created a seamless spatial framework to allow assessments that integrate influences on fish habitat from headwaters to the estuary. This effort began when the Chesapeake Bay Program Fish Habitat Action Team expressed interest in a Baywide fish habitat assessment spanning tidal salt, tidal fresh, warm non-tidal, and cold non-tidal waters. However, the complexity of the myriad of implementation details to consider when developing such an assessment necessitated the need for a tributary-specific pilot assessment. To conduct this pilot assessment, a NOAA/USGS joint partnership was formed with cooperation and support from the Chesapeake Bay Agreement and Chesapeake Bay Fish Habitat Action Team (FHAT).

Maryland↗
USGS2025

Assessment of density pattern retention of generalized data for 1:100,000-scale United States topographic maps

Cartographic generalization reduces the complexity of geographic data to produce legible, smaller-scale displays that retain essential information and logical geographic patterns. Generalization is a vital process in topographic map production. An important challenge in this process is managing and evaluating consistency across scale in the density and spatial distribution of map features such as buildings, roads, streams, water bodies, and elevation contours. Density patterns in these features reflect underlying physiographic conditions, which include factors such as bedrock geology, tectonics, climate, and landforms. Assessments of an acceptable level of change in feature density patterns are critical to ensuring the readability, usability, and accuracy of generalized maps and data. Preserving realistic density patterns across mapping scales also supports sustainable development goals in cartography, by helping to prioritize and communicate the relative reliability of geospatial data at specific scales.

Conference Paper↗
USGS2024

Geologic input databases for the 2025 Puerto Rico – U.S. Virgin Islands National Seismic Hazard Model update: Crustal faults component

The last National Seismic Hazard Model (NSHM) for Puerto Rico and the U.S. Virgin Islands (PRVI) was published in 2003. In advance of the 2025 PRVI NSHM update, we created three geologic input databases to summarize new onshore and offshore fault source information in the northern Caribbean region between 62°–70° W and 16°–21° N. These databases, of fault sections, fault‐zone polygons, and geologic estimates of fault activity (fault‐slip rate and earthquake recurrence intervals) at specific sites, document updates to fault parameters used in prior seismic hazard models in PRVI. Fault sources were reviewed from published studies since 2003, which document substantial changes to the understanding of fault location, geometry, or activity. New fault section sources were added for features that meet the criteria of (1) length ≥7 km, (2) unequivocal evidence of recurrent tectonic Quaternary activity, and (3) documentation that is publicly available in a peer‐reviewed source. In addition, we revised several broad areal sources, such as the Mona and Anegada extensional zones. The 2003 model included three fault sections and two fault‐zone polygons (areal sources). These databases include 35 fault sections, 6 fault‐zone polygons, and 51 earthquake geology sites. To characterize fault activity rates, slip‐rate bins were assigned based on landscape expression and paleoseismic trench observations for faults without published slip‐rate sites. Additional fault sources were evaluated but not included in these databases due to a lack of published information about fault location, geometry, or recurrent Quaternary activity. The PRVI NSHM 2025 geologic input databases describe crustal faulting; the geometries and coupling of Puerto Rico subduction zone and Muertos Trough models are considered in a separate database. Updates to the fault sections, fault‐zone polygons, and earthquake geology databases can help inform the location and recurrence rate of damaging earthquakes in the PRVI NSHM implementation.

Puerto Rico, U.S. Virgin Islands↗
USGS2026

Cnidarian–algal partnerships structure bacterial communities during strobilation in Cassiopea xamachana

Cnidarian–algal (Symbiodiniaceae) symbioses rely on complex interactions among the cnidarian host, algal symbionts, and associated bacterial communities. In the upside-down jellyfish Cassiopea xamachana , the polyp-to-medusa transition (strobilation) requires the establishment of symbiosis with Symbiodiniaceae algal partners, yet bacterial community dynamics during this developmental process remain unknown. Here, we experimentally induced symbiosis in aposymbiotic polyps using four algal treatments: xenic Symbiodinium microadriaticum (native symbiont), xenic Breviolum minutum , antibiotic-treated B. minutum , and a photosynthetically impaired B. minutum mutant. We combined 16S rRNA gene sequencing with measurements of photosynthetic efficiency, asexual budding, and algal surface N-glycan profiles to characterize holobiont assembly during symbiosis onset and strobilation. Algal treatment structured bacterial communities in both algal cultures and polyp tissues. Our analyses identified a set of amplicon sequence variants that consistently distinguished strobilating polyps from non-strobilating aposymbiotic and mutant polyps, in addition to potential bacterial biomarkers associated with successful metamorphosis. Strobilation was associated with the enrichment of bacterial communities putatively involved in sulfur and nitrogen cycling, whereas non-strobilating aposymbiotic and mutant polyps were characterized by opportunistic bacteria and increased community variability. Together, these results reveal coordinated changes in algal physiology, surface glycan profiles, and bacterial community structure associated with successful strobilation in C. xamachana and support a model in which tripartite host–alga–bacteria interactions influence cnidarian life stage transitions.

ISME Communications↗
USGS2024

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

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

Michigan, Wisconsin↗
USGS2025

Smectite-rich horizons in Inceptisols trigger shallow landslides in tropical granitic terranes

Puerto Rico was affected by >70,000 landslides in the wake of 2017 Hurricane Maria, and landslide prevalence was especially high in the Utuado region in the Cordillera Central. Landslide density was highest where soil parent material is granodiorite; landslide slip surfaces tended to be shallow (<60 cm), and often were mobilized rapidly and with long runout distances. This study combines field observations with soil mineralogy (bulk and clay fractions), soil geochemistry (bulk fraction), and soil strength as determined by field cone penetrometer testing (CPT) to test the hypothesis that clay-rich subsoil horizons function as slip planes when water-logged. Soil pits were excavated to depths of ∼200 cm in Ultisols on an undulating plateau and to ∼100 cm in Inceptisols on steep slopes (36-43 o ) that flank the plateau and cone penetrometer tests (CPT) were done within 2 m of the soil pit. Six pits were located adjacent to scarps from previous landslides, enabling analysis of soil profiles downward through extrapolated slip surfaces. Results from X-ray diffraction (XRD), X-ray fluorescence (XRF) and thermogravimetric analysis (TGA) indicate that soils are heterogeneous, often with subsoil horizons enriched in clay minerals and immobile elements (Al, Fe, Ti). Inceptisols on steep slopes often contain smectite-rich horizons at 30–60 cm depth that appear to function as slip surfaces; in other Inceptisols, such horizons are not present and landslide susceptibility is potentially lower. In Ultisols, soil mineralogy is dominated by kaolinite with minor halloysite, and increased kaolinite content at soil depths ≥80 cm at some sites suggests potential slip surfaces enhancing probability of landslides. The origin of clay-rich horizons appears to be (1) fractures in granodiorite that facilitate water flow and leaching, accelerating mineral dissolution during early weathering stages, and (2) smectite-rich buried soils under permeable colluvium likely deposited by a prior mass wasting event. Where clay-rich layers occur beneath more-permeable horizons, rapid infiltration then absorption of water in clay-rich subsoil horizons causes decreased shear strength and increased landslide susceptibility.

Puerto Rico↗
USGS2000

Water quality of selected springs and public-supply wells, Pine Ridge Indian Reservation, South Dakota, 1992-97

This report presents results of a water-quality study for the Pine Ridge Indian Reservation, South Dakota. The study was a cooperative effort between the U.S. Geological Survey and the Water Resources Department of the Oglala Sioux Tribe. Discharge and water-quality data were collected during 1992-97 for 14 contact springs located in the northwestern part of the Reservation. Data were collected to evaluate potential alternative sources of water supply for the village of Red Shirt, which currently obtains water of marginal quality from a well completed in the Inyan Kara aquifer. During 1995-97, water-quality data also were collected for 44 public-supply wells that serve about one-half of the Reservation's population. Quality-assurance sampling was used to evaluate the precision and accuracy of environmental samples. Ten of the springs sampled contact the White River Group, and four contact the Pierre Shale. Springs contacting the White River Group range from calcium bicarbonate to sodium bicarbonate water types. Two springs contacting the Pierre Shale have water types similar to this; however, sulfate is the dominant anion for the other two springs. In general, springs contacting the White River Group are shown to have better potential as alternative sources of water supply for the village of Red Shirt than springs contacting the Pierre Shale. Nine of the springs with better water quality were sampled repeatedly; however, only minor variability in water quality was identified. Six of these nine springs, of which five contact the White River Group, probably have the best potential for use as water supplies. Discharge from any of these six springs probably would provide adequate water supply for Red Shirt during most periods, based on a limited number of discharge measurements collected. Concentrations of lead exceeded the U.S. Environmental Protection Agency (USEPA) action level of 15 µg/L for three of these six springs. Five of these six springs also had arsenic concentrations that exceeded 10 µg/L, which could be problematic if the current maximum contaminant level (MCL) is lowered. Blending of water from one or more springs with water from the existing Inyan Kara well may be an option to address concerns regarding both quantity and quality of existing and potential sources. All nine springs that were sampled for indicator bacteria had positive detections on one or more occasions during presumptive tests. Although USEPA standards for bacteria apply only to public-water supplies, local residents using spring water for domestic purposes need to be aware of the potential health risks associated with consuming untreated water. One spring contacting the White River Group and two springs contacting the Pierre Shale exceeded 15 pCi/L for gross alpha; these values do not necessarily constitute exceedances of the MCL, which excludes radioactivity contributed by uranium and radon. Additional sampling using different analysis techniques would be needed to conclusively determine if any samples exceeded this MCL. Nine springs were sampled for selected pesticides and tritium. The pesticides atrazine, carbaryl, and 2,4-D were not detected in any of the samples. The nine springs were analyzed for tritium in order to generally assess the age of the water and to determine if concentrations exceeded the MCL established for gross beta-particle activity. Tritium results indicated two springs are composed primarily of water recharged prior to atmospheric testing of nuclear bombs and two other springs have a relatively large percentage of test-era water. The remaining five springs had tritium values that indicated some percentage of test-era water; however, additional sampling would be needed to determine whether water is predominantly pre- or post-bomb age. Of the 44 public-supply wells sampled, 42 are completed in the Arikaree aquifer, one is completed in an alluvial aquifer, and one is completed in the Inyan Kara aquifer. Water from the alluvial well is a sodium bicarbonate water type, water from Arikaree aquifer ranges from calcium bicarbonate to sodium bicarbonate types, and water from the Inyan Kara well is a calcium sulfate bicarbonate type. Of the 44 wells sampled, 28 (64 percent) tested positive for indicator bacteria in presumptive tests. Because these were single samples that generally were collected upstream from chemical treatment feeders, positive detections do not necessarily constitute exceedances of drinking-water standards. A single sample from an Arikaree well exceeded the MCL for arsenic of 50 µg/L. Arsenic exceeded 10 µg/L for six additional Arikaree wells and for the alluvial well and the Inyan Kara well, which could be problematic if the current MCL is lowered. The alluvial well also exceeded the secondary maximum contaminant level (SMCL) for dissolved solids, which is non-enforceable, and the action level for lead. The Inyan Kara well exceeded the SMCL's for iron and for manganese and the MCL of 5 pCi/L for radium-226 and 228 combined. Several Arikaree wells exceeded SMCL's for either pH, sulfate, dissolved solids, iron, or manganese. One Arikaree well exceeded the MCL of 4.0 mg/L for fluoride and another exceeded the MCL of 10 mg/L for nitrite plus nitrate. Ten Arikaree wells equalled or exceeded 15 pCi/L for gross alpha; however, these values do not necessarily constitute exceedances of the MCL, which excludes radioactivity contributed by uranium and radon. Additional sampling using different analysis techniques would be needed to conclusively determine if any samples exceeded this MCL. Eight wells, all from the Arikaree aquifer, equalled or exceeded the proposed MCL of 20 µg/L for uranium and 33 wells (75 percent) equalled or exceeded one-half of the proposed MCL. Although this standard has only been proposed, additional information regarding the extent of elevated uranium concentrations in the Arikaree aquifer, and the geochemical processes involved, may be beneficial. It was determined from analyses of uranium isotope data for five wells that the source of elevated uranium concentrations is naturally occurring, rather than anthropogenic.

South Dakota↗
USGS2025

The 3D National Topography Model Call for Action—Part 2: The Next Generation 3D Elevation Program

The three-dimensional (3D) National Topography Model initiative to integrate elevation and hydrography data includes the next generation of hydrography data from the 3D Hydrography Program and the next generation of elevation data from the 3D Elevation Program (3DEP). The first-ever collection of light detection and ranging (lidar) data for the nation (IfSAR for Alaska) provides a critical baseline reference, and the addition of multiple repeat elevation mapping projects as part of the next generation of 3DEP would substantially expand analysis capabilities. As the U.S. Geological Survey (USGS) is closing in on our goal of complete coverage with 3DEP data available or in progress for 98.3 percent of the Nation at the end of fiscal year 2024, the USGS is already transitioning to the next generation of 3DEP. Based on the 3D Nation Study results and input from a broad range of stakeholders, the USGS National Geospatial Program has finalized a new design for 3DEP that provides increased lidar quality levels and refresh rates. The new program is designed with more flexibility to meet changing user needs and take advantage of improvements in mapping technologies. The program will aim to expand the level of interagency coordination for topobathymetric lidar acquisition for inland rivers. The next generation of 3DEP will also aim to emphasize research, including advancing program design, products, and services and engaging and leveraging the evolving 3D industry. Research goals also include becoming more flexible in meeting user needs and taking advantage of evolving remote-sensing technologies. The program also plans to move from focusing on producing standard products to producing a concept of a 3D Nation Ecosystem with a variety of inputs, products, and services.

Circular↗
USGS2026

Decadal shifts in groundwater age detected by environmental tracers across California, USA

Groundwater age offers important insight into recharge, storage, and contamination risk. Although models predict age changes can be driven by pumping and climate variability, direct observational evidence remains limited. Here, we analyzed paired environmental tracer suites (tritium, carbon-14, and tritiogenic helium-3) collected a decade apart from 268 wells across California to assess the prevalence of groundwater age transience. Travel-time distribution models and statistical tests indicated age transience at 29% of sites, occurring most often in agricultural regions, such as the San Joaquin Valley and Southern Coast Ranges, where large carbon-14 changes coincided with substantial nitrate and chloride shifts. Sites with tritiogenic helium-3 data showed more frequent age transience, underscoring the value of multi-tracer data sets. These results provide the first regional evidence of widespread groundwater age change and a method for detecting changing water balances with implications for groundwater sustainability and water quality.

California↗
USGS2026

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

Ungulate migrations of the Western United States, volume 5

Many ungulates migrate between distinct summer and winter ranges to take advantage of spatially and temporally variable food sources and avoid threats such as predators and deep snow. In 2018, the U.S. Department of the Interior established Secretarial Order 3362, which provided Federal support to expand existing research efforts to study ungulate populations and conserve their migrations by enhancing habitat quality for ungulates across the Western United States. In response to the order, the U.S. Geological Survey (USGS) created the Corridor Mapping Team, which is a collaboration among 11 State agencies, regional and Federal partners, and an expanding number of Tribal wildlife agencies. Together, the Corridor Mapping Team maps ungulate migrations throughout the Western United States and publishes them in the USGS “Ungulate Migrations of the Western United States” report series. This report details migrations and seasonal ranges from 36 additional herds and includes 2 herd updates detailed in previous reports. The Corridor Mapping Team has mapped the migrations and seasonal ranges of 218 unique herds for the report series, including this report. The report series serves as a map-based inventory of the ungulate migrations across the Western United States for biologists, managers, policymakers, and conservation practitioners. Building on the previous report volumes in the series, volume 5 additionally describes some of the local and national initiatives that are incorporating the products, tools, and information from this growing USGS report series.

Navajo Nation↗
USGS2025

Evidence for offset of Cretaceous plutons by the Tintina fault in eastern Alaska: Implications for regional metallogeny

Cretaceous magmatism in eastern interior Alaska is voluminous, but temporally and spatially diverse – suggestive of varying sources and drivers. More than 150 new U-Pb zircon and more than 500 geochemical analyses of Cretaceous plutonic units allow for the grouping of distinct plutonic suites. Magmatism was continuous from 120-66 Ma but can be grouped into temporally distinct pulses from ca. 115-100 Ma, 100-90 Ma, and 75-66 Ma. Geochemical diversity occurs during each pulse, further distinguishing multiple suites. Diverse metallogenic epochs are strongly correlated to pluton chemistry. Mineralization is largely absent prior to 108 Ma. From 108-100 Ma, plutonism is coeval with sparse, but notable Au-quartz veins with variable Bi, As, W, and Mo. From 100-90 Ma, intrusion-related mineralization zones from Au-Cu(-Bi) and U-Th in the northwest to central Au- Bi-As-Te(-W), and Mo-W to the southeast. Porphyry style Cu-Mo(-Au) occurrences occur with the latest Cretaceous plutons emplaced from 75-66 Ma . Restoration of ~450 km of dextral movement on the Tintina fault and comparison of metallogenic and geochemical characteristics of Alaska plutons suggest 100-90 Ma plutons may be the continuation of the metallogenically significant Tombstone, Mayo, and Tungsten suites from the Yukon.

Alaska↗