USGS Science⌕ Search

FIND YOUR NEXT DISCOVERY

Results for “Remote Sensing Applications: Society and Environment”

Original records, connected by a shared subject.

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

1,224 records · Page 3Linked to original sources
USGS2026

Groundwater drought in the United States: Spatial and temporal variability

Many communities and ecosystems in the United States that are dependent on groundwater are potentially adversely affected by groundwater drought. We computed yearly groundwater-drought metrics and mean groundwater levels at well locations across the conterminous United States (CONUS), using data from wells and remotely sensed and modeled Gravity Recovery and Climate Experiment Drought Monitor Data Assimilation (GRACE-DADM). We also modeled the probability of low or high human impact at each well location. The spatial distribution of groundwater-drought duration and severity from 2001 to 2020 for 1,510 wells shows longer maximum duration and higher maximum severity events in drier regions like the Southwest than in wetter regions like the Northeast. Based on 613 wells in CONUS from 1981 to 2020, there are many significant decreases in drought duration and severity in the Northeast and many significant increases in annual-mean groundwater levels. In contrast, there are many significant increases in drought metrics and decreases in mean water levels in parts of the Southeast. There are major differences in trends from 2001 to 2020 between well-based and GRACE-DADM-based groundwater metrics in some CONUS regions and a very low correlation between trends at individual locations across CONUS. A potential reason for this disparity is the low GRACE-DADM resolution (∼12 km) and the potential for a large amount of groundwater variation at the local scale. Also, GRACE-DADM represents shallow, unconfined aquifers which may not match the screened interval of the monitoring wells we evaluated. Large spatial gaps in long-term, high frequency, and quality-assured groundwater-well monitoring data present a challenge for understanding groundwater-drought variability across CONUS. Remote sensing tools such as GRACE can help but cannot fully replace well monitoring, as highlighted by our study results. Substantially more long-term monitoring wells would more accurately represent groundwater-drought trends and spatial variability across CONUS, particularly in western regions.

conterminous United States↗
USGS2026

The U.S. Geological Survey 2025 Puerto Rico and U.S. Virgin Islands time-independent earthquake rupture forecast

We present the 2025 U.S. Geological Survey Puerto Rico and U.S. Virgin Islands (PRVI) time‐independent earthquake rupture forecast (ERF), developed for the 2025 update to the National Seismic Hazard Model (NSHM) for PRVI. The updated ERF improves upon a prior model from 2003, including an expanded fault inventory with slip‐rate estimates, updated seismicity catalogs, and refined subduction zone geometries and deformation models. It applies the fault‐system inversion methodology to solve for rates of ruptures on modeled faults, adapted from the 2023 NSHM (NSHM23) for the western United States, including the first application of the inversion to model rates on a U.S. subduction interface. Off‐fault and intraslab seismicity are constrained by observed seismicity and use updated methods developed for NSHM23. Uncertainties in model components are substantial, and the ERF represents epistemic uncertainties through a comprehensive logic tree consisting of 1.7 billion logic‐tree branches combined across all sources.

Puerto Rico, U.S. Virgin Islands↗
USGS2026

Restoring dryland water cycles for precipitation feedback and climate stability; a review

Drylands across the globe are experiencing intensifying water scarcity, land degradation, and hydroclimatic extremes. This review integrates evidence from multidecadal field studies, hydrologic monitoring, geomorphic and ecological assessments, remote sensing, and land–atmosphere science to evaluate how restoration influences key components of the terrestrial water cycle. Low-tech natural infrastructure in dryland streams (NIDS)—including check dams, leaky weirs, one-rock dams, and gabions—has emerged as a promising but under-synthesized nature-based solution for restoring hydrologic function in these environments. We describe the mechanisms through which these interventions modify runoff detention, infiltration, sediment and alluvial storage, shallow-groundwater recharge, vegetation recovery, and surface-energy partitioning, and we summarize outcomes across diverse dryland settings. Findings consistently show increased water residence time, enhanced soil-moisture storage, expanded riparian vegetation, extended flow duration, and shifts toward greater latent-heat flux—producing localized cooling and strengthened ecohydrological feedbacks. Building on these localized effects, we articulate a hypothesis that links the spatial extent of restoration, the density of NIDS per unit drainage area, and the magnitude of the latent-to-sensible-heat contrast generated by wetter post-rainfall conditions. Specifically, we hypothesize that when NIDS are implemented at densities permitted by topography and across areas large enough to maintain elevated soil moisture after storm events, the resulting increases in latent heat flux, surface cooling, and boundary-layer moistening may enhance moisture convergence and boundary-layer development, potentially increasing the likelihood or stability of convective precipitation, analogous to how reductions in these processes have contributed to regional drought intensification. These land–atmosphere feedbacks remain untested at scale but represent an important research Frontier. By integrating hydrologic, geomorphic, ecological, and atmospheric perspectives, this review provides a comprehensive framework for considering how low-tech, landscape-scale interventions can strengthen watershed resilience and contribute to climate-relevant nature-based solutions.

