USGS Science⌕ Search

SEARCH · USGS Science

Results for “Soil Dynamics and Earthquake Engineering”

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,696 records · Page 30Linked to original sources
USGS2025

Concentration-discharge relations and transient metal loads reveal spatiotemporal variability in solute-generation mechanisms in a mine-affected watershed

Concentration-discharge (CQ) relations are commonly used to understand geochemical and hydrologic controls on the generation of solutes in watersheds. Despite the widespread application of CQ relations, this technique has been infrequently applied to acid mine drainage (AMD) sites, but the CQ framework may allow mechanistic understanding of remedial outcomes such as impoundment of water within underground mines. Results of CQ analyses and changes in metal loads in an AMD affected watershed in Colorado, USA indicate that dissolved loads increased at many individual locations following water impoundment within mine workings. Although increased loads were observed at most individual locations, these increases were offset by a large decrease in loading from the largest mine. A loading analysis that included data from an instream monitoring location showed a statistically significant decrease in Fe and Zn after bulkhead emplacement, indicating a net positive effect of bulkheads. Streams generally displayed dilution CQ patterns whereas mines and springs showed either flushing or chemostatic patterns prior to bulkheading, which transitioned to chemostatic patterns following bulkheading, indicating a transition from dynamic to equilibrium geochemical processes. Saturation indices for sulfide and secondary minerals indicated that mines and springs were near equilibrium for phases including schwertmannite, fluorite, and gypsum. Saturation indices vary through time for mines suggesting progressive leaching of sulfide minerals as the mass of available minerals in the mine workings decreases. Together, these diverse analyses provide an integrated understanding of the variability in solute generating processes in this watershed and may inform remediation plans for similarly affected sites by indicating the nature of mineralogic controls on water quality.

Colorado↗
USGS2025

The influence of postfire root strength decay on shallow landslide susceptibility in western Oregon

Severe wildfires can increase the potential for postfire landslides, partly due to the loss of vegetation and root reinforcement. This study investigates the multi-year dynamics of root strength following wildfire in coniferous forests of western Oregon and evaluates implications for postfire slope stability. Root tensile strength was measured through laboratory testing of 969 roots collected from 26 test pits over four postfire timeframes: unburned conditions, 4–12 months postfire (Cedar Creek, 2022), 24–35 months postfire (Holiday Farm, 2020), and 59–67 months postfire (Eagle Creek, 2017). Variables analyzed include root diameter, tensile thread strength, stiffness, progressive tensile strength with displacement, time-dependent strength and ductility, root area ratio, and cohesion. Results indicate a 50% reduction in root strength after wildfire, with the lowest values occurring around 4 years postfire. Projections indicate it could take 10–22 years following fire for root strength to recover to 70% and 90% of the prefire original root strength, respectively. This timeline indicates a multi-year window of vulnerability for shallow landslides that reaches a maximum approximately 4 years following fire. Brittle failure behavior was observed in burned root systems compared to unburned root systems, signaling a more abrupt loss of strength at yield. Simple slope stability analyses show that reduced root strength can lead to instability on more gentle slopes and under less saturation in comparison to unburned conditions. These findings highlight the critical role of root reinforcement in postfire slope stability and the long-term implications of wildfire disturbance on landslide susceptibility.

Oregon↗
USGS2026

Reconstructing the Quaternary depositional history using geologic mapping and three-dimensional modeling of the subsurface near Fort Morgan, northeastern Colorado

Centered on Fort Morgan, Colorado, this study is intended to build from previous work by adding a three-dimensional (3D) view of the subsurface to better understand the depositional history of Quaternary deposits. A 1:100,000 scale geologic map was made by combining previous geologic maps, regional soil maps, and recent field investigations. In addition to the geologic mapping, drill hole lithologic data from water wells and oil and gas exploration were compiled and lithologic units simplified to best represent the stratigraphy of the Quaternary deposits. From these subsurface data, a 3D subsurface model was constructed, trimmed at the surface by a digital elevation model, and a bedrock surface foundation gridded from drill hole data was added. The surface of the 3D model was then compared visually to the surficial geologic map. Cross sections were constructed from the 3D model and compared to site-specific drilling that was done as part of this project. Finally, the model was examined in detail to reconstruct the depositional history of the subsurface alluvial and eolian units. Alluvial and fluvial drainage basins exposed in the subsurface have a greater areal extent than the present-day narrow drainages. Older eolian sand in the subsurface tends to be interbedded with loess indicating coeval deposition. Holocene sand, both eroded from bedrock exposed at the surface north of the study area and reworked from the South Platte River, buries most of the interbedded older sand and loess.

