FORECASTING DRY SEASON STREAMFLOW ON THE PEACE RIVER AT ARCADIA,FLORIDA, USA1

ABSTRACT: The Peace River at Arcadia, Florida, is a municipal water supply supplement for southwestern Florida. Consequently, probabilities of encountering low flows during the dry season are of critical importance. Since the association between Pacific Ocean sea surface temperatures (SSTs) and seasonal streamflow variability in the southeastern United States is well documented, it is reasonable to generate forecasts based on this information. Here, employing historic records of minimum, mean, and maximum flows during winter (JFM) and spring (AMJ), upper and lower terciles define “above normal,”“normal,” and “below normal” levels of each variable. A probability distribution model describes the likelihood of these seasonal variables conditioned upon Pacific SSTs from the previous summer (JAS). Model calibration is based upon 40 (of 50) years of record employing stratified random sampling to ensure equal representation from each decade. The model is validated against the remaining 10 samples and the process repeated 100 times. Each conditional probability distribution yields varying probabilities of observing flow variables within defined categories. Generally, a warm (cold) Pacific is associated with higher (lower) flows. To test model skill, the forecast is constrained to be the most probable category in each calibration year, with significance tested by chi‐square frequency tables. For all variables, the tables indicate high levels of association between forecast and observed terciles and forecast skill, particularly during winter. During spring the pattern is less clear, possibly due to the variable starting date of the summer rainy season. This simple technique suggests that Pacific SSTs provide a good forecast of low flows.     

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.     

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.     

Morphology Characteristics of Southern Appalachian Wilderness Streams1

Zink, Jason M., Gregory D. Jennings, and G. Alexander Price, 2012. Morphology Characteristics of Southern Appalachian Wilderness Streams. Journal of the American Water Resources Association (JAWRA) 48(4): 762‐773. DOI: 10.1111/j.1752‐1688.2012.00647.x Abstract: Watersheds without urbanization or impacts from logging are rare in the southern Appalachian Mountains. The Joyce Kilmer/Slickrock Wilderness of North Carolina and Tennessee contains 24 km² of old‐growth forest, with the balance of the wilderness in a mature second‐growth forest. The watersheds of Little Santeetlah and Slickrock Creek are located within the wilderness. Morphological information, including channel dimensions and longitudinal profiles, was gathered from 14 alluvial stream reaches in these watersheds. The study sites had drainage areas from 0.25 to 41.6 km² and stream slopes from 0.014 to 0.104 m/m. Bankfull cross‐section dimensions of the study stream reaches were strongly correlated to drainage area across the observed range of slopes and bed morphology. Cross‐section area and width relationships for the streams in this study did not differ significantly from regional curves for the mountain physiographic region of North Carolina. Observations of these reaches did not suggest a definitive rule regarding the proportion of steps and riffles in streams. Pools occupied greater than 50% of the length in all stream reaches with slopes less than 0.07 m/m. Significant correlation existed between step height ratio and slope, suggesting that step height can be approximated as the product of channel width and slope. Riffle length and riffle slope ratios were also significantly correlated with slope, though pool spacing was not.     

INFLUENCES OF A FOREST ON THE HYDRAULIC GEOMETRY OF TWO MOUNTAIN STREAMS1

ABSTRACT The influence of a forest on the formation of steps in two small streams of the Colorado Rocky Mountains was studied. Steps provided by logs fallen across the channel added to flow energy reduction. The streams required additional gravel bars to adjust to slope. Average step length between logs and gravel bars was strongly related to channel gradient and median bed material size. Based on the average number of log steps per 50 feet of channel, an average of 116 percent of gravel bars were added at Fool Creek and 60 percent at Deadhorse Creek. The latter had 52 percent more logs in the channel and therefore required less bed material movement than the former. Although these are “rushing mountain streams,” most flow velocities ranged between 0.5 and 2.5 f.p.s. Exponents of a function relating rate of change of depth or velocity to discharge indicated that dynamic stream equilibrium was attained. Implications for forest management are that sanitation cuts (removal of dead and dying trees) would not be permissible where a stream is in dynamic equilibrium and bed material movement should be minimized.     

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.     

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.     

Improved Bankfull Channel Geometry Prediction Using Two‐Year Return‐Period Discharge1

