This work was also the first to comprehensively demonstrate wildfire-associated changes in DOC character (by measuring HPO %, UV254, SUVA, FI, and FEEMs) and related DBP-FPs, at the watershed-scale and over multiple flow regimes. The disturbance impacts indicated by all of these quantitative, DOC-associated metrics were all statistically significant, except for FI. Qualitative FEEM results were consistent with these significant shifts. Notably, despite the continued development and promotion of various proxy indicators, UV254 offered the most precise linear correlation with THM-FP, with a coefficient of determination (R2) of 0.6 (in contrast to values of 0.47, 0.42, and 0.39 for DOC, SUVA, and HPO %). Thus, changes in the proxy indicators were related to changes in THM-FP; however, they could not adequately explain response variability, thereby demonstrating the need to 1) better understand relationships between disturbance-associated changes in DOC and their implications to DOC reactivity and 2) advance modeling approaches for describing these relationships. While the mass of various DOC fractions obtained using LC-OCD and HAA-FPs was not analyzed in this manner because of the limited size of the data sets, similar relationships were suggested. Overall, these data suggest that severe wildfire may lead to significant DOC-associated drinking water treatability challenges and that post-fire salvage logging may further exacerbate them-notably, UV254 is unequivocally the best available tool for monitoring these potential impacts at present. THM-FP is generally understood to be linearly correlated with aromatic NOM as measured by UV254 and/or SUVA. In Phase 4, simple strategies for enhancing the prediction of THM-FPs using NOM-associated proxy indicators were investigated. Specifically, the relationship between NOM aromaticity (HPO %, HS, UV254, and SUVA) and THM-FP was examined. Then, HPO and HS were re-analyzed after weighting by mass (DOC concentration)-this appreciably enhanced their prediction performance. This improvement was especially evident for HS, for which the coefficients of determination (R2) increased from 0.10 and 0.26, to 0.85 and 0.88 (Phase 2 and 3 data, respectively). Thus, data processing and reporting are critical to anticipating NOM reactivity; absolute quantities have superior prediction performance. Notably, regardless of these improvements, the relationships between DBP-FP and NOM proxy indicators can be quite variable spatially and temporally, and frequently site specific. More work is required to link source water quality to DBP-FP and drinking water treatability more broadly. In Phase 5, changes in DOC concentration and character and their relationships to regulated DBP-FPs were comprehensively characterized using multiple NOM characterization techniques in the two years during and immediately after forest harvesting in the eastern slopes of the Rocky Mountains in south-western Alberta. Several NOM fractions also were characterized by LC-OCD to inform the relative potential for membrane fouling and microbial regrowth in distribution systems. Like Phase 3, this work was conducted as part of the ongoing SRWP in which two watersheds that served as unburned-reference watersheds in Phase 3 were studied. They were fully calibrated for climate, streamflow, and water quality for 11 years [2004-2014]). Three sub-watersheds within one watershed were harvested using clear-cut with patch retention, strip-shelterwood cut, and partial cut. All possible best management practices (BMPs) were followed to minimize disturbance impacts on water quality. Samples were collected during the dominant regional streamflow regimes. Notably, no substantial impacts of forest harvesting on water quality and treatability were observed during the harvest and first post-harvest years. Thus, this work suggests that forest harvesting with careful implementation of BMPs for erosion control may mitigate the potentially catastrophic impacts of wildfire on drinking water treatability without significantly compromising it.