Event Scale Analysis of Streamflow Response to Event Scale Analysis of Streamflow Response to Wildfire in Oregon, 2020 Wildfire in Oregon, 2020

: Wildﬁre increases the magnitude of runoff in catchments, leading to the degradation of ecosystems, risk to infrastructure, and loss of life. The Labor Day Fires of 2020 provided an opportunity to compare multiple large and severe wildﬁres with the objective of determining potential changes to hydrologic processes in Oregon Cascades watersheds. Geographic information systems (GIS) were implemented to determine the total percentage burned and percentage of high burn severity class of six watersheds on the west slope of the Oregon Cascade Range. In addition, two control watersheds were included to contrast the inﬂuence of climatic effects. Spatial arrangements of burned patches were investigated for correlation to streamﬂow response by utilizing landscape metrics algorithms, including Largest Patch Index (LPI), mean gyration (GYRATE), Contiguity Index (CONTIG), Patch Cohesion Index (COHESION), and Clumpiness Index (CLUMPY). Results of the ﬁrst-year post-ﬁre response were consistent with other studies of ﬁre effects in the Paciﬁc Northwest (PNW) and indicated changes to runoff dynamics were difﬁcult to detect with inferential statistics, but the largest changes in runoff coefﬁcients occurred in watersheds having the greatest percentage burned. Correlation analysis indicated relationships between event runoff coefﬁcients and percentage burned during the 2020 ﬁre season. Control watersheds show confounding runoff coefﬁcients, point to the inﬂuence of ongoing drought, and complicate conclusions about the role of spatial burn severity patterns. These results could guide future post-ﬁre studies of spatial patterns of burn severity and could assist watershed managers to prioritize at-risk PNW catchments to minimize harm to ecological and societal values.


Introduction
In the western United States, over 60% of municipal water supplies originate in forested areas, and wildfire reduces the water resiliency of communities that rely on these sources [1].Stresses to water resources are increasingly common in burned watersheds due to changes to runoff dynamics that result from enhanced surface runoff and soil erosion [2].Implications of these shifts include changes in annual water yield [3] and degradation of water quality [4,5], which challenges the flexibility of socio-economic systems to adapt.Additionally, severely burned watersheds can have a reduced capacity to hold water in the soil, which increases the occurrence of extreme hydrological events such as floods and landslides, threatening expanding populations near the urban-wildland interface and leading to loss of life and property [6].
Classes of burn severity are determined by the assessment of remotely sensed imagery and/or field reconnaissance of watershed disturbances that are driven by fire intensity, the physical release of energy from combusting organic material over a period of time, which can cause a gradient of alterations to biomass and soil properties occurring during a wildfire event [7,8].Measures of burn severity have been shown to be correlated with increased runoff in plot-scale studies and attributed to alterations of soil and vegetative cover [9,10].First, the removal of vegetative canopy and forest litter increases throughfall to the soil surface [11].Secondly, alterations of soil structure include the reduction in aggregation and water repellency, which can drastically reduce the hydraulic conductivity of soils [12].Fire disturbance is dependent on the fire intensity, duration [13], and soil properties [14] and varies spatially [15].These disturbances can be categorically mapped in a landscape using remotely sensed imagery including Landsat data.Landsat data can be adapted to represent burn severity by first developing a Normalized Difference Vegetation Index (NDVI), which can be transformed into the Normalized Burn Ratio (NBR) to represent burn severity [16].These map coverages are commonly used as proxies for the level of consumption of ground cover and alteration of the soil, assisting land managers to prioritize responses in burned areas, and have been increasingly used in fire research [17].
Streamflow metrics can quantify the influence of fire severity on hydrological response at the catchment and watershed scale.Higher burn severity often results in increased runoff and erosion rates on catchment hillslopes due to a reduction in soil infiltration [10].At the watershed scale, these alterations can lead to changes in streamflow characteristics, such as peak flow, timing of runoff, base flow, and water yield [18].While comparing burned and unburned catchments using the unit peak discharge method, [6] found that peak flow is more sensitive to wildfire disturbance than annual water yield.Short-term increases in base flow have been observed because of the transitory reduction of evapotranspiration and increased groundwater availability [19,20].Timing of runoff has been shown to be correlated with percentage of burned area within watersheds [21,22].
The influence of burn patch connectivity on hydrologic response has been little studied, but it has been proposed that there is a link between the distribution and location of fire disturbances in a basin and increased runoff delivery to stream channels [23].Additionally, the spatial heterogeneity of burn severity can form a patchy mosaic, further increasing the spatial variability of soil infiltration rates, which enhances overland flow in some areas [24] or acts as a sink for runoff [25,26].Preferential flow paths have been described with the term functional connectivity, which describes the tendency of landscape patterns to regulate the movement of overland flow [27].Thereby, the hydrological response of catchments to precipitation is a balance between functional connectivity of runoff processes and storage [28].In this sense, the removal of ground cover and alteration of hydraulic soil properties can further enhance the functional connectivity in areas of high disturbance [29].Conversely, normal infiltration rates remain in unburned areas within the fire perimeter, effectively reducing connectivity and increasing soil storage.The destabilization of this balance can be rapid and highly variable due to the spatial heterogeneity of precipitation and catchment characteristics [30].Additional variability in fire extent, severity, and overlapping disturbances cause difficulties in comparing the hydrological response of different watersheds to fire disturbance [11].
Understanding how the spatial distribution of landcover might influence physical processes between different locations is a primary interest of geographers.One method of describing these patterns is the use of landscape metrics-algorithms that allow the quantification of landcover pattern and arrangement [31].Additionally, landscape metrics can quantify subjective descriptions of landcover aggregation and density observed in spatial, categorical data to produce comparable indices.There are many metrics that can be calculated, but often, metrics describe similar aspects of patch characteristics.Furthermore, it has been suggested that metrics describing patch compaction, the number of patches, contagion, texture, patch shape, and average patch composition are adequate in describing the spatial arrangement of most landscape patterns [32].Fire ecologists have modified this rationale to fit project objectives, having described significant relationships between forest composition and burn severity patterns [33,34].Landscape metrics have been used to quantify spatial patterns in the post-fire environment representing two drivers of streamflow response: vegetation mortality [35] and soil hydrophobicity [36].The configuration of burn patches has been quantified with landscape metrics across modeled hillslopes, and the spatial structure of vegetative patterns has been found to be significant in describing runoff and discharge [37,38].
In urban environments, the spatial arrangement of landcover on hydrologic response has been well studied.For example, the fragmentation of impervious urban surfaces was found to be negatively correlated with flashiness and streamflow variability in North Carolina and Utah, with large, continuous urban patches producing flashier streamflow responses [39].Landscape metrics have been utilized in quantifying the effect of urban land cover change on aquatic health, indicating a correlation of spatial patterns with water quality and degradation at accelerating rates [40].Similarly, landscape fragmentation metrics have been employed to examine changes in total suspended solids in an urbanizing Oregon watershed [41].The use of landscape metrics in determining runoff response in forested catchments is less common, but modeling efforts directed at the influence of hydrological connectivity on runoff response demonstrate that landscape metrics are effective as a proxy for heterogeneity that influences runoff [42].Studies using experimental plots in burned areas indicate that unburnt patches effectively reduce hydrologic connectivity and produce less post-fire runoff [26].Additionally, modeling experiments have tested the ability of patch aggregation metrics such as the Contiguity Index (CONTAG) to represent hydrological connectivity in an idealized snow-dominated catchment in Idaho, USA, finding spatial patterns of burned and unburned patches that are important in explaining variations in discharge and soil condition, with less aggregated burn mosaics showing a greater departure from normal [37].Limited investigation at increased scales shows that there is a need to examine the ability of landscape metrics to represent spatial patterns of wildfire and possible effects on watershed runoff, with the caveat that landscape patterns do not always represent physical hydrological processes [43].
Specifically, this work seeks to answer the following questions: 1.
To what degree does the arrangement of burn severity patches influence streamflow response at event scales? 2.
Are changes to event runoff ratios more detectable using daily values or using instantaneous values for individual events?

