POLLUTANT RUNOFF MODELS: SELECTION AND USE IN DECISION MAKING1

Models for pollutant runoff can be useful in water quality management planning if appropriately structured for the problem at hand. Accordingly, a “top-down” approach is proposed for the examination of extant pollutant runoff models. The approach consists of the identification of objectives and attributes that reflect the needs of planners and decision makers when these models are used for water quality management planning. Ideally, the attributes should concern the effect of model information on improved decision making and the cost of model application. Practical difficulties with the first attribute necessitates substitution of surrogate attributes reflecting model appropriateness, resolution, and uncertainty. Common pollutant runoff models, in particular export coefficients and hydrology-driven simulation models, are found to have serious weaknesses on some of the attribute scales. The “top-down” approach leads to a set of desirable pollutant runoff model attributes; alternate modeling techniques are thus examined in order to identify promising future directions for model development. The focus of this examination is phosphorus, due to its importance in the eutrophication of surface waters. Models for both sediment-attached and dissolved phosphorus are considered. Among the conclusions is the belief that the partial contributing area concept can yield an effective yet simple simulation despite the variable and complex nature of runoff.     

MODELING DEVELOPED COASTAL WATERSHEDS WITH THE AGRICULTURAL NON-POINT SOURCE MODEL1

ABSTRACT: Many coastal states are facing increasing urban growth along their coast lines. The growth has caused urban non-point source nitrogen runoff to be a major contributor to coastal and estuarine enrichment. Water resource managers are responsible for evaluating the impacts from point and non-point sources in developed watersheds and developing strategies to manage future growth. Non-point source models provide an effective approach to these management challenges. The Agricultural Non-Point Source Model (AGNPS) permits the incorporation of important spatial information (soils, landuse, topography, hydrology) in simulating surface hydrology and nitrogen non-point source runoff. The AGNPS model was adapted for developed coastal watersheds by deriving urban coefficients that reflect urban landuse classes and the amount of impervious surface area. Popperdam Creek watershed was used for model parameter development and model calibration. Four additional watersheds were simulated to validate the model. The model predictions of the peak flow and total nitrogen concentrations were close to the field measurements for the five sub-basins simulated. Measured peak flow varied by 30 fold among the sub-basins. The average simulated peak flow was within 14 percent of the average measured peak flow. Measured total nitrogen loads varied over an order of magnitude among the sub-basins yet error between the measured and simulated loads for a given sub-basin averaged 5 percent. The AGNPS model provided better estimates of nitrogen loads than widely used regression methods. The spatial distribution of important watershed characteristics influenced the impacts of urban landuse and projecting future residential expansion on runoff, sediment and nitrogen yields. The AGNPS model provides a useful tool to incorporate these characteristics, evaluate their importance, and evaluate fieldscale to watershed-scale urban impacts.     

Effects of calibration on L-THIA GIS runoff and pollutant estimation

Urbanization can result in alteration of a watershed's hydrologic response and water quality. To simulate hydrologic and water quality impacts of land use changes, the Long-Term Hydrologic Impact Assessment (L-THIA) system has been used. The L-THIA system estimates pollutant loading based on direct runoff quantity and land use based pollutant coefficients. The accurate estimation of direct runoff is important in assessing water quality impacts of land use changes. An automated program was developed to calibrate the L-THIA model using the millions of curve number (CN) combinations associated with land uses and hydrologic soil groups. L-THIA calibration for the Little Eagle Creek (LEC) watershed near Indianapolis, Indiana was performed using land use data for 1991 and daily rainfall data for six months of 1991 (January 1-June 30) to minimize errors associated with use of different temporal land use data and rainfall data. For the calibration period, the Nash-Sutcliffe coefficient was 0.60 for estimated and observed direct runoff. The calibrated CN values were used for validation of the model for the same year (July 1-December 31), and the Nash-Sutcliffe coefficient was 0.60 for estimated and observed direct runoff. The Nash-Sutcliffe coefficient was 0.52 for January 1, 1991 to December 31, 1991 using uncalibrated CN values. As shown in this study, the use of better input parameters for the L-THIA model can improve accuracy. The effects on direct runoff and pollutant estimation of the calibrated CN values in the L-THIA model were investigated for the LEC. Following calibration, the estimated average annual direct runoff for the LEC watershed increased by 34%, total nitrogen by 24%, total phosphorus by 22%, and total lead by 43%. This study demonstrates that the L-THIA model should be calibrated and validated prior to application in a particular watershed to more accurately assess the effects of land use changes on hydrology and water quality.     

