USGS Science⌕ Search

SEARCH · USGS Science

Results for “Agricultural Water Management”

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.

754 records · Page 35Linked to original sources
USGS2026

Increasing artificial light at night enhances salmonid predator foraging efficiency in an urbanized lake

Artificial light at night (ALAN) poses a threat to ecosystems globally. It includes both direct and indirect light, or skyglow, which occurs when ALAN scatters in the atmosphere, extending beyond its original source. We analyzed ALAN trends in Lake Washington, WA, from 2014 to 2023 using Visible Infrared Imaging Radiometer Suite (VIIRS) nighttime light measurements, evaluated the relationship between in situ and satellite measurements, and modeled predator search volumes for a juvenile salmon predator, cutthroat trout ( Oncorhynchus clarki ), as a function of ambient light and turbidity conditions. Open water regions experienced significant increases in ALAN, while nearshore areas primarily showed no or negative trends, revealing the role of skyglow in shaping open water light environments. Using a visual foraging model for cutthroat trout, we found that juvenile salmon at the shallow southern pelagic site experienced light 28 times brighter, resulting in a 168% greater predation vulnerability than those at the northern site. In the Ship Canal, a narrow corridor for outmigrating salmon, predator search volumes were 249% higher than at the southern site. These contrasts in predation vulnerability demonstrate how local conditions influence predator–prey dynamics and provide critical insight for targeting mitigation of both nearshore and distant light sources.

Washington↗
USGS2025

Properties of solid bitumen formed during hydrous, anhydrous, and brine pyrolysis of oil shale: Implications for solid bitumen reflectance in source-rock reservoirs

Solid bitumen reflectance (BR o in %) is widely used as a thermal maturity proxy in source-rock reservoirs, yet solid bitumen texture and physical attributes may be affected by varying environmental constraints during its formation, e.g., water saturation, mineral catalysis, or salinity. Herein we investigated the development of solid bitumen properties during artificial maturation using three diverse (lacustrine to marine) oil shale samples containing abundant amorphous H-rich sedimentary organic matter (bituminite). The oil shales were treated via pyrolysis (320 °C, 72 h) using hydrous, anhydrous, and brine conditions, causing the development of a newly formed solid bitumen in the experiment residues from the thermal conversion of bituminite. The properties of the newly formed solid bitumen then were evaluated via geochemical screening tests, optical and electron microscopy, and infrared spectroscopy. Experimental residues also were treated via solvent extraction, allowing characterization of the effects of extraction on solid bitumen. Results showed that bituminite with higher resin and asphaltene components evolved to a solid bitumen with higher reflectance (as a percentage of the original value) when exposed to the same heat treatment, relative to bituminite with higher saturate and aromatic components. Aromatization of solid bitumen also was generally more pronounced in the presence of deionized water relative to anhydrous conditions, supporting prior observations. These results suggest the compositions of primary sedimentary organic matter, and the local concentration of water, affect the development of solid bitumen during thermal advance, potentially explaining the origin of multiple solid bitumen populations with varying reflectance in source-rock reservoir samples or in a single microscope field. Experiments using brine were inconclusive regarding enhanced/decreased aromatization, which could be attributed to the salinities of the brines used in the tests. Extraction of residues caused a consistent increase in BR o values, suggesting that migrating or expelled oils could cause an increase in BR o via natural solvation and absent an increase in temperature. This work provides insights into the development and evolution of BR o in source-rock reservoirs as a function of the composition of the original bituminite and changing environmental conditions, with potentially broad application for petroleum prospecting and resource estimation.

Colorado↗
USGS2026

Estimating paleotemperature using stable isotopes of soil-formed phyllosilicates from paleosols: A review

Fossilized soils, or paleosols, contain soil-formed phyllosilicates whose stable isotopic compositions may be used to calculate paleotemperature and thus reconstruct ancient terrestrial environments. Though paleosols are common in the geologic record, the use of phyllosilicates as paleotemperature proxies is limited in the literature owing to difficulties with selecting optimal paleosols, isolation from non-clay minerals and organic materials, mixtures of phyllosilicates in natural samples, wide variations of chemical compositions for phyllosilicates, and limited to undefined equilibrium fractionation factors between phyllosilicates-water. Here, we address these challenges by examining and comparing methods used for sample selection, mineral isolation, pretreatments, mineral identification, conventional and developing methods for oxygen and hydrogen isotopic analyses, and determination of phyllosilicate-water equilibrium fractionation factors, concluding with recommendations for best approaches for paleotemperature estimation. Additionally, we discuss how to identify and avoid detrital phyllosilicates, the impacts of diagenesis, comparison of stable isotope and non-isotope paleosol paleotemperature proxies, and challenges and opportunities for broadly using paleosols as paleoclimate archives. With ongoing efforts to refine this multi-faceted paleotemperature approach, the stable isotope geochemistry of soil-formed phyllosilicates continues to be an invaluable proxy system, enhancing our understanding of terrestrial paleoenvironments and paleoclimate.

