Application of Tabu Search Algorithm With a Coupled AnnAGNPS‐CCHE1D Model to Optimize Agricultural Land Use1

Abstract: A principal contributor to soil erosion and nonpoint source pollution, agricultural activities have a major influence on the environmental quality of a watershed. Impact of agricultural activities on the quality of water resources can be minimized by implementing suitable agriculture land‐use types. Currently, land uses are designed (location, type, and operational schedule) based on field study results, and do not involve a science‐based approach to ensure their efficiency under particular regional, climatic, geological, and economical conditions. At present, there is a real need for new methodologies that can optimize the selection, design, and operation of agricultural land uses at the watershed scale by taking into account environmental, technical, and economical considerations, based on realistic simulations of watershed response. In this respect, the present study proposes a new approach, which integrates computational modeling of watershed processes, fluvial processes in the drainage network, and modern heuristic optimization techniques to design cost effective land‐use plans. The watershed model AnnAGNPS and the channel network model CCHE1D are linked together to simulate the sediment and pollutant transport processes. Based on the computational results, a multi‐objective function is set up to minimize soil losses, nutrient yields, and total associated costs, while the production profits from agriculture are maximized. The selected iterative optimization algorithm uses adaptive Tabu Search heuristic to flip (switching from one alternative to another) land‐change variables. USDA’s Goodwin Creek experimental watershed, located in Northern Mississippi, is used to demonstrate the capabilities of the proposed approach. The results show that the optimized land‐use design with BMPs using an integrated approach at the watershed level can provide efficient and cost‐effective conservation of the environmental quality by taking into account both productivity and profitability.     

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.     

Evaluation of the MIKE SHE Model for Application in the Loess Plateau,China1

Abstract: Quantifying the hydrologic responses to land use/land cover change and climate variability is essential for integrated sustainable watershed management in water limited regions such as the Loess Plateau in Northwestern China where an adaptive watershed management approach is being implemented. Traditional empirical modeling approach to quantifying the accumulated hydrologic effects of watershed management is limited due to its complex nature of soil and water conservation practices (e.g., biological, structural, and agricultural measures) in the region. Therefore, the objective of this study was to evaluate the ability of the distributed hydrologic model, MIKE SHE to simulate basin runoff. Streamflow data measured from an overland flow‐dominant watershed (12 km²) in northwestern China were used for model evaluation. Model calibration and validation suggested that the model could capture the dominant runoff process of the small watershed. We found that the physically based model required calibration at appropriate scales and estimated model parameters were influenced by both temporal and spatial scales of input data. We concluded that the model was useful for understanding the rainfall‐runoff mechanisms. However, more measured data with higher temporal resolution are needed to further test the model for regional applications.     

The Watershed Flow and Allocation Model: An NHDPlus‐Based Watershed Modeling Approach for Multiple Scales and Conditions

The Watershed Flow and Allocation model (WaterFALL^®) provides segment‐specific, daily streamflow at both gaged and ungaged locations to generate the hydrologic foundation for a variety of water resources management applications. The model is designed to apply across the spatially explicit and enhanced National Hydrography Dataset (NHDPlus) stream and catchment network. To facilitate modeling at the NHDPlus catchment scale, we use an intermediate‐level rainfall‐runoff model rather than a complex process‐based model. The hydrologic model within WaterFALL simulates rainfall‐runoff processes for each catchment within a watershed and routes streamflow between catchments, while accounting for withdrawals, discharges, and onstream reservoirs within the network. The model is therefore distributed among each NHDPlus catchment within the larger selected watershed. Input parameters including climate, land use, soils, and water withdrawals and discharges are georeferenced to each catchment. The WaterFALL system includes a centralized database and server‐based environment for storing all model code, input parameters, and results in a single instance for all simulations allowing for rapid comparison between multiple scenarios. We demonstrate and validate WaterFALL within North Carolina at a variety of scales using observed streamflows to inform quantitative and qualitative measures, including hydrologic flow metrics relevant to the study of ecological flow management decisions.     