Context of 2020 Labor Day Fires
The 2020 fire season in western Oregon was remarkably severe and resulted from extreme atmospheric aridity that coincided with a strong east to north-east katabatic wind event, which drove the rapid expansion of existing lightning ignitions and human-caused fires [44].Between September 7 and September 15, these strong down-slope winds drove multiple large wildfires into populated areas, forcing the evacuation of about 40,000 people and placing almost 500,000 people under evacuation readiness notices [45].Ultimately, 400,000 hectares burned across Oregon in 2020, representing the most substantial fire season in about 120 years.Despite causing extensive harm to life and property, the 2020 fires were not unusual for Cascade west-slope mesic forests.Several large fire events, greater than 300,000 hectares, have occurred post-European settlement, which include the Great Fire of 1845 and Yacolt Burn of 1902 [46].The most unprecedented aspect of the Labor Day Fires of 2020 was the population centers that were impacted, both directly through loss of socioeconomic values and indirectly due to cumulative threats of poor air quality and interrupted water delivery.It is unclear what implications these fires will have on the forest ecology and what the ongoing impacts will be for communities that may suffer degraded water sources.The Oregon Labor Day Fires provide a unique opportunity to compare the hydrologic response of multiple watersheds in a relatively small geographic area with the aim of determining the effect of burn severity patterns on the rate and timing of rainfall-runoff at the event scale.

Study Area Overview
The study area is located in the Western Cascades ecoregion, with two distinct geologic provinces: the High Cascade and the Western Cascade (Figure 1).These provinces have contrasting influences on the hydrology of the watersheds that drain them due to the unique underlying geology [47].The High Cascade province contains Quaternary and Pliocene volcanics, related to the subduction of the Juan de Fuca plate.These areas have primarily produced basaltic and andesitic lava flows from young shield and cinder volcanos, resulting in a broad north to south running plateau.These volcanic terrains provide strong subsurface pathways for groundwater interaction in the headwaters of west-slope Cascade watersheds.At lower elevations, the geology is dominated by Tertiary aged rocks of the Western Cascade Province, which have much lower permeability and provide greater runoff rates (Figure 2).The dominance of the High Cascades province in the upper portion of the watersheds contributes to the low intra-annual variability of watersheds in the Oregon cascades [52].Most precipitation falls from November to May, and the amount and type are dependent on the elevation of the watershed, with lower elevation watersheds receiving a rain-snow mix and higher elevations accumulating dense snowpacks [53].Many watersheds in this region are utilized for timber extraction, which occurs on private, federal, and state-managed lands.It is worth mentioning that extensive salvage logging, the commercial harvesting of burned trees, was undertaken shortly after the fire, primarily concentrated on state managed lands, but specific information about the extent was not available at the time of study and not accounted for herein.

