Linking Hydrologic and Hydraulic Data with Models to Assess Flow and Channel Alteration at Hog Park, Wyoming USA

: Transbasin diversions and dams allow for water uses when and where there is high demand and low supply, but can come with an expense to the environment. This paper presents a linkage of hydrologic and hydraulic modeling and datasets to assess the hydrologic and hydraulic stability within a transbasin watershed as an approach for meeting water use targets and safeguarding environmental sustainability. The approach used a Prediction in Ungauged Basin (PUB) regionalization technique that completed the parameterization of a study watershed hydrologic model by transferring calibrated parameters from a reference watershed hydrologic model. This resulted in a long-term, simulated natural ﬂow record that was compared to the measured modiﬁed ﬂow record for the same time period to assess ﬂow alteration. In the sensitive reach, hydraulic modeling results tracked channel response from before hydrologic modiﬁcation to baseline using repeated survey years during the hydrologic modiﬁcation. The combined assessment of hydrology and hydraulics highlighted the relation between ﬂow regime and channel form. deﬁned for peak ﬂows (no ﬂood, small ﬂood, and large ﬂood), intermediate ﬂows, and low ﬂows. The x -axis provides an indication of the ﬂow attribute timing and the y -axis provides an indication of the ﬂow attribute magnitude. Following a stable, modiﬁed ﬂow regime for a long-term period would be one tool to promote future channel stabilization.