Establishing aquatic restoration priorities using a watershed approach

Since the passage of the Clean Water Act in 1972, the United States has made great strides to reduce the threats to its rivers, lakes, and wetlands from pollution. However, despite our obvious successes, nearly half of the nation's surface water resources remain incapable of supporting basic aquatic values or maintaining water quality adequate for recreational swimming. The Clean Water Act established a significant federal presence in water quality regulation by controlling point and non-point sources of pollution. Point-sources of pollution were the major emphasis of the Act, but Section 208 specifically addressed non-point sources of pollution and designated silviculture and livestock grazing as sources of non-point pollution. Non-point source pollutants include runoff from agriculture, municipalities, timber harvesting, mining, and livestock grazing. Non-point source pollution now accounts for more than half of the United States water quality impairments. To successfully improve water quality, restoration practitioners must start with an understanding of what ecosystem processes are operating in the watershed and how they have been affected by outside variables. A watershed-based analysis template developed in the Pacific Northwest can be a valuable aid in developing that level of understanding. The watershed analysis technique identifies four ecosystem scales useful to identify stream restoration priorities: region, basin, watershed, and site. The watershed analysis technique is based on a set of technically rigorous and defensible procedures designed to provide information on what processes are active at the watershed scale, how those processes are distributed in time and space. They help describe what the current upland and riparian conditions of the watershed are and how these conditions in turn influence aquatic habitat and other beneficial uses. The analysis is organized as a set of six steps that direct an interdisciplinary team of specialists to examine the biotic and abiotic processes influencing aquatic habitat and species abundance. This process helps develop an understanding of the watershed within the context of the larger ecosystem. The understanding gained can then be used to identify and prioritize aquatic restoration activities at the appropriate temporal and spatial scale. The watershed approach prevents relying solely on site-level information, a common problem with historic restoration efforts. When the watershed analysis process was used in the Whitefish Mountains of northwest Montana, natural resource professionals were able to determine the dominant habitat forming processes important for native fishes and use that information to prioritize, plan, and implement the appropriate restoration activities at the watershed scale. Despite considerable investments of time and resources needed to complete an analysis at the watershed scale, the results can prevent the misdiagnosis of aquatic problems and help ensure that the objectives of aquatic restoration will be met.     

HYSTAR Sediment Model: Distributed Two‐Dimensional Simulation of Watershed Erosion and Sediment Transport Using Time‐Area Routing

An erosion and sediment transport component incorporated in the HYdrology Simulation using Time‐ARea method (HYSTAR) upland watershed model provides grid‐based prediction of erosion, transport and deposition of sediment in a dynamic, continuous, and fully distributed framework. The model represents the spatiotemporally varied flow in sediment transport simulation by coupling the time‐area routing method and sediment transport capacity approach within a grid‐based spatial data model. This avoids the common, and simplistic, approach of using the Universal Soil Loss Equation (USLE) to estimate erosion rates with a delivery ratio to relate gross soil erosion to sediment yield of a watershed, while enabling us to simulate two‐dimensional sediment transport processes without the complexity of numerical solution of the partial differential governing equations. In using the time‐area method for routing sediment, the model offers a novel alternative to watershed‐scale sediment transport simulation that provides detailed spatial representation. In predicting four‐year sediment hydrographs of a watershed in Virginia, the model provided good performance with R² of 0.82 and 0.78 and relative error of ?35% and 11% using the Yalin and Yang's sediment transport capacity equations, respectively. Prediction of spatiotemporal variation in sediment transport processes was evaluated using maps of sediment transport rates, concentrations, and erosion and deposition mass, which compare well with expected behavior of flow hydraulics and sediment transport processes.     

Baseflow nitrate in relation to stream order and agricultural land use

