# Patent

## Title
Method and system for water quality trading

## Abstract

A method and system for calculating trading ratios between a buyer and a seller of water quality pollutant load credits. The water quality trading system (“WQTS”) inputs a definition of watershed and land uses of the watershed. The WQTS then calculates the load credits (e.g., the difference between total maximum daily load of the pollutant and the actual load needed) of applying various controls (e.g., best management practices and treatments) to the point sources and non-point sources of the watershed. The load credits take into consideration the uncertainty of the calculated loads.

## Background
TECHNICAL FIELD
The described technology relates generally to analysis of pollutant loading in a watershed and particularly accounting for the uncertainty in both the pollutant loading and pollutant removal efficiency of best management practices implemented as controls and in determining trading ratios.
BACKGROUND
Many watersheds throughout the country are approaching or already exceeding their allowable pollutant loading capacities referred to as total maximum daily loads (“TMDLs”). As a result, it is becoming increasingly important for watershed management to provide an economical way for ensuring that the TMDLs are not exceeded. Watershed management is important for watersheds not only because TMDLs for many watersheds are currently defined, but also because many existing TMDLs are expected to be revised and TMDLs for many other watersheds are expected to be promulgated in the near future. The changes to the TMDLs will likely have a significant impact on the treatment costs for both point and non-point sources of pollutants. The United States Environmental Protection Agency published a Water Quality Trading Policy on Jan. 13, 2003, which is hereby incorporated by reference, stating that a market-based approach to water quality trading can provide flexibility and achieve water quality and environmental benefits greater than would otherwise be achieved under more traditional regulatory approaches. A market-based, water quality trading approach can be used by utilities, industries, county governments, and so on to achieve economically sound development and meet the EPA's water quality goals.
Any such water quality trading should factor in the high degree of uncertainty that results from various watershed processes such as rainfall runoff, soil erosion, sediment delivery, pollutant generation, and so on. It is desirable to have a tool to assist in analyzing water quality trading potentials that factor in this uncertainty.



BRIEF DESCRIPTION OF THE DRAWINGS
FIG. 1 illustrates an initial display page of the WQTS in one embodiment.
FIG. 2 illustrates a sample report listing the current and future pollutant loads for a selected confluence point in one embodiment.
FIG. 3 illustrates a display page for specifying the BMPs and treatments in one embodiment.
FIG. 4 illustrates a report listing the local credits for a subwatershed in one embodiment.
FIG. 5A illustrates a graph of the earned credits for the subwatersheds contributing to a confluence point in one embodiment.
FIG. 5B illustrates a report of the earned credits for a subwatershed contributing to a confluence point in one embodiment.
FIG. 6 illustrates a display page of a GIS interface for the WQTS in one embodiment.
FIG. 7 illustrates a map overlaid with the loading information of the subwatersheds in one embodiment.
FIG. 8 illustrates a map overlaid with the local credit information of the subwatersheds in one embodiment.
FIG. 9 illustrates a map overlaid with the earned credit information of the subwatersheds in one embodiment.
FIG. 10 is a block diagram illustrating the components of the WQTS in one embodiment.
FIG. 11 is a flow diagram illustrating the processing of a component that identifies subwatersheds in one embodiment.
FIG. 12 is a flow diagram illustrating the processing of a component that calculates the loading for the subwatersheds that are upstream from a confluence point in one embodiment.
FIG. 13 is a flow diagram posted in the processing of a component that calculates the local credits for the subwatersheds that are upstream from a confluence point in one embodiment.
FIG. 14 is a flow diagram illustrating the processing of a component that calculates the earned credits for subwatersheds that are upstream from a confluence point in one embodiment.



