Bankfull Regional Curves for the Inner and Outer Bluegrass Regions of Kentucky1

Brockman, Ruth R., Carmen T. Agouridis, Stephen R. Workman, Lindell E. Ormsbee, and Alex W. Fogle, 2012. Bankfull Regional Curves for the Inner and Outer Bluegrass Regions of Kentucky. Journal of the American Water Resources Association (JAWRA) 48(2): 391‐406. DOI: 10.1111/j.1752‐1688.2011.00621.x Abstract: Bankfull regional curves that relate channel dimensions and discharge to watershed drainage area are useful tools for assisting in the correct identification of bankfull elevation and in stream restoration and reconstruction. This study assessed 28 stable streams located in two physiographic regions of Kentucky: the Inner Bluegrass and the Outer Bluegrass. Bankfull channel dimensions, discharge, and return period as well as average channel slope, median bed material size, sinuosity, Rosgen stream classification, and percent impervious area were determined. Significant relationships were found between drainage area and the bankfull characteristics of cross‐sectional area, width, mean depth, and discharge for both the Inner Bluegrass and Outer Bluegrass regions (α = 0.05). It was also found that the percent impervious area in a watershed had minimal effect on bankfull dimensions, which is attributed to the well‐vegetated nature of the streambanks, cohesive streambank materials, and bedrock control. No significant differences between any of the Inner Bluegrass and Outer Bluegrass regional curves were found (α = 0.05). Comparisons were made between the Inner Bluegrass and Outer Bluegrass curves and others developed in karst‐influenced areas in the Eastern United States. Although few significant differences were found between the regional curves for bankfull discharge and width, a number of the curves differed with regards to bankfull cross‐sectional area and mean depth.     

Bankfull Regional Curves for North and Northwest Florida Streams1

Abstract: Regional curves, which relate bankfull channel dimensions and discharge to watershed drainage area, are developed to aid in identifying the bankfull stage in ungaged watersheds, and estimating the bankfull discharge and dimensions for river studies and natural channel design applications. This study assessed 26 stable stream reaches in two hydro‐physiographic regions of the Florida Coastal Plain: the Northwest Florida Coastal Plain (NWFCP) and the North Florida Coastal Plain (NFCP). Data from stream reaches in Georgia and Alabama were also used to develop the Florida regional curves, since they are located in the same hydro‐physiographic region. Reaches were selected based on the presence of U.S. Geological Survey gage stations and indicators of limited watershed development (e.g., <10% impervious surface). Analyses were conducted to determine bankfull channel dimensions, bankfull discharge, average channel slope, and Rosgen stream classification. Based on these data, significant relationships were found between bankfull cross‐sectional area, width, mean depth, and discharge as a function of drainage area for both regions. Data from this study suggested that bankfull discharges and channel dimensions were larger from NWFCP streams than from Coastal Plain streams in North Carolina and Maryland. Bankfull discharges were similar between NFCP and Georgia coastal plain streams; therefore, the data were combined into one regional curve. In addition, the data were stratified by Rosgen stream type. This stratification strengthened the relationships of bankfull width and mean depth as a function of drainage area.     

FLUVIAL GEOMORPHOLOGICAL METHODOLOGY FOR NATURAL STABLE CHANNEL DESIGN1

A fluvial geomorphological methodology for designing natural stable channels is being widely applied for river restoration. It is an analogue procedure, as the W/d ratio and sinuosity from a reference reach are scaled to determine the restoration design. The choice of reference reach is crucial and published criteria specify that it should be stable, correspond to the stream type at the restoration site, have the same valley type, and be from the same hydrophysiographic region. For stable, meandering gravel cobble bed rivers flowing through alluvial flood plains (C3 and C4 stream types), UK regime equations are used to evaluate the procedure. Successful design requires particular combinations of the ratios of bankfull discharge, bed material size and load, valley slope, and bank vegetation category between the reference and restoration sites. These critical ratios, which are confirmed by U.S. field data, provide guidelines for selecting a suitable reference reach for C3‐C4 stream types. They also indicate that the reference reach can be in any valley type or hydrophysiographic region. The geomorphological procedure will apply to all stable stream types, provided the reference reach is correctly identified. Specific guidelines for each stream type await the development of additional regime equations.     

Bankfull Regional Curves for the Alleghany Plateau/Valley and Ridge,Piedmont, and Coastal Plain Regions of Maryland

