Cropland Riparian Buffers throughout Chesapeake Bay Watershed: Spatial Patterns and Effects on Nitrate Loads Delivered to Streams

We used statistical models to provide the first empirical estimates of riparian buffer effects on the cropland nitrate load to streams throughout the Chesapeake Bay watershed. For each of 1,964 subbasins, we quantified the 1990 prevalence of cropland and riparian buffers. Cropland was considered buffered if the topographic flow path connecting it to a stream traversed a streamside forest or wetland. We applied a model that predicts stream nitrate concentration based on physiographic province and the watershed proportions of unbuffered and buffered cropland. We used another model to predict annual streamflow based on precipitation and temperature, and then multiplied the predicted flows and concentrations to estimate 1990 annual nitrate loads. Across the entire Chesapeake watershed, croplands released 92.3 Gg of nitrate nitrogen, but 19.8 Gg of that was removed by riparian buffers. At most, 29.4 Gg more might have been removed if buffer gaps were restored so that all cropland was buffered. The other 43.1 Gg of cropland load cannot be addressed with riparian buffers. The Coastal Plain physiographic province provided 52% of the existing buffer reduction of Bay‐wide nitrate loads and 36% of potential additional removal from buffer restoration in cropland buffer gaps. Existing and restorable nitrate removal in buffers were lower in the other three major provinces because of less cropland, lower buffer prevalence, and lower average buffer nitrate removal efficiency.     

PHYSICAL DETERMINANTS OF ECONOMIC VALUE OF RIPARIAN BUFFERS IN AN AGRICULTURAL WATERSHED1

ABSTRACT: Economic values of riparian buffers in a watershed are evaluated by the changes in the net economic return for farming with and without riparian buffers when achieving the same water quality objectives. The simulated water quality impacts of alternative farming systems using SWAT and experimental data for riparian buffers are used in a mathematical optimization model to estimate net economic return for farming subject to a water quality objective. Physical characteristics such as stream length, channel slope, average land slope, cropland percentage and several soil attributes are identified in the watershed using ARC/INFO GIS. A regression model is then used to evaluate the impacts of these physical characteristics on the estimated economic values of buffers. The study is conducted in Goodwater Creek watershed, Missouri. The results show the estimated economic value of buffers is significantly affected by some soil properties, stream length, and cropland percentage in watershed and can be used to improve the effectiveness of riparian buffers at watershed and regional levels.     

Nonpoint Source Pollution Responses Simulation for Conversion Cropland to Forest in Mountains by SWAT in China

Several environmental protection policies have been implemented to prevent soil erosion and nonpoint source (NPS) pollutions in China. After severe Yangtze River floods, the “conversion cropland to forest policy” (CCFP) was carried out throughout China, especially in the middle and upper reaches of Yangtze River. The research area of the current study is located in Bazhong City, Sichuan Province in Yangtze River watershed, where soil erosion and NPS pollution are serious concerns. Major NPS pollutants include nitrogen (N) and phosphorus (P). The objective of this study is to evaluate the long-term impact of implementation of the CCFP on stream flow, sediment yields, and the main NPS pollutant loading at watershed level. The Soil and Water Assessment Tool (SWAT) is a watershed environmental model and is applied here to simulate and quantify the impacts. Four scenarios are constructed representing different patterns of conversion from cropland to forest under various conditions set by the CCFP. Scenario A represented the baseline, i.e., the cropland and forest area conditions before the implementation of CCFP. Scenario B represents the condition under which all hillside cropland with slope larger than 25° was converted into forest. In scenario C and D, hillside cropland with slope larger than 15° and 7.5° was substituted by forest, respectively. Under the various scenarios, the NPS pollution reduction due to CCFP implementation from 1996–2005 is estimated by SWAT. The results are presented as percentage change of water flow, sediment, organic N, and organic P at watershed level. Furthermore, a regression analysis is conducted between forest area ratio and ten years’ average NPS load estimations, which confirmed the benefits of implementing CCFP in reducing nonpoint source pollution by increasing forest area in mountainous areas. The reduction of organic N and organic P is significant (decrease 42.1% and 62.7%, respectively) at watershed level. In addition, this study also proves that SWAT modeling approach can be used to estimate NPS pollutants’ impacts of land use conversions in large watershed.     