Arizona↗
USGS2026

Orientation dependence of probabilistic seismic hazard estimates from CyberShake physics-based simulations

Earthquake ground‐motion intensities, such as pseudospectral accelerations (SAs), can vary significantly with horizontal orientation. However, conducting probabilistic seismic hazard analysis (PSHA) for each horizontal orientation is challenging because current ground‐motion models used in PSHA consider only a single horizontal intensity value, usually the median across all orientations, known as RotD50. To address this limitation, we employ physics‐based simulations for PSHA, which contain full waveforms from which ground‐motion intensities can be computed for all horizontal orientations to study directional seismic hazard. We apply our approach to the latest CyberShake study of the Greater Los Angeles metropolitan area, developed by the Statewide California Earthquake Center, finding that seismic hazard at a 2475‐yr return period, a common value used for earthquake‐resistant design, varies significantly with horizontal orientation. For instance, for SAs at 3 s, the maximum seismic hazard across all horizontal orientations is, on average, 15% higher than the median RotD50 hazard, with these differences becoming more pronounced at longer periods. These observed variations can generally be attributed to physical mechanisms that polarize seismic waves, such as the radiation pattern of the earthquake source and the influence of the subsurface structure. These results may have important implications for earthquake engineering applications, particularly for long‐period structures in areas with substantial horizontal variations in seismic hazard.

California↗
USGS2024

Investigation of land cover within wetland complexes at Dixie Meadows, Churchill County, Nevada, from October 2015 to January 2022

The U.S. Geological Survey investigated land cover at subannual time steps within six wetland areas in Dixie Valley, Churchill County, Nevada, from October 2015 to January 2022. As requested by the U.S. Fish and Wildlife Service, we used aerial photography and satellite remote sensing data to map surface water and other land cover types within the wetland complexes. We identified five land cover classes using the green normalized difference vegetation index (gNDVI) and its inverse relationship to the normalized difference water index (NDWI) within three U.S. Department of Agriculture National Agriculture Imagery Program aerial images (acquired in 2015, 2017, and 2019) and 110 European Space Agency Sentinel-2 satellite images (acquired 2015–2022). The relative wetness of soil conditions within each land cover class is estimated by comparison to previously published observations of relative conductivity measured by 79 field-based sensors within the wetlands from 2019 to 2021. We mapped the areal coverage of the five land cover classes for approximately 385 acres (1,559,000 square meters [m²]) comprising six individual wetland complexes as well as a larger 1,298- acre (5,254,000-m 2 ) area of interest inclusive of the wetland complexes and adjacent landscape. Land cover of open water (Class 5) primarily within ponds at one of the wetland complexes comprised 8,333 m 2 , on average, of the wetland complexes. Land cover of mixed shallow surface water, saturated soil, and vegetation (Class 4) comprised 111,723 m 2 on average of the wetland complexes. Land cover of dense green vegetation canopy cover (Class 3) that often (46 percent of observations) had underlying surface water or saturated soil conditions comprised 592,522 m 2 on average of the wetland complexes. The remaining areas of the wetland complexes not mapped as these three land cover types (Classes 2 and 1) had sparse vegetation or bare soil cover and commonly (greater than or equal to 67 percent of observations) had dry soil conditions. The investigation of land cover detailed in this report could inform future efforts to map land cover more precisely via higher resolution remote sensing or ground-based surveying or could be incorporated with other environmental monitoring data to characterize habitat and hydrology of the wetland complexes at Dixie Meadows.

Nevada↗
USGS2026

Perspectives and insights gained from applying landscape ecological approaches to the Upper Mississippi riverscape