CALIBRATION OF STORM LOADS IN THE SOUTH PRONG WATERSHED,FLORIDA, USING BASINS/HSPF1

ABSTRACT: The South Prong watershed is a major tributary system of the Sebastian River and adjacent Indian River Lagoon. Continued urbanization of the Sebastian River drainage basin and other watersheds of the Indian River Lagoon is expected to increase runoff and nonpoint source pollutant loads. The St. Johns River Water Management District developed watershed simulation models to estimate potential impacts on the ecological systems of receiving waters and to assist planners in devising strategies to prevent further degradation of water resources. In the South Prong system, a storm water sampling program was carried out to calibrate the water quality components of the watershed model for total suspended solids (TSS), total phosphorous (TP), and total nitrogen (TN). During the period of May to November 1999, water quality and flow data were collected at three locations within the watershed. Two of the sampling stations were located at the downstream end of major watercourses. The third station was located at the watershed outlet. Five storm events were sampled and measured at each station. Sampling was conducted at appropriate intervals to represent the rising limb, peak, and recession limb of each storm event. The simulations were handled by HSPF (Hydrologic Simulation Program‐Fortran). Results include calibration of the hydrology and calibration of the individual storm loads. The hydrologic calibration was continuous over the period 1994 through 1999. Simulated storm runoff, storm loads, and event mean concentrations were compared with their corresponding observed values. The hydrologic calibration showed good results. The outcome of the individual storm calibrations was mixed. Overall, however, the simulated storm loads agreed reasonably well with measured loads for a majority of the storms.     

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.     

NUTRIENT CONTRIBUTION OF NONPOINT SOURCE RUNOFF IN THE LAS VEGAS VALLEY1

ABSTRACT: A Geographic Information System (GIS) based non‐point source runoff model is developed for the Las Vegas Valley, Nevada, to estimate the nutrient loads during the years 2000 and 2001. The estimated nonpoint source loads are compared with current wastewater treatment facilities loads to determine the non‐point source contribution of total phosphorus (TP), total nitrogen (TN), and total suspended solids (TSS) on a monthly and annual time scale. An innovative calibration procedure is used to estimate the pollutant concentrations for different land uses based on available water quality data at the outlet. Results indicate that the pollutant concentrations are higher for the Las Vegas Valley than previous published values for semi‐arid and arid regions. The total TP and TN loads from nonpoint sources are approximately 15 percent and 4 percent, respectively, of the total load to the receiving water body, Lake Mead. The TP loads during wet periods approach the permitted loads from the wastewater treatment plants that discharge into Las Vegas Wash. In addition, the GIS model is used to track pollutant loads in the stream channels for one of the subwatersheds. This is useful for planning the location of Best Management Practices to control nonpoint pollutant loads.     

WATERSHED OPTIMIZATION OF AGRICULTURAL BEST MANAGEMENT PRACTICES: CONTINUOUS SIMULATION VERSUS DESIGN STORMS1