Nitrous oxide accumulation in soils from riparian buffers of a coastal plain watershed carbon/nitrogen ratio control

Riparian buffers are used throughout the world for the protection of water bodies from nonpoint-source nitrogen pollution. Few studies of riparian or treatment wetland denitrification consider the production of nitrous oxide (N2O). The objectives of this research were to ascertain the level of potential N2O production in riparian buffers and identify controlling factors for N2O accumulations within riparian soils of an agricultural watershed in the southeastern Coastal Plain of the USA. Soil samples were obtained from ten sites (site types) with different agronomic management and landscape position. Denitrification enzyme activity (DEA) was measured by the acetylene inhibition method. Nitrous oxide accumulations were measured after incubation with and without acetylene (baseline N2O production). The mean DEA (with acetylene) was 59 microg N2O-N kg(-1) soil h(-1) for all soil samples from the watershed. If no acetylene was added to block conversion of N2O to N2, only 15 microg N2O-N kg(-1) soil h(-1) were accumulated. Half of the samples accumulated no N2O. The highest level of denitrification was found in the soil surface layers and in buffers impacted by either livestock waste or nitrogen from legume production. Nitrous oxide accumulations (with acetylene inhibition) were correlated to soil nitrogen (r2=0.59). Without acetylene inhibition, correlations with soil and site characteristics were lower. Nitrous oxide accumulations were found to be essentially zero, if the soil C/N ratios>25. Soil C/N ratios may be an easily measured and widely applicable parameter for identification of potential hot spots of N2O productions from riparian buffers.     

A Geospatial Methodology to Identify Locations of Concentrated Runoff from Agricultural Fields

A geospatial methodology has been developed that utilizes high resolution lidar‐derived DEMs to help track runoff from agricultural fields and identify areas of potential concentrated flow through vegetated riparian areas in the Coastal Plain of Virginia. Points of concentrated flow are identified across 74 agricultural fields within the Virginia portion of the Chesapeake Bay watershed. On average, 70% of the surface area of the agricultural fields analyzed drains through less than 20 m of the field margin, and on average 81% of the field surface area drains through 1% or less of the field margin. Within the riparian buffer, locations that were predicted by the geospatial model to have high levels of concentrated flow were found to exhibit evidence of channelization. Results indicate that flow concentration and channelized flow through vegetated riparian areas may be common along the margin of agricultural fields, resulting in vegetated riparian areas that are less effective at sediment trapping than assumed. Additional results suggest that the regulations governing the location and width of vegetated riparian may not be sufficient to achieve goals for reducing sediment and nutrient runoff from nonpoint agricultural sources. Combined with the increasing availability of lidar‐derived DEMs, the geospatial model presented has the potential to advance management practices aimed at reducing nonpoint source pollution leaving agricultural fields.     

Stream Nitrogen Sources Apportionment and Pollution Control Scheme Development in an Agricultural Watershed in Eastern China

A modeling system that couples a land-use-based export coefficient model, a stream nutrient transport equation, and Bayesian statistics was developed for stream nitrogen source apportionment. It divides a watershed into several sub-catchments, and then considers the major land-use categories as stream nitrogen sources in each sub-catchment. The runoff depth and stream water depth are considered as the major factors influencing delivery of nitrogen from land to downstream stream node within each sub-catchment. The nitrogen sources and delivery processes are lumped into several constant parameters that were calibrated using Bayesian statistics from commonly available stream monitoring and land-use datasets. This modeling system was successfully applied to total nitrogen (TN) pollution control scheme development for the ChangLe River watershed containing six sub-catchments and four land-use categories. The temporal (across months and years) and spatial (across sub-catchments and land-use categories) variability of nonpoint source (NPS) TN export to stream channels and delivery to the watershed outlet were assessed. After adjustment for in-stream TN retention, the time periods and watershed areas with disproportionately high-TN contributions to the stream were identified. Aimed at a target stream TN level of 2 mg L^?1, a quantitative TN pollution control scheme was further developed to determine which sub-catchments, which land-use categories in a sub-catchment, which time periods, and how large of NPS TN export reduction were required. This modeling system provides a powerful tool for stream nitrogen source apportionment and pollution control scheme development at the watershed scale and has only limited data requirements.     