He, Laien and Gregory V. Wilkerson, 2011. Improved Bankfull Channel Geometry Prediction Using Two‐Year Return‐Period Discharge. Journal of the American Water Resources Association (JAWRA) 47(6):1298–1316. DOI: 10.1111/j.1752‐1688.2011.00567.x Abstract:  Bankfull discharge (Q_bf) and bankfull channel geometry (i.e., width, W_bf; mean depth, D_bf; and cross‐section area, A_bf) are important design parameters in stream restoration, habitat creation, mined land reclamation, and related projects. The selection of values for these parameters is facilitated by regional curves (regression models in which Q_bf, W_bf, D_bf, and A_bf are predicted as a function of drainage area, A_da). This paper explores the potential for the two‐year return‐period discharge (Q₂) to improve predictions of W_bf, D_bf, and A_bf. Improved predictions are expected because Q₂ estimates integrate the effects of basin drainage area, climate, and geology. For conducting this study, 29 datasets (each representing one hydrologic region) spanning 14 states in the United States were analyzed. We assessed the utility of using Q₂ by comparing statistical measures of regression model performance (e.g., coefficient of determination and Akaike’s information criterion). Compared to using A_da, Q₂ is shown to be a “clearly superior” predictor of W_bf, D_bf, and A_bf, respectively, for 21, 13, and 25% of the datasets. By contrast, A_da yielded a clearly superior model for predicting W_bf, D_bf, and A_bf, respectively, for 0, 0, and 14% of the datasets. Our conclusion is that it alongside with developing conventional regional curves using A_da it is prudent to develop regional curves that use Q₂ as an independent variable because in some cases the resulting model will be superior.     

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.     

ANALYSIS OF STREAMFLOW TRENDS AND THE EFFECTS OF CLIMATE IN PENNSYLVANIA, 1971 TO 20011

An understanding of temporal trends in total stream‐flow (TSF), base flow (BF), and storm runoff (RO) can help in the development of water management plans for watersheds and local communities. In this study, 47 streams across Pennsylvania that were unregulated and unaffected by karst environments or coal mining were studied for flow trends and their relationships to selected climate parameters for the period 1971 to 2001. LOWESS curves for annual flow showed that almost all of the selected streams in Pennsylvania had downward trends in total TSF, BF, and RO. Using a seasonal Mann‐Kendall analysis, downward trends were significant at an α= 0.05 level for 68, percent 70 percent, and 62 percent of the streams and at an α= 0.10 level for another 19, 17, and 13 percent of the streams for TSF, BF, and RO, respectively. The ratio of BF to TSF (RBS) had significant upward trends for 34 percent of the streams at an α= 0.05 level and for another 9 percent of the streams at an α= 0.10 level, indicating that TSF decreased relative to BF for more than 40 percent of the streams during the previous 30 years. Downward trends in TSF, BF, and RO were most common for the months of June, July, and December. Trend analyses using monthly and annual total precipitation and mean temperature showed some association between climate and the streamflow trends, but Spearman's correlation and partial Mann‐Kendall analyses revealed that the trends in TSF, BF, and RO could not be explained by trends in precipitation and temperature alone, and thus urbanization and development may have played a role.     

Legacy Effects of Colonial Millponds on Floodplain Sedimentation,Bank Erosion,and Channel Morphology,Mid‐Atlantic,USA1

Abstract: Many rivers and streams of the Mid‐Atlantic Region, United States (U.S.) have been altered by postcolonial floodplain sedimentation (legacy sediment) associated with numerous milldams. Little Conestoga Creek, Pennsylvania, a tributary to the Susquehanna River and the Chesapeake Bay, is one of these streams. Floodplain sedimentation rates, bank erosion rates, and channel morphology were measured annually during 2004‐2007 at five sites along a 28‐km length of Little Conestoga Creek with nine colonial era milldams (one dam was still in place in 2007). This study was part of a larger cooperative effort to quantify floodplain sedimentation, bank erosion, and channel morphology in a high sediment yielding region of the Chesapeake Bay watershed. Data from the five sites were used to estimate the annual volume and mass of sediment stored on the floodplain and eroded from the banks for 14 segments along the 28‐km length of creek. A bank and floodplain reach based sediment budget (sediment budget) was constructed for the 28 km by summing the net volume of sediment deposited and eroded from each segment. Mean floodplain sedimentation rates for Little Conestoga Creek were variable, with erosion at one upstream site (?5 mm/year) to deposition at the other four sites (highest = 11 mm/year) despite over a meter of floodplain aggradation from postcolonial sedimentation. Mean bank erosion rates range between 29 and 163 mm/year among the five sites. Bank height increased 1 m for every 10.6 m of channel width, from upstream to downstream (R² = 0.79, p < 0.0001) resulting in progressively lowered hydraulic connectivity between the channel and the floodplain. Floodplain sedimentation and bank erosion rates also appear to be affected by the proximity of the segments to one existing milldam, which promotes deposition upstream and scouring downstream. The floodplain and bank along the 28‐km reach produced a net mean sediment loss of 5,634 Mg/year for 2004‐2007, indicating that bank erosion was exceeding floodplain sedimentation. In particular, the three segments between the existing dam and the confluence with the Conestoga River (32% of the studied reach) account for 97% of the measured net sediment budget. Future research directed at understanding channel equilibria should facilitate efforts to reduce the sediment impacts of dam removal and legacy sediment.     

