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.     

MINIMIZING THE IMPACT OF URBANIZATION ON LONG TERM RUNOFF1

Increasing concern about the problems caused by urban sprawl has encouraged development and implementation of smart growth approaches to land use management. One of the goals of smart growth is water resources protection, in particular minimizing the runoff impact of urbanization. To investigate the magnitude of the potential benefits of land use planning for water resources protection, possible runoff impacts of historical and projected urbanization were estimated for two watersheds in Indiana and Michigan using a long term hydrological impact analysis model. An optimization component allowed selection of land use change placements that minimize runoff increase. Optimizing land use change placement would have reduced runoff increase by as much as 4.9 percent from 1973 to 1997 in the Indiana study watershed. For nonsprawl and sprawl scenarios in the Michigan watershed for 1978 to 2040, optimizing land use change placement would have reduced runoff increase by 12.3 percent and 20.5 percent, respectively. The work presented here illustrates both an approach to assessing the magnitude of the impact of smart growth and the significant potential scale of smart growth in moderating runoff changes that result from urbanization. The results of this study have significant implications for urban planning.     

The Impact of Future Land Use Scenarios on Runoff Volumes in the Muskegon River Watershed

In this article we compared the response of surface water runoff to a storm event for different rates of urbanization, reforestation and riparian buffer setbacks across forty subwatersheds of the Muskegon River Watershed located in Michigan, USA. We also made these comparisons for several forecasted and one historical land use scenarios (over 140 years). Future land use scenarios to 2040 for forest regrowth, urbanization rates and stream setbacks were developed using the Land Transformation Model (LTM). Historical land use information, from 1900 at 5-year time step intervals, was created using a Backcast land use change model configured using artificial neural network and driven by agriculture and housing census information. We show that (1) controlling the rate of development is the most effective policy option to reduce runoff; (2) establishing setbacks along the mainstem are not as effective as controlling urban growth; (3) reforestation can abate some of the runoff effects from urban growth but not all; (4) land use patterns of the 1970s produced the least amount of runoff in most cases in the Muskegon River Watershed when compared to land use maps from 1900 to 2040; and, (5) future land use patterns here not always lead to increased (worse) runoff than the past. We found that while ten of the subwatersheds contained futures that were worse than any past land use configuration, twenty-five (62.5%) of the subwatersheds produced the greatest amount of runoff in 1900, shortly after the entire watershed was clear-cut. One third (14/40) of the subwatersheds contained the minimum amount of runoff in the 1960s and 1970s, a period when forest amounts were greatest and urban amounts relatively small.     

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.     

ANALYSIS AND PREDICTION OF SURFACE RUNOFF IN AN URBANIZING WATERSHED USING SATELLITE IMAGERY1

ABSTRACT: This paper demonstrates how satellite image data [e.g., from Landsat 5 Thematic Mapper (TM)], in conjunction with an urban growth model and simple runoff calculations, can be used to estimate future surface runoff and, by implication, water quality within a watershed. To illustrate the method, predictions of land use change and surface runoff are shown for Spring Creek Watershed, a medium sized urbanizing watershed in Central Pennsylvania. Land cover classifications for this watershed were created from images for summertime 1986 and 1996 and subsequently used as input to the Clarke urban growth model, called SLEUTH, to predict land use changes to the year 2025. Simulations with this model show a progressive growth in the percentage of urban pixels and in impervious surface area in the watershed but also an increase in woodland, primarily in previously clear‐cut areas. Given that woodland area will continue to increase in area, surface runoff into Spring Creek is predicted to remain only slightly above present level. However, should the woodland amount fail to increase, surface runoff is then predicted to increase more significantly during the next 25 years. Finally, the concept of urban sprawl is addressed within the context of predicted increases in urbanization by relating the implied increase in impervious surface area to population density within the watershed.     

MODELING THE EFFECTS OF URBANIZATION ON BASIN WATER YIELD AND RESERVOIR SEDIMENTATION1

