This resource contains the data and scripts used for:
Goeking, S. A. and D. G. Tarboton, (2022). Variable streamflow response to forest disturbance in the western US: A large-sample hydrology approach. Water Resources Research, 58, e2021WR031575. https://doi.org/10.1029/2021WR031575.

Abstract from the paper:
Forest cover and streamflow are generally expected to vary inversely because reduced forest cover typically leads to less transpiration and interception. However, recent studies in the western US have found no change or even decreased streamflow following forest disturbance due to drought and insect epidemics. We investigated streamflow response to forest cover change using hydrologic, climatic, and forest data for 159 watersheds in the western US from the CAMELS dataset for the period 2000-2019. Forest change and disturbance were quantified in terms of net tree growth (total growth volume minus mortality volume) and mean annual mortality rates, respectively, from the US Forest Service’s Forest Inventory and Analysis database. Annual streamflow was analyzed using multiple methods: Mann-Kendall trend analysis, time trend analysis to quantify change not attributable to annual precipitation and temperature, and multiple regression to quantify contributions of climate, mortality, and aridity. Many watersheds exhibited decreased annual streamflow even as forest cover decreased. Time trend analysis identified decreased streamflow not attributable to precipitation and temperature changes in many disturbed watersheds, yet streamflow change was not consistently related to disturbance, suggesting drivers other than disturbance, precipitation, and temperature. Multiple regression analysis indicated that although change in streamflow is significantly related to tree mortality, the direction of this effect depends on aridity. Specifically, forest disturbances in wet, energy-limited watersheds (i.e., where annual potential evapotranspiration is less than annual precipitation) tended to increase streamflow, while post-disturbance streamflow more frequently decreased in dry water-limited watersheds (where the potential evapotranspiration to precipitation ratio exceeds 2.35).

The following scripts (R language and environment for statistical computing) produce the results, figures, and tables in this paper (in the order in which they appear in the paper; requires either running data compilation/aggregation scripts first OR using provided data files watersheds.csv and wb_filtered.csv):
1. Map_watersheds.r
2. Analysis_M-K_trend_test.r
3. analysis_M-K_quadrant_figure.r
4. analysis_timetrend_linear.r
5. analysis_regressn_w-veg.r

The following scripts (R) compile the data, aggregated from other sources prior to the analyses in the scripts listed above:
1. compilation_CAMELS.r
2. compilation_Daymet.r
3. compilation_USGS.r
4. compilation_FIA.r
5. compilation_CAMELS_Daymet_USGS.r (must run scripts #1-3 first)
6. watershed_compilation.r (must run scripts #1-5 first)