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.     

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.     

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.     

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.     

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 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.     

A Framework to Develop Nationwide Flooding Extents Using Climate Models and Assess Forecast Potential for Flood Resilience

The methods used to simulate flood inundation extents can be significantly improved by high‐resolution spatial data captured over a large area. This paper presents a hydraulic analysis methodology and framework to estimate national‐level floodplain changes likely to be generated by climate change. The hydraulic analysis was performed using existing published Federal Emergency Management Agency 100‐year floodplains and estimated 100‐ and 10‐year return period peak flow discharges. The discharges were estimated using climate variables from global climate models for two future growth scenarios: Representative Concentration Pathways 2.6 and 8.5. River channel dimensions were developed based on existing regional United States Geological Survey publications relating bankfull discharges with channel characteristics. Mathematic relationships for channel bankfull topwidth, depth, and side slope to contributing drainage area measured at model cross sections were developed. The proposed framework can be utilized at a national level to identify critical areas for flood risk assessment. Existing hydraulic models at these “hot spots” could be repurposed for near–real‐time flood forecasting operations. Revitalizing these models for use in simulating flood scenarios in near–real time through the use of meteorological forecasts could provide useful information for first responders of flood emergencies.     

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.     

ESTIMATING STREAMFLOW AND ASSOCIATED HYDRAULIC GEOMETRY,THE MID‐ATLANTIC REGION,USA1

ABSTRACT: Methods to estimate streamflow and channel hydraulic geometry were developed for unpaged streams in the Mid‐Atlantic Region. Observed mean annual streamflow and associated hydraulic geometry data from 75 gaging stations in the Appalachian Plateau, the Ridge and Valley, and the Piedmont Physiographic Provinces of the Mid‐Atlantic Region were used to develop a set of power functions that relate streamflow to drainage area and hydraulic geometry to streamflow. For all three physiographic provinces, drainage area explained 95 to 98 percent of the variance in mean annual streamflow. Relationships between mean annual streamflow and water surface width and mean flow depth had coefficients of determination that ranged from R²= 0.55 to R²= 0.91, but the coefficient of determination between mean flow velocity and mean annual streamflow was lower (R²= 0.44 to R²= 0.54). The advantages of using the regional regression models to estimate streamflow over a conceptual model or a water balance model are its ease of application and reduced input data needs. The prediction of the regression equations were tested with data collected as part of the U.S. Environmental Protection Agency (USEPA) Environmental Monitoring and Assessment Program (EMAP). In addition, equations to transfer streamflow from gaged to ungaged streams are presented.     

Optimizing Bankfull Discharge and Hydraulic Geometry Relations for Streams in New York State1

Mulvihill, Christiane I. and Barry P. Baldigo, 2012. Optimizing Bankfull Discharge and Hydraulic Geometry Relations for Streams in New York State. Journal of the American Water Resources Association (JAWRA) 48(3): 449-463. DOI: 10.1111/j.1752-1688.2011.00623.x Abstract: This study analyzes how various data stratification schemes can be used to optimize the accuracy and utility of regional hydraulic geometry (HG) models of bankfull discharge, width, depth, and cross-sectional area for streams in New York. Topographic surveys and discharge records from 281 cross sections at 82 gaging stations with drainage areas of 0.52-396 square miles were used to create log-log regressions of region-based relations between bankfull HG metrics and drainage area. The success with which regional models distinguished unique bankfull discharge and HG patterns was assessed by comparing each regional model to those for all other regions and a pooled statewide model. Gages were also stratified (grouped) by mean annual runoff (MAR), Rosgen stream type, and water-surface slope to test if these models were better predictors of HG to drainage area relations. Bankfull discharge models for Regions 4 and 7 were outside the 95% confidence interval bands of the statewide model, and bankfull width, depth, and cross-sectional area models for Region 3 differed significantly (p < 0.05) from those of other regions. This study found that statewide relations between drainage area and HG were strongest when data were stratified by hydrologic region, but that co-variable models could yield more accurate HG estimates in some local regional curve applications.     

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.     

River Channel Geometry and Rating Curve Estimation Using Height above the Nearest Drainage