Regional curves are empirical relationships that can help identify the bankfull stage in ungaged watersheds and aid in designing the riffle dimension in stream restoration projects. Bankfull regional curves were developed from gage stations with drainage areas less than 102 mi² (264.2 km²) for the Alleghany Plateau/Valley and Ridge (AP/VR), Piedmont, and Coastal Plain regions of Maryland. The AP/VR regions were combined into one region for this project. These curves relate bankfull discharge, cross‐sectional area, width, and mean depth to drainage area within the same hydro‐physiographic region (region with similar rainfall/runoff relationship). The bankfull discharge curve for the Coastal Plain region was further subdivided into the Western Coastal Plain (WCP) and Eastern Coastal Plain (ECP) region due to differences in topography and runoff. Results show that the Maryland Piedmont yields the highest bankfull discharge rate per unit drainage area, followed by the AP/VR, WCP, and ECP. Likewise, the Coastal Plain and AP/VR streams have less bankfull cross‐sectional area per unit drainage area than the Piedmont. The average bankfull discharge return interval across the three hydro‐physiographic regions was 1.4 years. The Maryland regional curves were compared to other curves in the eastern United States. The average bankfull discharge return interval for the other studies ranged from 1.1 to 1.8 years.     

Development of Regional Curves for Hydrologic Landscape Regions (HLR) in the Contiguous United States

Regional curves relate drainage area to the bankfull channel characteristics discharge, cross‐sectional area, width, and mean depth. These curves are used for a variety of purposes, including aiding in the field identification of bankfull elevation and in the natural channel design process. When developing regional curves, the degree to which landform, geology, climate, and vegetation influence stream systems within a single physiographic province may not be fully considered. This study examined the use of the U.S. Geological Survey's Hydrologic Landscape Regions (HLR), as well as data from 2,856 independent sites throughout the contiguous United States (U.S.), to develop a set of regional curves (bankfull discharge, cross‐sectional area, width, and mean depth) for (1) the contiguous U.S., (2) each of the 20 HLRs, (3) each of the eight physiographic divisions, (4) 22 of the 25 physiographic provinces, and (5) individual HLRs within the physiographic provinces. These regional curves were then compared to each other, as well as those from the literature. Regional curves developed for individual HLRs, physiographic divisions, and physiographic provinces tended to outperform the contiguous U.S. indicating increased stratification was beneficial. Further stratifying physiographic provinces by HLR markedly improved regional curve reliability. Use of HLR as a basis of regional curve development, rather than physiographic region alone, may allow for the development of more robust regional curves.     

Estimating the Frequency and Extent of Overbank Flow of a Baseflow‐Dominated Stream Flowing through a Wetland

River floodplains provide critical habitat for a wide range of animal and plant species and reduce phosphorus and nitrogen loads in streams. It has been observed that baseflow‐dominated streams flowing through wetlands are commonly at or near bankfull and overflow their banks much more frequently than other streams. However, there is very little published quantitative support for this observation. The study focuses on a 1‐km reach of Black Earth Creek, a stream in the Midwestern United States (U.S.). We used one‐dimensional hydraulic modeling to estimate bankfull discharge at evenly spaced stream cross sections, and two‐dimensional modeling to quantitate the extent of wetland inundation as a function of discharge. We then used historical streamflow data from two U.S. Geological Survey gaging stations to quantitate the frequency of wetland inundation. For the with‐sediment case, the frequency of overbank conditions at the 38 cross sections in the wetland ranged from 3 to 85 days per year and averaged 43 days per year. Ten percent of the wetland was inundated for an average of 35 days per year. For the without‐sediment case, the frequency of overbank conditions ranged from 2.6 to 48 days per year and averaged 14 days per year. Also, 10% of the wetland was inundated for an average of 25 days per year. These unusually high rates of floodplain inundation are likely due in part to the very low stream gradient and shallow depths of overbank flow.     

HYDRAULIC GEOMETRY RELATIONSHIPS FOR URBAN STREAMS THROUGHOUT THE PIEDMONT OF NORTH CAROLINA1

