Data for Variable streamflow response to forest disturbance in the western US: A large-sample hydrology approach

This resource contains the data and scripts used for:
Goeking, S. A. and D. G. Tarboton, (2021), \"Variable streamflow response to forest disturbance in the western US: A large-sample hydrology approach,\" Submitted to Water Resources Research, Earth and Space Science Open Archive: 32, http://doi.org/10.1002/essoar.10508683.1.

Abstract from the paper:
Forest cover and streamflow are generally expected to vary inversely because reduced forest cover typically leads to less transpiration. 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 evidence for cases where forest cover loss leads to decreased streamflow using hydrologic, climatic, and forest data for 159 watersheds in the western US from the CAMELS dataset. Forest change and disturbance were quantified in terms of net tree growth (total growth volume minus mortality volume) and mean annual rate of tree mortality, respectively, from the US Forest Service’s Forest Inventory and Analysis database. Annual water budget components were analyzed using multiple methods: Mann-Kendall trend analysis, time trend analysis to quantify change not attributable to precipitation and temperature, and multiple regression. 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 that factors other than disturbance, precipitation, and temperature are driving streamflow changes. Finally, multiple regression analysis indicated that although change in streamflow is positively related to tree mortality, the direction of this effect is dependent upon aridity. Specifically, forest disturbances in wet, energy-limited watersheds (i.e., where potential evapotranspiration is less than 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)