DETAILED DESCRIPTION
A method and system for calculating trading ratios between a buyer and a seller of water quality pollutant load credits is provided. In one embodiment, the water quality trading system (“WQTS”) inputs a definition of watershed and land uses of the watershed. The WQTS then calculates the load credits (e.g., the difference between the total maximum daily load (“TMDL”) of the pollutant and the actual load) of applying various controls (e.g., best management practices (“BMPs”) and treatments) to the point sources and non-point sources of the watershed. The load credits take into consideration the uncertainty of the calculated loads. The uncertainty can result from uncertainty in the export coefficient of a land use, in the removal efficiency of a BMP, and so on. The WQTS calculates an uncertainty ratio that will reduce the number of transferable credits generated by calculated non-point source reduction to account for variability in the effectiveness of the controlling practice (e.g., BMPs) implemented. For example, a seller may have a daily load credit of 500 lbs/day and the uncertainty is 0.7, then the available credits for sale are 500×0.7=350 lbs/day. To account for the difference between the uncertainty ratios of the available credits of a seller and the needed credits of a buyer, the WQTS calculates a trading ratio to indicate the equivalency between the available and needed credits. In one embodiment, the trading ratio is the exchange ratio of the seller divided by the exchange ratio of the buyer. When the exchange ratio of both the buyer and the seller are the same (e.g., both 0.7), then the trading ratio is 1, which means the available credits can be traded on a one-for-one basis with the needed credits. This happens mostly when the sellers and buyers are located in the same vicinity. When the trading ratio is greater than one, then an available credit is worth more than a needed credit. For example, if the trading ratio is 2, then the buyer would need to purchase 2 available credits from the seller for every needed credit. In contrast, when the trading ratio is less than one, then an available credit is worth less than a needed credit. For example, if the trading ratio is 0.5, then the buyer would need to purchase 1 seller credit for every 2 needed credits. The WQTS thus provides a way to calculate the equivalency between available credits and needed credits so that potential credit trades can be evaluated.
The WQTS allows a user to evaluate the feasibility of conducting pollutant trades while accounting for the economics of the trade. The WQTS models the pollutants of a watershed taking into account the contributing pollutants and their loads, the contributing sources and their locations, the uncertainty in the pollutant loads, and BMPs for non-point sources and associated uncertainty. The WQTS calculates credits and cost differentials among various sources and estimates trading ratios for buyers and sellers within a watershed.
The WQTS initially calculates the pollutant loading of the subwatersheds of a designated watershed for use when analyzing the effects of various controls. In one embodiment, the WQTS uses a database that contains a description of river basins that specifies each confluence point of a river basin and each watershed and each subwatershed. After a user selects a confluence point at which the trading analysis is to be performed, the WQTS identifies upstream subwatersheds that contribute to that confluence point. The WQTS then calculates the pollutant load of each contributing subwatershed so that the effect of the controls can be analyzed. The database includes for each subwatershed a description of each current and planned land use (including the area and export coefficient of the pollutant) and the load of each current and planned point source. In one embodiment, the WQTS calculates the total current and future pollutant loads for the point and non-point sources of a subwatershed. The current and future pollutant loads for point sources take into consideration the current point sources and the planned point sources (e.g., add a water treatment plant or remove the secondary treatment of an existing water treatment plant) for a subwatershed. The current and future pollutant loads for non-point sources take into consideration the current and planned land uses in the subwatershed. The WQTS also calculates an estimated future load at 50% and 95% confidence levels.
The WQTS allows a user to specify the modifications (e.g., adding, removing, or changing) of the controls such as the BMPs of the land uses and the treatment of the point sources of a subwatershed. The WQTS allows a user to specify the land uses of a subwatershed and the percent of the subwatershed affected and to specify the BMPs that are to be applied to each land use, the year of installation, and whether the BMPs are to be in serial or in parallel. The WQTS also allows a user to specify the types of treatments along with the year of installation that are to be applied to point sources of a subwatershed.
After the BMPs and treatments are defined, the WQTS calculates the local load credit of each subwatershed. The WQTS calculates the load of each land use of a watershed taking into consideration the percent of the subwatershed affected by the land use, the area of the land use, the export coefficient of the land use, and the pollutant removal efficiency of each BMP. The total load for the subwatershed is the sum of the load of its land uses and its point sources. The “local credit” of a subwatershed is the difference between the load prior to and after application of the controls and represents the cumulative effect of the controls of a subwatershed.
When credits are to be traded, the WQTS calculates the “earned credit” of a subwatershed that represents the portion of the local credits that are realized at the confluence point taking into consideration the uncertainty. Various factors influence what percent of a local load contributes to the load at the confluence point. The aggregate of the factors is referred to as the “exchange ratio.” The load of a subwatershed that contributes to a confluence point is the exchange ratio times the local load, and the earned credit of a subwatershed at a confluence point is the difference between the load prior to and after application of the controls multiplied by the exchange ratio. The exchange ratio may aggregate the effects of the “uncertainty ratio,” “delivery ratio,” and so on. The uncertainty ratio represents the uncertainty in the calculated load of a subwatershed resulting from the variation in various parameters including the export coefficient of a land use and the removal efficiency of a BMP. The uncertainty ratio represents the percent of the load that will be generated with a certain confidence level. For example, if the load is calculated to be 100, then the load with a 50% confidence level may be 90 and the load with a 95% confidence level may be 75. The delivery ratio represents the fraction of the load of a subwatershed that is actually delivered to the confluence point. Each downstream subwatershed on the way to the confluence point may only deliver a certain percent of its inflow to its outflow, which is represented by its delivery ratio. Thus, the effective delivery ratio of a subwatershed to a confluence point is calculated based on the delivery ratios of each downstream subwatershed. As discussed above, the WQTS calculates a trading ratio based on the exchange ratio for the available earned credits of a seller and the needed earned credits of a buyer at a certain confidence level. The exchange ratio can be used to determine the relative value of the available and needed credits.
In one embodiment, the WQTS provides the cost of various BMPs for land uses and treatments for point sources. A user can specify various combinations of BMPs and treatments and analyze their cost-effectiveness. The WQTS in one embodiment may automatically identify an optimum combination of BMPs and treatments based on various constraints specified by a user. The optimum combination can be identified using various well-known techniques such as linear programming.
In one embodiment, the WQTS provides a geographic information system (“GIS”) that allows a user to identify a river basin and a confluence point of the river basin that is to be analyzed. The GIS database describes the river basins, the confluence points, the subwatersheds, the land uses (current and future), and the point sources. A river basin table provides geographic information of the river basin. The confluence point table provides for each river basin the location of each confluence point along with the “from” subwatersheds whose outflows are to the confluence point and the “to” subwatersheds whose inflows are from the confluence point. The subwatershed table provides geographic information for each subwatershed. The current land use table and the future land use table describe each current and future land use of a subwatershed and include its area and export coefficient. The point source table provides the current and future loads of each point source.
The WQTS includes other database tables that describe the BMPs, the local delivery ratios, the confidence levels, the treatment types, the export coefficient and removal efficiency variations, and the TMDL at each confluence point. The BMP table provides the removal efficiency, the life, the cost, and the coefficient of variation of the removal efficiency for each BMP. The delivery ratio table provides the delivery ratio for each subwatershed. The confidence level table provides the z value of the standard normal distribution at various confidence levels. One skilled in the art will appreciate that different types of distribution can be used. The treatment type table provides the concentration of the pollutant, the life, and the cost of various treatment types. The export coefficient variations table provides the coefficient of variation of the export coefficient of each land use. The TMDL table provides the TMDL of the pollutant at each confluence point.
FIG. 1 illustrates an initial display page of the WQTS in one embodiment. The display page 100 includes a river basin field 101, a confluence point field 102, and a calculation level field 103. The river basin field provides a drop-down list of the river basins defined in the database and allows a user to select the river basin to be analyzed. The confluence point field provides a drop-down list of the confluence points within the selected river basin as defined by the database. The user selects the confluence point at which the pollutant load is to be analyzed. The calculation level field specifies the level (e.g., watershed or subwatershed) at which the analysis is to be performed. (This description assumes that the user has selected the subwatershed calculation level.) The display page also includes a loading button 104, a controls button 105, a trading button 106, and a requirements button 107.
After a user selects the river basin, the confluence point, and the calculation level, the user selects the loading button to direct the WQTS to calculate the current and future pollutant loads for subwatersheds that contribute to the selected confluence point. The WQTS generates a report that lists for each subwatershed the total current load and total future load for the point and non-point sources. FIG. 2 illustrates a sample report listing the current and future pollutant loads for a selected confluence point in one embodiment. The report also lists the total future loads with a 50% and 95% confidence level for non-point sources. The report lists the loads for each subwatershed and then provides the totals at the confluence points.
After a user selects the river basin, the confluence point, the calculation level, and then the loading button, the user selects the controls button to specify the BMP to be applied to the land uses of a subwatershed and the treatments to be applied to point sources of the subwatershed. FIG. 3 illustrates a display page for specifying the BMPs and treatments in one embodiment. The display page 300 includes a subwatershed field 301, non-point source area 310, and point source area 320. The subwatershed field identifies the subwatershed currently selected. Scroll buttons 330 allow a user to select different subwatersheds that are upstream from the selected confluence point. The non-point source area includes a land use area 311 with a land use field 312, a proportion controlled field 313, and land use scroll buttons 314. The land use scroll buttons allow the user to scroll through and add and delete land uses. The land use field provides a drop-down list through which a user can select a land use of the selected subwatershed. The proportion controlled field 313 specifies the portion of the subwatershed that is controlled by the land use. A BMP field 317 lists the BMPs that have been applied to the selected land use and provides a drop-down list of BMPs that can be applied. An initial year field 318 specifies the starting year of the associated BMP. A BMP area 316 lists the BMPs of the selected land use. The BMPs, shown in the BMP area at the same time, represent BMPs that are applied in series to the land use. BMP scroll buttons 319 allow the user to scroll through sets of serial BMPs that are applied in parallel to the land use. The point source area includes a point source field 321, a treatment type field 322, and an installation year field 323. The point source field provides a drop-down list of the point sources of the selected subwatershed so that the user can select a point source to which a treatment is to be supplied. The treatment type field provides a drop-down list of the possible treatment types of the point source. The installation year field allows the user to specify when the selected treatment is installed. Point source buttons 324 allow a user to select, add, and delete different point sources. A view results button 331 allows a user to direct the WQTS to calculate the local credits of each subwatershed based on the specified controls (e.g., BMP and treatments).
FIG. 4 illustrates a report listing the local credits for a subwatershed in one embodiment. The report is for the currently selected subwatershed. The report lists future load with and without the controls applied, with their difference being the local credit. The report also summarizes the effect of the point source and non-point source controls. The summary of the non-point source controls includes an identification of the BMP, the install year, the life, the efficiency, the load reduction, the cost per pound of reduction, and the total cost of the control. The summary of the point source controls includes an identification of the point source, the treatment type, the effluent concentration, the install year, and the total cost of the control. This report can be generated at the subwatershed, watershed, and river basin levels depending on the selected calculation level.
After a user selects the river basin, the confluence point, the calculation level, and then the loading and controls buttons, the user selects the trading button to direct the WQTS to perform the trading calculations. The WQTS requests the user to select a confidence level and then calculates and displays the earned credits (aka the “watershed credits”) at the selected confluence point for each subwatershed, factoring in the uncertainties and the selected confidence level. FIG. 5A illustrates a graph of the earned credits for the subwatersheds contributing to a confluence point in one embodiment. The graph provides an indication of the credits versus year for each subwatershed. Although not illustrated in the figures the WQTS also allows a user to select pairs of subwatersheds for which trading ratios are to be calculated. FIG. 5B illustrates a report of the earned credits for a subwatershed contributing to a confluence point in one embodiment. The report includes the local credits, the earned credits, and the uncertainty ratios for each year.
FIG. 6 illustrates a display page of a GIS interface for the WQTS in one embodiment. The display page 600 displays a map 601 of the subwatersheds contributing to a user-selected confluence point in a user-selected river basin. A layers area 602 allows a user to specify the information that is to be overlaid onto the map by selecting the appropriate boxes. For example, the user may select to overlay onto the map an outline of each subwatershed. After the trading information has been generated using the display page 100 of FIG. 1, the WQTT toolbar 603 allows the user to display the results. FIG. 7 illustrates a map overlaid with the loading information of the subwatersheds in one embodiment. FIG. 8 illustrates a map overlaid with the local credit information of the subwatersheds in one embodiment. FIG. 9 illustrates a map overlaid with the earned credit information of the subwatersheds in one embodiment.
The following scenario provides an overview of a typical use of the WQTS. A person may need to offset an increase in the pollutant load of a point source. That person may have the option of adding certain BMPs to land uses under their control to offset the increase. Alternatively, it may be more cost-effective to apply BMPs to land uses under the control of others. That person would use WQTS to enter the BMP information for their land use to generate the required earned credit to offset the increased load. That person would enter the BMP information and select a confidence level and be presented with cost and the earned credits with an associated uncertainty ratio. That person can then repeat the process for the BMPs that might be applied to others. That person can then use the trading ratios of the BMPs applied to their land uses and those of others to determine the equivalency of the earned credits. Assuming that equivalent credits can be obtained at a lower cost by applying BMPs to the land uses of others, then that person would attempt to trade credits with the other person who controls that land use.
Loading Calculations
The loading calculations performed by the WQTS in one embodiment are described in the following.
The total current load of a pollutant for a subwatershed is calculated by summing the loads generated from each individual land use of the subwatershed. The load of a land use is calculated by multiplying the area of the land use in the subwatershed by the export coefficient of the land use. The export coefficient represents the amount of a pollutant generated by the land use in mass per unit area per time duration (e.g., pounds/acre/year). Thus, the total current load for a subwatershed is represented by the following equation:




  
    
      
        
          L
          i
        
        =
        
          
            ∑
            
              j
              =
              1
            
            
              n
              i
            
          
          ⁢
          
            
          
          ⁢
          
            (
            
              
                A
                ij
              
              *
              
                C
                j
              
            
            )
          
        
      
    
    
      
        (
        1
        )
      
    
  