Kansas↗
USGS2026

Inference of pattern-based geological CO2 sequestration and oil recovery potential in a commingled main pay and residual oil zone CO2-EOR flood

Several detailed studies have shown that residual oil zones (ROZs) can present significant resources for additional hydrocarbon recovery as well as subsurface carbon dioxide (CO 2 ) sequestration via enhanced oil recovery by injecting CO 2 (CO 2 -EOR). Field development strategies included new wells drilled dedicated to main pay zones (MPZ) and ROZs, or existing wells in MPZs deepened to ROZs for commingled injection-production using different well patterns. The latter presented a challenge when discerning the injection and production from each of the zones, and for subsequent quantification of CO 2 sequestration and EOR potential from different patterns and from the field. In this paper, an innovative method for analyzing commingled injections and productions from MPZs and ROZs, with application to pattern-based data from four staggered line drive patterns in Wasson Field's Denver Unit, Texas, USA, was developed. Decline curve and ratio-trend methods were used as means of history-matching and forecasting. Cumulative production-time and cumulative production-rate data for oil, gas, and water, as well as water-oil ratio (WOR) and gas-oil ratio (GOR), were analyzed along with injection data for time intervals covering major injection events in MPZ, or MPZ and ROZ combined. A combined analysis enabled inference of allocation of fluids into different zones during WAG (water alternating gas) injection and thereby estimation of CO 2 storage, utilization, and retention in different zones as a function of total injection. Results show that ROZs generally present higher CO 2 sequestration potential compared to MPZs, and a comparable incremental oil recovery factor of ∼20%, on average. Results based on ratio analysis further show that while the WOR trend of the pattern production is mostly dominated and controlled by ROZ, GOR is controlled by both intervals. Although the method relying on decline curves and the approach used in zonal fluid allocations are subject to their limitations, this study presents a practical and innovative well-pattern-based method to infer and forecast CO 2 sequestration and oil recovery quantities and fluid ratios from MPZs and ROZs in commingled operations and highlight the added potential offered by ROZs.

Texas↗
USGS2026

Microplastics undergo fragmentation during pressurized membrane filtration

Low-micrometer microplastics (<10 μm) have been detected in drinking water, driving growing interest in using pressure-driven membrane filtration to remove these particles and ensure drinking water safety. However, little attention is paid to the fate of microplastics concentrated in the filtration byproduct. In this study, using well-defined polystyrene (PS) and poly(methyl methacrylate) microspheres as model particles, we observed microplastic fragmentation during nanofiltration and the subsequent release of smaller fragments. After operating for 3 h at 20 bar, 67.9 ± 8.0% of the PS spheres in the concentrate, based on particles counted in selected fields of view, were transformed into fractured particles. The characteristic Raman band signals of microplastics were significantly weakened after fragmentation, leading to detection challenges. To address this, a data processing algorithm was developed to identify PS fragments as small as 1 μm. Preliminary statistical analysis reveals that within the tested pressure range, operation time has a stronger influence on fragmentation than pressure magnitude alone and that fragmentation is further governed by the intrinsic mechanical properties of the tested model polymers. This work provides direct evidence that pressure-driven membrane processes induce microplastic fragmentation and highlights the environmental risks associated with the discharge of fragmented microplastics into the concentrate.

Environmental Science & Technology Letters↗
USGS2025

Groundwater budget for the surficial aquifer surrounding Lake Nokomis, Minneapolis, Minnesota