Management of agricultural nonpoint-source pollution continues to be a challenge because of spatial and temporal variability. Using stream order as an index, we explored the distribution of nitrate concentration and load along the stream network of a large agricultural watershed in Pennsylvania-the East Mahantango Creek Watershed and two of its sub-watersheds. To understand nitrate concentration variation in the stream water contributed from ground water, this study focused on baseflow. Impacts of agricultural land use area on baseflow nitrate in the stream network were investigated. Nitrate concentration showed a general decreasing trend with increasing stream order based on stream order averaged values; however, considerable spatial and temporal variability existed within each snapshot sampling. Nitrate loads increased with stream order in a power function because of the dominant effect of stream flow rate over the nitrate concentration. Within delineated sub-watersheds based on stream orders, positive linear functions were found between agricultural land use area percentage and the baseflow nitrate concentration and between agricultural drainage area and the nitrate load. The slope of the positive linear regression between the baseflow nitrate concentration and percent agricultural land area seems to be a valuable indicator of a watershed's water quality as influenced by agricultural practices, watershed size, and specific physiographic setting. Stream order seems to integrate, to a certain degree, the source and transport aspects of nonpoint-source pollution on a yearly averaged basis and thus might provide a quick estimate of the overall trend in baseflow nitrate concentration and load distribution along complex stream networks in agricultural watersheds.     

THREE‐DIMENSIONAL KRIGING ANALYSIS OF SEDIMENT MERCURY DISTRIBUTION: A CASE STUDY1

ABSTRACT: This study presents an application of a three‐dimensional kriging technique to estimate spatial distribution of total mercury (Hg) in the Cedar‐Ortega Rivers watershed located in the lower St. Johns River basin, Florida. The kriging procedures, including preliminary data analysis, structural data analysis and the log kriging estimation, are presented. Results show watershed wide Hg contamination of river sediment to a depth of 1.0 m. A three‐dimensional plot of Hg against the Florida Sediment Assessment Guidelines (i.e., the probable effect level or PEL) demonstrates that the Cedar River is more contaminated with Hg than the rest of the watershed. The maximum sediment depth with Hg concentrations above PEL value (0.696 mg/kg) is 1.5 m. Hg concentrations at or above this level could pose a significant hazard to aquatic organisms. Analysis of the spatial distribution of Hg in the watershed finds multiple input sources. This study suggests that there is a need to identify the major sources of Hg in the watershed, and to determine the pathways that allow Hg to enter the river.     

A Gis‐Based Approach to Watershed Classification for Nebraska Reservoirs1

Abstract: The U.S. Environmental Protection Agency is charged with establishing standards and criteria for assessing lake water quality. It is, however, increasingly evident that a single set of national water quality standards that do not take into account regional hydrogeologic and ecological differences will not be viable as lakes clearly have different inherent capacities to meet such standards. We demonstrate a GIS‐based watershed classification strategy for identifying groups of Nebraska reservoirs that have similar potential capacity to attain a certain level of water quality standard. A preliminary cluster analysis of 78 reservoirs was performed to determine the potential number of Nebraska reservoir groups. Subsequently, a Classification Trees method was used to refine number of classes, describe the structure of reservoir watershed classes, and to develop a predictive model that relates watershed conditions to reservoir classes. Results suggest that Nebraska reservoirs can be represented by nine classes and that soil organic matter content in the watershed is the most important single variable for segregating the reservoirs. The cross‐validation prediction error rate of the Classification Tree model was 26.3%. Because all geospatial data used in this work are available nationally, the method could be adopted throughout the U.S. Hence, this GIS‐based watershed classification approach could provide water resources managers an effective decision‐support tool in managing reservoir water quality.     

Development and Application of an Agricultural Intensity Index to Invertebrate and Algal Metrics from Streams at Two Scales