ABSTRACT: Nonpoint source (NPS) models and expert opinions are often used to prescribe best management practices (BMPs) for controlling NPS pollution. An optimization algorithm (e.g., a genetic algorithm, or GA) linked with a NPS model (e.g., Annualized AGricultural Nonpoint Source pollution model, or AnnAGNPS), can be used to more objectively prescribe BMPs and to optimize NPS pollution control measures by maximizing pollutant reduction and net monetary return from a watershed. Pollutant loads from design storms and annual loads from a continuous simulation can both be used for optimizing BMP schemes. However, which strategy results in a better solution (in terms of providing water quality protection) for a watershed is not clear. The specific objective of the study was to determine the differences in watershed pollutant loads, in an experimental watershed in Pennsylvania, resulting from optimization analyses performed using pollutant loads from a series of five 2‐yr 24‐hr storm events, a series of five 5‐yr 24‐hr storm events, and cumulative pollutant loads from a continuous simulation of five years of weather data. For each of these three different event alternatives, 100 near optimal solutions (BMP schemes) were generated. Sediment (Sed), sediment nitrogen (SedN), dissolved N (SolN), sediment organic carbon (SedOC), and sediment phosphorus (SedP) loads from a different five‐year period (an evaluation period) suggest that the optimal BMP schemes resulting from the use of annual cumulative pollutant loads from a continuous simulation of five years of weather data provide smaller cumulative NPS pollutant loads at the watershed outlet.     

APPLICATION OF THE CONTINUOUS STORMWATER POLLUTION SIMULATION SYSTEM (CSPSS): PHILADELLPHIA CASE STUDY1

ABSTRACT: This paper describes the Continuous Stormwater Pollution Simulation System (CSPSS) as well as a site-specific application of CSPSS to the Philadelphia urban area and its receiving water, the Delaware Estuary. Conceptually, CSPSS simulates the quantity and quality or urban stormwater runoff, combined sewer overflow, municipal and industrial waste water effuent, and upstream flow on a continuous basis for each time step in the simulation period. In addition, receiving water dissolved oxygen, suspended solids, and lead concentrations resulting from these pollutant sources may be simulated. However, only rceiving water dissolved oxygen (DO) response is considered in this paper. The continuous Do receiving water response model was calibrated to existing conditions usinv observed data at Chester, Pennsylvnia, located on the Delaware Estuary approximately 10 miles down stream from the study area. Average annual pollutant loads to the receiving water were estimated for all major sources and receiving water quality improvements resulting from removal of various portions of these pollutant loads were estimated by application of the calibrated stimultion model. It was found that the removal of oxygen-demanding pollutants from combined sewer overflow and urban stormwater runoff would result in relatively minor improvements in the overall dissoved oxygen resources of the Delaware Estuary; whereas. removal of oxygen demanding pollutants from waste water treatment plant effluent would result in greater improvemens. The results of this investigation can be used along with appropriate economic techniques to identify the most cost-effective mix of point and nonpoint source pollution control measures.     

ALAWAT: A spatially allocated watershed model for approximating stream,sediment, and pollutant flows in Hawaii,USA

The Ala Wai Canal Watershed Model (ALAWAT) is a planning-level watershed model for approximating direct runoff, streamflow, sediment loads, and loads for up to five pollutants. ALAWAT uses raster GIS data layers including land use, SCS soil hydrologic groups, annual rainfall, and subwatershed delineations as direct model parameter inputs and can use daily total rainfall from up to ten rain gauges and streamflow from up to ten stream gauges. ALAWAT uses a daily time step and can simulate flows for up to ten-year periods and for up to 50 subwatersheds. Pollutant loads are approximated using a user-defined combination of rating curve relationships, mean event concentrations, and loading/washoff parameters for specific subwatersheds, land uses, and times of year. Using ALAWAT, annual average streamflow and baseflow relationships and urban suspended sediment loads were approximated for the Ala Wai Canal watershed (about 10,400 acres) on the island of Oahu, Hawaii. Annual average urban suspended sediments were approximated using two methods: mean event concentrations and pollutant loading and washoff. Parameters for the pollutant loading and washoff method were then modified to simulate the effect of various street sweeping intervals on sediment loads.     

MODELING THE BUILDUP AND WASHOFF OF POLLUTANTS ON URBAN WATERSHEDS1