where Li is the total current load of subwatershed i, Aij is the area of land use j in subwatershed i, Cj is the export coefficient of land use j, and ni is the number of land uses in subwatershed i.

The total future load of a pollutant for a subwatershed is calculated in a similar manner as for the total current load, except using future land uses. The total future load for a subwatershed is represented by the following equation:




  
    
      
        
          FL
          i
        
        =
        
          
            ∑
            
              j
              =
              1
            
            
              n
              i
            
          
          ⁢
          
            
          
          ⁢
          
            (
            
              
                A
                ij
              
              *
              
                C
                j
              
            
            )
          
        
      
    
    
      
        (
        2
        )
      
    
  




where FLi is the total future load of subwatershed i, Aij is the area of future land use j in subwatershed i, Cj is the export coefficient of land use j, and ni is the number of future land uses in subwatershed i.

The standard deviation of the future load of a pollutant of a subwatershed is used to calculate future load at certain confidence levels. The standard deviation is the square root of the variance, which is calculated by summing the square of the product of future average load of each land use with the coefficient of variation for that land use. Thus, the standard deviation of the future load is represented by the following equation:




  
    
      
        
          σ
          i
        
        =
        
          
            [
            
              
                ∑
                
                  i
                  =
                  1
                
                
                  n
                  i
                
              
              ⁢
              
                
              
              ⁢
              
                
                  (
                  
                    
                      A
                      ij
                    
                    *
                    
                      
                        C
                        _
                      
                      j
                    
                    *
                    
                      CVC
                      j
                    
                  
                  )
                
                2
              
            
            ]
          
          
            1
            /
            2
          
        
      
    
    
      
        (
        3
        )
      
    
  




where σi is the standard deviation of the future load of subwatershed i, Aij is the area of future land use j within subwatershed i, {overscore (C)}j is the mean export coefficient of land use j, ni is the number of future land uses in subwatershed i, and CVCj is the coefficient of variability for land use j.

The future load of a subwatershed at a certain confidence level is calculated by multiplying the future load by the standard deviation of the subwatershed by the z value of the confidence level. The future load at a certain confidence level is represented by the following equation:

FLi(x)=FLi[1+z(x)σi]  (4)

where FLi(x) is the future load of subwatershed i at confidence level x, FLi is the future load of subwatershed i, σi is the standard deviation of the future load of subwatershed i, and z(x) is the z value for confidence level x.

The current point source load of a subwatershed is calculated by summing the loads of each point source in the subwatershed. The point source load is represented by the following equation:




  
    
      
        
          PSL
          i
        
        =
        
          
            ∑
            
              j
              =
              1
            
            
              m
              i
            
          
          ⁢
          
            
          
          ⁢
          
            L
            ij
          
        
      
    
    
      
        (
        5
        )
      
    
  




where PSLi is the current point source load of subwatershed i, Lij is the current load of current point source j in subwatershed i, and mi is the number of current point sources in subwatershed i.

