FTABLE Generation Method Effects on Instream Fecal Bacteria Concentrations Simulated With HSPF1

Abstract: Computer simulation models are used extensively for the development of total maximum daily loads (TMDLs). Specifically, the Hydrological Simulation Program‐FORTRAN (HSPF) is used in Virginia for the development of TMDLs for bacteria impairments. HSPF estimates discharge from a reach using function tables (FTABLES). The FTABLE relates stream stage, surface area, and volume to discharge from a reach. In this study, five FTABLE estimation methods were assessed by comparing their effect on various simulation outputs. Four “field‐based” methods used detailed cross‐sectional data collected via site surveys. A fifth “digital‐based” method used digital elevation data in combination with the Natural Resources Conservation Service Regional Hydraulic Geometry Curves. Sets of FTABLEs created using each method were used in simulations of instream bacteria concentration for a Virginia watershed. Several statistics relating to instream bacteria including long‐term average concentration, die‐off, and the violation rate of Virginia’s bacteria criterion were compared. The pair‐wise Student’s t‐test was used for the comparison. The HSPF simulations that used FTABLES estimated from digitally based data consistently produced significantly higher long‐term average instream fecal bacteria concentrations, significantly lower instream fecal bacteria die‐off, which is related to differences in residence time in the streams, and significantly higher water quality criterion violation rates.     

Enhancing Hyrdological Simulation Program – FORTRAN Model Channel Hydraulic Representation1

Abstract: The hydrological simulation program – FORTRAN (HSPF) is a comprehensive watershed model that employs depth‐area‐volume‐flow relationships known as the hydraulic function table (FTABLE) to represent the hydraulic characteristics of stream channel cross‐sections and reservoirs. An accurate FTABLE determination for a stream cross‐section site requires an accurate determination of mean flow depth, mean flow width, roughness coefficient, longitudinal bed slope, and length of stream reach. A method that uses regional regression equations to estimate mean flow depth, mean flow width, and roughness coefficient is presented herein. FTABLES generated by the proposed method (Alternative Method) and FTABLES generated by Better Assessment Science Integrating Point and Nonpoint Sources (BASINS) were compared. As a result, the Alternative Method was judged to be an enhancement over the BASINS method. First, the Alternative Method employs a spatially variable roughness coefficient, whereas BASINS employs an arbitrarily selected spatially uniform roughness coefficient. Second, the Alternative Method uses mean flow width and mean flow depth estimated from regional regression equations whereas BASINS uses mean flow width and depth extracted from the National Hydrography Dataset (NHD). Third, the Alternative Method offers an option to use separate roughness coefficients for the in‐channel and floodplain sections of compound channels. Fourth, the Alternative Method has higher resolution in the sense that area, volume, and flow data are calculated at smaller depth intervals than the BASINS method. To test whether the Alternative Method enhances channel hydraulic representation over the BASINS method, comparisons of observed and simulated streamflow, flow velocity, and suspended sediment were made for four test watersheds. These comparisons revealed that the method used to estimate the FTABLE has little influence on hydrologic calibration, but greatly influences hydraulic and suspended sediment calibration. The hydrologic calibration results showed that observed versus simulated daily streamflow comparisons had Nash‐Sutcliffe efficiencies ranging from 0.50 to 0.61 and monthly comparisons had efficiencies ranging from 0.61 to 0.84. Comparisons of observed and simulated suspended sediments concentrations had model efficiencies ranging from 0.48 to 0.56 for the daily, and 0.28 to 0.70 for the monthly comparisons. The overall results of the hydrological, hydraulic, and suspended sediment concentration comparisons show that the Alternative Method yielded a relatively more accurate FTABLE than the BASINS method. This study concludes that hydraulic calibration enhances suspended sediment simulation performance, but even greater improvement in suspended sediment calibration can be achieved when hydrological simulation performance is improved. Any improvements in hydrological simulation performance are subject to improvements in the temporal and spatial representation of the precipitation data.     

Simulation of Surface Water for Un‐Gauged Areas with Storage‐Attenuation Wetlands1