River channel geometry is an important input to hydraulic and hydrologic models. Traditional approaches to quantify river geometry have involved surveyed river cross sections, which cannot be extended to ungaged basins. In this paper, we describe a method for developing a synthetic rating curve to relate flow to water level in a stream reach based on reach‐averaged channel geometry properties developed using the Height above Nearest Drainage (HAND) method. HAND uses a digital elevation model (DEM) of the terrain and computes the elevation difference between each land surface cell and the stream bed cell to which it drains. Taking increments in water level in the stream, HAND defines the inundation zone and a water depth grid within this zone, and the channel characteristics are defined from this water depth grid. We apply our method to the Blanco River (Texas) and the Tar River (North Carolina) using 10‐m terrain data from the United States Geological Survey (USGS) 3D Elevation Program (3DEP) dataset. We evaluate the method's performance by comparing the reach‐average stage‐river geometry relationships and rating curves to those from calibrated Hydrologic Engineering Center's River Analysis System (HEC‐RAS) models and USGS gage observations. The results demonstrate that after some adjustment, the river geometry information and rating curves derived from HAND using national‐coverage datasets are comparable to those obtained from hydraulic models or gage measurements. We evaluate the inundation extent and show our approach is able to capture the majority of the Federal Emergency Management Agency (FEMA) 100‐year floodplain.     

Evaluating Pre‐ and Post‐Fire Peak Discharge Predictions across Western U.S. Watersheds

This study reviews five models commonly used in post‐fire hydrologic assessments: the Rowe Countryman and Storey (RCS), United States Geological Survey (USGS) Linear Regression Equations, USDA Windows Technical Release 55 (USDA TR‐55), Wildcat5, and U.S. Army Corps of Engineers (USACE) Hydrologic Modeling System (HEC‐HMS). The models are applied to eight diverse basins in the western United States (U.S.) (Arizona, California, Colorado, Montana, and Washington) affected by wildfires and assessed by input parameters, calibration methods, model constraints, and performance. No one model is versatile enough for application to all study sites. Results show inconsistency between model predictions for events across the sites and less confidence with larger return periods (25‐ and 50‐year events) and post‐fire predictions. The RCS method performs well, but application is limited to southern California. The USGS linear regression model has wider regional application, but performance is less reliable at the large recurrence intervals and post‐fire predictions are reliant on a subjective modifier. Of the three curve number‐based models, Wildcat5 performs best overall without calibration, whereas the calibrated TR‐55 and HEC‐HMS models show significant improvement in pre‐fire predictions. Results from our study provide information and guidance to ultimately improve model selection for post‐fire prediction and encourage uniform parameter acquisition and calibration across the western U.S.     

Stage‐Discharge Relationships for U‐, A‐, and W‐Weirs in Un‐submerged Flow Conditions1

Thornton, Christopher I., Anthony M. Meneghetti, Kent Collins, Steven R. Abt, and S. Michael Scurlock, 2011. Stage‐Discharge Relationships for U‐, A‐, and W‐Weirs in Un‐submerged Flow Conditions. Journal of the American Water Resources Association (JAWRA) 47(1):169‐178. DOI: 10.1111/j.1752‐1688.2010.00501.x Abstract: Instream rock weirs are routinely placed into stream systems to provide grade control, reduce streambank erosion, provide energy dissipation, and allow fish passage. However, design and performance criteria for site specific applications are often anecdotal or qualitative in nature, and based upon the experience of the design team. A study was conducted to develop generic state‐discharge relationships for U‐, A‐, and W‐weirs. A laboratory testing program was performed in which scaled, near‐prototype U‐, A‐, and W‐rock weir structures were constructed in 11 configurations. Each configuration encompassed a unique weir shape, bed material, and/or bed slope. Thirty‐one tests were conducted in which each structure was subjected to a sequence of predetermined discharges that minimally included the equivalent of 1/3 bankfull, 2/3 bankfull, and bankfull conditions. All tests were performed in subcritical, un‐submerged flow conditions. Stage‐discharge relationships were developed using multivariant, power regression techniques for each of the U‐, A‐, and W‐rock weirs as a function of the effective weir length, flow depth, mean weir height, rock size, and discharge coefficient. Unique coefficient expressions were developed for each weir shape, and a single discharge coefficient was proposed applicable to the weirs for determining the channel stage‐discharge rating.     

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.     

Channel Roughness in North Carolina Mountain Streams

Channel roughness, often described by Manning's n, is used to represent the amount of resistance that flow encounters, and has direct implications on velocity and discharge. Ideally, n is calculated from a long‐term record of channel discharge and hydraulic geometry. In the absence of these data, a combination of photo references and a validated qualitative method is preferable to simply choosing n arbitrarily or from a table. The purpose of this study was to use United States Geological Survey (USGS) streamflow data to calculate roughness coefficients for streams in the mountains of North Carolina. Five USGS gage stations were selected for this study, representing drainage areas between 71.5 and 337 km². Photo references of the study sites are presented. Measured discharges were combined with hydraulic geometry at a cross‐section to calculate roughness coefficients for flows of interest. At bankfull flow, n ranged between 0.039 and 0.064 for the five study sites. Roughness coefficients were not constant for all flows in a channel, and fluctuated over a large range. At all sites, roughness was highest during low‐flow conditions, then quickly decreased as flow increased, up to the bankfull elevation.     