ABSTRACT: A model for urban stormwater quality was developed in this study. The basis for the model is the process by which pollutants build up on the watershed surface. For the wet climate of the study site, it was assumed that there exists an interval of time over which the pollutant buildup equals the pollutant washoff (no accumulation of pollutant). The buildup model was represented by a linear function of the antecedent dry time. The buildup function was then linked with a pollutant washoff model represented by a power function of the storm runoff volume. Various time intervals for no net accumulation were tested to calibrate the model. The model was calibrated to observed data for two small urban basins in Baton Rouge, Louisiana, and model results were used to analyze the behavior of phosphorus concentrations in storm runoff from these basins over a long period of time.     

Runoff Curve Numbers at the Agricultural Field‐Scale and Implications for Continuous Simulation Modeling

This study used monitoring in the waterways of agricultural fields to understand the use of the runoff curve number (CN) in continuous simulation models. The CN has a long history as a design tool for estimating runoff volumes for large, single storms on small watersheds, but its use in continuous simulation models to describe runoff from smaller storms and relatively small areas is more recent and controversial. We examined 788 nonwinter rainfall events on four agricultural fields over five years (2004‐2008) during which runoff was generated in 87 events. The largest 20 runoff events on each field generated approximately 90% of the total runoff volume. The runoff event CNs showed an inverse correlation with storm depth that could not consistently be explained by previous precipitation. We review how small areas of higher runoff generation within larger areas will systematically increase the apparent CN of the larger area as the storm size decreases. If this variation is not incorporated into a model explicitly, continuous simulation modelers must understand that when source areas are aggregated or when runoff generation is spatially variable, the overall CN is not unique when smaller storms are included in the calibration set.     

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.     

Climate Sensitivity of Phosphorus Loadings to an Urban Stream

We investigate the sensitivity of phosphorus loading (mass/time) in an urban stream to variations in climate using nondimensional sensitivity, known as elasticity, methods commonly used by economists and hydrologists. Previous analyses have used bivariate elasticity methods to represent the general relationship between nutrient loading and a variable of interest, but such bivariate relations cannot reflect the complex multivariate nonlinear relationships inherent among nutrients, precipitation, temperature, and streamflow. Using fixed‐effect multivariate regression methods, we obtain two phosphorus models (nonparametric and parametric) for an urban stream with high explanatory power that can both estimate phosphorus loads and the elasticity of phosphorus loading to changes in precipitation, temperature, and streamflow. A case study demonstrates total phosphorus loading depends significantly on season, rainfall, combined sewer overflow events, and flow rate, yet the elasticity of total phosphorus to all these factors remains relatively constant throughout the year. The elasticity estimates reported here can be used to examine how nutrient loads may change under future climate conditions.     

Dry Weather Water Quality Loadings in Arid,Urban Watersheds of the Los Angeles Basin,California, USA1

Abstract: Dry weather runoff in arid, urban watersheds may consist entirely of treated wastewater effluent and/or urban nonpoint source runoff, which can be a source of bacteria, nutrients, and metals to receiving waters. Most studies of urban runoff focus on stormwater, and few have evaluated the relative contribution and sources of dry weather pollutant loading for a range of constituents across multiple watersheds. This study assessed dry weather loading of nutrients, metals, and bacteria in six urban watersheds in the Los Angeles region of southern California to estimate relative sources of each constituent class and the proportion of total annual load that can be attributed to dry weather discharge. In each watershed, flow and water quality were sampled from storm drain and treated wastewater inputs, as well as from in‐stream locations during at least two time periods. Data were used to calculate mean concentrations and loads for various sources. Dry weather loads were compared with modeled wet weather loads under a range of annual rainfall volumes to estimate the relative contribution of dry weather load. Mean storm drain flows were comparable between all watersheds, and in all cases, approximately 20% of the flowing storm drains accounted for 80% of the daily volume. Wastewater reclamation plants (WRP) were the main source of nutrients, storm drains accounted for almost all the bacteria, and metals sources varied by constituent. In‐stream concentrations reflected major sources, for example nutrient concentrations were highest downstream of WRP discharges, while in‐stream metals concentrations were highest downstream of the storm drains with high metals loads. Comparison of wet vs. dry weather loading indicates that dry weather loading can be a significant source of metals, ranging from less than 20% during wet years to greater than 50% during dry years.     