Abstract: A nontraditional application of the Hydrological Simulation Program – FORTRAN (HSPF) model to simulate freshwater discharge to upper Charlotte Harbor along Florida’s west coast was performed. This application was different from traditional HSPF applications in three ways. First, the domain of the model was defined based on the hydraulic characteristics of the landforms using small distributed parameter discretization. Second, broad wetland land forms, representing more than 20% of this area, were simulated as reaches with storage‐attenuation characteristics and not as pervious land segments (PERLNDs). Finally, the reach flow‐tables (F‐Tables) were configured in a unique way to be calibrated representing the uncertainty of the storage‐attenuation process. Characterizing wetlands as hydrography elements allows flow from the wetlands to be treated as a stage‐dependent flux. The study was conducted for the un‐gauged portion of the Peace and Myakka rivers in west‐central Florida. Due to low gradient tidal influences, a large portion of the basin is un‐gauged. The objective of this study was to simulate stream flow discharges and to estimate freshwater inflow from these un‐gauged areas to upper Charlotte Harbor. Two local gauging stations were located within the model domain and were used for calibration. Another gauge with a shorter period of record was used for verification. A set of global hydrologic parameters were selected and tested using the parameter optimization software (PEST) during the calibration. Model results were evaluated using PEST and well‐known statistical indices. The correlation coefficients were very high (0.899 and 0.825) for the two calibration stations. Further testing of this approach appears warranted for watersheds with significant wetlands coverage.     

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.     

CONSTRUCTING A PROBLEM SOLVING ENVIRONMENT TOOL FOR HYDROLOGIC ASSESSMENT OF LAND USE CHANGE1

ABSTRACT: This paper describes the integration of a comprehensive hydrological model known as the Hydrological Simulation Program Fortran (HSPF) into a problem solving environment (PSE) for watershed management. The original PSE concept was a structure providing web‐based access to a suite of models, including HSPF and other models of in‐stream hydrodynamics, biological impacts and economic effects, for the watershed‐wide assessment of alternative land use scenarios. The present paper describes only the HSPF integration into the PSE program. Example applications to the 148 square kilometer (57 square mile) Back Creek subwatershed in the upper Roanoke River system (1,479 square kilometers or 571 square miles) in southwest Virginia are used to illustrate important concepts and linkages between land development and hydrological change using hypothetical' what if'scenarios. The features of HSPF and its limitations in this context are discussed. The paper as such is a proof‐of‐concept paper and not a completion report. It is intended to describe the PSE tool building process rather than analysis of the many possible simulation outcomes. However, the dominance of raw imperviousness as a contributor to hydrograph response is apparent in all the PSE simulations described in this paper.     

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.     

Estimating Evapotranspiration for Dryland Cropping Systems in the Semiarid Texas High Plains Using SWAT

The Soil and Water Assessment Tool (SWAT) is one of the most widely used watershed models for simulating hydrology in response to agricultural management practices. However, limited studies have been performed to evaluate the SWAT model's ability to estimate daily and monthly evapotranspiration (ET) in semiarid regions. ET values were simulated using ArcSWAT 2012 for a lysimeter field managed under dryland conditions at the USDA‐ARS Conservation and Production Research Laboratory at Bushland, Texas, and compared with measured lysimeter values from 2000 to 2010. Two scenarios were performed to compare SWAT's performance: (1) use of default plant leaf area index (LAI) values in the embedded plant database and (2) adjusted LAI values. Scenario 1 resulted in an “unsatisfactory” Nash‐Sutcliffe efficiency (NSE) of 0.42 and 0.38 for the calibration and validation periods, respectively. Scenario 2 resulted in a “satisfactory” NSE value for the calibration period while achieving a “good” NSE of 0.70 for the validation period. SWAT generally underestimated ET at both the daily and monthly levels. Overestimation during fallow years may be due to the limitations of the pothole function used to simulate furrow diking. Users should be aware of potential errors associated with using default LAI parameters. Inaccuracies in ET estimation may also stem from errors in the plant stress functions, particularly when evaluating water management practices for dryland watersheds.     

Improvement in HSPF's Low‐Flow Predictions by Implementation of a Power Law Groundwater Storage‐Discharge Relationship