During prolonged periods of above-average precipitation, rising groundwater levels have the potential to cause damage to and interfere with underground infrastructure and building foundations. To understand the relations between precipitation and groundwater in the vicinity of Lake Nokomis, the U.S. Geological Survey, in collaboration with the University of Minnesota, quantified five components of the groundwater budget: groundwater recharge, change in surficial aquifer storage, surficial aquifer groundwater discharge to Lake Nokomis, groundwater evapotranspiration, and groundwater discharge to underlying bedrock aquifers. Field data, geologic records, and empirical calculation methods were used to quantify groundwater budget components for April 2023 through April 2024. Lake water budget data indicate that Lake Nokomis is a flowthrough system during periods with no outflow through the weir, with groundwater inputs equal to outputs. Roughly 40 percent of precipitation that fell in the study area was added to the surficial aquifer as recharge. Uncertainty in the vertical hydraulic conductivity resulted in wide-ranging estimates (spanning three orders of magnitude) of water discharging from the surficial aquifer to the underlying bedrock aquifer. Drought conditions persisted for the duration of this study and were not representative of the conditions that motivated this study. This study is a start towards understanding relations between precipitation, Lake Nokomis levels, and groundwater levels that could affect local underground infrastructure.

Minnesota↗
USGS2025

Examining the compositional selectivity of hydrocarbon oxidation products using liquid–liquid extraction and solid-phase extraction techniques

The effect of extraction methods on detecting hydrocarbon oxidation products (HOPs) in groundwater remains unclear. HOPs are polar, water-soluble byproducts of petroleum biodegradation. Our previous work showed that liquid–liquid extraction (LLE), a method commonly used in regulatory monitoring, has a significantly lower extraction efficiency for HOPs compared to solid-phase extraction (SPE). In this study, we evaluate the analytical limitations and compositional selectivity of LLE and SPE using groundwater samples from the Bemidji, MN, crude oil spill site. Optical properties were characterized using excitation–emission matrix spectroscopy (EEMs), and a three-component PARAFAC model was validated, showing consistent trends across both extracts and whole water samples. Ultrahigh-resolution mass spectrometry (UHR-MS) revealed that LLE selectively recovered aliphatic-like compounds but underrepresented more polar oxygenated HOPs. In contrast, SPE methods were more effective at isolating highly oxidized compound classes. These differences were consistent across a gradient of contamination. Overall, the LLE was less precise and less representative of polar HOPs, introducing bias in the characterization of HOPs. This study is the first to quantitatively demonstrate the compositional selectivity and analytical bias of LLE versus SPE for HOPs using combined EEM-PARAFAC and UHR-MS techniques, with implications for long-term monitoring and site assessment protocols.

Environmental Science and Technology↗
USGS2026

Evaluation of stream capture related to groundwater pumping, middle Humboldt River Basin, Nevada

Historical, future, and potential stream capture from groundwater pumping in the middle Humboldt River Basin (MHRB), Nevada, is estimated using a calibrated numerical groundwater flow model. The model was developed to estimate (1) stream capture, which is the change in flux between the groundwater system and the Humboldt River and tributaries, and (2) change in streamflow, which is the change in streamflow estimated for the Imlay gage on the Humboldt River (U.S. Geological Survey streamgage 10333000). Historical stream capture for water years (WYs) 1961–2015 is estimated using recorded and estimated groundwater pumping during that period. Future (predictive) stream capture was based on historical stresses (WYs 1961–2015) using a scenario that simulated non-mine pumping from WY 2015 at a uniform rate for 100 years into the future. Potential stream capture throughout the middle Humboldt River Basin from groundwater pumping during varying durations of time are presented in a series of capture maps. Maps also are presented that show the potential to capture from groundwater evapotranspiration, as well as the storage changes for pumping duration of 100 years. Estimates of historical stream capture from the mainstem Humboldt River during the early 1960s are less than 400 acre-feet per year (acre-ft/yr) when groundwater withdrawals and pumping rates were relatively small compared to more recent times. In the late 1980s and early 1990s, groundwater withdrawals increased and estimated historical stream capture also increased from about 4,000 acre-ft/yr in the late 1980s and early 1990s to as much as 18,800 acre-feet (acre-ft) in WY 1998. In WY 2015, estimated historical stream capture declined to about 13,000 acre-ft because of decreasing groundwater withdrawals and lower streamflow during the drought of WYs 2012–15, resulting in less stream water available for capture. Stream capture was estimated for 100 years into the future based on WY 2015 non-mine pumping rates and mine-dewatering activity through WY 2015. Stream capture is forecast to increase to about 23,000 acre-ft/yr, and streamflow in the Humboldt River could decrease by as much as 19,000 acre-ft/yr. Pumping for mine-dewatering and the associated discharge of that water affects streamflow in the Humboldt River at Imlay, Nevada (U.S. Geological Survey streamgage 10333000). Historically, from WYs 1991 to 2015, streamflow was greater at Imlay gage during active mine-dewatering from mine-water discharge operations and increased by as much as 105,000 acre-ft in WY 1998. The increase was attributed mostly to the discharge of groundwater from mine-related dewatering operations directly into the mainstem Humboldt River or its tributaries, with some of this increase associated with return flows from discharge to rapid infiltration basins. Results indicate that streamflow at Imlay gage is expected to decrease by as much as 1,600 acre-ft/yr 30 years after mine-related pumping and discharge are discontinued. The streamflow reductions at the Imlay gage are expected to then decrease to around 500 acre-ft/yr, 100 years after mine-related pumping and discharge are discontinued. Potential capture maps were produced for pumping durations of 10, 25, 50, and 100 years. Capture map results indicate that areas of greater potential stream capture occur adjacent to the Humboldt River and for upstream tributaries areas north of the Humboldt River.