Colorado↗
USGS2025

Airborne geophysical analysis to decipher salinization for coastal Louisiana

Coastal Louisiana is known for saltwater intrusion that threatens wetlands, aquifers, and rivers. However, the extent of saltwater intrusion is not well understood. This study develops an innovative framework with airborne electromagnetic (AEM) data to map chloride concentration distributions for wetlands in the Mississippi River deltaic plain and Chenier plain as well as for the Mississippi River Valley alluvial aquifer (MRVA) and Chicot aquifer. Moreover, the framework maps chloride concentrations along the Mississippi River and Atchafalaya River. Key components in the framework include the establishment of resistivity-to-chloride concentration transformation, 3D resistivity architecture building through geostatistics, and the employment of a lithologic model. The transformation functions correlate AEM resistivity data with porewater salinity measurements and groundwater and river chloride samples. The results show that AEM data reliably infers soil water chloride concentrations and correlates well with the distribution of various marsh types. AEM data reveals extensive saltwater presence at depth and near the coast, originating from salt domes and the Gulf of Mexico, respectively. The saltwater upconing pattern in the Chicot aquifer is likely due to excessive groundwater withdrawals. The AEM data also confirms a distinct tongue of saltwater intruding into the Atchafalaya Basin from the Gulf. The AEM data helps to identify faults that are obscured or eroded at the surface, which appear as leaky barriers in the subsurface where dramatic changes in chloride concentration are apparent. Finally, this study uses the AEM data to infer the presence of an extensive seawater wedge in the Mississippi River and Atchafalaya River.

Louisiana↗
USGS2026

The digital archivist: Automating legacy macroseismic data processing using large language models

Macroseismic data are a key resource to investigate shaking and damage from preinstrumental and early instrumental eras. However, data are often stored as inconsistently formatted reports describing observed shaking and damage, making manually parsing and interpreting accounts labor‐intensive. We introduce a novel workflow using Google’s Gemini 2.5 Pro large language model (LLM) to automate the extraction and structuring of macroseismic observations from summary reports. We apply this workflow to the 22 March 1957 M 5.3 Daly City, California, earthquake as a case study. We used Gemini to extract addresses, originally assigned modified Mercalli intensity values, and descriptions from each report. To address coordinate precision limits, addresses were geocoded via Google’s Geocoding application programming interface. This workflow yielded over 2300 geocoded intensity reports for the Daly City earthquake. We use the geocoded accounts, with the original report intensity assignments, to develop a shaking intensity map that in some respects rivals modern Did You Feel It? Maps. We also extract and present data for the 9 February 1971 M L 6.7 Sylmar, California, earthquake. Our results demonstrate the potential of LLMs for reliably extracting and analyzing large, unstructured macroseismic datasets. LLMs offer a scalable solution for rapidly digitizing macroseismic archives, enabling their broader use to constrain ground‐motion models in modern seismic hazard analysis and to improve our understanding of site effects in urban areas. The concepts explored here may also be applied to the handling of other legacy seismological and earth science data.

Seismological Research Letters↗
USGS2026

Only very strong shaking can break a tree

A handful of credible accounts, both recent and historical, describe trees being broken during earthquakes. Apart from landslides and tsunamis, such accounts are uncommon, suggesting that trees break only when ground acceleration is extreme or when trees are especially vulnerable to breakage. Here, I consider the question: what ground acceleration is required to break a living tree? I present simplified first‐principles calculations, considering the properties of common trees, focusing on tall, thin trees that can be approximated as uniform vertical cantilevers. The results suggest that, possibly excepting tall trees with heavy canopies, standing trees can only be broken by shaking approaching or exceeding 1 g , as anticipated by the rarity of documented accounts of snapped trees. I briefly consider several well‐documented instances of tree damage during strong earthquakes, drawing inferences about local ground motions.

The Seismic Record↗
USGS2026

A roadmap for identifying and interpreting physical processes and national water model prediction bias associated with baseflow index regimes across the contiguous United States