RIPARIAN COMMUNITIES ASSOCIATED WITH PACIFIC NORTHWEST HEADWATER STREAMS: ASSEMBLAGES,PROCESSES, AND UNIQUENESS1

Riparian areas of large streams provide important habitat to many species and control many instream processes — but is the same true for the margins of small streams? This review considers riparian areas alongside small streams in forested, mountainous areas of the Pacific Northwest and asks if there are fundamental ecological differences from larger streams and from other regions and if there are consequences for management from any differences. In the moist forests along many small streams of the Pacific Northwest, the contrast between the streamside and upslope forest is not as strong as that found in drier regions. Small streams typically lack floodplains, and the riparian area is often constrained by the hillslope. Nevertheless, riparian‐associated organisms, some unique to headwater areas, are found along small streams. Disturbance of hillslopes and stream channels and microclimatic effects of streams on the riparian area provide great heterogeneity in processes and diversity of habitats. The tight coupling of the terrestrial riparian area with the aquatic system results from the closed canopy and high edge‐to‐area ratio for small streams. Riparian areas of the temperate, conifer dominated forests of the Pacific Northwest provide a unique environment. Forest management guidelines for small streams vary widely, and there has been little evaluation of the local or downstream consequences of forest practices along small streams.     

CUMULATIVE EFFECTIVE STREAM POWER AND BANK EROSION ON THE SACRAMENTO RIVER,CALIFORNIA, USA1

ABSTRACT: Bank erosion along a river channel determines the pattern of channel migration. Lateral channel migration in large alluvial rivers creates new floodplain land that is essential for riparian vegetation to get established. Migration also erodes existing riparian, agricultural, and urban lands, sometimes damaging human infrastructure (e.g., scouring bridge foundations and endangering pumping facilities) in the process. Understanding what controls the rate of bank erosion and associated point bar deposition is necessary to manage large alluvial rivers effectively. In this study, bank erosion was proportionally related to the magnitude of stream power. Linear regressions were used to correlate the cumulative stream power, above a lower flow threshold, with rates of bank erosion at 13 sites on the middle Sacramento River in California. Two forms of data were used: aerial photography and field data. Each analysis showed that bank erosion and cumulative effective stream power were significantly correlated and that a lower flow threshold improves the statistical relationship in this system. These correlations demonstrate that land managers and others can relate rates of bank erosion to the daily flow rates of a river. Such relationships can provide information concerning ecological restoration of floodplains related to channel migration rates as well as planning that requires knowledge of the relationship between flow rates and bank erosion rates.     

A Multi‐Agency Nutrient Dataset Used to Estimate Loads,Improve Monitoring Design,and Calibrate Regional Nutrient SPARROW Models1

