Groundwater discharge generates streamflow and influences stream thermal regimes. However, the water quality and thermal buffering capacity of groundwater depends on the aquifer source-depth. Here, we pair multi-year air and stream temperature signals to categorize 1729 sites across the continental United States as having major dam influence, shallow or deep groundwater signatures, or lack of pronounced groundwater (atmospheric) signatures. Approximately 40% of non-dam stream sites have substantial groundwater contributions as indicated by characteristic paired air and stream temperature signal metrics. Streams with shallow groundwater signatures account for half of all groundwater signature sites and show reduced baseflow and a higher proportion of warming trends compared to sites with deep groundwater signatures. These findings align with theory that shallow groundwater is more vulnerable to temperature increase and depletion. Streams with atmospheric signatures tend to drain watersheds with low slope and greater human disturbance, indicating reduced stream-groundwater connectivity in populated valley settings.

Human activities and climate change threaten coldwater organisms in freshwater ecosystems by causing rivers and streams to warm, increasing the intensity and frequency of warm temperature events, and reducing thermal heterogeneity. Cold‐water refuges are discrete patches of relatively cool water that are used by coldwater organisms for thermal relief and short‐term survival. Globally, cohesive management approaches are needed that consider interlinked physical, biological, and social factors of cold‐water refuges. We review current understanding of cold‐water refuges, identify gaps between science and management, and evaluate policies aimed at protecting thermally sensitive species. Existing policies include designating cold‐water habitats, restricting fishing during warm periods, and implementing threshold temperature standards or guidelines. However, these policies are rare and uncoordinated across spatial scales and often do not consider input from Indigenous peoples. We propose that cold‐water refuges be managed as distinct operational landscape units, which provide a social and ecological context that is relevant at the watershed scale. These operational landscape units provide the foundation for an integrated framework that links science and management by (1) mapping and characterizing cold‐water refuges to prioritize management and conservation actions, (2) leveraging existing and new policies, (3) improving coordination across jurisdictions, and (4) implementing adaptive management practices across scales. Our findings show that while there are many opportunities for scientific advancement, the current state of the sciences is sufficient to inform policy and management. Our proposed framework provides a path forward for managing and protecting cold‐water refuges using existing and new policies to protect coldwater organisms in the face of global change.

Headwater stream responses to climate change will depend in part on groundwater‐surface water exchanges. We used linear modeling techniques to partition likely effects of shallow groundwater seepage and air temperature on stream temperatures for 79 sites in nine focal watersheds using hourly air and water temperature measurements collected during summer months from 2012 to 2015 in Shenandoah National Park, Virginia, USA. Shallow groundwater effects exhibited more variation within watersheds than between them, indicating the importance of reach‐scale assessments and the limited capacity to extrapolate upstream groundwater influences from downstream measurements. Boosted regression tree (BRT) models revealed intricate interactions among geomorphological landform features (stream slope, elevation, network length, contributing area, and channel confinement) and seasonal precipitation patterns (winter, spring, and summer months) that together were robust predictors of spatial and temporal variation in groundwater influence on stream temperatures. The final BRT model performed well for training data and cross‐validated samples (correlation = 0.984 and 0.760, respectively). Geomorphological and precipitation predictors of groundwater influence varied in their importance between watersheds, suggesting differences in spatial and temporal controls of recharge dynamics and the depth of the groundwater source. We demonstrate an application of the final BRT model to predict groundwater effects from landform and precipitation covariates at 1075 new sites distributed at 100 m increments within focal watersheds. Our study provides a framework to estimate effects of groundwater seepage on stream temperature in unsampled locations. We discuss applications for climate change research to account for groundwater‐surface water interactions when projecting future thermal thresholds for stream biota.