ABSTRACT: Hydraulic geometry relationships, or regional curves, relate bankfull stream channel dimensions to watershed drainage area. Hydraulic geometry relationships for streams throughout North Carolina vary with hydrology, soils, and extent of development within a watershed. An urban curve that is the focus of this study shows the bankfull features of streams in urban and suburban watersheds throughout the North Carolina Piedmont. Seventeen streams were surveyed in watersheds that had greater than 10 percent impervious cover. The watersheds had been developed long enough for the streams to redevelop bankfull features, and they had no major impoundments. The drainage areas for the streams ranged from 0.4 to 110.3 square kilometers. Cross‐sectional and longitudinal surveys were conducted to determine the channel dimension, pattern, and profile of each stream and power functions were fitted to the data. Comparisons were made with regional curves developed previously for the rural Piedmont, and enlargement ratios were produced. These enlargement ratios indicated a substantial increase in the hydraulic geometry for the urban streams in comparison to the rural streams. A comparison of flood frequency indicates a slight decrease in the bankfull discharge return interval for the gaged urban streams as compared to the gaged rural streams. The study data were collected by North Carolina State University (NCSU), the University of North Carolina at Charlotte (UNC), and Charlotte Storm Water Services. Urban regional curves are useful tools for applying natural channel design in developed watersheds. They do not, however, replace the need for field calibration and verification of bankfull stream channel dimensions.     

Improved Bankfull Discharge Prediction Using 2‐Year Recurrence‐Period Discharge1

Abstract:  Knowledge of bankfull discharge (Q_bf) is essential for planners, engineers, geomorphologists, environmentalists, agricultural interests, developments situated on flood prone lands, surface mining and reclamation activities, and others interested in floods and flooding. In conjunction with estimating Q_bf, regionalized bankfull hydraulic geometry relationships, which relate Q_bf and associated channel dimensions (i.e., width, depth, and cross‐section area) to drainage basin area (A_da), are often used. This study seeks to improve upon the common practice of predicting Q_bf using A_da exclusively. Specifically, we hypothesize that predictions of Q_bf can be improved by including estimates of the 2‐year recurrence‐period discharge (Q₂) in regression models for predicting Q_bf. For testing this hypothesis, we used Q_bf estimates from 30 reports containing data for streams that span 34 hydrologic regions in 16 states. Corresponding values of Q₂ and A_da were compiled from flood‐frequency reports and other sources. By comparing statistical measures (i.e., root mean squared error, coefficient of determination, and Akaike’s information criterion), we determined that predicting Q_bf from Q₂ rather than A_da yields consistently better estimates of Q_bf. Other principal findings are (1) data are needed for at least 12 sites in a region for reliable hydraulic geometry model selection and (2) an approximate range of values for Q_bf/Q₂ is 0.10‐3.0.     

UNCERTAINTY IN ESTIMATING BANKFULL CONDITIONS1

ABSTRACT: Bankfull depth and discharge are basic input parameters to stream planform, stream restoration, and highway crossing designs, as well as to the development of hydraulic geometry relationships and the classification of streams. Unfortunately, there are a wide variety of definitions for bankfull that provide a range of values, and the actual selection of bankfull is subjective. In this paper, the relative uncertainty in determining the bankfull depth and discharge is quantified, first by examining the variability in the estimates of bankfull and second by using fuzzy numbers to describe bankfull depth. Fuzzy numbers are used to incorporate uncertainty due to vagueness in the definition of bankfull and subjectivity in the selection of bankfull. Examples are provided that demonstrate the use of a fuzzy bankfull depth in sediment trans. port and in stream classification. Using fuzzy numbers to describe bankfull depth rather than a deterministic value allows the engineer to base designs and decisions on a range of possible values and associated degrees of belief that the bankfull depths take on each value in that range.     

PUMP SYSTEM FOR IMPROVING RUNOFF RECORDS DURING WINTER1

ABSTRACT: Freezing winter temperatures can cause icing of outdoor weirs used to measure surface runoff. Ice typically forms in the notch and on the crests of short-crested V-notch weirs, causing incorrect gage heights to be measured. A method for reducing the effects of ice formation on weirs and weir pools using a pump is presented and evaluated. Warmer water from the bottom of the weir pool is pumped to the surface, reducing the opportunity for the water surface to freeze. The pump is shown to work except under extremely cold conditions, improving runoff records from 27 percent to 60 percent. The pump system has no practical effect on measured gage height. Frequency distributions of flow rates and air temperatures under measured ice-free and other weir conditions are presented. Suggestions for use of the pump system under temperature conditions other than those in this study are given.     

Development and Evaluation of Bankfull Hydraulic Geometry Relationships for the Physiographic Regions of the United States