Effect of Watershed Subdivision and Filter Width on SWAT Simulation of a Coastal Plain Watershed1

Cho, Jaepil, Richard R. Lowrance, David D. Bosch, Timothy C. Strickland, Younggu Her, and George Vellidis, 2010. Effect of Watershed Subdivision and Filter Width on SWAT Simulation of a Coastal Plain Watershed. Journal of the American Water Resources Association (JAWRA) 46(3):586-602. DOI: 10.1111/j.1752-1688.2010.00436.x Abstract: The Soil and Water Assessment Tool (SWAT) does not fully simulate riparian buffers, but has a simple filter function that is responsive to filter strip width (FILTERW). The objectives of this study were to (1) evaluate SWAT hydrology and water quality response to changes in watershed subdivision levels and different FILTERW configurations and (2) provide guidance for selecting appropriate watershed subdivision for model runs that include the riparian buffer feature through the FILTERW parameter. Watershed subdivision level is controlled by the critical source area (CSA) which defines the minimum drainage area required to form the origin of a stream. SWAT was calibrated on a 15.7 km² subdrainage within the Little River Experimental Watershed, Georgia. The calibrated parameter set was applied to 32 watershed configurations consisting of four FILTERW representations for each of eight CSA levels. Streamflow predictions were stable regardless of watershed subdivision and FILTERW configuration. Predicted sediment and nutrient loads from upland areas decreased as CSA increased when spatial variations of riparian buffers are considered. Sediment and nutrient yield at the watershed outlet was responsive to different combinations of CSA and FILTERW depending on selected in-stream processes. CSA ranges which provide stable sediment and nutrient yields at the watershed outlet was suggested for avoiding significant modifications in selected parameter set.     

ECONOMIC EVALUATION OF RIPARIAN BUFFERS IN AN AGRICULTURAL WATERSHED1

ABSTRACT: This study determines the most cost effective spatial pattern of farming systems for improving water quality and evaluates the economic value of riparian buffers in reducing agricultural nonpoint source pollution in a Midwestern agricultural watershed. Economic and water quality impacts of alternative farming systems are evaluated using the CARE and SWAT models, respectively. The water quality benefits of riparian buffers are estimated by combining experimental data and simulated water quality impacts of fanning systems obtained using SWAT. The net economic value of riparian buffers in improving water quality is estimated by total watershed net return with riparian buffers minus total watershed net return without riparian buffers minus the opportunity cost of riparian buffers. Exclusive of maintenance cost, the net economic value of riparian buffers in reducing atrazine concentration from 45 to 24 ppb is $612,117 and the savings in government cost is $631,710. Results strongly support efforts that encourage farmers to develop or maintain riparian buffers adjacent to streams.     

Tracing nitrate transport and environmental impact from intensive swine farming using delta nitrogen-15

Natural-abundance delta15N showed that nitrate generated from commercial land application of swine (Sus scrofa domesticus) waste within a North Carolina Coastal Plain catchment was being discharged to surface waters by ground water passing beneath the sprayfields and adjacent riparian buffers. This was significant because intensive swine farms in North Carolina are considered non-discharge operations, and riparian buffers with minimum widths of 7.6 m (25 ft) are the primary regulatory control on ground water export of nitrate from these operations. This study shows that such buffers are not always adequate to prevent discharge of concentrated nitrate in ground water from commercial swine farms in the Mid-Atlantic Coastal Plain, and that additional measures are required to ensure non-discharge conditions. The median delta15N-total N of liquids in site swine waste lagoons was +15.4 +/- 0.2% vs. atmospheric nitrogen. The median delta15N-NO3 values of shallow ground water beneath and adjacent to site sprayfields, a stream draining sprayfields, and waters up to 1.5 km downstream were + 15.3 +/- 0.2 to + 15.4 +/- 0.2%. Seasonal and spatial isotopic variations in lagoons and well waters were greatly homogenized during ground water transport and discharge to streams. Neither denitrification nor losses of ammonia during spraying significantly altered the bulk ground water delta15N signal being delivered to streams. The lagoons were sources of chloride and potassium enrichment, and shallow ground water showed strong correlation between nitrate N, potassium, and chloride. The 15N-enriched nitrate in ground water beneath swine waste sprayfields can thus be successfully traced during transport and discharge into nearby surface waters.     