Research was conducted at 28‐30 sites within eight study areas across the United States along a gradient of nutrient enrichment/agricultural land use between 2003 and 2007. Objectives were to test the application of an agricultural intensity index (AG‐Index) and compare among various invertebrate and algal metrics to determine indicators of nutrient enrichment nationally and within three regions. The agricultural index was based on total nitrogen and phosphorus input to the watershed, percent watershed agriculture, and percent riparian agriculture. Among data sources, agriculture within riparian zone showed significant differences among values generated from remote sensing or from higher resolution orthophotography; median values dropped significantly when estimated by orthophotography. Percent agriculture in the watershed consistently had lower correlations to invertebrate and algal metrics than the developed AG‐Index across all regions. Percent agriculture showed fewer pairwise comparisons that were significant than the same comparisons using the AG‐Index. Highest correlations to the AG‐Index regionally were ?0.75 for Ephemeroptera, Plecoptera, and Trichoptera richness (EPTR) and ?0.70 for algae Observed/Expected (O/E), nationally the highest was ?0.43 for EPTR vs. total nitrogen and ?0.62 for algae O/E vs. AG‐Index. Results suggest that analysis of metrics at national scale can often detect large differences in disturbance, but more detail and specificity is obtained by analyzing data at regional scales.     

Instream Large Wood: Denitrification Hotspots with Low N2O Production

We examined the effect of instream large wood on denitrification capacity in two contrasting, lower order streams — one that drains an agricultural watershed with no riparian forest and minimal stores of instream large wood and another that drains a forested watershed with an extensive riparian forest and abundant instream large wood. We incubated two types of wood substrates (fresh wood blocks and extant streambed wood) and an artificial stone substrate for nine weeks in each stream. After in situ incubation, we collected the substrates and their attached biofilms and established laboratory‐based mesocosm assays with stream water amended with ¹⁵N‐labeled nitrate‐N. Wood substrates at the forested site had significantly higher denitrification than wood substrates from the agricultural site and artificial stone substrates from either site. Nitrate‐N removal rates were markedly higher on woody substrates compared to artificial stones at both sites. Nitrate‐N removal rates were significantly correlated with biofilm biomass. Denitrification capacity accounted for only a portion of nitrate‐N removal observed within the mesocosms in both the wood controls and instream substrates. N₂ accounted for 99.7% of total denitrification. Restoration practices that generate large wood in streams should be encouraged for N removal and do not appear to generate high risks of instream N₂O generation.     

Structured Decision Making to Meet a National Water Quality Mandate

Water quality criteria are necessary to ensure protection of ecological and human health conditions, but compliance can require complex decisions. We use structured decision making to consider multiple stakeholder objectives in a water quality management process, with a case study in the Three Bays watershed on Cape Cod, Massachusetts. We set a goal to meet or exceed a nitrogen load reduction target for the watershed and four key objectives: minimizing economic costs of implementing management actions, minimizing the complexity of permitting management actions, maximizing stakeholder acceptability of the management actions, and maximizing the provision of ecosystem services (recreational opportunity, erosion and flood control, socio‐cultural amenity). We used multi‐objective optimization and sensitivity analysis to generate many possible solutions that implement different combinations of nitrogen‐removing management actions and reflect tradeoffs between the objectives. Results show technological advances in controlling household nitrogen sources could provide lower cost solutions and positive impacts to ecosystem services. Although this approach is demonstrated with Cape Cod data, the decision‐making process is not specific to any watershed and could be easily applied elsewhere.     

Integrating Modeling and Field Experiments to Evaluate Impacts of Vegetative Practices on Ponderosa Pine Watersheds1

Abstract: Determining watershed response to vegetation treatment has been the subject of numerous hydrologic studies over the years. However, generalizing the information obtained from traditional paired‐watershed studies to other watersheds in a region is problematic because of the empirical nature of such studies and the context dependence of hydrologic responses. This paper addresses the issue of generalizing hydrologic information through integration of process‐based modeling and field observations from small‐scale watershed experiments. To this end, the results from application of a process‐based model were compared with the results from small‐scale watershed experiments in ponderosa pine forests of Arizona. The model simulated treatment impacts reasonably well when compared to the traditional paired‐watershed approach. However, the model tended to overestimate water yields during periods of low flow, and there was a significant difference between the two approaches in the estimation of treatment impacts during the first four years following treatment. The results indicate that the lumped‐parameter modeling approach used here may be limited in its ability to detect small changes, and tends to overestimate changes that occur immediately following treatment. It is concluded that watershed experiments can be highly informative due to their direct examination of cause‐effect relationships, while process‐based models are useful for their processing power and focus on functional relationships. The integrated use of both watershed experiments and process‐based models provides a way to generalize hydrologic information, illuminate the processes behind landscape treatment effects, and to generate and test hypotheses.     