We have enhanced the ability of a widely used watershed model, Hydrologic Simulation Program — FORTRAN (HSPF), to predict low flows by reconfiguring the algorithm that simulates groundwater discharge. During dry weather periods, flow in most streams consists primarily of base flow, that is, groundwater discharged from underlying aquifers. In this study, HSPF's groundwater storage‐discharge relationship is changed from a linear to a more general nonlinear relationship which takes the form of a power law. The nonlinear algorithm is capable of simulating streamflow recession curves that have been found in some studies to better match observed dry weather hydrographs. The altered version of HSPF is implemented in the Chesapeake Bay Program's Phase 5 Model, an HSPF‐based model that simulates nutrient and sediment loads to the Chesapeake Bay, and is tested in the upper Potomac River basin, a 29,950 km² drainage area that is part of the Bay watershed. The nonlinear relationship improved median Nash‐Sutcliffe efficiencies for log daily flows at the model's 45 calibration points. Mean absolute percent error on low‐flow days dropped in five major Potomac River tributaries by up to 12 percentage points, and in the Potomac River itself by four percentage points, where low‐flow days were defined as days when observed flows were in the lowest 5th percentile range. Percent bias on low‐flow days improved by eight percentage points in the Potomac River, from ?11 to ?3%.     

Daily Bank Erosion Rates in the Lower Yellow River Before and After Dam Construction

During the period of water impoundment and sediment detention of the Sanmenxia Reservoir, riverbank erosion processes played a key role in the channel evolution of the Lower Yellow River (LYR). However, research into bank erosion rates of the LYR has been neglected due to the lack of direct field monitoring. In this study, an indirect method is proposed to determine bank erosion rates at daily time scales by outlining a detailed calculation procedure using measured hydrological data. A total of 810 data points of daily bank erosion rates before and after the construction of Sanmenxia Dam was calculated at seven hydrometric sections along the LYR, with the corresponding values of the bank stability coefficient and the width‐to‐depth ratio also being calculated. Empirical relations were then developed to estimate the daily bank erosion rates, using these parameters at the sections. Temporal and spatial variability in daily bank erosion rates in the LYR before and after dam construction were also investigated, revealing that: (1) the bank erosion rates had a mean value of 16.7‐29.1 m/day in the braided reach, with a maximum value of 290.0 m/day, while they were relatively low in the meandering reach, with a mean value of 2.5 m/day; (2) the erosion rates before dam construction were slightly greater than those after dam construction, with the difference reaching 5‐10 m/day in the braided reach, decreasing in the transitional reach gradually, and being slight in the meandering reach.     

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.     

HYDROLOGICAL MODELING OF THE IROQUOIS RIVER WATERSHED USING HSPF AND SWAT1

ABSTRACT: The performance of two popular watershed scale simulation models — HSPF and SWAT — were evaluated for simulating the hydrology of the 5,568 km² Iroquois River watershed in Illinois and Indiana. This large, tile drained agricultural watershed provides distinctly different conditions for model comparison in contrast to previous studies. Both models were calibrated for a nine‐year period (1987 through 1995) and verified using an independent 15‐year period (1972 through 1986) by comparing simulated and observed daily, monthly, and annual streamflow. The characteristics of simulated flows from both models are mostly similar to each other and to observed flows, particularly for the calibration results. SWAT predicts flows slightly better than HSPF for the verification period, with the primary advantage being better simulation of low flows. A noticeable difference in the models' hydrologic simulation relates to the estimation of potential evapotranspiration (PET). Comparatively low PET values provided as input to HSPF from the BASINS 3.0 database may be a factor in HSPF's overestimation of low flows. Another factor affecting baseflow simulation is the presence of tile drains in the watershed. HSPF parameters can be adjusted to indirectly account for the faster subsurface flow associated with tile drains, but there is no specific tile drainage component in HSPF as there is in SWAT. Continued comparative studies such as this, under a variety of hydrologic conditions and watershed scales, provide needed guidance to potential users in model selection and application.     

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.     

Field Testing the Riparian Ecosystem Management Model on a Riparian Buffer in the North Carolina Upper Coastal Plain