Nevada↗
USGS2025

Ordovician stratigraphy, structure, and karst of the Falling Spring Valley, Alleghany County, Virginia, USA

This one-day trip highlights new findings on a preliminary bedrock geologic map that shows results from ongoing geologic mapping in the Falling Spring Valley of Alleghany County, Virginia, USA, which is the southern end of the larger Warm Springs Valley, an elongated anticlinal valley rimmed by Ordovician and Silurian siliciclastic rocks, and which is famous for its thermal springs. This mapping includes stratigraphic, structural, and karst field and lab research focused on the Ordovician strata exposed in the area, the oldest of which is the dolomitic upper part of the Beekmantown Formation (Lower Ordovician, Darriwilian), and the youngest of which is the Juniata Formation (Upper Ordovician, Katian), a sequence of siliciclastic redbeds. Warm Springs Valley is the location of the only known caves in the eastern United States—three at present—with thermal waters flowing in some of their passages. Stops on the trip will highlight key details from mapping efforts, primarily within the structurally deformed Ordovician carbonate sequence that is exposed in the core and limbs of the anticline, as well as the associated karst features that are developed in those carbonate rocks, including results of recent dye traces and water temperature monitoring that have improved our understanding of the karst hydrogeologic systems developed in these strata.

Virginia↗
USGS2025

Wetlands, groundwater and seasonality influence the spatial distribution of stream chemistry in a low-relief catchment

Evaluating stream water chemistry patterns provides insight into catchment ecosystem and hydrologic processes. Spatially distributed patterns and controls of stream solutes are well-established for high-relief catchments where solute flow paths align with surface topography. However, the controls on solute patterns are poorly constrained for low-relief catchments where hydrogeologic heterogeneities and river corridor features, like wetlands, may influence water and solute transport. Here, we provide a data set of solute patterns from 58 synoptic surveys across 28 sites and over 32 months in a low-relief wetland-rich catchment to determine the major surface and subsurface controls along with wetland influence across the catchment. In this low-relief catchment, the expected wetland storage, processing, and transport of solutes is only apparent in solute patterns of the smallest subcatchments. Meanwhile, downstream seasonal and wetland influence on observed chemistry can be masked by large groundwater contributions to the main stream channel. These findings highlight the importance of incorporating variable groundwater contributions into catchment-scale studies for low-relief catchments, and that understanding the overall influence of wetlands on stream chemistry requires sampling across various spatial and temporal scales. Therefore, in low-relief wetland-rich catchments, given the mosaic of above and below ground controls on stream solutes, modeling efforts may need to include both surface and subsurface hydrological data and processes.

Michigan↗
USGS2025

Detecting hidden sedimentary geothermal systems in the Upper Colorado River Basin

Geothermal resources exist in sedimentary rock where circulation of water for efficient extraction or storage of heat is possible. Except in rare instances where hot water is expressed at the land surface, sedimentary geothermal resources are hidden, so the identification of these systems is optimally accomplished using predictive subsurface modeling. An integrated approach using detailed paleogeographic interpretations, subsurface geologic mapping, and numerical modeling has produced regional geologic and temperature models for the Upper Colorado River Basin, a large watershed in central North America that contains many sedimentary basins. These models identify areas of hidden sedimentary geothermal resource potential in low temperature (<90°C), moderate temperature (90–150°C), and high temperature (>150°C) fairways across the study area. These models incorporate maps of key horizons in outcrop and the subsurface to create a robust structural framework that can be used to target favorable geology for natural or engineered permeability. This framework is populated with lithologies derived from detailed palaeogeographical maps and over 40,000 bottom hole temperature (BHT) values were used to create a calibrated three-dimensional (3D) temperature model across the region. The resulting maps serve as a regional sedimentary geothermal play fairway screening tool for evaluating different grades of sedimentary geothermal resources and for identifying areas of interest where more detailed, prospect-scale studies can be undertaken.