Individual Watershed Characteristics
To investigate the streamflow response to burn severity patterns, six watersheds in the Oregon Cascades that burned in the 2020 Labor Day fires were selected (Table 1, Figure 1, Figure 3).Criteria for selection included having greater than 5% burned area within the watershed boundary and containing a complete discharge record of one pre-fire and post-fire year.These watersheds were also limited to having either unregulated streamflow or minimum impacts by upstream reservoirs to reduce the influence of anthropogenic alterations.Additionally, two watersheds that were not burned by the fires were selected as control watersheds (Table 1, Figure 1, Figure 3).These control watersheds were selected based on having similar watershed areas, low flow alteration, and complete discharge records for the pre-and post-fire periods.The hydrology of these areas is characterized by high rates of erosion and sedimentation, punctuated by frequent mass wasting events [48].In addition, the hydrology of the region is heavily influenced by the percentage of permeable High Cascade terrain within the contributing area, where watersheds with greater mean elevations translate to greater mean permeability and lower elevation watersheds to low mean permeability (Table 1).The permeable geology obtained from USGS StreamStats was originally produced from a 1:2,500,000 map considering the percentage of the basin surface area that contained high permeability geologic units, which were Pliocene volcanic rocks, Upper Tertiary andesite, and Quaternary volcanic rocks [49].For soil permeability, we used watershed-wide average soil permeability (SOILPERM) in mm per hour, estimated with a data layer created from a 1:250,000 STATSGO map [50].Permeability values are assumed to be uniform across watersheds (Table 1), although the permeability of these formations can vary spatially due to local hydrothermal alterations and depth of soil formation [51].
The dominance of the High Cascades province in the upper portion of the watersheds contributes to the low intra-annual variability of watersheds in the Oregon cascades [52].Most precipitation falls from November to May, and the amount and type are dependent on the elevation of the watershed, with lower elevation watersheds receiving a rain-snow mix and higher elevations accumulating dense snowpacks [53].Many watersheds in this region are utilized for timber extraction, which occurs on private, federal, and state-managed lands.It is worth mentioning that extensive salvage logging, the commercial harvesting of burned trees, was undertaken shortly after the fire, primarily concentrated on state managed lands, but specific information about the extent was not available at the time of study and not accounted for herein.

Individual Watershed Characteristics
To investigate the streamflow response to burn severity patterns, six watersheds in the Oregon Cascades that burned in the 2020 Labor Day Fires were selected (Table 1, Figures 1 and 3).Criteria for selection included having greater than 5% burned area within the watershed boundary and containing a complete discharge record of one pre-fire and post-fire year.These watersheds were also limited to having either unregulated streamflow or minimum impacts by upstream reservoirs to reduce the influence of anthropogenic alterations.Additionally, two watersheds that were not burned by the fires were selected as control watersheds (Table 1, Figures 1 and 3).These control watersheds were selected based on having similar watershed areas, low flow alteration, and complete discharge records for the pre-and post-fire periods.Oregon study area including the location USGS streamgaging stations, watershed bo aries, and burn severity classes.Low burn severity includes areas where surface organics have been completely combusted, roots are intact, and some patches of grey ash are present.Mode burn severity includes areas where some surface organics and roots have been consumed, s needles and leaves may remain, and patches of grey/black ash are present.High burn severit cludes areas where almost all the organic groundcover has been destroyed, many surface root consumed, bare soil is exposed and soil structure has been disturbed, stumps and trees are c pletely burned out in places, and thick deposits of white and grey ash are present [54].

Figure 3.
Oregon study area including the location USGS streamgaging stations, watershed boundaries, and burn severity classes.Low burn severity includes areas where surface organics have not been completely combusted, roots are intact, and some patches of grey ash are present.Moderate burn severity includes areas where some surface organics and roots have been consumed, some needles and leaves may remain, and patches of grey/black ash are present.High burn severity includes areas where almost all the organic groundcover has been destroyed, many surface roots are consumed, bare soil is exposed and soil structure has been disturbed, stumps and trees are completely burned out in places, and thick deposits of white and grey ash are present [54].Beginning in the southern portion of the study area, the North Umpqua River at Winchester, OR (USGS site 14319500) has the largest contributing area (1005 km 2 ) and receives 19% less annual precipitation than the Clackamas watershed in the north.Overall, the North Umpqua has the third highest mean elevation (1005 m) and contains about 26% permeable geology (Table 1).The North Umpqua is primarily a snow-driven watershed that can respond rapidly to rain on snow events despite the porous High Cascades geology that facilitates ground water inputs.The Winchester Dam is located just upstream of the gaging station, one of eight dams in the watershed, but all dams have a small storage capacity and limited influence on the streamflow of the North Umpqua River.The North Umpqua was impacted by the Archie Creek fire, which originated 20 miles upstream of Glide, OR and spread to 53,233 hectares.The USGS gage North Umpqua River near Winchester sustained 12.9% of its drainage area burned, with 4.5% burning at high severity (Table 2).Due to fragmented landcover, the burn severity patterns were highly influenced by the availability of fuels in lands dissected by timber parcels [56].The northwestern portion of the watershed contains densely clustered high burn severity patches primarily located in the Rock Creek drainage.The North Santiam watershed drains 558 km 2 above the USGS streamgage (North Santiam River below Boulder Creek near Detroit, OR) before entering the Detroit Lake reservoir.The geology of the upper North Santiam watershed is typical of other westslope watersheds in Oregon, having soils dominated by easily mobilized, unconsolidated volcanic parent material containing abundant smectite clays and kaolinite deposits [57].The underlying geology is similar to the above-mentioned description of the North Umpqua but contains a greater proportion of High Cascade province than Western Cascade province, resulting in a much higher percentage of permeable geology (Table 1).Additionally, the hydrology is primarily influenced by snow accumulation in high elevation areas and rain on snow events that produce the largest peak flows.The northwestern section above the USGS streamgage was impacted the Lionshead fire, and 24.9% was within the fire perimeter, although much of the upper watershed was not disturbed.The investigation of the burn severity classes within the drainage area shows a majority low burn severity with only 1.7% high burn severity (Table 2).Areas of high burn severity are concentrated in the Marys Creek and Boulder Creek drainages near the streamgaging location.
The Breitenbush River is a subbasin in the North Santiam watershed that drains 274 km 2 into the Detroit Reservoir near Detroit, OR (Table 1).The geology of the Breiten-bush river differs slightly from nearby watersheds as it is dominated by the Little Butte Terrain containing late Oligocene to early Miocene volcanics with highly varying lithologies.These ignimbrites and breccia are the oldest exposed rocks in the north central cascades resulting from pyroclastic density currents [58].The upper watershed contains rocks typical of the Western Cascade and High Cascade terrains.The Breitenbush watershed was heavily burned (79%) by the Lionshead fire in the mid to upper watershed and, to a lesser extent, the Beachie Creek fire in the lower watershed (Table 2).
The Little North Santiam River flows into the North Santiam downstream and west of the Breitenbush River.Although the neighboring watersheds have moderate permeabilities, the Little North Santiam watershed has the lowest permeability among the burned watersheds (Table 1).Additionally, high concentrations of suspended sediment in the North Santiam River often originate from the Little North Santiam due to aquitard forming, clay-rich soils, which cause frequent mass movement events [57].The geology is dominated by Oligocene and Miocene volcanics at the valley bottom and is capped by younger andesitic rocks of the Late High Cascade Volcanic event.This follows the same trend of less permeable geologies at low elevations and highly permeable geologies at high elevation.Annual precipitation averages 2235 mm, which is the highest of the study watersheds.The largest percentages of high severity burn class and burned area are found in the Little North Santiam River (Table 2).
In the northernmost portion of the area of interest, the upper Clackamas River (USGS site 14209500, Clackamas River above Three Lynx Creek) has a contributing area of 1266 km 2 (Table 1).The geology of the Clackamas River watershed is consistent with the descriptions of other west-slope Cascade watersheds, containing proportional amounts of High Cascade and Western Cascade provinces.The contributing area was affected by the Riverside Fire, burning approximately 11.3% of the watershed.The watershed holds the lowest percentage of high soil burn severity class out of the six burned watersheds (Table 2).
The Molalla River has only one active stream gage near the town of Canby, OR, (14200000) with an upstream area of 841 km 2 (Table 1).Almost all of the headwaters of the Molalla River originate in less permeable Western Cascade terrain, and the low elevation areas are comprised of a mixture of course and fine-grained fluvial deposits from the Missoula Floods [59].Most of the burned areas are found in the headwaters of the Molalla River (Table 2).
Two unburned watersheds were selected as controls to check the influence of climate on streamflow.The first control was the South Santiam River, which has a contributing area 451 km 2 above the USGS streamgage (14185000) and flows into the Foster Reservoir near Sweet Home, OR.Approximately 5% of the watershed area is comprised of permeable geology, with the majority being Western Cascade terrain.Second, Thomas Creek was also unburned by the Labor Day Fires of 2020 and has the lowest mean elevation of all sites at 557 m.Thomas Creek has the lowest percentage of permeable geology and the lowest average soil permeability (Table 1).