Stream buffer effectiveness in an agriculturally influenced area, southwestern Georgia: responses of water quality, macroinvertebrates, and amphibians

To determine useful metrics for assessing stream water quality in the Southeastern Coastal Plain, we examined differences among two buffered and three unbuffered streams in an agricultural landscape in southwestern Georgia. Potential indicators included amphibian diversity and abundance, aquatic macroinvertebrate populations, riparian vegetative structure, water quality, and stream physical parameters. Variability among sites and treatments (buffered vs. unbuffered) existed, with sites in the same treatment as most similar, and disturbances from a nearby eroding gully strongly affecting one unbuffered site. Of the invertebrate metrics examined, percentages of clingers, Ephemeroptera-Plecoptera-Trichoptera (EPT), Elmidae (Coleoptera), Crustacea (Decapoda and Amphipoda), and dipterans were found to be possible indicators of stream health for perennial streams within this region. Overall, buffered sites showed higher percentages of sensitive invertebrate groups and showed lower and more stable concentrations of nitrate N, suspended solids, and fecal coliforms (FCs). Percent canopy cover was similar among sites; however, riparian vegetative coverage and percent leaf litter were greatest at buffered sites. No differences in amphibian abundance, presence, and absence within the riparian area were apparent between sites; however, instream larval salamanders were more abundant at buffered streams. In this study, stream buffers appeared to decrease nutrient and sediment loads to adjacent streams, enhancing overall water quality. Selected benthic macroinvertebrate metrics and amphibian abundance also appeared sensitive to agricultural influences. Amphibians show potential as indicator candidates, however further information is needed on their responses and tolerances to disturbances from the microhabitat to landscape levels.     

Riparian buffers and potentially unstable ground

The spatial coincidence between riparian buffers of various widths and extents and potentially unstable ground was quantified using a physically based model for shallow landslide initiation and GIS for two watersheds on the Olympic Peninsula, Washington, USA. The proportion of the potentially unstable ground in each watershed within riparian buffers is a function of both buffer width and the extent of the stream channel network being buffered. While current buffers required by Washington State cover less than 5% of the potentially unstable ground, buffering all stream channels in these watersheds with 100-m buffers covered 75%–90% of the potentially unstable areas. Our analyses further show that: (1) riparian buffers are not efficient mechanisms for protecting potentially unstable ground, and (2) identifying potentially unstable ground using a physically based model should prove more effective for designing methods to reduce shallow landsliding hazards than relying on extensive buffer zones along stream channels.     

EFFECTS OF LAND COVER AND GEOLOGY ON STREAM CHEMISTRY IN WATERSHEDS OF CHESAPEAKE BAY1

ABSTRACT: We measured the base‐flow stream chemistry in all the major physiographic provinces of the Chesapeake Bay drainage basin. The spatial variation of stream chemistry was closely related to differences in geology and land cover among the sampled watersheds. Some stream chemistry variables were strongly affected by geological settings in the watersheds while others were more influenced by land cover. The effects of land cover differed among chemical constituents and regions. Concentrations of Ca²⁺, Mg²⁺, pH, total alkalinity, and conductivity were mainly functions of carbonate bedrock, especially in the Great Valley. Nitrate‐N and total dissolved N were closely related to cropland and increased as the percentage of cropland increased. The rate of increase varied from region to region with the highest in the Piedmont. Na⁺ and Cl^? were mainly affected by the percentage of developed area in a watershed, especially in the Coastal Plain and Piedmont. We observed no significant effects of region or land cover on species of phosphorus because samples were collected under base flow conditions and only dissolved forms were measured. Dissolved silicate (DSi) was not related to any other water chemistry variables. DSi increased as developed area decreased and cropland increased in the Coastal Plain, but these patterns were reversed in the Piedmont. There was no consistent pattern in the spatial variation of land cover effects on the reduced forms of N, dissolved organic P, dissolved organic matter, and K⁺.     

Assessing Watershed-Scale, Long-Term Hydrologic Impacts of Land-Use Change Using a GIS-NPS Model