Introduction
Across the world, transbasin diversions move large volumes of water from one river basin to another, and major dams impound and store water for later use [1]. These large-scale water collection systems allow for the beneficial use of water when and where there is high demand and low supply. However, the benefits may come at the expense of the environment, and the extent is not fully understood [2]. Though environmental analyses for such projects may have been completed to minimize potential adverse effects, some were difficult to predict and persisted many decades later [3][4][5]. As the supply and demand for water changes, whether due to climate change, land use change, or population growth, there is a pressing need to better understand the balance of vital water resource objectives and environmental sustainability.
The environmental effects from extensive water collections systems are influenced by the changes in hydrology. The changes in hydrology occur as a result of water being augmented to or diverted from a basin, stored behind dams, and released at different rates, magnitudes, and times throughout the year. This study analyzes the scenario at Hog Park, Wyoming ( Figure 1) where a major dam stores the water received from a transbasin diversion. One result is downstream channel enlargement [6].
Transbasin diversions have been found to increase intermediate and low flows in the receiving basin, such as the Kemano River in central British Columbia that receives water from the Nechako River for hydroelectric power generation [7], and the Milk River in Montana, which receives water from the Saint Mary River in Alberta [8]. Similarly, a 10-fold increase in low flows was observed at the River Ter, UK, which receives water from groundwater pumped to Leighs Reservoir [9]. These studies indicated that increased low and intermediate flows can play a key role in channel enlargement. On smaller rivers, peak flows can play a role in channel enlargement due to flow alteration from transbasin diversion and dam storage, such as diversion from the Colorado River headwaters to La Poudre Pass Creek of the Poudre River where a combination of increased peak Transbasin diversions have been found to increase intermediate and low flows in the receiving basin, such as the Kemano River in central British Columbia that receives water from the Nechako River for hydroelectric power generation [7], and the Milk River in Montana, which receives water from the Saint Mary River in Alberta [8]. Similarly, a 10-fold increase in low flows was observed at the River Ter, UK, which receives water from groundwater pumped to Leighs Reservoir [9]. These studies indicated that increased low and intermediate flows can play a key role in channel enlargement. On smaller rivers, peak flows can play a role in channel enlargement due to flow alteration from transbasin diversion and dam storage, such as diversion from the Colorado River headwaters to La Poudre Pass Creek of the Poudre River where a combination of increased peak flows and decreased low flows contribute to channel enlargement [10]. Peak flows were double, with longer duration on the Lake Creek and the Lake Fork of the Arkansas River, which receive water from the Roaring Fork of the Colorado River, resulting in more than fifty years of channel enlargement, as well as substantially decreased riparian cover area [11]. In the Owens River of California, which receives water from the nearby Mono Basin, increased flood magnitude, frequency, and duration resulted in a decrease in annual growth rates of willows, with the best willow growth rates spurred from floods in the high range of natural flooding and a low range of the modified flooding [12]. Together, these studies indicate that changes in hydrology play a key role in channel enlargement, as well as the associated loss of riparian ecosystems.
The changes in hydrology have been estimated using a framework known as predictions in ungauged basins (PUB). PUB utilizes a regionalization of hydrologic process knowledge to understand hydrologic response in ungauged basins [13]. Regional regression methods can be used to estimate ungauged, natural flow characteristics. The regional regression equations derived from the USGS flood flow frequency analysis method [14] are published for most regions in the US, including the mountainous regions of Wyoming, but have a high standard error [15]. Locally developed multiple regression models have demonstrated the potential for accurate estimates of natural flow regime metrics [16]. Reference sites that have spatially and physically similar proxies and have a long-term, measured, natural flow record can be used to estimate natural flow attributes [17]. Furthermore, the attributes of natural hydrologic processes from a reference watershed can be transferred to a hydrologically modified study watershed based on the close spatial proximity and physical similarity [18]. One application of the latter regionalization technique is to parameterize a hydrologic model of a hydrologically modified study watershed to simulate a natural flow record [19,20].
In response to a change in hydrology, a channel changes cross-sectional, plan, and profile geometry. The progression of this channel change can be assessed using measured channel geometry data for a long-term period of time and hydraulic modeling. A hydraulic model simulates flow characteristics in a channel or through a network of multiple channels [21]. For hydraulic properties, such as water surface elevations, wetted perimeters, and flood inundation area, the hydraulic model uses the measured channel geometry data [22]. Hydraulic models of sequential years show the progression of hydraulic changes. An improved method to concurrently assess hydrologic and hydraulic stability could provide a tool to inform managers on how to simultaneously meet water resource targets and safeguard environmental sustainability. Towards this purpose, the study objectives are: (1) use PUB methods to simulate long-term natural flow in a hydrologically modified watershed without a measured, natural flow record, but with a measured, modified flow record; (2) use channel geometry data from sequential years to model the long-term hydraulic response of a sensitive reach within the hydrologically modified watershed; and (3) assess relations between the flow regime and channel response in a hydrologically modified watershed. The unique constraints that lead to this type of assessment are: (1) limited long-term flow data to understand flow regime alteration in hydrologically modified watersheds, and (2) limited long-term channel geometry data to understand channel response in hydrologically modified watersheds. This approach recognizes that establishing a new long-term, modified flow regime and channel form that are both in harmony may be more practical than returning hydrologically modified systems to their natural state.

Study Site
The Hog Park study area is in the Sierra Madre Range in the Southern Rocky Mountains (Figure 1). Glaciation, fluvial processes, climate, and dam management helped shape the area. Moderate to deep loamy-skeletal soils and steep, forested hillslopes are underlain by a Precambrian igneous and metamorphic geology [23,24]. The climate consists of cold winters and cool summers where the bulk of precipitation falls as snow, with a median peak snow water equivalent of about 760 mm and an average annual precipitation of about 1080 mm, as measured at the nearby USDA NRCS Whiskey Creek snow telemetry (SNOTEL) site [25].
The Hog Park study site encompasses a reference watershed, study watershed, and study watershed sensitive reach ( Figure 1). The reference watershed is the Encampment River drainage above its confluence with Hog Park Creek and is considered mostly pristine with minor anthropogenic influences due to its location in designated wilderness [26]. The study and reference watersheds are in close spatial proximity and consequently share similar climate, geology, and land cover. Notable differences in the study watershed are a large reservoir with a storage capacity of 3.08 × 10 7 m 3 and large transbasin diversion which receives water piped through the Continental Divide from 126 diversion structures on the North Fork of the Little Snake River. Other differences include drainage area, aspect, elevation, and percent land cover types (Table 1). A sensitive reach is located within a large, unconfined valley downstream of the main dam on Hog Park Creek ( Figure 1).