Precipitation and Flow
Daily precipitation time-series from the PRISM Climate Group [60] were retrieved for watershed centroid coordinates using the PRISM R package [61].PRISM time-series for single locations were preferred over precipitation observations due to the limited number and inconsistent distribution of meteorological stations in the study area.Daily and instantaneous discharge values were obtained from National Water Information System (NWIS) using the R package dataRetrieval [62].

Streamflow Metrics
At the event scale, two streamflow metrics were selected to capture potential changes to the rainfall-runoff pathways caused by varying spatial arrangements of patch size, aggregation, and fragmentation.The first was the event runoff ratio (Q r ), which is the quotient of runoff depth to precipitation amount or the amount of precipitation that is transferred to the catchment outlet as stormflow.This metric was selected based on the findings of other studies, including [63], who found that runoff ratio increased most at the event scale compared with daily and monthly time steps.Some have found runoff coefficients to be most altered for individual events as opposed to longer time scales, noting that preceding soiling moisture influenced some storms [64].
A second runoff metric, event peak runoff ratio (Q p ), was selected to account for the influence of antecedent soil moisture conditions on surface runoff.Several studies have used the metric for the long-term comparison of runoff ratios for annual peak flow events [65,66].Here, it is adapted to any runoff event to aid intra-storm comparison without the compounding effects of antecedent soil conditions.In the context of individual storms, the event peak runoff ratio is defined as the net storm flow divided by the 3 day precipitation amount leading up to the peak flow event.Net storm flow is the difference between the peak-flow and minimum daily average discharge during the 3 day period.The peak runoff ratio is then normalized to the watershed area to aid comparison between catchments of different sizes.The benefit of this method is that soils can be considered for rainfall-runoff events, where runoff can be highly dependent on antecedent moisture conditions-an important factor in explaining non-linear streamflow response [66,67].The event peak runoff ratio was included in the analysis because the metric can account for antecedent soil conditions preceding a precipitation event.Considering the influence of soil properties is an important part of any comparative watershed study, it was expected that this metric would be paramount for parsing fire effects from drought impacts.
It was hypothesized that Q r and Q p would increase most in watersheds with a greater coverage of high burn severity patches and show little to no change for the control watersheds.Because high burn severity patches can be associated with hydrophobic soil conditions, reduced rainfall interception, and increased hillslope connectivity, these areas can influence rainfall-runoff pathways, leading to a more efficient transfer of precipitation to the catchment outlet.