The future point source load of a subwatershed is calculated in a similar manner to the current point source load and is represented by the following equation:




  
    
      
        
          FPSL
          i
        
        =
        
          
            ∑
            
              j
              =
              1
            
            
              m
              i
            
          
          ⁢
          
            
          
          ⁢
          
            L
            ij
          
        
      
    
    
      
        (
        6
        )
      
    
  




where FPSLi is the future point source load of subwatershed i, Lij is the future load of future point source j in subwatershed i, and mi is the number of future point sources in subwatershed i.

The delivery ratio represents the ratio of the load of a subwatershed that makes it downstream to a confluence point. Each subwatershed has a local delivery ratio that represents the portion of the inflow load that makes it to the outlets of the individual watersheds. The delivery ratio is calculated by multiplying the local delivery ratio of each downstream subwatershed together and is represented by the following equation:




  
    
      
        
          DR
          i
        
        =
        
          
            ∏
            
              j
              =
              1
            
            
              n
              i
            
          
          ⁢
          
            
          
          ⁢
          
            LDR
            j
          
        
      
    
    
      
        (
        7
        )
      
    
  




where DRi is the delivery ratio of subwatershed i, LDRj is the local delivery ratio of downstream subwatershed j, and ni is the number of downstream subwatersheds of subwatershed j, including subwatershed i itself.

The point source load of a subwatershed that makes it to the confluence point is the point source load of the subwatershed multiplied by the delivery ratio. The point source load of a subwatershed contributing to a confluence point is represented by the following equation:

PSi=PSLi*DRi  (8)

where PSi is the point source load contribution of subwatershed i at the confluence point, where PSLi is the current point source load of subwatershed i, and DRi is the delivery ratio of subwatershed i.

The future point source load of a subwatershed that makes it to the confluence point is calculated in a manner similar to the way the current point source load is calculated and is represented by the following equation:

FPSi=FPSLi*DRi  (9)

where FPSi is the future point source load contribution of subwatershed i at the confluence point, FPSLi is the future point source load of subwatershed i, and DRi is the delivery ratio of subwatershed i to the confluence point.

Local and Earned Credit Calculations

The local and earned credit calculations performed by the WQTS in one embodiment are described in the following. The non-point source load from a subwatershed having at most one BMP per land use is calculated based on the removal efficiency of the BMP and the portion of the land use affected by the BMP and is represented by the following equation:

Lij=Cj*Aij*(1−ηj)*fij+Cj*Aij*(1−fij)  (10)

where Lij is the load of subwatershed i from land use j, Cj is the export coefficient of land use j, Aij is the area of land use j in subwatershed i, and fij is the fraction of land use j in subwatershed i treated with a BMP having a pollutant removal efficiency of ηj.

The mean (or average) load per interval (e.g., season) of a subwatershed over a period (e.g., one year) is calculated based on the mean of the export coefficient and the removal efficiency. The mean load of a subwatershed is represented by the following equation:

{overscore (L)}ij={overscore (C)}j*Aij*(1−{overscore (η)}j)*fij+{overscore (C)}j*Aij*(1−fij)  (11)

where {overscore (L)}ij is the mean load of subwatershed i from land use j, {overscore (C)}j is the mean export coefficient of land use j, Aij is the area of land use j in subwatershed i, and fij is the fraction of land use j in subwatershed i treated with a BMP having a mean pollutant removal efficiency of {overscore (η)}j.

The coefficient of variation of the load of a subwatershed is derived based on the coefficients of variation of the export efficiency and the removal efficiency and is represented by the following equation:




  
    
      
        
          CVL
          ij
        
        =
        
          
            (
            
              
                
                  (
                  
                    1
                    +
                    
                      CVC
                      j
                      2
                    
                  
                  )
                
                ⁢
                
                  (
                  
                    1
                    +
                    
                      
                        
                          f
                          ij
                          2
                        
                        *
                        
                          
                            (
                            
                              1
                              -
                              
                                
                                  η
                                  _
                                
                                j
                              
                            
                            )
                          
                          2
                        
                        *
                        CV
                        ⁢
                        
                          
                        
                        ⁢
                        
                          η
                          j
                          2
                        
                      
                      
                        
                          (
                          
                            1
                            -
                            
                              f
                              ij
                            
                            +
                            
                              
                                (
                                
                                  1
                                  -
                                  
                                    
                                      η
                                      _
                                    
                                    j
                                  
                                
                                )
                              
                              *
                              
                                f
                                ij
                              
                            
                          
                          )
                        
                        2
                      
                    
                  
                  )
                
              
              -
              1
            
            )
          
          
            1
            /
            2
          
        
      
    
    
      
        (
        12
        )
      
    
  




where CVLij is the coefficient of variation of the load generated from land use j within subwatershed i, CVCj is the coefficient of variation of the export coefficient of land use j, CVηj is the coefficient of variation of the removal efficiency of the BMP, and fij is the fraction of land use j within subwatershed i treated with the BMP having a mean pollutant removal efficiency of {overscore (η)}j.

Equations (11) and (12) for the mean and coefficient of variation of the land use assume implementation of only one BMP per land use as a pollution controlling unit. When multiple BMPs are used in series on a land use, then the removal efficiency is calculated as one minus each of the inefficiencies (i.e., 1−ηj) multiplied together. The mean removal efficiency of the group of multiple BMPs in series combination is represented by the following equation:




  
    
      
        
          
            η
            _
          
          ij
        
        =
        
          1
          -
          
            
              ∏
              
                k
                =
                1
              
              
                η
                ij
              
            
            ⁢
            
              
            
            ⁢
            
              (
              
                1
                -
                
                  
                    η
                    _
                  
                  jk
                
              
              )
            
          
        
      
    
    
      
        (
        
          13
          ⁢
          A
        
        )
      
    
  




where {overscore (η)}ij is the mean removal efficiency of the serial BMPs of land use j in subwatershed i, {overscore (η)}jk is the mean removal efficiency of BMP k in land use j, and nij is the number of BMPs of land use j in subwatershed i. Using equations (13A) and (11), the mean pollutant loading from a land use having a group of multiple BMPs in series is represented by the following equation:





  
    
      
        
          L
          ij
        
        =
        
          
            
              C
              j
            
            *
            
              A
              ij
            
            *
            
              f
              ij
            
            ⁢
            
              
                ∏
                
                  k
                  =
                  1
                
                
                  n
                  ij
                
              
              ⁢
              
                
              
              ⁢
              
                (
                
                  1
                  -
                  
                    
                      η
                      _
                    
                    jk
                  
                
                )
              
            
          
          +
          
            
              C
              j
            
            *
            
              A
              ij
            
            *
            
              (
              
                1
                -
                
                  f
                  ij
                
              
              )
            
          
        
      
    
    
      
        (
        
          13
          ⁢
          B
        
        )
      
    
  