Forecasting the Probability of Future Groundwater Levels Declining Below Specified Low Thresholds in the Conterminous U.S.

We present a logistic regression approach for forecasting the probability of future groundwater levels declining or maintaining below specific groundwater‐level thresholds. We tested our approach on 102 groundwater wells in different climatic regions and aquifers of the United States that are part of the U.S. Geological Survey Groundwater Climate Response Network. We evaluated the importance of current groundwater levels, precipitation, streamflow, seasonal variability, Palmer Drought Severity Index, and atmosphere/ocean indices for developing the logistic regression equations. Several diagnostics of model fit were used to evaluate the regression equations, including testing of autocorrelation of residuals, goodness‐of‐fit metrics, and bootstrap validation testing. The probabilistic predictions were most successful at wells with high persistence (low month‐to‐month variability) in their groundwater records and at wells where the groundwater level remained below the defined low threshold for sustained periods (generally three months or longer). The model fit was weakest at wells with strong seasonal variability in levels and with shorter duration low‐threshold events. We identified challenges in deriving probabilistic‐forecasting models and possible approaches for addressing those challenges.     

Improvements in Fluvial Stability Associated with Two‐Stage Ditch Construction in Mower County,Minnesota

Water quality and stream habitat in agricultural watersheds are under greater scrutiny as hydrologic pathways are altered to increase crop production. Agricultural drainage ditches function to remove water quickly from farmed landscapes. Conventional ditch designs lack the form and function of natural stream systems and tend to be unstable and provide inadequate habitat. In October of 2009, 1.89 km of a conventional drainage ditch in Mower County, Minnesota, was converted to an alternative system with a two‐stage channel to investigate the improvements in water quality, stability, and habitat. Longitudinal surveys show a 12‐fold increase in the pool‐riffle formation. Cross‐sectional surveys show an average increase in bankfull width of approximately 10% and may be associated to an increased frequency in large storm events. The average increase in bankfull depth was estimated as 18% but is largely influenced by pool formation. Rosgen Stability Analyses show the channel to be highly stable and the banks at a low risk of erosion. The average bankfull recurrence interval was estimated to be approximately 0.30 years. Overall, the two‐stage ditch design demonstrates an increase in fluvial stability, creating a more consistent sediment budget, and increasing the frequency of important instream habitat features, making this best management practice a viable option for addressing issues of erosion, sediment imbalance, and poor habitat in agricultural drainage systems.     

A Power Law Model for River Flow Velocity in Iowa Basins

This study explores power law relationships to estimate water flow velocity as a function of discharge and drainage area across river networks. We test the model using empirical data from 214 United States (U.S.) Geological Survey gauging stations distributed over the state of Iowa in the U.S. The empirical data are the measurements of the mean cross‐sectional velocity and concurrent discharge. The data are used to estimate parameters for a state‐wide model and to test for spatial variability for 15 large river basins contained within the state. Spatial differences among the basins are small but some parameters significantly differ from the state‐wide model. Using individual station data, the authors also explore a simpler power law model that disregards dependence on the drainage area. Overall, the study shows that including drainage area improves the model. Our study provides parameter values that can be directly incorporated into a regional scale routing model, and provides a framework for developing flow velocity models for hydraulically similar rivers in the U.S. and the world.     

COMPARISON OF MODELS FOR ESTIMATING FLOOD QUANTILES IN NEW HAMPSHIRE AND VERMONT1

We develop and compare three regression models for estimating flood quantiles at ungaged stream reaches in New Hampshire and Vermont. These models emerge from systematic analysis and validation of relations between flood magnitude and six candidate predictors reflecting basin size, topography, and climate and channel size at 36 gaging stations with record lengths exceeding 20 years. Of the candidate predictors, bank full width is most highly correlated with flood magnitude and the best prediction equation is based on width. Thus channel geometry is closely related to the current hydrologic regime in spite of geologically recent glaciation and apparently non-alluvial bank materials. We also develop models that use information obtainable from maps or GIS. The best of these uses drainage area and drainage-basin elevation as predictors, but it is substantially less precise than the width-based relation. A third relation using only drainage area as a predictor is even less precise but may be useful for some purposes. No other single predictors or combinations yielded useful predictions, although some had been included in previously-established models for the region. Model comparison included examination of residuals generated by regression using one-at-a-time suppression of data points and comparison with precision obtainable with gaging records of varying lengths.