Applying High-Resolution Imagery to Evaluate Restoration-Induced Changes in Stream Condition , Missouri River Headwaters Basin , Montana

Degradation of streams and associated riparian habitat across the Missouri River Headwaters Basin has motivated several stream restoration projects across the watershed. Many of these projects install a series of beaver dam analogues (BDAs) to aggrade incised streams, elevate local water tables, and create natural surface water storage by reconnecting streams with their floodplains. Satellite imagery can provide a spatially continuous mechanism to monitor the effects of these in-stream structures on stream surface area. However, remote sensing-based approaches to map narrow (e.g., <5 m wide) linear features such as streams have been under-developed relative to efforts to map other types of aquatic systems, such as wetlands or lakes. We mapped preand post-restoration (one to three years post-restoration) stream surface area and riparian greenness at four stream restoration sites using Worldview-2 and 3 images as well as a QuickBird-2 image. We found that panchromatic brightness and eCognition-based outputs (0.5 m resolution) provided high-accuracy maps of stream surface area (overall accuracy ranged from 91% to 99%) for streams as narrow as 1.5 m wide. Using image pairs, we were able to document increases in stream surface area immediately upstream of BDAs as well as increases in stream surface area along the restoration reach at Robb Creek, Alkali Creek and Long Creek (South). Although Long Creek (North) did not show a net increase in stream surface area along the restoration reach, we did observe an increase in riparian greenness, suggesting increased water retention adjacent to the stream. As high-resolution imagery becomes more widely collected and available, improvements in our ability to provide spatially continuous monitoring of stream systems can effectively complement more traditional field-based and gage-based datasets to inform watershed management.


Introduction
Remotely sensed imagery has been widely applied to characterize variability in surface-water extent across space and time [1,2].The spatial resolution (≥30 m) of commonly used sources of imagery (e.g., Landsat, MODIS, AVHRR), however, has limited our ability to remotely monitor river systems, except for large rivers (e.g., >40 m wide) [3] or rivers under flood conditions [4][5][6][7].Yet remote monitoring of the spatial distribution of river stage and condition has several applications including enhancing our ability to predict and monitor flood events, informing the source and distribution of flow to downstream gaged points, helping monitor ungaged watersheds, predicting carbon dioxide emissions, and informing river management [8][9][10][11].The rapidly increasing availability of multispectral, high-resolution imagery (≤5 m resolution, Dove, RapidEye (Planet, San Francisco, CA, USA), Worldview-2, 3 (DigitalGlobe, Westminster, CO, USA)) provides increased opportunity to potentially monitor river systems across diverse watershed sizes and flow conditions.
Multiple sources of fine resolution imagery have been applied to aquatic systems.LiDAR [12,13] and synthetic aperture radar (SAR) imagery have been successfully used to map surface water and can be preferable in forested environments or during storm events under cloud cover [14][15][16].Multispectral, high-resolution imagery has also been used effectively to map surface water [17,18].Riverscape units including the active channel have primarily been mapped by digitizing very high-resolution multispectral satellite imagery or aerial imagery [19][20][21] or by applying geographical object-based image analysis (GEOBIA) methods [11,22,23].A GEOBIA approach segments an image into homogenous objects prior to object classification.Such an approach can help account for the greater within-class spectral variability that can occur with high-resolution imagery, relative to moderate-resolution imagery [18].In general, however, efforts to remotely monitor narrow, linear water features, such as rivers and streams, have lagged behind efforts to remotely monitor lakes and wetlands [24,25].
The Upper Missouri River Headwaters Basin in southwestern Montana faces increasingly uncertain water supplies attributable to high water demand for agricultural irrigation [26,27] and public water supply [28,29].In addition, shifts in the timing of runoff and peak streamflow are predicted with increasing amounts of winter precipitation and a declining snowpack related to climate change [30][31][32].Societal water demands as well as climate-induced shifts in streamflow can threaten habitat critical for fish and aquatic species [33,34].These risks have raised interest in increasing the capacity of streams to respond to extreme events [35][36][37].One approach that is growing in popularity is to slow runoff, absorb excess floodwater, and encourage groundwater recharge by enhancing natural water storage in stream channels, riparian areas, and floodplains [38,39].One way to create natural water storage is using in-stream, channel spanning structures called beaver dam analogues (BDAs) [40][41][42].Over time BDAs have been shown to slow water flow, encourage channel stability and riparian vegetation, activate side channels, and improve water quality and fish habitat [40,41,43].Installing BDAs along a reach of stream can potentially increase spring overbank flow and elevate riparian water tables [43].If water is a limiting factor in the riverscape, elevating near-surface soil moisture along stream channels can encourage riparian vegetation [44,45].Depending on channel shape and river stage, slowing the flow of water can increase the stream surface area through ponding and temporary flooding [38,45].Downstream from beaver dams, channels are more likely to be stable, with lower sediment loads and a lower range of stream discharge [45].Although the installation of BDAs is becoming relatively common, analyses that evaluate their impact are limited.Most studies to date have focused on the impact of actual beaver dams that differ from BDAs in the magnitude of change and maintenance activity.In this study, we tested how multispectral high-resolution imagery can be used to monitor stream condition along four stream reaches in the Upper Missouri Headwaters Basin, and how pairs of images can potentially be used to monitor BDA stream restoration projects.Our research questions included: 1.
What methodological approaches are most effective to map stream surface area using multispectral high-resolution imagery?And, 2.
How can image pairs (e.g., pre-and post-restoration) be used to monitor changes in stream surface area and riparian greenness?