where {overscore (L)}ij is the mean load of subwatershed i from land use j, Cj is the export coefficient of land use j, Aij is the area of land use j in subwatershed i, fij is the fraction of land use j in subwatershed i treated with a BMP having a mean pollutant removal efficiency of {overscore (η)}j, and nij is the number of BMPs of land use j in subwatershed i.

The coefficient of variation of the removal efficiency for serial BMPs is the product of one plus the variance of the removal efficiency divided by the inefficiency for each BMP. That coefficient of variation is represented by the following equation:




  
    
      
        
          CV
          ⁢
          
            
          
          ⁢
          
            η
            ij
          
        
        =
        
          
            ∏
            
              k
              =
              1
            
            
              η
              ij
            
          
          ⁢
          
            
          
          ⁢
          
            
              (
              
                1
                +
                
                  
                    Var
                    ⁡
                    
                      (
                      
                        η
                        jk
                      
                      )
                    
                  
                  
                    
                      (
                      
                        1
                        -
                        
                          
                            η
                            _
                          
                          jk
                        
                      
                      )
                    
                    2
                  
                
              
              )
            
            
              1
              /
              2
            
          
        
      
    
    
      
        (
        
          14
          ⁢
          A
        
        )
      
    
  




where CVηij is the coefficient of variation of the removal efficiency of land use j in subwatershed i, Var(ηjk) is the variance of the removal efficiency for BMP k in land use j, {overscore (η)}jk is the mean removal efficiency for BMP k in land use j, and nij is the number of BMPs of land use j in subwatershed i. Using equations (14A) and (12), the coefficients of variation of pollutant loading from land use j having a group of multiple BMPs in series are represented by the following equation:





  
    
      
        
          CVL
          ij
        
        =
        
          
            {
            
              
                
                  (
                  
                    1
                    +
                    
                      CVC
                      j
                      2
                    
                  
                  )
                
                ⁡
                
                  [
                  
                    1
                    +
                    
                      
                        
                          f
                          ij
                          2
                        
                        ⁢
                        
                          
                            ∏
                            
                              k
                              =
                              1
                            
                            
                              n
                              ij
                            
                          
                          ⁢
                          
                            
                          
                          ⁢
                          
                            
                              
                                (
                                
                                  1
                                  -
                                  
                                    
                                      η
                                      _
                                    
                                    jk
                                  
                                
                                )
                              
                              2
                            
                            ⁢
                            
                              
                                ∏
                                
                                  k
                                  =
                                  1
                                
                                
                                  n
                                  ij
                                
                              
                              ⁢
                              
                                
                              
                              ⁢
                              
                                (
                                
                                  1
                                  +
                                  
                                    
                                      Var
                                      ⁡
                                      
                                        (
                                        
                                          η
                                          jk
                                        
                                        )
                                      
                                    
                                    
                                      
                                        (
                                        
                                          1
                                          -
                                          
                                            
                                              η
                                              _
                                            
                                            jk
                                          
                                        
                                        )
                                      
                                      2
                                    
                                  
                                
                                )
                              
                            
                          
                        
                      
                      
                        
                          (
                          
                            1
                            -
                            
                              f
                              ij
                            
                            +
                            
                              
                                f
                                ij
                              
                              ⁢
                              
                                
                                  ∏
                                  
                                    k
                                    =
                                    1
                                  
                                  
                                    n
                                    ij
                                  
                                
                                ⁢
                                
                                  
                                
                                ⁢
                                
                                  (
                                  
                                    1
                                    -
                                    
                                      
                                        η
                                        _
                                      
                                      jk
                                    
                                  
                                  )
                                
                              
                            
                          
                          )
                        
                        2
                      
                    
                  
                  ]
                
              
              -
              1
            
            }
          
          
            1
            /
            2
          
        
      
    
    
      
        (
        
          14
          ⁢
          B
        
        )
      
    
  




where CVLij is the coefficient of variation of the load of subwatershed i, CVCj is the coefficient of variation of the export coefficient of land use j, Aij is the area of land use j in subwatershed i, {overscore (η)}jk is the mean removal efficiency of BMP k in land use j, and nij is the number of BMPs of land use j in subwatershed i.

When multiple BMPs are in parallel on a land use, then the mean load for a land use in the subwatershed is the product of the mean export coefficient of that land use for the subwatershed times the sum of the areas of the land use adjusted by the mean removal inefficiency. The load with multiple BMPs in parallel is represented by the following equation:




  
    
      
        
          
            L
            _
          
          ij
        
        =
        
          
            
              C
              _
            
            j
          
          *
          
            
              ∑
              
                k
                =
                1
              
              
                n
                ij
              
            
            ⁢
            
              
            
            ⁢
            
              
                A
                ij
              
              *
              
                (
                
                  1
                  -
                  
                    
                      η
                      _
                    
                    jk
                  
                
                )
              
            
          
        
      
    
    
      
        (
        15
        )
      
    
  




where {overscore (L)}ij is the mean load of subwatershed i from land use j, {overscore (C)}j is the mean export coefficient of land use j, Aij is the area of land use j in subwatershed i, {overscore (η)}jk is the mean removal efficiency for BMP k in land use j, and nij is the number of BMPs in parallel of land use j in subwatershed i.