Advances in remote sensing, spatial data collection, hydraulic and geospatial models, and the maturation of landscape ecology during the early 2000s spurred the concept of rivers as landscapes or ‘riverscapes’. Since that time, the Upper Mississippi River Restoration (UMRR) Program - one of the world’s longest running large-river ecosystem restoration and monitoring programs - has applied landscape ecological principles to the Upper Mississippi River System (UMRS). I examine spatial patterns of aquatic nutrient concentrations, submersed aquatic vegetation, and fish and mussel communities, finding that within-river patterns are often patchy in nature, likely resulting from how nutrient processes and community traits interact with spatial variability in hydraulic exchange patterns. Local-scale patchy spatial patterns in the UMRS are embedded in system-scale longitudinal gradients and a hierarchical view of the river has provided a foundation for strategic river restoration planning. Key Words: Aquatic Habitat, hydraulic connectivity, patch, restoration, riverscape

BioScience↗
USGS2025

A review of post-wildfire adaptations of surface-water-quality models: Synthesis, gaps, and opportunities

As wildfires increasingly affect water-supply watersheds, the demand for models to predict water-quality responses is increasing. This work reviews and synthesizes existing post-wildfire applications of water-quality models in the context of geographic and ecohydrological distribution, hydrologic and water-quality response process representation, model parameterization, model and input data scales, model calibration data availability, as well as calibration and performance evaluation approaches. Emphasis is placed on models that simulate water-quality output, rather than sediment and erosional response as the primary focus. Here, identified gaps and opportunities to advance the post-wildfire application of water-quality models include: 1. applying models in under-represented geographic and ecohydrologic regions, 2. simulating multiple streamflow generation mechanisms, including groundwater, with an emphasis on shifting dominant flow pathways as the landscape recovers following wildfire, 3. adding studies that include the simulation of metals, 4. incorporating more biogeochemical and in-stream processes to model applications, 5. applying finer spatial and temporal resolution of precipitation data input as well as finer spatial resolution hydrologic response units, 6. implementing fully distributed grid or element models or finer resolution response units to capture burn severity heterogeneity, 7. collecting enhanced water-quality data for model calibration and validation, 8. conducting model-intercomparison studies, and 9. developing model parameter value guidance in post-wildfire applications. These identified gaps and opportunities may assist users in deciding on key processes and approaches to consider in modeling post-wildfire water-quality conditions.

Science of the Total Environment↗
USGS2025

Relationship between peak and cumulative ground motions from 49 Mw 3-6 earthquakes in the San Francisco Bay Area

We examine the relationship between peak ground velocity (PGV) and cumulative absolute displacement (CAD) for 49 M w 3 – 6 earthquakes in the San Francisco Bay Area (SFBA) and gain insight into the spatiotemporal partitioning of seismic energy in ground motion records with respect to source, path, and site effects. PGV and CAD are positively correlated, but there can be large deviations from the average trend. For example, ground motion records with either very long duration resonance or short pulse-like motions will have higher or lower CAD, respectively, but could have very similar PGV. We perform principal component analysis (PCA) on PGV-CAD for >7000 records in the SFBA with the goal of investigating what influences positive or negative anomalies in cumulative motions. PCA rotates the PGV-CAD datapoints into two principal components, where the one with the larger variance, which we call the “primary intensity” represents mostly the distance-dependence of ground motion amplitudes. The other principal component, which we call the “excess motion”, represents the deviation from cumulative motions that would be typical for a ground motion record with a given PGV. The excess motion will be positive in the case of records with long duration ringing and will be negative for short duration pulse-like ground motions. We find that excess motion is generally positive at sites in sedimentary basins and in soft sediments around the SF Bay. Excess motion is generally negative in the very near field, as well as at sites on hard bedrock. We discuss the findings here in terms of implications for seismic hazard applications and other wave propagation phenomena.

California↗
USGS2025

Seismic moment and local magnitude scales in Ridgecrest, CA from the SCEC/USGS Community Stress Drop Validation Study