The riparian ecosystem management model (REMM) was field tested using five years (2005‐2009) of measured hydrologic and water quality data on a riparian buffer located in the Tar‐Pamlico River Basin, North Carolina. The buffer site received NO₃‐N loading from an agricultural field that was fertilized with inorganic fertilizer. Field results showed the buffer reduced groundwater NO₃‐N concentration moving to the stream over a five‐year period. REMM was calibrated hydrologically using daily field‐measured water table depths (WTDs), and with monthly NO₃‐N concentrations in groundwater wells. Results showed simulated WTDs and NO₃‐N concentrations in good agreement with measured values. The mean absolute error and Willmott's index of agreement for WTDs varied from 13‐45 cm and 0.72‐0.92, respectively, while the root mean square error and Willmott's index of agreement for NO₃‐N concentrations ranged from 1.04‐5.92 mg/l and 0.1‐0.86, respectively, over the five‐year period. REMM predicted plant nitrogen (N) uptake and denitrification were within ranges reported in other riparian buffer field studies. The calibrated and validated REMM was used to simulate 33 years of buffer performance at the site. Results showed that on average the buffer reduced NO₃‐N concentrations from 12 mg/l at the field edge to 0.7 mg/l at the stream edge over the simulation period, while the total N and NO₃‐N load reductions from the field edge to the stream were 77 and 82%, respectively.     

MODELING NITROGEN TRANSPORT IN THE IPSWICH RIVER BASIN,MASSACHUSETTS, USING A HYDROLOGICAL SIMULATION PROGRAM IN FORTRAN (HSPF)1

ABSTRACT: Increased riverine nitrogen (N) fluxes have been strongly correlated with land use changes and are now one of the largest pollution problems in the coastal region of the United States. In the present study, the Hydrological Simulation Program‐FORTRAN (HSPF) is used to simulate transport of N in the Ipswich River basin in Massachusetts and to evaluate the effect of future land use scenarios on the water quality of the river. Model results show that under a land use change scenario constructed with restrictions from environmental protection laws, where 44 percent of the forest in the basin was converted to urban land, stream nitrate concentrations increased by about 30 percent of the present values. When an extreme land use scenario was used, and 100 percent of the forest was converted to urban land, concentrations doubled in comparison to present values. Model simulations also showed that present stream nitrate concentrations might be four times greater than they were prior to urbanization. While pervious lands with high density residential land use generated runoff with the highest N concentrations in HSPF simulations, the results suggested that denitrification in the riparian zone and wetlands coupled with the hydrology of the basin are likely to control the magnitude of nitrate loads to the aquatic system. The simulation results showed that HSPF can predict the general patterns of inorganic N concentrations in the Ipswich River and tributaries. Nevertheless, HSPF has some difficulty simulating the extreme variability of the observed data throughout the main stem and tributaries, probably because of limitations in the representation of wetlands and riparian zones in the model, where N processes such as denitrification seem to play a major role in controlling the transport of N from the terrestrial system to the river reaches.     

Impact of Length of Dataset on Streamflow Calibration Parameters and Performance of APEX Model

Due to resource constraints, long‐term monitoring data for calibration and validation of hydrologic and water quality models are rare. As a result, most models are calibrated and, if possible, validated using limited measured data. However, little research has been done to determine the impact of length of available calibration data on model parameterization and performance. The main objective of this study was to evaluate the impact of length of calibration data (LCD) on parameterization and performance of the Agricultural Policy Environmental eXtender model for predicting daily, monthly, and annual streamflow. Long‐term (1984‐2015) measured daily streamflow data from Rock Creek watershed, an agricultural watershed in northern Ohio, were used for this study. Data were divided into five Short (5‐year), two Medium (15‐year), and one Long (25‐year) streamflow calibration data scenarios. All LCD scenarios were calibrated and validated at three time steps: daily, monthly, and annual. Results showed LCD affected the ability of the model to accurately capture temporal variability in simulated streamflow. However, overall average streamflow, water budgets, and crop yields were simulated reasonably well for all LCD scenarios.     

ASSESSMENT OF RESTORED RWERINE HABITAT USING RCHARC1

ABSTRACT: The Riverine Community Habitat Assessment and Restoration Concept (RCHARC) was developed to integrate habitat enhancement into the stream restoration process. RCHARC assumes that aquatic habitat quality is closely related to hydraulic diversity based upon a “comparison standard” reach approach to stream restoration. A Beta test was performed by applying the RCHARC process to Rapid Creek in Rapid City, South Dakota. Standard and restored stream reaches were selected and data were collected. A comparison of field data and velocity-depth distributions indicated that the restored stream closely replicated the standard reach. The RCHARC methodology has the potential to assess habitat quality for planned comparison reaches and indicate the level of success resulting from restoration.     