Study Area and Restoration Activities
The four stream sites examined in this study occur within the Upper Missouri Headwaters Basin in southwest Montana (Figure 1).Annual precipitation across the Headwaters Basin averages 565 mm yr −1 , while the annual temperature maximum and minimum average 10 • C and −3 • C, respectively   [46].Across the basin herbaceous vegetation (35%) and shrub/scrub (20%) dominate the large river valleys while evergreen forest dominates the higher elevations (35%) [47].Two restoration sites occurred along reaches separated by 4.8 stream km in Long Creek (Figure 1), which flows south into the Red Rock River in the Red Rock River Hydrological Unit (Red Rock HUC8).Land cover adjacent to Long Creek is dominated by herbaceous vegetation, shrub/scrub, and emergent herbaceous wetlands [47].The third restoration site occurred along Alkali Creek, which flows northwest into Blacktail Deer Creek in the Beaverhead HUC8.This site showed evidence of beaver activity just upstream from the restoration site.The fourth restoration site occurred along Robb Creek that flows north into the Ruby River (Ruby River HUC8).Land cover adjacent to both Alkali Creek and Robb Creek is dominated by herbaceous vegetation and shrub/scrub habitat [47].
Remote Sens. 2018, 10, x FOR PEER REVIEW 3 of 29 while the annual temperature maximum and minimum average 10 °C and −3 °C, respectively (1981-2010) [46].Across the basin herbaceous vegetation (35%) and shrub/scrub (20%) dominate the large river valleys while evergreen forest dominates the higher elevations (35%) [47].Two restoration sites occurred along reaches separated by 4.8 stream km in Long Creek (Figure 1), which flows south into the Red Rock River in the Red Rock River Hydrological Unit (Red Rock HUC8).Land cover adjacent to Long Creek is dominated by herbaceous vegetation, shrub/scrub, and emergent herbaceous wetlands [47].The third restoration site occurred along Alkali Creek, which flows northwest into Blacktail Deer Creek in the Beaverhead HUC8.This site showed evidence of beaver activity just upstream from the restoration site.
The fourth restoration site occurred along Robb Creek that flows north into the Ruby River (Ruby River HUC8).Land cover adjacent to both Alkali Creek and Robb Creek is dominated by herbaceous vegetation and shrub/scrub habitat [47].All restoration activities were developed and completed by the Nature Conservancy.A series of BDAs were installed in stream reaches at each restoration site.The structures were created from wooden posts installed vertically into the streambed across the channel with willow branches woven between posts.The structures collect organic material and sediment behind them, building up the stream bed height, ponding water upstream from the structures, stabilizing the channel and increasing connectivity with its floodplain [48].The BDAs were accompanied by willow plantings along the stream to stabilize banks and cattle exclusions at most of the sites [48].The design is cost effective as no heavy equipment is used and the in-stream structures are designed to be temporary [49].BDAs were installed in two reaches of Long Creek (9 BDAs on the north reach and 7 BDAs on the south reach), a reach of Alkali Creek (6 BDAs) and a reach of Robb Creek (12 BDAs) (Table 1).The goal of the restorations at the Long Creek and Alkali Creek sites was to aggrade the streambed, improving hydrologic connectivity between the stream channel and associated floodplains.Along Robb Creek, the BDAs were designed to encourage reactivation of abandoned side channels.The time since restoration ranged from one to three years across the sites (     1. Characteristics of the restoration sites and a summary of the restoration activities performed at each site.A list of the high-resolution images representing preand post-restoration conditions across the four sites is also shown.Length refers to the stream distance from the most upstream to most downstream beaver dam analog (BDA).

Image Acquisition and Preprocessing
A total of six high-resolution images (2 m resolution) were acquired from DigitalGlobe (Westminster, CO, USA) via the NextView license for this analysis.These images included one QuickBird-2 image, three Worldview-2 images, and two Worldview-3 images (Table 1)."Pre-restoration" conditions were represented by images acquired during summer 2014, while "post-restoration" conditions were represented by images acquired during summer 2017.Using historical (1895-2017) Palmer Hydrological Drought Index (PHDI) values, we found that the preand post-restoration image dates represented similar historical wetness conditions (37.6% relative to 41.8% PHDI); however, 2013, the year prior to the pre-restoration images experienced a drought which may have influenced stream conditions in 2014 (Figure 2).We also compared the stream discharge values on the date the image was collected using a USGS stream gage downstream from the four restoration sites (Jefferson River, USGS Gage #06026500) (Figure 1).Discharge was reasonably similar between the image dates (10% higher post-restoration for the Long Creek image dates and 8% lower post-restoration for the Robb Creek image dates), for all sites except Alkali Stream.The post-restoration image at this site was collected in August when discharge was much lower relative to the early summer period (Figure 2).We converted the image (processing Level 1) pixel values from Digital Numbers to top-of-atmosphere reflectance in PCI Geomatica.For each image the panchromatic and multispectral bands were orthorectified together using PCI Geomatica's 2014 OrthoEngine.National Agricultural Imagery Program (NAIP) images (1 m resolution) were used as reference images (Long Creek and Robb Creek-22 October 2015, Alkali Creek-3 August 2013) together with the U.S. Geological Survey's 10 m National Elevation Dataset (NED) [50].Images were pan-sharpened to 50 cm resolution using PCI Geomatica's PANSHARP2 tool [51].The panchromatic band width for QuickBird-2 extends across all four of the multispectral bands (blue, green, red, NIR); however, the panchromatic band for Worldview-2 and Worldview-3 extends across only five of the eight spectral bands (excludes the coastal band, NIR1 and NIR2 bands).Prior work has shown that applying pan-sharpening methods to bands outside of the panchromatic range can distort the values in these bands [52].To account for this, we pan-sharpened only the bands overlapping the panchromatic band using the PANSHARP2 tool.Bands outside of the panchromatic band range were resampled to 50 cm resolution using cubic convolution.

Image Acquisition and Preprocessing
A total of six high-resolution images (2 m resolution) were acquired from DigitalGlobe (Westminster, CO, USA) via the NextView license for this analysis.These images included one QuickBird-2 image, three Worldview-2 images, and two Worldview-3 images (Table 1)."Prerestoration" conditions were represented by images acquired during summer 2014, while "postrestoration" conditions were represented by images acquired during summer 2017.Using historical (1895-2017) Palmer Hydrological Drought Index (PHDI) values, we found that the pre-and postrestoration image dates represented similar historical wetness conditions (37.6% relative to 41.8% PHDI); however, 2013, the year prior to the pre-restoration images experienced a drought which may have influenced stream conditions in 2014 (Figure 2).We also compared the stream discharge values on the date the image was collected using a USGS stream gage downstream from the four restoration sites (Jefferson River, USGS Gage #06026500) (Figure 1).Discharge was reasonably similar between the image dates (10% higher post-restoration for the Long Creek image dates and 8% lower postrestoration for the Robb Creek image dates), for all sites except Alkali Stream.The post-restoration image at this site was collected in August when discharge was much lower relative to the early summer period (Figure 2).We converted the image (processing Level 1) pixel values from Digital Numbers to top-of-atmosphere reflectance in PCI Geomatica.For each image the panchromatic and multispectral bands were orthorectified together using PCI Geomatica's 2014 OrthoEngine.National Agricultural Imagery Program (NAIP) images (1 m resolution) were used as reference images (Long Creek and Robb Creek-October 22, 2015, Alkali Creek-August 3, 2013) together with the U.S. Geological Survey's 10 m National Elevation Dataset (NED) [50].Images were pan-sharpened to 50 cm resolution using PCI Geomatica's PANSHARP2 tool [51].The panchromatic band width for QuickBird-2 extends across all four of the multispectral bands (blue, green, red, NIR); however, the panchromatic band for Worldview-2 and Worldview-3 extends across only five of the eight spectral bands (excludes the coastal band, NIR1 and NIR2 bands).Prior work has shown that applying pansharpening methods to bands outside of the panchromatic range can distort the values in these bands [52].To account for this, we pan-sharpened only the bands overlapping the panchromatic band using the PANSHARP2 tool.Bands outside of the panchromatic band range were resampled to 50 cm resolution using cubic convolution.

Object-Based Water Classification
We used the software eCognition (version 9.2.1, Trimble, Westminster, CO, USA) to process the high-resolution images into maps of surface-water extent.This software uses an object-oriented approach where an image is first segmented into objects representing meaningful features of the physical landscape, and the objects are then classified using user-defined rules and algorithms.Rules can be set in a hierarchical order as child rules under a parent process so that the ruleset can be automatically run in sequence.In this case, our objective was to segment, then classify objects into water and non-water where the objects of interest were the stream channels.Each pan-sharpened image was first clipped to the spatial area of interest, which included a minimum of 300 m stream length upstream from the restoration reach, the restoration reach and approximately 1 stream km downstream of the restored site (Figure 3).The stream length of the restoration reach (from the upstream to the downstream BDA) ranged from 830 m at Alkali Creek to 3.8 km at Long Creek (North) (Table 1).
To segment each image into objects we focused the segmentation along edges or sharp contrasts in the image.To do this we first modified the panchromatic band.Within eCognition, the Edge Extraction Lee Sigma filter was applied to the panchromatic band to create a (1) bright edge layer and (2) dark edge layer from the original image.The dark edge layer was then added to and the bright edge layer subtracted from an inverted version of the panchromatic band to enhance the edge contrasts along streams.The edge-enhanced panchromatic band was then used with the pansharpened blue, green, red, and near-infrared bands to guide the initial image segmentation (scale = 50-100).This approach focused segmentation along stream boundaries while segmenting the image at larger scales, where scale refers to the maximum allowable heterogeneity within an object.The larger scale was desired so that individual trees and tree shadows, which are spectrally similar to water, were segmented into larger objects that contained multiple trees, tree shadows, and intervening vegetation (Figure 4A).After the initial image segmentation, the objects were classified using spectral indices applied in a hierarchical rule-based approach.Our goal was to classify all objects as (1) water; (2) vegetation; (3) soil; or (4) shadow.Objects that were not spectrally similar enough to fit in any of these categories were classified temporarily as (1) water candidates (i.e., potentially water) or (2) unclassified.Objects classified as water candidates or unclassified using the rule-based approach were then re-segmented at a finer scale (scale = 15) to create smaller and more spectrally pure objects (Figure 4B).

Object-Based Water Classification
We used the software eCognition (version 9.2.1, Trimble, Westminster, CO, USA) to process the high-resolution images into maps of surface-water extent.This software uses an object-oriented approach where an image is first segmented into objects representing meaningful features of the physical landscape, and the objects are then classified using user-defined rules and algorithms.Rules can be set in a hierarchical order as child rules under a parent process so that the ruleset can be automatically run in sequence.In this case, our objective was to segment, then classify objects into water and non-water where the objects of interest were the stream channels.Each pan-sharpened image was first clipped to the spatial area of interest, which included a minimum of 300 m stream length upstream from the restoration reach, the restoration reach and approximately 1 stream km downstream of the restored site (Figure 3).The stream length of the restoration reach (from the upstream to the downstream BDA) ranged from 830 m at Alkali Creek to 3.8 km at Long Creek (North) (Table 1).
To segment each image into objects we focused the segmentation along edges or sharp contrasts in the image.To do this we first modified the panchromatic band.Within eCognition, the Edge Extraction Lee Sigma filter was applied to the panchromatic band to create a (1) bright edge layer and (2) dark edge layer from the original image.The dark edge layer was then added to and the bright edge layer subtracted from an inverted version of the panchromatic band to enhance the edge contrasts along streams.The edge-enhanced panchromatic band was then used with the pan-sharpened blue, green, red, and near-infrared bands to guide the initial image segmentation (scale = 50-100).This approach focused segmentation along stream boundaries while segmenting the image at larger scales, where scale refers to the maximum allowable heterogeneity within an object.The larger scale was desired so that individual trees and tree shadows, which are spectrally similar to water, were segmented into larger objects that contained multiple trees, tree shadows, and intervening vegetation (Figure 4A).After the initial image segmentation, the objects were classified using spectral indices applied in a hierarchical rule-based approach.Our goal was to classify all objects as (1) water; (2) vegetation; (3) soil; or (4) shadow.Objects that were not spectrally similar enough to fit in any of these categories were classified temporarily as (1) water candidates (i.e., potentially water) or (2) unclassified.Objects classified as water candidates or unclassified using the rule-based approach were then re-segmented at a finer scale (scale = 15) to create smaller and more spectrally pure objects (Figure 4B).Several spectral indices were included in our hierarchical rule-based approach (Table 2).Worldview-2 and 3 images (8 spectral bands) provide data from several bands not available in Landsat TM, ETM+ or QuickBird-2 including a coastal band, red edge band and two separate near infrared bands (NIR1 and NIR2), offering opportunities for unique band combinations [53].The Normalized Difference Water Index (NDWI) [54] and the Worldview Water Index (WWI) [55] (Table 2) were used as the primary means to identify water objects.Objects were classified as water when either the NDWI or WWI object values were greater than zero.Objects were identified as vegetation when the objects had high Enhanced Vegetation Index (EVI) [56] values, using the coefficients generally adopted [57], or when they showed both a high Normalized Difference Vegetation Index (NDVI) [58] as well as a minimal difference between NDVI and EVI.Objects were classified as soil when the Green-Red Vegetation Index (GRVI) [59] value of objects was less than zero or where brightness values, derived from the panchromatic band, were high (Table 2).Shadow objects were identified using the normalized difference between the coastal band and blue band.For QuickBird-2 Several spectral indices were included in our hierarchical rule-based approach (Table 2).Worldview-2 and 3 images (8 spectral bands) provide data from several bands not available in Landsat TM, ETM+ or QuickBird-2 including a coastal band, red edge band and two separate near infrared bands (NIR1 and NIR2), offering opportunities for unique band combinations [53].The Normalized Difference Water Index (NDWI) [54] and the Worldview Water Index (WWI) [55] (Table 2) were used as the primary means to identify water objects.Objects were classified as water when either the NDWI or WWI object values were greater than zero.Objects were identified as vegetation when the objects had high Enhanced Vegetation Index (EVI) [56] values, using the coefficients generally adopted [57], or when they showed both a high Normalized Difference Vegetation Index (NDVI) [58] as well as a minimal difference between NDVI and EVI.Objects were classified as soil when the Green-Red Vegetation Index (GRVI) [59] value of objects was less than zero or where brightness values, derived from the panchromatic band, were high (Table 2).Shadow objects were identified using the normalized difference between the coastal band and blue band.For QuickBird-2 images in which the coastal band was not available this index was adapted to the normalized difference between the blue band and green band.A low brightness threshold, derived from the panchromatic band, was also used to identify shadowed areas.Finally, water candidates (i.e., potentially water) were identified using a series of indices.The rules to identify objects as water candidates were applied as child rules below the rules identifying objects as soil or vegetation.Indices used to identify water candidates included the NDWI using the coastal band instead of the green band [60] and panchromatic brightness, both of which were effective at identifying deeper water, as well as several novel band combinations including the Red Edge NDWI and Red Edge WWI both of which were helpful in mapping shallow water, particularly where sandy soils were visible below the water, and the Normalized Difference Coastal Red Edge Index (NDCREI) which was helpful in identifying turbid water.A list of the indices used, and the band combinations are shown in Table 2.
images in which the coastal band was not available this index was adapted to the normalized difference between the blue band and green band.A low brightness threshold, derived from the panchromatic band, was also used to identify shadowed areas.Finally, water candidates (i.e., potentially water) were identified using a series of indices.The rules to identify objects as water candidates were applied as child rules below the rules identifying objects as soil or vegetation.Indices used to identify water candidates included the NDWI using the coastal band instead of the green band [60] and panchromatic brightness, both of which were effective at identifying deeper water, as well as several novel band combinations including the Red Edge NDWI and Red Edge WWI both of which were helpful in mapping shallow water, particularly where sandy soils were visible below the water, and the Normalized Difference Coastal Red Edge Index (NDCREI) which was helpful in identifying turbid water.A list of the indices used, and the band combinations are shown in Table 2.The hierarchical rule-based approach of classification was initially applied to all objects (scale = 50-100).For objects that were classified as water candidates or unclassified by the initial rule-based approach, these objects were re-segmented to create smaller, more homogenous objects (scale = 15) and the rule-based approach was re-applied to this subset of smaller objects.For objects classified as water candidates or unclassified after both rounds of rule-based classifications, we applied a Random Forest classifier to determine if these remaining objects were water, vegetation, soil, or shadow.The Random Forest classifier was trained using the objects already classified (either using scale = 50-100 or scale = 15) by the rule-based approach.Bootstrap iterations (n = 500) were run using all indices shown in Table 2, the individual band values, and the band standard deviations of the objects as independent variables.
Following the classification of all objects, in images that contained dense riparian vegetation with shadows, we applied the "grow region" algorithm to the shadow class in eCognition.This step reclassified the objects neighboring shadow objects as shadow candidates.We repeated the process as needed.More inclusive shadow thresholds were then applied to the shadow candidates and the objects were converted to shadow if they were within the thresholds.The object-based elliptic fit shape attribute was also used to classify individual tree shadows for trees that occurred near the stream.Although a similar segmentation and image classification approach was applied to all high-resolution images, as is common in the eCognition environment, a trial-and-error approach to segmentation and image classification was used [61] so that segmentation scale, index thresholds and rulesets were not identical across images.A flowchart showing our eCognition methods from image segmentation to an output of water and non-water is shown in Figure 5.

Pixel-Based Water Classification
Because a GEOBIA approach can be time-intensive and site specific, we were interested in comparing the performance of eCognition outputs with outputs producing using simplistic, single, spectral index thresholds.Using the pan-sharpened TOA reflectance values, we calculated the (1) NDWI [54]; (2) WWI [55]; and (3) the panchromatic brightness value.Brightness was calculated as the pixel value of the panchromatic band, a grayscale image of portions of the electromagnetic spectrum (Worldview-2, 3 (450-800 nm) and QuickBird-2 (450-900 nm)) (e.g., [62]).These three indices were selected as they were most prominently used to identify water objects in the eCognition image processing approach.The spectral index values of the validation points, described in Section 2.5, were used to guide the threshold selection.The rasters were thresholded to water and non-water using the maximum Youden's index, which maximizes the difference between the true positive rate and the false positive rate from the ROC curve and provides an optimal threshold independent from class prevalence [63,64].The Youden's Index optimal threshold and corresponding AUC were calculated for each of the pixel-based outputs.

Stream Surface Area Validation
At each of the four sites, a stream line was manually delineated along each stream reach.The stream line was buffered (above and below the streamline, 1 m, 1.5 m, 1.8 m, 0.75 m for Alkali Creek,

Pixel-Based Water Classification
Because a GEOBIA approach can be time-intensive and site specific, we were interested in comparing the performance of eCognition outputs with outputs producing using simplistic, single, spectral index thresholds.Using the pan-sharpened TOA reflectance values, we calculated the (1) NDWI [54]; (2) WWI [55]; and (3) the panchromatic brightness value.Brightness was calculated as the pixel value of the panchromatic band, a grayscale image of portions of the electromagnetic spectrum (Worldview-2, 3 (450-800 nm) and QuickBird-2 (450-900 nm)) (e.g., [62]).These three indices were selected as they were most prominently used to identify water objects in the eCognition image processing approach.The spectral index values of the validation points, described in Section 2.5, were used to guide the threshold selection.The rasters were thresholded to water and non-water using the maximum Youden's index, which maximizes the difference between the true positive rate and the false positive rate from the ROC curve and provides an optimal threshold independent from class prevalence [63,64].The Youden's Index optimal threshold and corresponding AUC were calculated for each of the pixel-based outputs.

Stream Surface Area Validation
At each of the four sites, a stream line was manually delineated along each stream reach.The stream line was buffered (above and below the streamline, 1 m, 1.5 m, 1.8 m, 0.75 m for Alkali Creek, Long Creek (North), Long Creek (South), and Robb Creek, respectively) so that the total buffered area represented the average stream width.Points were randomly selected within the buffered stream area to represent water points (n = 200).The stream line was then buffered by 50 m and 200 points were randomly selected within the buffered area to represent non-water points.The 400 validation points per image were visually inspected using the raw pan-sharpened image to confirm status (water or non-water).Accuracy metrics calculated included overall accuracy, omission error, commission error, Dice coefficient, and relative bias.Omission and commission errors were calculated for the category of water.The Dice coefficient is the conditional probability that if one classifier (product or reference data) identifies a pixel as water, the other one will as well, integrating errors of omission and commission [65,66].The relative bias provides the proportion that water is underestimated (negative bias) or overestimated (positive bias).Accuracy metrics were calculated for each of the pixel-based and object-based stream surface area outputs and presented by site, year, and methodology.

Changes in Stream Surface Area
To evaluate changes in stream surface area, our goal was to select the most accurate pair of stream surface area maps per site.This was determined using both the accuracy statistics as well as a visual assessment of quality.For Long Creek (North) and Long Creek (South), we used the panchromatic brightness output (97.4 and 98.4% overall accuracy, respectively when averaged across the two years, Table 3).For Robb Creek, we used the eCognition outputs (96.6% overall accuracy, averaged across the two years, Table 3), and for Alkali Creek, we merged the eCognition outputs with the panchromatic brightness outputs so that if water was identified by either output it was included as water.This was necessary only for Alkali Creek as the stream segment immediately downstream from the restoration was narrow (~1 to 1.5 m wide) and not adequately mapped using the eCognition output alone.Each output was edited manually to remove errors of commission.All outputs were converted to polygons and projected to WGS 1984 UTM zone 12 N prior to calculating area (m 2 ) in ArcGIS 10.3 (ESRI, Redlands, CA, USA).
All surface water continuous with the stream centerline was included in stream surface area calculations, while waterbodies that were disconnected from the stream centerline were excluded from the stream surface area calculations.The total stream surface area was calculated for (1) the restoration reach, which extended from the upstream BDA structure to the most downstream BDA structure; (2) a reach extending upstream from the restoration reach (stream length of 300 m), and three reaches extending downstream from the restoration reach including; (3) a stream length of 0 m to 250 m downstream; (4) a stream length of 250 m to 500 m downstream; and (5) a stream length of 500 m to 1 km downstream from the restored reach.A schematic showing the distribution of these reaches at each site is shown in Figure 3.The stream length evaluated upstream and downstream from the restoration reach was limited by the extent of the pre-and post-restoration images.The same five reaches were used in analysis of the 2014 and 2017 images.At Robb Creek, the BDAs were aimed at reactivating side channels so each of the side channels were considered separately from the main channel (Figure 3).At Alkali Creek, beaver activity upstream from the restoration site influenced conditions in the upstream reach, so that the reach containing extensive beaver activity was analyzed separately (Figure 3).In addition, dense riparian vegetation approximately 300 m downstream from the restoration site in Alkali Creek limited our analysis of downstream area, so that results were presented only for a single downstream reach (0 m to 300 m) (Figure 3).Table 3.A comparison of the accuracy of methods to map stream extent across sites and image dates.Errors are presented for the accuracy of mapping stream surface area extent.All methods listed are pixel-based except for eCognition which is object-based.WWI: Worldview Water Index; NDWI: Normalized Difference Water Index; WV2: Worldview-2; WV3: Worldview-2; QB2: QuickBird-2; OE: omission error; CE: commission error; OA: overall accuracy; DC: Dice coefficient; RB: relative bias; AUC: area under curve; NIR: near-infrared.The challenge in evaluating the effects of restoration is to separate change attributable to natural variation (interannual, seasonal, event) from change attributable to the restoration activities.

Site and
To accomplish this, we assumed that variability in hydro-climatic conditions would propagate similarly at the scale of adjacent or nearby stream reaches.Therefore, the ratio between the condition (e.g., amount of water, greenness) of an upstream reach and downstream reach at T 1 should be equivalent to the ratio between the condition of the same upstream reach and same downstream reach at T 2 .

O upstream reach
where O refers to the observed value and E refers to the expected value.If we let the downstream reach represent the restoration reach, we can provide the values for the pre-restoration conditions (upstream and downstream reach values at T 1 ) and the value for the upstream reach at T 2 and solve for the "expected" value of the downstream reach (i.e., restoration reach) at T 2 .We can then use the observed downstream reach value at T 2 , to calculate the percent change from the expected value that we can attribute to the restoration activities: Change attributed to restoration (%) = (O downstream reach T2 − E downstream reach T2 ) E downstream reach T2 * 100 (2) Calculating change as a function of the expected value can help account for variability between the images due to natural variability in climate and stream discharge.Although the outputs used to analyze changes in stream surface area tended to show high accuracy, if we assume that error induced by the methodology is random or consistent across the image extent, this approach also allowed us to take into account between-image variability in the accuracy of the mapped stream extent.We also recognize that changes to the restoration reach can potentially propagate upstream (e.g., [67]), possibly influencing the relationship between the upstream and downstream reaches; however, due to image extents we were limited in how far upstream we could document stream surface area necessitating the above assumptions.
In addition to calculating changes in stream surface area at a reach scale, we also calculated changes in stream width, inundated length, and changes to stream surface area just upstream from each BDA.To calculate changes in stream width, 20 points were randomly selected along the stream centerline within each site and reach and stream width was measured manually using the pan-sharpened raw imagery and averaged to obtain a single mean stream width per reach.The same points were used in both years.Inundated stream length was calculated as the percent of the stream centerline mapped as water.To quantify local changes in stream surface area induced by the structures, we calculated the stream surface area immediately upstream of each of the structures that showed a visually evident change (from 2014 to 2017) in stream width.Stream surface area was calculated from the classified (edited) stream surface area used to quantify reach-scale changes.The local change in stream surface area was presented as the change relative to pre-restoration stream surface area.The localized change was observed over a variable stream length distance, but the observed effect averaged 26 m of stream length upstream from the installed BDAs.

Changes in Riparian Condition
Riparian greenness was evaluated using three vegetation indices, NDVI, EVI and the Soil-Adjusted Vegetation Index (SAVI) [65].Multiple indices were included because vegetation indices are sensitive to several conditions including canopy geometry (trees versus herbaceous vegetation), soil properties, sun position, and cloudiness [68][69][70].While NDVI is the most commonly used vegetation index, EVI can help eliminate atmospheric noise, and SAVI reduces the influence of soil by including a soil adjustment factor (L) [56,68].Changes in riparian greenness, averaged across the three indices, were then evaluated, (1) along the restoration reach (from the upstream BDA to the downstream BDA); (2) 0 m to 250 m downstream from the restoration reach; (3) 250 m to 500 m downstream from the restoration reach; and (4) 500 m to 1 km downstream from the restoration reach.At Robb Creek, the BDAs were designed to reactivate side channels so the side channels were considered separately from the main channel.We also tested how the effect changed as the buffer from the main channel increased from 10 m to 20 m from the channel centerline.
It was critical to control between-image differences in greenness not related to the restoration, therefore changes in greenness between the pre-and post-restoration riparian corridors were corrected using the difference in greenness across reference areas.Three reference polygons were selected at each site (ranging from 0.2 ha to 2.5 ha in size) representing: (1) herbaceous photosynthetic riparian vegetation upstream from the restoration site; (2) an upland patch dominated by photosynthetic grasses; and (3) an upland patch dominated by non-photosynthetic grasses.We avoided areas that appeared to show a difference in grazing intensity between the two image dates.The greenness values for each of the reference polygons were averaged to obtain a reference greenness value for each site and date.Because herbaceous vegetation is more sensitive to interannual change than riparian tree species, the reference polygons included a mix of riparian and upland vegetation samples.Additionally, because changes may have occurred for upstream vegetation in response to the restoration measures, we used a mix of riparian and non-riparian patches.Even at sites with willow plantings, short-term changes (one to three years) in riparian condition can be expected to primarily result from growth of herbaceous species, a change attributable to increased water availability in the shallow subsurface areas adjacent to the stream.The "reference greenness" values were seen as equivalent to the role of the upstream reach values when analyzing changes to stream surface area.The same analysis used to evaluate restoration effects on stream surface area above (Equations ( 1) and ( 2)) were applied here to evaluate restoration-induced changes to riparian condition.

Accuracy of Stream Delineation Approaches
The accuracy of our estimates of stream surface area depended on both the stream width and the classification approach.Across all sites and classification approaches, the relative bias tended to be negative, indicating that the stream surface area, on average, was underestimated.Long Creek (South), which showed a stream width averaging 3 m to 3.5 m and Long Creek (North), which showed a stream width averaging 2.5 to 3.5 m wide, showed more consistent accuracy statistics across the approaches tested relative to the other two sites.The eCognition output and Panchromatic brightness consistently performed the best with errors of omission for water ranging from 2% to 6% and errors of commission ranging from 0.5% to 7%.In contrast, the NDWI tended to show higher errors of omission, with omission errors ranging from 7% to 37% and commission errors ranging from 1% to 5% (Table 3).Overall accuracy and dice coefficients were >95% except for NDWI outputs and one of the four WWI outputs (Table 3).
As stream width decreased at the Alkali Creek (averaged 1 m to 2.5 m) and Robb Creek (averaged 1.5 m to 2 m) sites, the accuracy of published indices (NDWI and WWI) was relatively poor.Using these indices errors of omission, for instance, ranged from 16% to 69% along the two streams.However, panchromatic brightness and eCognition outputs maintained relatively strong accuracy even as stream width decreased with errors of omission and commission across the two sites and years ranging from 3% to 11% and 2% to 13%, respectively, and the corresponding overall accuracy and Dice coefficient ranging from 91% to 97%.
It was also evident that classification accuracy for a given method can be inconsistent over time.Examples of this are shown in Figure 6, in which we can compare the outputs for each classification approach and year along the Alkali Creek site.Figure 6 also demonstrated that when surface water is a minority cover type across the image extent, an output can visually appear to be relatively noisy but statistically show a relatively low amount of calculated commission error.An example of this is the panchromatic brightness output for 2014 in which an error of commission of 13% created a visually "noisy" output.The visualized variability in the surface-water extent across methods justifies the need to analyze changes in stream surface area using the most accurate method possible and including a manual editing component as time allows so that uncertainty in surface-water extent does not obscure "true changes" in stream condition.

Changes to Stream Condition
The installation of BDAs resulted in proximal changes (e.g., increases in surface water immediately upstream from the in-stream structures) as well as changes at the scale of the restoration reach and downstream reaches.Increases in stream surface area immediately upstream from structures or in reactivated side channels were observed at all four sites.Increases in total stream surface area along the restoration reach were observed at three of the four sites, while a decrease in stream surface area downstream from the restoration reach was observed at all four sites (Table 4).

Changes to Stream Condition
The installation of BDAs resulted in proximal changes (e.g., increases in surface water immediately upstream from the in-stream structures) as well as changes at the scale of the restoration reach and downstream reaches.Increases in stream surface area immediately upstream from structures or in reactivated side channels were observed at all four sites.Increases in total stream surface area along the restoration reach were observed at three of the four sites, while a decrease in stream surface area downstream from the restoration reach was observed at all four sites (Table 4).At Alkali Creek approximately one year post-restoration, we observed a net increase of 303 m 2 in stream surface area immediately upstream from installed BDAs (Figure 7).This amount of water represented 19% of the total stream surface area along the restoration reach in 2017.After controlling for differences in stream surface area attributable to interannual variability, we observed a 25% increase in stream surface area attributable to the restoration activities, and a 30% decrease in stream surface area downstream from the restoration reach (Table 4).These changes in stream surface area were matched by corresponding changes in stream width (Table 4).A complicating factor at this site was that beaver activity was present from 200 m to 800 m upstream from the restoration reach.Upstream ponding resulting from natural beaver dams was substantial in the 2014 image but total stream surface water along the beaver impacted reach (200 m to 800 m upstream) decreased 74% by 2017 (Table 4, Figure 8A).At Alkali Creek approximately one year post-restoration, we observed a net increase of 303 m 2 in stream surface area immediately upstream from installed BDAs (Figure 7).This amount of water represented 19% of the total stream surface area along the restoration reach in 2017.After controlling for differences in stream surface area attributable to interannual variability, we observed a 25% increase in stream surface area attributable to the restoration activities, and a 30% decrease in stream surface area downstream from the restoration reach (Table 4).These changes in stream surface area were matched by corresponding changes in stream width (Table 4).A complicating factor at this site was that beaver activity was present from 200 m to 800 m upstream from the restoration reach.Upstream ponding resulting from natural beaver dams was substantial in the 2014 image but total stream surface water along the beaver impacted reach (200 m to 800 m upstream) decreased 74% by 2017 (Table 4, Figure 8A).The goal of the restoration along Robb Creek was to reactivate two side channels.At approximately two years post-restoration, we observed a reactivation of the side channels that included a 78% increase in stream surface area along the northern side channel and a 356% increase in stream surface area along the southern side channel (Table 4, Figure 8).Correspondingly, we observed a 27% decrease in stream surface area along the main stream stem, and a decrease in stream surface area downstream from the restoration reach (0 m to 500 m downstream) (Table 4).
Long Creek (South) is 4.8 stream km downstream from Long Creek (North) and was the first restored of the four sites.At the time of the post-restoration image the site was approximately three years post-restoration.We observed a substantial amount of water stored upstream of many of the structures (net increase of water surface area of 747 m 2 ) (Figure 9), which represented 6% of the total stream surface area along the restoration reach in 2017.We observed a 14% increase in stream surface area along the restoration reach and a corresponding decrease in stream surface area (−17% to −20%) from 0 m to 500 m downstream from the restoration reach (Table 4).Long Creek (North) at one-year post-restoration was the only site where we did not observe an increase in stream surface area along the restoration reach, but instead observed a minor decrease of 5% in stream surface area.We also observed the smallest increase in stream surface area upstream from the BDA structures (net increase of proximal water surface area of 171 m 2 ), relative to the other sites evaluated (Figure 9).Similar to other sites, however, Long Creek (North) showed a decrease in stream surface area (−9% to −16%) through 1 km downstream from the restoration reach and an associated decrease in mean stream width (Table 4, Figure 8).The goal of the restoration along Robb Creek was to reactivate two side channels.At approximately two years post-restoration, we observed a reactivation of the side channels that included a 78% increase in stream surface area along the northern side channel and a 356% increase in stream surface area along the southern side channel (Table 4, Figure 8).Correspondingly, we observed a 27% decrease in stream surface area along the main stream stem, and a decrease in stream surface area downstream from the restoration reach (0 m to 500 m downstream) (Table 4).
Long Creek (South) is 4.8 stream km downstream from Long Creek (North) and was the first restored of the four sites.At the time of the post-restoration image the site was approximately three years post-restoration.We observed a substantial amount of water stored upstream of many of the structures (net increase of water surface area of 747 m 2 ) (Figure 9), which represented 6% of the total stream surface area along the restoration reach in 2017.We observed a 14% increase in stream surface area along the restoration reach and a corresponding decrease in stream surface area (−17% to −20%) from 0 m to 500 m downstream from the restoration reach (Table 4).Long Creek (North) at one-year post-restoration was the only site where we did not observe an increase in stream surface area along the restoration reach, but instead observed a minor decrease of 5% in stream surface area.We also observed the smallest increase in stream surface area upstream from the BDA structures (net increase of proximal water surface area of 171 m 2 ), relative to the other sites evaluated (Figure 9).Similar to other sites, however, Long Creek (North) showed a decrease in stream surface area (−9% to −16%) through 1 km downstream from the restoration reach and an associated decrease in mean stream width (Table 4, Figure 8).

Changes to Riparian Condition
At each site, the spectral greenness index used (NDVI, EVI or SAVI) showed a relatively minor but inconsistent influence on the reported change in greenness, after controlling for between-year differences in greenness.EVI, for example, showed a higher percent change in greenness relative to NDVI and SAVI along Alkali Creek and Robb Creek, but not along the Long Creek sites (Table 5).However, we also found that the directionality of change to riparian condition post-restoration at each of the sites did not depend on the greenness index used.Increases in riparian greenness along the restoration reach were observed at three of the four sites, while changes to riparian greenness downstream from the restoration reach were less consistent.

Changes to Riparian Condition
At each site, the spectral greenness index used (NDVI, EVI or SAVI) showed a relatively minor but inconsistent influence on the reported change in greenness, after controlling for between-year differences in greenness.EVI, for example, showed a higher percent change in greenness relative to NDVI and SAVI along Alkali Creek and Robb Creek, but not along the Long Creek sites (Table 5).However, we also found that the directionality of change to riparian condition post-restoration at each of the sites did not depend on the greenness index used.Increases in riparian greenness along the restoration reach were observed at three of the four sites, while changes to riparian greenness downstream from the restoration reach were less consistent.
Table 5. Percent change in riparian greenness between pre-(2014) and post-restoration (2017) images after controlling for between-image differences in greenness not related to the restoration action.Greenness was evaluated using the Soil-Adjusted Vegetation Index (SAVI), Enhanced Vegetation Index (EVI) and Normalized Difference Vegetation Index (NDVI).Indices were evaluated using a 10 m, 15 m, and 20 m buffer from the stream.Changes along the restored reaches are in bold.DS: downstream; N: North; S: South.Along Alkali Creek, we observed an average increase in greenness of 20% using a 10 m stream buffer, declining to a 15% increase using a 20 m stream buffer (Table 5).Although we had observed a decrease in stream surface area downstream from this restoration site, we found an increase in greenness that persisted through 1 km downstream from the restoration reach.Defining reference conditions required particular attention along Alkali Creek because the post-restoration image was late summer (August 2, 2017) when water is more limited.This is evident in the contrast observed in the NDVI between the riparian area and uplands in the 2017 Alkali Creek image (Figure 10).Along Robb Creek, minor increases in greenness were observed along the main stem and northern side channel (<5%), while a substantial green-up was observed along the southern side channel that showed an increase in greenness of 21% (Figure 10).The decrease in stream surface area observed downstream of this restoration site was found to co-occur with a decrease in greenness, which ranged from a 4% to 12% decrease in greenness through 1 km downstream of the restoration (Table 4).Although Long Creek (North) showed a minor decrease in stream surface area within the restoration reach, we observed an increase in greenness of 7 to 8% along the restoration reach and an increase in greenness of 13 to 15% just downstream of the restoration reach (0 to 250 m) (Figure 11).In contrast, while we observed clear changes in stream surface area along Long Creek (South), changes in greenness were minimal (<5% change) (Table 5, Figure 11).
Along Alkali Creek, we observed an average increase in greenness of 20% using a 10 m stream buffer, declining to a 15% increase using a 20 m stream buffer (Table 5).Although we had observed a decrease in stream surface area downstream from this restoration site, we found an increase in greenness that persisted through 1 km downstream from the restoration reach.Defining reference conditions required particular attention along Alkali Creek because the post-restoration image was late summer (August 2, 2017) when water is more limited.This is evident in the contrast observed in the NDVI between the riparian area and uplands in the 2017 Alkali Creek image (Figure 10).Along Robb Creek, minor increases in greenness were observed along the main stem and northern side channel (<5%), while a substantial green-up was observed along the southern side channel that showed an increase in greenness of 21% (Figure 10).The decrease in stream surface area observed downstream of this restoration site was found to co-occur with a decrease in greenness, which ranged from a 4% to 12% decrease in greenness through 1 km downstream of the restoration (Table 4).Although Long Creek (North) showed a minor decrease in stream surface area within the restoration reach, we observed an increase in greenness of 7 to 8% along the restoration reach and an increase in greenness of 13 to 15% just downstream of the restoration reach (0 to 250 m) (Figure 11).In contrast, while we observed clear changes in stream surface area along Long Creek (South), changes in greenness were minimal (<5% change) (Table 5, Figure 11).

Discussion
Long-term trends in the degradation of riparian and stream habitat are common across the western United States [71,72].Satellite imagery has the potential to provide spatially continuous monitoring of stream extent and condition, which can complement point-based field efforts and stream gage data, and better inform stream management in response to degradation.However, only a limited number of studies have attempted to apply satellite imagery to streams, particularly smaller streams (<5 m wide).We found that pan-sharpened high-resolution imagery can be used to effectively monitor streams as narrow as 1.5 m wide.We tested sites where stream width was <1 m and found the results too poor to include.Panchromatic brightness consistently outperformed more established indices such as the NDWI.However, in experimenting with the classification of stream surface area we found that spectrally mixed portions of a stream can be more challenging to identify using this approach, for instance, portions of the stream showing high turbidity, high chlorophyll levels or bright sands can be missed using brightness alone.In addition, vegetation shadows, common in riparian areas, can be erroneously mapped as inundated areas, creating substantial errors of commission.Alternatively, while far more time consuming, a GEOBIA approach was able to greatly reduce errors of commission outside of the stream area by increasing object size with distance from stream and adding region growing to the shadow class within heavily shadowed riparian areas.While the index-and object-based approaches showed distinct advantages, we found that regardless of efforts accuracy results could be uneven across years.This is a major challenge in change detection analysis in which uneven accuracy over time, potentially attributable to differences in sensors, the off-nadir view angle, as well as the time of day that the image was collected or variability in local hydro-climatic conditions, can obscure change attributable to human-caused degradation or restoration [73].We sought to minimize this source of error by (1) using only our most accurate outputs and further manually editing these outputs prior to analyzing changes in stream surface area; and (2) calculating the change as a function of the change from the expected value, therefore controlling, to the extent possible, for change due to image quality or variability in hydro-climatic conditions.This aspect of remote sensing change analysis, however, remains a challenge.

Discussion
Long-term trends in the degradation of riparian and stream habitat are common across the western United States [71,72].Satellite imagery has the potential to provide spatially continuous monitoring of stream extent and condition, which can complement point-based field efforts and stream gage data, and better inform stream management in response to degradation.However, only a limited number of studies have attempted to apply satellite imagery to streams, particularly smaller streams (<5 m wide).We found that pan-sharpened high-resolution imagery can be used to effectively monitor streams as narrow as 1.5 m wide.We tested sites where stream width was <1 m and found the results too poor to include.Panchromatic brightness consistently outperformed more established indices such as the NDWI.However, in experimenting with the classification of stream surface area we found that spectrally mixed portions of a stream can be more challenging to identify using this approach, for instance, portions of the stream showing high turbidity, high chlorophyll levels or bright sands can be missed using brightness alone.In addition, vegetation shadows, common in riparian areas, can be erroneously mapped as inundated areas, creating substantial errors of commission.Alternatively, while far more time consuming, a GEOBIA approach was able to greatly reduce errors of commission outside of the stream area by increasing object size with distance from stream and adding region growing to the shadow class within heavily shadowed riparian areas.While the index-and object-based approaches showed distinct advantages, we found that regardless of efforts accuracy results could be uneven across years.This is a major challenge in change detection analysis in which uneven accuracy over time, potentially attributable to differences in sensors, the off-nadir view angle, as well as the time of day that the image was collected or variability in local hydro-climatic conditions, can obscure change attributable to human-caused degradation or restoration [73].We sought to minimize this source of error by (1) using only our most accurate outputs and further manually editing these outputs prior to analyzing changes in stream surface area; and (2) calculating the change as a function of the change from the expected value, therefore controlling, to the extent possible, for change due to image quality or variability in hydro-climatic conditions.This aspect of remote sensing change analysis, however, remains a challenge.Additionally, as the method that produced the highest overall accuracy varied also across sites, it was evident that the appropriate processing approach to minimize uneven error over time will vary across sites depending on the amount and type of riparian vegetation (herbaceous or tree), stream width, as well as water depth and clarity, a finding supported by others [74,75].
Improving techniques to map stream surface area with commercial high-resolution imagery offers opportunities to remotely monitor changes in key aspects of stream condition induced by flood or drought events, shifts in local land uses, or in-stream restoration activities.However, it is important to clarify that the stream data gleaned from a remote sensing analysis is intrinsically different from the data a hydrologist typically uses.For example, at the Alkali Stream site, differences in the seasonality of the image pair (June vs. August) meant that the downstream discharge in the Jefferson River was much lower at the August date relative to the June date (15 m 3 s −1 compared to 127 m 3 s −1 ); however, the stream surface area upstream from the restoration reach along Alkali Creek was only 3% less at the August date relative to the June date.This contrast clarifies that a remote sensing analysis is not necessarily capturing changes in stream discharge, which are better measured with stream gages, but instead providing a spatially continuous dataset of changes in stream surface area, specifically stream width and the creation or change to riparian wetlands, which could in turn, impact downstream stream discharge.
Relying on image pairs, however, or only two points in time, can limit our understanding regarding seasonally specific effects.For example, local stakeholders are interested in the effects of BDAs on streams not just after snowmelt, but in particular during the late summer period in which water availability can be limited [27,31].Because very few DigitalGlobe images have been collected and archived across southwestern Montana we were restricted to the early summer period and were therefore unable to evaluate the impact of the BDAs during this late summer period.This limitation in image timing means we were unable to observe how changes in stream surface area or riparian condition documented near the start of the growing season influenced conditions near the end of the growing season.However, this limitation will likely be reduced in the near future.Sentinel-2 (10-20 m resolution), launched in June 2015, was too coarse for the streams evaluated, but in the future, could be used to regularly monitor the condition of rivers >10 m wide.CubeSats, such as those launched by Planet (San Francisco, CA, USA), also show high potential for improved monitoring of stream condition at more frequent intervals.Obstacles to the widespread application of CubeSats for monitoring stream condition, however, include the limited number of spectral bands (blue, green, red, near-infrared), the cost of Planet imagery, as well as challenges in calibrating reflectance and georeferencing between satellites e.g., [76].As these technical obstacles are overcome, satellite imagery can be more commonly used to monitor streams in a spatially continuous manner.
We found separating the influence of weather relative to the influence of human-induced change particularly challenging in the riparian areas.Trends in riparian condition could be very sensitive to how reference conditions were defined.In part, this is because we might expect that the magnitude of interannual variability may be inconsistent across areas dominated by trees versus herbaceous vegetation.Additionally, reference areas can be influenced by forces of change independent of weather patterns such as variability in grazing intensity.The use of image pairs, and the time-frame at which the projects were considered (one to three years post-restoration), also limited our ability to evaluate the success of the activities with a lagged response, such as the tree planting effort, which takes more time to establish and can act to further influence riparian-stream interactions over time [44,77].Riparian and wetland herbaceous species, however, can respond very quickly to changes in riparian soil moisture [71,77].Therefore, changes in riparian greenness can be used to indirectly assess short-term changes to the stream hydrologic processes [77][78][79].Because changes in riparian greenness were detected even when using a 20 m stream buffer, it may be possible to monitor riparian condition for streams <5 m wide using moderate resolution satellites, such as Landsat or Sentinel-2.As riparian trees grow larger and become more established, they can also impact our ability to monitor stream surface area.For instance, stream length at most sites was found to be <100% inundated, a finding that can be attributed to overhanging vegetation or hyporheic flow masking stream water.
Because our analysis relied on just two high-resolution images per site, we view this analysis as an example of the capabilities and potential of using high-resolution imagery to monitor stream condition and not a conclusion on the impact of specific restoration methods (e.g., BDAs, riparian planting) on streams.However, despite these concerns, we saw proximal evidence of increases in stream surface area upstream from structures at all four sites and evidence of either increases in stream surface area or riparian greenness along the restoration reach at all four sites.We also saw a decrease in stream surface area downstream from the restoration reach at all four sites.These findings suggest that the restoration activities have induced increases in instream water storage, at least during the early summer period.Uneven findings regarding the change induced by BDAs could be due to several factors including stream size, pre-restoration conditions [41], restoration goals, restoration installation (including extent and age), and site specifics, all of which can influence how a stream responds to change through time.However, these findings suggest that high-resolution imagery can provide a spatially continuous understanding of how narrow (<5 m wide) streams respond to restoration projects.

Conclusions
As high-resolution, multispectral imagery becomes more frequently collected and available, regular monitoring of stream surface area and condition in response to local or watershed-based changes will become increasingly feasible.However, mapping narrow, linear water features that are subject to rapid changes in water depth and turbidity may require approaches independent from those widely applied to map wetlands and lakes [1,2,80].We found that utilizing bands unique to Worldview-2 and 3 in an eCognition framework can produce accurate results and minimal errors of commission down to a stream width of approximately 1.5 m.By comparing images pre-and postrestoration across multiple sites in the Upper Missouri River Headwaters Basin, we were able to quantify proximal and reach-scale changes in stream surface area and riparian greenness in response to the installation of multiple BDAs at each site.In general, the installation of BDAs appeared to create an increase in stream surface area immediately upstream from many of the structures.These proximal changes tended to result in reach-scale increases in stream surface area and riparian greenness along the restoration reach as well as decreases in stream surface area for reaches just downstream (through 500 m) from the restoration reach at most of the sites.The consistency of the directional changes to stream surface area across the sites, despite differing patterns in discharge lower in the watershed between-image pairs, suggests that we were able to account for hydro-climatic variability.Restoring degraded streams can positively influence in-stream habitat, water quality as well as water quantity across a watershed [34,37,38].Monitoring the impacts of stream restoration projects, including resource-efficient structures such as BDAs, can improve and inform site selection and expectations for future stream restoration efforts.

Figure 1 .
Figure 1.Distribution of the stream sites within the Upper Missouri Headwaters Basin.Background image is a Landsat 8 image (path 39, row 39, June 9, 2016).Location of the U.S. Geological Survey stream gage (Jefferson River, #06026500) is also shown relative to the restoration sites.NHD: National Hydrography Dataset, HUC8s: 8-digit Hydrological Units.

Figure 1 .
Figure 1.Distribution of the stream sites within the Upper Missouri Headwaters Basin.Background image is a Landsat 8 image (path 39, row 39, 9 June 2016).Location of the U.S. Geological Survey stream gage (Jefferson River, #06026500) is also shown relative to the restoration sites.NHD: National Hydrography Dataset, HUC8s: 8-digit Hydrological Units.

Figure 2 .
Figure 2. (A) The monthly Palmer Hydrological Drought Index (PHDI) values for southwestern Montana over the past 10 years, converted to percentages based on the historical record (1895-2017).Stars indicate the PHDI value at the time of the pre-and post-restoration images; (B) The daily mean discharge for the Jefferson River near Twin Bridges, Montana (USGS Gage #06026500) downstream of the four restoration sites.Diamonds indicate the discharge value at the time of the pre-and postrestoration images.

Figure 2 .
Figure 2. (A) The monthly Palmer Hydrological Drought Index (PHDI) values for southwestern Montana over the past 10 years, converted to percentages based on the historical record (1895-2017).Stars indicate the PHDI value at the time of the pre-and post-restoration images; (B) The daily mean discharge for the Jefferson River near Twin Bridges, Montana (USGS Gage #06026500) downstream of the four restoration sites.Diamonds indicate the discharge value at the time of the pre-and post-restoration images.

29 Figure 3 .
Figure 3.A schematic of (A) Alkali Creek, (B) Long Creek (North and South), and (C) Robb Creek, indicating the reaches along each of the stream sites that were analyzed, relative to the distribution of the beaver dam analogues.Flow direction can be determined from the relative location of the upstream and downstream reaches.

Figure 3 .
Figure 3.A schematic of (A) Alkali Creek, (B) Long Creek (North and South), and (C) Robb Creek, indicating the reaches along each of the stream sites that were analyzed, relative to the distribution of the beaver dam analogues.Flow direction can be determined from the relative location of the upstream and downstream reaches.

Figure 4 .
Figure 4. Examples from Robb Creek, Montana, of segmenting an image into objects of interest.Segmentation approaches that focus on detecting edges as well as iterative segmentation at different scales can allow objects to vary in size from large, upland objects (A) to narrow, small objects that follow the stream (B); DigitalGlobe Copyright 2017.

Figure 4 .
Figure 4. Examples from Robb Creek, Montana, of segmenting an image into objects of interest.Segmentation approaches that focus on detecting edges as well as iterative segmentation at different scales can allow objects to vary in size from large, upland objects (A) to narrow, small objects that follow the stream (B); DigitalGlobe Copyright 2017.

Figure 5 .
Figure 5.A flowchart showing the order of steps taken to process each image to water and non-water using an object-based approach.Inputs and outputs are shown in ovals while processing steps are shown in rectangles.

Figure 5 .
Figure 5.A flowchart showing the order of steps taken to process each image to water and non-water using an object-based approach.Inputs and outputs are shown in ovals while processing steps are shown in rectangles.

Figure 6 .
Figure 6.A comparison of mapped surface water extent for a reach of Alkali Creek in which challenges included shaded riparian vegetation (left) and narrow stretches of stream (center).Comparisons include, natural color images (A,F); the Worldview Water Index (WWI, B,G); the Normalized Difference Wetness Index (NDWI, C,H); the panchromatic brightness (D,I); and the eCognition output (E,J).O: omission error; C: commission error Copyright DigitalGlobe, 2014, 2017.

Figure 6 .
Figure 6.A comparison of mapped surface water extent for a reach of Alkali Creek in which challenges included shaded riparian vegetation (left) and narrow stretches of stream (center).Comparisons include, natural color images (A,F); the Worldview Water Index (WWI, B,G); the Normalized Difference Wetness Index (NDWI, C,H); the panchromatic brightness (D,I); and the eCognition output (E,J).O: omission error; C: commission error Copyright DigitalGlobe, 2014, 2017.

Figure 7 .
Figure 7. Retention of water upstream of beaver dam analogues along Alkali Creek are visible by comparing pre-restoration stream reaches (A,C) with post-restoration stream reaches (B,D).Copyright DigitalGlobe 2014, 2017.

Figure 7 .
Figure 7. Retention of water upstream of beaver dam analogues along Alkali Creek are visible by comparing pre-restoration stream reaches (A,C) with post-restoration stream reaches (B,D).Copyright DigitalGlobe 2014, 2017.

Figure 8 .
Figure 8. Change in stream surface area between the most upstream and most downstream beaver dam analogues (BDAs) along (A) Alkali Creek; (B) Long Creek (North); (C) Long Creek (South); and (D) Robb Creek.Copyright DigitalGlobe, 2014.

Figure 8 .
Figure 8. Change in stream surface area between the most upstream and most downstream beaver dam analogues (BDAs) along (A) Alkali Creek; (B) Long Creek (North); (C) Long Creek (South); and (D) Robb Creek.Copyright DigitalGlobe, 2014.

29 Figure 9 .
Figure 9. Visually apparent changes with the installation of the beaver dam analogues varied along Long Creek.Within Long Creek (North) changes from pre-to post-restoration included retention of water upstream from structures as well as the reactivation of abandoned side channels (A,B); Along Long Creek (South) changes included widening of the stream as well as retention of water upstream from structures (C,D).Copyright DigitalGlobe 2014, 2017.

Figure 9 .
Figure 9. Visually apparent changes with the installation of the beaver dam analogues varied along Long Creek.Within Long Creek (North) changes from pre-to post-restoration included retention of water upstream from structures as well as the reactivation of abandoned side channels (A,B); Along Long Creek (South) changes included widening of the stream as well as retention of water upstream from structures (C,D).Copyright DigitalGlobe 2014, 2017.

Table 2 .
Spectral indices used in object-based and pixel-based classification of surface water, as well as the characterization of riparian greenness.NIR: near infrared.

Table 4 .
Change to stream surface area.Local storage increase refers to a net change in the surface area of water upstream from individual beaver dam analogues (BDAs).Because the BDAs at Robb Creek served a primary purpose of redirecting water toward side channels, local storage is not shown.