Landscape Metrics
The R package Landscape Metrics [68] was implemented to calculate landscape metrics for the Largest Patch Index (LPI), Clumpiness Index (CLUMPY), mean gyration (GYRATE), Contiguity Index (CONTIG), and Cohesion Index (COHESION).Metrics were selected out of the suite of shape metrics compiled by [31] that pertain to landscape aggregation, complexity, and density.Some caution during the selection of the landscape metrics was considered because landscape patterns do not always equate to landscape processes [42].Thus, landscape metrics were carefully evaluated for their relevance to the hypothesis that the level of aggregation and continuity of severe burn patches explain some variability of post-fire streamflow response.Using this logic, CLUMPY was selected due to linear behavior when measuring class aggregation and complexity without problems of autocorrelation between class areas [69].CLUMPY measures the proportion of like adjacencies that match a landscape with a random distribution.It ranges from −1 (fragmented) to 1 (aggregated), with 0 representing complete spatial randomness.The high burn severity class area was assessed using LPI, which has linear behavior to allow for simple interpretation.GYRATE was included because it can be considered a "shape" index and provides information about the elongation of a patch.Burned areas with higher levels of GYRATION would translate to large and elongated burn patches that could provide efficient runoff pathways.Small values of GYRATION would result in discontinuous burned patches, which would be more contained by other burn classes that could potentially buffer runoff.CONTIG measures the connectivity between classes from 0 (less connected) to 1 (highly connected) using a 9-by-9 moving window focal filter to determine the mean of a class belonging to that class.Finally, COHESION offers information about the dispersion of burn classes, with 0 being the most isolated and 100 being the most aggregated.It is important to note that landscape metrics applied are an attempt to represent the pattern and arrangement of burned patches and were not intended to quantify the hydrological connectivity of the burned area because the suite of metrics were not designed within a hydrologic context.Table 3 summarizes the values of the total burned area and the selected landscape metrics for the high severity patches in each burned watershed.Often, land and water resource managers are concerned with the immediate response to wildfire, which is challenged by uncertainty in unique watershed attributes that produce variable responses to disturbance.To approach this problem, burn severity maps were quantified using landscape metrics to determine if there were differences in burned area characteristics between locations.Then, event-based runoff metrics were produced to determine if there was a significant difference between pre-and post-fire responses to precipitation events.
Post-fire runoff events were investigated for changes using hydrologic data obtained at the daily resolution for water year 2021.Pre-fire years were compared based on precipitation totals and drought indices, with the period of interest set as 2010 to 2020 to reduce influence from long-term climatic shifts.Extreme drought conditions during the 2021 water year challenged the selection of a similar pre-fire year.Ultimately, water year 2015 was selected as the most similar pre-fire year, receiving slightly greater precipitation totals in the Oregon Cascades.
The event scale was further investigated using fine-resolution hydrologic data for individual precipitation events following wildfire.Due to the well-defined seasonality of the Oregon Cascades, two distinct periods were investigated following wildfire disturbance.The early-wet season was defined as the first major rainfall-runoff event occurring from September through November-a period which is characterized as having low antecedent soil moisture conditions.The mid-wet season, from December through February, was typical as having high antecedent soil moisture conditions, making the occurrence of saturation excess runoff likely.

Hydrograph Separation
To detect the quickflow component from each storm event, the baseflow separation method by Lyne and Hollick was applied, using a low-pass numerical filtering that splits the hydrograph into multiple signals, analogous to the filtering of electronic frequencies [70].This method has been adapted as an automated function in the EcoHydRology R package [71], which was used to derive baseflow and stormflow hydrograph components consistently between sites.Filter parameters were set to 0.925 for time-series having a daily interval [72] and increased to 0.98 for instantaneous time-series [73].A recursive filter was applied forward, backward, and forward again, for a total of three passes.

Selection of Storm Events
Storm event delineation was performed using the Eventseq function from the Hydro-Mad R package [74].This function allows for multiple hyetograph and hydrograph time series to be utilized to set thresholds for the start and termination of rainfall-runoff events.Because of differing watershed characteristics and variable antecedent moisture conditions, the Eventseq function parameters were varied for each event.The event delineation results were verified using hyetograph plots to determine appropriate start and end times.
For determination of events using daily streamflow and precipitation data, the Eventseq function was set as follows: 1.
Starting threshold was set to 7 mm of rainfall depth, 2.
A minimum event gap of 1 day, 3.
A minimum event duration of 2 days, 4.
Ending threshold was varied between 2 mm and 4 mm of runoff depth depending on watershed area and remained below the threshold for at least 1 day.
To delineate events using instantaneous 15 minute streamflow and precipitation data, the Eventseq function was set as follows: 1.
Starting threshold was set to 0.1 mm of rainfall depth, 2.
A minimum gap of 6 hours, 3.
A minimum duration of 6 hours, 4.
The ending threshold was varied using precipitation or runoff depth for each event based on the visual inspection of the recession limb of the hydrograph.

GIS Analysis
The USGS StreamStats web tool [54] was used to investigate the catchment attributes of each candidate watershed to determine the relative similarities between watersheds.Before landscape metrics could be applied to the six burned watersheds, geospatial data were compiled in a GIS to determine the burned area falling within watershed boundaries and individual burn classes were isolated.The preprocessing of geospatial data was conducted in ArcGIS Pro version 2.8.0 (ESRI, 2021) and output for further landscape metric analysis.Watershed boundaries and streamgaging locations were obtained as shapefile and point feature classes from the USGS NHDplus v2.1 database and input into ArcGIS Pro.Then, soil burn severity maps were sourced from BAER Imagery Support and overlain in the GIS to determine the percentage of each burn severity class within the watersheds.These burn severity maps were downloaded in raster format but were converted to polygon feature classes to calculate the percentage coverage of each burn severity class within the watershed boundaries.Burn severity rasters were input into the R environment as simple features using the raster package [75].Then, watershed area shapefiles were read using the sf package, which preserves the geospatial attributes of the features by considering both the vector geometry and spatial projection of each object.Finally, the burn severity rasters were cropped to the respective watershed area, and each burn severity class was isolated.

Statistical Analysis
Descriptive statistics were used to summarize the distribution of runoff coefficients between pre-and post-fire periods for the group of sites.Inferential statistics were used to test for differences between pre-and post-fire periods.All statistical computations and visualizations were created using R version 4.1.2[76].Specifically, box-whisker plots were created to characterize the distribution of runoff coefficients, identify outliers, and compare median values.Due to the small sample size, a non-parametric Mann-Whitney test for paired groups was used to evaluate statistical differences between distributions of runoff coefficients between pre-and post-fire years.Lastly, Spearman's rank correlation was used to investigate the relationship between burn severity derived landscape metrics and change in runoff coefficients.