Methodology
To understand the effects of hydrologic modification in the study watershed, hydrologic and hydraulic data were linked with models. The first step was collating the overlapping long-term hydrologic and hydraulic data that were used to parameterize hydrologic and hydraulic models. The Precipitation-Runoff Modeling System (PRMS) hydrologic model [27,28] was selected to simulate natural flows. The Hydrologic Engineering Center River Analysis System (HEC-RAS) hydraulic model [29] was selected to simulate the effects of channel change over time. The parallel results of flow and channel change indicated the cause-and-effect relation between hydrologic modification and channel form. This relation was assessed to identify the desired, stable hydrologically modified flow regime and channel form. The individual components of the approach are detailed below and are summarized in Figure 2. To achieve this, these in situ data are necessary: (1) measured modified flows at the study watershed, (2) measured natural flows at the reference watershed, and (3) baseline and repeat stream surveys of the sensitive reach at the study watershed ( Figure 1).

of 18
Hydrology 2020, 7, x FOR PEER REVIEW 5 of 17 Figure 2. Schematic outlining the method used to assess an existing, unstable, hydrologically modified system. Flow and channel alteration were assessed using a combination of hydrologic and hydraulic data with models to identify a desired, future, stable, hydrologically modified system.
Measured natural flows and stream surveys of the sensitive reach at the study watershed do not exist prior to hydrologic modification, which occurred in the 1960s, following the construction of the water collection system and again in the 1980s following enlargement of that system. The current system became fully operational by 1987 (Table A1). Expanding upon the unique constraints of this assessment, two explicit assumptions were made: (1) natural flows simulated for the study watershed for the period 1995-2015 will also be considered representative for the previous years, and (2) in the absence of geometric surveys for a period before the hydrologic modification, remnant features will be defined and used. The reliability of these remnant features as representative of pre-modification discharges shall also be discussed in light of the modification date.
The specific data selected begins with availability of the modified flow record at the study watershed, which were measured by the Wyoming State Engineer and Cheyenne Board of Public Utilities [30]. Daily average flow for the period 1995 to 2015 were used. The beginning of the water year (WY) was set on October 1st, following the standard used in the U.S. Measured natural flows at the reference watershed were measured as part of the USGS Hydrologic Benchmark Network, due to its long-term, unregulated flow record [31]. Daily average flow for the same time period as the study watershed measured modified flow dataset were used. Sensitive reach stream survey data from 2006 were repeated in 2015 after a long, continuous period of hydrologic modification. The baseline and repeat survey data included cross-sections (XSs), reach lengths between XSs, longitudinal Schematic outlining the method used to assess an existing, unstable, hydrologically modified system. Flow and channel alteration were assessed using a combination of hydrologic and hydraulic data with models to identify a desired, future, stable, hydrologically modified system. Measured natural flows and stream surveys of the sensitive reach at the study watershed do not exist prior to hydrologic modification, which occurred in the 1960s, following the construction of the water collection system and again in the 1980s following enlargement of that system. The current system became fully operational by 1987 (Table A1). Expanding upon the unique constraints of this assessment, two explicit assumptions were made: (1) natural flows simulated for the study watershed for the period 1995-2015 will also be considered representative for the previous years, and (2) in the absence of geometric surveys for a period before the hydrologic modification, remnant features will be defined and used. The reliability of these remnant features as representative of pre-modification discharges shall also be discussed in light of the modification date.
The specific data selected begins with availability of the modified flow record at the study watershed, which were measured by the Wyoming State Engineer and Cheyenne Board of Public Utilities [30]. Daily average flow for the period 1995 to 2015 were used. The beginning of the water year (WY) was set on October 1st, following the standard used in the U.S. Measured natural flows at the reference watershed were measured as part of the USGS Hydrologic Benchmark Network, due to its long-term, unregulated flow record [31]. Daily average flow for the same time period as the study watershed measured modified flow dataset were used. Sensitive reach stream survey data from 2006 were repeated in 2015 after a long, continuous period of hydrologic modification. The baseline and repeat survey data included cross-sections (XSs), reach lengths between XSs, longitudinal profiles, and pebble counts. Approximately a dozen XSs were spaced out to capture two meander wavelengths of the sensitive reach. A longitudinal profile identified the hydraulic grade controls and slopes throughout the sensitive reach.
Two applications of the hydrologic model were required to simulate the reference watershed and then the study watershed ( Figure 2). PRMS is a deterministic, distributed-parameter, physically based model that uses distributed parameters over sub-watersheds called hydrologic response units (HRUs) [27,28]. For each HRU, PRMS balanced energy and water budgets of the snowpack, plant canopy, and soil zone to simulate hydrologic processes including snowmelt, sublimation, interception, infiltration, evapotranspiration (ET), interflow, groundwater (GW) flow, and surface runoff [28].
The key input variables for each application of PRMS included temperature (T), precipitation (P), solar radiation (SR), and ET. The T and P data were obtained from NRCS SNOTEL stations [25] surrounding the study area ( Figure 1). The T and P data were distributed using three-dimensional multiple-linear regression based on latitude, longitude, and elevation [32]. The SR data were estimated by a degree-day relation [33]. These SR data estimates were then adjusted using daylight hours, days with P, and potential SR [34]. The ET data were estimated using the modified Jensen-Haise method and then distributed based on plant cover, soil properties, and potential ET (PET) [33,34]. In addition to distributing climate input variables, unique geospatial parameter values were extracted from stream network, land cover, canopy density, soils, elevation, slope, aspect, and radiation plane data with the GIS Weasel program [35].
The parameterization process for the reference watershed was automated by the Let Us Calibrate (LUCA) multiple-objective, stepwise calibration method. This method calibrated model variables (e.g., SR; ET; and monthly, daily, low, and high flows) against measured datasets over a ten-year calibration period (2006-2015) by adjusting the model parameter values [32]. These measured datasets used for calibration included global horizontal irradiance [36], actual ET [37], and daily average flow [31]. The parameter values used in each step were selected based on the sensitivity analysis conducted by the USGS on the adjacent Yampa River basin for snowmelt dominated basins [32]. LUCA calculated one parameter value for the watershed by averaging across all HRUs as opposed to for each HRU. LUCA calculated the accuracy of the modeled variables compared to the measured datasets using two objective functions, sum of the absolute difference in logarithms and normalized root mean square error [32]. The final parameter values were set after six rounds of six steps (Table 2). Each step continued until the objective functions made no improvement in accuracy following a series of parameter value adjustment iterations. With the final parameters (Table A2), the model was run for a separate ten-year evaluation period (1996)(1997)(1998)(1999)(2000)(2001)(2002)(2003)(2004)(2005) to test the parameterization. Nash-Sutcliffe Efficiency (NSE) [38] and Percent Bias (PBIAS) were used as metrics for determining satisfactory overall model performance [38,39].
The last step of the parameterization process adjusted all the parameters that could be recalculated using values derived from the unique characteristics from the study watershed. These adjusted parameters included those relating to site-specific land cover, elevation, aspect, stream network, and climate characteristics. There were 21 sensitive parameters obtained from a reference watershed parameterization that could not be recalculated, such as coefficients for routing groundwater and snowmelt ( Table 2). Utilizing both the adjusted parameter values from the study watershed parameterization and the 21 unadjusted parameters obtained from reference watershed parameterization, the study watershed natural flow dataset was simulated. The HEC-RAS hydraulic model [39] was parameterized for using both the baseline and repeat stream surveys of the sensitive reach ( Figure 1). The inputs included slope, XS geometry, reach lengths between XSs, and Manning's n estimates. The average bed slope, which approximated the energy slope, was used as the downstream boundary condition for the model. The active channel Manning's n value was initially estimated using flow depth and grain size from pebble count data [40]. Additional roughness from bank irregularities along the reach was added to the initial active channel Manning's n value [41]. The floodplain Manning's n value estimates were obtained from the table values of winter and summer, medium-to-dense brush [42].
A range of representative high flows were selected for the hydraulic models. The range of representative high flows were the 1-year to 20-year recurrence intervals calculated from the study watershed simulated natural daily average flows. For each flow between the 1-and 20-year recurrence intervals, the hydraulic model computed a corresponding water surface elevation as well as other hydraulic data such as wetted perimeter and flood inundation area.
Using the hydrologic and hydraulic data with models allowed the concurrent assessment of the flow and channel alteration in response to hydrologic modification (Figure 2). Flow alteration was indicated by comparing the attributes of simulated natural flows to measured modified flows using the Indicators of Hydrologic Alteration (IHA) tool [43]. IHA calculates a deviation factor to express the magnitude of change between pre-and post-impact periods. Deviation factor is calculated by taking the absolute value of the difference between the post-impact and pre-impact values and then dividing by the pre-impact value. To measure the significance of the deviation, IHA calculates a significance count for the deviation values by randomly shuffling all years of input data and recalculating fictitious pre-and post-impact medians 1000 times. The significance count is the number of fictitious deviation values that are greater than for the actual deviation factor out of 1000. If this fraction is close to zero, it means that the deviation is highly significant [43]. As a way to characterize the connection between flow and channel alteration, bankfull flows were analyzed. Specifically, the flood recurrence interval estimates of bankfull flow from pre-and post-hydrologic modification were compared to geomorphic estimates of bankfull capacity. This method is described in more detail below.
The key attributes of flow alteration were the monthly magnitudes, magnitude and duration of annual extremes, the timing of annual extremes, frequency and duration of pulse flows, rate and frequency changes, recurrence intervals, and snowmelt period flood characteristics. The snowmelt period was identified as the subset of the annual hydrograph from the onset of the rising limb to the terminus of the falling limb.
Since the channel was in a state of disequilibrium, the visual estimates of bankfull elevation at each XS were inconsistent within the sensitive reach ( Figure 1). While surveying the sensitive reach, reliable indicators of bankfull stage were difficult to find due to the amount of channel change. Instead, the geomorphic bankfull for each XS was estimated using wetted perimeter-flow curves [22]. A curve for each XS, for both 2006 and 2015, was plotted using the simulated wetted perimeter against flow. The geomorphic bankfull was indicated on the curve by a break in slope, called a major breakpoint. Conceptually, a major breakpoint represents the threshold where the water within the active channel reaches the top of the stream bank and spills into the floodplain, thereby substantially increasing the wetted perimeter for negligible increases in flow. A subtler break in slope, called a minor breakpoint, likely indicated a remnant geomorphic bankfull. While one could mathematically approximate the breakpoints by curve fitting, we visually determined them, which yielded a possible error of 0.1 to 0.2 m 3 /s.
Drawn from the evidence that channel form responds to flow, the existing channel may be brought into a quasi-equilibrium near its current geomorphic capacity by defining a corresponding hydrologically modified flow regime that is stable over a long-term period. The attributes for this new, stable, modified flow regime incorporated the flow attributes of the existing geomorphic channel capacity and the simulated natural flow regime where possible. The attributes that were defined included the frequency, magnitude, duration, rise rate, fall rate, and timing for low, intermediate, and peak flows. A range of variability was applied to the attributes of the new, stable, modified flow regime to accommodate dry, average, and wet years. Namely, the attributes of the new, stable, modified flow regime bounded the hydrologically modified flow regime to promote stabilization of the existing channel over a long-term period of time. In doing so, this provides a means to both meet water resource objectives and protect environmental sustainability.

