UNITED STATES ENVIRONMENTAL PROTECTION AGENCY
                           WASHINGTON, DiC. 20460
                                        "'
                                   281984
                                                                 OFFICE OF
                                                                  WATER
 MEMORANDUM
 SUBJECT:  Technical Guidance Manual for Performing Waste Load
           Allocations Book VII, Permit Averaging Periods

 IDs       Regional Water Management Division Directors
           Regional Environmental Services Division Directors
           Regional Wasteload Allocation Coordinators                      .
 •        : • '       ''.     -      '      -..-••          -         .  • .      •  .'
      Attached, for national use, is the final version  of the  Technical
 Guidance.Manual for Performing Waste Load Allocations, Book VTI, Permit
 Averaging Periods;  We are sending extra copies of this manual  to  the
 Regional Wasteload Allocation Coordinators for distribution to  the States
 to use in conducting waste load allocations. ,.
           '  •                  =                    '  '
 '                             '  -         '
      Modifications to the February 1984 draft include:                     "

      o  The method to calculate, the Reductions7Factor  in Chapter 2 has been
      ,   elaborated to include the use of 95% cut-offs  for frequency of permit
         violations.  • •  •     '  ..   .
-       .     .        ..''"'     '     ,  •  '  .     , ,        ..'•••    .   •_•('."•_
      o  The example calculation in Chapter 3 has been  expanded.. Step 1 has
  ,       been added to the step-procedure to show how permit limits can. be
         specified using 95% cut-offs for frequency of  permit  violations.

      o  The document recommends that advanced treatment facilities should  be
         built to meet the long-term average and the selected  effluent
       -  variability.    .        '..''•"                          ',     '

      o  A flow diagram and an IBM PC-compatible program have  been  added to
         Appendix D.

      If you have any questions or comments or desire additional information
 please contact Tim S. Stuart, Chief, Monitoring Branch, Monitoring and Data
 Support Division (WH-553) on (FTS) 382-7074,       .     '                  .;
                                   Edwin L.  Johnson,  Director
                                   Office of Water Regulations
                                    and Standards (WH-551)
 Attachment

-------

-------
  Technical Guidance Manual for
Performing Waste Load Allocations

Book VII: Permit Averaging Periods
              September 1984
                Final Report

                 •.  .for
                                  **
      Office of Water Regulations and Standards
       Monitoring and Data Support Division,
              Monitoring Branch
       U.S. Environmental Protection Agency
      401 M Street, S.W. Washington, D.C. 20460

-------

-------
                                         VII (0)
                                         Revision No. 0
       TECHNICAL  GUIDANCE MANUAL FOR-.

     PERFORMING WASTE  LOAD ALLOCATIONS


                 Book  VII

         Permit Averaging Periods
     Contract Number S8-03-3131-WA9



             Project Officer

             Hiranmay Biswas
Office of Water Regulations and Standards
  Monitoring and Data .Support Division
            Monitoring Branch


  U.S. Environmental  Protection Agency
          401 M Street, S.W.    '     '.
        Washington,  D.C. 20460
            September 1984  .

-------
_l

-------
                                 FOREWORD

      - This guidance document is a product of several years of research
on many complex water quality issues.  Although much progress has been made,
some  issues still remain.  User participation will be needed to develop
answers to these unresolved issues and will be key to future revisions of
this  document.
       Selection of permit averaging periods, as presented in this manual, is
based on an assumed-,exceedance frequency of an acute violation in the stream
no more than 1 day in 10 years-;  The EPA is currently considering the issue
of allowable duration and frequency of exposure to acute as well  as chronic
toxicity.  Based on this study, the choice of duration and frequency used in
this document as examples may have to be changed.

-------
                                                          VII  (0)
                                                          Revision  No.  0
                                 CONTENTS              .