Event Runoff Ratio (Q r )
Results of the Mann-Whitney U test could not reject the null hypothesis that no difference existed between the median Q r values of the pre-and post-fire groups (Table 4).Although only minor shifts in median Q r were observed at burned sites, some showed slight increases in median Q r (Figure 4).Breitenbush River shows the largest post-fire increase in Q r of 0.09 while also having several outliers corresponding to wet season storm events (0.85 for 12 December 2020 event and 0.89 for 12 January 2021 event).Other sites having increases in Q r include Little North Santiam, with an increase of 0.07, and one outlier of 0.79 occurring during an event on 16 December 2020.The Molalla River saw a 0.04 increase in median Q r and an extreme value of 0.82 occurring during the post-fire period on 12 January 2021.Minimal change in Q r occurred in the North Santiam watershed, but several wet season post-fire events show elevated Q r over pre-fire response.The Clackamas River watershed shows a slight decrease in Q r values from pre-fire median Q r.Thomas Creek and South Santiam control watersheds show increases in Q r during the post-fire period.Both sites show highest Q r values occurring during wet season events, and outliers are less pronounced than in the burned watersheds.Percentage of area burned was significantly correlated to Q r (r s = 0.94, p = 0.005), and percentage of high severity class area was marginally correlated to Q r (r s = 0.77, p = 0.07).Computed landscape metrics did not show a relationship to changes in Q r .

Event Peak Runoff Ratio (Qp)
Like Qr, the null hypothesis of the Mann-Whitney u-test could not be rejected for burned and unburned watersheds, even at a 0.10 significance level (Table 5).The Qp median value was slightly less in the post-fire year than in pre-fire years for control sites

Event Peak Runoff Ratio (Q p )
Like Q r , the null hypothesis of the Mann-Whitney u-test could not be rejected for burned and unburned watersheds, even at a 0.10 significance level (Table 5).The Q p median value was slightly less in the post-fire year than in pre-fire years for control sites (Thomas Creek and South Santiam) and in sites that had the least amount of burned area (North Santiam, North Umpqua, and Clackamas) (Figure 5).The largest change in median Q p values was found in the Breitenbush watershed (43% increase) and North Santiam watershed (33% increase).All other sites showed negative median change in Q p values, with the unburned Thomas Creek showing the largest decrease in median Q p (−55%) and Clackamas River having the largest decrease among burned sites (−60%), percentage area burned and Q p (r s = 0.89, p-value = 0.02), and percentage high burn severity and Q p (r s = 0.89, p-value = 0.02).Due to the short temporal scale of this study, runoff ratios at the event scale were chosen to compare one-year pre-fire and one-year post-fire.It was assumed that change

Event Runoff Ratio (Q r )
Due to the short temporal scale of this study, runoff ratios at the event scale were chosen to compare one-year pre-fire and one-year post-fire.It was assumed that change would be most detectable at the event scale because runoff ratios have been shown to increase most for individual rainfall events over monthly and annual scales [63].The comparison of median Q r and landscape metric values shows that sites with a greater burned area and greater area of high severity patches have a larger magnitude of change in Q r values (Table 2, Figure 4).As expected, the heavily burned Breitenbush and Little North Santiam watersheds show the greatest change in Q r .Both control watersheds have increased median Q r , with a decrease in the number of post-fire events (n = 14) from the pre-fire year (n = 21).Surprisingly, several of the burned watersheds show decreases in Q r values: Clackamas and North Umpqua.These two watersheds hold the largest contributing areas, and it is possible that undetermined precipitation thresholds were not reached to produce enhanced runoff [77,78].In order to determine these thresholds, finerresolution precipitation data are needed to determine the exact duration and intensities of precipitation events.
It is clear that the heavily burned Breitenbush and Little North Santiam watersheds show differences in runoff coefficients that are likely attributed to post-fire effects, but in those watersheds with less fire disturbance, the influence of disturbance is less certain.It is possible that the drought conditions coinciding with water year 2021 may attenuate the hydrologic response to wildfire by reducing the duration, intensity, and number of wet-season precipitation events [79].The influence of climatic effects can be difficult to disentangle from watershed disturbance, and ongoing drought in the PNW might complicate runoff change detection.In addition, drought can decrease soil moisture by increasing evaporative demand, thereby encouraging infiltration and reducing runoff rates [80].