Arizona, Colorado, New Mexico, Utah, Wyoming↗
USGS2026

Assessment of groundwater quantity and quality contributions to Lake Huron

Lake Huron, one of the five Great Lakes, borders the United States and Canada, with Michigan as the only U.S. State on its shoreline. Like other freshwater lakes, it faces water-quality challenges from nutrients and chemicals applied across its drainage basin. Although past studies focused on surface-water sources, groundwater contributions remain less understood. To address this gap, the U.S. Geological Survey, as part of the Cooperative Science and Monitoring Initiative, classified drainage basins to Lake Huron into eight hydrogeologic zones based on bedrock rock type and glacial sediment transmissivity. Utilizing existing data and empirical field data, we quantified groundwater discharge and identified areas of concern for loading of chloride and nitrate to Lake Huron. Groundwater contributions, including indirect and shoreline discharge, ranged from 5.8 to 11.5 inches annually, totaling 1.9 cubic miles and 0.09 cubic mile, respectively. Hydrogeologic zones with higher glacial sediment transmissivity yielded greater indirect groundwater discharge. Chloride levels above the U.S. Environmental Protection Agency’s 250-mg/L recommendation were mainly in the Saginaw lowlands, whereas nitrate above the 10-mg/L standard was rare—found in only 11 wells. Together, the analysis of where groundwater discharge is occurring in the Lake Huron Basin and the identification of areas with potential groundwater-quality concerns can help prioritize areas that are critical to protecting the long-term health of Lake Huron.

Michigan↗
USGS2025

Linking permafrost to the abundance, biomass, and energy density of fish in Arctic headwater streams

Permafrost thaw alters groundwater flow, river hydrology, stream-catchment interactions, and the availability of carbon and nutrients in headwater streams. The impact of permafrost on watershed hydrology and biogeochemistry of headwater streams has been demonstrated, but there is little understanding of how permafrost influences fish in these ecosystems. We examined relations among permafrost characteristics, the resulting changes in water temperature, stream hydrology (e.g., discharge flashiness), and macroinvertebrates, with the abundance, biomass, and energy density of juvenile Dolly Varden ( Salvelinus malma ) and Arctic Grayling ( Thymallus arcticus ) across 10 headwater streams in northwestern Alaska. Macroinvertebrate density was driven by concentrations of dissolved carbon and nutrients supporting stream food webs. Dolly Varden abundance was primarily related to water temperature with fewer fish in warmer streams, whereas Dolly Varden energy density decreased with the flashiness of the headwater streams. Dolly Varden biomass was related to both temperature and bottom-up food web effects. The energy density of Arctic Grayling decreased with warmer temperatures and discharge flashiness. These relations demonstrate the importance of terrestrial–aquatic connections in permafrost landscapes and indicate the complexity of landscape effects on fish. Because permafrost thaw is one of the most impactful changes occurring as the Arctic warms, an improved understanding of how stream temperature, hydrology, and bottom-up food web processes influence fish populations can aid forecasting of future conditions across the Arctic.

Alaska↗
USGS2024

Hydrogeologic framework and extent of saltwater intrusion in Kings, Queens, and Nassau Counties, Long Island, New York