We illustrate the systematic difference between moment magnitude and local magnitude caused by underlying earthquake source physics, using seismic moments submitted to the Statewide California Earthquake Center/United States Geological Survey Community Stress Drop Validation Study 2019 Ridgecrest data set. While the relationship between seismic moment and moment magnitude ( M or M w ) of log 10 ( M 0 ) ~ 1.5* M is uniformly valid for all earthquake sizes by definition (Hanks and Kanamori, 1979), the relationship between local magnitude M L and moment is itself magnitude dependent. For moderate events, ~3< M < ~6, M and M L are coincident; for earthquakes smaller than ~3, M L ~ 1.0 log 10 M 0 (Hanks and Boore, 1984). This is a physical consequence of the corner frequency fc becoming larger than the upper frequency of observation and implies that M L and M differ systematically by a factor of 1.5 for these small events. While this idea is not new, we propose a new, continuous relationship between local magnitude and moment, for magnitudes 2 to 6 which extrapolates to smaller and larger magnitudes, applicable to southern California specific to the Ridgecrest region. We make use of the plethora of seismic moments as submitted by many participants of the Community Stress Drop study, compared to the Southern California Seismic Network (SCSN) catalog magnitudes. Overall, the seismic moments in the Community Study recover moment magnitude well, so we use our new M L - M 0 to convert M L to M , refining the SCSN operational M Lr scale. This systematic difference of 50% in slope between local and moment magnitude at small magnitudes has implications for spectral stress drop estimates, earthquake ground motion modeling, as well as other magnitude scales and earthquake occurrence statistics.

California↗
USGS2025

Combining multisite tsunami and deformation modeling to constrain slip distributions for the 1700 C.E. Cascadia earthquake

A major earthquake ruptured the Cascadia subduction zone (CSZ) on 26 January 1700. Key paleoseismic evidence associated with this event include tsunami deposits, stratigraphic evidence of coastal coseismic subsidence, written Japanese records of a tsunami unaccompanied by earthquake shaking, and margin‐wide turbidites found offshore and in lacustrine environments. Despite this wealth of independent clues, important details about this event remain unresolved. Dating uncertainties do not conclusively establish whether the proxies are from one earthquake or a sequence of them, and we have limited knowledge of the likely slip distributions of the event or events. Here, we use a catalog of 37,500 candidate synthetic ruptures between M w 7.8 and 9.2 and simulate their resulting coseismic deformation and tsunami inundation. Each model is then compared against estimated Japan tsunami arrivals, regional coastal subsidence records, and local paleotsunami deposits mapped at six different coastal marshes and one coastal lake along the CSZ. We find that seven full‐margin ruptures with a median magnitude of M w 9.1 satisfy all three constraints. We favor one M w 9.11 model that best matches all site paleoseismic observations and suggests that the Cascadia megathrust slipped up to ∼30 m and must have shallow geodetic coupling. We also find that some sequences composed of three or four ruptures can still satisfy the observations, yet no sequences of two ruptures can. Sequences are differentiated into three groups based on whether they contain a mainshock rupture located in the south (>44° N) or further north. All sequences contain unruptured portions of the megathrust and most contain mainshocks with peak slip above 40 m. The fit of the geologic evidence from sequences is poor in comparison to single‐event models. Therefore, sequences are generally less favored compared to full‐margin events.

British Columbia, California, Oregon, Washington↗
USGS2026

Simulation-based scenario ShakeMaps for large magnitude (MW6.5+) crustal earthquakes on the Seattle, Tacoma, and southern Whidbey Island faults, Washington, USA

Scenario ground‐motion maps based on empirical ground‐motion models (GMMs) provide a rapid and generally reliable means of estimating the amplitude and distribution of earthquake shaking. However, because GMMs are designed for broad applicability, they often rely on simplified representations of Earth structure, which can limit their accuracy in regions with complex source, path, and site effects. This can substantially impact the accuracy of predicted shaking in areas like western Washington State, where deep, interconnected basin structure exerts a strong influence on seismic‐wave propagation. In this study, we present a new suite of simulation‐based scenario ShakeMaps that characterize ground shaking from large‐magnitude ( ⁠ M W 6.5–7.5) earthquakes on the Seattle, Tacoma, and southern Whidbey Island faults. These maps are developed using results from recent 3D wave propagation simulations ( Stone et al. , 2022 , 2023 , 2025 ) that incorporate realistic rupture geometries, variable slip distributions, and a regional 3D seismic velocity model with shallow soils. Broadband ground motions are estimated by combining the low‐frequency (<1 Hz) deterministic seismograms from these studies with high‐frequency (1–10 Hz) stochastic seismograms. Simulated ground motions are corrected to account for the enforced minimum shear‐wave velocity and nonlinear site response. The resulting ShakeMaps represent median ground‐shaking estimates derived from multiple rupture scenarios with varying slip distributions and hypocenter locations for each fault. To extend ShakeMap coverage beyond the simulation domain (i.e., into eastern Washington, northern Oregon, and southwestern British Columbia), we scale GMM‐based ground‐motion estimates using amplification patterns observed in the simulations. These new ShakeMaps reveal the substantial influence of deep basin structure on shaking intensity, underscoring the importance of considering crustal structure complexity in regional hazard assessments for the Pacific Northwest.