SPATIALLY EXPLICIT HYDROLOGIC MODELING OF LAND USE CHANGE1

ABSTRACT: Using a Geographic Information System (GIS), a method is presented to develop a spatially explicit time series of land use in an urbanizing watershed. The method is prefaced on the existence of independent observations of land use at different times and data that describes the spatial‐temporal land use transition characteristics of the watershed between these two points in time. A method is then presented to generalize the TR‐55 graphical method, a common lumped hydrologic model for estimating peak discharge, for use in a spatially explicit scheme. This scheme predicts peak discharge throughout a watershed, rather than at a single selected watershed outlet. Coupling these two methods allows the engineer to model both the temporal and spatial evolution of peak discharge for the watershed. An illustrative watershed in a suburban area of Washington, DC is selected to demonstrate the methods. The model results from these analyses are presented graphically to highlight the complex features in peak discharge behavior that exist both spatially, as a function of position within the watershed drainage network, and temporally, as the watershed undergoes urbanization. These features are not commonly noted in most hydrologic analyses but are captured in these analyses because of the high spatial and temporal resolution of the methods presented. The physical implications of the modeled results are discussed in the context of the information content of a stream gauge located at the overall outlet of the illustrative watershed. This work shows that the common practice of transposition of gauge information to locations internal to the watershed would neglect internal variability in peak discharge behavior, and could potentially lead to the determination of inappropriate design discharges.     

Factors influencing scaling up of agroforestry-based spatial land-use integration for soil fertility management in Arsamma Watershed,Southwestern Ethiopian Highlands

This study was set to examine factors influencing agroforestry upscaling, inter-plot natural fertilizer transfer and inter-plot income flow in Arsamma watershed. A semi-structured questionnaire was used to gather necessary information. Contingency table, chi-square, Phi and Cramer's V were used to analyze the data. Access to seedlings was the most important determinant of agroforestry upscaling; and farmers' production orientation, farm size and wealth status ranked, respectively, second, third and fourth. Inter-plot natural fertilizer transfer was primarily influenced by participation in agricultural extension. Wealth status and livestock size ranked, respectively, second and third in influencing inter-plot natural fertilizer transfer. The study indicated a geographic concept of spatial land-use integration for soil fertility management and key factors influencing agroforestry-based land-use integration. Agroforestry-centered diversified small-scale agricultural commercialization, tree-crop-livestock integration, agricultural extension services promotion and multi-purpose tree species supply are the way out to upscale agroforestry and agroforestry-based spatial land-use integration.     

Identification of Priority Areas for Integrated Management of Semiarid Watersheds in the Ecuadorian Andes

Integrated watershed management (IWM) is a priority, especially in semiarid regions that are concurrently affected by population growth, land use change, soil erosion, and poor governance. In developing countries, IWM is often done without any support tool, scientific data, or deep knowledge of territory characteristics. The aim of this study was to present a case study to apply a decision support tool to prioritize areas for territory management. A simple, quantitative multi‐criteria analysis was applied in a semiarid basin of the Ecuadorian Andes to identify the zones of greatest concern for implementation of resource conservation and management practices at a local and regional scale. In addition to describing the current state of the conditions of this basin, our results suggest scenarios of change in relation to official population projections based on spatial analysis of land use change. Analysis resulted in a scattered distribution of priority values within the watershed, so a hierarchical rule was incorporated to define priorities at the subwatershed (SW) scale. Our analysis identified four SW of very high priority and urgent need to implement management practices. Based on projections of future change due to population growth and land cover change, the number of subbasins that require more attention was doubled. Finally, this study includes zones for management or conservation of the land, according to the Sustainable Development Goals.     