Chapter      •

        FOREWORD.	.  .  .  ....'.  .  .  .  ,  . i

        LIST OF TABLES. . ^ ..............  .  . iv

        LIST OF FIGURES '. .	v-vi

        LIST OF ABBREVIATIONS AND SYMBOLS  ......... vii-viii

        ACKNOWLEDGMENTS ...............  /*•  .  .  , ix

        EXECUTIVE SUMMARY	  .  a ..'.'*......!

  1     INTRODUCTION                                         •

        1.1   Background.  . . .	 1-1
        1.2   Objectives.	 1-2
        1.3   Approach	 1-3
        1.4   Organization. . . . .	 1-5

  2     METHODOLOGY                                       /

     1   2.1-   Description  of Probabilistic Dilution Model  . 2-1
        2.2   Choice of the Permit  Averaging- Period .  .  .  .2-9

  3    .EXAMPLE COMPUTATION

        3.1   Hypothetical  Site-Specific Conditions .... 3-2
        3.2   Example Computation - Hand Calculation.  .  .  . 3-8
        3.3   Example Computation - Computer Program.  .  .  . 3-28

  4     RANGE  OF EXPECTED  VALUES FOR STREAMS IN U.S.                .

        4.1  Analysis for Conservative Substances. .... 4-1
        4.2   Use  as a Screening Tool	 4-9
        4.3   Preliminary  Analysis  for Dissolved Oxygen .  . 4-10
        4.4  Analysis for Conservative Substances  in
             Effluent-Dominated Streams	 4-22

  5     USES AMD LIMITATIONS.  ....... e  ....... 5-1

  6    REFERENCES.  . .  ...  ........  ..'. .. . . , ..6-1
                                     11

-------
                                                          VII  (0)
                                                          Revision  No.  0
                               ;  , CONTENTS-   .          = .    ;

                1     ,            (Continued)
                                     \             •        .         '

Appendix                                                  .       page


   A  .STATISTICAL PROPERTIES OF LOG-NORMAL DISTRIBUTIONS

     .' A-l.  General  Considerations. . . ....,: . v  ........ A-l
       A-2.  Probability Distributions . . . '. .  . . . . . . . . A-3
       A-3.  Relationship Between .Distributions*  . . ... . . . A-6
       A-4.  Properties of Log-Normal Distributions	 . A-6
       A-5.  Standard Normal  Tables. . . ....  . . . . j . . . A-10
       A-6.  References. . ....  .  . . . . . ."-..  . . . .. . . . . A-12

   B   FIELD VALIDATION OF LOG-NORMAL .DISTRIBUTION AND
       RELATED ASSUMPTIONS                                           .

'.'      8-1.  Use of the Log-Normal  Distribution.  ........ B-l
       B-2.  Verification of  the  Probabilistic Dilution  Model.. B-3
       B-3.  Appropriateness  of Assumptions	 B-8
       B-4.  References. . .  ...  . . . .  . ......... .- B-ll

   C   CHARACTERISTIC VALUES  FOR  INPUT PARAMETERS            '

       C-l.  Treatment  Plant  Effluent Flows. .......... C-l
       C-2.  Treatment  Plant  Effluent Concentrations  ...... C-l
       C-3.  Stream Flow	 C-3
       C-4.  References. . ...'....'•....'......'.,.  . C-6

   D   COMPUTER PROGRAM FOR THE PROBABILISTIC  DILUTION
       MODEL -  POINT  SOURCE (PDM-PS)

    •.  D-l.  Formulation and  Normalization  .- .  .  . .  ......  D-l
       D-2v  Description of Program  Use.  .  ...........  D-5
                                    iii

-------
                                                         VII (0)
                                                         Revision No. 0
                            .  LIST OF TABLES
Table   .               -                                   .               page

 2-1  Reduction factors for various coefficients of variation. . . .  «, . 2-13
 4-1  Averaging period selection matrix for conservative substances:
      effluent dilution ratio  - 7Q10/OJ =50 	 ......... 4-5
 4-2  Averaging period selection matrix for conservative substances:
      effluent dilution ratio  - 7Q10/OT = 5. .............. 4-6
 4-3  Averaging period selection matrix for conservative substances:
      effluent dilution ratio  - 7Q10/TJT = 3. . ...	, .... 4-7
 4-4  Averaging period selection matrix for conservative substances:
      effluent dilution ratio  - 7Q10/TJF = 1.	4-8
 4-5  Conditional.moments for  the low flow sub-population.  ...... .4-16
 4-6  Permit averaging period  selection matrix for BOD/DO:   effluent
      dilution ratio - 7Q10/OJ' = 5	 4-19
 4-7  Permit averaging peri Deselection matrix for BOD/DO:   effluent
      dilution ratio - 7Q10/QE = 3 ..................... 4-20
 4-8  Permit averaging period  selection matrix for BOD/DO:   effluent
      dilution ratio - 7Q10/QE" =, 1 ..... s .............. 4-21
 4-9  Averaging period selection matrix for effluent  dominated streams  . 4-25
 A-l  Probabilities  for the  standard  normal  distribution  ....  .....  .  . A-ll
 B-l  Comparison of  observed and computed  downstream  concentrations.  .  . B-6
 B-2  Approximate  overestimation of 10 year return period stream
      concentration  by ignoring  serial  correlation-.  	  ....  8-10
 C-l  Coefficient  of variation of  daily effluent  flows, vg£.  ...... c-2
 C-2  Summary of secondary treatment  plant performance - median
      coefficients of variation, vr,E	C-4
 C-3  Effluent concentration variability  for trickling filters  ....  .-C-6
 C-4  Summary of stream flow characteristics  .............  a  C-8
                                    IV

-------
                                                         VII (0)
                                                         Revision No. 0
                              LIST OF FIGURES

Figure     .             •              .                                   Page


 1-1  Schematic-outline of probabilistic method. . .	 1-4 ^

 2-1  Simple dilution model.	 . 2-2

 2-2  Illustration of analysis results:  stream concentration versus     .
      return period for three permit averaging periods . .	2-16

 3-1 "'Step procedure to select optimal  permit averaging period . . . . . 3-1

 3-2  Sample stream concentration versus probability plot for 30-day
      averaging period .	 . ... 3-15

 3-3  Sample stream concentration versus mean recurrence interval
      for 30-day averaging period. . ,...•.....•...•.»	 . 3-17

 3-4  Concentration versus probability  plot for 1-,  7-, and 30-day
      averaging periods, i ......,...,'........... .3-19

 3-5  Concentration versus-mean  recurrence interval  plot for IT,   •
      7-,  and 30-day averaging periods  . . . ............... 3-19

 3-6 • Concentration versus probability  for PDM-PS computation. .  .- . . .3-33

 3-7  Concentration versus mean  recurrence interval  for PDM-PS
      computation.  ...........*..,,.  ...  . .  ...  .  . . . 3-33

 4-1  Effect  of permit averaging period  on stream concentrations  for
      conservative  substances:   general  analysis  .	  4-4

 4-2 .Effect  of permit averaging period  on stream concentrations  for
      BOD/DO  ....... .,..  . . .  .  . . . ...-'.  .  .  .  . ,.  ...... .4-18

 4-3  Effect  of permit averaging period  on stream concentrations  for  .
      conservative  substances  in effluent-dominated  streams.  ......  4-24

 A-l  Probability distributions. .. .........'.•..'-.;.,..  .  A-5

 A-2  Effect  of coefficient  of variation on  frequency distribution  .  .  .  A-7

 A-3  Pertinent relationships  for  log-normal  distribution. -.  .  ....  .  A-8

 A-4  Cumulative log-normal  distribution .................  A-9

 B-l  Evaluations of  Tog-normaT  distribution  for  stream  flows..  ....  B-2

-------
                                                         VII. (0)
                                                         Revision  No.  0
                              LIST OF FIGURES
                                (continued)
Figure
 B-2  Probability distribution of treatment  plant  effluent  concen-
      trations - conventional  pollutants ................ B-4
 B-3  Probability distribution of treatment  plant  effluent  concen-
      trations - heavy metals	 8-5
 C-l  Typical  low flow characteristics  of  U.S.  streams  .......... C-7
                       /                            •                    .
 D-l  CRT displays	,	 D-9
 D-2  Example  of printed  output.  ...................  .0-10
 D-3  Flow chart for PDM-PS program.  .................. D-ll
 D-4  PDM-PS program listing - HP85 compatible  ............. D-12
 0-5  PDM-PS program listing - IBM-PC and MS-DOS compatible.,  .  	 D-15
                                    v1

-------
                                                          vii
                                                         .Revision Wo.  0
                      LIST OF  ABBREVIATIONS AND SYMBOLS,

 BASIC     Computer  language
 BOD      Biochemical  oxygen demand
 BODs     The  amount of dissolved  oxygen  consumed  in  five days by biological
          oxidation of organic matter         .
 CE        Treatment plant  effluent concentration
 CFS      Cubic  feet per second, unit of  flow
 CL        Concentration equal  to a wa.ter  quality standard
 CO        Downstream concentration, after complete mixing   ,
 CRT      Cathode ray  tube                                   .
 Csat      Saturation concentration of dissolved oxygen
 CS    '    Stream concentration upstream of  discharge
 D         Flow ratio,  equal to QS/QE
 DC        Critical  (or maximum) dissolved oxygen deficit
 DO        Dissolved oxygen
 EL        Effluent  limit.  A maximum effluent concentration determined from
          a waste load allocation  analysis, and specified by an NPDES
          permit
 FAV      Final acute  value
 FCV      Final chronic value                      '
 K         Stream purification factor.
 Ka        Stream reaeration rate constant
 Kd    .    BOD oxidation rate constant
MRI       Mean recurrence interval, expressed in years
 WPDES     National pollutant discharge elimination system
f         Pollutant
                                     vn

-------
                                                          VII (0)
                                                          .Revision No. 0
                      LIST OF ABBREVIATIONS AND SYMBOLS
                                (Continued)

 PDM-PS   Probabilistic dilution  model:   point source
 P01W     Publicly-owned treatment works
 Pr        Probability          .   '
 7Q10     The  lowest 7-day average stream flow with  a  recurrence  interval
          of 10 years
 Q£        Treatment pi ant  eff 1 uent flow
 QS        Stream  flow
 QT        Total downstream flow,  equal  to QS + QE
 R         Reduction factor, equal  to  the  ratio of the  mean  CE  for which  a
          treatment'pi ant  is  designed to  the EL
 TSS       Total suspended  solids
 WLA       Waste load allocation
 WQ        Water quality
 a         Exceedence probability                        .
 P         Dimensionless  unit  of concentration  equal to  CO/CL
 i*x        Mean value of  x
 «t»        Dilution factor
 °x        Standard deviation  of x
 vx        Coefficient of variation of x   .
Za        Value of statistical parameter  Z  for a probability of a

-------
                                                           VII  (0)
                                                           Revision No. 0
                   .             ACKNOWLEDGMENTS  .-,

         This report was developed by JACA Corporation under contract to the
  U.S. Environmental Protection Agency (Contract No. 68-03-3131), and is-the
  product of the contributions of the following individuals:

         Eugene D. Driscoll            (E.D. Driscoll and Associates)
         Dominic M. DiToro             (Hydro Quai, Inc.)

  Editorial assistance has been provided by:

         Patrick J. Rafferty           (JACA Corporation)
         Virginia R. Hathaway          (JACA Corporation)

       .  Hi ranmay Biswas, the EPA Project Officer,  USEPA Washington, -DC,
  provided  direction on the basic content and emphasis, and coordinated the
  review process with EPA.Regional  Offices.  Tim Stuart, Chief of the Monitor-
  ing Branch, and Mark Morris, Chief of the Waste!oad Allocation Section,
  provided  overall  guidance and support in developing this manual. :
         Water quality personnel  in EPA headquarters, the Regions, States  and
.  ORD laboratories  and the private  sector provided  review, comments and sugges-
  tions that contributed significantly to this  effort.  .
  Special thanks to:

         Charles App,  EPA - Region  III
         Thomas  Barnwell, EPA, ERL  - Athens,  GA
         John Hall, OWPO, EPA-HQ
         Norbert Huang, OWPO,  EPA-HQ

-------
                                                   VII (0)
                                                   Revision  No.  0
Henry Kahn,  EPA-HQ
Narendar M.  Kumar, DER,  State of  Florida
Noel Kohl, EPA-Region V
Maurice Owens, EPA-HQ
James J. McKeown, Regional' Manager, NCASI
Bruce Newton, OWEP, EPA-HQ
Lewis A. Rossman, EPA-MERL/Cincinnati, OH
Donald R. Schregardus, EPA-Region V
                                                                            V
                                                                             I

-------
                                                           VII  (0)
                                                           Revision No, 0
  '                            EXECUTIVE SUMMARY   '-'

 Background                                                              _   .

        The conventional approach to developing Waste Load Allocations
 (WLAs) is based on a steady state analysis qf.. stream conditions, using a
 design stream flow (usually the 7Q10) and a receiving water concentration   ,  '
 (usually a water quality standard based on chronic  criteria) for the
 pollutant to be allocated.  An  effluent concentration limit is computed for
 these conditions, and is used  to establish the .NPDES permit conditions.

        The water quality based, permit conditions  apply,  in  addition  to
 technology based requirements  (e.g.,  BAT,  BCT, and  secondary treatment).
 This  effluent requirement may be incorporated  into  the permit'as the daily   •
 maximum limit,  the average limit over a week (for POP^s)  or the  .average limit
 over  a month'(for industrial as  well  as municipal source) .1, Typical
 practice for toxic pollutants is to incorporate the  wasteload  allocation
 result as  the daily maximum  permit limit.   This document  provides an  inno-
 vative approach to  determining which  types  of  permit  limits  (dai-ly maximum,
 weekly, or monthly  average)  should be specified for the steady-state.model
 output based  on  the frequency of acute  criteria violations.

 Approach

        The method used to evaluate the. effect of, permit averaging periods'
                                                           •T '."• .••'•'•
 is based on a probabilistic dilution model  (PDM)  in which it is assumed that
the stream flows, effluent flows and concentration are Iqg-pnormally distributed
     40 CSR 122.45(d).
                                     -1-

-------
                                                          VII (0)
                                                          Revision No. 0
 and tin correlated. The log-normal distribution is known to be representative
 of effluent behavior and to almost always under-estimate the lowest stream
 flows somewhat.  Thus, the analysis is generally conservative (overprotec-
 tlve) to some extent.  However, a verification of the probabilistic dilution
 model Indicates that, for the cases tested, it correctly estimates observed
 downstream concentration probability distributions-to within the confidence
 limits of the data.                    .                •     \    '•

        The method applied in using  this model  to evaluate permit averaging
 period choices-is based  on the following  observation.  If chronic  criteria
 and 7-day, 10-year Tow flow, or any other state-specified low flow,  are used
 in the WLA analysis  to develop the  maximum effluent  concentration, the use
 Of monthly or  weekly permit limits  for  specifying this  effluent  requirement
 presents the possibility  that  simultaneous occurrences  of high effluent
 concentrations and low stream  flows may result in stream  concentrations
 which  exceed the  acute criteria  for a pollutant  without violating maximum
 average discharge permit  conditions..
 •      •                                                               •••
       The  analysis  consists of  computing  the level of• treatment required for
 the three averaging  period  options for specifying the WLA results as permit
 limits.  The analysis computes the frequency at which acute stream criteria
 concentrations are violated under each of the permit averaging period options,
taking Into account the likely range of stream and effluent variability.
 Computation results are normalized so that summary results can be applied  to
a variety of pollutants based on their ratio of acute-to-chronie  criteria
concentrations.
                                     -2-

-------
                                                          VII (0)
                                                          Revision No. 0
 Uses
        The'primary use of this methodology will  be specifying the  required
 level of treatment and deriving permit limits based ,on  water quality  require-
 ments.  Care must be taken in the  assumptions related to  the permit limits
 and assumptions used in the methodology.   For example,  throughout  this
 document, reference is made to 7-day  and  30-day  averages.  These averages
 are equivalent to weekly and monthly  permit limits where  the assumption can
 be made that the monitoring data is adequate (i.e., that  the data  collected
 In a month adequately reflects the SCMay  average). Where  this requirement
 is not valid,  alternative  limits may  be calculated which  incorporate monitoring
 frequency, or  monitoring frequency may be  adjusted so that  these conditions
 are met.  ,
                                          ''
        In  addition  to the  usefulness  of this method for permit writers in
 selecting  the  averaging period  for discharge  permits, the method has been
 used to calculate suitable  averaging,  periods  for the range of stream and
 effluent conditions typified  in  the U.S.   The results have been summarized in
 convenient  graphic and tabular displays, and can be. used as a "screening
tool" that  provides, a guide for water quality decisions.  These summaries   •
 show,  for  instance, that for toxic pollutants with acute-to-chrohic ratios  of
 10 or  greater, 30-day permit averages will virtually always meet the criteria
that  have been adopted; that is, that acute criteria violations in  the stream
will recur with a frequency that averages  less than 1  day in 10 years.1
1
 The EPA is presently considering the issue of allowable duration and
 frequency of exposure to toxicity.  Based upon this workv duration and
 frequencies used as the decision criteria may change.   This  guidance does
 not recommend any particular minimum acceptable duration or  frequency.
                                    -3-

-------
                                                          vii (9)
                                                          Revision No. 0
 For pollutants with acute-to-chronic ratios of between 5 and 10, monthly
 permit averages will be appropriate in most cases, although there will be
' some site-specific conditions that would call for the use of weekly averages.
 For pollutants with acute-to-chronic ratios of less than 5,site specific
 conditions must be considered, and no general rule is possible.  In these
 cases, site-specific analyses of the effects of different permit averaging
 periods can be performed using the methods outlined in the text.

 Limitations                         ,
        Several  technical  refinements to the probabilistic  model would  be
 required  to  more accurately  reflect  the deviation  of  lowest  stream  flow
 from  log-normality,  and to ac'count for  serial  and  cross-correlation, of
 stream  flows and effluent loads.  For  coupled reactions,  such as BOD/DO,
 the procedures  would have to  be extended to  provide a seasonal approach
 and results  should be verified against  field data.  .The analysis
jnethod  would have to be extended to  incorporate the variability of  sec-
 ondary  water quality parameters such  as pH,  hardness and temperature,
 since these  affect the toxicity of a  number  of pollutants.  Finally, the
 chronic exposure  event, as defined by the state design.flow conditions, was
used throughout the  document to estimate the maximum effluent concentration.
Further analyses  to  determine the possible underprotection or overprotection
of chronic criteria  based on the state design flow* were not done.
 The EPA is considering studying the impact of uncertainties  involving
 the low flow estimating techniques on the selection of stream design  flow.

-------
                                                          Revision No. 0
                                  CHAPTER 1

                          ;      INTRODUCTION                     v



 1.1    Background              "                            ,



        The conventional procedure for establishing a point source  effluent

 limit using a waste load allocation (WLA) analysis begins  by  specifying  a

 target concentration of the pollutant in the stream, such  as  a  state  water

 quality standard based oil chronic criteria.   This  stream concentration is

 converted to a maximum effluent concentration using a mass balance calcula-

,tion (for conservative substances)  or. a  steady-state analysis (for reactive

 substances).  The inputs to these analyses are a design stream  flow (repre-

 senting  low stream-flow conditions)1  and  a measure of the  effluent flow,

 typically.the mean  effluent flow.  Although  this technique  is presumed to
 • *'                '            /                        '                 '
 provide  adequate protection for receiving water quality, it fails to account.

 for  random and other  fluctuations in  the  flow rate  and concentration that

 naturally occur in  both  the stream  and effluent.   Thus, the degree to which a

 given  limit protects  against exceedances of  acutely toxic concentrations
                                    ' .  - -   (  -   -     . •        • '  . ' .
 Is not quantified.    ,            ,



       Effluent  permit limitations are currently specified  as maximum  concen-

 trations  for one day  or  averaged  over a week or month.  The number of  obser-

 vations from .which the average is computed depends on the frequency of    ' ; ••,
         '                  '                    - '                   •
 The design stream flow most commonly used is the  7Q10 flow, which
 represents the low-flow condition with a recurrence interval' of 10 years
 cased on a 7-day averaging period.  Other flows,  suefi as the 30Q10 or
 Jugs are occasionally used as the design stream flow.  Wherever the use
 or stream design flow is called for, these or other stream design flows
 can be substituted throughout this document.
                                     i-i

-------
                                                         VII (1)
                                                         Revision No. 0
monitoring.  Although there is no generally accepted rational basis for
selecting permit averaging periods, the effluent requirement derived from a
      •                   ••*-'.
WLA is typically expressed as a monthly average for conventional pollutants
and as the daily maximum for toxic pollutants.  A set of conversion factors
is then'used to convert these concentrations to other averaging periods.
                                                 •»•                >
In this document the maximum daily, weekly, and monthly permit limits are
referred to as 1-day, 7-day, and 30-day permit levels, respectively.

       .The permit limit used to incorporate a WLA effluent requirement can
have a substantial influence on the degree (and cost) of treatment re-
quired and on the quality of the receiving water.  It is clear that a permit
limit imposed as a daily maximum requirement is more restrictive than when
the same permit limit is used as a 30-day average requirement, since in the
latter case the effluent concentration can fluctuate above"the effluent limit
for days at a time and still meet the 30-day average requirement.  Such
fluctuations may or may not be significant in terms of receiving water
quality.  The appropriate choice of the averaging period, then, is one which
ensures acceptable receiving water quality without imposing unnecessarily
restrictive treatment requirements.

1,2    Objectives

       This guidance document is intended to achieve the following:
       (1)  Present a rational method for selecting the level of treatment
            required based on considerations of water quality;
       (2)  Present a rational method to incorporate the water quality based
            treatment requirements as permit-limits;   •
                                     1-2

-------
                                                           VII  (1)
                                                           Revision Wo.  0
         (3)   Provide specific'Information, including detailed examples, so

              that the method can be applied to site-specific cases;

:   ;,     (4)   Use the method to provide an overall analysis of a broad
                         /        •    .           •          "         •' ' -
             'range of conditions likely to be .encountered, so as to provide

              a'screening tool for the rapid'assessment of a wide variety of

              cases;

         (5)   Discuss the uses and limitations of the method.



 1.3    Approach                v  .'.••'.



        The basis of the method is  an evaluation  of the extent  and frequency

 of acute criteria violations to  be expected  in the stream receiving the

 discharge as a  result  of imposing  the  effluent concentration,  computed  from a
                    •                   r               "         * '       '    '
 steady state wasteload allocation, as  a  daily, weekly, or monthly average

 permit.  A.probabilistic framework 'is  adopted to  account  for the  inherent

 variability of  flows and concentrations.   Acute criteria  violations are

 assumed to be associated with  random simultaneous  occurrences of  high

 effluent loadings  and  low stream flows.1   The analysis  is based on  an

 examination  of  the probability distributions  involved  and how they  combine  to

 Influence the concentration  downstream.   The  probabilistic dilution model

 provides the analysis  framework.  •   .•  "  -              •       •



        The probabilistic  dilution model is summarized  in Figure 1-1.  The
                ' .                                     '                ' " • •  .- i

 Inputs  to the model include the flow and concentration histories (or pro-

jections) of both the effluent and the receiving stream.  Each of these is
 While it is apparent that effluent loadings and stream flows experience '
 both random and nonrandom (e.g-., seasonal) variations, the problem is
 analyzed here in purely random terms to, limit the complexity of the analysis
                                     1-3

-------
                                                     VII  (1),
                                                     Revision No,
                                                                       ft
                                                                       as
                                                                       ui
                                                                       u
                                                                       o
                                                                       UI
                                                                       03

                                                                       CE
                                                                       o
                                                                       OT
                                                                       Uj
                                                                       o

                                                                       OT
                                                                      o
                                                                      H
                                                                      a.

                                                                      o
                                                                      u
                                                                      CO
                                                                      til
                                                                      O
                                                                      5
                                                                      z

                                                                     o
                                                                      uj
                                                                      o
                     A/                           —V
                    MOId                   NO!lVdiN30NOO
Figure  1-1 - Schematic outline of probabilistic method.
                               1-4

-------
                                                           VII (1)   ,
                                                           Revision No. 0
  expressed as a probability distribution; that is, in terms of the probability
  that a given value is. exceeded.  Next, the effluent and stream  flows are  '.
  combined, to yield the probability distribution of the dilution factor; then
  the dilution factor and concentrations are combined to provide, the  probability
  distribution for the resulting stream concentration.  The  stream  concentra^
.  tion probability distribution is  then converted  to  a plot  showing the.re-
  currence interval  to be associated  with each  stream concentration so  that the
  frequency  of occurrence of a  given  (high)  stream concentration  can  be com-
  pared to water  quality  objectives.  '                          •

         The probabilistic dilution model is used  to  guide the choice of the
  permit  averaging period as  follows.   Given an effluent  requirement  from a
  WLA analysis, the mean effluent required to meet that WLA requirement is
  calculated for each of the  three averaging periods-, based on an assumed
  allowable frequency of effluent limit violation.  This provides three levels
 of treatment for the plant  in question.  Each-mean effluent concentration is
 then used, together with the parameters that characterize the stream varia-
 bility, in the probabilistic dilution model.  The result is a probability
 distribution of resulting stream concentration for each of  the  three treat-
 ment plant options, which can  be compared to  daily concentration/frequency
 water quality goals.  The use  of daily concentration frequencies allows the
 use of acute  criteria in establishing  water quality  goals.

 1-4    Organization
                                                                        -•.   •   r
        This document  is  organized  as follows.  Chapter  2 provides a
detailed description  of the  methodology  for finding  an optimum averaging
                                      1-5

-------
                                                         VII (1)
                                                         Revision No. 0
period based on a probabilistic dilution method.  Chapter 3 presents an
annotated example of the .method performed first as a hand calculation and
then using the computer program provided in Appendix D. " Chapter 4 uses
the model in several representative applications, and Chapter 5 discusses
the uses of the method.  Several appendices .to this document, provide
detailed additional  material, including a review of relationships for
log-normal distributions (Appendix A) and a discussion of technical  issues
and assumptions employed in the analysis (Appendix B).
                                 1-6

-------
                                                          VII  (2)  .
                                                          Revision No. 0
                           .     CHAPTER' 2
';-.'"                    METHOD OF ANALYSIS
                       ' *                '        ,            '
       This  chapter  lays the theoretical groundwork for.the.application of
the  probabilistic dilution model to the problem of permit averaging period
selection. /This discussion is  presented in two parts.  Section 2.1
describes the probabilistic dilution model.  Section 2.2 develops the
method whereby the probabilistic dilution model is employed to predict the
water quality effects of the selection of different averaging periods.   x

2.1    Description of the Probabilistic Dilution Model

       The probabilistic dilution model is based on a simple stream
dilution calculation.  The complexity of the model arises from the proba-
bilistic framework that is superimposed upon the dilution equation..  This
section is intended to provide a description of the derivation of the
model, and to reduce it to a manageable set of equations.  While a strict
mathematical derivation of the model is available [1],  a rigorous treat-
ment is considered beyond the scope of this manual .

       Figure 2-1 illustrates a treatment plant discharge entering a
stream. The effluent discharge flow (QE), having a concentration ('CE) of
the pollutant of interest, mixes with the stream flow (QS),  which may have
a background concentration (CS).  The receiving water concentration (CO)
is the concentration that results after complete mixing of the effluent
and stream flows.  It is  the  cross-sectional  average concentration down-
stream of the discharge,  and  is given by:
                                     Z-l

-------
                                                 VII  (2)
                                                 Revision No. 0
  •r-EFFLUENT
  \  FLOW=QE
    \ CONCENTRATION =CE
               t
STREAM

                      DISTANCE FOR COMPLETE MIXING
  -UPSTREAM
     FLOWsQS
     CONCENTRATION=CS
                                                          7
DOWNSTREAM
  FLOW^QS-I-QE
  CONCENTRATION s
Figure 2-1 - Simple dilution model.
                              2-2

-------
                                                          •VII (2)
                                                           Revision  No.  0
                         CO  = (QE  •' CE)  + (QS  •  CS)       ,              (2-1)
                          '."'       QE  + QS .        .     '
  If the dilution  factor,  0,  is  defined  as:
                                                   1                    (2-2)
     .         .          ..             QE +  QS     1 +..D

 where   ,                                                           .

             D = QS/QE, the ratio of stream flow to effluent flow.

                       *                                  •   '               \
 CO may be calculated in terms of 


-------
                                                         VII (2)'
                                                         Revision No. 0
effluent  flows  and concentrations in a-concise and realistic fashion.
       The probabilistic dilution calculation procedure used in this
report permits  the probability distribution of downstream concentrations
(CO) to be computed directly from the probability distributions of the
flows and concentrations.
                                                  ' ,    ' •     .
       The first step in the use of the probabilistic dilution model is to"
develop the statistics of the concentration and flow of both the stream
and effluent.1  These statistics include'both the arithmetic and logarith-
mic forms of the mean ( \i ), standard deviation ( 
-------
                                                            VII (2)
                                                            Revision No. 0
          The probability distribution of the dilution factor, 

approaching 0 are of no practical significance to the calculations being performed since they occur at high dilutions. .For smaller streams relative to the size of the discharge, deviations from a log-normal approximation can.be appreciable*. They are large enough to ..Introduce significant error into the calculated recurrence interval of higher . stream concentrations. The error introduced is .almost always conservative; .that is, it projects high concentrations to recur more frequently than they . actually would. The appropriateness of this assumption is discussed in detail in Appendix B* . . . A procedure is provided in this report for accurately calculating the probability distribution of. the dilution factor (-) and stream concentra- tion (CO}. This numerical method uses quadratures and would be prohibitively tedious to perform manually. It has, therefore, been provided in'the form of a computer program which can be utilized on a microcomputer (Appendix D). For purposes of .presenting the approach in a form which can be solved manually, and thereby better illustrate the basic procedure employed, the-methodology description which follows In this section develops a log-normal approximation for the dilution function and then proceeds •with the calculations for stream concentration. Whether the log-normal approximation or'the quadrature calculation is used, the subsequent steps in determining the appropriate averaging period are the same. • . " , •' ' ' ' . • ' 2-5


-------
                                                          VII (2)
                                                          Revision No. 0
        The manual procedure (moments method) estimates the mean  and
 standard deviation of a log-normal approximation of dilution  by  first
 calculating, and then interpolating, between the 5% and 95% probability
 values.  The value of the dilution factor (o) for any probability  per-
 centile (a) is given by:
                                                                      (2-5)
                              (QE -f QS)  exp  (Za«rinD)
 where the value of Za is taken  from  a standard  normal  probability ta'ble
 for the corresponding value  of  a  (see Appendix  A).
         For example,  where    o  =  95%; 395  =  1.65
                              a  =   5%; 15   =  -1.65
                              a  =  50%; Z50  =  0
                              a  =  84J13%; Z84 =  1.0

        The  log  mean dilution  factor  is estimated by interpolating between
the 5% and  95%  values, calculated  above.

                •  mn  -s 1/2  [In (95) + In  (05)]                (2-6a)
            •                            ;       *                         *
        The  log  standard deviation  is  determined by the following formula
which,  in effect, determines the slope of the straight line on the log-prob-
ability plot:
                         —   __	2

From the log mean and log standard deviation of the dilution  factor
the arithmetic statistics are computed using

                                     2-6

-------
                                                           VII  (2)
                                                           Revision No. 0
        The arithmetic mean of the receiving water contaminant 'concentra-
 tion  (CO) downstream of the discharge after complete mixing, then, can be
 found by:        '                           .

                MCO= C"CE < %)! + *"cs (1 •-%  )?.''               (2-8)
       'The arithmetic standard deviation of stream concentration  is:
                                                                      (2-9)

        The. coefficient  of variationgf stream concentration (CO) is:

                                VCO  =  "CO/*1 CO                       (2-16)

        The arithmetic statistics  used -to derive the  log statistics will -be
 used .to develop  the  desired probability of exceedence.
           log  standard deviation =olnCO = Vln (1 +  vc§)         '•(2-lij'
          ,          .             .'.'• \  •      •            .
           log  mean               =  ^InCO =^n /  /C° , A        (2'12)
          :          ;                      \ V  I'+.-McS'-/

       The probability (or expected frequency) at which a value of CO Will'-
occur is determined by constructing a probability distribution plot on
log-probability paper.  This is accomplished by computing the 50th percen-
tlle and 84th percentile Concentrations and connecting  them with a straight
line:
                      '.'•'•    2-7. '

-------
                                                         VII (2)
                                                         Revision No. 0
             50% concentration » CO = exp ( M]nco)
             84% concentration = exp ( winco+  o"lnCO)
 Using this procedure, any concentration of interest can. be identified and
 its probability of occurrence scaled directly from the plot.

       Alternatively, the concentration that will not be exceeded at some
 specific frequency (or probability) can be calculated from:
C0
exp
                                         (Za
-------
                                                           VII  (2)
                                                           Revision No. 0
  probability  that  the  concentration,  will  be exceeded is obtained by sub-
  tracting the value  obtained  from  1.0.

  2.2    Choice of  the  Permit  Averaging Period
                                      *               *  »                       .
         In order to  examine the comparative effects of different choices of
  permit  averaging  periods on  water quality, it is necessary to define the
  relationships between the established effluent limit (EL) from the steady
  state WLA, the permit averaging period, the treatment plant performance that
  results, in particular the mean effluent (£E>, the downstream concentra-
 tion (CO), and a stream target concentration (CL).

        The objective of this section is to examine the  relationships  among
 these parameters in  order to be able to  predict the probability of an (adverse)
.waterDuality outcome based on known or estimated  stream  and  efflueht charac-
 teristics and the  choice of permit averaging  period.  .The approach is  based
 on the  assumption  that the  EL will be violated  with a particular frequency.
 The mean effluent.required  to meet this  level  of compliance, with EL is then
 calculated for each  of the  three permit  averaging  periods,  and  the probabi-
 listic  dilution  model  is  then used to develop a  probability distribution of-
 the downstream concentration  (CO)  for the  three  cases.  A level of acceptable
 adverse water quality  (a  decision  expressed.in terms of the probability or
 frequency of  experiencing a selected  high  value of  CO, such as the acute
 criteria concentration) °is then compared with the probability distribu-
 tions to determine the longest permit averaging period that meets the water
 quality  goals.                                                   '."  ~
                                     2-9

-------
                                                           VII  (2)
                                                           Revision No. 0
         The first  step  in this  sequence  is  to establish the relationship
  between the mean  effluent  (TIE), the effluent limit  (EL), and the permit aver-.
  aging  period.   In fact, what is required is the relationship between the
  treatment  plant performance necessary to meet the effluent limit as either a
  daily, weekly, or monthly maximum permit.  The reason for this is that the
  daily  variation of stream quality is governed,'not by the effluent limit
  which  is a regulatory  upper limit, but by the probability distribution
  of the daily effluent  concentrations which results from the design of the
  treatment  plant consistent with the effluent limit and the permit averaging
  period.  For log-normally distributed random variables, this.distribution is
  specified  by the  mean  effluent concentration, "CF, and its coefficient of
  variation, VCE.                                   %

        A particular effluent limit (say 30mg/l)  established  by permit as a
•maximum daily value would require a higher level  of plant performance (a
 lower mean effluent concentration) to avoid permit violations, than would the
          •   *     •                         f                    '     •'   i
 same limit specified as a maximum monthly average.  In the  latter case,
 excursions above  the effluent limit could be tolerated oh individual  days,
 without causing a violation of permit conditions.  The reason  for this  is
 that a monthly average of 30 individual  daily effluent concentrations  is less
 variable than the daily concentrations themselves.  Occasional  high daily
 concentrations  are averaged together with lower concentrations  to produce a
 less variable monthly average.   Hence,  treatment  pVant performance  is directly
 related to the  averaging period specified in  the  permit.

        In order to proceed  with the  analysis  a quantification of  this relation-
 ship is required.   Daily treatment plant effluent  concentration variations
                                      2-10

-------
                                                           VII  (2)
                                                           Revision  No. 0
  are well  described  by  a  log-normal  distribution parameterized by a long term
  average'concentration, lEVand a coefficient of variation, vCE.  Thus, -a   '
  relationship between these parameters and the permit effluent limit and
  averaging period is required.

        A  method to be employed is based upon an interpretation of what is
 meant, in practice, by specifying permit effluent limits as maximum values
 which may never be exceeded for the specified averaging period without
 causing a violation.  AsHaugh, et al. [2] observe, fixed upper limits
 which are never to be exceeded are conceptually inconsistent  with the
 stochastic nature of wastewater treatment processes and the effluent concen-
 trations  they produce.   Realistically,  some  exceedence  frequency  must  be
.acknowledged* regardless  of the averaging period  assigned.   For the present
 analysis,  it will  be assumed  that  the effluent  limit .specified by a permit is
 not  to" be  exceeded more frequently than  5 percent or 1  percent of the
 time.  Of  course, any other choice is possible.

        Once  a specific  choice  is made, say 1 percent, then the probability  '
 of compliance is  a= 99 percent and  that  establishes the fact that EL is the
,a-Percentile effluent concentration:  CEa.  This procedure, then,  gives a   .  '
specific probabilistic interpretation to the effluent limit.  It is the
effluent concentration that is exceeded with no greater frequency than  (1-q)
percent of the time.  If the permit is specified as  a daily maximum value,
then EL is the a-percentile of daily effluent concentrations.;  If  the permit
is specified as a weekly, (or monthly) maximum value, then EL is the a-percentile
of 7-day (or 30-day)  average effluent concentrations.
                                     2-11

-------
                                                         VII (2)
                                                         Revision No.  0
        In order to compute the long term average effluent concentration,
"CE, that would insure that CEa = EL as a daily, weekly, or monthly permit
the coefficients of variation are required for 1-day and 7-day or 30-day
averages of effluent concentrationse  Table C-2 presents representative
                                      t
values.

       Thus, the requirement that:                    .
                                  *             '           .
                                CEa=EL                             (2-15)

and for a coefficient of variation vCE» tne average effluent  concentration
EEcan be computed from

                             "CF =.Ra . EL                           (2-16)

•where the reduction factor relating CEa = EL to UE", that is,  Ra = UF/CEg is
                 Ra =   1 + V    exP £- Zaln  (1 + V)1           (2-17)
the ratio of the arithmetic average to the a-percentile  of  a  log-normal
random variable with coefficient of variation,  vCE-   Table  2-1 gives the
values of RQ for various coefficients of variation.

       The derivation'of this formula follows  from the expression for the
a-percentile of a log-normal  random variable:

                      CEa = exb (MIOCE.+ Za °lnCE)     ,             ,(2-18)
and the arithmetic average of a log-normal  random variable:
                                     2-12
                                                                                      1

-------
                     VII  (2)
                     Revision  No*  0
TABLE 2-1 - Reduction
Coefficient of
Variation
VCE
0.1
0.2
0.3
• ' , 0.4-
>0.5
0.6
0.7
0.8
' . » ,
0.9
i.o
1.1
1.2
1.3
1.4
i.5
factors for various coefficients of. variation.
•\
* ' >'
Reduction Factor •
. ' . R
-------
                                                           VII (2)
                                                           Revision No. 0
                                                                      (2-19)
  Thus:
                                                  lnCEJ
                                                                      (2-20)
and since exp (l/2anCE) .^i + V£ and alnCE =
(appendix A, page A-8)  equation (2-17)  follows.
                                                       (i
        At this point the effect of the choice of permit averaging period on
 treatment plant design can be illustrated.  If the permit averaging period is
 1-day, and the daily effluent coefficient of variation is VCE = 0.7 (for
 example, extended -aeration activated sludges, Table C-2), then for a 1 per- •
 cent violation frequency „= 99 percent, Ra = 0.281, which indicates that
 the long term average 'effluent concentration must be 28.1 percent  of the
 daily maximum permit limit.

        However,  if the  permit  averaging period. is  7 days,  then  the  coefficient
 of variation  of  7-day averages is vCE  = 0.6  and Ra  = 0.321.   Now the
 treatment plant  can be  designed  to produce a long term average effluent
 concentration of 32.1 percent  of the weekly  permit  limit.  For a 30-day
 average permit limit VCE .  0.45  and R0  = 0.404.  Hence, if EL = 10 mg/1 ,
the treatment plant average effluent concentration  must be 2.81, 3.21, or
4.04 mg/1 for a daily, weekly, or monthly permit specification, respectively.

       Hence  the selection of the .permit averaging period is  related to
the -CF required for each of the three averaging, periods' in order to
                                     2-14

-------
                                                           VII (2)
                                                           Revision  No.  0
  avoid exceeding, the  EL more  often  than  the  selected  frequency.   These average
  values are then used in the  probabilistic dilution model  (with other input
.  parameters such as "OT and IpT) to develop the probability distribu-
  tion of CO for  each  of the three permit averaging periods.    >           '

         The value of  CO in the probability distribution can be normalized in
  terms of a stream target concentration  (such as the chronic criteria concen-
  tration, CL) so that  the calculation can be.used for a wide variety of
  pollutants.  Stream  concentration  is therefore expressed in terms of
  $=  CO/CL, p being a  dimensionless unit of concentration.

        A convenient presentation of the resulting probability distribution
 jfiakes use of the  concept of return period.  For daily stream concentrations
 the  1 percent exceedence value has an average recurrence rate of one day
 every 100 days  so that  its average return period is 100 days.  Thus the
 return period for daily values is defined as:

             Return Period (days) = I/Probability of Exceedence       (2-28)
 The basic assumption in the  use  of return period as  defined above is that the
 event whose probability is being examined has a characteristic time  associated
  -.''.'.'             •            •         -                       •
 with it, in this case, one day for daily concentrations.   Thus,  it  is assumed
 that daily stream concentrations are  of concern,  and  each  event  corresponds
 to one day.                                                       .

        Figure 2-2 illustrates how the results of such an analysis can be
 expressed in a  plot  of concentration  versus  return  period.
                                      2-15

-------
                                                          VII (2)
                                                          Revision  No.  0
~  10
ea.
<
uj   5
   tr
   i
   o

      0.5
   S
   —  O.2
        O.I
                                            PERMIT AVERAGING PERIOD-*
                                                      30 DAY-—"
         	ACUTE	•^-**sr	^
                                            	CL= CHRONIC —
                                I   9   filltfl     1    I  lllttl
                   0.5    I     2      5     10    20
                        REf URN  PERIOD (YEARS)
50   !OO
    Figure 2-2  - Illustration of analysis results:  stream concentration
                 versus return period for three permit averaging periods.

       The stream target concentration (CL) for a typical WLA is the chronic
criteria concentration of the pollutant under consideration.  The use of the
chronic criteria as the stream target concentration is convenient for the
comparison of permit averaging periods because it represents a specific and
frequently used procedure.  The analysis that follows does not attempt to
quantify the frequency with which chronic criteria concentrations are met by
either the conventional WLA procedure or the guidance provided for selecting
permit averaging period.  Instead, the analysis is designed to relate the
choice of the permit averaging period to the frequency with which severe,
short term water quality impacts are expected as a result of an effluent
limit.                             :
                                     2-16

-------
                                                          VII (2)
                                                          Revision No.  0
        These short-term impacts are perhaps most effectively  evaluated
 with respect to acute criteria concentrations.   If the  stream concentra-
 tion exceeds the acute.criteria as  a result of  an occasional  high  daily
 effluent loading, the result  is presumed  to be  an undesirable impact.
 Hence,  there is a direct  connection between the permit  averaging period
 and the probability  of acute  criteria violations.   Specifying that the      '
 WLA requirement be met. as a daily maximum permit limit  significantly reduces
                                       i             • ''   ''.•'''    '  •     ,  '
 the possibility of acute  criteria violation since the effluent limit is
 specified using the  chronic criteria,  which is  always a smaller concentration.

        The  frequency with which daily  stream concentrations are allowed to
 exceed  acute criteria is  a regulatory  decision.    The analyses presented
 herein  employ a frequency that  corresponds  to a  1-day in 10-year recurrence,
 on  average.   The choice of 10 years  is, of  course,  used for example
 purpose only but it  i.s consistent with the  10 year  return period that is
 conventionally  used  for the design  stream flow.         "-       '

        the  results of  the permit averaging  period analysis are presented in
terms of  CO/CL which  is exceeded with a.particular  frequency,  such as once in
 10 years.  This  ratio can  then  be compared  to the acute-to-chronic criteria
concentration ratio  for the pollutant of  concern.   For pollutants with large
acute-to-chronic  ratios,  occasional  large daily  fluctuations can be tolerated;
and a 30-day  permit  averaging period provides protection from acute criteria.
 violations.   Conversely,  pollutants with  small acute-to-chronic ratios are
more likely to require shorter day permit averaging periods.   Site specific
 This is currently under EPA study.
                                     9-1 7

-------
                                                         VII (2)
                                                         Revision No. 0
considerations, primarily the ratio of effluent.to stream flow and stream
flow variability, become significant in these cases.

       The final translation of the selected averaging period option to
permit limits requires consideration of the monitoring frequency.  The
method assumes either daily monitoring or other monitoring adequate to
describe the performance of the plant on a monthly basis.  If such conditions
are not met, alternate limits may be calculated which incorporate monitoring
frequency, or monitoring frequency may be adjusted so that these conditions
are met.                                                              -•.,.•'
                                    2-18

-------
                                                             VII (3)
                                                             Revision No.  0
                        .  .          CHAPTER 3                      r   /   .

                              fEXAMPLE COMPUTATION           •

            • (      .'    ' • •      .                   • • •.    '               "  "
         This  chapter presents  an  example-problem,  showing step by step

 computations  using  the methodology.described in the previous chapter.  A

 set of  hypothetical conditions that apply to a site-specific situation is

 assumed, and  an analysis is performed to determine the. effect on receiving

 water quality resulting from the  assignment of different  permit averaging

 periods to the steady-state model  output.  The steps used.to conduct this

 analysis are  summarized below in  Figure 3-1.  The format  used in this

 chapter presents  data and computations  on the left-hand page, and pertinent

 commentary and  supporting discussion  on the facing page immediately  opposite

 those computations.  The manual computation using the moments approximation

 is described  first, followed by an  analysis using the computer program

 (POM-PS) in Appendix D.   Both examples  use  the  same set of hypothetical

 site-specific conditions.                                   ,
         GIVEN:
  •acute and chronic
   toxicfty
  •design flows
  •flow and concentration
   variability for specific
   averaging periods
       STEP "1:
Compute statistical
parameters of stream
and effluent flow and
concentration
         STEP 5:
       Repeat steps
  1 through 4 for remaining
      averaging periods
Compute
frequency distributions
       STEP 2:
Compute statistical
parameters of dilution
factor
      STEPS:
Compute statistical
parameters of the
resulting stream
concentration
, STEP 6:
Select optimal permit averaging period
»-

STEP 7:
Translate into daily, weekly, and
monthly permit limits
Figure 3-1  - Step procedure to  select optimal permit  averaging period*

-------
VII (3)
Revision
                                                                  No. 0
                   .HYPOTHETICAL SITE-SPECIFIC CONDITIONS
This section provides an example of the type and amount of information
required to perform the analysis.  It also establishes the basis for the
example computations and assumes that pertinent site-specific conditions
are as follows:                          ,

A.  Site-Specific Waste Load Allocation (WLA) Results ;   -.

   .The pollutant (P) to be allocated has a chronic toxicity  concentration
    (CL) of 2.5, and an acute toxicity concentration of 6.25.
    WLA policy for the agency performing the analysis  is  to  use  7Q10 as
    stream design flow, to use the design capacity of  the treatment  plant
    as the effluent flow, and to compute (e.g., using  a water quality
    model) the effluent concentration of pollutant (P) that  will  result  in
    a stream concentration after dilution less  than or equal  to  the
    chronic value (2.5 = the stream, target concentration, CL) .  For  this"
    example, it was assumed that:
               Design Effluent Flow (QE)  = 5 MGD = 7.77 cfs

               Design Stream Flow (7Q10)  = 23.3 cfs
    The stream target concentration  (CL = 2.5) will  be met  under these
    design flow conditions,  when  the effluent concentration  is  CE =  10.
    Therefore,, based on the  WLA analysis, the effluent limit  (EL) for
    pollutant (P)  is specified  by the permit as:
               EL = 10
                                                                                 ;j
                                     3-2

-------
                                                         :vn (3)
                                                         Revision No. 0.
                                COMMENTARY
	from EPA Criteria

       State water quality standards do not usually specify both values;
       they are usually based on chronic values.

       (Any concentration units may be assigned; stream concentrations  will
       have to be in the same units.)
       7Q10 (the lowest 7-day average stream flow with a recurrence
       interval-of 10 years^  is the most common  "design stream flow".
       Some states .use other  values (e.g.,  30Q5).  This analysis  uses  the
       numerical  value of the "design flow".  However,  although the.
       example terminology uses "7Q10",  it  should be  interpreted  as
       "design stream flow" and the appropriate  value substituted,  regard-
       less of the averaging  period or the  recurrence interval  on which  it
       is  based.   (For example, if design flow in a state were 30Q5,
       assume that 30Q5 =23.3 cfs)v

           NOTE;   The only exception to  this is  in Figure C-l,  in which
           -       the ratio of 7Q10 to average stream flow  is  used  to
                  estimate the variability  of  daily flows in the  absence
                  of  a specific local  analysis.   The  use of this  figure  is
                  not requisite to  either-the  analysis  methodoloay  or  the
                  computations.     • •                              -..''-
            (QE> CE),+  (QS  . CS)

                  QE +  QS


                        + (23.3
      25 -       "
                   7.77 + 23.3
      CE = 10 = EL


       ...     '.  ' '     '    '    •••."' 3-3

-------
                                                          VII  (3)
                                                          Revision  No.
                   HYPOTHETICAL SITE-SPECIFIC  CONDITIONS
                                (continued)

B.  Site-Specific Conditions         „

    Stream Flow                             Mean Flow  (T|S~)  =  467 cfs
                            Coefficient of Variation (VQS)  =  1.5
    Upstream Concentration                       Mean  (UF) .= 0
                            Coefficient of Variation (VGS) - 0
    Effluent Flow          .     '                 Mean ("DJ) = 7.77 cfs
                            Coefficient of Variation (V) = 0.20
                                     3-4

-------
                                                   VII (3)
                                                   Revision No. 0
                          COMMENTARY
 Stream flow data are obtained from analysis of flow gaging records
 for the stream in question; where the stream reach is ungaged,  it
 is obtained by extrapolation from an appropriate record.

 At present, records are not normally analyzed for the coefficient
 of variation, although the computation is straight forward and  can
 be readily incorporated into a routine statistical  analysis of
 daily stream flows.  In the absence of specific analysis results,
 the coefficient of variation of daily stream flows can be estimated
 using the material  presented in Figure C-l.
 Upstream concentration can be assumed to be zero if the stream
 concentration of the pollutant is very low compared to the dis-
 charge, or if the effect of the discharge only is to be examined.
 Site-specific values for upstream concentration statistics would be
 obtained from analysis of an appropriate STORE! station, or from
 local monitoring records.  If upstream concentrations are assigned,
 enter data here and in the equations when called for.
The design effluent flow is assumed to be the mean effluent flow.
The variability of daily effluent flows for a new facility must be
estimated on the basis of available data for existing treatment
facilities (such as Table C-l).  For an existing facility being
expanded, or simply re-permitted, variability could be based on an
analysis of past plant records.  For many industrial  dischargers,
this data will be available in Book. VI (Design Conditions) of the'
waste load allocation technical guidance document series (specific-
ally, in Chapter 4:  Effluent Design Conditions).         .'"•'•
                              3-5

-------
                                                          VII (3)
                                                          Revision  No.  0
                    HYPOTHETICAL SITErSPECIF.IC CONDITIONS  .
                    '        '.     (continued)       _ .       V •


     Effluent Concentration                      Mean  ("CT)   = (*)

                             Coefficient  of Variation  (VCE)  = .7


        The mean concentration is a function 'of the  permit averaging
        period and is  that concentration  required to avoid exceeding the
        effluent limit concentration (EL)  more often than the compliance
        probability.

        The coefficient of variation, for  the hypothetical treatment plant.
        is  not known because the plant  has yet to be constructed.  Assuming
        that the plant will produce an  effluent with a variability similar
        to  the values  given in Table. C-2,  the  following values are used:  ,


                 Permit Averaging      Coeff.  of Var.
                    Period
                     Daily                 0.70
                     7-Day                 0.40
                     30-Day                0.20
       Equation 2-17 is then used to determine the mean effluent concentra-
       tion of (P) Which is required to avoid a violation of EL more often
       than the compliance probability.  For this example, assume that the
       exceedence probability is 1 percent.  For  0 = 0.99 percent,
       Za- 2.327.  For VCE = 0.70, Ra * CF/EL is:
                             exp [
                 = V 1 + 0.49 exp [-2.327   Vln (1 + 0.49)]
                 = 1.221 exp [-2. 327  . 0.6315]
                 = 0.281

       The reduction factor for 7-day and 30-day averages are computed
       similarly with VCE (7-day) = 0.40 and VCE (30-day) = 0.20.  The
       results are:

                    Coeff. of Var.
                     of Averaged
                       Effluent         Reduction          Required Mean
     Permit         Concentrations     - Factor            Effluent Cone
Averaging Period        (vrF)           Rg = CT/EL         ("CT = Rg EL)
                          0.70           0.281           0.281 • 10 = 2.81
     7-Day     .           0.40           0.439           0.439 - 10 = 4.39
     30-Day               0.20           0.643           0.643 • 10 = 6.43
                                     3-6

-------
                                                       VII  (3)
                                                       Revision  No.  0
                             COMMENTARY
    The mean effluent concentration that a treatment facility, is
    capable of producing is influenced significantly by process selec-
    tion.  FOP this example, it will be .assumed that process selection
    will be made following the issuance ,of a permit, and influenced by
    its provisions.                  .'"'•

    The mean effluent concentration that a facility is required to
    produce is influenced by the permit averaging period and the vari-
    ability of effluent concentrations of the pollutant in question.

    The analysis employed.here, which bases permit averaging period
    selection on receiving water impacts, is based on, exceedance of the
    acute criteria on a daily basis.  Therefore, all subsequent stream
    impact computations (Step 4) are based on the coefficient of variation
    of daily effluent concentrations, or 0.7, as shown.

    The mean concentration is shown by (*), because a different value
    is used for each permit averaging period.                   .
 .-  .The recommended exceedence probability for the effluent limit is
.. r  either 5 percent or 1 percent.  For 5 percent, Za would be Z'g$ =
    1.645.
    Longer averaging periods reduce the variability of effluent concen-
    trations, and allow permit exceedance limits to be met with higher
    effluent means.  Computation of the required mean (tlT) uses the
    values of VCE for the corresponding permit averging period.
                                  3-7

-------
                                                          VII  (3)
                                                          Revision  No,  0
                   EXAMPLE  COMPUTATION - -HAND CALCULATION
 This  section  illustrates the  hand  computation  using  the moments approxima-
 tion  to evaluate  the  stream concentration  probability distribution.


 STEP  1;   Compute  statistical  parameters  (arithmetic  and logarithmic) of
          inputs using relationships  for  log-normal distributions  (see   .
          notes on page 3-9 or Appendix A for equations).

          o  For the mean effluent  concentration  ("CF) for a  30-day permit
            averaging period  with  X  = CE,  that is for the variable CE:
ARITHMETIC
i can
Coef. Vaj
Std. Dev,
Median
LOGARTHMIC
Log Mean
Log Std.
0
x
Stream
How: QS
Effluent
Flow:- QE
v MX/ — 	 	 	 - v
r( 11 * \ — _«»_«_>_^
• v vx) --------(
• ( CTx) = MX * vx = (6.43
(x) = nv/ Jl + v2 =
• y w v
** ~ .A
( Minx) = 1" (^<) = I" (5.
Dev. ( 
-------
                                                          vii (3)   ;
                                                          Revision No. 0
                                 COMMENTARY
 The following parameters  are used subsequently:
                                                                        i
                                Arithmetic  _                Logarithmic
Input
Parameter

Stream Flow
Stream Cone .
Effluent Flow
Effluent Cone.


QS
CS
QE
CE
Median
x
QS
CS
QE
.CE
Mean
"x
"QS .
MCS
"QE
"CE
Std.
Dev.-
<*x
ffQS
'. CTCS
°QE
CTCE „
Coef.
Var.
vx
VQS
VCS
VQE
VCE
Log Log
Mean S.D.
Mln x CTln x
MlnQS CTlnQS
MlnCS nnCS
plnQE alnQE
plnCE alnCE
The following definitions and equations summarize the relationships among
the statistical parameters of log-normal random variables.           "


       Arithmetic              Terms           '    Logarithmic


           x         Random Variable                  In x


            nx       Mean                              Mln x


            o        Variance                        .  -of
             x.                                          In x

            ox       Standard Deviation     _  ..       o-jn x   .


         • '  vx       Coefficient of Variation         (not used)
                    "    •                  •                           %
          rt_i
          x         Median                            (not used)



                                            = In
         = exp  [Mlnx]
                                      'inx•
     °x =
                                     3-9

-------
                                                         VII  (3)
                                                         Revision  No. J)

                  EXAMPLE COMPUTATION..- HAND CALCULATION
                          :'(continued)


STEP 2;  (a)  Compute the log standard deviation  of the flow  ratio  QS/QE = D,
                 alnD = V CT?nQS + c fnQE + 2P-°inQS .• CTlnQE
               The first two  terms  are  taken  from the table in Step 1
               (and squared).  Since, for this  example,  flows are not
               correlated (P= 0), the third term drops out.  Therefore,
                °lnD --v/Cl-0857)2 + (0.1980)2  =  1.1036
         (b)   Compute the  5th  and  95th  pereentiles of the actual distribu-
              tion  of the  dilution 'factor  (0).
                            QE
                   0,
                    'a
                        (QE  + QS)'  exp(Za ="^5	$	1-0546
                                     3-10

-------
                                                  VII (3)
                                                  Revision  No.  0
                         COMMENTARY
This equation accounts for any correlation  that may exist between
stream flow and effluent flow; e.g.,  where  higher effluent flows
tend to occur during periods of high  stream flow.               .
Ordinarily, ttoere is no reason  to  expect  any  such correlation;
therefore  P= 0, and the computation  in step  (a) is simplified as
shown.
     *95 =            QE
          (QE + QS)exp (ZaolnD)

         =       .      7.62
          (7.62 + 259) exp  [(1.645), (1.1036)]

              7.62           ,
           7,62 + 1591

         = 0.004766        .

-------
                                                          VII  (3)
                                                          Revision No. JO

                 ,   EXAMPLE COMPUTATION  -'HAND CALCULATION
                             !(continued)'~~



           (d)  Compute arithmetic.statistical parameters (using equations on
              . Page 3-9 and tabulate for convenience.


                     	  Arithmetic                  Logarithmic

                     Mean   Median   Std Dev   Coef.Var "    Mean    Std- Dev

      Dilution (0)  0.0471  0.0270   0.0673      1.43'    -3.6115   1.0546
      Factor
 STEP 3:  Compute the statistical  parameters of the resulting in-stream
          concentration (CO).


          (a)  Compute the arithmetic mean concentration using previously
               tabulated values, using Equation 2-8.


                HCO = [MCE • M0  ] + [MCS • (1 -  *"«, )3

                    = [6.43 • 0.0471] + [0] = 0.303

                        '             '              ,


          (b)  Compute the standard deviation,  using Equation  2-9.
                      (  MCE-  MCS)2        /0.0673)2  ..  (6.43  - 0)2   /0.187



                                             (4.50)2  .                * 0,137
                    =  0.569
         (c)  Compute and tabulate for use in subsequent graphical or
              other summaries, the other statistical parameters-of stream
              concentration*


                    __	Arithmetic                  Logarithmic

                    Mean •  Medi an   Std Dev   Coef Var     Mean    Std Dev

Stream              0.303  0.142    0.569       1.88      -1.95     1.23
Concentration (CO)                            .    •     -     '
                                     3-12

-------
                                                         VII (3),
                                                         Revision No. 0
                                COMMENTARY
— -  The equations are as follows:
           = exp  [-3.6115 +   (1,0546)2]
           =0.0471
        V
     : Jexp (o-2
               l
                    [(1.0546)2]-i
          = 1.429                     '                                    •  .
       
-------
                                                         VII (3)
                                                         Revision No. 0
                  EXAMPLE COMPUTATION - HAND CALCULATION
                                (continued)      '.       ~ .     .

                                         h    /'      i    '      -     f
STEP 4:  Use the statistical  parameters of stream concentration computed
       "  in the previous step to construct graphical .or tabular displays
       .  summarizing the, frequency distribution.        .

         (a)  To construct a probability plot using log-probability graph
              paper:       '   ,                                     .

              o  The median concentration is plotted/at the 5.0th percentile
                 position.                                 ;

                                CO = C050%
              '                  '   = exp (-1.95).

                                   = 0.142                /


              o  Any other plotting position is determined as follows:

                 (1)  From Table A-l, select a probability (  a )  and
                      determine the corresponding  value of Za.  For
                      example,                                   •

                      Probability = 0.841 (84%) --- ' - -   Z84.1% = 1-00

                      Probability =. 0.159 (16%) --- - -   Z15.9% =.-1.00
                 (2)   Compute the concentration at probability  ( a )  from
                      the log mean and  log  standard deviation of  stream
                      concentration (CO)..

                         C0fl  = exp(  M
                      84% plotting  position  .

                         C084% = exp(-1.95 + 1.00 '•  1.23)  = 0.487  '


                      16% plotting  position

                        "C016% = exp(-1.95 - 1.00 •  1.23)  = 0.0416
                 (3)   Plot  these  concentrations  on  log-normal  probability
                      paper and connect  with  a straight  line,
                                     3-14

-------
                                                           VII (3)
                                                           Revision No.  0
                                  COMMENTARY
5  10   3O  50  7O    90'
% EQUAL OR LESS THAN
                                                         99
99.99
 Figure 3-2 - Sample stream concentration  versus  probability plot  for
 30-day averaging period.
 The probability plot indicates,  for example, -that the  stream concentration
 ?  P°3^upanMp) will  exceed  a concentration of  1.0, at a frequency
 (probability)  of about 5%.. Since  the  analysis is based on daily values,
.-cms is interpreted  as:   5% of all  days wit! have stream concentrations
 greater than 1.  .                                      ,
                                      3-15

-------
                                                         VII (3)
                                                         Revision No. 0
                 • EXAMPLE COMPUTATION.- HAND CALCULATION
                             ~'  (continued)
STEP 4 (continued)
         (b)  To construct a recurrence interval  (return period)  plot
              using log-log graph paper:,

              o  The formula used in the previous step
                   can  be rearranged:

                   ,    In (CO.)  - p
              The log mean  and, log  standard  deviation were determined in
              Step 3:
                 
-------
                                                          VII (3)
                                                          Revision No. 0
                                 COMMENTARY
        Probability results can be misleading for the water quality issues
        being considered here,  unless interpreted very carefully.   FOP
        example,  a-1% probability of exceeding a significant stream concen-
        tration jneans that this occurs nearly 4-times in  1 year, and for
        more than a month of individual  days  over a 10 year period.
        Expressi-R-g results as recurrence intervals is believed  to  provide a
        more useful  expression  of analysis  results.  ,   •   '> .
       100
      o
      o

      o
UJ
o
g;
o
      UJ I
      o: '
        .1
                                                    10 YEAR
                                             RECURRENCE INTERVAL
           .	ACUTE TOXIC1TY CONC. = 6.25 —	—^a?-. _ J	'.		
           	CL= 2.5
                                                      I
                                   - ' '
                                                    JL
         .OI  .02   .05   .1   .2,.   .51    2    5    10   20
                      MEAN RECURRENCE  INTERVAL' - YEARS
                                                         50
Figure 3-3 - Sample stream concentration versus mean recurrence interval
             for 30-day averaging period.


       Note that the acute concentration assumed for the pollutant (6.25).
       is exceeded an average, of once every 2.6 years.  If the exceedance
       criteria to be met is an average of 1 acute toxicity exceedance
       every 10 years, then the assignment of a 30-day permit averaging
       period is insufficient; shorter averaging periods must be examined.

       However, if the pollutant had an acute concentration of 12.5 (or an
       acute-to-chronic ratio of 5), the recurrence interval  of 20 years
       would-be sufficiently protective for acute events.
                                     3-17

-------
                                                         VII (3)
                                                         Revision No. 0
                  EXAMPLE COMPUTATION - HAND CALCULATION
                      ~          (continued)               '
STEP 5:  Compute the receiving water quality impact that would result
         from assigning other permit averaging periods.

         Repeat Steps 1-4 using the values for UE" that have been
         calculated for weekly and daily permit assignment.
              7-day permit average ------- cT = 4.39

              Daily maximum permit average - - - CT = 2.81

         All  other imputs remain unchanged.
         When the computations are repeated using these values,  the
         statistical  parameters for stream concentration (Step 3)  that  are
         developed are as follows:
                            STREAM CONCENTRATION (CO)  STATISTICS
Permit
Averaging
Period
30-Day
7 -Day
• 1-Day
Mean
MCO
0.303
0.207
0.132
Median
••w
CO
0.142
• 0.0971
0.0622
Std. Dev.
""CO
0.570
0.389
0.248
*
Coef . Var.
VCO
1.88
1.88
1.88
Mean
-1.95
-2.33
-2.78
Std. Dev.
""InCO
1.23
1.23
1.23
         Probability and  recurrence  interval  plots  are  then constructed as
         described  in Step  4 to  provide  a  graphical  comparison of the
         influence  of alternative  choices  for averaging period on the
         frequency  of exceeding  acutely  toxic concentrations of pollutant
         (P)  in  the receiving  system.  ,
                                     3-18

-------
                                                           VI I.(3)
                                                           Revision  No.  0
                                 COMMENTARY"
                                            PERMIT AVERAGING PERIOD
            .Ol   .1
               I  2  5  IO     30 5O 70    90
                    % EQUAL OR LESS THAN
                                                            99   99.9 99.99
Figure 3-4 - Concentration versus probability plot for 1-, 7-, and
            ,  30-day averaging periods.             '
 100



8

o
           10
         I
         111
         o
         o
         o
         <  I
         iu
         o:
                                                          I
                                                        10 YEAR
                                                 RECURRENCE INTERVAL
    ;	ACUTE TOXICITYCONC



      —CL = 2.5	-,— -
         30 DAY
            7 DAY
                     PERMIT AVERAGING PERIOD
            .1
            .01  .02
                                                          I
                      * ' •
                               LUU
JL
             .05   .1   .2    .5   I     2     5   10  20
                 MEAN RECURRENCE  INTERVAL - YEARS
      '50
Figure 3-5 - Concentration versus mean recurrence  interval  plot  for  1-,
             7-, and 30-day averaging periods.

                                      3-19

-------
                                                         VII (3)
                                                         Revision No. 0
                  EXAMPLE COMPUTATION - HAND CALCULATION
                                (continued)'"
STEP 6:  Select the appropriate permit averaging period,
         The appropriate permit averaging period is chosen to provide.
         an acceptable level  of receiving water quality.   The decision  is
         based on the assumption that an unacceptable exceedence of the  acute'
         criteria in the receiving stream is more than once every 10 years,
         on average.

         Therefore, the permit averaging period selected is the highest
         one that does not result in  a mean recurrence interval  for acute
         criteria violations  that is  less than 10 years.  For this  example,
         recurrence intervals for a stream concentration of 6.25 are
         approximately

              30-Day Avg. Period = 2.6 years

               7-Day Avg. Period = 7.7 years

               1-Day Avg. Period ='31 years

         For the  site specific conditions assumed  for this  example,  a
         1-day permit averaging period could" be assigned to the  effluent
         limit of 10.  However, as shown below using  more exact  calculations,
         a  7-day  permit averaging period is sufficiently protective  for  acute
         events.   Thus  a 7-day permit averaging period is assigned to the
         effluent limit of 10.
                                    3-20

-------
                                                   VII (3)
                                                   Revision No. 0
                          COMMENTARY
 For marginal  case's,  it should be recognized  that  the  projections made
 using the moments  approximation  tend  to  be conservative.   As  shown
 below the more  exact recurrence  intervals are  6.4,  32,  and 280 years.
 The  acceptable  frequency of. acute  criteria violation  is, of course,
 a  policy  decision.   Alternate  levels are evaluated directly from
 Figures  3-3.and 3-4.             ,                               '   .
The moments approximation used for the foregoing computations
(because it approximates" the distribution of dilution factor
with a log-normal distribution) provides an approximation of the
probability distribution and recurrence interval of the stream   ,
concentrations.                                  •        ,

An exact computation that avoids the necessity of this approxima-
tion, is provided by use of the computer program detailed in the
next section and in Appendix D. .In this case,its use is warranted
since a 7-day permit averaging period is sufficiently protective.

Based on the selection of the 7-day permit averaging period, the
maximum 7-day average permit limits -EL = 10 mg/1.  Thi-s permit
limit-is equivalent to a long-term average effluent concentration
CE = Ra . EL = (0.439)(10) = 4.39, with coefficient of variation
daily effluent concentration (VCE) = 0.7.  Thus, the design of the
treatment facility, and the selection of treatment process should be
made to meet these specifications  of UE = 4.39 mg/1  with coefficient
of variation of daily effluent concentrations VCE- * 0*7.
                              3-21

-------
                                                         VII (3)
                                                         Revision No, 0
                  EXAMPLE COMPUTATION .--.HAND CALCULATION
                            ~~(continued)
STEP 7;  Compute permit limits for other averaging periods (daily maximum
         and monthly) and exceedence percentiles (1 percent, and 5 percent)
         that are consistent with the treatment performance level established
         in Step 6.


         At this point in the analysis, it has been determined that  assigning
         the effluent limit of EL = 10 as a weekly permit, applicable  to  7
         day averages of the daily concentrations, is sufficiently protective.
         .This choice is based upon an effluent limit violation frequency  of
         one percent.  The mean effluent concentration for these choices  is
         "CT=4.39.     '                               -

         If it is assumed that the same violation  frequencies  apply  to the
         other permit concentrations, then they can be computed directly:

                 Permit Limit ="CT/Ra

         since Ra  ~ "CT/CEa and the permit limits.are assumed  to be  the
         o-percentile concentrations for each  averaging period.

         If other violation frequencies are desired, for example,  5
         percent, then permit limits of this frequency can also be-calculated
         using the appropriate Ra for* 1-a = 5  percent.  The table below
         presents the results for the example  considered  above.
          Permit
        Averaging
          Period

           1-day .
           7-day
          30-day
Coeff. of Var. of
Avg.'ed Effluent
 Concentration3
     VCE

     0.70
     0.40
     0*20,
Reduction Factors0
	w

0.281
0.439
0.643
 5% '

0.432
0.571
0.736
                      Permit Limitsc
                       1%

                      15.6
                      10.0
                       6.83
 5%

10.2
 7.69
 5.96
         It should be pointed out that any or all of these permits are  '
         equivalent in the sense that a treatment plant meeting any of these
         requirements will also meet the desired water quality goal.  Of
         course, this is true only if the actual coefficients of variation
         for daily values and 7 and 30 day average plant effluent concentra-
         tions are as specified-.
        aThese are assumed to be representative of the treatment plant'
         effluent behavior.
        bTable 2-1, equation 2-17.
        cpermit limit = "OT/R • "CT = 4.39.
                                    3-22

-------
                                 COMMENTARY
                                                          VII (3)
                                                          Revision No. 0
 Permit Limits

 Reduction factors  (see p.  3-6)

 Choice of averaging period
 (from step 6).

 Value for the selected averaging
 period (from  step  6 -  steady
 state model output).

 Permit limits using'reduction
 factors,  RQ'S

 Long-term average effluent
 concentration, cF  (see p.  3-6)

 Coefficient of variation of
 daily, weekly, and  monthly
permit limits  (see  p.  3-6)
    Daily
    Maximum

    0,281

    no
10
0.281

 4^39


 0.7
        • 15.6
Weekly

0.439

yes


10




10.0


4.39


0.4 ,
                            Monthly

                            0.643

                            no
                                     =6.8
                               4.39


                               0.2
       The long term average effluent concentration for the required level  of
       treatment is equal to 4.39 mg/1 with the coefficient of variation of
       daily effluent concentrations equal to 0.7.  To meet the water quality
       standard at .the state specified design flow and to meet the acute
       criteria at all times except for 1 day once in 10 years, the treatment
        *1 l«ies need to be built to meet the long term average concentration
       of 4,39 mg/1 with coefficient of variation of daily effluent concen-
       tration VCE =0.7. -The permit limits derived above are based on
       daily, weekly, and monthly reporting procedures.  If less than
       adequate monitoring is required, the appropriate permit limits
       must be derived using the long term average and equivalent coefficient
       of variation.                                                       .
                                     3-23

-------
                                                   VII (3)
                                                   Revision No. 0
            EXAMPLE COMPUTATION - HAND CALCULATION
                   ''       (continued)       "   ~~
 Recapitulation-

 In order to aid in the understanding of the suggested procedure,  the
 sequence is reviewed below in outline form.
 1.   Establish streamflow characteristics.
 2.   Establish effluent flow characteristics.
 3.    Establish effluent  concentration  variability  characteristics
      (VCE)  for daily  values-and  7  and  30 day  averages.

                                   Coefficient of Variation
             Averaging Period       •        VQ^

                   1-day                     0.7
                   7-day                     0.4       .
                  30-day                     0.2


.4.    Establish effluent  limit  from steady state waste!oad al location.

                EL = 10


 5.,   Establish violation  frequency of  EL.       -


                1-0=1%
                  «=  99%   ' '     '                    '

      and assume CEQ -  EL
                              3-24

-------
                                                   VII  (3)
                                                   Revision  No.  0
                         COMMENTARY
 1.   These should .be site specific since the computation is usually
     sensitive to the values*
2.   Mean effluent flow i's important, but the coefficient of variation,
     since it is usually small, is usually .not significant if VQE « VQS-
3.   These coefficients of variations specify the behavior of the
     'daily values and temporal averages of effluent concentrations.
     More detailed evaluations for industry specific or pollutant
     specific .situations are required to be more definitive.  The
     values used are not suggested as universal.       .    .
4.   The analysis presented in this manual  does not evaluate  the
     degree of protection afforded by this  choice.   That  is,  the
     probability of violation of the chronic criteria is  not  calculated,
     It is assumed to be sufficiently protective.

5.   The choice of violation frequency is necessary in order  to  give  a
     specific probabilistic meaning to EL.   Reasonable values appear
     to be one or five percent.  However, a problem may arise if too
     frequent a violation frequency is chosen.   It  may turn out  that
     even specifying the permit as a daily  maximum  does not insure
     that acute criteria violations 'are sufficiently rare.  In this
     case, a lower probability of violation must be specified.
                              3-25

-------
                                                  VII (3)
                                                  Revision No. 0
           EXAMPLE COMPUTATION - HAND CALCULATION
                         (continued):,
6.   For a(step 5) and coefficients of variation (step 3) compute
     ratio of mean effluent to effluent limit, Ra = EF/CEa and the
     resulting mean effluent concentration rF for each averaging
     period.
            Averaging Period

                  1-day
                  7-day
                 30-day   •
     Reduction Factor
Mean Effluent Concentration
  R~
 0.281
 0.439
 0.643
2.81
4.39
6.43
7.   Evaluate each mean effluent concentration using POM'to compute
     the return period of acute criteria violation.  Choose the
     appropriate averaging period.
Averaging
Period
1-day
7-day
30 -day
CT
2.81
4.39
"5733.
Moments
Approximation
31
7.7
2.6
                            Return,Period (years) for
                                     CO = 6.25
                                             Quadrature
                                               Method

                                               281
                                                31.8
                    > 10 years
8.   Establish appropriate permit limits for other averaging periods.
     "CE » 4.39, 1 -a =1%.
     Averaging Period
        Permit Limita
1-day
7-day
30-day
0.70
0.40
0.20
0.281
0.439
0.643
15.6
10.0
6.83
     aPermit Limit = "CT/Ra ; 1% violation frequency.
                              3-26

-------
                                                   VII (3)   ,
                                                   Revision  No.  0
                          COMMENTARY
      This .calculation  makes  the  connection between the effluent limit.
      and  the  mean  effluent concentration  required to meet the effluent
      limit if it is  assigned to  daily values or 7 or 30 day averages.
      A treatment plant designed  to produce UF and whose variability
      is as specified in  (3)  will meet the effluent limit with one
      percent  violation frequency.
7.   The three treatment plant .designs (the three mean effluent
    •concentrations) and the daily effluent variability are used in
,  .   POM to compute the return period of an acute criteria violation.
  -  The moments approximation is sufficient if the return periods are
     significantly less than or greater than the 10 year criteria
     violation frequency being examined.  In this case, the 7-day
     averaging period result is close to 10 years and the more accurate
     computer method is used to improve the accuracy of the calculation,
     The calculation indicates that a mean effluent concentration of
     Ct = 4.39 and a daily vr,g = 0.7 is sufficiently protective
     for acute criteria violations.  This, then, is the base's  for the
     treatment plant design.
8.   The permit limits for the other averaging  periods  are now calcu-
     lated to be consistent with  the treatment  plant design.  That is,
     these permit limits  are consistent  with  effluent mean and coeffic-
     ients of variation as indicated,  and  specify the same performance.
     Thus, they are equivalent requirements.
                              3-27

-------
                                                          VII (3)
                                                          Revision No. 0
                   EXAMPLE COMPUTATION - COMPUTER PROGRAM


 This section illustrates the use of the PDM-PS computer program (included
 and described in Appendix D) to the solution of the example presented  in
 the previous section.  The site-specific conditions used to define input
 values in the .previous section are used in this section as  well.

 The PEW-PS is structured to accept inputs in the form  of statistical
 parameters and ratios, determined readily from the  data. The  following
 ratios are entered for this example computation:


        Stream Flow Ratio,         7Q10/33" = 23.3/467 =  0.05

        Effluent Dilution Ratio   7Q10/3T = 23.3/7.77 = 3.0.
        Effluent. Concentration       /EL    =       (*)
          Reduction  Factor


          (*) Reduction  factor assigned depends on permit averaging  :
           •  period.  As determined earlier,


                                 30  Day -  --- R = 0.643      -

                 '   "CF/EL        7  Day.-  - - - R = 0.439

                                 1  Day -•  - - - R = 0.281


The only other inputs called for are the coefficients of variation of
stream flow, effluent flow, and effluent concentration, which have already
been determined.


The facing page illustrates the input prompts that are displayed when
the program is run, and the values entered in response to the prompts,
in this case for evaluating the 30-day permit averaging period.
                                     3-28

-------
                                 COMMENTARY
                                                          VII (3)
                                                          Revision No. 0
      DISPLAY AND PROMPTS

 POINT SOURCE - RECEIVING WATER
   CONCENTRATION ANALYSIS
                                          RESPONSE ENTRIES
 INPUTS:  COEF  VAR OF QS,QE,CE
         RATIO... 7Q10/avgQS
         RATIO... 7Q10/avgQE
         RATIO... avg CE/EL

   BACKGROUND  STREAM CONC  (CS)
     IS ASSUMED TO BE ZERO
 ENTER  COEF  VAR  OF  QS,QE,CE?

 ENTER  FOLLOWING RATIOS:
  	,7Q10/avg OS ? -. - -
        -7Q10/avg QE 7

        .avg CE/EL?  -
ENTER LOWEST,HIGHEST,AND INCREM-
ENT OF MULT OF TARGET FOR WHICH
% EXCEED IS DESIRED
7              •  •


ENTER LOWEST,HIGHEST,AND INCREM-
ENT OF MULT OF TARGET FOR WHICH
% EXCEED IS DESIRED
                                           1.5,  0.2,  0.7


                                                0.05

                                              •  3.0

                                                0.643
                                   This prompt repeats after the
                                   selected range of values has
                                   been computed'and displayed.
                                   It allows the user to be guided
                                   by output in selecting values
                                   and ranges for subsequent comp-
                                   utations.
NOTE:
                                          0.01, 0.06, 0.01.   ,

               •                         ' 0.08, 0.36, 0.04

                                          0.40, 4.0, 0.2


The manual analysis, presented earlier, computed the exceedance
probability and recurrence interval for specific stream concentra-
tion values.  The computerized computation generates these results
for stream concentrations expressed as multiples of the target
concentration (CL) that is explicitly assumed to result when

Effluent Concentration CE = EL   (the effluent limit)

Effluent Flow          ,QE = TJE"   (average QE)

Stream Flow            QS = 7Q10.(the design stream flow)
                                     3-29

-------
                                       VII (3) /.
                                       Revision No. 0
EXAMPLE COMPUTATION.- COMPUTER PROGRAM
              (continued)       '.    ~"~"
            PROGRAM OUTPUT

   ********************************
      RECEIVING WATER CONC (CO)
       PROBABILITY DISTRIBUTION
          AND RETURN PERIOD
    FOR MULTIPLES OF TARGET CONC.
     DUE TO POINT SOURCE LOADS  .
   ********************************
COEF VAR 	 QS =
COEF VAR 	 QE =
COEF VAR 	 .CE =
7Q10/ave QS =
7Q10/ave QE »
ave CE/ EL =
1.50
0.20
0.70
0.05
3.00
0.64
       STREAM  CONCENTRATION  (CO)
MULT OF
TARGET
(CO/CL)
0.01
0.02
0.03
0.04'
0.05
0.08
0.12
0.16
0.20
0.24
0.28
0.32
0.36
*»0.40
0.60
0.80
1.00
1.20
1.40
1.60
1.80
2.00
PERCENT
OF TIME
EXCEEDED
92.699
80.916
71,039
62.788
55.862
40.808
28.659
21.170
16.201
12.728-
10.206
8.320
6.875
5.746
2.650
1.399
0.804
0.490
0.312
0.206
0.140
0.097
RETURN
PERIOD
.(YEARS)
0.003
0.003
0.004
0.004
0.005
0.007
0.010
0.013
0.017
0.022
0.027
0.033
0.040
0.048
0.103
0.196
0.341
0.559
0.878
1.331
1.961
2.821
                  3-30

-------
                                                             (3)
                                                         Revision No. 0
                                COMMENTARY
- - -  This output is for a 30-day permit average period (Ra =  0.643)

       The range of values selected here is broad enough to facilitate
       construction of probability and recurrence interval  plots.

       Stream concentrations listed are in terms  of a ratio to  the'target
       concentration (CL).  In  this example, the  target  stream  concentpa-
      . tion is:        ,

              CL = 2.5                             ,  ;    •'•.'..•','•

       Actual  stream concentration is  this value  multip!ied by  the  listed
       value:  e.g., the  multiple  of Target (CO/CL)  =0.4

       Corresponding stream concentration is:

              0.4 X 2.5 =1.0

       Since the  acute-to-chronic  ratio  for pollutant  (P) is  6.25/2.50 = 2.5,
       acute exceedences  are reflected  by multiple  2.5.     .

       Probability or recurrence interval  plots can be constructed, simply
       by  plotting the  values listed in  the computer printout.

       Note that  the probability distribution of  stream concentrations
       deviates from log-normal (a  straight line) at the higher exceedance
       percent!les.        ~
                                    3-31

-------
                                       VII (3)
                                       Revision No. 0
EXAMPLE COMPUTATION - COMPUTER PROGRAM
              (continued)
   STREAM CONCENTRATION (CO) '(cont.)
    MULT OF
    TARGET
    (CO/CL)
PERCENT
OF TIME
EXCEEDED
RETURN
PERIOD
(YEARS)
2.20
2.40
2.60
2.80
3.00
3.20
3.40
3.60
3.80
4.00
0.069
0.050
0.036
0.027
0.020
0.016
0.012
0.009
0.007
0.006
3.977
5.507
7.509,
10.098
13.41.1
17.612
- 22.894
29.482
37.640
47.674
                  3-32

-------
                                                          VII  (3)

                                                          Revision No. 0
                                 COMMENTARY
     o
     o
     a
     ca
     tr


     UJ
     u

     o
     o
     ac,

     fe
          I 2   5 10    30  50  70    90

                % EQUAL OR LESS THAN
                                                          99   99.9 99-99
Figure 3-6-  Concentration  versus  probability for PDM-PS  computation
        10
     8
                                                    10 YEAR

                                              RECURRENCE INTERVAL
	ACUTE TOXICITY CONC.  2.5 — — ——•—^^n^.—, —	;	

    Ul
    u
    z
    o
    o
    UJ
    tr

    fc
       O.l
       .01
J—i  i I ittii.    i  ,  .
         -Ol  .02   ,05   .1
                          nl
J
                                    i MIII
                                                          J	•  • '  •
                .2    .5    I    2     5    10

                  RETURN PERIOD-YEARS
               20    50
Figure 3-7- Concentration versus mean recurrence interval for PDM-PS computation.


                                     3-33

-------
                                                         VII (3)
                                                         Revision No. 0
                  EXAMPLE COMPUTATION - COMPUTER PROGRAM
                                (continued)
To examine stream concentration effects for other permit averaging periods,
repeat the analysis, substituting the appropriate value .for the reduction
factor (R = CF/EL).                     ,                               .

The return period curves provide a useful  summary and perspective; however,
the evaluation can be performed without constructing the graph.  In this
case, the range of concentrations specified might (as shown below) simply
bracket those of principal interest.  In this case, a range of CO/CL from
0,5 to Sis selected, because the chronic limit (CL = 1), and the acute
limit to be exceeded no more than once every 10 years is CO/CL = 2.5.

The relevant portions of the output for the three permit averaging periods
are shown below:

                                     STREAM CONCENTRATION (CO)
       30-Day Average

        EF/EL = 0.643
MULT OF
TARGET
(CO/CL)

  0.50
  1.00
  1.50
  2.00
  2.50
  3.00
                                              PERCENT
                                              OF TIME
                                              EXCEEDED
          3.
          0.
 .818
 .804
0.252
0.097
0.043
0.020
                    RETURN
                    PERIOD
                    (YEARS)
            0.072
            0.341
             .085
             .821
            6.443
           13.411
                      1.
                      2.
       7-Day Average
               0.439
  0.50
  1.00
  1.50
  2.00
  2.50
  3.00
.1.717
0.272
0.069
0.023
0.009
0.004
0.160
1.008
3.957
12.149
31.819
74.364
       1-Day Average

       "CT/EL - 0.281
  0,
  1.
  1.
  2.
 .50
 .00
 ,50
 ,00
2.50
3.00
0.560
0.060
0.011
0.003
0.001
0.000
            0.489
            4.601
           23.866
           90.571
          281.076
          756.249
                                     3-34

-------
                                                         VII. (3)
                                                         Revision No. 0
                                COMMENTARY
 In this case a different averaging period would be selected than that
 based upon the manual computation.  Acute criteria exceedences have a mean
 recurrence interval shorter than 10 years for a 30-day permit average, so it
 would be rejected in favor of a 7-day average, which meets the guideline.  •

 Note that the exact computation using the computer program indicates a    .
 6.4 year return period for acute violations, compared with a 2.6 year
 return period estimated by the manual  approximation.  The manual  approxi-
mation tends to give conservative projections for the longer return periods
that are of interest, though differences vary depending on specific input
conditions.

Hence, there will be marginal cases where the approximate computation may
reject a 30-day average inappropriately.

On the other hand, wherever the manual  approximation accepts a 30-day
permit average as appropriate, it is safe to assume that the more exact
computation will not modify the choice.

For the site specific conditions assumed for the example analysis:

     o  Any pollutant with ah acute-to-chronic ratio of 9.5 or greater
        would, based on the manual  approximation, always be assigned a
        30-day permit average.

     o  The PDM-PS computation extends  this to pollutants with acute-to-
        chronic ratios of 3 or more.                            .
i   •  '   . •  ,\         ,      •                           .         .

NOTE:  EPA interprets any return period greater than 25 years as  being
       highly improbable.                 '  , .•
                                     3-35

-------

-------
                                                          Revision No. 0
 , -                        .  ...    CHAPTER 4
              .  RANGE OF EXPECTED VALUES FOR STREAMS IN U.S.
        As illustrated in Chapter .3,  the method can  be applied  to  any  site
 specific evaluation for which the relevant statistical  parameters are
 available or can  be estimated.  The  purpose of this  section  is to present
 a concise summary of the results of  such computations for the  range of
 site conditions that are likely  to be  encountered in  practice.  This
 chapter provides  such'a compilation  along three lines.   Section 4.1
 describes the basis for the input values selected to  provide a representa-
 tive range of site conditions, and presents  the results  of an analysis
 using these typical  ranges in  the methodology  described  previously.  The
 stream  flow characteristics .were determined  from an analysis of 180
 streams  and rivers;  treatment  plant  effluent characteristics are  based on
 analysis  of data  from over 400 POTWs.   The  results in,this section apply
 for  conservative  (nonreacting) pollutants.   Section 4.2 describes how the
 Information  provided  by  such an  analysis  can be used  as a screening .tool
 for  selecting  permit  averaging periods.   Section 4.3  presents results of a
 similar analysis,  except that it is  specific to oxygen depletion by
biochemical  oxygen demand (BOD)  loadings.  Section 4.4  extends the
analysis  for conservative pollutants to the special  case of streams that
are  highly effluent dominated, including those.with significant zero-flow
periods.

4.1    Analysis for Conservative Substances        •
                1         '         . •               >           •     .- ' •
       The review of stream flow and effluent statistics presented in-
•Appendix 8 indicates that the following ranges  are reasonable.   Effluent

-------
                                                          VII  (4)
                                                          Revision No. 0
 concentration variability,  (VQE)»  is  in  the  range  of  VQE  =  °*3  -  1.1."
 Effluent  flow variability,  (VQE),  is  generally  small  relative , to stream flow
 variability and,  therefore,  does not  greatly influence-, the computation.
 VQ£ s 0.2 is consistently used.  Stream  flow variability follows from the
 empirical  relationship  ofvqs and 7Q10/T45V For  a specified ratio,  the
 range of  VQS,  as indicated  by the data  discussed  in  Appendix B, is used.
 The ratio  7Q10/TJS" varies considerably.   A  representative range is  7Q1D/W
 * 0.01 -  0.25.  Finally, the magnitude of  the effluent flow  relative to the
 stream flow is  specified by  the effluent dilution. ratio: 7Q10/OJE\  A range
 from 7Q10/TJF =1-50 is chosen to represent  effluent dominated streams and
 large streams with small discharges.  A  10 year return period has  been
 selected as  the acute criteria violation frequency.

       In  order to compute the ratio  of the mean effluent  concentration to
the  effluent limit Ra = TJF/EL,' it  is  assumed  that  the permit -violation
 frequency  is one  percent.  The final  specification required  is the relationship
 of 7 and 30 day average effluent concentrations to the daily effluent concen-
tration coefficient of variation,  vcg.  Based upon the data presented in
Table C-2,  it appears' reasonable to expect that the 7-day  averages have a
 coefficient  of variation that is 0.8  of the daily values,  and that" 30 day
 averages have a coefficient  of variation of 0.6 of the daily values.  Thus,
the  reduction factors used are:
       Coefficient of Variation
           of Daily Values
Reduction 'Factor,  Ra
    a = 99 Percent
VCE
0.3
0.7
1.1
1-day
0.527
0.281
0.187
7-day
0.593
0.340
0.229
' 30-day
0.671
0.425
0.296
                                     4-2

-------
                                                           VII (4)
                                                           Revision No.  0
         The results of these computations are summarized  in  Figure  4-1  and.,
  given in detail  in Tables  4-1 to 4-4.   The  three  choices  for  permit
  average are shown.  Each group of bars  represents the  range in effluent
  concentration  variability,  vCE'   Each individual  bar represents a
  particular effluent  dilution,  7Q10/T4F.   Finally,  the length of each bar
,  represents the range that results  from  the  range  of stream  flow variability
  (7Q10/Q5"= 0.01  -  0.25) and  the  associated  coefficient of variation,
  VQS.  The  ordinate is the downstream concentration (in multiples of the
  chronic  criteria)  which has  a  10 year, return period.1

      ..  A number of features  are  immediately apparent.  For pollutants with
 an acute to -chronic  ratio of greater than 10, no acute criteria violations
 are projected over the ranges investigated, and 30-day average permit
 specifications appear to be sufficiently protective.  For. acute-to-chronic
  ' »                 "          ,        .                                  . - -•   f
 ratios of less than  10, site specific considerations are important.

        The results are most sensitive to. the stream flow parameter  7Q10/T33T,
 as can be seen from .the range covered by each bar.  For example,  the  last
 bar in the figure, 30-day permit averaging period, 7Q10/QE" = 50,  vCE  =  1.1.
 covers the range  from p= 0.9 to 4.6, corresponding to  7Q10/Q3" =  0.01 and
  VQS = 2-4.

        Following, in  order  of decreasing sensitivity, is the effluent
 dilution ratio:   7Q10/T3F.   A significant distinction can be  .found between
  The EPA is presently  considering  the  issue  of  allowable duration  and
•  frequency of exposure to  toxicity.  Based upon this work, duration and
  Trequencies used  as the decision  criteria may  change.  This guidance does
  not recommend any particular minimum  acceptable duration or frequency.
                                      4-3

-------
                                                          VII  (4)
                                                          Revision  No.  0
  8
       LEGEND:
                  NOTE'            .
                  HEIGHT OF BAR INDICATES STREAM
                  FLOW VARIABILITY (7QIO/QS)
          rr  r
         135  SO
o
o
*
C3L
 EFFLUENT
 DILUTION
(7QIO/AV6.QE)
in
o
  2 -T
      *CE'
             2^0.7  Z/CE*UO
                              HE30'7   ye.-I.IO
                                                                Ol M
                                                      ycE»0.3  VCE = 0.7
                                         ^CE=I-10^
    EFFLUENT LIMIT FROM WLA
    SPECIFIED-AS    DAYAV6.
EFFLUENT LIMIT FROM WLA
SPECIFIED AS  7 DAY AVG.
                                               EFFLUENT LIMIT FROM WLA
                                               SPECIFIED AS 30  DAY AVG-
   "INDICATES THE STREAM CONCENTRATION (CO) WHICH WILL BE EXCEEDED WITH A
    FREQUENCY OF ONCE IN TEN YEARS, EXPRESSED AS A MULTIPLE OF THE CHRONIC
    CRITERIA (CL).
Figure 4-1 - Effect of  permit  averaging period on stream concentrations
             for conservative  substances:  general analysis.
                                     4-4

-------
01
o-

CO
 •a


 c
 0)



 «*-
 01




 o

' c
 

-t->
 » en


             4J
             C •   •
             CO'
            <4-     I

            <4-    O
            UJ    CO
                     r-..    i ,_• x

                     o   ^ g'S

                     It

                     UJ
                     U          «  i
                     *    "  >»cn
                                        tn un co    «3 ^r ,   en ir> r-.
                                          ***"   «-••     •••     ••••
                                        OrHCVI    Or->cn
            t~    o «j >
            UJ    CO O «C
           <»>     I  r- X

           o    ^5^
           ii
           UJ
           u         •
           i    . I  >» C7>
                    (U
          UJ
                                                                                     CM CO
 cntncr>   cncMCM   i—< co CM •   CM CM

• o CM '«a- •   o CM «*   <-* r* *f -', ^4 CNJ .
                                                                                       CM
                                                                                     CM
                                       ^•^oo   f^-cor^   co-*«a-.  COVOCOOO.CM

                                       OCM^   o'r-Tco'   d^'co*    o^'oo*    ^^'co*
                                       CTi P~. ^
                                        »   a        ,
                                       o CM in   o CM
                                        COOCM    cnincn   cncorv    ,-(^^10

                                                   O i—I CO   Or-ICO    <-HCMCO
cr>


•ffili



a


i
        ca
        
            +j u. cr >
                                o
                                 •
                                o
             m
             o
                                                                         in
                                             in
                                             CM
                                                      4-5

-------
01
                as
                03
                VI
                0)
                o
               CD


               03
               Oi
               tn

               o
               u
               I
               o
               tu
               (U i
               in
               Q.t-1

                 O-
 cn
 'O CD
 C-t-
 CU 4->

< 2

 1

04



IX!

CQ
                                     i —  X
                             LU
                             O
                            (U
                            u-     i >, en

                            LU    co a «c
                            LU
cu
3
                                   1  >, 0)
                           LU    CO
                                 Q«S.t"
                           o

                            II

                            LU
                            O
                            > en
      f"- 
                                  i  >, en
                           LU    CO
                                 03
        •4-J  CO
   CO    T-  O*
  •(->    I—  ^
   (O    »i—

  •r-  O to  CO
  •P    -r-  cn
   «>    t_  c
  LU    to  ra
        > C£.
 0}  O t—I   «
 s_ ^-o- en
.<_> a. r~  >
co       »COO     CO CM CM    CTttOCM   'COCO^-

                    C3 t-l CM    O t-1 CM     O t—t CM    O* I-H CNJ    t-H I-H CVJ





                    !**• tO ^    CO VO tO     O to IO'    i—I CO P^    IT) CM O

                    O t—I CM    O i—I CM     t—( i—I CM   .r-t i—I CM    t—I CM CO





                    cnoi^i    O o CM    cocn«3-    tocoto    ococn

                    o CM co    t—i CM co    t-i r^ en    i—i CM co    CMCMCO








                    *o co t—i    vocoo    r^ t—4 t-H    cocot-H   •—IVOCM

                    O»--ioto    C0to«^-    cn'sj-to    ovoto   coc?>r*.

                    Oi-HCM    O»-HCM    O^-ICM    i—I »4 CM   ^Ht-HCM





                    o o t—i     cncno    t—< r^« •—•    CMOCM   tO'Q'^'

                    OCMCO     Ot-HCO    ^Hi-HCO,   t-J.CMCO   r-H CM  CO








                    CO r^ to     (»-. LO CO    CO CM CM    CO «* CM    O VO  CM

                    O ^H CM     Ot-HCM    O t-1 CM    O I-H CVJ    I-H ^J  CM





                    cncneo     coc».io    cn^to    envoto    t—i co  m

                   Ot-HCM     O t—i CM    Ot-HCM    O«-HCM    »-i,-4cM





                   o'cMt—i     cr>cT»cn  .otoco    o co co    cooco

                   •-4CMCO    O«-HCM    r-tt-HCM    i—I »-« CM    i-HCMCM








                   222    ooo    tooo    ooto    oioo
                   ooo    otoo    i~«oto    vocncM    tor^o

                   c^co^a-    t-Hi-HCM    o r-« t—i    oot~<    oo>~J


                                  8          8          8          8 •
                                Soi^H   Oe;'->   OCH'-H    OC£M
                                •1°-—   -J Q. n:   _j o_ IE    _j a. ;c
                                                o o: —•
                                                _! Q- 1C
                                                   o

                                                   o*
                                                               to
                                                               o
                                                                                        to
                                                                        to
                                                                        CM
                                                                4-6

-------




c
o
3

^
• c.

1
• £_
OJ
. O
o
. 'fc.
*•—
•X
t,
-M
'I '
s=
Q
*f»
•i^
(J
0)-
CU CO
M
II
 '
•^O ^3
0) Z>
•a: «-
i
CO
w
ca
*""•
*"? rH .•? *
rH Q1 ^
ll
LU
CJ .
P*. ns >
•S 9*
LU CO O «C

„
.''
f^ ' 1 r— X
• rH -r- a

II

LU
, o . .
* ~ I*". «* 2*
-M
' C
(U
- 13
r~ •/
<*- • >> o»
<«- . O n> > .
LU CO Q , o> ,.•
^J ° **
c
O)
3

V- O ro >
LU CO Q «C


co £• cr
•U r— ^

10 t. C
UJ fO (Q
> ce


TO s o"°"
 O <— 1 •
«- r- O- 01
4-> u. r*^ >
oo <;-
IO CM 00
O rH rH


r>« tn CM
CD rH CM


cy» cy> oo
O rH CM



«3 CM CO


•


r^. in CM
O rH CM



OS OS ^^
O rH CM




,co in o
O rH CM


en r^ co
O rH CM
•
.


C3 OS ttj^



o o o
000
* • •
CM co «d-
8
32 =



o
C5
^^ co os
O rH rH


co in co
"O rH CM


rH CM CO



VO CM 00





oo in CM
• '• A
O rH CM



O CO P"*
rH rH CM




r>. oo os
• « a
O <— 1 rH


oo m CM
O t-H CM




cr> r»* in



o o o
o in o
rH rH.CM
08M



in
o
C3
CT( CM O
O rH CM


o in in
rHr4cM
•

CO OS CM
rH rH CO



00 rH OS
O rH rH




CT> CO CO
O rH CM


r.
rH rH CM




CO CM OS
o • e
O rH rH


cr> co CM
O *-t CM




o in «a-
rH rH CM


in o o
P*. o in
« a e
O rH rH
of-



o
rH
0
O in rH CO OO CO
rH rH CM rH rH CM


CM 00 IO VO CM CO
rH rH CM rH CM CM
1

m co co I-H co P-.
rH CM CO CM CN4 CO



oo co os I-H m o
O r- 1 rH t— 1 ^H CM



1 1
o — i o ce. >— •
— J Q- 3C _1 Q, Z



in in • ,
rH CM
0 CO
4-7

-------
pi
                cr

                3

                C
                V)
                <0
                

                0)
                w
                c=
                O
                U
                t_
               •4J
                ta
                m.
                O
                0]
                OJ
                  O
                  T—I
                  cr
 *o o
 QJ I *
«c 2
 t
               UI

               ca
               >  zn
                             O)

                             r—            •
                             u-     i  >»  cn
                             <+_    O m  >
                             LU    CO Q «C
               LU
               U
               >    J  >>
                             > 0>
                            u-    O  ta  >
                            UJ    CO Q «C
                            CO     ,
              UJ
              o
                             (U
                     I  >, cn
                                    i  >, cn
                                   o 
                                   co a <;
                                      •u  on
                                
 r^. i—< 10    OOCOP^.    ocooo    cvji«o>     m co i—<

 C3 i—< rH    O t-l t—I.    i—1 •—I 5—I    i-H f—I i—I     r-< r—I INI





 00«*GO    O(OO:CMVOCvl    «iS-Cr>CO     COOJVO

 O i—Is—<    ^M>-HCsl     r-l «—I C\J    <-H i—I CSJ     •—I CM CM





 r-ICO«^    COOVO     VOOCft    CO«*O-COCf»PO

 t-li-HCVJ    ^HCSJCVI     t—I CM CVJ    i—I CM CO     CMCMCO








 r*. t—i <*    t~«.CMto     cn i—i vo.    ocovo     cgtor-*

 C3t—lt-4    O V-1 «—I     O «—I «—'     •—I i—I «-«     t—(  T-l »-4





 cocof»*cr>«5rco     C3?rer>     CM«OO     ^co*-*

 O i—< t—i     O «-H »—I     ^H^HJ-H     ?-l »—I CM     _H  ,-) CM





 O r--«—i     i-rir««CM     co^^^-     «^-ef»LO    COCMVO

 <-H<—3-    Q PO UO

 O «—< t—I     O i—I t-H     O t—( »-H    O >—< 5—)    t—I i—I rH





 COCO^O     COCOVO     CTICMVO    OCOVO,f-H-!3-f^.

 Ot—I •—I     O»—If—IO<—ii—I.   «—»a-co    ococo    I-^LOCO    covocr*

 «—< I—4 i—I     O »-H «—I    t—I 1-H »—«    t—I r—I i—I    r-4 r-< t-H








 ooo     ooo    unoo    oom     OLOO
 O o o     omo    r»«oio    VOCT>CM     tor»«C3
o o; >—i    o CK •—•
_i o- =c    _j S: ac

                                          8           8
                                                                                                      o en »~«
                                                                                                      _J a. 3=
   o
     •
   o
                LO
                o
LO
r—I
  •
o
LO
CM
                                                                   4-8

-------
                                                          VII (4)
                                                          Revision No. 0
 the effluent dominated streams, 7Q10/HE" <_ 5, and the large stream case,
 7Q1D/QF = 50. For the latter cases', the stream flow variability is a.more
 important determinant of the normalized downstream concentration.  Finally,
 the effluent variability,  VCE, affects the results by approximately  a
 factor of 2, all  other things being equal.
         \           •                    "           .                '
 4.2    Use As a Screening Too,!

        It is suggested that Figure  4-1  may  be  used as  a screening tool  to
 separate the cases  which  can be dealt with  immediately  from those for
 which  more site specific  information .is required.   For  the latter cases,
 the flow ratios,  7Q10/OT and 7Q10/QS" can  usually be  found quite  easily  so
 that a more specific  answer can be  found  in  Tables 4-1  to 4-4;   The final
 determinant,  vgs,  requires  a log-normal  analysis  of the stream  flow
 record.'  Since  this is  reasonably straightforward, a more refined  analysis
 Is  not excessively  burdensome and would serve  to reduce the range  of
 possible  values of  p,  from which the permit  averaging decision can be
made,        ,           '     '   .   .                  •     -.-     --..-'•

       As  ah example of such a  screening  analysis, consider the  hypothet-
ical case  of a..state'establishing permit  averaging periods for phenol.
Phenol has an acute-to-chronic  ratio of 4, so that stream concentrations
which  exceed a multiple of 4 times the chronic concentration  will not  be
accepted  (assuming that the acute criteria is not to be exceeded on a  daily
basis more often than once every 10 years).                       "

       Comparing the bars on Figure 4-1 with the multiple of.  P= 4, the
following conclusions relative to the permit averaging period can be.
                                     4-9

-------
                                                          VII  (4)
                                                          Revision  No. 0
 drawn.  For situations  with  an  effluent  dilution  ratio of -5 or, less
        a.    A 30-day permit averaging period will be selected whenever
             the VCE is  0.7 or  less.                            .
        b.    Where  VCE  s  !•!»  a. 7-day permit averaging period will meet.
             requirements  under all reasonable possibilities of stream flow
             variability (VQS). (The upper ends of the bars correspond
             to high  values of VQS.)
        c.    Even for effluent  variability as high as VCE =1.1, there
             will be  ma'ny  s.t reams where it would be appropriate to select  a
             30-day permit average, since only the upper end of the bars
             exceeds, a multiple of 4.

For an effluent dilution ratio 7Q10/3F = 5, the third column from the
right (VCE = l-.l; 30-day permit average} in Table 4-2 indicates that
only the highly variable stream flows approach violations using a 30-day
pejTnit average.  State  records could be examined to  determine if. the set
of streams under consideration (or a representative  set from Appendix C)
experiences  VQS in this range.

       A conservative decision, then, would be to select  a 7-day  permit
averaging period, although a site-specific assessment of  stream flow
variability or a restriction of VQS values could be  expected (in  most
cases) to support selection of a 30-day permit averaging  period.   -
                                     4-10

-------
                                                          VII (4)
                                                          Revision No. 0
 4.3    Preliminary Analysis for Dissolved Oxygen

        The Choice of permit averaging periods for effluent  limits of
 oxygen-consuming pollutants., such as BOD or ammonia,  is  a more complex
 problem than that addressed in  the previous sections. The  variations of
 the minimum' or critical  DO are  caused not only by effluent  concentration
 and dilution fluctuations, which  are addressed by the probabilistic
 dilution model, but also  by fluctuations  in reaction  rates  and other
 sources and sinks ,of DO,  such as  algal  production, respiration, and
 sediment oxygen demand.   Stream flow and  temperature  variations affect
 these parameters, the  latter also  determining  the  DO  saturation.  A  .
 comprehensive  probabilistic  analysis that would  include these effects as"
 well  is beyond  the scope  of  this  report. •

        It  is desirable, however, to  provide at least  a preliminary analy-"..
 sis for suitably  restricted, cases that are  amenable to analysis using the
 probabilistic dilution model.  The method to be employed makes use of
 the similarity  of the formula for critical DO deficit for those streams
 for which the simple Streeter-Phelps formulation is adequate, and  the
 dilution equation.  The principal  assumptions are (1)  a  single point
 source  of BOD is the only DO sink; (2) the stream flow,  geometry and
 reaction rates are spatially constant; and (3) the reaction  rates  are
                                           /           '   '   • .           ' -
temporally constant.  For this restrictive situation,  the critical or
maximum dissolved oxygen deficit (DG) is a function of the reaeration
rate (Ka), the BOD oxidation rate  (Kd), and the ultimate-to-5-day.  BOD
ratio.                                             .
                                     4-11

-------
                                                           VII  (4)
                                                           Revision  No.
                                                                           .
         The  Streeter-Phelps equation  can  be  solved for the critical or
                '       -               -                 '
 dissolved oxygen deficit , (Dc) :
                                    •
                            Dc = CE e F • , «J»  « P-                     (4-1)
 where:
          CE = treatment plant effluent BODs concentration.
           F = ratio of ultimate/ 5-day BOD.  Stream calculations are based
               on ultimate BOD; effluent criteria on  5-day BOD.
           «p = stream dilution factor = QE/(QS + 1QE).
           P = stream purification factor; for a BOD  'oxidation rate (Kd)
               and  stream reaeration rate  (Ka),            .
                        A
               P =  (A,)1"*; where  A  = Ka/iq

 tNote that  if the purification factor were constant  then  Equation  4-1 '
                                      ^                 '
 would be formally equivalent  to  the dilution equation analyzed previously.)
 One  remaining  difficulty 'is that  it is  not the critical DO deficit  (Dc)
 that is  of concern  but rather the  critical dissolved oxygen  (DOc) 'itself:

                               D0c  - Csat  - Dc                        (4-2)

which  is  a function of stream temperature  through the DO saturation
concentration, Csat.'  Hence,  the applicability of probabilistic, dilution
to the dissolved oxygen problem requires that the analysis be  restricted
to periods for which temperatures are essentially constant and fluctuations
in the purification factor (P) are small.
                                     4-12

-------
                                                           VII (4)
                                                           Revision No. 0
         An evaluation of this, latter effect can be made as  follows.   A
  relationship  between P and  stream depth,  H,  which .follows  from  Ka and Kd
  versus  depth  relationships  is  [3]:

                                P  -.H0.8'-    I"      •-.                  (4-3)

  and for many  streams,  depth is proportional  to  total stream flow,. QT, to
  a power H *0* with m = 0.4 - 0.6.  Thus,   -
                            Q"     n = 0.3 - 0.5-
(4-4)
 Consider Equation 4-1 for critical deficit.  Taking natural logs and
 applying the formula for the variance of a sum of independent random  '
 variables yields: .'
                        -   2   '     2        ->
                        -.°"lnCE   ""In*   + n
where  QT = QS + QE;   This  equation,  of course,  ignores  the fact  that'®
and  QT are correlated,  but the  point is that  n2  =  0.09 - 0.25 so  that
If the log variance  of  QT  is  comparable to  the  effluent concentration.
log  variance,  then. the  n2  term  is  not  a major contribution  to critical
deficit log variance; hence,  it; can  be neglected.   The  fact that dilution
(«p)  and total  stream flow  are negatively' correlated would  further  reduce
the  effect.

        Hence,' the key observation  is that if  it were possible to restrict
consideration to those flows .for which  VQS  *  VCE, then purification
factor fluctuations would not be very significant and probabilistic
dilution can be applied.  If these flows also correspond to periods of
                                     4-13

-------
                                                         .VII  (4)
                                                         Revision No. 0
 approximately constant  temperature,  then the two  requirements for applying
 probabilistic dilution  to  critical dissolved oxygen have been met.  For a
 site-specific analysis,  the  obvious  solution is to seasonally analyze the
 stream flow and  temperature  data and apply probabilistic dilution, making
 any necessary corrections  for  purification.factor variations.  However,
 for the general  case considered here, an alternate approach is.required.

        Consider, instead,  restricting consideration to that period of the
year during which  flows  are  low.  This period corresponds, presumably, to
the period  of time during  which 7Q10 occurs, and.includes the conditions
for which the WLA was performed.  Considering this period alone  signifi-
cantly  reduces the variability of the stream flows to be considered.  If,
in  addition,  it can be argued that these low flows tend to occur during
the same season each year, then the  temperature variation is less than the
annual  variability and will be less  significant as well. Hence,  for these
low flow periods, the assumption of  constant"? is much more realistic.

        The technical problem to 'be solved is to compute the.reduction  in
the average stream flow and coefficient of variation  when flows  are
restricted to  the low values for this restricted period.  We  restrict
consideration  to the lowest one-sixth of the total  population.   This
corresponds to an average of 2 out of 12 months in each year, and the
presumption is that this period recurs during the same months each year  so
that the temperature variation during this  restricted  period  is  small.
This simplification also assumes that the lower one-sixth'of the  daily
stream flows occur only in  the two  month period when  temperature  and
reaction rates are assumed  to be approximately  constant.
                                     4-14

-------
                                                           VII  (4)
                                                           Revision No. 0
         As  indicated  earlier, a  statistical analysis of actual stream data,
  stratified by month  or critical season, could be performed to provide
  actual  results and avoid the need for. this type of estimate.  However,
  data of this type are not presently available.  The estimation described
  here is performed in order to allow a preliminary analysis for BOD/DO to
  be made.

        The computation of the required stati stical  parameters , the  stream
 flow average and coefficient of variation  for flows restricted to the •'
 lowera-quantile of the  total  population*  is  straightforward.  For  log-
 normal  random variables,  it  can  be. shown that these conditional moments,
 denoted by  primes,  are:

                      US"1                                      '
                          =  Q( alnQS>-Za  )/Q(Za )     '              (4-6)
                                                                          .
                                                      --1            (4-7)
where Q(Z*) = Pr  Z > Z*  for Z, a standard normal random variable, and '
2a are the Z scores for the a-quantile which is the upper bound for the
flows being considered.  .For a = 1/6, Za > 0.967.  Table 4-5 presents  the
results.  These corrections, when applied to 7Q10/QT and VQS in the first
two columns of. Tables 4-1 to 4-4 adjust these parameters to  represent  the
low flow periods.  For highly variable streams, vQS and  therefore 
-------
                                                          VII  (4)
                                                          Revision No. 0
 TABLE 4-5 - Conditional  moments  for  the  low  flow subpopulation.

                               (a= 16.7%)
Coefficient of Variation Reduction
for .in Reduced
Entire Record Mean Coefficient of Variation
yQS W/TjS Z'QS'
, , .0.50
0.60
,0.75
0.90
1.00
1.25
• 1.50
2.00
3.00
4.00
0.450
0.384
• 0.306
0.247
0.216
0.158
0.120 "
0.0761
0.0389
0.0241
0.188
0.216
0.254 '
0.287
0.306
0.348
. 0.381
0.431
0.500
0.547
Th'is table provides a basis for a preliminary estimate of the average
stream flow and flow variability during critical  low flow periods,  rel-
ative to overall long-term characteristics.  For site-specific cases, the
actual values can be determined readily from a statistical  analysis  of
stream flows during the selected critical  period of the year.
                                     4-16

-------
                                                           VII (4)
                                                           Revision No. 0
  compressed from VQS = 0.5 - 4.0 to  VQS'  = 0.19 - 0,55,  so  that  the  sub-
  population of low flows fluctuates much less  violently than the-entire
- population, which includes the  annual  cyclical  variation as well.

         A 10 year return period  was selected for consistency with, the
  general  analysis,  but  since only one-sixth of the flow-population is being
  considered,  and  we  assume  that  no  DO acute criteria violations occur
  during  the  remaining higher flows, the exceedence probability to be
  applied  in the probabilistic dilution calculation is a 10/6 = 1.67 year
 .returnperiod.   Figure 4-2 and Tables 4-6 to 4-8 present  the results.

         In order to properly evaluate the computations, it is necessary  to
 realize that they apply to 10 year return period critical deficit ratios.
 To convert critical  DO concentrations to the deficit ratio  (p)  shown  by .
 the tables,  the DO standard (CL),  the DO saturation- (Csat) used in  the
 WLA, and the DO concentration  taken to  represent an  acute criteria  value
 are required. For most  reasonable combinations of these values, the ratio
 will  be  between  approximately  2.0'Vnd 2.5.  For  example, if  CS =  8,
 CL = 5,  and  acute DO =2, then  p=  2.0.  Alternatively, if these concen-
 trations  are  CS  = 9.0,,  CL = 6.0, acute DO = 1.5, then  (the acute-to-
 chronic deficit ratio)    p= 2.5.

       Appropriate permit averaging periods'are seen in Tables 4-6 to  4-8
to be strongly influenced by local  conditions of effluent  load and stream
flow variability.  Because of this, a general  statement on permit  averag-
ing period for effluent BOD/DO is not possible;  it must be selected  on the
basis of site conditions.
                                     4-17

-------
                                                         VII (4)   '
                                                         Revision  No. 0
8
7
_i
o
o
S5
z
o
p
CONCENT
H 4

1


". LEGL
" *-\ '^
1 1
•A

>0i
I NOTES
HEIGHT OF B/
> FLOW VARIAE
1 35
EFFLUENT
DILUTION
(7QIO/AVG.QE)
•m
••
.
f ''
•
•* - .
•»
.?; ;
• /
fLr ipl'
EFFLUENT LIMIT FROM WLA
SPECIFIED AS 1 DAY AVG.
\R INDICATES STREAM
MLITY C7QIO/QS)








f
>
\_
/;
••
4
*
y

fr
f
•>'
bu
«
^
V
I ;
V V
H
V
ES
T
i^
EFFLUENT LIMIT FROM WLA
SPECIFIED AS 7 DAY AVG.
*
.
. -
**
rTJ "
7^ .
IM , • '•" S **
x _, y ,
^ . jf '
* $ *
" ' - • * ' ^ ', "
tt i : 4 —
'.fi £_, |j
EFFLUENT LIMIT FROM WLA
SPECIFIED AS 30 DAY AVG.
  INDICATES THE STREAM CONCENTRATION (co) WHICH WILL BE EXCEEDED WITH A
    FREQUENCY OF ONCE IN TEN YEARS, EXPRESSED AS A MULTIPLE OF THE CHRONIC
    CRITERIA (CL).
Figure 4-2 -
Effect of permit  averaging period on stream concentrations
for BOD/DO.
                                   4-]

-------






to
II
' £
°,
• o-
1
0
• «J

, en
i-l . Q <
II .
UJ
O •
^ i ;>> cn
^^ (H ^
« «•
•W
CO
3
i- i >> en
<*— O (« >
UJ CO Q «C



'- i*- • i >> en
/**** cs ^c
II

UJ
<-> •

•M
C
, en
<4— O ra >
^kJ CO (CZ^ «^f





t •
co i >, cn
^^^ c3 *ic
n
UJ
o .
^ i >> en
!**• (O >
c5 .5
4->
. C
O)
3
C i >, cn
<4— O 
uj co a <;




CO
TJ
o oo
•r- ' O"
co aJ *
u a_

i * ^ w
cn o (c/>
O O*
oq 0-
 O*
U- 0 >
p- » '
*a m t/%
OJ Q. ^
-|__* p_v £3
CO r— r-4
< o-


«!f CO '.CM




tO Q 10


vo CM cn




«* r-. CM





.»*• cn «*
O O i-H




CO .-H CO
O «— 1 1— 4

'

,


•>f CO CO




^ cn *^*



to o vo
O «-H r- 4





CO O tO
«* to to





CO VO t— 4
«-H CM "3-




O O O
O O O
CM CO «3-


• i— I
O
' c
0

vo O vo




^^ CO ^^
O i— 1 f— I

CD VO IO
O t-1 CM



to o to
O i-H «-4




VO CM CO
O r-4 r-*




CO ^- CM
O r-4 CM






to o to
O 1—4 .—4



VO 1-H (>.
0 rH ,-4


vo CM cn
O 1— 4 t— 1





«-H CO CO
CO CO ^J-
.000




CO CM VO
CM «* VO
o o o

/.

o o o
o to o
i-H ^H CM


to
o
*
o
4-19
. fr* ^^ CO
O i—4 r-4



cn co CM
O i—4 CM

CM r-. cn
»-H i-H CM



^^» C^3 ^^
O i-H i— 1


\

CO CM 0
O r-4 CM




o «* to
i— 1 i-H CM






vb cn vo
O O 1-4



r*^ f-s co
O 1—4 t— 4


CO CM O
O i— 1 CM





tO r-4 CO
CM CO CO
o o o




CO VO CO
CO SJ- CO
o o o



to o o
I~«. O to
C3 r— 4 i— 4


o
o

; co co o
- O »-< CM

''

O VO «*•
•— I I-H CM

CO r- 1 r-H
i-H CM CO



CO CM CO
O r— 4.1— 4




CT> «3- CM
O i-4 CM




r-H CO" r**»
«— 4 r— 1 CM






r^r-i r^.
O i— 1 i— 1



co co cn
O r-H rt


cn «g- CM
O r-H CM





CM O% tO
CM CM CO
000




cn I-H to
co vo cn
000



o o to
VO CT\ CM
o o *-<


to
1—4
C3

. r-4 VO CM
r-4 t-H CM



^^ ^D i^^
r-4 CM CM

CO VO tO
rH CM CO



0 U? 0 .
r-l r-l CM




CM CO «5T
. i-H 1-4 CM




LO CM O
I-H CM CO






cn co cn
O 1— 1 1— *



O tO 1-4
1-4 1-4 CM


CM !-•» «3-
i-H r-H CM





cn to j-H-
i-H CM CO
O Q O



' i
tO CM VO
tO CO «-H
O O 1-H



o to'-
10 r-» o
O O i-H


to
CM
• o













•
3
• ' " c
•r»

^CJ
O)
CO
• 3
-u
• fHf
u
1 •!—
»^
J
as
T3
(U
c
o

T3
QJ
•a
(U

•r"
o


-------
   OJ
*-.  >
                   CO

                    II
4J
ea
                   o
                   CO
                  «*-
                   , 01
                             <*-    o  «  >
                             UJ    CO O «C
          o

          It

          UJ
          o



          4-3

          cu
                            UJ
         co

         o

          II

          UJ
          o
          5k

         4-9
          c
          cu
                            UJ
                                       >> at
                                         "
                                          cn
                                   o
                                   CO i
                                          0>
    >> Ol

  '  c§ .^
                                   I  >, 01
                                  2  
                                   CO
                                  •§
                                         to
                 LU


                 CQ
                            CO
                            o
                            4->
                            to
         to     o    Ico
         T    rr    tp-
                            o
                            us
 o

iZ

 o
         s
         o    co

        LZ    ID

         E    "tl

         cu    a.

        4->    r—
        CO    i—
               ,  co CM    too^    r~»ovo    cocMrv.     o^-cn

  o o i—i   . o T—11—i    o 7-4 »-H    o r-I '^4     r4^4r4



                                                                    •a
                                                                     cu

 tocn^-    P*.CMI^«.    cocMcn    cn^s-t—4    CMI*~.CO        "*
   **  *  *     •••     •••     e  •  A      •   •  •
 O O i—I    Or—I «-H    O i—I i—I     OrHCM    t—I r-4 CM

                                                                     O



 "•ocMr>.    co tri CM     oto«a-     CMCOVO    tocsicn

 O 1—I r-l    Of—I CM     rHrHCM     rHrHCM   'rHCMCM       4J
                                                                     cu
                                                                     Ol

                                                                     "3


                                                                     cu .

    CO CM  .tO O ««-     P««CT>tO    f^rHtO,   CTlCOl*^        o

 OOrH    OrHrH     OOi—4    OrHrH    O


                                                                    i»


    oico    vorHto     r^or^    co CM co    o to cn

 OOrH    OrHrH     OrHrH    OrH»-H    ^4,-J^       *"'




                                                                    (O

 tooto     P^ CM co     co  CM cn    cn «* o     CM r*». CM        ^

 OrHrH     OrHrH    O  r-H r-l    OrHCM     rHrHCM        O








 COOtO'    rHCOCO    tOrHCO     CMCntO     CTVtOrH
 -^toto    coco«a-    CMCOCO     CMCMCO     ,-HCMCO        cu

 OOO    OOO    cJcDO     OOO     O O G3        o

                                                                    •a
                                                                    cu
                                                                    •o

COtOrH    COCMtO    COtOCO     CT>rHtO    tOCMtO        Sj
rHCM«*    CM«srtO    CO«3-CO     COtOCntOCOrH        o

OOO    OOO    CDOO     OOO    OOrH.      gj



                                                                    o



OOO    000    tOOO     OOtO    OtO          -S
ooo    otoo    t^oto     tocncM    tor^-o

CMCO«3-    rHrHCM    OrHrH     OOrH    OOrH       Q
                                                     o
                                                      *
                                                     o
                to
                o
o
r-H

o
to
«—I
  e
o
to
CM
                                                                                                                    ' 
-------
<=>!



•
t—4
II
O
• r-H
Or
l"
o
•i—

cr
cu

i
c
_0
0
cu
r-«
cu
CO
.'•o
o
•r—
c.
en
c
•r—
en
03
cu
03
-u
•r—
1
-co
UJ
— J
CO
•
«— i i >» en
it
UJ
t^ t = •
>, i >» en
*j
• c
cu
C « >> en
•
>*• i >) en
f^ ^3 *C
II
UJ
0 •
». • >> en
cu
£ ci>»£
UJ CO C3 «$
•
co • i >, en
II
UJ
O •
> i >> en
4->
C
cu
£ o ^s1
CO
•o -
o co
co eu ^
0 Q.
•r—
4J 2 -
_U> O JCO
cu o
-U 3 t-H
o o o-
,t_
I
2
o co co
f— -O O-
U- 0 >
•r—
CU 'Q- JQ*
5 .£' 2
 o
U
CU
"cu
en
4-)
<"
'I
CO
ns
CU
>>
f-H
•r— '
"^3
QJ •'
«=•:,
CU
deficit ex
C3
O
0
•r"
4_>
E
                                                4-21

-------
                                                         VII  (4)
                                                         Revision No. 0
       *A table  for  the  effluent dilution ratio  (7Q10/Tir) equal to 50 has
not  been prepared for BOD/DO.  For small discharges entering larger
streams, it  is  likely that an effluent BOD limit determined from a steady-
state  WLA analysis  would be greater than the technology-based limit which
would  be used in the permit.  The use of the standard matrix table, which
would  show a higher pattern of violations, would tend to be misleading,
since  the computations  and the tables assume that the allowable effluent
concentration determined from a WLA becomes the effluent limit (EL)
specified by the permit.

       It should be emphasized at this point that the dissolved oxygen
analysis  presented  in this section is meant only as a preliminary applica-
tion.  There are, as yet, no verification examples that support the
applicability of a  probabilistic dilution/critical  deficit analysis.  It
has not been shown  that actual stream DO data conform to' the probabilistic
assumptions and simplifications used in this preliminary analysis.
Further,  it is well  known that th'e DO distribution in streams  cannot
always be described  by the simplest (Streeter-Phelps) model.  Upstream
sources of BOD and deficit are common, as are nitrification, algal  effects,
and sediment oxygen  demand.  A more comprehensive analysis would  b.e  •
required to incorporate these effects into a calculation of the effect of
selecting a permit averaging period.

4.4    Analysis for  Conservative  Substances  in  Effluent-Dominated Streams

       An effluent-dominated stream.is defined,  for the.purpose of this
analysis, as one in  which the effluent flow  exceeds the  design  stream  flow
                '            '         4-22

-------
                                                          VII (4)
                                                          Revision No. 0
                 \.
 (e.g., thef7Q10).  There are then two bounds to this  analysis.   The  upper
 bound is the effluent dilution ratio 7Q10/avg QE = 1,  which  was  the  lowest
Dilution ratio examined  in Section 4.1.   The lower bound  is  provided by
 the case where the design stream  flow is  zero (7Q10 =0).         , ;

        It should be recognized that as the degree  of dilution decreases, a
 WLA-based EL becomes increasingly restrictive.   When the design  stream flow
 is zero,  the effluent limit  must  equal the stream  target concentration (CL).

       While the degree  of effluent  domination  has a subsequent  influence
 on the magnitude of an EL assigned  in a permit, the screening analysis
 results presented  below  suggest that  in most  situations, a 30-day permit
 averaging period  will be  adequate for effluent dominated streams.

       The'results  of a broad hypothetical analysis of effluent dominated
 streams are  summarized in  Figure 4-3 and Table 4-9, using the format  used
 earlier to illustrate the  influence of permit averaging period,  effluent  "
 variability  and dilution ratio.

       o  The bars on the right provide the  upper bound;  i.e., the condition
          where 7QIO/avg QE = 1 (these results were also  shown  in Figure
          4-D.    •      • -     ;     ;          •            ..     ;   .,--  :   •
     .  o  The. bars on the left represent an effluent dilution ratio of    '
          7Q10/avg QE = 0.1,  that  is, where effluent flow  is  ten  times
         .greater than .design stream flow.  High variability  of daily flow is
          expected for such streams, together with  a very small ratio of
          stream design flow  to average stream flow.  The screening analysis
          assumes that the coefficient of variation ranges between VQS•-
                                    4-23

-------
                                                           VII (4)
                                                           Revision No. 0
      8
             LEGEND:
                           NOTE:
                           HEIGHT OF BAR INDICATES STREAM
                           FLOW VARIABILITY (7010/OS)
   O
   o
   M
   G3.
   |4

   ui
   o
   §3
 ! Q2 I 1.0
0.1  0.8
  EFFLU-ENT
  DILUTION
  (70IO/AV6.QE)
        I/CE'0.3   yeg=0.7
        EFFLUENT LIMIT FROM WLA
        SPECIFIED AS     DAY AV6.
               EFFLUENT LIMIT FROM WLA
               SPECIFIED AS  7 DAY AV6.
EFFLUENT LIMIT FROM WLA
SPECIFIED AS 30 DAY AVG.
      * INDICATES THE STREAM CONCENTRATION (CO ) WHICH WILL BE. EXCEEDED WITH A
        FREQUENCY OF ONCE IN TEN YEARS, EXPRESSED AS A MULTIPLE OF THE CHRONIC
        CRITERIA (CL).
Figure 4-3 - Effect  of. permit averaging period on stream concentrations for
             conservative  substances"in effluent-dominated streams.
                                      4-24

-------
 to


 cu
•o
 0)
 o
•o

4->

 0)
O)


•r—
 0>



 



 I



 I



UJ



.«*
                       0>
                 I  >> 0>
 uj    eo
                O  (T3  >
                   o <
 O

 II

 UJ
        I  >> 0>
<•-   •   i  >> en-

oj    co c2 *^c
CO    - I _' x




 II

 UJ
 CO  •        .
 >.     I  >> 01
         c
         OJ
         <«-     «  >> CT>
         <4-    O  «  >
         UJ    CO Q  00
             O •!-  en
            w    t_  c
                »    t
            CO    LU
            a> f- o&
            =3 +J .1-	
           i— 3 •!-» O
           t- r— OS i—I
           ». en o     co en o

 O O i—t    O O •—I    O O'«-4     C3 C3 i-4






 tOOCO    tOOCM    CO O i—I     Cnr-HCM

 O<—I-H    O~H<-4    C3^4r-H.     Or-5t-4



















g%g\ggg    ggg    S§§

CM«fl-LO    CM«a-lT)    CNJ,«^lrt    CM^IO
                     — i a- a:  -•  _ 10.2:
                                                CLO:    _i
                                                                                    
-------
                                                         VII '(4)
                                                         Revision No. 0
           2 and vgs s  5, and estimates a stream flow ratio 7Q10/avg QS =
           0.005, for this condition near the lower bound for effl uent-dom-
           inated streams.

       The conditions  under which the design stream flow is greater than
zero are listed in more detail in Table 4-9.  Results for several additional
intermediate effluent  dilution ratios (7Q10/QE = 0.2 and 0.5) are also
presented.  A comparison of results for an effluent ratio of 1.0 presented
here as an upper bound, and previously (Table 4-4 and Figure 4-1) as a .
lower bound will indicate that results are similar but not exactly the
same.  Ihe differences are due to different assumed values for 7Q10/QS and
the range of coefficients of variations used as inputs for the PDM-PS
model.                                                         •   .

       For the case where the design stream flow is zero, 7Q10 is zero and
there appears to be a  problem since 7Q10/Q5" and 7Q10/DT are both zero.
However, what actually matters is "OT and "QF.  Thus, in order to
evaluate these cases, the use of the actual "DT, TJE" and a small
7Q10 suffices since the computation depends only on QSVTjF and 7Q10 cancels
                                             i
out {Equation D-14).   Finally, the use of a small  7Q10/OTT correctly indicates
that the WLA is. done with QS = 0 (Equation D-15).  Thus, no problems arise.  .

       Screening analysis results indicate that in the case of effluent-
dominated streams, a 30-day permit averaging-period provides adequate  protec-
tion for pollutants with the acute-to-chronic ratios summarized below:
                                     4-26

-------
                                            VII (4)
                                            Revision No. 0
                                  When
                         30-Day Permit Average
Acute-to-Chronic          .  Is Adequate for
      Ratio            ,     Acute Protection

    3 or more                   Al ways

    2 to 3               Effluent variability is
                         relatively high, but
                         less than VCE - 1.1
                       4-27

-------

-------
                                                 VII (A)
                                                 Revision Wo..-0
                    APPENDIX A

Statistical Properties of Log-Normal Distributions

-------
•\

-------
                                                          VII (5)
                                                          Revision  No.  0
           .                .       CHAPTER 5     .        -   •
          .           . .       USES AND LIMITATIONS

        The probabilistic  dilution model  has  been demonstrated to be useful in
 selecting the appropriate averaging  period for discharge permits.  The method
 is easily adaptable  to  situations which  vary widely in  terms of stream and
 effluent  characteristics, data availability, and'policy-level assumptions    •
 used in the  analysis.   Although  the  example  in Chapter  3 of how to use
 the method is based  on  the typical WLA assumptions of 7Q10 as the design-flow
 and chronic  criteria as the effluent limit, the method  is easily adjusted to
 accommodate  other assumptions.

        The method is intended to  apply to pollutants for which the regulatory
 concern is at  the point of complete mixing and for which the toxicity can be
 evaluated  in  terms of the  total pollutant concentration.  The method has been
 applied to a  range of stream and effluent characteristics which typify the
 characteristics of streams and effluents in the United States.  The'results
 of this application are useful as a screening tool, by which the appropriate
 averaging  period for many  field situations can be  readily identified.
 However, pollutants whose toxicity is a function  of pH,  temperature, and
 hardness require site-specific evaluations incorporating these parameters.,

       There are also several  limitations on'the use  of  the method.  .One of
the technical limitations is  that the level of chronic protection  is based  on.
                        »                      '           •.'•-'";•
state-specified design flow,  e.g., 7Q10,  7Q2, etc., which may  be overprotec-
tive or underprotectiye for many  site-specific  conditions.   The  EPA is
presently considering the issue of allowable  duration  and frequency of
                                     5-1

-------
                                                        .VII (5)
                                                         Revision No. 0
exposure to acute as well as chronic toxicity.  Users of this manual are
advised to refer to Part A, Stream .Design Flow, of Book VI, Selecting Design
Conditions, when considering the choice of an appropriate, chronic exposure
event.  Book VI is currently under peer review and will be issued by the
Office of Water Regulations and Standards once the peer review process is
completed.   '                                         ,'•-.-'
       Modifications are required to compute the probability distribution of
30-day average-concentrations, as-required for chronic criteria compliance;
these would have to be investigated and verified in the field.         .
       The major shortcoming, of the log-normal probabilistic dilution model
is Its misrepresentation of the lowest stream flows, thus tending to overesti-
mate the probability of high stream concentrations.  The use of a seasonally
segmented approach could be investigated.
       The effect of serial correlation on. the return period specification
would also need to be investigated, particularly with regard to the duration
of criteria violations.  For example, a knowledge of the return period for
n-day successive violations could be compared'to the time scales of the
criteria themselves.  This would provide a direct link to the toxicity data.
At a less sophisticated level of analysis, the tendency of criteria
violations to cluster on successive days could be investigated to provide a
basis for modifications to the method.
       For pollutants whose toxicity is a function of such secondary vari-
ables as pH, temperature and hardness, probabilistic methods are essential  in
that it is not possible to rationally choose "critical" or "sufficiently
                                     5-2

-------
                                                          VII  (5)     :
                                                          Revision  Wo.  0
protective" values for these variables.  Arbitrary choices cannot be defended
in terms of the probability of criteria violations.  Methods for analyzing
these situations could be developed, following the logic of probabilistic
dilution and incorporating the additional random variations of the variable..

      ..The application of this method to dissolved oxygen has indicated that
the probabilistic method provides a useful  approach to the problem of
DO deficit.  However this work has only been a first step.  Probabilistic.
methods can be further developed to assess  the effects of DO fluctuations on
fishery resources and to provide a more rational  approach to advanced waste
treatment decisions.
                                    5-3

-------
"'.I

-------
                                                           VII (6)
                                                           Revision No. 0
                                   CHAPTER 6 .
                .  . '     . ':  .      REFERENCES    .

  1..   DIToro,  D. M.t Probability Model  of Stream  Quality Due to Runoff.  J.
       Environmental  Engineering, American Society of Chemical  Engineers,  Vol
       110,  No. 3, June 1984,  p.  607-628.

  2.    DIToro,  D. M.  and,Fitzpatrick,  J. J.,  Verification  Analysis of the
     •  Probabilistic  Dilution  Model., report prepared  for  EPA  Contract No.
       68-01-6275, U. S.,Environmental.Protection  Agency,  Washington, O.C.,
       1982.

  3.    Driscoll  & Associates,  Combined Sewer  Overflow  Analysis  Handbook for
       Use in 201 Facility Planning, report prepared for EPA  Contract No.
       68-01-6148, U.  S; Environmental Protection Agency, Washington, D.C.
       (1981).     '                             •    . _                '  ,

 4.    Hazen and  Sawyer, Review of Performance of Secondary Municipal
    ,   Treatment  Works. Draft  Final Report  for Contract 68-01-6275, Work
   ;   Assignment  No.  5, U. S. Environmental Protection Agency, Washington,
       D^C., December  1982.

 5.    Niku, Shroeder, and Samaniego,  Performance of Activated Sludge Process
      and Reliability Related  Design. JWPCF, Vol.  51, No. 12, December  1979.
.6.   Niku, et al., Performance of Activated Sludge Processes:  Reliability.
      Stability and Variability.  EPA  600/52-81-227, December  1981.
                                     6-1

-------
                                                          VII  (6)
                                                          Revision No. 0
7.   Haugh, et al., Performance of Trickling Filter Plants:  'Reliability

     Stability and Variability, EPA 600/52-81-228, December  1981.

                                                                           •       * _»
8.   H/droscience, Inc., Simplified Mathematical Modeling of Water Quality.        "'

     for U.'S. Environmental  Protection Agency, March 1971.
                                    6-2

-------
                                                          VII (A)
                                                          Revision No. 0
      . This appendix is intended to present a brief, simplified review of
 the statistical properties of log-normal distributions which characterize
 the important variables in the water quality analysis procedures used for
 this'-report.  It is designed to help the user without a formal  background
 in statistics to appreciate the physical significance of the statistical
 properties employed.  It is not the intent of this appendix  to  present  a
 theoretical  discussion or to provide technical  support for developing
 relationships or equations used in the development of the methods  employed.

 A-l, ,  General  Considerations                    '

        The factors  which  influence the concentration  of ,a pollutant in  a
 receiving  water  body are  subject to a  significant  degree  of  variability.
 This  variability results  in  fluctuations in the resulting  stream concen-
 tration, which  is compared with target  concentrations  such" as criteria
 or standards, and which provides a  basis  for decisions on  treatment ..
 requirements. The approach adopted  in this report  for examining the
 effects-of different averaging periods  on treatment plant discharges uses
 the concept "how much  — how often" as  a basis for such decisions.  It is,
 therefore, essential that statistical aspects be incorporated into the
 methodology even though they may add complexity.

       The standard statistical parameters of a population of values for
 a random variable which are used as a concise means of describing central
tendency and spread are:  .                   .

       Mean:       ( nx pr x) the arithmetic  average,   x" defines  the
                   average of the available (usually limited)  data  set;
                                     A-l

-------
 Variance:
Standard
Deviation:
                                                   VII  (A)
                                                   Revision  No. 0
             MX denotes the true mean of the total population of
            variable x.  Twill be an increasingly better approx-
            imation of  MX as the size of the sample (the number
            of data points) increases.

            («Tx) *>y definition, the average of the square of the
            differences between individual  values of x and the
            mean (x").  The greater the variation in the data, the
            higher the variance:
                   2     (xi-x)2 + (xe-lT)2 +
                      — - - '
                                         N
            (.ox) another measure of the spread of a population of
           'random variables; by definitioni the square root of
            the variance:
Coefficient of                           ,   "
Variation:       ( vx)  is defined as  the  ratio  of the  standard
                 deviation (  ox)  to the mean (  Mx):

                             vx  = ax/. MX
                 It  is  the principal  measure of variation used in
                 the analyses described in this report.  The
                 coefficient  of variation is a  dimension!ess
                 quantity and is  thus freed from any dependence on
                             A-2

-------
                                                          VII  (A)
                                                          Revision No. 0
   •_    .• <                the specific dimensions used to describe the

                         variable (e.g., flow rate, concentrations, etc.).

                         High coefficients of variation reflect greater

                         variability in the random variable x.



        Median:          (~)  This is the value in a data set for which


                         half the values are greater .and half are lesser.



        Mode:            The "most probable value" -- more of the individual
/         •                                                       •    ..     ^
       .  .                data points are at this value (or are within  this

                        ^interval) than at  other values,or ranges.  On a

                         frequency histogram,  this is the highest point on
   •            /                   '                           •
    ;              .       the graph.   The mode  has no real  significance in

                         the calculations in the methodology  employed.



        Comparing the statistical  properties of  different  data sets  pro-

 vides a  convenient,  concise way of  recognizing  similarities  and  differ-

 ences.  This could not  be  accomplished  simply by "looking  at  the data"

 where reasonably large  data sets  are involved.   These statistical proper-

 ties  convey ho information concerning  frequency, or"the probability at '  '

 which any  particular value or  range of  values in the total popul ation.will
                              ,              •          '        ,- t       •
 occur. . This essential  item of  information  is provided by  a knowledge of

the type of distribution,  technically,  the  probability distribution

 function (PDF).




A-2.   Probability Distributions                          ^  -



       There are  several different" patterns which characterize the distri-


bution of  individual  values  in a large  population of variable events.



r  /••'":'•-.'•'-•'        A-3          .           •,

-------
                                                          VII (A)
                                                          Revision No. 0
 Most analysts are familiar with the normal  distribution,  in which a
 histogram of the frequency of occurrence of various values  describes
 the familiar bell-shaped curve (Figure A-l(a)). When the  cumulative
 frequency! is plotted'on probability paper, a straight  line is  generated
 as in Figure A-l(b)*
                                                   ,                   "L
        Many variables,  particularly those which are  important in  water
                      '                             •                         *.
 quality applications, have been shown  by a  rapidly accumulating body of
 data to be represented  by or adequately approximated by a log-normal
 distribution.  A log-normal  distribution has  a  skewed frequency histogram
 (Figure A-l(c))  which indicates an  asymmetrical  distribution of values
 about an axis defining  the central  tendency of  the data set.  There is a
 constraining  limit to lower values  (sometimes zero)  and a relatively small
 number of rather large  values  but ho upper  constraint.  Point source
 effluent concentrations  [1,2],  and  pollutant concentrations in combined
 sewer overflows  and  separate storm  runoff [3,4], are parameters which are
 usually well  characterized  by log-normal distributions.  In general, daily
 stream flows  are  satisfactorily approximated by log-normal distributions
 E5,6].   Scattered data from a number of  unpublished sources suggest that
 receiving water concentrations are also  log-normally distributed.   Stream
 flows  and concentrations are currently being examined from this  perspec-
tive.   A log-normal distribution appears as  a straight line on  log/proba-
bility  paper  (using cumulative frequency) as shown  in Figure A-l(d).   In
this report natural (base "e") logs are used throughout.
 Cumulative frequency is the relative frequency (or probability)  of
 values being less than or equal  to a specific value.
                                     A-4

-------
                                                  Vll  (A)
                                                  Revision No. 0
  tu
  <

 UJ
 u.
 O

 UJ
 O
 e>
 <
       CG)
      MAGNITUDE OF VALUE (X)


              NORMAL
I   '10      50      9O
      PROBABILITY  .

 % LESS THAN  OR EQUAL
                              99
                                    O
                                    I-
                              a:
                              UJ
                              en
                              ca
                              O
                              a.
                              O
                                .
                             ox

                             u.^E
                             o-r
                             UJ
                             a
                                   MAG-NITUDE OF VALUE  (X)
                                         LOG-NORMAL
I    IO      5O      9O   9.9
      PROBABILITY
% LESS THAN OR EQUAL
Figure A-1 - Probability distribution. .
                                A-5

-------
                                                          .VII (A)
                                                          Revision No. 0
 A-3.   Relationship Between Distributions                 .

        There are circumstances when two different types of distribution
 can begin to look similar — so that either one will  provide  a  reasonably
 good approximation of the probability distribution of a particular data
 set. For example, as the coefficient of variation becomes  smaller and
 smaller, approaching .zero, log-normal  distributions begin  to  look more and
 more like a normal  distribution.  Figure A-2 shows a  series of  histograms
 for log-normally distributed populations, all  having  (arithmetic) pop-
 ulation means  of 100, but with different coefficients  of variation (v )
 as  shown.  As  discussed  above, smaller values  of  v approach  a  normal
 distribution.
                                                                        , .,   •
 A-4.    Properties of Log-Normal  Distributions                         -   •

      •  Figure  A-3 summarizes the pertinent  statistical relationships for
 log-normal  probability distributions.   The mathematical  formulas shown
 are based on statistical  theory, and permit back-and-forth conversions
 between  arithetic properties (in which  concentrations,.flows,  and loads
 are reported)  and the log of the'variable (In which probability and  frequency
 characteristics are  defined).

       Normalized plots of probability versus the magnitude of a variable
 expressed as a multiple of the mean are presented in Figure A-4 for
 log-normal distributions. These plots present a family of curves reflecting
the effect of coefficient of variation on probability  of .occurrence of
events of specific magnitude. These, plots can be used  directly in the
                                     A-6

-------
                                                     VII  (A)
                                                     Revision Mo. 0
 o
 z
 UJ
 ID
 o
 UJ
 or
 u.
MEDIAN 7


MEAN /xx
COEFFICIENT OF

VARIATION

    Z>x=,'0.25
UJ
•3
a
UJ
or
u.
o
•z.
UJ
13.
O
UJ
cc.
u.
                                       X
                                       ,  I


                                   COEFFICIENT OF ;

                                   VARIATION
                                  COEFFICIENT OF

                                  VARIATION

                                      l/ = 2.5
.        .     -••  •        •   .   •       ;X       • •.--•  -    •

                            RANDOM VARIABLE   .  .

Figure A-2 - Effect, of coefficient of variation on frequency, distribution.
                                 A-7

-------
 Frequency
                  Arithmetic Space
  —Mode
     — Median x
           Mean P.
        X
x is a random variable
                 Definition or' Terms
•x	 Random Variable .•	
 u... 	-. - Mean	
                cr 	 Variance
                 /\
                                                           VII (A)
                                                           Revision No. 0
                                              Log  Space
                                          Natural  Logs  (base e)
                                                                 Pr
                o"x	  Standard Deviation ......... o^
                vx 		  Coefficient of Variation ... (not used)
                x  ..'	Median
                 Relationships  Between Statistical  Properties
                 	In  Arithmetic and Log Space
   °x
     x « exp [nlnx]
         >/exp (
  lnxj
                                      Tnx s ln
Figure A-3 - Pertinent relationships for log-normal distribution.
                                    •   A-8  .

-------
                                                         VII (A)  •
                                                         Revision No.  0
               10
UJ
2
UJ
X

u.
o
UJ
_J
a.
»-
_j
o
               1.0
                                    COEFFICIENT
                                         OF   ..
                                    VARIATION,
              O.I
                0.1  0.51 2  5 10 2O

               iOcr	'	——	
                                                           i  t
                            50  70   90    98   99.8 99.99
          UJ
          X
          Hv
         UJ
         _J
         a.
                                                COEFFICIENT OF
                                                  VARIATION,
                                                        1.0
                                                        O.8
               O.I 0.5 I  2  5 10 20    50  7O   90   98    99.8  99.99
                      PERCENT LESS THAN OR EQUAL  TO
Figure A-4 -  Cumulative log-normal distribution-
                                     A-9

-------
                                                          VII  (A)  .
                                                          Revision No. 0
 analysis methodology  and  permit  direct  determination of frequency for
 events of any specified magnitude  with  a  known or estimated coefficient of
 variation.

 A-5.    Standard  Normal Tables   .

        For normal  (or log-normal)  distributions, probabilities can be
 defined in terms of the magnitude  of a  value, normalized by the standard
 deviation. 'This technique is used in the calculations of the probability
 of exceeding  specified receiving water  concentrations in this analysis.
 Standard normal tables can be obtained  from any statistics textbook [8,9].
 Table A-l  presents the standard normal  table to provide a convenient source
 for the analyses used,jn this report.   Table A-l lists the probability for
the interval betw.een 0 and the ;value of Z listed.  Thus, it represents the
probability that a value will be less than or equal  to the selected value
of Z.
                                     A-10

-------
TABLE  A-l - Probabilities for the standard normal  distribution.


  Each entry in the table  indicates the proportion of the total area under the
  normal curve to the left of a perpendicular raised  at a distance of Z
  standard deviation units.
  Example:  88.69 percent of the area under a normal curve lies to the left
  of a point 1.21 standard deviation units to the right of the mean.    .   -
Z
0.0
0.1
0.2
0.3
0.4
0.5
0.6
0.7
0.8
0.9
1.0
1.1
1.2
1.3
1.4
1.5
1.6
1.7
1.8
1.9
2.0
2.1
2.2
2.3
2-4
2.5
2.6
2.7
2.8
2.9
3.0
3.1
3.2
3.3
3.4
3.5
3.6
3.7
3.8
3.9
0.00
0.5000
0.5398
. 0.5793
0.6179
0.6554
0.6915
0.7257
0.7580
0.7881
0.8159
0.8413
0.8643
0.8849
0.9032
0.9192
0.9332
0.9452
0.9554
0.9641
0.9713
0.9772
0.9821
0.9861
0.9893
0.9918
0.9938
0.9953
0.9965
0.9974
0.9981
0.9986
0.9990
0.3993
0.9995
0.9997
0.9998
0.9998
0.9999
0.9999
1.0000
0.01
0.5040
0.5438
0.5832
0.6217
. 0.6591
.0.6950
0.7291
0.7612
0.7910
0.8186
0.8438
0.8665
0.8869
0.9049
0.9207
0.9345
0.9463
0.9564
0.9649
0.9719
0.9778
0.9826
0.9864
0.9896
0.9920
0.9940
0.9955
0.9966
0.9975
0.9982
0.9987
0.9991
0.9993
0.9995
0.9997
0.9998
0.9998
0.9999
0.9999
1.0000
0.02
0.5080
0.5478
0.5871
0.6255
0.6628
0.6985
0.7324
0.7642
0.7939
0.8212
0.8461
0.8686
0.8883
0.9066
0.9222
0.9357
0.9474
0.9573
0.9655.
0.9726
0.9783
0.9830
0.9868
0.9898
0.9922
0.9941
0.9956
0.9967
0.9975
0.9982
0.9987
0.9991
0.9994
0.9995
0.9997
0.9998
0.9999
0.9999
0.9999
1.0000
0.03
0.5120
0.5517
0.5910
0.6293
0.6664
0.7019
0.7357
0.7673
0.7967
0.8238
0.8485
0.8708
0.8907
0.9082
0.9236
0.9370
0.9484
0.9582
0.9664
0.9732
0.9788
0.9834
0.9871
0.9901
0.9925
0.9943
0.9957
0.9968
0.9977
, 0.9983
0.9988
0.9991
0.9994
0.9996
0.9997
0.9998
0.9999
0.9999
'0.9999
l.OCOO
0.04
0.5160
0.5557
0.5948
0.6331
0.6700
0.7054
0.7389
0.7704
0.7995
0.8264
0.8508
0.8729
0,8925
0.9099
0.9251'
0.9382
0.9495
0.9591
0.9671
0.9738
0.9793
0.9838
0.9875
0.9904
0.9927
0.9945
O.S949
0.9969
0.9977
0.9984
0.9988
0.9992
0.9994
0.9996
0.9997
0.9998
0.9999
0.9999
0.9999
l.OOQO
0.05
0.5199
0.5596
0.5987
0.6368
0.6736
0.7088
0.7422
0.7734
0.8023
0.8289
0.8531
0.8749
0.8944
0.9115
0.9265
0.9394
0.9505
0.9599
0.9678
0.9744
0.9798
0.9842
0.9878
0.9906
0.9929
0.9946
0.9960
0.9970
0.9978
0.9984
0.9989
0.9992
0.9994
0.9996
0.9997
0.9998
0.9999
0.9999
0.9999
1.0000
0.06
0.5239
0.5636
0.6026
0.6406
0.6772
0,7123
0.7454
0.7764
0.8051
0.8315'
0.8554
0.8770
0.8962
0".9131
0.9279
0.9406
0.9515
0.9608
0.9686
0.9750
0.9803
0.9846
0.9881
0.9909
0.9931
0.9948
0.9961
0.9971
0.9979
0.9985
0.9989
0.9992
0.9994
0.9996
0.9997
0.9998
0.9999
0.9999
0.9999
l.OCOO
0.07
0.5279
0.5675
0.6064
0.6443
0.6808
0.7157
-0.7486
0.7794
0.8078
0.8340
0.8577
0.8790
0.8980
0.9147
0.9292
0.9418
0.9525
0.9616
0^9693
0.9756
0.9808
0.9850
0.9884
0.9911
0.9932
0.9949
0.9962
0.9972
0.9979
0.9985
0.9989
0.9992
0.9995
0.9996
6.99.97
0.9998
0.9999
0.9999
1.0000
1.0000
0.08
0.5319
^0.5714
0.6103
, 0.6480
5 0.6844
0.7190
0.7518
0.7823
0.8106
0.8365
0.8599
0.8810
0.8997
0.9162
0.9306
0.9429
0.9535
0.9625
0.9699
0:9761-
0.9812
0.9854
0.9887
0.9913
0.9934
0.9951
0.9963
0.9973
0.9980
0.9986
0.9990
0.9993
0.9995
0.9996
0.9998
0.9998
0.9999
0.9999
1.0000
1.0000
0.09
0."5359
0.5753
0.6141
0.6517
0.6879
0.7224
0.7549
0.7852
0.8133
0.8389
0.8621
0.8830
0.9015
0.9177
0.9319 ,
0.9441
0.9545
0.9633
0.9706
0.9767
0.9817
0.9857
0.9890
0.9916
0.9936
0.9952
0.9964
0.9974
0.9981
0.9986
0.9990
0.9993
0.9995
0.9997
0.9998
0.9993
0.9999
0.9999
1.0000
1.0000
                                        A-ll

-------
                                                         VII (A)
                                                         Revision No. 0
A-6.   References
1.     Niku, et al ./'Performance of Activated Sludge-Processes and Reli-
       ability Based Design."  Journal WPCF, Vol. 51, No. 12, (December,
       1979).-   •

2.     McCarty, et al ., "Reliability of Advanced Wastewater Treatment."

3.     EPA Water Planning Division, "Preliminary Results of the Nationwide
       Urban Runoff Program," (March 1982).

4.     Mancini, J. L., "Methpds for Developing Wet Weather Water Quality
       Criteria."  Progress Report, June 1981; EPA ORD Grant No. R806828010,
       Cincinnati.

5.     Chow, V.T.  "Handbook of Applied Hydrology."  Mc-Graw Hill, New York
       (1964).          '                 -

6.     Linsley, et al., "Hydrology for'Engineers." Mc-Graw Hill, 2nd
       Edition, (1975).                     ,

7.     Hydroscience, In., "A Statistical  Method for'the Assessment of
       Urban Stormwater."  USEPA, EPA 440/3-79-023,  (May 1979).

8.     Benjamin,  J. R.  and C. A. Cornell, "Probability, Statistics and
       Decision for Civil  Engineers."  McGraw-Hill,  New York, (1970).

9.     Johnson, R. R.,  "Elementary Statistics."  Duxbury Press, North
       Scituate,  Massachusetts, (1980).
                                     A-12

-------
                                          VII (B)
                                          Rev'fsfon No. 0
                APPENDIX B    ,      .

Field Validation of Log-Normal Distribution

          and Related Assumptions  .

-------

-------
                                                           VII  (BJ
                                                          .Revision No. 0
         This appendix presents a discussion of several technical issues and
  assumptions which are necessary to the use of the probabilistic dilution  .
  model to guide selection of permit averaging periods.  This discussion is
  organized  in two sections:  the first provides a justification for the use
  of the  probabilistic dilution model in the method; the second, provides a
  discussion of several key assumptions.
                  ,        •                                   - •     «•* i
  B-l.    Use of the Log-Normal Distribution

       • A relatively simple and straightforward analysis is made possible
 by the'assumption that each of the input, variables is log-normally dis-
 tributed and independent.  The appropriateness of these assumptions  and
-.their implications are discussed below.

        A basic feature of any random time series  of numerical  values
 is its probability distribution  function, which  specifies  the distribution
 of values and  their frequency of occurrence.   More detailed characteriza-
 tions which account for seasonal  trends  and day-to-day correlations  are
 also possible,  but at  minimum the  univariate  probability density function
 is required.  An  examination onflow data from a  number of streams indi-
 cates that  the  data  can  be  reasonably well represented  by  a log-normal
 distribution.   Figure  B-l summarizes  an  examination  of the adequacy  of  a
 log-normal  distribution  for  daily  flows  of 60 streams with long periods of
 record.  The actually observed 10th  and  1st percentile low flows are
.compared with the  flow estimated by  a  log-normal  distribution.  The major
 important discrepancy occurs  at the  lowest flows  where the predicted
 distribution is lower than thai actually  observed.   The most likely cause
                                      B-l

-------
                                                      VII  (B)
                                                      Revision  No. 0
      UJ
      .J
      UJ
      o
      a:
      UJ
      CL
      ID
      H
      O
       UI
      «*•
       o



      UJ
      UJ
      o
      or
      13
      l-
      O
                 =  10 PERCENTILE
             10'
             10'
             10'
              -I  s  t  titiMt?   t.titntri   i  tfttttri	;  tit
10
!0°
!0
                                                  I0
                    LOG  NORMAL APPROXIMATION (cfs )
                 «=  i PERCENTIL£
               10           !0W          10'        ^lO*1         10"
                   LOG NORMAL  APPROXIMATION  (cfs)
Figure B-l:   Evaluations of log-normal distribution   for stream flows.
                                 '8-2

-------
                                                     VII (B)

                                                     Revision No. o
the prince of a base strea™ f,ow whi
                                    which does
  Log probabiHty plots of treataent


    are .,,«,t«tid in Figure B 2 f
                                             in sane cases.
                                      p,ant efHuent fl
         ,3

  Uted do»stream
          r: rr
  for upstrea™ and eff,uent  f,ow

      concentrations  I

               T         9"
                                                    as .„ as
PercentiUs
                f
                              B-3

-------
                              in
                                                                  m
  °l
    *
   o
  =

   c

2 •*-
•»* ta
  " T"
*-« >
»—« 03
£» £C
§;
it  ;
ZL
                              3



                              a



                              8

                              Q
o
       21
I      -
           SSI
                              01

                              p-



                              B
      9



      a-



      a
      a
                       en
                              or
                               »
                              -H
                                                            CT

                                                            a
                                                  HN
                                     a
                                     a
                                                             cu
                                                             in
                                                             •rl
                                                               '  .'P
                                                             H
                                                             li

                                                             r   ;
                                                                                 —a
                                                                  B"
                                                                                      CO
                                                                                      o
                                                          tu
                                                                                  ..Si
ri



a



a

9
                                                                      v-
                                                                      GO

                                                                      O
                                                                      (T
                                                           tn
                                          a
                                          tu
   (1/6UJ) 009
                                                       (I /6tu) QOO iN3nidJ3
                                                                                                 I

                                                                                                 w
                                                                                                 c
                                                                                                 o
                                                                                                  1^
                                                                                v> o
                                                                                i- O
                                                                                                  (O 0)
                                                                                                 ja >
                                                                                                  o c
                                                                                                  s_ c
                                                                                                . a. o
                                                                                CM

                                                                                 I

                                                                                CO
                                  B-4

-------
o.
u»
>

^








- ,



'
' •(
Q • *"

**
• • . *
• ' ' •
: j
• ' • ; -
* . . •
a*
* . • a
a «

°: - '.
D • B
* * - • a
• *
' * »•
•
*t ' "
»•
f * • "
•s :
i ** i
utti » » i mm f t* f suit tit
«l •* C
> 0 O C
•1 »-• *H • T-
(I/fin) no



,
•













*i
> c
•« *•

• . • : - ' ;. •'• . ;
: •
• '* «
•: J
' : ' •
• *
, **
• •
' •
' *
.
":.
•
.
•• .
'-. • • = •

' - -
* • • ' *
• > •
i • .'•••' !
UTt i ( f nitti i i \ inn f i ",
01 -* o
> O 0 C
•< »-< •,_!.•.,.

-.


0
o»
i

x™'
^ :
o 3
m
<
CD
O
cc
«4


>
^ '







"•• i
M
O
^
•|
l
4J
C
03
"a..
+j
i










^
c
• , •
; "
• . ^
* • •»
*
- • . *
1 « • - " **
• - . ' 'a
• •• '
• _
• . '•
5 . . ' •

• ^
• •
. ' ;^ . " * «
:: -
mi»i « i Wto « « t Vntt iii "
> «< «« 0
> o o c











1 A
> c
• «
' : ' • - :
. ' 4
* ' 4
e , ' «
*
. •
.... ,
• ' a
*
* ' " ^
•
*
•
* , ^
* . ~
•• •
1. j
itiil f 1 l mtt« i i f bllM ill
» ei >• a
3 O 0 C
! • • j
. r
o •£•
°* - 'U
*^ M- 4->
• ^- ^
.1 • .2 >•
Si II
03 j
03" i
' « |
2 c


>
(I/On)
(I/fin)  jo
                   B-5

-------
*
TABLE B-l - Comparison of
* *
Median
Location


North Buffalo Creek, NC
Jackson River, VA
Haw River, NC
Pigeon River, NC
Mississippi River, MN
VII
: Revi
(B)
sion No, £ - '
observed and computed downstream concentrations^2) .
(50th Percentile) Concentrations
Model . Observed
Variable- Prediction Quantile


BOD (rag/1) 9.7 10.0
COD (mg/1) 51.0 59.0
TSS (rag/1) 16.0 15.0
BOD (rag/1)' 6.0 5.3
TSS (mg/1) 15.8 13.6
Color (PCU) 110.0 100.0
BOD (mg/1) 2.0 1.7
COD (mg/1) 23.8 22.0
BOD (mg/1) 3.7 3.8
COD (mg/1) 85.0 78.0
"NHs (mg/1) 1.0 1.1
»
Confidence
Limit of
Observed
Qjjantile
8.5 - 11. 0
47.0 - 66.0
12.0 - 22.0
4.2 - 5.0
10.0 - 17.0
90.0 - 130.0
1.5 - 1.7 .
19.0 - 26.0
3.0 - 5.1
65.0 - 87.0 •
1.0 - 1.2
95th Percent! le Concentrations

North Buffalo Creek, NC
Jackson River, VA
Haw River, NC
Pigeon River, NC
BOD (mg/1) 31.0 22.0
COD (mg/1) 120.0 97.0
TSS (mg/1) 15.8 13.6
BOD (mg/1) 18.1 15.6
TSS (mg/1) 41.6 32.0
Color (PCU) 324.0 330.0.
BOD (mg/1) 4.5 4.7
COD (mg/1) 43.0 46.0
BOD (mg/1) 8.7 7.6
COD (mg/1) 186.0 229.0
20.0 - 33.0
82.0 - 129.0
10.0 - 17.0
13.0 - 20.0
30.0 - 40.0 !
300.0 - 410.0
3.2 - 5.6
33.0 - 53.0
6.4- 9.4
188.0 - 233.0
Mississippi River, MN
(mg/1)
3.5
4.3
3.2 -   5.0
                                     B-6

-------
                                                          .VII (B)
                                                          Revision No". 0
 distribution of quantiles, are also listed.  In all  but one case, the ""
 computed quantiles are within the confidence limits. '

        Thus, there is no statistical  evidence to reject the computed  quan-
 tiles as not being the true quantiles of the observed concentration distribu-
 tion.  This is strong statistical  evidence that indeed  the  log-normal
 probabilistic dilution model  is  representative of actually  observed down-
 stream concentration  distributions for the 95th percentile  at  least.

        The  11 data sets used  in  the verification  analysis were examined
 for cross correlations between effluent, flows  and  concentrations.  The
 observed ranges  in correlation coefficients  have  no  significant impact on
 the computation.   -Correlations between  stream  flow and  effluent load  for
 a point  source are not expected.   Upstream concentrations are not employed
 in  the comparison  of  permit averaging  period effects, so that any correla-
 tion  between  stream flow and concentration is not  relevant to this analysis.
 Modifications to the  probabilistic  dilution model  computations are avail-
 able  for use  in situations .where cross correlations must be considered ,[!].

       The influence  of, possible deviations from the assumed log-normality
 of the upstream and effluent flows and concentrations upon more extreme
 quantiles is unknown  at present due to lack of larger data sets that  encompass
these extreme quantiles.  However, the quality of the alternatives to  and  the '
simplicity of this model argue strongly for its use in the present context  of
describing comparative differences in  water quality impacts.
                                     B-7

-------
                                                          VII (B)
                                                          Revision No. 0
 B-3.   Appropriateness of Assumptions       -            «        .
                                                         .,  '       '
        We have chosen to ignore the' seasonal and day-to-day correlation
                               ••'.'.•        '••-.-
 structure of both'stream flow and effluent behavior in order to simplify
 the characterization of each variable.  TTie consequences of this simpli-'
 fication are discussed below in more detail, but it should be pointed out
 that trends and correlations do not invalidate the use of the log-normal
 probability distribution function to characterize the frequency of  occur-
                          »,'...
 rence of flows and concentrations.  Trends and day-to-day correlations
 affect the time sequences with which certain values occur, but not  their
 long term frequency of occurrence.  This is judged to be an  acceptable
 penalty to be endured when compared to the simplification achieved.  If a
•more refined, site specific analysis is required, then a seasonal breakdown
 of the data, with the appropriate means and standard deviations for each time
 period, can  be generated  and the analysis performed as described below.

        The consequence of a possible serial  correlation can  be approxi-
 mately quantified as follows.  If, in  fact,  the  serial  correlation  is  such
 that 10 consecutive daily violations always  occur when  one violation
 occurs, then the  proper percent! le to  consider is not  0.0274 (10 years)
 but rather 0.274  (1 year  return  period).  The degree to  which  the 10 year
 return period concentration  is overestimated can  be estimated  by comparing
 the ratio of the  10 year  to  the  1  year stream concentrations which  are  .
 computed  without  regard to  serial  correlation.

        The ratio  of the 10 -year  return period concentration to  that for
                                     B-8

-------
                                                          VII (B)
                                                          Revision No. 0
 some other return period can be computed for log-normally distributed
 random concentrations by:                   •
            x
 where
                    EXP[(Z10yr - Zxyr)
          lnc          = log  standard  deviation  of  stream concentrations  (C)
                                                                 \
            yr>  C10.yr  = z  score  and concentration  corresponding to a
                         10 year  return  period

        Zx yr, Cx yr    = Z  score  and concentration  corresponding to an
                         x  year return period   • '   .                    .

        Table B-2 summarizes  results for' a range of values for coefficient
of variation of stream  concentrations.  Clustering tendencies of 5 and 10
are examined as approximations of the degree of serial correlation which
might exist.  If clusters  of 10 occur, the comparison is between 10 and 1
year return periods as discussed above; for clusters of 5S  the comparison
is between  10 and 2 year return periods.  On the basis of this analysis,
the water quality effects presented in Chapter 4 for various permit
averaging periods may overstate the 10 year stream  concentrations  by
approximately a factor of  1.5 to 2.0.            •

       Until stream and effluent  data  can  be analyzed to define  the serial
correlation  structure and the methodology  modified  to incorporate  it,, the
results presented  in  Table  B-2 should  be interpreted  to  indicate with the
following possibilities:
                                     B-9

-------
                                                        .VII (B)
                                                         Revision No. 0
TABLE B-2-- Approximate overestimation of 10 year return period  stream
            concentration by ignoring serial  correlation.
Variability of
Stream Concentration
Coefficient
of Variation
(*c)
0.5
1.0
1.5
2.0
Log
Sigma .
( °lnc)
0.4724
0.8326 .
1.0857
1.2686
Ratio of Stream Concentration
At Indicated Average Return Periods
10 Year
. to 1 Year
(CIQ/CI)
1.4
1.8
2.1
2.4
10 Year
to 2 Year
(Clo/C2)
'1.25
1.50
1.65
1.80
               Cl,2
                     =  EXP
                            lncl
    (10 year Return Period)  = 3.456
"Li  (lyear Return Period)   = 2.778
7-2  (2 year Return Period)   =2.
                                               996
                                    B-10

-------
                                                         VII  (B)
                                                         Revision No. 0
       o  Stream concentrations indicated by the methodology used in the

          report to recur on average for 1 day every 10 years would, if.

          they actually never occur except in clusters of 5 to 10 days,

          have return periods of 50 to 100 years.


       o  Conversely, for the same clustering assumptions, the stream

          concentrations that occur at 10-year intervals should be; 50 to

       ,   70% (1/2 to 1/1.5) of the 10-year concentrations projected by

         . the report methodology.


B-4.  References


1.    DiToro, D.M., "Probability Model  of Stream Quality Due to Runoff."
      J.  Environmental  Engr. ASCE,  Vol. 110,  #3, June 1984 p. 6.07-628.

2.    DiTqro, D.M.  and  Fitzpatrick, J.J.,  "Verification Analysis of the
      Probabilistic Dilution Model"  Report prepared for EPA Contract No.
      68-01-6275, U.S.  Environmental  Protection Agency, Washington,  D.C.,
      (1982),   . ,
                                     B-ll

-------

-------
                                          VII (C)
                                          Revision  Ndi  0
                 APPENDIX C

Characteristic Values for Input Parameters

-------
••1

-------
                                                          VII (C)
                                                          Revision No, 0
        The results  reported here represent an attempt  to .develop character- •

 istic values  and  ranges  for stream flow and effluent variability.  These


 values and ranges have been extracted  from the results of  published


 analyses,  and are used in  Chapter 4 to evaluate the influence'of the


 permit averaging  period  on typical  receiving water conditions.   These
         '                '                '            •   '            ' !'  ' •

 values are provided  for,  effluent flows (Section 1),, effluent coneentra-  .   -

 tions (Section<2), and stream  flow (Section 3).




 C-l.   Treatment'Plant Effluent  Flows
                                         • s     '  '   •              '      '


 .  '     A recent study [1]  analyzed  several  years of performance  data from


 approximately 400 secondary treatment  plants  in  8 different process  ,


 categories.   Average plant  effluent  flows  ranged from .0.002 to 82 MGD.     .


 Table  C-l  summarizes the coefficient of variation of treatment plant

 effluent flows.                .

                     • i                      i      .       • .     • ' .


 C-2.   Treatment Plant Effluent Concentrations          "'



        Data on the variability of effluent  BOD5 and total  suspended

 solids  (TSS)  from municipal biological  treatment plants are available •from

 several sources. . Niku, et  al. £2] provide analysis results for 37 acti-

vated  sludge plants  which show the coefficient of variation of effluent


 BOD5 concentrations  to  range between 0.34 and 1.11 for individual

plants.  The median  of  the individual plant.values was  0.635.   The EPA

research report [3]  on  which the foregoing was based  reported  a  mean


coefficient of variation  for 43 activated sludge plants using  a  variety of

processes.   Daily  effluent concentrations were found to be  well  represented
                                     C-l

-------



VII (C)
Revision No. £
TABLE C-l - Coefficient of variation of daily effluent flows, vgE.
• n
Process Category
Trickling Filter
Rock
Trickling Filter
Plastic
Conventional Activated
Sludge
Contact Stabilization
Activated Sludge
Extended Aeration
Activated Sludge
Rotating Biological
Contact
Oxidation Ditch-
Stabilization Pond
Number of
Plants
64
17
66
57
28
27
28
37
Range
Individual
0.06-
0.16 -
0,04 -
0.06 -
0.11 -
• 0-12 -
0.09 -
. 0.00,-
For Median of -
Plants All Plants
0.97 0.27
0.88 0.38
1.04 . 0.24
1.35 0.34
1.32 0.34
1.19 0.31
1.16 • 0.31
0.83 0.31
C-2

-------
                                                           VII (C)
                                                           Revision No. 0
  by a log-normal  distribution.  The mean of all  plants analyzed had co-
  efficients of variation of 0.7 for BQD5 and 0.84 for TSS.

         Two recent studies  have extended the analysis of effluent  concentra-
  tion variability, and, report  coefficients  of variation  of  BOD5  and TSS
  for 7-and  30-day  averages  as  well  as  for daily  values.   Results reported
  by Hazen and  Sawyer  [1]  provide the .basis  for the summary  presented.in
  Table C-2  as  well  as the two  other sources  cited in  the  table.  An analysis
  of  the performance of  11 trickling filter plants by  Haugh, et al.  [4] produced
  the results summarized by Table C-3.

        Based on available data, a single'representative value for coeffi-
  cient of variation of effluent concentrations cannot be defined.  The most
Appropriate characteristic value will  be influenced by process category,
 effluent concentration averaging period, .and the pollutant  in question
  (e.g., BOD, TSS,  etc.), as  well as  individual plant  differences.  The
 computations in this  report are performed using  a range  of  values  esti-
 mated, to encompass most of  the conditions of interest.

 C-3.  Stream Flow              ;     '   ,               •      ,. '

       Figure  C-l  provides  a basis  for estimating the  coefficient of
 variation of daily stream flows  on the basis  of the ratio of 7Q10 to -
 average  (7JS) stream flow.  These flow values  are usual!/  readily avail-
 able.  The  relationship shown  is .derived from a set of flow measurements and
 statistics, which has been developed for a sample of 130 streams in  various  "•'
 areas of the country [5] and is summarized in Table C-4,  along with addi-
tional details  on the location of the stream gages used.   The  ranges
                                     C-3

-------
of
J

rf—*-*1 O
On- •>
"— » «/> C
*-«**" °
— • QJ 4J
s» cc n
•r-

ra
u-
o
•«J
c
a
•r-
O
51
u-
o
u

ra
•5
cu
s
1
to
0
ra
&n
0
u-
a>
a.
4->.
c°"
ra
a.
c
fll
rj£
^^
ra
CD
t_ .
>* *""^
f*
•° S
"c *-
8 ?!
Si"
" 2
0 E
« ^- O
»*^ S™
In **~
§ LU

J ^
1
CM
O
LU

<:


^•^
>^.
cr
E
^_*
(—

[ <
C
cu
3

LU


ro •
a (/>
4— 1 O»
o o >
* CO «£
-r^ ^"
c o
cu ••->,•
•r- 4-) ra CO
O ro O O>
•f- .r- 1 >
<4- ro
CU > V)
o >> 
^
O
ca

4->

0 0 >
* CO «C
4-3 C
•r- 4-> ra to
o ra a o>
.^- .r- I >
<>- ro
cu > to
o >, cu
0 ,->3
ra ra


c
S




c- to

-Q <4- C
S O ra
S3^ ^
—









t"
o
o*
ra
0
to
/»!
cu
o
o.
a.
in
CM

o

v
0
CO
• 0


0
in
CD

CO
in
CM






CM
•
0

o
CO
o


o •
CO


o
vo
CM





to







^
o
o
az
c_
a;
4-J
P>_
uZ
a>
c
r"~
^£
U
'l—

.0

o


in
0


in
to
o

•*
i-H






O
co
•
o

in
CO
0
0


o
in
O


o
•
i—4





*£






• o
•r™
CO
^O
^•"
n.
&.
4-*
»^^
Lu
C31
C
^^
^^
u
1—
'• "-'
in

o


O
. to
0


in
03
O

CO
•— I






o
9
o

. in
m
d


in
to
p


CO
•
*— i





VO
to







•o
(U
4^
ra

o
^c

ra
c
o
4J CJ>
c -o
(U 3
> 1 —
c to
o
o

0
VO

o


to
0


o
o

03
co
«— 4






o
0
• o

o
m
0


0
VO
C3*


VO
CM
f— 4





• in







c
o

j_> cj
ra C7>
N T3
•«- 3

^ co
ro -a
4J CU
C/) 4J
ro
4J >
0 ••-
ra 4->
4-9 U
C «=C
o
C_3

0
CO

C3


m
0


in
VO
o

03
a\






in
^
o

0
to
e'
0


o
o


CM
•^





CO
CM










 •—
ra CO
Ig
"O ra
O) >
"O •!—
C 4->
cu u
4J «=C
X
LU •

in
co

o


0
m
o


o
0

CM
m
l-H






m
CO ^

o

in
C3


O
VO
0
O


O
1— 1





CM









t^«
ro
(J
cn
o

*o
•i- C_
ca cu
•M
O) (J
C ra
•r™ ), y
4-> C
 -Q
0 §
o
c_
r~- oj
CM ro
CM
C
3
0
.c

to
CO 3
r™"
ra








"O
1

c
0
N
,H.
.^K
ra
4-5
OO














-












,
*"
.c
ra
01
c
•o
•r-
•p
ra
4- i—
o a.

c c
O ra
•r- '^
4-> T3
rO CO
^
CU 4—
0 0
•a c
t- ra
ra QJ
^3 ^^
ra
CO
II

LU
o


in
1J-
to
ra
23 '



C-4

-------
                                                           VII (C)
                                                           Revision No. 0
  Table  C-2 (Cont.)
                       Chemical  Precipitation/Settling*

                 Pollutant                    Coefficient, of Variation

                    ,Cr                                  .99
                    £u                               .   .60
                    M6                '                  «57
                    Mn                                  .84
                    Nl                          .        .81
                    Zn                                  .84
                    Tss    .                             .66
                          Pharmaceutical. Industry2

                                Coefficient of Variation
        Plant Number           BOD     (n)TSS"
12015
- 12072
12026 • '
12036
12097
12098
12117
12160
12161
12186
12187
12136
12248
12257
12294
12307
1.01
.97
.95
.74
1.08
1.37
.70
.92
.55
.71
.21
1.02
.58
.64
.93
1.55
»•• '
' 46
392
44
366
222
24
39
. 34
249
54
12
110
50
56
56
39
.85
.63
.49
1.12
1.21
1.52
.81 .
1.11
.99
.50
.26
1.16
.55
.92
1.25
1.34
195
395
53
364. >
249
25
51 '
32
355
54
12
111
52
56 •
50
'38
                    14'of 10-18-83 memorandum .from H. Kahn to:E
                         t0 Data and Analys1s of the Combined".
2From preliminary descriptive statistics generated on pharmaceu-
 tical  data by SRI International, 11-12-82.
                                     C-5

-------
                                                         VII (C)
                                                 .:       Revision No. J3


TABLE C-3 - Effluent concentration variability for trickling filters
            (from reference 4).   .
,
Mean for 11 plants (mg/1)
Coefficient of Variation (median of
individual plant values):
Daily Values
7-Day Averages .
30-Day Averages
BOD5
29.6
•
0.39
0.35
0.31
TSS
29.3

0.55
0.31
0.26
shown reflect the bulk of the data in the sample of stream records which

were used.  However, a relatively small percentage of streams will have

coefficients of variation which fall outside the indicated ranges.  The

statistical analysis was performed for the entire period of record.

Results in some-cases may be distorted, if flow regulation works were

installed on the stream sometime during the period of record.



C.4.  References                ,
1  .  Hazen a*Hl Sawyer, "Review of Performance of Secondary Municipal
      Treatment Works."  Draft Final Report for Contract 68-01-6275, Work
      Assignment No. 5, U.S. Environmental Protection Agency, Washington,
      D.C., (December 1982).

2.    Niku, Shroeder, and Samaniego, "Performance of Activated Sludge
      Process and Reliability Related Design."  JWPCF, Vol. 51, No. 12,
      (December 1979).  •                                                              J

3.    Niku, et al., "Performance of Activated Sludge Processes:  Reliability;
      Stability and Variability."  EPA 600/52-81-227, (December 1981).                 i

4.    Haugh, etal.  "Performance of Trickling Filter Plants:  Reliability,
      Stability and Variability."  EPA 600/52-81-228.  (December 1981).            .'-   r

5.    Driscoll & Associates, "Combined Sewer Overflow Analysis Handbook
     , for Use in 201 Facility Planning."  Report prepared for EPA Contract
      No. 68-01-6148, U.S. Environmental Protection Agency, Washington,                I
      D.C. (1981).                                                                     4


                                     C-6                     .                          I

-------
                                                       VII  (C)
                                                       Revision No* 0
o
r-
O

u.
UJ

-------



ol
« *
3.
CJ 4-
*•— ' t/J

•— *>
4—4 CO






'-






* .









CO
0
•JM*
•*•>
CO
r
(U
5 * ^J
U
re
t_
*tj
j::
o
5
o
r"~
^4—
•


CU
I—
•4J
CO
O
*
^>
§
ys
3
«O

1

*s3*
i
o

UJ
CO
J2











3
O ^"^
r- CM
ft )^^

E *-.
re co
cu <«-
r- u


,
c

re
• h














c
o
•r—
4_d
re
u
2

CO
CT
re
CIS























2 CM
Or or
r-. t-4
sU
M


CD-
CM
CD-

CD
i— i
r-.




c"o-

lo-

re CM
Qj >•— i
•* 5.



^— ^ .
&.
re
zs
£_
O

^J*
S



















CO
*f™
0£
CO

re
4_5
CO


•
£
to
C3 CO
=D re1


in co cr> o r*.
en co in vo
CM t-4 i—4 t-4
.
. CO VO «*• CO
to O VO *— 4
O O O O CD


VO CO V-.' CM r-<
^- in i— i o cr>
•-I CM •— i. CM *-4
«3- co ^* vo r-»
CO O CO «-4 .— 1
o o o o o
CM »— 4 CD en
O i— 4 CO CM
l-H O 1-4 O O




^" ^J ^? ^O CO
CO vo co vo r*.
O i-l
en CM o en co
** ^ ° "* "2

o co !••» CM «a*
LO i-» m ft o
CM f-4 «S- •— I
1— 1
UJ
• LU
•2: A re
CO 2='
LU A, —
^*^ CHr^* "^ 'A
re *f— 21 c
A CO > 3
JZ^tf >> « 0
CO CO CO E 4J
re 3 c re co
OTO 4-> JZ "—
.re c -r- e_  £_ 4->
 > CO 0
•r- .OJ-r— > CO
r^ re /^^ •!<— >p—
co 01 a.
J= ^ CO
co co re c_ •
re 3 i- co t_
CDT3 J= 4-> CO
re c o co
r— co re >•> •
. o en 1-4
in o co o CM
r^ vo in en co

o en •— i <3* i— i
I-H »— i in en en
t— 4 i-H '.-
^ •
«c t_ -f-
-1— -^ CO "4— t-
Q£ O > 4J CO
O -r- CO >.
re t_ ce co •>-
3 CO 3 O£.
<— O^
co 4-> re • jc
re co o t- u
2: . . .
21
A ""3
co >- 2:
C712T
•a •
.,- A CO
i— i— CQ i— re
o c; T- >
CO 3= t-
A «,— QJ
* • r^T\ vx ^
• • re re •!-
• • LU O O£








c.
c- co
CO -^ >
> CO -r—
.,- t- C CO 0£.
C£. CO  > C_J ^i
cn-r- -i- o
3 o; as r— re
re r— co
^3 O O -i- C
CO T- •»- ^tf CO
C 4-> CO CO J^
•i- e o 4-» o
"•• .M g^ - as o z

!— 1— >- >- -D
o o 2: 2: 2:



c o o o o
O O 0 O 0
o m in in o

CM CM CO VO r*.


I^^HIJ^

co en in co CM
1—41—4 t— I

^r vo vo vo vo.
O O O O ir-4


co in r^ i-^ ^^
CD co 1—4 en P**.
in co CM en in
en vo O ^j" r^
O O t— 4 i— 1 1— 1
vo co CM en r*-
r»» co in t— i <— t
O i— 1 i— 4 r-4 CM




o «=f «a- en t-«
• CM co r-» «a- 1-4
CM VO CO CM CO
r-» CM vo o co
i—4 CM CM CM t-4
vo •— 4 vo i** ^r
CM -sj- VO r-4 —4
CM CO

>- >- UJ
-a 2: 2: a

A A *"* A
* to o c
CO r— f— • O
•— re •!- t- en
•r— U. > CO C
«-r~ ^.^ ,— •
•*£ t- O £=
t_ o re 4J ^=
re Or— re •»-
U_ 0 0 = 3





C_ '
CO

•r—
C£ 3£
c. -a
c at co
re > c_
4-> -i- C_5
•f— r— O£
g r— • >x Q|
re t- -s«i co c
cc ^^ c co i—
t- t_ 3
• C_ CO 0 >»
CQ > CO -C C
re > i — re

2:03 2: 3 ca

••3 >- >- «C uj
2: 2: Z 0_ £2



0 O O O O
O CD O O O
LO in o in in
co o LO en t-4

en CM co «a- co
CO «=1- -3- -3- «3-

—^——

-------
o|
o
.;:
VI

>


















.














T
c .
a
•
LU
CQ
^C
1 —
0
i— I CM
o- o*
^^ «™l
o

cr
. ' '

o •—
p— CM
ii i— i
(O CO
« en to co
«3- co r^ CM in

P^ CO CM CM in
CO i— 1 •— 1 •— I CM
. O O O 0 O
»-4-p>. in 1-4 in
CO P-- -!S- •— I CM
O O O O CD


en co in to en
CO P"^ <3" to * "5f CM i-r4


£^Sc?
••2; •> •> «
re *••• en re
r— *» p— C «r»
r— re ai -i— +J
•O co ex t- E
re re s cv r—
3 i— i O OO Q


N>
0)
 > -K: o
•r- .r- c_ o en
as os  t_ re
re re T- ea 4->
c cnos c
co °a re
re -i- c: o .c
. £Z G Q. O a
.3 C7) U 4->
cr c o 4- +J
CO O -C -I- 01
<^i=<3c%i3

2= 2= z: CL 2

o o o o o
o o o o o "
in in m o in
o 1—4, en .co in
O i— 1 CM •!!• in
in in in in un
fl— 4
C
CMO 0
i— i •— i • -: CM
• - :
CM en co
r-* O O O O
CO O CO CO CM
p** co en p** co -
. O O i— i «— I
to to co i-i en
O CO i— 4 P»> O
*- ' O O O O
«3- to P-»
CM CO r-l
i-4 O O O O.


CM in in co «—<
co p^» p**> vo co
«* co co o -=r
o en to co to
in co i— 1 1— i in

><
o en •
0 3 £•
-— «S * **** i— •
si c co •!—
« O C >
* >•> •*-> ••- >>
re -!-> C 4-> t_'
8S&£-£


i.

CO <1>  OJ -i£
Q. > O) t_ O)
re -T~ c. o or
+j QS o • e_
nj • as O
Q- j_> (1) t_
C OJ O C
•  re 3
CO 3.J3 0 0-
-t-> re to a)
• «3 £_ 3 a.
2= 0- 0 1— 0

Q O Q > cC
s s s 5: ^

o o o o o
O 0 O O O
o o o o o
to »-i i — r--.in
co en en •— * I-H
in in in to to


'.t»j ^j-
p

VO r"-<
C3 O 0
f"H ^^ d^
en to 'S-
CO i-l
to -« 1-1
to o co
o o o
0 «d-
in 1-4
000


to co r^.
to CM to
en co in
CO CO i-H
"— i CO P^
0 «q- CO
^H r-( CM

a
«<:
O ^ ' A
tn. « ^_
•r- CO p—
C CO >
O CO >,
co re tu
3 c x
o. s 5










a; , «.
0) . (U
3 as
re as
u i —
 §

o o o
C3 O 0
o o in
in i — co
«?!• m to
to to to


i— 1 CO »-4 CO
•* r-< CO

CO CO «?T CO
, — < 0 — • 0 0
CM CM O CM OV
CO «fl- i— 1 in r-H
i— 1 i— 1 i— i CO i—t
to rv, to •— t if
CO CM <* O «3-
o o o o o
i— 1 to CM ^~
in co •=!• co
•-H O i-i O O


o CM en r-i p^
r**» to to co r*^
tO CO CM CM O
i-H O O i— 1 i-H
«-< «3- in r~» o
f— < to «— i in co
^" to ^1" to
CM
•-§
A
en o
S- cC < 0 2=
CO •
, A A A \>*
CD re a) a> 2
c t. c T- re
i— so co re
!— ' 	 O O =
LU a. ca «c i—
£_
ai
^
c^


c
Q

C '
t. i- 3 ^£
^J ^J 1^ 't^
> > j_j OJ C-
.,_.,_ co £_ aj
OS OS -^- C_3 >
O C -iS -i-
co c o -^ -i:
j^ re c co o
o > «s3" to .in
co to en en
.rt •
O r- < «^ >— »
•— t co co en
O CM 1-1 O
CO to CO CO
«;r r»» CM in
O CM CM CM


in co CM CM
»«• co in o
I— 1 i— 1 f-H
to en o i— «
en in ,-1 -3-
I-H CM i— 1
^* en p>» co
- CM CM to en
r-l i— 1


2= 2: 0
t_ .C «
«a o o o
C CO 2= _
CD c. re
r— . oj «-a
•a -M o s
cz *j ja s




v


£_
a-

£- -r-
«- >
<1> -r- -0
t- > as .re
QJ T— O
> OS GJ «_
•f— F— CO
as c i—
•t— V- 4->
' Q.JK: > co
aj -a c t_
aj re -i— -P-
Q >- _I Ix.

0 O O 0

o o o o
O O O 0
in o in in
cr> •— i co CM
en i—i co in
O i— 4 r— * i— I


                                                 C-9

-------



°J
^s

e
•— . * m
•-* ">
*-H m
•

*















6





•





















•


T*
4J
C
O
0
^r
i
CJ

UJ
1
CO
«-f
O
rH CM
,o* o*

o
rH
0 |0*

cr
^




o *—*
r- CM
U. HH
ta vt
CO 4-
t- o




tr

CM
o-
rH
O
rH
O*


-


lo-

.,- _>
c_

•r-
CO
4J

I <
co
,'
O

CO
CD CO
co o

CD
j—
0 0 O
• • • C?
rH in CO
CO

in rH ^f
o o co o
^O ^3* ^M ^^
O CO O rH
rH co CM in
CO rH in
in o o
o o o o
rH VO «*• CO
vo CD in I-H
O O rH O O


vo o o r>. co
r^ CM CM I-H cn

co o in en co
rH oo «* oo cn
rH
co cn co CM o
r-. CM «a- PO I-H
CO CO VO rH i— 1

s'
,,-J ^J
« —1 U. U.
C 
s v> o •<-> jz
§t~* ta cn
in o ••-
8r5l]SS


£.
cu
> f
•r~
Q£

CO CO
CO CO
cn , t_
t- r— J«J ^£ O
CO 3 CO CO
> S CO CO £-
•r- O C_ C_ CO
oi o o cj j->
(O
3 CO  rtJ 3
O r- CO 3 -S£
r— •»-> ex: o

co -I— ta o r—
>i— I 3= -^ CD

^C ^C i i t
Co CD LL. LL. LL.


. O O 0, 0 O
0 O 0 0 O
in o CD I-H m
r-» vo co r». CM

o I-H co cn CD
CM CM CM CM CO

CM
O '

co CM P>. I-H r»»
rH in rH «* CM

. -..
* ^* *3j" LO CO GO
*T O rH CM f-*
^O CO CM ^ O
CO CO GO CAJ O^
. -*
oo rH m vo vo
m i— i CM in I-H
«*• O rH i-H rH
o r^ co m CM
o m o co rH
vo o CM vo co


en I-H co o vo
CM CO O CM CM
<-H rH CM rH
i*- m cn r^ o
co co co CM r*.

m vo «a- •53- vo
co «a- «=»• r«~ in
•in CM I-H «a- oo
_.s
_J U_ •>
— J  CO
« 3 Q.
A M cn co ta
4J r-* S~ *r- <•*
t- r- 3 > CJ
O CO t^> JJ
'D.4-> CO VI CO
3 I/I r- CO C
co 3 
CO CO -i-
> t_ - cc
•r- O
Q£ J£ t- CO
t_ CO CO CO
in CO CO > -M
.Si -M (_ >r- 4->

t. g o cn  co r— ta
• O) -r—  03 O J* CO
•r- > CO O
0 C- C_
.C 0. C. O >>
,0 ITS CO CO
§N^
t_
££^£,2

j^ 1 I C/} t
CD ^C ^ 21 <


OO O O 0
O O O O O
o o m o o
CM CM CM ««• VO

en rH CM co in
CO ^* ^" ^" ^"




oo
VO O
CM

CMrH
O OO
vo PO «a-
CM rH P>
co com
^
o oo
,0 00-
CM f^-
CO rH
o oo


CM O rH
m vo CM

cb coin
o cn CM

CM CM «T
in cnco


co
CO S
y
co •
•> s o

CO-r-

co u
CO >>
C. CO C
(_}•!— ITS
|SJ -^x •
CO C O.
CO  0

CO «r^
t. > Wl O
O 3 ^ O
0 Q.
t. •••- t_
•  CO
•t— CO
rv t
-^ C_3
O C
CO O t-
t_ C^f CO
O C C- >
'O ^2 'O
<- J= U. CO
(0 J* CQ
• en o >%
3 3 t_ •
CO CQ O. _1

^C ^* ^* ~~
0- 3 3 O


o o o o
0 0 0 O
o in o in
in co in cri

CM in vo o
O O O rH

CO
o
••••

-------
 OJ
  ta
. •!•—
  >
  CO

































J
c
0
LU
_J
CQ
^
1—
S CM
o- o-
0

o-


3'
o ^->
r— CM
U_ r-l

113 VI
, CO 
*os
CD ,
2
•S
CO
CD 0}
C0.cn
.CD

VO CO O CO
«*• «3- CM
CO
O O rH O O
I-H en rH ' in en
co, vo co en o
1-4,1-4 rH CO
CO
1—4 rH i-H CO O
O O rH O O
CO CO OO
ro. vo CM o
O O CM O O
• « « e



O C7» CO CM VO
LO -sj- en I-H co
rH
en LO o O ro,
en en CM LO rH
r-l CM rH

ro. cn o co i-H
CO LO O CVJ CO
LO «3- CO rH •
3: «=c
0 3» .
A A ^X
Z CO C • 3
O to O
•r- +J > «
« c. e 3 t_
Jtt 0. O. CO
C- t- 10 «JS
ro co ,:>> c. +j
^ *4J 10 cj C' .
^ScDC^cS





t- t- CD -1£
Ol .CD
> > O) -i- O!
•r- -r- >, Oi £_
as as -r- cj
a: to
en en = «_
c: c  «C
0 0 > 3 >
O O O O O
• o o o o .0
LO LO O LO LO
vo ro. o .vo co
^i- LO ro, co i— i
I— 1 rH rH rH CM
co
o
1—4 CO VO LO
CM «* r-l rH
• -
O O O O O
•4* VO O 00 CO
O CM CM CM rH
CO rH CO r-l CM
CO CO CO
O CM O CO O
O O O O O

LO I— 1 ** rH
o o o o o



CO VO LO CO rH
CM «d- CM «* «3-
en «^ ^f co co
co fs. co ro. en


en en co co vo
en CM vo en CM
rH CM CM CO
s: -SS.
rH' A *— i
— ®
S O C r3 CO
O O •!- 1 	
" c en to •'m
>» 3 c t. >
O) O -i— CO -to
i — -M -U r— i—
J— iwi C r-— CD
•gssss




^<:
co e~ CD
CO CO t.
c. > c. cj
CJ -r- CJ
a; > ji +j
CD -r- CO 3
C -r- OS CO C '
OS S_ r- .
•M IO CO CJ .-O
tO »r™ r— 35
+•> S 4-> r-r
CO JJ r— C31
.C • 'r- (O 'r-
3 — I — I U_ CQ
3C 3C 3= Z Z
oo-S55
o o o o o
o o o o o
LO O O LO LO
^S- O «3- CM fo.
CM ^f CM LO LO
.CM CM CO CO CO


OO CO
•a- co
- f
LO rH
.00
rH CO
CTl rH
rH CO
CO
CM O
0 0
'•
,cn I-H
O'O



CO CO
CO «5T
CO CO
ro. LO
rH 1—4

O CM
VO CM
1—

CD "ZZ.
r- 1—
•r- ^
>. to
c to
C IO
•^ °
S in
£3


£_
CD
>
£
to
t— CO
CD c
> 0
as oo
to -^
•r- C3
r^ ^
O -•
CJ LU

1— 1—
o o
o o
O LO
CM CM


VO *^ CO ^* OO
CO CM r-4 i-H CM

LO co vo CM en
«3- CM rH 0 0
CTI'CO LO rH •*
t^ ^^ OO »^^ OO
CM rH
LO
co en ro. o CM
rH O O O O

	
.•^ CM rH O O
-


LO LO rH CO «*
OO fo. VO CM CO
«=}• co en I-H vo
O en ro LO to
rH

en co co I-H rH
(o, CM CM CM CO
CM LO rH CM

^^ t*^
CX3
»»-4 <_t
3—i2 «
S a- en>— i
•O ' ' • « 3 .
C JJ OS J2 '•>
(O CD t- CD'
r— C CD 03 r-
f— t» c "o en *
tO fO -i— OJ (O
 > CO CD CD
•1- .1- > t_
a; as -r- cj en
OS C
CD ^^ £- -r-
.-Jf ^ c -S "o
-C i— -i— CD O
3 co a. o _i

3 S 3 3 S
0 O. O O 0
o o o o o
LO O LO LO LO
ro. vo «* vo «^-
CM «!j- vo oo i— i .
O O O O rH
^.
o .
l/> «5 <=>  1—

to its O O *-
CD i — en e_ to

,ifc iE £ S S




. •
CO .
OS C. SU Jbi
t- cj en CD
CD CO > C CD '
r- T- as ess cj
re os '
00 ^4 t- S-
• CD U O) (O
en c  -a
•i— •!— r— -c— CO
CQ a. co as cj

SSSSH-
o o o o o
O O O O CD
O LO LO LO O
co LO en vo o
CM LO LO VO CO
rH rH rH rH rH


                                                                                                     r.  11

-------
3
IO
at


























*
























T*
c
a
CJ

LU
«J
CD
*^C
1—








I.
3
0 "^
i— CM
1 f hurt
U"* rT;
E "«*.

•1—.
0£
CO
+J
1 *
s
c/3
CD CO
co cn
3D (13
CO

r^r-
CM CM


fH LO
O O
cn vo
P-. VO
CM i-H
co en
0 ^H
0 0
#
CO
O LO
o o



•*f CO
CM 10

CM CM
P*1* r~i
1— <

f«H r^
P-. i— i
CO ^T
i-H
A
CO
OJ
[ "CJ
5 *j"
CO
A
C t/)

CO T-
> ^"^
C£ CO
CO
c 
O CO
c- c
5-
03^
*•
0 0
0 0
o to

cn CM
•—i CM

^"
o
o cn cn co
CM CO CM

• " -"
•3- CO LO CO
O «* CM O «3-
LO en en en en
to o I-H
CO LO .-H O
cn co T-H co
o o o o co

VO CO
cn co co vo
O «—< r-H O CM


10
CM O CM «=J- O
O «3- 'S- CO VO

^" LO cn *^* *™H
i-> *y to vo vo


co LO cn •— < vo
cn P^ vo vo »!r
CO CM CM

3
*™^ _^T
' *s  4-> SO
CJ 3 t- C rtJ
S5 CO
> T- CO
•r- o; t-
Q£ O
-U
^ o t. ^:
CO CO ^ 4->
(B OC > CO C_
CO -r- CO  J^
i— U_ c «- .O
r-* - CO CD CO CQ

Z |E I-H eC w

O O O O O
O O O O O
O IO LO f- LO
., co to co r*-*»_ vo

en co •— i •— i o
CM CO ^3- -sf «3-

LO
o
cn «*• »-H CM
CO CO LO CM
"

co r^ CM I-H
1— t 1— < l-H O O
(— * P^ o cn ^°
i-l O IO CM LO
i— 1 i-H i-H CO LO
O O O •— < VO
CO CO ^H O O
o o o o o

r>. co
i— I O LO O '
r-4 r-l O O O



*sj* co co vo en
«3- CO CM I-H O

vo vo i— i r-. cn •
VO LO «3- LO *3-

-
co vo LO co en
P** ^* CM P^* ^"
CM i-H «3- LO CO

I-H
3
A .^J ^^
C ^""^ ^T ^C ^"^
0 S i-«
"cn co" « «jT'
e to c ns i—
•i- o •<- c: -> o 3
t_ c to i— t_
cSc^^




c^
CO

•r"" "^^ £••
o; co co
CJ U > «_
 0£ >
c e £ _c S
O t- Ifl
4J O t- ••- J=
n3 -C -• -J ^ <: «c

O O O O 0
o o o o o
LO O O LO O
CM *3r r~» LO — i _1
A r\ I-H I— 1
, to «~l * 'O **
J3 CO IO i —
•r- J= CO •— . 	
C J-J C -u O)
C t- t_ c 3
= s<5<£li



c_
CJ C.
> CO
C-r- >
CO OC T-
> A^
•I- WJ
*^ r** CO CT ^^
CO C O CO
CO CO •!- -i- CO
£_ CO «3 i— t-
O J^ r- ^- O
03 Q. .i—
t. ^tf g ^J
» LO
CM o <— « «a- «a-
O i-H O -^ «3-

vo o «a-
«S- CO CM «*• CM,
•— t CO i-H LO CO



LO •— 1 CM P«- CM
 O -C O - CO > CO
v- Ol •-- t_
oc t_ o;  >, Or-
o -i- c- j= -u
•»~ o; nj o 4-J
j^* ^C? ^ * ^3
O C C T3 4->
« re 3 c co
•— 3 O CO 4->
sa co ea 3 t/5

1— 1— Q < C£.

O O 0 O O
o o o o o
CD CD LO O LO
LO CD «— « LO P»»

co r~> CM LO p^.
CO CO CO ^f •— «

CM
•— 1
                                                                                               -1
                                                 r-1?  .

-------
















'



























T
•*->
w
o
c_>
u
LU
—1
<=Q
2




.



• 3'
' O •*"-*
«— CM
1 I !_(

<5 co
c. u
_i_\ ^~- .-•
oo



c
Q













o
•p—
^rt
O
o
— 1

CT
(O
C3











-
<

01
i — 11 Osl

O
P-I
cr 1 cr
cr


CM
o-
v— t
o
i-H
r*.

20-

Icy
iv"'
<0 C\P

» t— i o «3-
' P*» en ^j* en en
^ .«3- CM «S* •— I
O o r^ CM in
en o CM en «s-
vo o «a- in CM
f—4
in o vo co vo
in CM en co in
CO I-H f-4





• «c
3 " cr
c fa O) IB a)
»-H CJ — 1 2: OO

 •
/C*
£_
-C 
•f~ *f™ Cu &_
£ OS QJ OJ
O > >
t. 4J T- T-
>,r- {_ OS OS
-^ O QJ '
c/> .f— > >> ^
•r^ f^ trt
^£ JX Q^ r»« rtj
t. C « C
O O C 3 -H ^f


CM 00
O CO
CM CM
o r^^
in «*
O CM
co «a-
i— i







3 3
M A
1 :
p*™" ^^
r— ' (XJ
•SOS

£_
di
>
OS'
- s
*"^
C- 3
 .Z
•p- O)
OS £=

fO -^
Q. t.
is oo

< «C
3 3'

O O
o o •
in o
CO «5l-
i-H CM
O C3

i en co co r^ in
. co in «* o I-H

CO OO CO CO i-H
CO PH i-H r-4 O
CO 00 O O CO
oo CM co en o
i-H p-4 i— t CO
in co «*> en I-H
o «3- CM I-H o
CM O O O O
o in o o o
CO CM I-H CM O


r^^ r^ en in r**
co, co *^* r** o
i— i
^- i— 1 O «— I «3-
^f ;^r co vo CM
CM l-H pH
VO O CM CO O
CM CO CM T-H O
co co vo CM in
CM


a
I-H Q
A ^
Q co "
i— i cj> aj
C OS C  a. * •>
<- OO 19 3 t_
t- -COO)
•p- C «O i — Q.
3 -P- C •— O
.<^r?^'^5







^, c
t- d) OJ
S. U 
£_•!-
> > -i- O OS
•r- -r- OS
OS OS C . in CM co I-H
«a- cc en I-H in
O i-l CM CO CO
CO
1—4
»3- VO CO i-H CM
CM p-4 CO CO CM

in ^ vo in ^^
O VO r-4 O i-H
in 1-4 en CM r^
in co o r*» en
rH i— 4 i— 4
•'•*
p-« co »sr en co
O CO i-H O CO
CO CM
co co in o co
O i—4 ^ CO 00
CM

in r^ r^, oo CM
co CM en o -H
co I-H co ^*
in I-H o CM in
vo ^ en I-H r*^
,co CM vo in
I-H m CM in P*.
vo *3- en CM oo
CO CO i—4 CM

OS
o

,- «
f^
0 • eC
• « oo c o
co a; o o
 (/T NX
:, — c. o -i- o ,
O (O <** QJ (X3
CO — 1 (O = O-
• s~
 p— 3 >
•1- C_ -r- O) •!-
OS « P*. m CM vo
i—4 in «3"CM CM
O O i-H CM CM
1—1 • '
' in
VO CM CM i-H «tf- O
in CM i— 4 vo "ff vo

p-4 CO i-4 CM P*. CO
o o o o o o
I-H I-H oo vo o CM'
CO pi CM CM .-4 CM

O OO i-H O •— 4 CM
000000
VO VO O CM
co oo 1— i co in «a-'
O I-H O O O I-H


vo oo ob vo en CM
co I-H o in co vo
CO I-H O f«
P-^ in oo ^* co co
en co I-H p^. p^ co
CM VO CO i-H O «*•
en CM in o co CM
in P-I.CM P«. co «^-
p-4 i— 1 r-H VO i-H
1-H


i

t_ •
o os es •
0 0 OS
. "OS OS O
JZ O O f » .

^3 * * fl3 * ^3
S • >> * • £
O • fl.» ^ O —
i— . • r«. {J £J >PH.

Jd QJ *f— t_ fO ^O
0- 2 Q CQ _J H-
• .
cu o t_
> 
•i- 
OS • +J -r-
_ 4J OS
13 eo co
•i— C- -r—
4-> U Ol 0) O
. C OJ > r— C
t— IB > *p- 4-> *r—
'03 OO .f— QS J-> i—
> * OS -I— r—
•p- • fl3 _l i— i
OS SI C 3

CO OJ 4-> SX £_ £_
- •— « j= o o
>>•!-> i — = U_ U-
£- +J (O
f3 «p- 3 * * •
S -J p— OO OO LU

es as es os os ' es'-
f-^ ^^ <"» CO ^"^ fj>

CO O O O C3 O
O 0 O O 0 O
o in in -o in in
i-H CM CO CM i-4 CM
r**» co o r-H 
-------

-------
                                          VII  (0)
                                          Revision  No.  0
                APPENDIX D     '

         Computer Program for the
Probabilistic Dilution Model  - Point Source
                . (PDM-PS) .   "-

-------

-------
                                                            VII(D)
                                                            Revision No. 0
         This appendix describes a computer program (PDM-PS) which performs
  the computations of the Probabilistic Dilution Model for Point .Source
  discharges using numerical methods based on quadratures.  The program is
  written in BASIC for the HP-85 and the IBM-PC, and should be readily applicable
  to other personal  computers with perhaps minor modifications to reflect
  individual  machine characteristics.

         The  program is  structured around  a slightly different input  format  -
  than  that used  for the manual  calculation using the  moments  approximation.
  A  series of normalizations  (ratios) of certain of  the  input  data  items  is
  used  to provide  a  computation  framework  that provides  a more  generalized
  perspective.                                     ,                         .

        The  appendix is  organized as follows.   Section  1 describes the
  basis for the formulation and  normalization of the input data, as .used in
  the program.  Section 2 provides an annotated  description of the CRT and
  printer functions, as well as the nature-of the user's response.  Figures
  D-l and D-2 provide the results of running the  PDM-PS through the example
 described in Section 3.2 of this report.   Finally,  Figure D-3 provides a
 Tisting of the PDM-PS program for entry into a personal computer.

 Drl•  Formulation and  Normalization                        •           .

        The analysis can be made more useful  in  a general way  if the
 normalization described  below is applied  to  reduce  certain  of the.inputs
;to  readily recognized ratios,  and to express results  (stream  concentra-'
 tions) as a  multiple or  fraction of  the target  stream concentration  (CL).    .
                                     D-l

-------
                                                          VII (D)
                                                          Revision No.  0
        The explicit assumptions in the normalization scheme  that  is  used
               . .         „                         , .                 .
 are that:

        o  The stream target concentration (CL)  is  produced when the
           discharge flow is the mean effluent  flow (TJE"), the discharge
           pollutant concentration is equal  to  the  permit effluent limit
           (EL), and the stream flow is equal to the  design value  (here
           designated 7Q10 - though  any other basis may be used for desig-
           nating the numerical  value of stream design flow, e.g., 30Q5, '
           30Q10, etc.).

        o  The'reduction factor  (R = "CF/EL) determines the mean effluent   •
          concentration of  the  pollutant  being evaluated.  It could be
          selected  arbitrarily;  however,  as applied in this manual for
          evaluating the permit  averaging period, the value selected  will
          be  dictated by the variability of effluent concentrations and
          the permit averaging period.

 In the  usual  case, where the stream target concentration  (CL)  is  set  at
 the chronic toxicity level, the multiples .of the, target  - in  which stream
 concentrations are expressed (CO/CL) - correspond with the  acute toxicity
 level.  The basis for the normalization scheme  adopted is as  follows.

       The downstream concentration, CO, is  given by  the  dilution  equation:
                                                         "         -     *'
                             co = QTTir*CE.                     ^
                                            •         ''      "   '           .
For a chronic criteria concentration, CL,  the effluent limit concentration,
                                     D-2

-------
                                                          vii  (d)
                                                          Revision No. 0
 EL, is computed using QS = 7Q10 and  an  average effluent flow, TJF:


                  '-    CL =             EWSTO                       (D-2)
 where STO         .                  (p.3)



 However,  the  choice of permit averaging period forces a reduction of


"CF of magnitude, R, so that permit violations occur only 5 percent


 or 1 percent  of the time. Thus the actual long term average, effluent


 concentration  is:



                         ?F = R EL = R CL/0STD       .             .     (D-4)



        The problem is to compute the probability  that the  downstream
                        \       ' .                        ,             .

concentration exceeds a multiple,  3, of the  chronic  concentration, CL.   In


particular, if the acute criteria  concentration is  selected", then  p  is the


acute to chronic criteria ratio for the pollutant being  regulated.   Hence

It  is necessary to compute:



                 Pr [CO >   PCL] = Pr [CO >  p
-------
                                                          VII (D)
                                                          Revision  No. 0
 and CE/CITis the normalized effluent concentration.   The  probability
 distribution of this  random variable no longer depends  upon the mean
 effluent concentration,  but only on the coefficient  of  variation, VCE.
 This is easily seen from the following  representation of  a log-normal
 random variable:

                         InCE =  InCE + Zo-lnCE                            (D-7)
t       *^^               '                  •          '
 where CE is  the median,  
-------
                                                             VII (D)
                                                             Revision  No.  0
    Note that QS/QE is  log-normally distributed  since  both QS and QE are
    assumed to be log-normal.  Thus, only the ratio of the average flows,
   WOE, is required.  A convenient normalization using ratios that are
   more, readily available results if the average effluent and stream flows
   are standardized relative to design stream flow (here designated by 7Q1Q).
   Defining

                ..         ,        Fl  = 7Q10/QS     '                    (D-12)
                                  F2  = 7Q10/QF           ". _.             {D.13)
   Then
                 '          ;    \Wqr-  F2/F1                  ;         (
   and
  These ratios, Fl and F2. together with the coefficients of variation,
  VQS.VQE, and VCE, completely specify the characteristics of the random
  variables in the normalized dilution Equation D-ll.  R specifies the
  effe,ct of permit averaging period and p, the acute to chronic criteria   ,
  ra-tio, specifies the toxicity behavior of the substance being considered.
  This completes the normalization.

  D~2-  Description  of Program Use

        The program  is easy  to  use.   The  values of  the  input variables  are
  sequentially  requested on the  CRT.   Once  the  input  values are entered, a
  summary of the input data is printed out, as  is a tabular listing of the
,           : .                          D-s

-------
                                                         VII  (D)
                                                        . Revision No. 0
results of the calculations.  The user should be thoroughly familiar with
the theoretical and practical bases for the PDM-PS as described in Chap- -.
tens 2 and 3 before attempting to use the PDM-PS.
                                                         *  • •
       USER:     Initiates program execution.
       PRINTER:  Writes'title.
       CRT:      Displays title and general  descriptive material  shown in
                 Figure D-l.           •     ;

       CRT:      Question #1 is displayed:  "Enter coefficient of varia-
                 tion of  QS, QE, and CE:"
       USER:      Enters the values of VQS, VQE and VCE,  separated by
                 commas.                *-                       .        .

       CRT:      question #2 is  displayed:  "7Q10/avg  QS?"
       USER:      Enters the ratio  of the  7Q10  flow to  the average stream
                 flow (TJS).

       CRT:       Question #3 is  displayed:  "7Q10/avg  QE?"
       USER:      Enters the design dilution  ratio,  i.e., the ratio of  7Q10
                 flow rate  to the  average  effluent  flow rate (T|E").

       CRT:       Question  #4 is  displayed:   "avg CE/EL?"
       USER:      Enters the ratio of the average effluent concentration
                which  the  treatment plant will be designed to produce
                 (avg  CE),  to  the effluent concentration derived from the
                                    D-6

-------
                                                  VII (0)
                                                  Revision No. 0
          WLA analysis (EL).  This latter value is that concentra-
          tion in the effluent which will result in the stream
          target concentration being met, when the following flow
          conditions prevail:                  ,                 ,.

             Stream flow (QS) is at the 7Q10 flow rate.

             Effluent flow (QE) is at the average discharge rate of
             f 1 ow.

PRINTER:  Prints a tabular summary of the input data selected.

CRT:      Question #5 is displayed:  "Enter lowest, highest and
          increment of multiple of target for which % exceedence -is
          desired."
USER:     Decides on a range of stream concentrations (expressed as
          multiples of the target concentration1,  CL)  for  which the
          probability of occurrence and the recurrence interval  are
          desired.   The user enters,(1) the lowest  value,  (2)  the
          highest value and (3)  the incremental step  desired  for -
          values  between  the highest  and lowest.
PRINTER:  Prints  tabular  listing of results.  For each multiple  of
          CL,  the exceedence frequency and  return period are
          listed.  When the'printing  is completed,  a  tone  sounds
          and  Question  5  is repeated.

USER:     Enters  a  new  set  of  values  for multiples  of  CL,  if
                              D-7

-------
                                        VII (D)
                                        Revision No. 0
desired.  This allows the user to conveniently search out
              • •.,.       .
the ranges of interest and select the most appropriate

levels of incremental.detail.  When the desired amount of

output has been  obtained9 the program is interrupted, and

begun again at Question #1 to examine another set  of

conditions.  The user can formally "end" the  program by

entering 0,0,0 in response to Question 5»
                    D-8

-------
 POINT SOURCE - RECEIVING WATER
   CONCENTRATION ANALYSIS
 INPUTS:  COEF VAR OF QS,QE,C£
          RATIO... 7Q10/avgQS
          RATIO... 7Q10/avgQE
          RATIp...ayg CE/EL

   BACKGROUND STREAM CONC (CS)
     IS ASSUMED TO BE ZERO
                                                           VII  (D)
                                                           Revision No. 0
 GENERAL  DESCRIPTIVE  MATERIAL
 ENTER  COEF VAR  OF  QS,QE,C£?
QUESTION
ENTER  FOLLOWING  RATIOS:
        7QIO/avg QS ?
          '  '
 .05
        .7Q10/avg QE ?
        .avg CE/ EL? -
•.67
ENTER LOWEST, HIGHEST, AND INCREM
ENT OF MULT OF TARGET FOR WHICH
 % EXCEED IS DESIRED
?
ENTER LOWEST,HIGHEST,AND INCREM
ENT OF MULT OF TARGET FOR WHICH
 % EXCEED IS DESIRED
7
2.5, 3, .05  .        ,  ••..'..•
QUESTION #2

QUESTION #3

QUESTION #4  ,

QUESTION #5 (CONTINUES TO REPEAT
             AS NEEDED)
                        Figure  D-l  -  CRT displays,
                                    'D-9

-------
                                                         VII  (D)
                                                         Revision
               No. 0
•***********************.**********
    RECEIVING WATER CONC (CO)
     PROBABILITY DISTRIBUTION
        AND RETURN PERIOD
 - FOR MULTIPLES OF TARGET CONC
   DUE TO POINT SOURCE LOADS
*********************************
COEF VAR..»..QS =
COEF VAR.....QE =
COEF VAR.....CE =•
7Q10/avg QS =
7Q10/avg QE -
avg CE/ EL *
VIOLATION PERCENT
MULT OF OF TIME
TARGET EXCEEDED
1.00 0.894
2.00 0.112
3.00 0.024
• 4.00 0.007 •
5.00 0.002
2.50 0.050
2.55 0.046
2.60 0.043
2.65 0.040
2.70 0.037
2.75 0.034
2.80 0.032
2.85 0.030
• 2.90 0.028
2.95 0.026
3.00 0.024
1.50
0.20
0.70
0.05
3.00
0.67
RETURN
PERIOD
( YEARS)
0.3
2.4
11.3
39.4
114.4
5.5
5.9
6.4
6.9
7.4
8.0
8.6
9.2
9.9
10.6
11.3
TITLE
                                                   SUMMARY OF INPUT DATA
                                                   CALCULATED RESULTS .
                  Figure D-2 - Example of printed output
                                     D-10

-------
  t
Start
     Clear screen
     Print header
       messages
    Prompt for and
    input coef. of
    variations for
      QS, QE.-CE
    Prompt for and
     input ratios
      of 7Q10/QS,

      7Q10/QE,
    and  avg.  CE/CL
   Compute normal
 and reverse normal
    coefficients
Prompt for and'input
lowest, highest,  and
 delta increment  of
 multiples of CO/CL
       to use
                                  .  Print input values
                                     and table header
                                           1
                                         Iterate on CO/CL values
                                      Evaluate Q(x)
                                     Compute return
                                     .  . period
                                 Print CO/CL,  % of time.
                                      exceeded, and
                                   •,  return period
                                       Next  CO/CL
         Figure D-3 - Flow chart for PDM-PS program.

-------
                                                     VII  (D)
                                                     Revision No.  0
10
20
30
40
50
60
70
80
90
100
110
128
13©
140
150
170
180
190
200
210
220
233
240
259
260
270
280
290
380
1
i
!
i
i
i
i
i
i
i
i
t
t
I
!
t
i
i
i.
i"
t
i
i
;*
j
i
i
i
D
                POM-PS
            PROBABALI3TIC
           DILUTION  MODEL
       FOR POINT SOURCE DISCHARGE
     DEFINITION - INPUT TERM*

    QS  «  STREAM FLOW
    GE  »  EFFLUENT FLOW
    CE  -  EFFLUENT CONCENTR.

  7Qi0/av-gi3S =       RATIO
     SPECIFIED STREAM FLOW'i

  7Q10/av-3i?E -      DESIGN
    EFFLUENT DILUTION RATIO

  av^CE/EL    =    RATIO OF
    THE SPECIFIED AVERAGE
    PLANT EFFLUENT CONCENTR.
    -TO THE  EFFLUENT LIMIT
      CONCENTRATION.
             EL IS THE EFFL
    CONC  THAT PRODUCES THE
    STREAM TARGET CONC WHEN-
              AND QE=av-9i3E
DIM R5<32),Z5t32>
DIM R<8>,S8<8>
DIM PC
PRINT   .

PRINT  "    RECEIVING  WATER
NC          PROBABILIT
ISTRIBUTION  "
PRINT  "   •    AND  RETURN
IGD"
PRINT  " FOR-,MULTIPLES  OF
GET CONC.    DUE TO  POINT
URGE LOADS"
PRINT
320
330

340
350

360


370
380 D.ISP  "POINT  SOURCE - RECIEVI
    NG WATER"
390 DISP  "  CONCENTRATION ANALtS
    IS"
430 DISP
410 DISP  "•
                               CO
                              Y D

                              PER-

                              TAR-
                               SO
420 DISP  "INPUTS--  COEF  VAR OF QS

                   RATIO...7(

                   RATIO...7C

450 DISP  "         RATIO... av<3 CE
•*30 DISP
    v^QS"
440 DISP
                                     470
480

490.
560

510
520

530
549
550
560
570
580
530
680
610

62>3

630

640
658

660

670

680
690
                                     DISP
                                     DISP  "
                                     C  
                                     BE ZERO
                                     DISP
                                                 BACKROUHD STREAM
                                                       IS ASSUMED TO
                                         DISP
                                         DISP
                                               ENTER .COEF VAR OF OS, ft
                                         INPUT  V1,V2.'..'3
                                         DISP  "ENTER FOLLOWING RATIOS
                                                                QS ".

                                                                i3E -,

                                                              EL " ?



                                                               COEF ;

                                                               COEF .

                                                              .COEF 5
                                                                   Cli
                                                           • f *• + -r 4. 4.
     DISP  "  ..
     INPUT Fl
     DISP  "
     INPUT F2
     DISP  "  . i .  . . av;-^
    ERCENT" ; TAB <25 > .= " RETURN "
    PRINT " TARGET  "; TAB< 13) .= "OF
     TIME"; TAB (. 25 > ; " PER IOD "
    PRINT "  "jTAB<13>;"EX
    CEEDED";TAB(25>;" >
                                        W3=SQR
-------
                                                 VII (0)
                                                 Revision No. 0
 850
 860
 870
 890
 920

 930
 940
 950
   ' .
 960
 970
 980
 990
 1910
 1 020
 1030
 1040
 1 850
 1070
 1080
 1090
 1100
 1119
 1120
.1130
 1140
 1150

 1160
 1170
 1180
 1190
 1200
 1210
 1220
 1230
1240
1250
1260
1270
1230
,1290
1300
1310
1320
1330
1340
1350
1360
     INPUT. B1,B2,B3
     IF B1+B2+B3-9 THEN 1199
     !  - LOAD QUAD. WGTS.  & ROOTS
     GOSU8 1488
     !  - COMPUTE PORTION OF' GKH)
     ARGUMENT IHDEP OF C©
     DIM Z9<32>
     FOR 1=1  TO NO
     !  - EVALUATE USING INV PROB
     TRANSFORMATION
     P9=R5
     GOSUB 1380
     Z9 < I > =LOG < 1 + EXP < U9-W9#X9 > > -U
     3    '      • . .   " ' ;  '  .
     NEXT  I
     !  - CONCENTRATION LOOP
     FOR C0=B1 TO BZ STEP B3
     15=0
      !  — QUAD LOOP—  EVALUATE. CK
     X> = F  AND SUM
     FOR  1=1 TO N0
     X= CLOG < C0 > +Z9 < I » /M3
     X0=SGN
     X=ABS-;X>
     F= 1 +X* < 0 1 +X-* -, D2+X* < 03+X* < 04
     +X* (. 05+X*06 ;• > > > ;
     F=.-'5*FA-lb      .     •  •   .
     IF X0>0 THEN 1090
     F=1-F  '          .     '
     I5=I5+F*Z5
     NEXT I-
     !  — COMPUTE RETURN PERIOD
15=100*15
PRINT USING 1150'.: C8.I5.I0'-
IMAGE 2DZ.DD,5X,2DZ.3D, 5X,3
DZ.3D
NEXT C8
PRINT 8 BEEP
GOTO 840    -
FOR L=l TO 7
PRINT
NEXT L               '
END
!  -SUBROUTINE TO LOAD NORMA
L AND REVERSE NORMAL COEFFI
CIENTS
01=. 0,49867347
D2=;.0211410061 ,
03=.0032776263
D4=.0000380036
D5=.0000483986
D6=.000005383
!  -n-t-++++++-
El=2.515517
E2=.802853
E3=.010328
E4=l.432788
E5=.189269
F6=.001388
1370 RETURN
1380 !  -SUBROUTINE  TO  COMPUTE  IN
     VERSE NORMAL TRANSF
1390 !  POLYNOMIAL APPP.OX  TO  INVE
     RSE NORMAL  TABLE
1400 DEF FNC*S9
 RETURN
 ! -QUADRATURE SUBROUTINE -
'COMPUTE  ROOTS AND WEIGHTS
 !   15 =  INTEGRAL
 !   R5  =  N0 ROOTS  C+- GR
 USSIAN ROOTS «. N9/-2 LAGRR R
 OOTS>   '          .-•'.'
 !   Z5CN0.)  =  H0 WEIGHTS
 !
 ! LOAD ROOTS AND  WEIGHTS C0
 R 32nd ORDER QUADS
 ! FIRST THE  GAUSSIAN  AND TK
 EN  THE LAGUERRE TERMS
 ! -QUAD ROOTS & WEIGHTS
 16th ORDER  GAUSSIAN
 Rl=8
       989400935  .
       944575023
       86563 12024
       7554044084
     R < 1 > =- .
     R<2>=-.
     R<3>=-.
     R<4>=-.
     R<5>=-.
     R<6^=-.
                                                6178762444
                                                4530 1677.76
RC8>=-.
88. =.
S8<2>=.
S8<3>=.
S8<4>=.
S8<5>=. '
S8<6>=.
S8<7>=.
S8<8>=.
N0=4#P1
!  CONVERT
WEIGHTS FOR
NTVL
                                                09501250984
                                                02 7. 1-5245942
                                                06225352394
                                                09515851168
                                                1246289713
                                                1495959888
                                                1691565194
                                                1826034154
                                                1894506.105
                 (0,1''
                                                            ROOTS
                                                           INTEGP
                                                        .
                                    1750 !  AND DIVIDE BY THO FOR COM
                                         POSITE FORMULA
                                    1760 FOR K2=l TO Rl
                                    1770 R5=.5+.5*RCK2>
                                    1780 R5=.5-.5*R
                                    1790 Z5=S8X4
                                    I860 Z5=Z5'
                                    1810 NEXT K2             ,
                                    1820 !  -LOAD THE LAGUERRE PfOOTS
                                       - AND WEIGHTS, PROPERLY CON'VE
                                         RTED
                    Figure D-4 (cont'd.)
                            0-13

-------
                                                    VII (D)
                                                    Revision No,
1S30  ! -loth  ORDER  LfiGUERRE  ROOT
      S & WEIGHTS
184Q  P<1>=51.7011603395
1850  P<2>=41.9494526477
1860  P<3>-34.5333987023
1870  P<4>-28.5787297429
1880  P<5>=23.515905694
1898  P<6>=19.1301568563
19013  P<7>=15.4415273683
1919  P<8>=12.2142233689
1920  P<9>=9.43831433639
1938  P<10>^7.07033353505
1940  P-5.07801861455
1950  P<12>s3.43708663389
I960  P<13>=2.1292836451
1970  P<14>=1.14105777433
1980  P<15>s.462696328915
1990  P<16)=8.76494l04739E-2
2000  Q-4.16146237E-22
£810  Q<2>=5 0504737E-'13
2020  Q<3>-6.297967003E-15
2030  Q<4>=2.127079033E-12
2040  Q<5>=2.862350243E-18
2Q50  Q<6>sl.8S1024841E-8
2060  €K7>=6.828319331E-7
2076  Q<8>-1.4S4458687E-5
2030  Q<9^2.042719153Er4
2090  Q<10>-l.S4907094353E-3
2106  6<11>=1.12999000803E-2
2110  Q<12)=4.732S9286941E-2
2120  Q<13>=.136296934296
£130  Q<14>-.265795777644
2140  Q<15>=.331057854951
2150  Q<16>=.206151714953
2160  FOR K2-1 TO N0/2
2170
2180
2130  NEXT K2
2200  RETURN
                        Figure D-4 (cont'd.)
                                 n-u

-------
                                                                     VII (0)    ,
                                                                     Revfsfon Wo. 0
  A>TYPE B:DILMOB.BAS -
  10 REM  -H.+-M.-M.   PDM-PS
  20 REK  / •   PROEABALISTIC
  30 EEt--       DILUTION   MODEL
  40 REM FOB POINT SOURCE DISCHARGE
  50 EEH.   '  '     '      . '  •   .
 '55' EEK-   AUGUST,  1984
  60 REJa  ,  IBK-PC AMD MS-DOS  COMPATIBLE VERSION
  70 REM    HORIZON SYSTEl-S  CORPORATION
  80 EEK       ,  (703) 471-0460
  85 REM                           .
  90 REIi •+*+-M
 300 DIM R5#(32),Z5#(32)
 310 DIM R#(8),S8>(8)
 320 DIM P#(16),Q#(16),Z9#(32)
 321 CLS
 322 KEY OFF            ;.'.•
 330 PHINT n**
 340 PRINT » BECEIVING WATER CQNC (CO) PROBABILITY DISTRIBUTION  «
 350 PRINT »               AND RETURN PERIOD"
 360 PRINT »   •       FOR MULTIPLES OF TARGET CONG"
 370 PRINT «           DUE TO POINT SOURCE LOADS"
 380 PRINT f******«*****«**5*»***»s.*«»**«K«
 390 PRINT "POINT SOURCE - RECEIVING WATER"
 ^00 PRINT "CONCENTRATION ANALYSIS"
 •410 PRINT               '             .
 •420 FRUIT "H •! •! i i .
 •430 PRINT "INPUT COEF OF VAR OF QS,QE,CE«
 •440 PRINT "       RATI0...7Q10/AVGQS"
 450 PRINT n       RATIO... 7Q10/AVGQS"
 •460 PRINT "•       RATIO... AVG CE/CL"
 470 PRINT " BACKGROUND STREAM CONC (CS) IS ASSOKED' TO BE ZERO"
 480 PRINT
 490 PRINT
 500 PRINT "ENTER COEF OF VAR OF QS.QE.CE"
 510 'INPUT V1,V2,V3
 520 PRINT "EHTEB THE FOLLOWING RATIOS:"
 530 INPUT " ..... ..7Q10/AVG QS ";F1
 540 INPUT "... ---- 7Q10/AVG QE tt;F2
 550 INPUT "... ---- ..AVG CE/EL ":F3
 560 PRINT
 565 CLS
 570 PRINT "  COEF OF  VAR.....QS =  ";V1
580 PRINT »  COEF OF  VAR ..... QE, =  «;V2
 581  PRINT "  COEF OF  VAR, ____ CE =  "-V3
 590 PRINT             '
 600 PRINT «        7Q10/AVG  QS = ";F1
610 PRINT "        7Q10/AVG  QE =  "jF2
620 PRINT "      .    AVG CE/EL =  »:F3
630 PRINT ' '.          '              .    •
640  PBINT "-S-H
         Figure D-5 - PDM-PS program listing - IBM-PC and MS-DOS compatible.

-------
                                                                   VII (D)
                                                                   Revision No. 0
720 i
730 V2cSQR(LOG(1+V2"2))
740 H3-SQR(LOG(1-f-V3~2})
750 W9=SQR(in"2+W2~2)
760 U9sLOGCF2/F1)+LOG(SQR(1+V2''2)/SGR(1+Vl''2))
770 L3=LOC-(F3*(1+F2)/SQR(1-t-V3~2))
7SO GOSDB 1160
790 PRINT "ENTER LOWEST,  HIGHEST, AND IHCEEKEKT OF KULT OF TARGET FOR"
795 INPUT "WHICH $ EXCEED IS DESIRED";B1,E2,B3
796 IF B1+B2+B3=0 THEN GOTO 1120
797 CLS
BOB PRINT w COEF OF VAR	QS = ";VT
804 PRINT « COEF OF VAR	QE = ";V2
805 PRINT * COEF'OF'VAH	CE = n;V3                          .
806 PRINT
807 PRINT «       '7Q10/AVG QS = ";F1
808 PRINT «        7010/AVG QE = ";F2
809 PRINT •          AVG  CE/EL = ";F3
810 PRIIJT                      ,
811 PRINT "•! I •! I -I 'I •! '1
812 PRINT                                     ,
813 PRINT «     STREAK COMC (CO)"
EU PRINT
.815 -PRINT n MULT OF";TAB( 13) ;nPERCEKT"-;TAB(25); "RETURN",
816 PRINT " TARGET ";TAE(13);"OF TIME";TAB(25);"PERIOD"
817 PRINT "(CO/CL) ";TAB(13);"EXCEEDED";TAE(25);"(YEARS)"
818 PRINT "	«;TAB(13);"	";TA3(25);"	"
820 EEM - LOAD QUAD.  WGTS & ROOTS
83t) GOSUB 1410                         ,
S^K) REI-1 COiffDT PORTION OF Q(X) ARGUMENT IIIDEP OF CO
850 FOR 1=1  TO NO
860 REK - EVALUATE USING  IHV PROB TRANSFORMATION
870 P9£=R5£(I)
880 GOSUE 1310
S90 Z9f(I)sLOG(1-i-EXP(U9-W9sX9))-U3
900 NEXT I
910 EEK - CONG LOOP
520 F.OR CO=B1  TO  E2  STEP  33
930 15=0
SlJO REK - QUAD LOOP  - EVALUATE Q(X)  = F AND SUM
950 FOR 1=1  TO KO
960 X={LOG(CO
970 XO=SGN(2)
980 X=
990 F=1
1000 F=.5*F"(-16)
1010 IF X0>0 THEii GOTO 1030
T020 Fs1-F
1030 I5=I5+F«Z5#(I)
1040 EEXT I
1050 REI« COHPUTE  RETURN PERIOD
1060 10=1/365/15
                                 '  Figure  D-5  (cont'd.)                                       -j
                                            D-16

-------
                                                                      VII (D).
                                                                      Revision  No.  0
  1070 15=100*15         -
  1080 PRINT USING  "###.###      "; CO, 15, 10
  1090 NEXT CO     ••'.-.
  1100 P?:INT CER$(7)
  1101 INPUT "ENTER  TO  CONTINUE,  OR  'STOP'  ";A$     '     x
  1102 IF A$OnSTOPw THEN GOTO 560          ' '"
  1110 £EM GOTO 790                                      ,           .
  1120 FOR L=1 TO 7                     )
  1130 PRINT                                .    '•
  1140 NEXT L                  •
  1 1 45 KEY OH
  1150 END
  1160 HEM SUBROUTINE TO LOAD. NORMAL AND REVERSE NOBK4L COEFFIC^rTS
  11~r LI =.049867347*                                 ?
  1180 D2=..021U10001-'-    •
  'ISO D3a.0032776263#
  1200 D4=3.80036E-05
  1210 D5=4.88906E-05   .  -
  1220 D6s5.383E-06                                    .
  1230 REl'i
  1240  El =2.51 5517  •
  1250  E2=.302S53              "    .
  1260  E3=. 010328  .-.••-,
  1270  £4=1.432788
  1280  S5=. 189269                 '         ,  .
  1290  E6 = . 00.1308
  1300  RETURN      '                 •          ,
  1310  RET: SUBROUTINE  TO  COMPUTE  INVERSE iiOFj-iAL
  1320  BEM POLYNOMIAL  APPROX  TO INVERSE TABLE
  1330  EEF FlICa#)= X-J-(E1-i-E2*X#-hE
  1340  S9=1
  1349  IF P9#<1E-18 THEN  P9#=1E-18
'1350  IF P9#<.5 TKBT  GOTO  1380
  1360  P9,?=i
  1370  S9=-1
  1380  P9i
  1390  X9=FKC(P9#)*S9
  1400 .RETURN
  1410  REM QUADRATURE  SUBROUTINE  - COMPUTE ROOTS AND WEIGHTS
  1420'REM I5=INTEGRAL
  1430  REM R5(KO)= NO  ROOTS
  1440  REM Z5(KO)= NO WEIGHTS
 1450  REM LOAD ROOTS AND WEIGHTS FOR 32ND ORDER QUADS
 1460  REM, FIRST THE GAUSSIAN, THEN THE LAGUERRE TERt-S
 1470  REK QUAD ROOTS & WEIGHTS FOR 16TH ORDER GAUSSIAN
 T480 R1=8                              •
 1490 R#(1) =-.989400935*
 1500 R#(2)=-. 944575023*
 1510 R#(3)=-.8656312024#
 1520 R#(4).=-.7554C44084#          ;
 1530 R#(5)=-.617S762444#  '            ...'•-
'1540 R#(6)=-.458oi67776#
 1550 E#(7) =-.2816035508*                   - -
                                   Figure D-5  (cont'd.)
                                            n_i7

-------
                                                                   VII (D)
                                                                   Revision No.
1560 M( 8)=-:. 09501250984*
1570 S8£(1)=.02715245942*                '              ''
1580 S8£<2)=.06225352394*               •             '"
1590 S8£(3)=.Q9515851168*         '      .                        '   '
1600 S8#(4)=.1246289713*
1610 S8£(5)=.1495959888*     '           .         .        .
1620 S8#(6)=,1691565194*                               •
1630 S8#(7)=.1826034154*
1640 SS*(8)=.1894506105*
1650 KO=4*R1
1660 HEM CONVERT GAUSSIAN HOOTS & WEIGHTS FOR (0,1) INTEGR. INTERVAL
1670 REH AED DIVIDE BY TUO FOR COMPOSITE FORMULA
1680 FOR K2=1 TO R1
1690 R5#(K2)=.5+.5*R*(K2)

1710 Z5*(K2)=S8*(K2)/4                    '  .    .

1730 IffiXT K2
1740 REM LOAD THE LAGUERRE ROOTS AND WEIGHTS, PROPERLY CONVERTED
1750 'REH 16TH ORDER LAGUERRE ROOTS AND WEIGHTS
1760 P#(1)=51.7011603395*
1770 P#(2)=41,9404526477*
1780 P*C3)=34.5833987023*
1790 P£(4)=28.5787297429*
1800 ?£(5)=23.515905694*
1-310 P*(6)=19.1301568568*                                       '   .
1820 P*(7)=15.4415273688*                                '
1830 P#(8)=12.2142233689*
1840 Pf(9)=9.43831433639*   •                              -
1850 P#(10)=7.07033853505*
1860 P*(11)=5.07801861455*
1870 P#(12)=3-43708663389*
1880 P*{13)=2.1292836451*
1390 P*( 1-4) =1.14105777483*
1900 P#{15)=.462696328915*
1910 ?#(16)=.0876494104789*
1920 Q£(1)=4.l6l46237D-22
1930 Q£(2)=5.05C4737D-18
1940 p£(3)=6.297967003D-15             '
1950 Q^C-4) =2.1270790330-12
I960 G£t5)=2.862350243D-10                                     -
1970 Qf-(6) =1.8810248410-08
1980 Q#(7)=.0000006828319331*
1990 Ojf(8)=.0000l484458687*
2000 Q*(9)=.0002042719153*
2010 Q*(10)=,00184907094353*
2020 Q*(11)=.0112999000803*
2030 Q*(12)=.0473289286941*
2040 q*(133=.136296934296*
2050 QS( 1-4) =.265795777644*
2060 Q$(15)s.331057854951*                   '
2070 ^(163=.206151714958*                                    '

                                  Figure D-5 (cont'd.)
                                           D-18

-------
                                                                  VII (D)
                                                                  Revision Wo. 0
2080 FOB K2=1 TO
2090 P.5£(X2+NO/2) =EXP(-P#(K2);
2100 Z5#(K2+NO/2)=Q#(K2)/2
2110 NEXT IC2
2120 RETURN •

A> -      •'   '          •   ' '  •
                    RECEIVING" WATER  CONG  (co)  PROBABILITY DISTRIBUTION
                                   AND  RETURN PERIOD
                              FOR MULTIPLES  OF  TARGET CONG
                               DUE  TO POINT  SOURCE  LOADS
                   *«******««»«««*«$*«*»*»«*»*««****«***«*«***«*«*«£«**
                   POINT SOURCE -  RECEIVING WATER
                   CONCENTRATION ANALYSIS
                   INPUT COEF OF' VAR OF QS,QE,CE
                          RATIO...7Q10/AVGQS
                          RATIO...7Q10/AVGQE
                         .RATIO...AVG CE/CL
                   'BACKGROUND .STREAM CONG (CS) IS ASSUMED TO BE ZERO
                   ENTER COEF OF VAR OF QS,QE,CE
                   71.5.2.7         '
                   ENTER. THE FOLLOWING RATIOS:
                    ... ----- 7Q10/AVG QS ? .05  '
                    ----- '...7Q10/AVG QE ? 3.0
                    .. ------ ..AVG CE/EL ? .67
                   COEF  OF VAR.....QS  =1.5
                   COEF  OF VAR ..... QE  =   .2
                   COEF  OF VAR ---- .CE  =   .7

                          7Q10/AVG QS  =   .05
                          7Q10/AVG QE  =   3
                            :AVG CE/EL  =   .67
                  ENTER LOWEST, HIGHEST, AND INCREMENT  OF MULT OF TARGET FOR
                  WHICH % EXCEED IS DESIRED? 1,5,1
                    COEF OF VAR	QS' =  1.5
                    COEF OF VAR	.QE =  .2
                    COEF OF VAR	CE s  .7

                           7Q10/AVG QS =  .05
                           7Q10/AVG QS a  3
                             AVG CE/EL =  .67.
                                 Figure D-5 (cont'd.)

                                          D-19

-------
VII (D)
Revision No. 0^
STREAM COHC (CO)
MULT OF PERCENT
TARGET OF TIME
(CO/CL) EXCEEDED
1.000 0.894
2.000 0.112
3.000 0.024
4.000 0.007
5.000 0.002
ENTER  TO CONTINUE
COEF OF VAR 	 QS =
COEF OF VAR 	 QE =
COEF OF VAR 	 CE =
7Q10/AVG QS =
7Q10/AVG QE =
AVG CE/EL =

RETURN
PERIOD
(YEARS)
0.306
2.443
11.313
39.429
114.356
, OR 'STOP' ?
1.5.
.2
.7
.05
3
.67
ENTER LOWEST, HIGHEST, 'AND INCREMENT OF MULT OF TARGET FOR
WHICH % EXCEED IS DESIRED? 2.5,3,..1
COEF OF VAE 	 QS' =
COEF OF VAR 	 QE =
COEF OF VAR 	 CE =
7Q10/AVG QS =
7Q10/AVG QE =
AVG CE/EL '=
Jlllllllllllllll'l I1 'I11! TTH
STREAM CONG (CO)
MULT OF PERCENT
TARGET OF TIME
(CO/CL), EXCEEDED
2.500 0.050
2.600 0.043
2.700 0.037 •
2.800 0.032
2.900 0.028
3.000 0.024
1.5
.2
.7
.05
3 •
.67

.RETURN
PERIOD
(YEARS)
5.501 -
6.395
7.410
8.558
9.854 .
11.313
ENTER  TO CONTINUE,  OR 'STOP' ? STOP
            Figure D-5 (cont'd.)

-------