Bankfull hydraulic geometry relationships are used to estimate channel dimensions for streamflow simulation models, which require channel geometry data as input parameters. Often, one nationwide curve is used across the entire United States (U.S.) (e.g., in Soil and Water Assessment Tool), even though studies have shown that the use of regional curves can improve the reliability of predictions considerably. In this study, regional regression equations predicting bankfull width, depth, and cross‐sectional area as a function of drainage area are developed for the Physiographic Divisions and Provinces of the U.S. and compared to a nationwide equation. Results show that the regional curves at division level are more reliable than the nationwide curve. Reliability of the curves depends largely on the number of observations per region and how well the sample represents the population. Regional regression equations at province level yield even better results than the division‐level models, but because of small sample sizes, the development of meaningful regression models is not possible in some provinces. Results also show that drainage area is a less reliable predictor of bankfull channel dimensions than bankfull discharge. It is likely that the regional curves can be improved using multiple regression models to incorporate additional explanatory variables.     

Development and Comparison of Multiple Regression Models to Predict Bankfull Channel Dimensions for Use in Hydrologic Models

Channel dimensions are important input variables for many hydrologic models. As measurements of channel geometry are not available in most watersheds, they are often predicted using bankfull hydraulic geometry relationships. This study aims at improving existing equations that relate bankfull width, depth, and cross‐sectional area to drainage area (DA) without limiting their use to well‐gauged watersheds. We included seven additional variables in the equations that can be derived from data that are generally required by hydrologic models anyway and conducted several multiple regression analyses to identify the ideal combination of additional variables for nationwide and regional models for each Physiographic Division of the United States (U.S.). Results indicate that including the additional variables in the regression equations generally improves predictions considerably. The selection of relevant variables varies by Physiographic Division, but average annual precipitation (PCP) and temperature (TMP) were generally found to improve the models the most. Therefore, we recommend using regression equations with three independent variables (DA, PCP, and TMP) to predict bankfull channel dimensions for hydrologic models. Furthermore, we recommend using the regional equations for watersheds within regions from which data were used for model development, whereas in all other parts of the U.S. and the rest of the world, the nationwide equations should be given preference.     

A GIS‐ASSISTED DISTRIBUTED WATERSHED MODEL FOR SIMULATING FLOODING AND INUNDATION1

ABSTRACT: A distributed watershed model combining kinematic wave routing, 1‐D dynamic channel‐flow routing, and 2‐D diffusive overland‐flow routing has been developed to simulate flooding and inundation levels of large watersheds. The study watershed was linked to a GIS database and was divided into an upstream mountainous area and a downstream alluvial plain. A kinematic wave routing was adopted at the mountainous area to compute the discharge flowing into the alluvial plain. A 1‐D dynamic channel routing solving the St. Venant equations by the Preissmann method was performed for the main channel of the alluvial plain, whereas a 2‐D overland‐flow routing solving the diffusion wave equation with the Alternating Direction Explicit scheme was used for floodplains. The above two routings were connected by weir‐link discharge formula. The parameters in the model were calibrated and independently verified by single‐event storms. An example application of flooding/inundation analysis was conducted for the Taichung station and the Woozi depot (Taiwan High Speed Rail). Suggested inundation‐proofing measures ‐ including raising ground surface elevation of the station and depot and building a waterproofing exterior wall and their combination ‐ were investigated. It was concluded that building the waterproofing exterior wall had a strong tendency to decrease peak inundation depth.     

TRANSIENT STORAGE ASSESSMENTS OF DYE‐TRACER INJECTIONS IN RIVERS OF THE WILLAMETTE BASIN,OREGON1

ABSTRACT: Rhodamine WT dye‐tracer injections in rivers of the Willamette Basin yield concentration‐time curves with characteristically long recession times suggestive of active transient storage processes. The scale of drainage areas contributing to the stream reaches studied in the Willamette Basin ranges from 10 to 12,000 km². A transient storage assessment of the tracer studies has been completed using the U.S. Geological Survey's One‐dimensional Transport with Inflow and Storage (OTIS) model, which incorporates storage exchange and decay functions along with the traditional dispersion and advection transport equation. The analysis estimates solute transport of the dye. It identifies first‐order decay coefficients to be on the order of 10^?5/sec for the nonconservative Rhodamine W.T. On an individual subreach basis, the first‐order decay is slower (typically by an order of magnitude) than the transient storage process, indicating that nonconservative tracers may be used to evaluate transient storage in rivers. In the transient storage analysis, a dimensionless parameter (A_s/A) expresses the spatial extent of storage zone area relative to stream cross section. In certain reaches of Willamette Basin pool‐and‐riffle, gravel‐bed rivers, this parameter was as large as 0.5. A measure of the storage exchange flux was calculated for each stream subreach in the simulation analysis. This storage exchange is shown subjectively to be higher at higher stream discharges. Hyporheic linkage between streams and subsurface flows is the probable physical mechanism contributing to a significant part of this inferred active transient storage. Hyporheic linkages are further suggested by detailed measurements of river discharge with an Acoustic Doppler Current Profiler system delineating zones in two large rivers where water alternately enters and leaves the surface channels through gravel‐and‐cobble riverbeds. Measurements show patterns of hyporheic exchange that are highly variable in time and space.     