ABSTRACT: The dam impounding White Rock Lake was completed in 1910 to provide water for the City of Dallas. Since then, land use on the watershed has changed from entirely rural to over 77 percent urban. A model called SWRRB (Simulator for Water Resources in Rural Basins) was utilized to determine the effect of urbanization on water and sediment entering the lake. The simulation results show that, if urbanization had not occurred, then the annual surface runoff would be 135 mm rather than 151 mm and the annual sediment yield would be 4.4 t/ha rather than 4.1 t/ha. Also, the effect of urbanization on delivery ratios was shown and a positive linear correlation was found. Finally, the weather generator in SWRRB was utilized to estimate the loss of reservoir capacity until 2050 for three different land use management scenarios.     

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.     

Weighting Nitrogen and Phosphorus Pixel Pollutant Loads to Represent Runoff and Buffering Likelihoods

Watershed models often estimate annual nitrogen (N) or phosphorus (P) pollutant loads in rural areas with export coefficient (EC) (kg/ha/yr) values based on land cover, and in urban areas as the product of spatially uniform event mean concentration (EMC) (mg/L) values and runoff volume. Actual N and P nonpoint source (NPS) pollutant loading has more spatial complexity due to watershed variation in runoff likelihood and buffering likelihood along surface and subsurface pathways, which can be represented in a contributing area dispersal area (CADA) NPS model. This research develops a CADA NPS model to simulate how watershed properties of elevation, land cover, and soils upslope and downslope of each watershed pixel influence nutrient loading. The model uses both surface and subsurface runoff indices (RI), and surface and subsurface buffer indices (BI), to quantify the runoff and buffering likelihood for each watershed pixel, and generate maps of weighted EC and EMC values that identify NPS pollutant loading hotspots. The research illustrates how CADA NPS model maps and pixel loading values are sensitive to the spatial resolution and accuracy of elevation and land cover data, and model predictions can represent the lower and upper bounds of NPS loading. The model provides managers with a tool to rapidly visualize, rank, and investigate likely areas of high nutrient export.     

GIS BASED LONG TERM HYDROLOGIC IMPACT EVALUATION FOR WATERSHED URBANIZATION1

ABSTRACT: To adequately manage impacts of ongoing or future land use changes in a watershed, the magnitude of their hydrologic impacts needs to be assessed. A grid based daily streamflow model was calibrated with two years of observed streamflow data, using time periods when land use data are available and verified by comparison of model predictions with observed streamflow data. Streamflow data were separated into direct runoff and baseflow to estimate the impacts of urbanization on each hydrologic component. Analysis of the ratio between direct runoff and total runoff from 30 years of simulation results and the change in these ratios with urbanization shows that estimated annual direct runoff increased from 49.2 percent (1973) to 63.1 percent (1984) and 65.0 percent (1991), indicating the effects of urbanization are greater on direct runoff than on total runoff. The direct runoff ratio also varies with annual rainfall, with dry year ratios larger than those for wet years. This suggests that the impact of urbanization on areas that are sensitive to runoff ratios, such as stream ecosystems, might be more serious during drier years than in wetter years in terms of water quality and water yield. This indicates that sustainable base‐flow is important to maintaining sound stream ecosystems.     

A Water Quality Model for Regional Stream Assessment and Conservation Strategy Development

Non-point-source (NPS) pollution remains the primary source of stream impairment in the United States. Many problems such as eutrophication, sedimentation, and hypoxia are linked with NPS pollution which reduces the water quality for aquatic and terrestrial organisms. Increasingly, NPS pollution models have been used for landscape-scale pollution assessment and conservation strategy development. Our modeling approach functions at a scale between simple landscape-level assessments and complex, data-intensive modeling by providing a rapid, landscape-scale geographic information system (GIS) model with minimal data requirements and widespread applicability. Our model relies on curve numbers, literature-derived pollution concentrations, and land status to evaluate total phosphorus (TP), total nitrogen (TN), and suspended solids (SS) at the reach scale. Model testing in the Chesapeake Bay watershed indicated that predicted distributions of water quality classes were realistic at the reach scale, but precise estimates of pollution concentrations at the local scale can have errors. Application of our model in the tributary watersheds along Lake Ontario suggested that it is useful to managers in watershed planning by rapidly providing important information about NPS pollution conditions in areas where large data gaps exist, comparisons among stream reaches across numerous watersheds are required, or regional assessments are sought.     