The coefficient of variation of pollutant loading for BMPs in parallel is represented by the following equation:




  
    
      
        
          CVL
          ij
        
        =
        
          
            (
            
              
                
                  (
                  
                    1
                    +
                    
                      CVC
                      j
                      2
                    
                  
                  )
                
                ⁢
                
                  (
                  
                    1
                    +
                    
                      
                        
                          ∑
                          
                            k
                            =
                            1
                          
                          
                            n
                            ij
                          
                        
                        ⁢
                        
                          
                        
                        ⁢
                        
                          
                            A
                            ij
                            2
                          
                          *
                          
                            
                              (
                              
                                1
                                -
                                
                                  
                                    η
                                    _
                                  
                                  jk
                                
                              
                              )
                            
                            2
                          
                          *
                          CV
                          ⁢
                          
                            
                          
                          ⁢
                          
                            η
                            jk
                            2
                          
                        
                      
                      
                        (
                        
                          
                            ∑
                            
                              k
                              =
                              1
                            
                            
                              n
                              ij
                            
                          
                          ⁢
                          
                            
                          
                          ⁢
                          
                            
                              A
                              ij
                            
                            *
                            
                              
                                (
                                
                                  1
                                  -
                                  
                                    
                                      η
                                      _
                                    
                                    jk
                                  
                                
                                )
                              
                              2
                            
                          
                        
                      
                    
                  
                  )
                
              
              -
              1
            
            )
          
          
            1
            /
            2
          
        
      
    
    
      
        (
        16
        )
      
    
  




where CVLij is the coefficient of variation of the load of subwatershed i from land use j, CVCj is the coefficient of variation of the export coefficient of land use j, Aij is the area of land use j in subwatershed i, {overscore (η)}jk is the mean removal efficiency of BMP k in land use j, CVηjk is the coefficient of variation of BMP k in land use j, and nij is the number of BMPs of land use j in subwatershed i.

The mean pollutant load of a subwatershed is the sum of the mean pollutant loads from each land use within the subwatershed and is represented by the following equation:




  
    
      
        
          
            L
            _
          
          i
        
        =
        
          
            ∑
            
              j
              =
              1
            
            
              m
              1
            
          
          ⁢
          
            
          
          ⁢
          
            
              L
              _
            
            ij
          
        
      
    
    
      
        (
        17
        )
      
    
  




where {overscore (L)}i is the mean pollutant load of subwatershed i, {overscore (L)}ij is the mean pollutant load of land use j in subwatershed i, and mi is the number of land uses in subwatershed i.

The variance of the total load from a subwatershed is represented by the following equation:




  
    
      
        
          Var
          ⁡
          
            (
            
              L
              i
            
            )
          
        
        =
        
          
            ∑
            
              j
              =
              1
            
            
              m
              i
            
          
          ⁢
          
            
          
          ⁢
          
            
              (
              
                
                  
                    L
                    _
                  
                  i
                
                *
                
                  CVL
                  ij
                
              
              )
            
            2
          
        
      
    
    
      
        (
        18
        )
      
    
  




where Var(Li) is the variance of the pollutant load of subwatershed i, {overscore (L)}i is the mean pollutant load of subwatershed i, CVLij is the coefficient of variation of the load of land use j in subwatershed i, and mi is the number of land uses in subwatershed i.

The coefficient of variation of the total load from a subwatershed is represented by the following equation:




  
    
      
        
          CVL
          i
        
        =
        
          
            
              (
              
                Var
                ⁡
                
                  (
                  
                    L
                    i
                  
                  )
                
              
              )
            
            
              1
              /
              2
            
          
          
            
              L
              _
            
            i
          
        
      
    
    
      
        (
        19
        )
      
    
  




where CVLi is the coefficient of variation of the load of subwatershed i, Var(Li) is the variance of the pollutant load of subwatershed i, and {overscore (L)}i is the mean load of subwatershed i.

The one-sided lower and upper confidence intervals for an assumed distribution of the load delivered from a subwatershed at its outlet is represented by the following equation:

CIi={overscore (L)}i±zα*{overscore (L)}i*[Var(Li)]1/2={overscore (L)}i[1±zα*[Var(Li)]1/2]  (20)

where CIi is the confidence interval for the load delivered from the outlet of subwatershed i, {overscore (L)}i is the mean load of subwatershed i, Var(Li) is the variance of the pollutant load of subwatershed i, and zα is the z value at significance level α.

The two-sided lower and upper confidence interval for an assumed distribution of the load delivered from a subwatershed at its outlet is represented by the following equation:

CIi={overscore (L)}i±zα/2*{overscore (L)}i*[Var(Li)]1/2={overscore (L)}i[1±zα/2*[Var(Li)]1/2]  (21)

where CIi is the confidence interval for the load delivered from the outlet of subwatershed i, {overscore (L)}i is the mean load of subwatershed i, Var(Li) is the variance of the pollutant load of subwatershed i, and zα/2 is the z value at significance level α/2.

The uncertainty ratio for a non-point source accounts for the variability in load delivered from a subwatershed due to variation in the export coefficient of a land use and the removal efficiency of a BMP. The uncertainty ratio is calculated by dividing the mean load by the mean load times one plus the z value at a certain confidence level (e.g., 95%) times the coefficient of variation of the load. When the value of zα corresponding to the 95th percentile of the total load distribution is approximately 2, the uncertainty ratio for the non-point source load delivered from a subwatershed is represented by the following equation:




  
    
      
        
          UR
          i
        
        =
        
          
            
              
                L
                _
              
              i
            
            
              
                
                  L
                  _
                
                i
              
              ⁡
              
                [
                
                  1
                  +
                  
                    
                      z
                      α
                    
                    *
                    
                      CVL
                      i
                    
                  
                
                ]
              
            
          
          =
          
            1
            
              1
              +
              
                2
                *
                
                  CVL
                  i
                
              
            
          
        
      
    
    
      
        (
        22
        )
      
    
  




where URi is the uncertainty ratio for subwatershed i, {overscore (L)}i is the mean load of subwatershed i, CVLi is the coefficient of variation of the load of subwatershed i, and zα is the z value at a certain confidence level.

The uncertainty ratio for a point source, such as a wastewater treatment plant, is assumed to have negligible variability in the load delivered to the outlet of a subwatershed (i.e., CV is equal to 0). Therefore, the uncertainty ratio due to the point sources of a subwatershed is represented by the following equation:

URi=1  (23)

where URi is the uncertainty ratio of a point source in subwatershed i.

The exchange ratio for a given point or non-point source in a watershed i is represented by the following equation:

ERi=DRi*URi*WRi*RRi  (24)

where ERi is the exchange ratio for subwatershed i, DRi is the delivery ratio for subwatershed i, URi is the uncertainty ratio for subwatershed i, WRi is the water quality ratio for subwatershed i, and RRi is the retirement ratio for subwatershed i. The water quality ratio is included to protect local water quality from adverse impacts of pollutant loading contributed by the subwatershed as it travels to the confluence point. The retirement ratio allows a certain percentage of the available credits to be retired or excluded from trading to increase the potential for water quality benefit or to provide a margin of environmental safety for the overall trading program.

Since the delivery ratio as calculated above does not take into consideration any variability, a delivery ratio correction factor may be used to adjust the delivery ratio. The correction factor is assumed to have a normal distribution with a unity mean and a certain coefficient of variation based on uncertainty in the delivery ratio. The adjusted exchange ratio is represented by the following equation:

Ê{circumflex over (R)}i=λ*ERi  (25)

where Ê{circumflex over (R)}i is the corrected exchange ratio for subwatershed i and λ is the correction factor.

The trading ratios are represented by the following equation:




  
    
      
        
          TR
          s
        
        =
        
          
            ER
            s
          
          
            ER
            b
          
        
      
    
    
      
        (
        26
        )
      
    
  




where TRs is the trading ratio for seller s, ERs is the exchange ratio for seller s, and ERb is the exchange ratio for buyer b.

The amount of credits earned by a subwatershed by applying a BMP is represented by the following equation:

{overscore (EC)}i=({overscore (L)}ipre−{overscore (L)}ipost)*Ê{circumflex over (R)}i  (27)

where {overscore (EC)}i is the mean earned credit for subwatershed i, {overscore (L)}ipre is the mean load of subwatershed i prior to applying the BMP, {overscore (L)}ipost is the mean load of subwatershed i after applying the BMP, and Ê{circumflex over (R)}i is the corrected exchange ratio for subwatershed i.

The variance of the earned credit is represented by the following equation:

Var({overscore (EC)}i)=(({overscore (L)}ipre*CVLipre)2+({overscore (L)}ipost*CVLipost)2)*Ê{circumflex over (R)}2  (28)

where Var({overscore (EC)}i) is the variance of mean earned credit for subwatershed i, {overscore (L)}ipre is the mean load of subwatershed i prior to applying the BMP, CVLipre is the coefficient of variation of subwatershed i prior to applying the BMP, {overscore (L)}ipost is the mean load of subwatershed i after applying the BMP, CVLipost is the coefficient of variation of subwatershed i after applying the BMP, and Ê{circumflex over (R)}i is the corrected exchange ratio for subwatershed i.

The lower and upper confidence intervals of the earned credit corresponding to a certain significance level are represented by the following equation:

CIi={overscore (EC)}i[1±zα/2[Var(ECi)]1/2]  (29)

where CIi is the confidence interval of subwatershed i, {overscore (EC)}i is the mean earned credit for subwatershed i, zα/2 is the z value at significance level α/2, and Var(ECi) is the earned credit for subwatershed i.

FIGS. 10–14 illustrate the processing of the WQTS in one embodiment.
FIG. 10 is a block diagram illustrating the components of the WQTS in one embodiment. The WQTS includes a user interface component 1001, a GIS component 1002, a loading component 1003, a controls component 1004, a trading component 1005, and a database 1006. The user interface component controls the user interface of the WQTS. The GIS component provides a geographic information system interface of the WQTS. The loading component performs the loading calculations of the WQTS including the current loads and future loads. The controls component controls the input of the control information and performs the controls calculations of the WQTS including the local credits. The trading component performs the trading calculations of the WQTS including the earned credits. The database contains the GIS information and the other information describing BMPs, coefficients of variation, and so on.
FIG. 11 is a flow diagram illustrating the processing of a component that identifies subwatersheds in one embodiment. After a river basin and a confluence point are selected, the component identifies the subwatersheds that are upstream from the confluence point. The illustrated processing is recursive. One skilled in the art will appreciate that non-recursive processing may alternatively be performed. The component is passed an identification of the confluence point. In blocks 1101–1104, the component loops selecting each subwatershed that directly contributes to the confluence point and recursively identifies each subwatershed that indirectly contributes to the confluence point. In block 1101, the component selects the next subwatershed that directly contributes to the passed confluence point. In decision block 1102, if all the subwatersheds have already been selected, then the component returns, else the component continues at block 1103. In block 1103, the component adds the selected subwatershed to a list of subwatersheds that are upstream from the original confluence point. In block 1104, the component recursively invokes the component passing the inflow of the selected subwatershed as a confluence point to identify the subwatersheds that contribute to the selected subwatershed. The component then loops to block 1101 to select the next subwatershed that directly contributes to the passed confluence point.
FIG. 12 is a flow diagram illustrating the processing of a component that calculates the loading for the subwatersheds that are upstream from a confluence point in one embodiment. The component calculates the current and future loads for the point and non-point sources and calculates the future loads of the non-point sources at various confidence levels. In block 1201, the component selects the next subwatershed. In decision block 1202, if all the subwatersheds have already been selected, then the component completes, else the component continues at block 1203. In block 1203, the component calculates the current and future loads for the non-point sources of the selected subwatershed as illustrated by equations (1) and (2). In block 1204, the component retrieves the future load of the selected subwatershed from the database. In block 1205, the component calculates the future load of the non-point sources of the selected subwatershed at various confidence levels as illustrated by equations (3) and (4). The component then loops to block 1201 to select the next subwatershed.
FIG. 13 is a flow diagram posted in the processing of a component that calculates the local credits for the subwatersheds that are upstream from a confluence point in one embodiment. The component loops selecting each subwatershed and calculating the local credits for point sources and non-point sources of that subwatershed. In block 1301, the component selects the next subwatershed. In decision block 1302, if all the subwatersheds have already been selected, then the component returns, else the component continues at block 1303. In block 1303, the component calculates the load of the point sources contributing to the local credits of the selected subwatershed. In decision blocks 1304–1306, the component loops calculating the load of each land use contributing to the local credits of the selected subwatershed. In block 1304, the component selects the next land use of the selected subwatershed. In decision block 1305, if all the land uses of the selected subwatershed have already been selected, then the component continues at block 1301 to select the next subwatershed, else the component continues at block 1306. In block 1306, the component evaluates the BMPs of the selected land use. The component calculates the reduction in load attributable to the serial BMPs according to equations (10) and (12) and the parallel BMPs according to equations (10), (14A), and (14B) to identify the portion of the load contributed by that land use. The component then loops to block 1304 to select the next land use. The component adds the load for each land use to determine a total non-point source load. The component then adds the load for the point sources and the non-point sources and subtracts the sum from the future load calculated without the controls to determine the local credits for the selected subwatershed.
FIG. 14 is a flow diagram illustrating the processing of a component that calculates the earned credits for subwatersheds that are upstream from a confluence point in one embodiment. In block 1401, the component selects the next subwatershed. In decision block 1402, if all the subwatersheds have already been selected, then the component returns, else the component continues at block 1403. In block 1403, the component calculates the coefficient of variation of the load for the selected subwatershed according to equation (18). In block 1404, the component calculates the uncertainty ratio for the selected subwatershed according to equation (22). In block 1405, the component calculates the exchange ratio for the selected subwatershed according to equations (24) and (25). In block 1406, the component calculates the earned credits for the selected subwatershed according to equation (27). The component then loops to block 1401 to select the next subwatershed.
From the foregoing, it will be appreciated that specific embodiments of the invention have been described herein for purposes of illustration, but that various modifications may be made without deviating from the spirit and scope of the invention. One skilled in the art will appreciate that the WQTS can be used to evaluate for various types of pollutants (e.g., phosphorus and sediment) and can be used at various levels of granularity (e.g., subwatersheds, watersheds, and point and non-point sources). Accordingly, the invention is not limited except as by the appended claims.

## Claims

1. A method in a computer system for determining a pollutant load of a source at a selected confidence level, the method comprising:
receiving an indication of a control to be applied to the source, the control having a removal efficiency;
calculating an uncertainty ratio for the source at the confidence level based on coefficients of variation of an export coefficient and the removal efficiency, the uncertainty ratio indicating the portion of the load that will be generated by the source with the confidence level; and
calculating a determined load for the source based on the load and the calculated uncertainty ratio.

2. The method of claim 1 wherein the source is a subwatershed that has multiple non-point sources of the pollutant and the control is applied to a specific non-point source.

3. The method of claim 2 wherein the controlled non-point source is a land use and the control is a best management practice.

4. The method of claim 1 wherein the calculating of the determined load is further based on a delivery ratio to a confluence point.

5. The method of claim 4 wherein the delivery ratio factors in local delivery ratios of sources downstream from the source.

6. The method of claim 4 wherein the calculating of the determined load is based on a water quality ratio and a retirement ratio.

7. The method of claim 1 including calculating a trading ratio based on uncertainty ratios of two sources.

8. The method of claim 1 wherein the uncertainty ratio is calculated by the following equation:




  
    UR
    i
  
  =
  
    
      
        
          L
          _
        
        i
      
      
        
          
            L
            _
          
          i
        
        ⁡
        
          [
          
            1
            +
            
              
                z
                α
              
              *
              
                CVL
                i
              
            
          
          ]
        
      
    
    =
    
      1
      
        1
        +
        
          2
          *
          
            CVL
            i
          
        
      
    
  




where URi is the uncertainty ratio for subwatershed i, {overscore (L)}i is the mean load of subwatershed i, CVLi is the coefficient of variation of the load of subwatershed i, and zα is the z value at a certain confidence level.

9. The method of claim 8 wherein the coefficient of variation of the load is calculated by the following equation:




  
    CVL
    i
  
  =
  
    
      
        (
        
          Var
          ⁡
          
            (
            
              L
              i
            
            )
          
        
        )
      
      
        1
        /
        2
      
    
    
      
        L
        _
      
      i
    
  




where CVLi is the coefficient of variation of the load of subwatershed i, Var(Li) is the variance of the pollutant load of subwatershed i, and {overscore (L)}i is the mean load of subwatershed i.

10. The method of claim 8 wherein the coefficient of variation of the removal efficiency for controls in series is represented by the following equation:




  
    CV
    ⁢
    
      
    
    ⁢
    
      η
      ij
    
  
  =
  
    
      ∏
      
        k
        =
        1
      
      
        n
        ij
      
    
    ⁢
    
      
    
    ⁢
    
      
        (
        
          1
          +
          
            
              Var
              ⁡
              
                (
                
                  η
                  jk
                
                )
              
            
            
              
                (
                
                  1
                  -
                  
                    
                      η
                      _
                    
                    jk
                  
                
                )
              
              2
            
          
        
        )
      
      
        1
        2
      
    
  




where CVηij is the coefficient of variation of the removal efficiency of land use j in subwatershed i, Var(ηjk) is the variance of the removal efficiency for BMP k in land use j, {overscore (η)}jk is the mean removal efficiency for BMP k in land use j, and nij is the number of BMPs of land use j in subwatershed i.

11. A method in a computer system for determining a load of a pollutant of a source within a confidence level, the source having a load, the method comprising:
calculating an uncertainty ratio for the source at the confidence level based on a coefficient of variation representing an uncertainty in the load of the source, the uncertainty ratio indicating a portion of the load that will be generated by the source with the confidence level; and
calculating a determined load for the source based on the load and the calculated uncertainty ratio.

12. The method of claim 11 wherein the coefficient of variation is based on an export coefficient.

13. The method of claim 12 wherein the export coefficient is based on a land use of the source.

14. The method of claim 12 wherein the calculating of the determined load is further based on a delivery ratio to a confluence point.

15. The method of claim 14 wherein the delivery ratio factors in local delivery ratios of sources downstream from the source.

16. The method of claim 14 wherein the calculating of the determined load is based on a water quality ratio and a retirement ratio.

17. The method of claim 12 including calculating a trading ratio based on uncertainty ratios of two sources.

18. The method of claim 12 wherein the uncertainty ratio is calculated by the following equation:




  
    UR
    i
  
  =
  
    
      
        
          L
          _
        
        i
      
      
        
          
            L
            _
          
          i
        
        ⁡
        
          [
          
            1
            +
            
              
                z
                α
              
              *
              
                CVL
                i
              
            
          
          ]
        
      
    
    =
    
      1
      
        1
        +
        
          2
          *
          
            CVL
            i
          
        
      
    
  




where URi is the uncertainty ratio for subwatershed i, {overscore (L)}i is the mean load of subwatershed i, CVLi is the coefficient of variation of the load of subwatershed i, and zα is the z value at a certain confidence level.

19. The method of claim 18 wherein the coefficient of variation of the load is calculated by the following equation:




  
    CVL
    i
  
  =
  
    
      
        (
        
          Var
          ⁡
          
            (
            
              L
              i
            
            )
          
        
        )
      
      
        1
        /
        2
      
    
    
      
        L
        _
      
      i
    
  




where CVLi is the coefficient of variation of the load of subwatershed i, Var(Li) is the variance of the pollutant load of subwatershed i, and {overscore (L)}i is the mean load of subwatershed i.

20. The method of claim 11 including receiving an indication of a control to be applied to the source, the control having a removal efficiency, wherein the coefficient of variation is based on the removal efficiency of the control.

21. The method of claim 20 wherein the source is a subwatershed that has multiple non-point sources of the pollutant and the control is applied to a non-point source.

22. The method of claim 21 wherein the controlled non-point source is a land use and the control is a best management practice.

23. The method of claim 20 wherein the calculating of the determined load is further based on a delivery ratio to a confluence point.

24. The method of claim 23 wherein the delivery ratio factors in local delivery ratios of sources downstream from the source.

25. The method of claim 23 wherein the calculating of the determined load is based on a water quality ratio and a retirement ratio.

26. The method of claim 20 including calculating a trading ratio based on uncertainty ratios of two sources.

27. The method of claim 20 wherein the uncertainty ratio is calculated by the following equation:




  
    UR
    i
  
  =
  
    
      
        
          L
          _
        
        i
      
      
        
          
            L
            _
          
          i
        
        ⁡
        
          [
          
            1
            +
            
              
                z
                α
              
              *
              
                CVL
                i
              
            
          
          ]
        
      
    
    =
    
      1
      
        1
        +
        
          2
          *
          
            CVL
            i
          
        
      
    
  




where URi is the uncertainty ratio for subwatershed i, {overscore (L)}i is the mean load of subwatershed i, CVLi is the coefficient of variation of the load of subwatershed i, and zα is the z value at a certain confidence level.

28. The method of claim 20 wherein the coefficient of variation of the load is calculated by the following equation:




  
    CVL
    i
  
  =
  
    
      
        (
        
          Var
          ⁡
          
            (
            
              L
              i
            
            )
          
        
        )
      
      
        1
        /
        2
      
    
    
      
        L
        _
      
      i
    
  




where CVLi is the coefficient of variation of the load of subwatershed i, Var(Li) is the variance of the pollutant load of subwatershed i, and {overscore (L)}i is the mean load of subwatershed i.

29. The method of claim 20 wherein the coefficient of variation of the removal efficiency for controls in series is represented by the following equation:




  
    CV
    ⁢
    
      
    
    ⁢
    
      η
      ij
    
  
  =
  
    
      ∏
      
        k
        =
        1
      
      
        n
        ij
      
    
    ⁢
    
      
    
    ⁢
    
      
        (
        
          1
          +
          
            
              Var
              ⁡
              
                (
                
                  η
                  jk
                
                )
              
            
            
              
                (
                
                  1
                  -
                  
                    
                      η
                      _
                    
                    jk
                  
                
                )
              
              2
            
          
        
        )
      
      
        1
        2
      
    
  




where CVηij is the coefficient of variation of the removal efficiency of land use j in subwatershed i, Var(ηjk) is the variance of the removal efficiency for BMP k in land use j, ηjk is the mean removal efficiency for BMP k in land use j, and nij is the number of BMPs of land use j in subwatershed i.