Effects of forest disturbance on water yield and peak flow in low‐relief glaciated catchments assessed with Bayesian parameter estimation

Empirical assessment of how forest disturbance affects streamflow has been traditionally limited to small, experimental catchments. However, larger catchments where landscape management occurs have emergent drivers of streamflow at scale, and thus may exhibit novel responses to land cover disturbance. We used statistical models of water yield and annual maximum peak streamflow for multiple large (> 50 km2) forested catchments in the low‐relief glaciated region of central North America to investigate how forest disturbance and climatic variability affect water yield and peak flows in similar landscapes. We used linear models, linear mixed effects models, and probabilistic flood‐frequency analysis with Bayesian parameter estimation in two case studies. These included: (1) a wildfire that burned ~30% of the 650 km2 wilderness Upper Kawishiwi catchment, and (2) 11 catchments within the St. Louis River Basin ranging from 56 to 8880 km2 with a patchwork disturbance regime wherein ~0.25% to 1% of the catchment is harvested or converted to non‐forest land use each year. We also assessed the most likely hydrological recovery year after forest disturbance, and the relative importance of stationary and nonstationary drivers of streamflow. We found forest disturbance correlated with declines in water yield for low‐level disturbance regimes in some catchments, but that water yield increased in response to the large‐scale wildfire. Positive and negative associations of forest disturbance with peak flows were observed. Hydrologic recovery time ranged from 5 to 13 years for water yield and peak flows following disturbance. Despite these effects of forest disturbance on streamflow, effects of climatic variability and stationary catchment size factors were more prominent streamflow drivers. Basins larger than ~50 km2 in low‐relief glaciated regions can be impacted by forest cover change even on