In 2016, the U.S. Geological Survey began a multiyear cooperative study with the New York State Department of Environmental Conservation to evaluate the sustainability of Long Island’s sole-source aquifer system through hydrogeologic mapping, compilation of groundwater chloride concentrations, and groundwater flow modeling. In the initial phase of the islandwide study, the hydrogeologic framework and extent of saltwater intrusion in aquifers in Kings, Queens, and Nassau Counties on western Long Island, N.Y., were investigated. The aquifer system underlying western Long Island has been under stress from pumping of public, irrigation (golf course), and industrial supply wells. Saltwater intrusion has occurred from surrounding embayments (East River, Long Island Sound, Jamaica Bay, and the Atlantic Ocean) due to pumping. Eighteen boreholes were drilled and cores taken during 2019–21 to collect hydrogeologic, geochemical, and geophysical data to delineate the complex subsurface hydrogeology and extent of saltwater intrusion within the study area. Evaluation of the new cores, reexamination of legacy core descriptions, and analysis of borehole geophysical logs was used to refine the previously published hydrogeologic framework of Pleistocene and Cretaceous unconsolidated sediments in the area, including delineation of a previously undefined hydrogeologic unit between the Magothy aquifer and the Raritan confining unit, herein named the “upper Raritan aquifer.” The upper Raritan aquifer was first recognized in southeastern Nassau County from an analysis of about 50 closely spaced boreholes with high-resolution core descriptions and gamma-ray (gamma) logs. Further analysis of borehole logs across the study area indicated that the upper Raritan aquifer was also present in Kings and Queens Counties. Nuclear magnetic resonance (NMR) logging was used for the first time on Long Island to provide estimates of the hydraulic properties of the major aquifer and confining units. Unlike other geophysical logs that record responses to the rock matrix and fluid properties and are strongly dependent on mineralogy, NMR logs record responses to the presence of hydrogen protons in the formation fluid to determine water fraction and pore-size distribution. NMR log analysis provided estimates of the clay-bound, capillary-bound, and mobile water fractions and hydraulic conductivity of aquifers and confining units penetrated by five wells in Nassau County. Pumpage for public-supply and industrial wells on Long Island began in the 1870s with small, localized suppliers of populated areas in Kings and Queens Counties. By 1904–16, pumpage for public water supply in Kings County averaged 21 million gallons per day, and averaged 37 million gallons per day in Queens County, mostly from the upper glacial aquifer. Saltwater intrusion was reported as early as the beginning of the 20th century and included the upper glacial-Jameco-Magothy and Lloyd-North Shore aquifer systems. By 1936, pumping in central Kings County created a major cone of depression in the water table extending to the south shore of much of Kings County and into southwestern Queens County. Saltwater intrusion has caused the shutdown of public-supply wells in Kings, Queens, and Nassau Counties. A large saltwater intrusion wedge in the Lloyd aquifer was indicated in southern Queens County in the early part of the 20th century, and the saltwater interface may have been onshore predevelopment. Most of Kings and Queens Counties are intruded with saltwater in both the upper glacial-Jameco-Magothy and Lloyd-North Shore aquifers systems. Saltwater increased during the 20th century and continues to increase to the present (2023) in the Lloyd-North Shore aquifer system in Great Neck and Manhasset Neck in northern Nassau County. A major wedge of saltwater intrusion in the upper glacial-Jameco-Magothy aquifer in southwestern Nassau County appears to be increasing.

New York↗
USGS2025

Phase 1 technical implementation plan for the expansion of the ShakeAlert earthquake early warning system to Alaska