SUSCEPTIBILITY OF INDIANA WATERSHEDS TO HERBICIDE CONTAMINATION1

ABSTRACT: An index of watershed susceptibility to surface water contamination by herbicides could be used to improve source water assessments for public drinking water supplies, prioritize watershed restoration projects, and direct funding and educational efforts to areas where the greatest environmental benefit can be realized. The goal of this study is to use streamflow and herbicide concentration data to develop and evaluate a method for estimating comparative watershed susceptibility to herbicide loss. United States Geological Survey (USGS) concentration data for five relatively water soluble herbicides (alachlor, atrazine, cyanazine, metolachlor, and simazine) were analyzed for 16 Indiana watersheds. Correlation was assessed between observed herbicide losses and: (1) a herbicide runoff index using GIS‐based land use, soil type, SCS runoff curve number, tillage practice, herbicide use estimates, and combinations of these factors; and (2) predicted herbicide losses from a non‐point source pollution model (NAPRA‐Web, an Internet‐based interface for GLEAMS). The highest adjusted R²value was found between herbicide concentration and the runoff curve number alone, ranging from 0.25 to 0.56. Predictions from the simulation model showed a poorer correlation with observed herbicide loss. This indicates potential for using the runoff curve number as a simple herbicide contamination susceptibility index.     

Modeling the Effects of Land Use Change on the Water Temperature in Unregulated Urban Streams

Streams, in their natural state, are typically diverse and biologically productive environments. Streams subject to urbanization often experience degradation brought about by the cumulative effects of flow alteration, unsanitary discharge and channelization. One of the water quality parameters affected by urbanization is stream temperature. This study offers a model for predicting the impact of land use change on the temperature of non-regulated streams during extreme events. A stream temperature model was created by considering the gains and losses of thermal energy resulting from radiation, convection, conduction, evaporation and advection. A sensitivity analysis showed that out of 14 variables, shade/transmissivity of riparian vegetation, groundwater discharge, and stream width had the greatest influence on stream temperature. These same three variables are highly influenced by land use. Individual component models were developed to predict how urbanization changes stream width and baseflow discharge. Using 3-D computer modeling, a model was also developed to illustrate the effects of altering the extent and composition of riparian vegetation on streams with different orientations. By modeling these three variables as a function of urbanization, the results became inputs into the stream temperature model. The critical urban stream temperature model (CrUSTe), an aggregation of these four models, allows the prediction of stream temperature change as a result of amount, type and location of urbanization within a watershed. It has the potential to become a valuable tool for environmental managers.     

Long-Term Impacts of Land-Use Change on Non-Point Source Pollutant Loads for the St. Louis Metropolitan Area, USA

A land-use-change simulation model (LEAM) and a non-point-source (NPS) water quality model (L-THIA) were closely coupled as LEAMwq in order to determine the long-term implications of various degree of urbanization on NPS total nitrogen (TN), total suspended particles (TSP), and total phosphorus (TP) loads. A future land-use projection in the St. Louis metropolitan area from 2005 to 2030 using three economic growth scenarios (base, low, and high) and a long-term precipitation dataset were used to predict the mean annual surface runoff and mean annual NPS pollutant loads in the region. Results show mean annual TN increases of 0.21%, 0.13%, and 0.14% by 2030 compared to 2000 under the base, high, and low scenarios, respectively. TSP and TP showed similar trends with different magnitudes. Corresponding changes in annual mean surface runoff were shown to be lower than expected, which might be attributed to the small-scale conversion pattern of land uses. In the most dramatic change (high growth) scenario, the runoff would increase across time but at varying rates, and temporal pollutant loads would result in a more complicated pattern than in the other scenarios. This is attributed to the complex interactions between event mean concentrations of pollutants and the magnitude of changes in land-use acreages. By integrating L-THIA with LEAM, LEAMwq was found to be a useful planning tool to illustrate in a quick and simple manner how future water quality is connected to decision-making on future land-use change.     