Understanding how groundwater–surface water interactions shape streamflow variability is critical for diagnosing low flow behavior and prediction bias in continental scale hydrologic models. We present a process informed framework that links observed baseflow (BF) dynamics, watershed attributes, and National Water Model (NWM) performance across the contiguous United States. Using daily observed streamflow from 797 reference quality streamgages, we developed monthly baseflow index (BFI) signatures using a streamgage specific, calibrated digital filter. Hierarchical clustering of these signatures identified seven distinct BFI regimes capturing regional and seasonal variability. We evaluated NWM v3.0 retrospective streamflow performance within each regime using multiple hydrograph and flow duration curve-based metrics. Model skill varied systematically across regimes: mixed flow systems were simulated most accurately, while predominantly BF dominated and quickflow dominated regimes exhibited substantially poorer performance. Across nearly all regimes, the NWM underestimated observed BFI magnitude and frequently failed to reproduce seasonal BF patterns, indicating systematic biases in simulated low flow contributions. To relate these regimes to potential process controls, we trained a Random Forest classifier using static watershed attributes and applied Shapley Additive Explanations to identify features most strongly associated with each regime. Results highlight regionally varying influences, including the dominant role of snow fraction and seasonal runoff timing in snow dominated basins and the importance of evapotranspiration and aridity in quickflow dominated systems. Collectively, these findings demonstrate how hydrologic signatures combined with interpretable machine learning can diagnose regime specific model biases and generate process-based hypotheses about limitations in large scale hydrologic prediction systems.

contiguous United States↗
USGS2025

Paleoseismology and paleogeodesy using coral microatolls

Establishing the rupture extent and slip distribution of individual paleo-earthquakes is vital for assessing fault behavior including the persistence of rupture segmentation, recurrence patterns, and similarity of successive events, key issues in both fault mechanics and hazard assessment. Techniques with high temporal and geodetic precision as well as a wide distribution of study sites are necessary to investigate past earthquakes in such detail. Coral microatoll growth is one of the best types of geologic record for paleoseismology and paleogeodesy given these needs, as it provides long, continuous, widely distributed records of centimeter-scale vertical tectonic motion with potentially annual-level temporal precision. This chapter describes the process of interpreting microatoll growth records to obtain time series of relative sea level, tectonic vertical deformation fields, and finally slip and coupling parameters on an underlying fault interface.

Book chapter↗
USGS2026

Surface-wave relocation and characterization of the 2024–2025 dyking episode along the Fentale–Dofen segment of the Ethiopian rift

From late-December 2024 to mid-March 2025, a 50-km-long dyke intrusion triggered over 300 earthquakes (magnitude 4 to 5.9) between Fentale and Dofen volcanoes along the Northern Main Ethiopian Rift. Dyke intrusions periodically occur along the Fentale–Dofen magmatic segment and are an expression of ongoing rift extension. Preliminary analyses using interferometric synthetic aperture radar revealed extensive ground deformation (up to 60 cm), which closely matched the temporal and spatial evolution of surface manifestations and earthquake locations from global catalogues. While global catalogues are critical for real-time monitoring, the precision of locations in remote and or sparsely instrumented regions can be low. In this investigation, we present surface-wave relocation results of the dyking episode that began near Fentale volcano in December 2024. We estimate relative locations using differential traveltimes measured from regional-to-teleseismic distance surface-wave observations of earthquakes reported by the U.S. Geological Survey. Relative relocations reduce the initial region of diffuse seismicity to a 50-km-long narrow band bounding the strike of surface manifestations and the zone of maximum surface deformation. We demonstrate the precision of surface-wave relocations over incremental time periods, capturing the progression of dyking from seismic onset through seismic migration and caldera subsidence. Results showcase the utility of surface-wave relocations in the characterization of dyking episodes and provide complementary insights into the current understanding of the Fentale–Dofen volcanic plumbing system.

Geophysical Journal International↗
USGS2024

Using surficial geologic maps, vegetation, and monitoring to address erosion impacts from grazing in Channel Islands National Park, California