Event Peak Runoff Ratio (Q p )
The event peak runoff ratio, Q p , was included in the study to account for soil moisture conditions leading up to runoff events.In the case of water year 2021, antecedent soil conditions were likely variable due to a reduced number of high-intensity precipitation events, abnormal seasonal variability, and as many regions of the state experienced all-time positive temperature anomalies that increased evaporative demand, effectively lowering soil moisture [81].For comparison, Oregon started the 2020 water year with no counties with drought, but conditions had shifted by the end of September 2020 to 34% of the state experiencing drought.By the end of water year 2021, 96% of Oregon was experiencing some form of drought [82].The severe drought conditions in water year 2021 caused difficulties in selecting a comparable reference water year.This method assumes that climatic variability between pre-and post-fire years can be controlled for by selecting a pre-fire year with similar levels of drought.Because of distinct differences between climatic zones east and west of the Cascade Range, drought indices were examined for the North Santiam watershed, which is centrally located in the study area (Figure 1).Water year 2015 was the most comparable year to 2021 based on the percentage area of the North Santiam watershed in U.S. drought monitoring categories [83].Differences in the percentage of drought levels can be observed between the two years: drought in water year 2015 progressively worsened from abnormally dry winter conditions to moderate and severe drought during the spring and summer months and finally developing into extreme drought near the end of the year.Post-fire water year 2021 began under extreme drought conditions but improved with normal levels of winter precipitation before progressing into extreme drought for much of the summer (Figure 6).drought levels can be observed between the two years: drought in water year 2015 progressively worsened from abnormally dry winter conditions to moderate and severe drought during the spring and summer months and finally developing into extreme drought near the end of the year.Post-fire water year 2021 began under extreme drought conditions but improved with normal levels of winter precipitation before progressing into extreme drought for much of the summer (Figure 6).Drought conditions likely dulled post-fire runoff effects across the study area, although alterations to runoff were still detected in the highly burned Breitenbush and Little North Santiam watersheds.One surprising finding was that the low elevation control site, Thomas Creek, had the greatest decrease in Qp.A possible explanation for this inconsistency was that high valley temperatures might have driven large vapor pressure deficits, thereby increasing transpiration and further lowering soil moisture conditions.Drought conditions likely dulled post-fire runoff effects across the study area, although alterations to runoff were still detected in the highly burned Breitenbush and Little North Santiam watersheds.One surprising finding was that the low elevation control site, Thomas Creek, had the greatest decrease in Q p .A possible explanation for this inconsistency was that high valley temperatures might have driven large vapor pressure deficits, thereby increasing transpiration and further lowering soil moisture conditions.
Although the pre-and post-fire periods do not show different median values, the largest changes in Q p values occurred in the sites with the greatest percentage of burn area and the greatest percentage of high burn area (Breitenbush and Little North Santiam).The median Q p value in the Mollala watershed remained unchanged, despite a widening of the distribution.The most notable difference is the decrease in Q p values for watersheds with less percentage burned including Clackamas, North Umpqua, and North Santiam.Both control watersheds (South Santiam and Thomas Creek) also displayed decreases in median Q p .Generally, these less burned watersheds had tighter distributions except for several large outliers.These results agree with other studies in the WCONUS that found consistent negative streamflow response in post-fire year one of eastern California, western Nevada, and western Oregon watersheds [84].They interpret these results as indicating that burned watersheds in these areas are likely show small to non-existent streamflow responses.Additional findings of regional post-fire research support the lack of runoff response observed in the present study.Runoff response is more likely in Mediterranean and semi-arid watersheds, while high elevation, snow-driven watersheds such as the Cascades and Rocky Mountains are less vulnerable to increases in runoff response [85].

Landscape Metrics
The six landscape metrics selected to represent patch density, elongation, and aggregation (Table 3) did not show correlations with percentage change in either runoff metric.Intuitively, burned area and high severity area were colinear (r s = 0.83, Figure 7).In addition, multicollinearity was present between landscape metrics, suggesting that the selected aggregation metrics were not independent and measured similar patch characteristics.