IMPACTS OF CALIFORNIA'S CLIMATIC REGIMES AND COASTAL LAND USE CHANGE ON STREAMFLOW CHARACTERISTICS1

ABSTRACT: To investigate the impacts of urbanization and climatic fluctuations on stream flow magnitude and variability in a Mediterranean climate, the HEC‐HMS rainfall/runoff model is used to simulate stream flow for a 14‐year period (October 1, 1988, to September 30, 2002) in the Atascadero Creek watershed located along the southern coast of California for 1929, 1998, and 2050 (estimated) land use conditions (8, 38 and 52 percent urban, respectively). The 14‐year period experienced a range of climatic conditions caused mainly by El Nino‐Southern Oscillation variations. A geographic information system is used to delineate the watershed and parameterize the model, which is calibrated using data from two stream flow and eight rainfall gauges. Urbanization is shown to increase peak discharges and runoff volume while decreasing stream flow variability. In all cases, the annual and 14‐year distributions of stream flow are shown to be highly skewed, with the annual maximum 24 hours of discharge accounting for 22 to 52 percent of the annual runoff and the maximum ten days of discharge from an average El Nino year producing 10 to 15 percent of the total 14‐year discharge. For the entire period of urbanization (1929 to 2050), the average increase in annual maximum discharges and runoff was 45 m³/s (300 percent) and 15 cm (350 percent), respectively. Additionally, the projected increase in urbanization from 1998 to 2050 is half the increase from 1929 to 1998; however, increases in runoff (22 m³/s and 7 cm) are similar for both scenarios because of the region's spatial development pattern.     

Hydrological Responses to Climate and Land‐Use Changes along the North American East Coast: A 110‐Year Historical Reconstruction

The North American east coast (NAEC) region experienced significant climate and land‐use changes in the past century. To explore how these changes have affected land water cycling, the Dynamic Land Ecosystem Model (DLEM 2.0) was used to investigate the spatial and temporal variability of runoff and river discharge during 1901‐2010 in the study area. Annual runoff over the NAEC was 420 ± 61 mm/yr (average ± standard deviation). Runoff increased in parts of the northern NAEC but decreased in some areas of the southern NAEC. Annual freshwater discharge from the study area was 378 ± 61 km³/yr (average ± standard deviation). Factorial simulation experiments suggested that climate change and variability explained 97.5% of the interannual variability of runoff and also resulted in the opposite changes in runoff in northern and southern regions of the NAEC. Land‐use change reduced runoff by 5‐22 mm/yr from 1931 to 2010, but the impacts were divergent over the Piedmont region and Coastal Plain areas of the southern NAEC. Land‐use change impacts were more significant at local and watershed spatial scales rather than at regional scales. Different responses of runoff to changing climate and land‐use should be noted in future water resource management. Hydrological impacts of afforestation and deforestation as well as urbanization should also be noted by land‐use policy makers.     

Combined monitoring and modeling indicate the most effective agricultural best management practices

Although water quality problems associated with agricultural nonpoint source (NPS) pollution have prompted the rapid and widespread adoption of a variety of so called "best management practices" (BMPs), there have been few realistic efforts to assess their combined effectiveness in reducing NPS pollution. This study used the Variable Source Loading Function (VSLF) model, a distributed watershed model, to simulate phosphorus (P) loading from an upstate New York dairy farm before and after the implementation of a suite of BMPs. With minimal calibration, the model calculates the dissolved P (DP) losses from impervious surfaces (e.g., barnyards), the plant/soil complex, field-applied manure, and loads associated with baseflow conditions. The simulated DP loads agreed well with measured loads for both the pre-BMP and post-BMP periods. More importantly, results showed that BMPs reduced DP loads by 35%, which is over half of the expected reduction if all manure was removed from the watershed, i.e., approximately 50% reduction. The model results indicate that had no BMPs been installed DP loads would be approximately 37% greater than observed at the watershed outlet. The most effective BMPs were those that disassociated pollutant loading areas from areas prone to generating runoff, i.e., hydrologically sensitive areas. By contrast, attempts to reduce P content in manure were somewhat less effective. This study demonstrates that a combination of distributed, mechanistic modeling and long-term monitoring provides better insights into the effectiveness of water quality protection efforts than either individually.     