Results
The reference watershed hydrologic model parameterization had a ten-year evaluation period NSE of 0.82% and PBIAS of −4.5%, as compared to the ten-year calibration period NSE of 0.89 and PBIAS of −1.9%. The simulated streamflow matched the observed streams at the reference well for average, high and low flow years (Figure 3). The individual years each had NSE above 60% and PBIAS within 20%. The lower NSE values in individual years indicated that snowmelt runoff was simulated early or late. The higher PBIAS values in individual years indicated the simulated total annual runoff volume was overestimated or underestimated. With satisfactory model performance [39], the final parameters values were selected for use in the study watershed hydrologic model application.
Hydrology 2020, 7, x FOR PEER REVIEW 9 of 17 The parameterization of the study watershed hydrologic model simulated daily average flows from 1995-2015. A comparison with the measured modified flows indicated hydrologic modifications, specifically the later timing of peak flows, lower peak flows, and lower low flows, especially in the winter (Figure 4a and Table 3); further cumulative volumes were less (Figure 4b).
The most altered attributes were low flows and monthly flows in the winter and spring (Table 3). To a slightly lesser degree, other attributes were lowered, including flood frequency, magnitude, duration, timing (later), and fall rates (Table 3)  The parameterization of the study watershed hydrologic model simulated daily average flows from 1995-2015. A comparison with the measured modified flows indicated hydrologic modifications, specifically the later timing of peak flows, lower peak flows, and lower low flows, especially in the winter (Figure 4a and Table 3); further cumulative volumes were less (Figure 4b). The most altered attributes were low flows and monthly flows in the winter and spring (Table 3). To a slightly lesser degree, other attributes were lowered, including flood frequency, magnitude, duration, timing (later), and fall rates (Table 3) The most altered attributes were low flows and monthly flows in the winter and spring (Table 3). To a slightly lesser degree, other attributes were lowered, including flood frequency, magnitude, duration, timing (later), and fall rates (Table 3). Channel change assessed by tracking geomorphic bankfull capacity was indicated by the minor and major breakpoints on the wetted perimeter-flow curves ( Figure 5 Table 3. Comparison of the simulated natural and measured modified flow attributes to assess the degree of alteration using a deviation factor and the significance of the alteration, ranging from highly significant (0) to not significant (1).

Flow Attribute
Meas. Modified Flow Sim. Natural Flow Deviation Significance  Figure 3. Table 3. Comparison of the simulated natural and measured modified flow attributes to assess the degree of alteration using a deviation factor and the significance of the alteration, ranging from highly significant (0) to not significant (1).   The range of bankfull flow estimates, from 2.8 to 7.6 m 3 /s, is large ( Figure 6). This variability indicated the potential for a future, unstable, hydrologically modified flow regime. Since the 1.5-year flood flow suggests the tendency of the geomorphic bankfull capacity, it was selected as the main criterion for assessing future channel stability. Assuming dam regulation of the upstream structures could be made to accommodate a desired hydrology, a magnitude of 5.7 m 3 /s was selected as the desired 1.5-year flood flow that would limit future channel change. Namely, this magnitude is between the 2015 geomorphic bankfull capacity indicated by the wetted perimeter-flow curve's major breakpoint of 5.5 m 3 /s ( Figure 5) and the 1987-2015 measured modified flow regime 1.5-year flood flow of 5.8 m 3 /s ( Figure 6). The 2-year, 5-year, and 20-year floods were selected by increasing the 1.5-year flood flow magnitude at similar percent increase intervals to the simulated natural flow flood recurrence interval percent increases (Table 4).  The range of bankfull flow estimates, from 2.8 to 7.6 m 3 /s, is large ( Figure 6). This variability indicated the potential for a future, unstable, hydrologically modified flow regime. Since the 1.5-year flood flow suggests the tendency of the geomorphic bankfull capacity, it was selected as the main criterion for assessing future channel stability. Assuming dam regulation of the upstream structures could be made to accommodate a desired hydrology, a magnitude of 5.7 m 3 /s was selected as the desired 1.5-year flood flow that would limit future channel change. Namely, this magnitude is between the 2015 geomorphic bankfull capacity indicated by the wetted perimeter-flow curve's major breakpoint of 5.5 m 3 /s ( Figure 5) and the 1987-2015 measured modified flow regime 1.5-year flood flow of 5.8 m 3 /s ( Figure 6). The 2-year, 5-year, and 20-year floods were selected by increasing the 1.5-year flood flow magnitude at similar percent increase intervals to the simulated natural flow flood recurrence interval percent increases (Table 4). The new, stable, modified flow regime incorporated two major components, the flow attributes of the existing geomorphic bankfull capacity of the channel (Table 4) and the individual attributes of the simulated natural flow regime ( Table 3). The first major component, the recurrence intervals (RIs), were grouped into peak flow size classes. The peak flow size classes included no overbank flows (no floods) with 1-to 1.5-year RI, small overbank flows (small floods) with 1.5-to 5-year RIs, and large overbank flows (large floods) with 5-to 20-year RIs (Figure 7). The magnitudes of each peak flow were based on the existing geomorphic bankfull capacity (Table 4). One indicator to predict the size class of peak flow is the peak Snow Water Equivalent (SWE) of the snowpack measured at nearby SNOTEL stations ( Figure 1). The duration, rise rates, fall rates, and timing were based on the simulated natural flow regime ( Table 3). The simulated natural flow regime indicated that the rising limb would begin no earlier than mid-April; the peak would occur, on average, near early June but no earlier than mid-May and no later than early July, and the falling limb would end no later than mid-August ( Figure 7).
Hydrology 2020, 7, x FOR PEER REVIEW 12 of 17 no earlier than mid-May and no later than early July, and the falling limb would end no later than mid-August ( Figure 7). In addition to peak flow attributes, the intermediate and low flows were estimated based on the simulated natural flow regime ( Table 3). The intermediate flows were restricted to the late spring and the late summer during the rising and falling limb periods (Figure 7). Moreover, these flows were simulated only as transition flows, meaning they were not sustained for long durations, but were instead either falling or rising. The magnitudes of the intermediate flows were between low flow and In addition to peak flow attributes, the intermediate and low flows were estimated based on the simulated natural flow regime ( Table 3). The intermediate flows were restricted to the late spring and the late summer during the rising and falling limb periods (Figure 7). Moreover, these flows were simulated only as transition flows, meaning they were not sustained for long durations, but were instead either falling or rising. The magnitudes of the intermediate flows were between low flow and the 1.1-year RI peak flow. The low flow magnitude was approximated to be 0.15 m 3 /s, which was slightly greater than the 7-day low flow of 0.13 m 3 /s (Table 3). Together the peak, intermediate, and low flow attributes provided a foundation for a new, stable, modified flow regime to promote future channel stability (Figure 7).

Discussion
To assess hydrologic response in a hydrologically modified watershed, PUB methods were utilized. This PUB method involved two applications of the PRMS hydrologic model, for reference and study watersheds. Since measured natural-flow data were unavailable at the study watershed, hydrologic model calibration and evaluation could not be performed in the typical way for the study watershed application of the PRMS hydrologic model. Instead, typical calibration and evaluation was conducted only for reference watershed hydrologic model parameterization. This allowed an overall model performance for the reference watershed application of the PRMS hydrologic model, which was an NSE of 0.89 for calibration and NSE of 0.82 for evaluation. For comparison, twelve applications of the PRMS hydrologic model in the Willamette River basin near the Cascade Range of Oregon, resulted in an NSE ranging between 0.62-0.93 for calibration and 0.58-0.90 for evaluation [19].
Furthermore, the Willamette River study treated five of the 12 applications of the PRMS hydrologic model as ungauged basins to test how well the regionalization method simulated flows in gauged basins, and found acceptable model performance, with an NSE of 0.60-0.76 [19]. Without any further calibration or evaluation, the study watershed application of the PRMS hydrologic model was parameterized using the regionalization method with the unique study watershed parameters and the missing sensitive parameters obtained from calibration of the reference watershed model. Although each PUB method includes an unknown degree of uncertainty, the main uncertainties in this approach were assumed to be minimized by using an adjacent watershed with similar climates, topography, geology, and land cover.
Recognizing the uncertainty, this particular PUB regionalization technique demonstrated its value in producing a long-term dataset of simulated daily average flow. An accurate estimation of the natural flow record is necessary in determining key components of the natural flow regime [17]. The long-term simulated natural daily average flow record at the study watershed did successfully identify the degree to which each attribute of the natural flow regime was altered (Table 3). While hydrologic modification did illustrate changes, such as flows that were lower with peaks occurring later (Figure 3), part of this difference could have come from the parameterization process that is based on the transfer of parameters calibrated on the reference watershed. As seen in this study, the specific type of flow alteration is dependent on the unique environmental setting and impoundment, diversion, and water collection system operation. This study also found many likenesses to similar studies. For example, the study watershed low flows increased three-fold similar to the 10-fold increase in low flows at the River Ter, UK [9]. The intermediate flows increased similarly to the Milk River, Alberta, as indicated by monthly average flow increases of 2-3 times [8]. The peak flows increased by nearly doubling, similar to the La Poudre Pass Creek of the Poudre River and Lake Creek of the Arkansas River in Colorado [10,11].
With the change in hydrology, in this study, as well as the many other similar cases, there was an alteration in channel form [7,[9][10][11]. This response is consistent with the qualitative response model and findings of similar trans-basin diversion case studies [5]. In quantifying the hydraulic response of the study watershed hydrologic modification, a key validation in the PUB regionalization technique used in this study was perhaps found, which may provide one more tool for understanding the uncertainty.
The agreement between the 1.5-year peak flow of the simulated natural flow regime and the flow associated with the minor breakpoint found consistently in the sensitive reach cross-sections ( Figure 5) appears to provide some validation for this PUB regionalization technique providing an adequate simulation of the natural flow regime.
In turn, characterizations of the simulated natural flow regime and the measured modified flow regime provided validation for the capability of the hydraulic simulation method. As noted above, there was agreement between simulated natural flow attributes (i.e., 1.5-year flood flow) and remnant minor breakpoints on the wetted-perimeter-flow curves ( Figure 6). There was agreement between measured modified flow attributes (i.e., 1-5-year flood flow) and major breakpoints in the sense that the channel bankfull capacity increase tracks the increasing 1.5-year flood flow ( Figure 6). With these consistencies, the hydraulic response method of collecting baseline and repeat geometry data and parameterization of the HEC-RAS model was deemed adequate.
Lastly, as indicated by the linkage between hydrology and hydraulics, a stable flow regime was found to be an integral component when considering future channel stabilization. Since a channel adjusts to both its water and sediment supplies, the other key factor for considering future channel stabilization is the sediment regime. In this case, it was assumed that the increased flows from augmentation outweighed the sediment deficit from the reservoir storage, dam, and outlet works. Thus, this approach focuses on the sediment transport relation to a stable flow regime, while also acknowledging the need for a better understanding of the changing sediment regime in hydrologically modified systems.

Conclusions
The linking of hydrologic and hydraulic data with models provided the foundation needed to understand the changes induced by the hydrologic modification of a transboundary watershed where collection systems export water. The Prediction in Ungauged Basin framework was used to simulate mean daily flows to evaluate the flow conditions of the natural flow regime and the corresponding flow alterations. Peak flows are less and occur later, with lower base flows and lower ramping (rise and fall) rates due to flow modification. Investigating a sensitive reach within a hydrologically modified system using hydraulic modeling demonstrated a change in the response of the system that would have occurred under natural flow conditions. The relation of hydrology and channel form is integral for improving environmental sustainability at complex, hydrologically modified systems. The design of a new, stable, modified flow regime demonstrated the potential for promoting a stable channel form. Lower flows are required at various recurrence intervals to maintain a new stable channel environment.

Appendix B
The appendix presents the optimized parameter values estimated for the reference watershed.