Assessing Surface Water Quality and Its Relation with Urban Land Cover Changes in the Lake Calumet Area, Greater Chicago

Urban land use and land cover change significantly affect spatial and temporal patterns of runoff, which in turn impacts surface water quality. With the exponential growth in urban areas over the past three decades, changes in land use and land cover to cater for the growth of cities has been a conspicuous spectacle in urban spaces. The main goal of this study was to assess the impacts of land cover change on runoff and surface water quality using a partial area hydrology framework. The study employed ArcHydro GIS extension and a modified version of Long-Term Hydrologic and Nonpoint Source Pollution model (L-THIA-NPS) in estimating runoff and nonpoint source pollutant concentration around Lake Calumet between 1992 and 2001. Data employed include National Land Cover Data set, rainfall data, digital elevation model (DEM), Soil Survey Geographic (SSURGO) data, and The United States Environmental Protection Agency’s STORET (storage and retrieval) water quality data. The model was able to predict surface water quality reasonably well over the study period. Sensitivity analysis facilitated a manual calibration of the model. Model validation was executed by comparing simulated results following calibration and observed water quality data for the study area. The study demonstrates that the level of concentration of nonpoint source pollutants in surface water within an urban watershed heavily depends on the spatiotemporal variations in areas that contribute towards runoff compared to the spatial extent of change in major land use/land cover.     

Coupling PCSWMM and WASP to Evaluate Green Stormwater Infrastructure Impacts to Storm Sediment Loads in an Urban Watershed

We coupled rainfall–runoff and instream water quality models to evaluate total suspended solids (TSS) in Wissahickon Creek, a mid‐sized urban stream near Philadelphia, Pennsylvania. Using stormwater runoff and instream field data, we calibrated the model at a subdaily scale and focused on storm responses. We demonstrate that treating event mean concentrations as a calibration parameter rather than a fixed input can substantially improve model performance. Urban stormwater TSS concentrations vary widely in time and space and are difficult to represent simply. Suspended and deposited sediment pose independent stressors to stream biota and model results suggest that both currently impair stream health in Wissahickon Creek. Retrofitting existing detention basins to prioritize infiltration reduced instream TSS loads by 20%, suggesting that infiltration mitigates sediment more effectively than detention. Infiltrating stormwater from 30% of the watershed reduced instream TSS loads by 47% and cut the frequency of TSS exceeding 100 mg/L by half. Settled loads and the frequency of high TSS values were reduced by a smaller fraction than suspended loads and duration at high TSS values. A widely distributed network of infiltration‐focused projects is an effective stormwater management strategy to mitigate sediment stress. Coupling rainfall–runoff and water quality models is an important way to integrate watershed‐wide impacts and evaluate how management directly affects urban stream health.     

URBAN STORMWATER QUANTITY/QUALITY MODELING USING THE SCS METHOD AND EMPIRICAL EQUATIONS1

ABSTRACT: Runoff depth and pollutant loading (Biological Oxygen Demand [BOD₅], Total Suspended Solids [TSS], Total Kjeldahl Nitrogen [TKN] and lead [Pb]) computations of urban stormwater runoff from four small sites (i.e., 14.7–58.3 ac) in South Florida were performed using the Soil Conservation Service (SCS) hydrology method and empirical equations developed by the U.S. Environmental Protection Agency (EPA). Each site had different predominant land uses (i.e., low density residential, high density residential, highway and commercial). Quantity and quality data from 95 storm events at these sites were measured by the U.S. Geological Survey (USGS), and used for calibration of the methodology to derive appropriate input parameters. Calibrated input parameters were developed for each land use to test the applicability of the methodology in small sub-tropical urban watersheds, and to provide hydrologists with a way to select appropriate parameter values for planning studies. A total of 16 independent rainfall events were used for verification of the methodology. Comparisons of predicted versus measured data for both hydrographs and pollutant loadings were performed.     