Employing a map-unit classification scheme based on geomorphic process and age, the U.S. Geological Survey (USGS) mapped Quaternary surficial deposits of the five islands comprising Channel Islands National Park (CHIS), California, as no such maps previously existed. Mapping was motivated through an agreement with the National Park Service (NPS) to aid natural resource assessments, including post-grazing disturbance recovery and identification of mass wasting and tectonic hazards. The resulting detailed (1:12,000 scale) maps portray areas of upland erosional transport processes and alluvial, fluvial, eolian, beach, marine terrace, mass wasting, and mixed depositional processes at a scale commensurate with the process signal on the landscape. Detailed, GPS-located observations of sedimentology, pedogenic development, landscape position, and geomorphology constrain map unit identification. In addition, map-unit boundary delineation was determined in a GIS using a variety of high-resolution (sub-meter) aerial imagery and lidar-based DEMs complemented by derivative raster products. The GIS geodatabase schema used is GeMS compliant. Resolvable map-unit widths are as small as 10-15 meter. Locational accuracy of unit boundaries is estimated to be as low as 10 m or better. Detailed mapping is at a scale allowing determination of geomorphic transport processes operating on the landscape and conveys geologic information to land managers at a resolution commensurate with park resources. Map unit ages are constrained by relative soil chronostratigraphy, vertebrate and invertebrate fossils, and materials collected for dating control using radiogenic, luminescence, and cosmogenic radionuclide techniques. Associated topical studies examined connections between geology, surface water, and vegetation distribution and restoration.

California↗
USGS2024

Seismic velocity changes from repetitive seismicity at Mauna Loa prior to and during its 2022 eruption

Mauna Loa’s short-lived eruption from late November to early December 2022 marked the culmination of nearly a decade of elevated seismic activity and geodetic inflation. The volcano has been monitored by a network of permanent, short period and broadband seismometers. I used the continuous waveform data from that network starting in 2012 to generate a catalog of seismicity that enhances the US Geological Survey Hawaiian Volcano Observatory’s public seismic catalog with four times the number of earthquakes, which were then grouped by waveform similarity. Analysis of subtle delays in the timing of arrivals of scattered waves between pairs of earthquakes in this catalog yields a history of small changes in the shallow seismic velocity structure of the volcano. Seismic velocities have been shown at other volcanoes to change during unrest and eruption. My results show a decrease in seismic velocity centered on the summit beginning in September 2022, corresponding to the onset of a vigorous precursory swarm of seismic activity and shallow inflation. During the eruption itself, I observe large changes due likely to dike opening along the northeast rift zone and deflation of the summit reservoir. However, seismic velocity changes associated with non-volcanic sources such as ground shaking from large earthquakes and meteorological influences at seasonal and diurnal time scales are also observed, and these dominate the velocity changes prior to the eruption. Proper accounting of these effects will be a requirement for use in real-time monitoring, and this work serves as a starting point in that endeavor for Mauna Loa.

Hawaii↗
USGS2024

Testing megathrust rupture models using tsunami deposits

The 26 January 1700 CE Cascadia subduction zone earthquake ruptured much of the plate boundary and generated a tsunami that deposited sand in coastal marshes from northern California to Vancouver Island. Although the depositional record of tsunami inundation is extensive in some of these marshes, few sites have been investigated in enough detail to map the inland extent of sand deposition and depict variability in tsunami deposit thickness and grain size. We collected 129 cores in marshes of the Salmon River estuary in Oregon and reanalyzed 114 core logs from a 1987–88 study that mapped the inland extent of circa 1700 CE sandy tsunami deposits. The ca. 1700 CE tsunami deposit in the Salmon River estuary is easily recognized in cores ≤1 m deep in which a buried marsh peat is overlain by a well sorted sand bed with a sharp lower contact that thins and fines inland. We use tsunami deposit data and models of sandy tsunami sediment transport (using Delft3D-FLOW) to test 15 rupture models that could represent a ca. 1700 CE earthquake. At least 12–16 m of slip offshore of the Salmon River, which results in 0.8–1.0 m of coastal coseismic subsidence, is required to match the ca. 1700 CE sand deposit's inland extent, which is consistent with models of heterogeneous megathrust slip in ca. 1700 CE. Our methods of detailed tsunami deposit mapping, combined with sediment transport modeling, can be used to test models of megathrust ruptures and their tsunamis to potentially improve earthquake and tsunami hazard assessments.

JGR Earth Surface↗
USGS2024

Seismicity zoning at Coso geothermal field and stress changes from fluid production and migration