Saad, David A., Gregory E. Schwarz, Dale M. Robertson, and Nathaniel L. Booth, 2011. A Multi‐Agency Nutrient Dataset Used to Estimate Loads, Improve Monitoring Design, and Calibrate Regional Nutrient SPARROW Models. Journal of the American Water Resources Association (JAWRA) 47(5):933‐949. DOI: 10.1111/j.1752‐1688. 2011.00575.x Abstract: Stream‐loading information was compiled from federal, state, and local agencies, and selected universities as part of an effort to develop regional SPAtially Referenced Regressions On Watershed attributes (SPARROW) models to help describe the distribution, sources, and transport of nutrients in streams throughout much of the United States. After screening, 2,739 sites, sampled by 73 agencies, were identified as having suitable data for calculating long‐term mean annual nutrient loads required for SPARROW model calibration. These sites had a wide range in nutrient concentrations, loads, and yields, and environmental characteristics in their basins. An analysis of the accuracy in load estimates relative to site attributes indicated that accuracy in loads improve with increases in the number of observations, the proportion of uncensored data, and the variability in flow on observation days, whereas accuracy declines with increases in the root mean square error of the water‐quality model, the flow‐bias ratio, the number of days between samples, the variability in daily streamflow for the prediction period, and if the load estimate has been detrended. Based on compiled data, all areas of the country had recent declines in the number of sites with sufficient water‐quality data to compute accurate annual loads and support regional modeling analyses. These declines were caused by decreases in the number of sites being sampled and data not being entered in readily accessible databases.     

The Impact of Asynchronicity on Event‐Flow Estimation in Basin‐Scale Hydrologic Model Calibration1

Abstract: The calibration of basin‐scale hydrologic models consists of adjusting parameters such that simulated values closely match observed values. However, due to inevitable inaccuracies in models and model inputs, simulated response hydrographs for multiyear calibrations will not be perfectly synchronized with observed response hydrographs at the daily time step. An analytically derived formula suggests that when timing errors are significant, traditional calibration approaches may generally underestimate the total event‐flow volume. An event‐adaptive time series is developed and incorporated into the Nash‐Sutcliffe Efficiency objective function to diagnose the potential impact of event‐flow synchronization errors. Test sites are the 50 km² Subwatershed I of the Little River Experimental Watershed (LREWswI) in southeastern Georgia, and the 610 km² Little Washita River Experimental Watershed (LWREW) in southwestern Oklahoma, with the Soil and Water Assessment Tool used as the hydrologic model. Results suggest that simulated surface runoff generation is 55% less for LREWswI when the daily time series is used compared with when the event‐adaptive technique is used. Event‐flow generation may also be underestimated for LWREW, but to a lesser extent than it may be for LREWswI, due to a larger portion of the event flow being lateral flow.     

Temporal and Spatial Trends in Nutrient and Sediment Loading to Lake Tahoe,California‐Nevada,USA

Since 1980, the Lake Tahoe Interagency Monitoring Program (LTIMP) has provided stream‐discharge and water quality data—nitrogen (N), phosphorus (P), and suspended sediment—at more than 20 stations in Lake Tahoe Basin streams. To characterize the temporal and spatial patterns in nutrient and sediment loading to the lake, and improve the usefulness of the program and the existing database, we have (1) identified and corrected for sources of bias in the water quality database; (2) generated synthetic datasets for sediments and nutrients, and resampled to compare the accuracy and precision of different load calculation models; (3) using the best models, recalculated total annual loads over the period of record; (4) regressed total loads against total annual and annual maximum daily discharge, and tested for time trends in the residuals; (5) compared loads for different forms of N and P; and (6) tested constituent loads against land use‐land cover (LULC) variables using multiple regression. The results show (1) N and P loads are dominated by organic N and particulate P; (2) there are significant long‐term downward trends in some constituent loads of some streams; and (3) anthropogenic impervious surface is the most important LULC variable influencing water quality in basin streams. Many of our recommendations for changes in water quality monitoring and load calculation methods have been adopted by the LTIMP.     

AN ANALYSIS OF REGION‐OF‐INFLUENCE METHODS FOR FLOOD REGIONALIZATION IN THE GULF‐ATLANTIC ROLLING PLAINS1

ABSTRACT: Region‐of‐influence (RoI) approaches for estimating stream flow characteristics at ungaged sites were applied and evaluated in a case study of the 50‐year peak discharge in the Gulf‐Atlantic Rolling Plains of the southeastern United States. Linear regression against basin characteristics was performed for each ungaged site considered based on data from a region of influence containing the n closest gages in predictor variable (PRoI) or geographic (GRoI) space. Augmentation of this count based cutoff by a distance based cutoff also was considered. Prediction errors were evaluated for an independent (split‐sampled) dataset. For the dataset and metrics considered here: (1) for either PRoI or GRoI, optimal results were found when the simpler count based cutoff, rather than the distance augmented cutoff, was used; (2) GRoI produced lower error than PRoI when applied indiscriminately over the entire study region; (3) PRoI performance improved considerably when Rol was restricted to predefined geographic subregions.     