Modeling Agricultural Nonpoint Source Pollution Using a Geographic Information System Approach

Agricultural non-point source (NPS) pollution, primarily sediment and nutrients, is the leading source of water-quality impacts to surface waters in North America. The overall goal of this study was to develop geographic information system (GIS) protocols to facilitate the spatial and temporal modeling of changes in soils, hydrology, and land-cover change at the watershed scale. In the first part of this article, we describe the use of GIS to spatially integrate watershed scale data on soil erodibility, land use, and runoff for the assessment of potential source areas within an intensively agricultural watershed. The agricultural non-point source pollution (AGNPS) model was used in the Muddy Creek, Ontario, watershed to evaluate the effectiveness of management strategies in decreasing sediment and nutrient [phosphorus (P)] pollution. This analysis was accompanied by the measurement of water-quality parameters (dissolved oxygen, pH, hardness, alkalinity, and turbidity) as well as sediment and P loadings to the creek. Practices aimed at increasing year-round soil cover would be most effective in decreasing sediment and P losses in this watershed. In the second part of this article, we describe a method for characterizing land-cover change in a dynamic urban fringe watershed. The GIS method we developed for the Blackberry Creek, Illinois, watershed will allow us to better account for temporal changes in land use, specifically corn and soybean cover, on an annual basis and to improve on the modeling of watershed processes shown for the Muddy Creek watershed. Our model can be used at different levels of planning with minimal data preprocessing, easily accessible data, and adjustable output scales.     

NONPOINT SOURCE POLLUTION POTENTIAL IN AN AGRICULTURAL WATERSHED IN NORTHWESTERN PENNSYLVANIA1

ABSTRACT: A 155,947 ha portion of the Shenango River watershed in western Pennsylvania was evaluated as to the potential impact of agriculture drainage on water quality. Approximately a third of the area is being used as either cropland or pasture with approximately an equal percentage in forest lands. Eleven subwatersheds were evaluated as to their potential for nonpoint source pollution according to the criteria established by the Pennsylvania Department of Environmental Resources for the Chesapeake Bay Pollution Abatement Program. The individual components and overall rating for each subwatershed were then evaluated as to their correlation with four water quality variables based on 104 samples collected at 26 sampling stations throughout the watershed. There was a significant correlation between the overall rating factor for each subwatershed and each of the four water quality variables. In general, the watershed delivery factor, animal nutrient factor, and management factors were correlated with fecal coliform and phosphorus in the receiving streams, whereas the ground water delivery factor appeared to be more important in determining nitrate concentrations in these streams. These results indicate that manure and nutrient management, along with the exclusion of livestock from streams and the enhancement and/or replacement of riparian wetlands, are important approaches in reducing agricultural impacts in fresh water ecosystems.     

AN EXPANDING ROLE FOR SUBARCTIC WATERSHED RESEARCH1

ABSTRACT: It is increasingly recognized that natural resources research should in many cases be broadened in scope and oriented toward more general “environmental” problems. Locales with a history of “watershed” research can be eminently suited for development of comprehensive, environmental research programs. This is recognized in many research efforts of the International Biological Program (IBP), where watershed research sites have been successfully utilized for intensive investigations of process and function of selected ecosystems or ecosystem components. In the North American Subarctic there is almost no history of “watershed” studies. Basic data on hydrometeorologic parameters such as precipitation amounts and areal and seasonal distribution of runoff are scarce; the data framework within which environmental understanding can be structured is exceedingly sketchy. Opportunity exists in the discontinuous-permafrost settings of central Alaska to begin rectifying this situation. A basic program of multi-agency, multi-discipline research and data acquisition for the most significant hydrologic subregions is being developed, based around several existing environmental research areas (chiefly the Bonanza Creek Experimental Forest, the Caribou-Poker Creeks Research Watershed, the Wickersham Dome Fire Study Area, and a series of outlying sites).