Characterizing Storm Hydrograph Rise and Fall Dynamics With Stream Stage Data1

Abstract: Storm‐flow transients (i.e., hydrograph rise and fall dynamics) may represent an important aspect of understanding streamflow dynamics. However, little is known about how temporal resolution of transient data and climate variability may color these potential indicators of hydrologic pattern or condition. Warm‐season stream stage and rainfall were monitored continuously (5 min) during the 2002 water year in eight tributaries of the Little Miami River (Ohio), which drain 17‐58 km² catchments. Rise rates generated using 5‐min data were different than those generated with mean daily data [calculated with the Indicators of Hydrologic Alteration (IHA) software], though fall rates were similar for fine and coarse temporal data. This result suggests that data with low temporal resolution may not be adequate to fully represent the dynamics of storm rise rates. Conversely, fall rates based on daily stage data (via IHA) were similar to those based on the 5‐min data, and so daily mean data may be appropriate for characterizing fall rates. We next analyzed the possible correlations between rainfall variability and storm‐flow stage dynamics. We derived rise and recession rates from storm stage hydrographs by assuming exponential rise and decay of a runoff peak. We found that raw rise rates (R_raw) were correlated with both the maximum rainfall rate and the time to the centroid of a rain event. We subsequently removed the trend based on these rainfall characteristics, which yielded new representations of rise rates abbreviated as R_rate and R_tcent, respectively, and that had lower variability than the uncorrected (raw) data. Fall rates were found to be independent of rainfall characteristics. Due to the predominant influence of stream hydrology upon aquatic biota and nutrient fluxes, our work suggests that these stage data analysis protocols can refine or otherwise reduce variability in these indices by accounting for relevant factors such as rainfall forcing. These protocols for derivation of transient indices should be tested for their potential to improve correlations between stream hydrology and temporally aligned biotic data and dissolved nutrient fluxes in streams.     

OPTIMIZING WELL PLACEMENT IN A COASTAL AQUIFER: OUTER CAPE COD,MASSACHUSETTS1

ABSTRACT: An inverse‐simulation approach is used to determine optimal strategies for developing public water‐supply systems in a shallow, coastal aquifer on the outermost arm of the Cape Cod peninsula in Massachusetts. Typically a forward simulation (or “trial and error”) approach is used to find best pumping strategies, but the chances of finding success with this tact diminish as the number of potential options grows large. Well locations and pumping rates are optimized with respect to: (1) providing sufficient water to areas of water‐quality impairment, (2) minimizing impacts to nearby surface waters, (3) preventing saltwater contamination due to overpumping, and (4) minimizing financial cost of well development. Potential well sites and water‐supply scenarios are separated into “politically‐based” and “resource‐based” categories to gain insight into the degree that pre‐existing political boundaries hinder best management practices. The approach provides a promising tool in transboundary water‐resources settings because it allows stakeholders to find solutions that best meet everyone's goals, as opposed to pursuing options that will create conflict, or are less than optimal.     

PREDICTING INFLUENCES OF URBAN DEVELOPMENT ON THERMAL HABITAT IN A WARM WATER STREAM1

ABSTRACT: Watershed and aquatic ecosystem management requires methods to predict and understand thermal impacts on stream habitat from urbanization. This study evaluates thermal effects of projected urbanization using a modeling framework and considers the biological implications to the fish community. The Stream Network Temperature Model (SNTEMP) was used in combination with the Hydrologic Simulation Program Fortran (HSPF) to assess changes in stream thermal habitat under altered stream‐ flow, shade, and channel width associated with low, medium, and high density urban developments in the Back Creek watershed (Roanoke County, Virginia). Flow alteration by the high density development scenario alone caused minimal heating of mean daily summer base flow (mean +0.1°C). However, when flow changes were modeled concurrently with reduced shade and increased channel width, mean daily temperature increased 1°C. Maximum daily temperatures exceeding the state standard (31°C) increased from 1.1 to 7.6 percent of the time using summer 2000 climatic conditions. Model results suggest that additional urban development will alter stream temperature, potentially limiting thermal habitat and shifting the fish community structure from intolerant to tolerant fish species in Back Creek. More research is needed on the sub‐lethal or chronic effects of increased stream temperature regimes on fish, particularly for those species already living in habitats near their upper limits.