Land-use change, dominated by an increase in urban/impervious areas, has a significant impact on water resources. This includes impacts on nonpoint source (NPS) pollution, which is the leading cause of degraded water quality in the United States. Traditional hydrologic models focus on estimating peak discharges and NPS pollution from high-magnitude, episodic storms and successfully address short-term, local-scale surface water management issues. However, runoff from small, low-frequency storms dominates long-term hydrologic impacts, and existing hydrologic models are usually of limited use in assessing the long-term impacts of land-use change. A long-term hydrologic impact assessment (L-THIA) model has been developed using the curve number (CN) method. Long-term climatic records are used in combination with soils and land-use information to calculate average annual runoff and NPS pollution at a watershed scale. The model is linked to a geographic information system (GIS) for convenient generation and management of model input and output data, and advanced visualization of model results. The L-THIA/NPS GIS model was applied to the Little Eagle Creek (LEC) watershed near Indianapolis, Indiana, USA. Historical land-use scenarios for 1973, 1984, and 1991 were analyzed to track land-use change in the watershed and to assess impacts on annual average runoff and NPS pollution from the watershed and its five subbasins. For the entire watershed between 1973 and 1991, an 18% increase in urban or impervious areas resulted in an estimated 80% increase in annual average runoff volume and estimated increases of more than 50% in annual average loads for lead, copper, and zinc. Estimated nutrient (nitrogen and phosphorus) loads decreased by 15% mainly because of loss of agricultural areas. The L-THIA/NPS GIS model is a powerful tool for identifying environmentally sensitive areas in terms of NPS pollution potential and for evaluating alternative land use scenarios for NPS pollution management.     

Role of wetlands in reducing phosphorus loading to surface water in eight watersheds in the Lake Champlain Basin

A landscape-level approach was applied to eight rural watersheds to assess the role that wetlands play in reducing phosphorus loading to surface waters in the Lake Champlain Basin. Variables summarizing various characteristics of wetlands within a watershed were calculated using a geographic information system and then compared to measured phosphorus loading through multiple regression analyses. The inclusion of a variable based on the area of riparian wetlands located along low- and medium-order streams in conjunction with the area of agricultural and nonwetland forested lands explained 88% of the variance in phosphorus loading to surface waters. The best fit model coefficients (P_load = 0.86Ag + 0.64For – 30Ripwet + 160) suggest that a hectare of riparian wetland may be many times more important in reducing phosphorus than an agricultural hectare is in producing phosphorus. These results provide additional support for the concept that protection of riparian wetlands is an important management strategy for controlling stream water quality in multiuse landscapes.     

Spatial-Scale Effects on Relative Importance of Physical Habitat Predictors of Stream Health

A common theme in recent landscape studies is the comparison of riparian and watershed land use as predictors of stream health. The objective of this study was to compare the performance of reach-scale habitat and remotely assessed watershed-scale habitat as predictors of stream health over varying spatial extents. Stream health was measured with scores on a fish index of biotic integrity (IBI) using data from 95 stream reaches in the Eastern Corn Belt Plain (ECBP) ecoregion of Indiana. Watersheds hierarchically nested within the ecoregion were used to regroup sampling locations to represent varying spatial extents. Reach habitat was represented by metrics of a qualitative habitat evaluation index, whereas watershed variables were represented by riparian forest, geomorphology, and hydrologic indices. The importance of reach- versus watershed-scale variables was measured by multiple regression model adjusted-R² and best subset comparisons in the general linear statistical framework. Watershed models had adjusted-R² ranging from 0.25 to 0.93 and reach models had adjusted-R² ranging from 0.09 to 0.86. Better-fitting models were associated with smaller spatial extents. Watershed models explained about 15% more variation in IBI scores than reach models on average. Variety of surficial geology contributed to decline in model predictive power. Results should be interpreted bearing in mind that reach habitat was qualitatively measured and only fish assemblages were used to measure stream health. Riparian forest and length-slope (LS) factor were the most important watershed-scale variables and mostly positively correlated with IBI scores, whereas substrate and riffle-pool quality were the important reach-scale variables in the ECBP.     

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.     

Landscape Planning for Agricultural Nonpoint Source Pollution Reduction III: Assessing Phosphorus and Sediment Reduction Potential