Executive Summary The conference report accompanying the fiscal year (FY) 2022 Consolidated Appropriations Act (Public Law 117–103) for the U.S. Department of the Interior and related agencies directed the U.S. Geological Survey (USGS) to “work with the State of Alaska to develop an implementation plan to be completed within two years in order to put ShakeAlert/Earthquake Early Warning in Alaska” (p. 29). Congress included $1 million in the FY 2022 appropriation to conduct this effort. The USGS Earthquake Hazards Program, along with partner organizations, has developed the ShakeAlert earthquake early warning (EEW) system for the West Coast, which currently operates in California, Oregon, and Washington. The purpose of the system and its alert delivery partners is to reduce the impact of earthquakes and save lives and property by delivering ShakeAlert-powered alerts that are transmitted to the public via mass notification technologies, and by providing more detailed data streams to institutional users and commercial service providers to trigger automated, user-specific, protective actions. ShakeAlert was designed in such a way that it could be expanded to other U.S. regions with high earthquake risk, after the build-out of seismic and geodetic networks to support ShakeAlert in a specified region is completed and the necessary funding is secured for long-term operation and maintenance. When an earthquake occurs, seismic waves radiate from the rupturing fault like waves on a pond. It is these waves that people feel as earthquake shaking and that can cause damage to structures. Using networks of ground-motion sensors and sophisticated computer algorithms, ShakeAlert can detect an earthquake seconds after it begins, calculate its location and magnitude, and estimate the resulting intensity of shaking. Early warnings of impending shaking are then sent to people and systems that may experience damaging shaking, allowing them to take appropriate protective actions. Depending on the user’s distance from the earthquake, alerts may be delivered before, during, or after the arrival of strong shaking. There will almost always be a region near the earthquake epicenter where alerts arrive after damaging shaking has begun. The ShakeAlert system updates its ground-motion estimates as an earthquake grows larger. In response to the FY 2022 congressional direction, the USGS worked with the State of Alaska to devise this implementation plan for ShakeAlert expansion to Alaska. The USGS engaged with the Alaska Division of Homeland Security and Emergency Management (DHS&EM) and the Alaska Division of Geological and Geophysical Surveys (DGGS). A cooperative agreement was awarded to the Alaska Earthquake Center (AEC) at the University of Alaska Fairbanks (UAF) for their contributions to the plan and their work coordinating with other networks in Alaska. The USGS engaged with the Alaska Seismic Hazards Safety Commission (ASHSC) throughout the process. The USGS also held a series of Alaska stakeholder engagements. The process of developing the implementation plan was facilitated by contracted staff from Corner Alliance, which is a government consulting firm. This implementation plan describes the details and estimates the costs for a Phase 1 expansion of the ShakeAlert system to Alaska. A geographically limited Phase 1 goal was chosen that covers the highest risk and most populated areas of Alaska. The areas proposed encompass the State’s main population centers and 90 percent of the State’s population. This Phase 1 design is considered very challenging and ambitious from the viewpoint of network operators. The lessons learned if this plan is implemented could be used to consider subsequent phases to expand EEW beyond Phase 1 in Alaska in the future. ShakeAlert is built on the foundation of the sensor networks and data processing infrastructure of the USGS-led Advanced National Seismic System (ANSS). This implementation plan calls for a total of 450 high-quality, real-time EEW-capable ANSS seismic stations in Alaska: 270 new stations, 160 upgraded stations, and 20 existing stations. These seismic station numbers are based on a station spacing of 10 kilometers (km) in urban areas, 20 km in seismic source areas that endanger population centers, and 40 km in other areas. The associated costs also include support for some EEW-capable global navigation satellite system (GNSS) stations, with a focus on improving warnings for large subduction zone earthquakes. For effective EEW, ShakeAlert requires low-latency, high-availability, robust telemetry links to deliver continuous, real-time data from field stations to the data centers. The Alaska data processing hardware infrastructure would follow the general design for fail-safe operation that is used for the ShakeAlert system on the West Coast. The ShakeAlert architecture uses two independent layers: the production layer for earthquake processing and the alert layer to make alerting decisions and serve alerts to users. This implementation plan includes two geographically separated data centers in Alaska, each with two fully independent production and alert layers using the same system design developed for the West Coast. As of March 2024, the ShakeAlert system is at version 3.0.1, with more advanced versions in the development and testing pipeline. ShakeAlert originally used two algorithms to determine the location and magnitude of earthquakes using seismic data. A third algorithm that can calculate very large magnitudes of very large earthquakes with geodetic data was added in March 2024. ShakeAlert publishes several data and alert products to meet the needs of different users. All messages include the location of the earthquake, either as a point or a line, and its magnitude. Ground-shaking estimates are published in two forms, as ground-motion contours and a map grid. Providing adequate warning time for strong shaking (the “target threshold”) requires sending alerts at a threshold lower than that strong shaking level (the “alert threshold”). The thresholds for public alerting in Alaska would be a joint USGS and State decision. To have the greatest benefit, ShakeAlert-powered alerts would be delivered to institutional users and individuals by all practical pathways. The USGS alert layer can support thousands of institutional users and alert redistributors, but the USGS does not have the mission nor the infrastructure and expertise to perform mass notifications to the public or implement automatic actions for end users of the alerts. To meet this need, ShakeAlert recruits private sector “technology enablers” that have the necessary expertise to develop end-user implementations using EEW alerts with the goal of stimulating an EEW industry. Earthquake early warning alerts are useless if people do not know how to respond to them. Although the alert messages include instructions about what to do (drop, cover, and hold on), alerts are more effective if people have been trained in advance. Messages about ShakeAlert’s capabilities, limitations, and benefits could be integrated with existing earthquake education programs, including State-run programs. Therefore, ShakeAlert would coordinate with both public and private partners and stakeholders through various partnerships and agreements to accomplish consistent and ongoing public earthquake hazard education. The estimated capital cost of completing the computing infrastructure and sensor networks for the Phase 1 ShakeAlert expansion to Alaska is approximately $66 million in 2024 dollars. The annual operation and maintenance cost of the completed system is estimated to be $12 million per year in 2024 dollars when fully built out.