WATER QUALITY IN AGRICULTURAL,URBAN, AND MIXED LAND USE WATERSHEDS1

ABSTRACT: Water quality and nonpoint source (NPS) pollution are important issues in many areas of the world, including the Inner Bluegrass Region of Kentucky where urban development is changing formerly rural watersheds into urban and mixed use watersheds. In watersheds where land use is mixed, the relative contributions of NPS pollution from rural and urban land uses can be difficult to separate. To better understand NPS pollution sources in mixed use watersheds, surface water samples were taken at three sites that varied in land use to examine the effect of land use on water quality. Within the group of three watersheds, one was predominately agriculture (Agricultural), one was predominately urban (Urban), and a third had relatively equal representation of both types of land uses (Mixed). Nitrogen (N), phosphorus (P), total suspended solids (TSS), turbidity, pH, temperature, and streamflow were measured for one year. Comparisons are made among watersheds for concentration and fluxes of water quality parameters. Nitrate and orthophosphate concentrations were found to be significantly higher in the Agricultural watershed. Total suspended solids, turbidity, temperature, and pH, were found to be generally higher in the Urban and Mixed watersheds. No differences were found for streamflow (per unit area), total phosphorus, and ammonium concentrations among watersheds. Fluxes of orthophosphate were greater in the Agricultural watershed that in the Urban watershed while fluxes of TSS were greater in the Mixed watershed when compared to the Agricultural watershed. Fluxes of nitrate, ammonium, and total phosphorus did not vary among watersheds. It is apparent from the data that Agricultural land uses are generally a greater source of nutrients than the Urban land uses while Urban land uses are generally a greater source of suspended sediment.     

A SPATIAL DECISION SUPPORT SYSTEM FOR ASSESSING AGRICULTURAL NONPOINT SOURCE POLLUTION1

ABSTRACT: A spatial decision support system (SDSS) was developed to assess agricultural nonpoint source (NPS) pollution using an NPS pollution model and geographic information systems (GIS). With minimal user interaction, the SDSS assists with extracting the input parameters for a distributed parameter NPS pollution model from user-supplied GIS base layers. Thus, significant amounts of time, labor, and expertise can be saved. Further, the SDSS assists with visualizing and analyzing the output of the NPS pollution simulations. Capabilities of the visualization component include displays of sediment, nutrient, and runoff movement from a watershed. The input and output interface techniques/algorithms used to develop the SDSS, along with an example application of the SDSS, are described.     

NONPARAMETRIC STATISTICAL METHODS IN URBAN HYDROLOGIC RESEARCH1

ABSTRACT. In urban hydrologic studies, it is often necessary to determine the effect of changes in urban land use patterns on such runoff characteristics as flood peaks and flow volumes. Nonparametric statistical methods have certain properties that make them a valuable tool for detecting hydrologic change caused by a treatment, such as urbanization, that changes watershed over a period of time. As many hydrologists do not have a working familiarity with nonparametric methods, a number of them are used for illustrative purposes to analyze the effect of urbanization on 24 years of annual flood peaks for a Louisville, Kentucky, watershed. In the example, urbanization was found to increase the central tendency, but not the dispersion of the peaks. Peak flows modeled by holding watershed parameters constant were also found to be increasing because of an upward trend in precipitation. By following the numerical examples in the paper and looking up test statistics in referenced sources, the reader can easily apply these methods to other situations.     

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.