Limitations
The lack of significant change to post-fire runoff coefficients could be affected by se eral aspects of the study design as well as uncertainties in the data used in the study.Wh water years 2015 and 2021 were closely matched according to their above-normal temp ature and precipitation totals, the timing and type of precipitation were different.In 20 the state of Oregon experienced record low snowpack despite precipitation totals stayi close to normal [89].The shift from winter snow accumulation to rain, as observed in 20 is considered the result of climate warming trends, as opposed to climate variability [9 In contrast, water year 2021 had normal snowpack, but exceptionally dry spring con tions led to extremely low soil moisture conditions across much of Oregon, which p gressed into the warmest summer on record [81].
Precipitation data used in calculating event-based runoff coefficients contain unc tainty due to the lack of precipitation observations across the study region.The 4 km r olution gridded PRISM daily time series do not include a measure of uncertainty, but i assumed that increased error of modeled values would be associated with greater distan from locations of observed precipitation.Another source of precipitation uncertainty the relatively coarse temporal resolution of the time-series, which limit precise delineati of storm events.
While differences between watershed characteristics were observed (Table 1), su differences were not directly accounted for between watersheds.Specifically, the wat sheds have widely ranging permeability values (Table 1), dependent on unique combin These findings are consistent with several post-fire studies that find percentage area burned to be significant in describing changes to hydrologic response [84,86].Most applicable to this study, an investigation into the influence of percentage area burned, vegetation recovery, and soil conditions in southern and central California shows that regional streamflow increases among watersheds with percentages burned above 20% and are most pronounced in wet years [87].
Several explanations exist for the lack of correlation between landscape metrics and runoff ratios.Although landscape metrics have been applied in watershed hydrology before, they were first developed to aid ecological studies of structural connectivity.They might fail to represent basic mechanics of the functional connectivity that govern runoff pathways.Some have attempted to quantify patch-based processes by integrating ecological and process-based connectivity metrics with node-based theory [88], which bypasses numeric descriptions of aggregation and could more realistically represent the directional connectivity present in a burned watershed.Second, processes that influence runoff generation may be present at different scales than what can be described with burned patch landscape metrics.Studies into the patchiness of hydrophobic soils have found evidence for small-scale hydrologic connectivity in shrub-dominated burned areas [29].The spatial variability of burned forest soils may display discontinuous areas of hydrophobicity, ash, and macro pores, likely complicating the quantification of hydrologic connectivity across broad catchments.The current study sought to increase the spatial scale of investigation to the catchment scale and assumes that high severity burn classes can accurately represent conditions that allow for increased runoff magnitude.

Limitations
The lack of significant change to post-fire runoff coefficients could be affected by several aspects of the study design as well as uncertainties in the data used in the study.While water years 2015 and 2021 were closely matched according to their above-normal temperature and precipitation totals, the timing and type of precipitation were different.In 2015, the state of Oregon experienced record low snowpack despite precipitation totals staying close to normal [89].The shift from winter snow accumulation to rain, as observed in 2015, is considered the result of climate warming trends, as opposed to climate variability [90].In contrast, water year 2021 had normal snowpack, but exceptionally dry spring conditions led to extremely low soil moisture conditions across much of Oregon, which progressed into the warmest summer on record [81].
Precipitation data used in calculating event-based runoff coefficients contain uncertainty due to the lack of precipitation observations across the study region.The 4 km resolution gridded PRISM daily time series do not include a measure of uncertainty, but it is assumed that increased error of modeled values would be associated with greater distance from locations of observed precipitation.Another source of precipitation uncertainty is the relatively coarse temporal resolution of the time-series, which limit precise delineation of storm events.
While differences between watershed characteristics were observed (Table 1), such differences were not directly accounted for between watersheds.Specifically, the watersheds have widely ranging permeability values (Table 1), dependent on unique combinations of geology and soils.Given that watersheds in the High Cascades (i.e., Clackamas and North Santiam) have buffering capacity due to groundwater input, post-fire stream response might have been delayed or have not been detected yet.Additionally, the size and location of burned areas, slope, and aspect are not consistent between watersheds, potentially challenging the comparability of post-fire runoff values.
Finally, this study assumes that all watersheds have similar pre-fire land management histories and that the post-fire response was consistent.Depending on land ownership, postfire salvage logging likely influenced these watersheds as private, federal, and state land managers had operations began immediately after the 2020 fire season due to categorical exclusions to permitting processes [91].Additional post-fire investigations should consider varying land management practices across these watersheds.

Conclusions
Event runoff ratios were computed for Oregon watersheds burned by the Labor Day Fires of 2020 to determine if post-fire runoff response was altered, whether the spatial arrangement of burn patterns was related to alterations, and at what temporal resolution the disturbance was most evident.Linking changes to burn severity patterns was difficult due to the non-significantly different runoff coefficient distributions, but climatic trends could play a role in concealing post-fire runoff characteristics.The role of drought on postfire runoff dynamics was evident in the two control watersheds that both had declining Q p .This finding suggests that the magnitude of runoff responses could have been larger in burned watersheds had the runoff enabling effect of reduced canopy interception on throughfall been dominant over evapotranspiration driven by high summer temperatures.Breitenbush and Little North Santiam watersheds, two of the most burned watersheds, had the most positive change in runoff coefficients.It is assumed that fire-inflicted disturbances were widespread enough to overtake the influence of increasing temperature trends.
Plot and small catchment scale studies [24,26,29,37] have indicated that the burned patch arrangement is important for determining the degree to which rainfall events are intercepted, allowed to infiltrate, or find efficient surface pathways.This study aimed to build upon these findings by increasing the spatial scale at which burned patchiness might influence watershed hydrology.However, the aforementioned issues of drought-induced soil conditions limit such conclusions.Another limitation is that the selected landscape metrics used in this study were not developed within a hydrologic framework and may

Hydrology 2022, 9 , 5 Figure 1 .
Figure 1.Surface geology of the study area displaying rocks of the Western Cascade Province the High Cascade Province.

Figure 1 .
Figure 1.Surface geology of the study area displaying rocks of the Western Cascade Province and the High Cascade Province.

Figure 2 .
Figure 2. Generalized geological cross section of the study area where (A) represents the High Cascade Province, which is more permeable, allowing for rapid groundwater exchange, while the Western Cascade Province (B) encourages surface runoff due to less permeable characteristics.

Figure 2 .
Figure 2. Generalized geological cross section of the study area where (A) represents the High Cascade Province, which is more permeable, allowing for rapid groundwater exchange, while the Western Cascade Province (B) encourages surface runoff due to less permeable characteristics.

Figure 3 .
Figure 3.Oregon study area including the location USGS streamgaging stations, watershed bo aries, and burn severity classes.Low burn severity includes areas where surface organics have been completely combusted, roots are intact, and some patches of grey ash are present.Mode burn severity includes areas where some surface organics and roots have been consumed, s needles and leaves may remain, and patches of grey/black ash are present.High burn severit cludes areas where almost all the organic groundcover has been destroyed, many surface root consumed, bare soil is exposed and soil structure has been disturbed, stumps and trees are c pletely burned out in places, and thick deposits of white and grey ash are present[54].

Figure 4 .
Figure 4. Box whisker plots of Qr values and for pre-fire and post-fire periods for eight Oregon Cascade watersheds.

Figure 4 .
Figure 4. Box whisker plots of Q r values and for pre-fire and post-fire periods for eight Oregon Cascade watersheds.

Figure 5 .
Figure 5. Box whisker plots of Qp values and for pre-fire and post-fire periods for eight Oregon Cascade watersheds.4. Discussion 4.1.Event Runoff Ratio (Qr)

Figure 5 .
Figure 5. Box whisker plots of Q p values and for pre-fire and post-fire periods for eight Oregon Cascade watersheds.

Figure 6 .
Figure 6.Comparison of pre-and post-fire year percentage area of North Santiam watershed (HUC 17090005) in U.S. drought monitoring categories [83].

Figure 6 .
Figure 6.Comparison of pre-and post-fire year percentage area of North Santiam watershed (HUC 17090005) in U.S. drought monitoring categories.The figure was created using the web-based tool available at U.S. Drought Monitor [83].

Figure 7 .
Figure 7. Spearman's rank correlation (rs) matrix for change in post-fire runoff coefficients to burned area landscape metrics.Only significant (p-value < 0.05) correlation values are shown.

Figure 7 .
Figure 7. Spearman's rank correlation (r s ) matrix for change in post-fire runoff coefficients to burned area landscape metrics.Only significant (p-value < 0.05) correlation values are shown.

Table 2 .
Characteristics of six burned watersheds in Oregon with the total percentage burned and percentage of each burn severity class (low, moderate, and high burned areas).

Table 3 .
Landscape metrics calculated for the high severity patches in each burned watershed.

Table 4 .
Q r results of Mann-Whitney U test and change in Q r values for six burned watersheds (B) and two controlled watersheds (C).The significance level was considered p-value < 0.10.

Table 5 .
Q p results of Mann-Whitney U test and change in median Q p values and interquartile range (IQR) of pre-and post-fire Q p values for six burned watersheds (B) and two controlled watersheds (C).The significance level was considered at a p value < 0.10.