Toward quantifying water pollution abatement in response to installing buffers on crop land

The scientific research literature is reviewed (i) for evidence of how much reduction in nonpoint source pollution can be achieved by installing buffers on crop land, (ii) to summarize important factors that can affect this response, and (iii) to identify remaining major information gaps that limit our ability to make probable estimates. This review is intended to clarify the current scientific foundation of the USDA and similar buffer programs designed in part for water pollution abatement and to highlight important research needs. At this time, research reports are lacking that quantify a change in pollutant amounts (concentration and/or load) in streams or lakes in response to converting portions of cropped land to buffers. Most evidence that such a change should occur is indirect, coming from site-scale studies of individual functions of buffers that act to retain pollutants from runoff: (1) reduce surface runoff from fields, (2) filter surface runoff from fields, (3) filter groundwater runoff from fields, (4) reduce bank erosion, and (5) filter stream water. The term filter is used here to encompass the range of specific processes that act to reduce pollutant amounts in runoff flow. A consensus of experimental research on functions of buffers clearly shows that they can substantially limit sediment runoff from fields, retain sediment and sediment-bound pollutants from surface runoff, and remove nitrate N from groundwater runoff. Less certain is the magnitude of these functions compared to the cultivated crop condition that buffers would replace within the context of buffer installation programs. Other evidence suggests that buffer installation can substantially reduce bank erosion sources of sediment under certain circumstances. Studies have yet to address the degree to which buffer installation can enhance channel processes that remove pollutants from stream flow. Mathematical models offer an alternative way to develop estimates for water quality changes in response to buffer installation. Numerous site conditions and buffer design factors have been identified that can determine the magnitude of each buffer function. Accurate models must be able to account for and integrate these functions and factors over whole watersheds. At this time, only pollutant runoff and surface filtration functions have been modeled to this extent. Capability is increasing as research data is produced, models become more comprehensive, and new techniques provide means to describe variable conditions across watersheds. A great deal of professional judgment is still required to extrapolate current knowledge of buffer functions into broadly accurate estimates of water pollution abatement in response to buffer installation on crop land. Much important research remains to be done to improve this capability. The greatest need is to produce direct quantitative evidence of this response. Such data would confirm the hypothesis and enable direct testing of watershed-scale prediction models as they become available. Further study of individual pollution control functions is also needed, particularly to generate comparative evidence for how much they can be manipulated through buffer installation and management.     

Integrated Modular Modeling of Water and Nutrients From Point and Nonpoint Sources in the Patuxent River Watershed1

Abstract: We present a simple modular landscape simulation model that is based on a watershed modeling framework in which different sets of processes occurring in a watershed can be simulated separately with different models. The model consists of three loosely coupled submodels: a rainfall‐runoff model (TOPMODEL) for runoff generation in a subwatershed, a nutrient model for estimation of nutrients from nonpoint sources in a subwatershed, and a stream network model for integration of point and nonpoint sources in the routing process. The model performance was evaluated using monitoring data in the watershed of the Patuxent River, a tributary to the Chesapeake Bay in Maryland, from July 1997 through August 1999. Despite its simplicity, the landscape model predictions of streamflow, and sediment and nutrient loads were as good as or better than those of the Hydrological Simulation Program‐Fortran model, one of the most widely used comprehensive watershed models. The landscape model was applied to predict discharges of water, sediment, silicate, organic carbon, nitrate, ammonium, organic nitrogen, total nitrogen, organic phosphorus, phosphate, and total phosphorus from the Patuxent watershed to its estuary. The predicted annual water discharge to the estuary was very close to the measured annual total in terms of percent errors for both years of the study period (≤2%). The model predictions for loads of nutrients were also good (20‐30%) or very good (<20%) with exceptions of sediment (40%), phosphate (36%), and organic carbon (53%) for Year 1.