Riparian buffers have the potential to improve stream water quality in agricultural landscapes. This potential may vary in response to landscape characteristics such as soils, topography, land use, and human activities, including legacies of historical land management. We built a predictive model to estimate the sediment and phosphorus load reduction that should be achievable following the implementation of riparian buffers; then we estimated load reduction potential for a set of 1598 watersheds (average 54 km²) in Wisconsin. Our results indicate that land cover is generally the most important driver of constituent loads in Wisconsin streams, but its influence varies among pollutants and according to the scale at which it is measured. Physiographic (drainage density) variation also influenced sediment and phosphorus loads. The effect of historical land use on present-day channel erosion and variation in soil texture are the most important sources of phosphorus and sediment that riparian buffers cannot attenuate. However, in most watersheds, a large proportion (approximately 70%) of these pollutants can be eliminated from streams with buffers. Cumulative frequency distributions of load reduction potential indicate that targeting pollution reduction in the highest 10% of Wisconsin watersheds would reduce total phosphorus and sediment loads in the entire state by approximately 20%. These results support our approach of geographically targeting nonpoint source pollution reduction at multiple scales, including the watershed scale.     

Assessing Critical Source Areas in Watersheds for Conservation Buffer Planning and Riparian Restoration

A science-based geographic information system (GIS) approach is presented to target critical source areas in watersheds for conservation buffer placement. Critical source areas are the intersection of hydrologically sensitive areas and pollutant source areas in watersheds. Hydrologically sensitive areas are areas that actively generate runoff in the watershed and are derived using a modified topographic index approach based on variable source area hydrology. Pollutant source areas are the areas in watersheds that are actively and intensively used for such activities as agricultural production. The method is applied to the Neshanic River watershed in Hunterdon County, New Jersey. The capacity of the topographic index in predicting the spatial pattern of runoff generation and the runoff contribution to stream flow in the watershed is evaluated. A simple cost-effectiveness assessment is conducted to compare the conservation buffer placement scenario based on this GIS method to conventional riparian buffer scenarios for placing conservation buffers in agricultural lands in the watershed. The results show that the topographic index reasonably predicts the runoff generation in the watershed. The GIS-based conservation buffer scenario appears to be more cost-effective than the conventional riparian buffer scenarios.     

ECONOMIC VALUATION OF RIPARIAN BUFFER AND OPEN SPACE IN A SUBURBAN WATERSHED1

Abstract: This study evaluates the economic value of riparian buffers and open space in a suburban watershed through two nonmarket valuation methods. A contingent valuation survey was implemented in the Dardenne Creek watershed, a suburban watershed of the St. Louis metropolitan area in Missouri, to evaluate the residents' perceptions of and willingness to pay (WTP) for adopting riparian buffers and preserving farmland in a hypothetical real estate market. A hedonic pricing model based on actual sale prices of homes in the watershed was applied to estimate the market value of open space and other environmental conditions such as flood zone and stream proximity in the study area. The results showed that residents' WTP was consistent with the economic values of open space and proximity to streams embedded in existing home prices. Through a better understanding of residents' perceptions and values, riparian buffer and open space programs can be designed and promoted to achieve greater implementation success and environmental benefit.     

THE USE OF BINARY OPTIMIZATION AND HYDROLOGIC MODELS TO FORM RIPARIAN BUFFERS1

ABSTRACT: Riparian buffers are considered important management options for protecting water quality. Land costs and buffer performance, which are functions of local environmental characteristics, are likely to be key attributes in the selection process, especially when budgets are limited. In this article we demonstrate how a framework involving hydrologic models and binary optimization can be used to find the optimal buffer subject to a budget constraint. Two hydrologic models, SWAT and REMM, were used to predict the loads from different source areas with and without riparian buffers. These loads provided inputs for a binary optimization model to select the most cost efficient parcels to form a riparian buffer. This methodology was applied in a watershed in Delaware County, New York. The models were parameterized using readily available digital databases and were later compared against observed flow and water quality data available for the site. As a result of the application of this method, the marginal utility of incremental increases in buffer widths along the stream channel and the set of parcels to form the best affordable riparian buffer were obtained.