Washington↗
USGS2025

Prospectivity modeling of the NASA VIPER landing site at Mons Mouton near the Lunar South Pole

We use a high-resolution digital elevation model and a numerical thermal model to produce a variety of inputs for a water-ice prospectivity model for the Volatiles Investigating Polar Exploration Rover (VIPER) landing site. These input data are maps of topography, surface slope, surface aspect, surface curvature, maximum temperature, depth to ice stability, permanently shadowed regions (PSRs), distance to PSRs, and PSR density. This model predicts where water ice is most likely within the top meter of regolith, assuming plausible relationships between ice concentration and the various inputs. The model is designed to be adjusted in near-real time as data are collected during the VIPER mission. As such, it is a tool for both analyzing data from the mission as well as planning operations. Since the current model, at this point, relies only on orbital remote sensing, the final version will also be a tool to extrapolate the VIPER mission results across the lunar poles.

Planetary Science Journal↗
USGS2025

Near-surface geophysics: Environmental applications

The field of geophysics encompasses a broad and diverse compilation of methodologies that employs principles of physics to characterize properties of earth materials within the subsurface. While geophysical methods have a long history in resource exploration and studies of Earth’s interior, the subdiscipline of “near-surface geophysics” has evolved in recent decades for examination of the shallow, near-surface environment for a range of purposes ranging from archaeological or forensic investigations to assessment of geologic, hydrologic, biologic, and geochemical properties and processes. “Environmental geophysics” are near-surface geophysical studies and methods that focus on understanding natural systems (e.g., watershed hydrology, groundwater–surface water connections, biophysical processes) as well as research pertaining to anthropogenic impacts and land management, (e.g., contamination and remediation, saltwater intrusion, agricultural practices). This field can be further subdivided into subdisciplines focused on specific topics and applications, such as water resources and hydrology (hydrogeophysics) or biologic and microbial processes (biogeophysics). Studies in environmental geophysics span a range of scales, from pore-scale laboratory tests to watershed-scale or regional field experiments. Methods vary by the nature of physics employed, the specific measurement acquired, and how that data is ultimately processed and analyzed to produce interpretable results. There exists further diversity in the acquisition logistics, geometry, and timing of data collection. Geophysical data can be collected in boreholes (one-dimensional, 1-D, vertical profiles), along survey lines (two-dimensional, 2-D, cross-sections), or in dense sensor arrays or gridded profiles (three-dimensional, 3-D, models). Regarding the temporal aspect, studies can conduct one-time geophysical surveys to obtain detailed imaging of subsurface structure or use timelapse and continuous monitoring to investigate variations in subsurface properties over time. The cumulation of all possible permutations of these factors (method, acquisition geometry, survey design, and target application) results in an immense diversity among environmental geophysical studies. Nevertheless, this field remains unified in the pursuit of understanding natural and human-impacted near-surface environments through geophysical investigations. Here we highlight some key references within environmental geophysics. Resources on geophysical theory, acquisition logistics, processing and inversion workflows, and example case studies are categorized into the most common geophysical classes within Geophysical Methods. Lastly, example references for the dominant types of applications in environmental geophysical studies are catalogued in Environmental Applications.

Book chapter↗
USGS2025

Protected from Pterygoplichthys? Predicting thermal habitat suitability for nonnative armored catfish in the Suwannee River