Evaluation of Overtopping Riprap Design Relationships

Rock riprap is one of the most widely used erosion control methods for protecting embankments, levees, spillways, and instream structures subjected to overtopping flow conditions. At least 21 stone‐sizing relationships exist to determine the median stone size of a protective riprap layer based on the results of 96 overtopping, laboratory experiments. Test parameters include median stone size, slope, unit discharge, coefficient of uniformity, and riprap layer thickness. A regression analysis was performed relating the observed median stone size to the predicted median stone size to each of the 21 relationships, yielding a coefficient of determination (R²) and percent error for the full spectrum of data. Zonal (partial spectrum of rock sizes) and complexity analyses were also conducted for each relationship. It was resolved that the Khan and Ahmad, and Chang relationships best aligned with the composite dataset. The predictive expressions by Olivier, Hartung and Scheuerlein, Knauss, Maynord, Abt and Johnson, and Siebel yield a noteworthy second tier of stone‐sizing relationships for overtopping conditions.     

AVERAGE DISCHARGE,PERENNIAL FLOW INITIATION,AND CHANNEL INITIATION ‐ SMALL SOUTHERN APPALACHIAN BASINS1

ABSTRACT: Regional average evapotranspiration estimates developed by water balance techniques are frequently used to estimate average discharge in ungaged streams. However, the lower stream size range for the validity of these techniques has not been explored. Flow records were collected and evaluated for 16 small streams in the Southern Appalachians to test whether the relationship between average discharge and drainage area in streams draining less than 200 acres was consistent with that of larger basins in the size range (> 10 square miles) typically gaged by the U.S. Geological Survey (USGS). This study was designed to evaluate predictors of average discharge in small ungaged streams for regulatory purposes, since many stream regulations, as well as recommendations for best management practices, are based on measures of stream size, including average discharge. The average discharge/drainage area relationship determined from gages on large streams held true down to the perennial flow initiation point. For the southern Appalachians, basin size corresponding to perennial flow is approximately 19 acres, ranging from 11 to 32 acres. There was a strong linear relationship (R²= 0.85) between average discharge and drainage area for all streams draining between 16 and 200 acres, and the average discharge for these streams was consistent with that predicted by the USGS Unit Area Runoff Map for Georgia. Drainage area was deemed an accurate predictor of average discharge, even in very small streams. Channel morphological features, such as active channel width, cross‐sectional area, and bankfull flow predicted from Manning's equation, were not accurate predictors of average discharge. Monthly baseflow statistics also were poor predictors of average discharge.     

ASSESSING STREAM CHANNEL STABILITY THRESHOLDS USING FLOW COMPETENCE ESTIMATES AT BANKFULL STAGE1

ABSTRACT: Stream channel stability is affected by peak flows rather than average annual water yield. Timber harvesting and other land management activities that contribute to soil compaction, vegetation removal, or increased drainage density can increase peak discharges and decrease the recurrence interval of bankfull discharges. Increased peak discharges can cause more frequent movement of large streambed materials, leading to more rapid stream channel change or instability. This study proposes a relationship between increased discharge and channel stability, and presents a methodology that can be used to evaluate stream channel stability thresholds on a stream reach basis. Detailed surveys of the channel cross section, water surface slope, streambed particle size distribution, and field identification of bankfull stage are used to estimate existing bankfull flow conditions. These site specific stream channel characteristics are used in bed load movement formulae to predict critical flow conditions for entrainment of coarse bed material (D₈₄ size fraction). The “relative bed stability” index, defined as the ratio of critical flow condition to the existing condition at bankfull discharge, can predict whether increased peak discharges will exceed stream channel thresholds.     

Comparison of Three Pebble Count Protocols (EMAP,PIBO, and SFT) in Two Mountain Gravel‐Bed Streams1

