Sediment-associated processes account for most of the spatial variation in stream ecosystem respiration in the Yakima River basin

Hyporheic zones (HZ) can contribute substantially to total stream ecosystem respiration (ERtot). HZ-focused process-based models may, therefore, effectively predict ERtot across sites, yet this remains untested under variable environmental conditions. Here we evaluate whether spatial variation in HZ respiration predicted via a process-based model explains spatial variation in field-estimates of ERtot across 33 sites in the Yakima River basin in Washington State, USA. We found that HZ respiration predictions did not explain spatial variation in field estimates of ERtot. To investigate further, we partitioned ERtot contributions into water column respiration (ERwc) and sediment-associated respiration (ERsed). ERsed contributed >50% of ERtot at 88% of sites, though relative contributions varied substantially. Despite this dominance, modeled HZ respiration explained neither spatial variation in ERtot nor in ERsed, suggesting that the HZ model alone does not capture the drivers of sediment-associated respiration across these sites. Instead, ERsed spatial variation was primarily explained by gross primary production, stream slope, velocity, and total dissolved nitrogen rather than median grain size, a primary control of HZ respiration predicted by the process-based model. Consistent with recent studies, our results indicate that improving basin-scale ERtot predictions requires integrating hydrologic and biogeochemical processes across hyporheic, benthic, and water column zones.