Objective Nonnative fishes can modify ecosystems and harm economies when they are introduced to new environments. Climate change is likely to assist the spread and establishment of some nonnative fishes (e.g., warmwater species), but spatiotemporal gaps in water temperature monitoring and modeling may prevent ecologists and managers from forecasting thermal habitat suitability for these taxa. The purpose of this study was to develop a predictive model of winter water temperatures and thermal habitat suitability for two priority nonnative armored catfish, Vermiculated Sailfin Catfish Pterygoplichthys disjunctivus and Orinoco Sailfin Catfish P. multiradiatus , in the Suwannee River, Florida and Georgia. Methods Precipitation- and groundwater-corrected air–water temperature models were developed and evaluated using a model selection procedure to predict water temperatures at four sites in the Suwannee River. These models were chosen because they blend the simplicity of air–water temperature models with the accuracy of hydrometeorological models to create an efficient, economical, management-relevant approach for analyzing and forecasting water temperature. Results Most of the top-performing water temperature models (92%) had precipitation or groundwater corrections to air–water temperature formulations. Projected mean and maximum water temperatures increased as simulated climate change intensified. All four Suwannee River sites studied were projected to be thermally hospitable to the survival of Vermiculated Sailfin Catfish. Lower river sites, noticeably warmer than upper river sites, were conducive to the survival of Orinoco Sailfin Catfish throughout the winter months. The upper river sites were too cold for Orinoco Sailfin Catfish survival in some climate-change scenarios, but the Suwannee River has an abundance of constant-temperature springs that are likely hospitable to Vermiculated Sailfin Catfish and Orinoco Sailfin Catfish throughout the year. Conclusions The findings suggest that winter water temperatures will likely not be a barrier to the survival of Pterygoplichthys catfish in the Suwannee River, amplifying the importance of conservation and management approaches to inhibit their spread and establishment. If the Pterygoplichthys population remains small and isolated and decision makers are able to devote required staff time and resources to managing these species, removal and eradication at local if not broader scales may be reasonable goals. This study provides a water temperature modeling approach that can aid ecologists and managers in prioritizing sites to prevent the introduction, slow the dispersal, eradicate, and control Pterygoplichthys catfish and other nonnative fishes in the Suwannee River and beyond.

Florida, Georgia↗
USGS2026

Understanding the resource potential of natural hydrogen on Earth: Scientific gaps, uncertainties and recommendations

A comprehensive scientific research roadmap is essential to bridge knowledge gaps and deepen the understanding of key geological, geochemical, and geophysical aspects of natural hydrogen (H 2 ) as a potential new energy resource. This paper reviews major scientific uncertainties on natural H 2 , suggesting research priorities, as a guide for defining exploration strategies, techniques, and data interpretation. The uncertainties concern all phases of the natural H 2 cycle, from generation (source rocks) through migration (advection and diffusion) and accumulation (reservoir and cap rocks) to the application and interpretation of subsurface and surface geochemical and geophysical exploration techniques. Understanding H 2 sources and generation rates (the amount of H 2 generated by a given volume of rock over time) is crucial for determining whether a geological H 2 system operates as a short-term dynamic system with rapid H 2 production and release, or as a conventional gas system with long-term accumulations, analogous to petroleum reservoirs. Preliminary estimates for serpentinisation, radiolysis, and organic matter degradation suggest that H 2 generation is not inherently fast, especially for non-hydrothermal continental systems (crystalline basement of shields, ophiolites, peridotite massifs, sedimentary basins), and long-term accumulations, like those of fossil natural gas systems, represent the most likely scenario. The mechanisms of H 2 migration through geological formations require application of fundamental principles of fluid-flow physics, distinguishing advection and diffusion, as well as their forms (from gas-phase, bubble flows to aqueous solutions). Additional studies of H 2 accumulation and retention in subsurface reservoirs could improve understanding of mechanisms of H 2 migration by focusing on the rock fluid-bearing properties and the factors affecting H 2 preservation, such as the presence of cap rocks impermeable to H 2 , pressure conditions, residence times, and microbial or abiotic consumption. Advanced techniques, including reservoir modelling, flow simulations, 3D imaging (micro-CT) of H 2 -bearing rocks, and extraction and analysis of gas occluded in rocks, can provide insights into the stability and potential recoverability of H 2 accumulations. The interpretation of surface exploration techniques, including gas geochemistry, geophysics, and remote sensing, long employed in mineral and energy resource exploration, is now being adapted for natural H 2 studies, but challenges remain in the data interpretation. Distinguishing H 2 seepage due to geological degassing from H 2 produced near the surface by modern microbial processes or artificial sources, such as hammering or drilling for soil-gas sampling, drilling into aquifers, and corrosion in boreholes, is an essential step in exploration. The simple detection of H 2 in soils, even in morphological structures like sub-circular depressions or “fairy circles”, cannot be cursorily interpreted as a signal of natural H 2 seepage from a deep source. A holistic geochemical approach, including isotopic analyses of gases associated with H 2 , is recommended to distinguish among the variety of possible H 2 origins. Observations of H 2 in wells should be interrogated to rule out possible artifacts such as corrosion and drill bit metamorphism. The integration of multiple geophysical methods, including seismic, gravimetric, magnetic, and electro-magnetic surveys, is recommended to mitigate interpretation ambiguities regarding the structure of a subsurface H 2 system (source and reservoir rocks, including fluid and gas storage), due to the non-uniqueness of rock-specific physical properties.