Abstract: Although the term ``pebble count'' is in widespread use, there is no standardized methodology used for the field application of this procedure. Each pebble count analysis is the product of several methodological choices, any of which are capable of influencing the final result. Because there are virtually countless variations on pebble count protocols, the question of how their results differ when applied to the same study reach is becoming increasingly important. This study compared three pebble count protocols: the reach‐averaged Environmental Monitoring and Assessment Program (EMAP) protocol named after the EMAP developed by the Environmental Protection Agency, the habitat‐unit specific U.S. Forest Service’s PACFISH/INFISH Biological Opinion (PIBO) Effectiveness Monitoring Program protocol, and a data‐intensive method developed by the authors named Sampling Frame and Template (SFT). When applied to the same study reaches, particle‐size distributions varied among the three pebble count protocols because of differences in sample locations within a stream reach and along a transect, in particle selection, and particle‐size determination. The EMAP protocol yielded considerably finer, and the PIBO protocol considerably coarser distributions than the SFT protocol in the pool‐riffle study streams, suggesting that the data cannot be used interchangeably. Approximately half of the difference was due to sampling at different areas within the study reach (i.e., wetted width, riffles, and bankfull width) and at different locations within a transect. The other half was attributed to using different methods for particle selection from the bed, particle‐size determination, and the use of wide, nonstandard size classes. Most of the differences in sampling outcomes could be eliminated by using simple field tools, by collecting a larger sample size, and by systematically sampling the entire bankfull channel and all geomorphic units within the reach.     

BANKFULL DISCHARGE RECURRENCE INTERVALS AND REGIONAL HYDRAULIC GEOMETRY RELATIONSHIPS: PATTERNS IN THE PACIFIC NORTHWEST,USA1

ABSTRACT: The model bankfull discharge recurrence interval (annual series) (T_a) in streams has been approximated at a 1.5‐year flow event. This study tests the linkage between regional factors (climate, physiography, and ecoregion) and the frequency of bank‐full discharge events in the Pacific Northwest (PNW). Patterns of Ta were found to be significant when stratified by EPA Ecoregion. The mean value for T_a in the PNW is 1.4 years; however, when the data is stratified by ecoregion, the humid areas of western Oregon and Washington have a mean value of 1.2 years, while the dryer areas of Idaho and eastern Oregon and Washington have a mean value of 1.4 to 1.5 years. Among the four factors evaluated, vegetation association and average annual precipitation are the primary factors related to channel form and T_a. Based on the results of the T_a analyses, regional hydraulic geometry relationships of streams were developed for the PNW, which relate variables, such as bank‐full cross‐sectional area, width, depth, and velocity, to bankfull discharge and drainage area. The verification of T_a values, combined with the development of regional hydraulic geometry relationships, provides geographically relevant information that will result in more accurate estimates of hydraulic geometry variables in the PNW.     

Stationarity of Streamflow Records and Their Influence on Bankfull Regional Curves1

Haucke, Jessica and Katherine A. Clancy, 2011. Stationarity of Streamflow Records and Their Influence on Bankfull Regional Curves. Journal of the American Water Resources Association (JAWRA) 47(6):1338–1347. DOI: 10.1111/j.1752‐1688.2011.00590.x Abstract: Bankfull regional curves, which are curves that establish relationships among channel morphology, discharge, drainage area, are used extensively for stream restoration. These curves are developed upon the assumption that streamflows maintain stationarity over the entire record. We examined this assumption in the Driftless Area of southwestern Wisconsin where agricultural soil retention practices have changed, and precipitation has increased since the 1970s. We developed a bankfull regional curve for this area using field surveys of bankfull channel performed during 2008‐2009 and annual series of peak streamflows for 10 rivers with streamflow records ranging from the 1930s to 2009. We found bankfull flows to correlate to a 1.1 return period. To evaluate gage data statistics, we used the sign test to compare our channel morphology to historic 1.5 return period discharge (Q_1.5) for five time periods: 1959‐1972, 1973‐1992, 1993‐2008, 1999‐2008, and the 1959‐2008 period of record. Analysis of the historic gage data indicated that there has been a more than 30% decline in Q_1.5 since 1959. Our research suggests that land conservation practices may have a larger impact on gaging station stationarity than annual precipitation changes do. Additionally, historic peak flow data from gages, which have records that span land conservation changes, may need to be truncated to represent current flow regimes.