The Coso geothermal field is a major geothermal power production site in the western United States. It has been observed that low-magnitude seismic events ( M < 3.71) are unevenly distributed in three distinct zones, namely, nearfield (<3 km), midfield (3–6 km), and farfield (> 6 km) from the Coso geothermal plant. These zones exhibit distinct changes in earthquake location before and during geothermal production episodes that began in 1986. After 1986, the midfield region of the main flank experiences a significantly lower seismicity rate than the surrounding areas before production episodes. During 2014–2019, the farfield earthquakes cluster in the eastern and western parts of the greater Coso area, which is discernably different from how those pre-production earthquake events were distributed along the conjugate NW-SE and SW-NW trending structures across the main flank. Here, we analyze the stage of stress with finite-element-based poroelastic simulations to illustrate how the spatiotemporal evolution of the seismicity is associated with the pattern of stress perturbations caused by fluid migration amid the operations of geothermal power plants. Generally, ∼70% of co-production seismicity is found in zones of increased Coulomb stress between 2014 and 2019 at >99% confidence. Meanwhile, the midfield zone of seismic paucity overlaps with the zone of decreasing pore-fluid pressure. Overall, the results provide a physical explanation of how decadal geothermal operations at Coso have perturbed stress-field changes and contributed to the evolving characteristic seismic pattern, shedding insights into assessing the seismic hazard in other geothermal settings.

California↗
USGS2025

U.S. Geological Survey geomagnetic variometer data: Capitalizing on seismic infrastructure

The U.S. Geological Survey’s Geomagnetism Program is collaborating with the Earthquake Hazards Program and Global Seismographic Network Program to densify magnetic field observations. This collaboration focuses on the installation of magnetometers, or magnetic variometers, at existing seismic stations. Along with improving the density of space weather observations for hazard monitoring, these data can be used to correct colocated magnetic field induced noise in seismic data. Such corrections are especially useful during time periods of large magnetic storms where the magnetic field‐induced instrument noise can be of similar amplitude to earthquake ground‐motion records.

contiguous United States↗
USGS2026

Magnitude conversion relations create substantial differences in seismic hazard models

Earthquake catalogs are essential data inputs for seismic hazard modeling. Because earthquake magnitudes are reported in a variety of types (e.g., local magnitudes and moment magnitudes), magnitude conversion relationships must be used to convert the different magnitude types present in a catalog to a uniform magnitude type to avoid biases in the hazard computation. However, these conversion relationships are often uncertain and have been shown to sometimes perform poorly. Here, we investigate the sensitivity of the gridded seismicity component of the National Seismic Hazard Model (NSHM) to the catalog conversion equations in the Eastern United States. In the 2023 NSHM, magnitudes of various types were converted to moment magnitudes using equations developed by the Central and Eastern United States Seismic Source Characterization for Nuclear Facilities (CEUS‐SSCn), based on least‐squares (LS) regressions made using data from a catalog containing events up through 2008. We recompute these equations using events in the Advanced National Seismic System Comprehensive Earthquake Catalog with multiple magnitudes from 2000 to 2023. Although we prefer the use of orthogonal regressions for our datasets, LS regressions produce broadly similar results, with both approaches exhibiting large deviations from the CEUS‐SSCn conversions, especially at smaller magnitudes. We compare the spatial distribution of annual rates using three different models: (1) the 2023 NSHM conversions, (2) our updated conversions, and (3) no conversions. We find that the choice of conversions leads to substantial differences in the rate forecasts, which can greatly impact the seismic hazard model, particularly in regions with low‐seismicity rates such as the Eastern United States, where the hazard is dominated by gridded seismicity rather than a fault model.

Seismological Research Letters↗
USGS2026

Interplay between tectonics and submarine mass transport deposits in Cortes Basin: New high-resolution geophysics in the Outer California Borderland