Alaska↗
USGS2024

Global geologic map of Europa

Discovered by Galileo Galilei more than 400 years ago and imaged in detail by the Voyager 2 Galileo spacecraft, Jupiter’s icy moon Europa has been a source of intrigue. A range of science investigations indicate that it contains the key ingredients for habitability, notably energy, chemistry, and liquid water. Europa’s surface is geologically complex and, based on the dearth of impact craters, interpreted to be as young as ~60 Ma. The array of geologic features that characterize the surface include extensive ridged plains, regions of broad disruption termed chaotic terrain, long, quasi-linear ridges that span thousands of kilometers, and bands as much as 60 kilometers wide and that extend 100s of kilometers. These features, along with other geophysical measurements, indicate the presence of a global briny liquid water ocean beneath the ice shell. It was not until the arrival of the Galileo spacecraft in 1995 that the true nature and level of complexity of the surface was revealed. Although image data returned by Galileo provided insight into the structure of a variety of regions, the entire satellite has yet to be observed at a regional scale (less than 250 meters per pixel) and the detailed geologic nature of much of its surface remains a mystery. Establishing the global context of the distribution and timing of Europan geologic units forms a basis to understand regional and local scale processes, serves as a tool for the planning of future missions, and most of all is essential to gaining insight into the potential habitability of this icy world.

Scientific Investigations Map↗
USGS2025

A new groundwater energy transport model for the MODFLOW hydrologic simulator

Heat transport in the subsurface is an important aspect of research related to the effects of a warming climate on ecological services (i.e., cold-water refugia); the development of geothermal resources for energy banking schemes (i.e., aquifer thermal energy storage [ATES]); and the effects of temperature on other aspects of groundwater quality, such as nutrient cycling. Historically, simulation of heat transport using the MODFLOW groundwater simulator and related codes was performed by scaling the input parameters of a solute-transport model to emulate heat transport. However, that approach required additional pre- and post-processing of input and output and could not account for the variation in effective thermal storage and transport properties during transient, unsaturated flow, for example. True heat-transport capabilities in the context of MODFLOW were first introduced in a variant called USG-Transport. More recently, a new groundwater energy-transport (GWE) model type has been added to MODFLOW 6, the core version of the MODFLOW hydrologic simulator. GWE supports the simulation of heat transport on structured or unstructured grids as well as within and between features of advanced packages that represent streams, lakes, multi-aquifer wells, and the unsaturated zone. GWE is integrated within MODFLOW 6 and is accessible through the FloPy Python package and the MODFLOW 6 application programming interface (API). An example simulation demonstrates conduction between grid cells through both the water and the solid aquifer material, including thermal bleeding from saturated overburden cells into a groundwater flow field.

Groundwater↗
USGS2025

Characterizing the scale of regional landslide triggering from storm hydrometeorology

Rainfall strongly affects landslide triggering; however, understanding how storm characteristics relate to the severity of landslides at the regional scale has thus far remained unclear, despite the societal benefits that would result from defining this relationship. As mapped landslide inventories typically cover a small region relative to a storm system, here we develop a dimensionless index for landslide-inducing rainfall, A * , based on extremes of modeled soil water relative to its local climatology. We calibrate A * using four landslide inventories, comprising over 11 000 individual landslides over four unique storm events, and find that a common threshold can be applied to estimate regional shallow-landslide-triggering potential across diverse climatic regimes in California (USA). We then use the spatial distribution of A * , along with topography, to calculate the landslide potential area (LPA) for nine landslide-inducing storm events over the past 20 years, and we test whether atmospheric metrics describing the strength of landfalling storms, such as integrated water vapor transport, correlate with the magnitude of hazardous landslide-inducing rainfall. We find that although the events with the largest LPA do occur during exceptional atmospheric river (AR) storms, the strength of landfalling atmospheric rivers does not scale neatly with landslide potential area, and even exceptionally strong ARs may yield minimal landslide impacts. Other factors, such as antecedent soil moisture driven by storm frequency and mesoscale precipitation features within storms, are instead more likely to dictate the patterns of landslide-generating rainfall throughout the state.

California↗