Earth-Science Reviews↗
USGS2025

Over, under, and through: Hydrologic connectivity and the future of coastal landscape salinization

Seawater intrusion (SWI) affects coastal landscapes worldwide. Here we describe the hydrologic pathways through which SWI occurs - over land via storm surge or tidal flooding, under land via groundwater transport, and through watersheds via natural and artificial surface water channels—and how human modifications to those pathways alter patterns of SWI. We present an approach to advance understanding of spatiotemporal patterns of salinization that integrates these hydrologic pathways, their interactions, and how humans modify them. We use examples across the East Coast of the United States that exemplify mechanisms of salinization that have been reported around the planet to illustrate how hydrologic connectivity and human modifications alter patterns of SWI. Finally, we suggest a path for advancing SWI science that includes (a) deploying standardized and well-distributed sensor networks at local to global scales that intentionally track SWI fronts, (b) employing remote sensing and geospatial imaging techniques targeted at integrating above and belowground patterns of SWI, and (c) continuing to develop data analysis and model-data fusion techniques to measure the extent, understand the effects, and predict the future of coastal salinization.

Water Resources Research↗
USGS2026

Groundwater dependency and hydroclimatic influences on riparian and upland vegetation productivity, Upper San Pedro, Arizona, United States

In arid and semi-arid regions, groundwater sustains vegetation through subsurface water access, yet the responses of groundwater-dependent ecosystems (GDEs) to changing hydroclimate and groundwater availability are relatively understudied. This study investigates seasonal and spatial patterns in vegetation greenness using Landsat Enhanced Vegetation Index (EVI) values across riparian and upland zones in the semi-arid Upper San Pedro (USP) watershed, southern Arizona, which experiences a bimodal precipitation regime. We paired 25 years (2000–2024) of EVI and depth to groundwater (DTG) data from 89 wells and climate metrics (precipitation and vapour pressure deficit) to quantify the sensitivity of vegetation to subsurface moisture as well as atmospheric moisture supply and demand. Vegetation at wells near the USP riparian area showed strong associations between EVI and DTG anomalies during the monsoon season, indicating sustained groundwater use even during this wet period when summer precipitation is abundant. In contrast, upland vegetation that lacked access to groundwater showed minimal sensitivity in EVI to DTG and was generally less responsive to vapour pressure deficit. Interestingly, the riparian GDEs were not decoupled from precipitation and climate variability. These results underscore the importance of groundwater for maintaining riparian productivity and highlight the utility of remote sensing in identifying vegetation-climate-groundwater linkages across heterogeneous dryland landscapes.

Arizona↗
USGS2025

State shifts in the deep Critical Zone drive landscape evolution in volcanic terrains

Understanding the near-surface environment where atmospheric and solid earth processes interact, often termed the “Critical Zone,” is important for assessing resources and building resilient societies. Here, we examine a volcanic landscape in the Oregon Cascade Range, an understudied Critical Zone setting that is host to major regional water resources, pervasive silicate weathering, and significant geohazards. We leverage a bedrock age chronosequence to show that the volcanic Critical Zone undergoes a structural shift, from depth extents of >1 km to meters, over timescales of ~1 My. We map an active groundwater volume comparable to major continental lakes, stored at the Cascade Range crest. This state shift makes volcanic landscape evolution a unique probe of deep coupling between Earth systems.

PNAS↗