The Outer California Borderland (OCB) is an active transform plate boundary offshore Southern California, where the relationship between faulting and submarine mass transport deposits (MTDs) remains poorly understood. Onshore paleoseismic data provide high-resolution earthquake records, whereas marine geophysical data capture longer-term histories. Offshore fault systems pose hazards to infrastructure and dense coastal populations, particularly when linked to submarine landslides. We present new high-resolution geophysical data set (cruise SR2303), including bathymetric and CHIRP sub-bottom data integrated with legacy seismic reflection data and chronostratigraphic constraints from ODP Site 1012 to examine Quaternary MTD recurrence and tectonic controls in the Cortes Basin, OCB. Bathymetry shows deformational features, including slide scarps and previously unmapped fault segments with evidence of Holocene activity. CHIRP profiles reveal 10 stacked MTDs in the East Cortes Basin and 8 in the West Cortes Basin, spanning ∼752 ka with an average recurrence of ∼83.6 ± 1 ka. Acoustic imaging shows 7 MTD intervals coinciding with fault offset increments and fault growth suggesting earthquake-triggered mass wasting. A strong association between MTD occurrences and sea-level extremes also supports glacio-eustatic contribution to slope failure. Stratigraphic correlations suggest quasi-synchronous MTDs across the eastern and western areas, likely triggered by larger eathquakes in the Quaternary. Although the identified MTDs occur relatively far from the Southern California coast, they still pose a potential tsunamigenic hazard requiring further assessment. Moreover, if linked to earthquakes along major strike-slip faults, for example, the Ferrelo fault, the MTDs may provide valuable proxies to constrain rupture scenarios and fault connectivity within the understudied OCB.

California↗
USGS2025

Effect of land cover type on 3D deformation recovery from synthetically deformed high-resolution satellite optical imagery

The limits of detection for earthquake surface deformation in the spatial domain have improved with advances in remote sensing imagery data availability, resolution, and analysis. Sub-pixel correlation and digital elevation model (DEM) differencing from sub-meter, earthquake-spanning satellite optical imagery has enhanced surface rupture mapping and deformation measurements. However, knowledge of measurement accuracy and uncertainty is limited. To address this, we construct orthophotos and digital elevation models (DEMs) from repeat high resolution (∼0.5 m) satellite optical imagery along two sections of the Garlock fault, California with clear fault geomorphology and differing land cover. We deform later sets of DEMs and images with synthetic earthquakes containing both diffuse and discrete horizontal and vertical displacements. Sub-pixel image correlation and DEM differencing demonstrate how vegetation degrades recovered displacement accuracy. In barren land cover, horizontal displacements are detectable to an expected ∼1/10th-pixel size. With shrubs, trees, and grass, detectable displacements increase to >1/2-pixel size, and filtering results by correlation score and using elevation values as input rather than image values improves accuracy. Vertical displacement detection thresholds remain lower in vegetation, at >1-pixel size. Higher slope angles degrade displacement recovery, worsened by vegetation. Diminishing seasonal separation improves accuracy over vegetated regions, though not to the level achieved in barren environments. These results will inform research and operational efforts on the utility of high resolution satellite optical imagery for detecting deformation in varied land cover. Furthermore, they reveal where alternative measurements, such as from LiDAR or radar interferometry, are required to mitigate the effects of vegetation and capture fine-scale crustal deformation.

California↗
USGS2025

A spatiotemporal deep learning approach for predicting daily air-water temperature signal coupling and identification of key watershed physical parameters in a montane watershed

Seasonal shifts from runoff to groundwater dominance influence daily headwater stream temperatures, especially where local groundwater input is strong. This input buffers temperature during hot periods, supporting cold-water habitats. Recent studies use air–water temperature signal metrics to identify zones of strong stream–groundwater connectivity. While Previous studies used air–water signal ratios as proxies for groundwater influence but were limited to specific sites and periods, without dynamic forecasting. This study is the first to forecast daily A r as a spatiotemporal signal using a Graph Convolutional Network–Long Short-Term Memory (GCN-LSTM) model. The model was trained using hydroclimate data (air temperature, precipitation, shortwave radiation, streamflow) and watershed physical features (e.g., sand content, slope). Results showed high predictive skill, achieving R 2 (NSE, RMSE) of 0.86 (0.73, 0.0004) for one-day-ahead to 0.52 (0.50, 0.0009) for seven-days ahead forecasts. Prior studies often have not explicitly incorporated spatial hydrogeologic drivers, but this model explicitly incorporates them to assess their impact on A r forecasting and stream-groundwater connectivity. Feature analysis identified mean sand, elevation, slope, clay, and TWI as key predictors of A r . Stronger groundwater signals appeared in hillslopes, elevations, and tributaries, highlighting watershed influence on streamflow. However, limitations include reliance on historical air–water temperature patterns for training and limited representation of extreme climate conditions. Despite these limitations, unlike previous studies relying on measured in-situ stream and air temperature, this study forecasts A r directly from climate and physiographic features after training, avoiding in-situ data requirements. Findings aiding predictions of stream ecosystem resilience.

New York↗