Modeling Pre‐ and Post‐Dam Removal Sediment Dynamics: The Kalamazoo River,Michigan1

Abstract:  The state of Michigan is interested in removing two low‐head dams in an 8.8 km reach of the Kalamazoo River between Plainwell and Otsego, Michigan, while minimizing impacts locally and to downstream reaches. The study was designed to evaluate the erosion, transport, and deposition of sediments over a 37.3‐year period using the channel evolution model CONCEPTS for three simulation scenarios: Dams In (DI), Dams Out (DO), and Design (D). The total mass of sediment emanating from the channel boundary, for the DI case, shows net deposition of 4,100 T/y for the study reach, with net transport (suspended and bed load) of 10,500 T/y passing the downstream boundary. For the DO case, net erosion is 19,200 T/y with net transport of 30,100 T/y (187% increase) passing the downstream boundary. For the D case, net deposition is 2,570 T/y (37% decrease) with transport of 14,200 T/y (35% increase) passing the downstream boundary. The most significant findings were: (1) removal of the low‐head dams will cause significant erosion of sediments stored behind the dams and increased sediment loads passing the downstream boundary and (2) sediment loads for the proposed channel design are similar to existing conditions and offer reduced fine‐sediment loadings.     

Response of Algal Biomass to Large‐Scale Nutrient Controls in the Clark Fork River,Montana, United States1

Suplee, Michael W., Vicki Watson, Walter K. Dodds, and Chris Shirley, 2012. Response of Algal Biomass to Large‐Scale Nutrient Controls in the Clark Fork River, Montana, United States. Journal of the American Water Resources Association (JAWRA) 48(5): 1008‐1021. DOI: 10.1111/j.1752‐1688.2012.00666.x Abstract: Nutrient pollution is an ongoing concern in rivers. Although nutrient targets have been proposed for rivers, little is known about long‐term success of programs to decrease river nutrients and algal biomass. Twelve years of summer data (1998‐2009) collected along 383 km of the Clark Fork River were analyzed to ascertain whether a basin‐wide nutrient reduction program lowered ambient total nitrogen (TN) and total phosphorus (TP) concentrations, and bottom‐attached algal biomass. Target nutrient and algal biomass levels were established for the program in 1998. Significant declines were observed in TP but not TN along the entire river. Downstream of the city of Missoula, TP declined below a literature‐derived TP saturation breakpoint and met program targets after 2005; TN was below targets since 2007. Algal biomass also declined significantly below Missoula. Trends there likely relate to the city’s wastewater facility upgrades, despite its 20% population increase. Upstream of Missoula, nutrient reductions were less substantial; still, TP and TN declined toward saturation breakpoints, but no significant reductions in algal biomass occurred, and program targets were not met. The largest P‐load reduction to the river was from a basin‐wide phosphate laundry detergent ban set 10 years before, in 1989. We document that nutrient reductions in rivers can be successful in controlling algal biomass, but require achievement of concentrations below saturation and likely close to natural background.     

POST‐HARVEST RIPARIAN BUFFER RESPONSE: IMPLICATIONS FOR WOOD RECRUITMENT MODELING AND BUFFER DESIGN1

Despite the importance of riparian buffers in providing aquatic functions to forested streams, few studies have sought to capture key differences in ecological and geomorphic processes between buffered sites and forested conditions. This study examines post‐harvest buffer conditions from 20 randomly selected harvest sites within a managed tree farm in the Cascade Mountains of western Washington. Post‐harvest wind derived treefall rates in buffers up to three years post‐harvest averaged 268 trees/km/year, 26 times greater than competition‐induced mortality rate estimates. Treefall rates and stem breakage were strongly tied to tree species and relatively unaffected by stream direction. Observed treefall direction is strongly biased toward the channel, irrespective of channel or buffer orientation. Fall direction bias can deliver significantly more wood recruitment relative to randomly directed treefall, suggesting that models that utilize the random fall assumption will significantly underpredict recruitment. A simple estimate of post‐harvest wood recruitment from buffers can be obtained from species specific treefall and breakage rates, combined with bias corrected recruitment probability as a function of source distance from the channel. Post‐harvest wind effects may reduce the standing density of trees enough to significantly reduce or eliminate competition mortality and thus indirectly alter bank erosion rates, resulting in substantially different wood recruitment dynamics from buffers as compared to unmanaged forests.