Hydrologic Frequency Analysis PDF
Hydrologic Frequency Analysis PDF
CECW-EH-Y
Engineer Manual
No. 1110-2-1415 5 March 1993
1. Purpose. This manual provides guidance and procedures for frequency analysis of: flood flows, low
flows, precipitation. water surface elevation, and flood damage.
2. Applicabilit}'. This manual applies to major subordinate commands. districts, and laboratories ha,·ing
responsibility for the design of civil works projects.
3. General. Frequency estimates of hydrologic, climatic and economic data are required for the
planning, design and evaluation of flood control and navigation projects. The text illustrates many of
the statistical techniques appropriate for hydrologic problems by example. The basic theory is usually
not provided, but references are provided for those who wish to research the techniques in more detail.
WILLIAM D. BROWN
Colonel, Corps of Engineers
Chief of Staff
DEPARTMENT OF THE ARMY EM 1110-2-1415
U.S. Army Corps of Engineers
Washington, D.C. 20314-1000
CECW-EH-Y
Engineer Manual
No. 1110-2-1415 5 March 1993
Table of Contents
CHAPTER INTRODUCTION
Purpose and Scope . . . . . . . . . . . . . . . . . . . . I- I 1-1
References ......................... . 1-2 1-1
Definitions ......................... . 1-3 1-1
Need for Hydrologic Frequency Estimates .. . 1-4 1-1
Need for Professional Judgement ........ . 1-5 1-2
11
EM 1110-2-1415
5 Mar 93
111
EM 1110-2-1415
5 Mar 93
List of Figures
iv
EM 1110-2-1415
5 Mar 93
List of Tables
v
EM I I 10-2-1415
5 Mar 93
INTRODUCTION
I -1. Purpose and Scope. This manual provides guidance in applying statistical principles
to the analysis of hydrologic data for Corps of Engineers Civil Works activities. The text
illustrates, by example, many of the statistical techniques appropriate for hydrologic
problems. The basic theory is usually not provided, but references are provided for those
who wish to research the techniques in more detail.
1-2. References. The techniques described herein are taken principally from "Guidelines
for Determining Flood Flow Frequency" (46) 1, "Statistical Methods in Hydrology" (1 ), and
"Hydrologic Frequency Analysis" (41 ). References cited in the text and a selected
bibliography of literature pertaining to frequency analysis techniques are contained in
Appendix A.
1
Numbered references refer to Appendix A, Selected Bibliography.
1-1
EM 1110-2 1415
S Mar 93
whether regulated frequency curves will be needed, etc. A frequency study should be
well coordinated with the hydrologist, the planner and the economist.
1-5. Need for Professional Judgment. It is not possible to define a set of procedures that
can be rigidly applied to each frequency determination. There may be applications where
more complex joint or conditional frequency methods, that were considered beyond the
scope of this guidance, will be required. Statistical analyses alone will not resolve all
frequency problems. The user of these techniques must insure proper application and
interpretation has been made. The judgment of a professional experienced in hydrologic
analysis should always be used in concert with these techniques.
1-2
EM 1110-2-1415
5 Mar 93
CHAPTER 2
FREQUENCY ANALYSIS
2-1. Definition.
a. Frequency. Many of the statistical techniques that are applied to hydrologic data
(to enable inferences to be made about particular attributes of the data) can be labeled
with the term "frequency analysis" techniques. The term "frequency" usually connotes a
count (number) of events of a certain magnitude. To have a perspective of the importance
of the count, the total number of events (sample size} must also be known. Sometimes the
number of events within a specified time is used to give meaning to the count, e.g., two
daily flows were this low in 43 years. The probability of a certain magnitude event
recurring again in the future, if the variable describing the events is continuous, (as are
most hydrologic variables), is near zero. Therefore, it is necessary to establish class
intervals (arbitrary subdivisions of the range) and define the frequency as the number of
events that occur within a class interval. A pictorial display of the frequencies within
each class interval is called a histogram (also known as a frequency polygon).
f. • n./N (2-1)
1 1
where:
2-1
EM I I 10-2 1415
5 Mar 93
J
~ -
E ~ -
z
:I
• I--
1- / -
-
- -
.. ~
I--·
~
~
1\ 1,m
-
-
-
-
~ -
0
.. ..... ....
Lopritbm of Annual Pew
1!1.2 1!1.8
•
•c:u o.e
.,•
••u
c3 o.e
c:'
0
...z
~
...
0
o ....
i
A
e
A.
O.R
2-2
EM 111 0-2- 1415
5 Mar 93
/
v
183 /
~
•
'"u.
3
... ~
v
v
D
18 2
II.
...
:I
-4
c•
181 /
v
c~
v
aa.aa aa.s as ae 68 18 1 .1 • 81
2-3
EM 1110-2 1415
S Mar 93
b. l.lm. Duration curves are useful in assessing the general low flow characteristics
of a stream. If the lower end drops rapidly to the probability scale, the stream has a low
ground-water storage and, therefore, a low or no sustained flow. The overall slope of the
flow-duration curve is an indication of the flow variability in the stream. Specific uses
that have been made of duration curves are: 1) assessing the hydropower potential of
run-of-river plants; 2) determining minimum flow release; 3) water quality studies; 4)
sediment yield studies; and 5) comparing yield potential of basins. It must be
remembered that the chronology of the flows is lost in the assembly of data for duration
curves. For some studies, the low-flow sequence, or persistence, may be more important
(see Chapter 4).
2-4
EM 1110-2-1415
5 Mar 93
FISHKILL CREEK
10000
/- MAX
en
""
u
....z
1000
---
.......-- 1%
5%
10%
30%
a
w
a
w
w
u
--
X
w
X
0
50% w
._....
.J ::E
100 70%
"">
._""
90%
.....J 0
0""' 95%
zw
99% u
a:
10 '. ·.· MIN w
Q.
MONTH
2-5
EM 1110-2-1415
5 Mar 93
b. Uniformity of Data.
(I) General Considerations. Data selected for a frequency study must measure the
same aspect of each event (such as peak flow, mean-daily flow, or flood volume for a
specified duration), and each event must result from a uniform set of hydrologic and
operational factors. For example, it would be improper to combine items from old records
that are reported as peak flows but are in fact only daily readings, with newer records
where the instantaneous peak was actually measured. Similarly, care should be exercised
when there has been significant change in upstream storage regulation during the period
of record to avoid combining unlike events into a single series. In such a case, the entire
record should be adjusted to a uniform condition (see Sections 2-3f and 3-9). Data should
always be screened for errors. Errors have been noted in published reports of annual
flood peaks. And, errors have been found in the computer files of annual flood peaks.
The transfer of data to either paper or a computer file always increases the probability
that errors have been accidentally introduced.
(2) Mixed Populations. Hydrologic factors and relationships during a general winter
rain flood are usually quite different from those during a spring snowmelt flood or during
a local summer cloudburst flood. Where two or more types of floods are distinct and do
not occur predominately in mutual combinations, they should not be combined into a
single series for frequency analysis. It is usually more reliable in such cases to segregate
the data in accordance with type and to combine only the final curves, if necessary. In
the Sierra Nevada region of California and Nevada, frequency studies are made separately
for rain floods which occur principally during the months of November through March,
2-6
EM 1110-2-1415
5 Mar 93
and for snowmelt floods, which occur during the months of April through July. Flows for
each of these two seasons are segregated strictly by cause - those predominantly caused by
snowmelt and those predominantly caused by rain. In desert regions, summer
thunderstorms should be excluded from frequency studies of winter rain flood or spring
snowmelt floods and should be considered separately. Along the Atlantic and Gulf Coasts,
it is often desirable to segregate hurricane floods from non hurricane events. Chapter I 0
describes how to combine the separate frequency curves into one relation.
c. Location Differences. Where data recorded at two different locations are to be
combined for construction of a single frequency curve, the data should be adjusted as
necessary to a single location, usually the location of the longer record. The differences in
drainage area, precipitation and, where appropriate, channel characteristics between the
two locations must be taken into account. When the stream-gage location is different
from the project location, the frequency curve can be constructed for the stream-gage
location and subsequently adjusted to the project location.
d. Estimating Missjng Events. Occasionally a runoff record may be interrupted by a
period of one or more years. If the interruption is caused by destruction of the gaging
station by a large flood, failure to fill in the record for that flood would result in a biased
data set and should be avoided. However, if the cause of the interruption is known to be
independent of flow magnitude, the record should be treated as a broken record as
discussed in Section 3-2b. In cases where no runoff records are available on the stream
concerned, it is usually best to estimate the frequency curve as a whole using regional
generalizations, discussed in Chapter 9, instead of attempting to estimate a complete series
of individual events. Where a longer or more complete record at a nearby station exists, it
can be used to extend the effective length of record at a location by adjusting frequency
statistics (Section 3-7) or estimating missing events through correlation (Chapter 12).
e. Climatic Cycles. Some hydrologic records suggest regular cyclic variations in
precipitation and runoff potential, and many attempts have been made to demonstrate that
precipitation or streamflows evidence variations that are in phase with various cycles,
particularly the well established 11-year sunspot cycle. There is no doubt that
long-duration cycles or irregular climatic changes are associated with general changes of
land masses and seas and with local changes in lakes and swamps. Also, large areas that
have been known to be fertile in the past are now arid deserts, and large temperate regions
have been covered with glaciers one or more times. Although the existence of climatic
changes is not questioned, their effect is ordinarily neglected, because the long-term
climatic changes have generally insignificant effect during the period concerned in water
development projections, and short-term climatic changes tend to be self-compensating.
For these reasons, and because of the difficulty in differentiating between stochastic
(random) and systematic changes, the effect of natural cycles or trends during the analysis
period is usually neglected in hydrologic frequency studies.
f. Effect of Basin Development on Frequency Relations.
(1) Hydrologic frequency estimates are often used for some purpose relating to
planning, design or operation of water resources control measures (structural and
nonstructural). The anticipated effects of these measures in changing the rate and volume
of flow is assessed by comparing the without project frequency curve with the
corresponding with project frequency curve. Also, projects that have existed in the past
have affected the rates and volumes of flows, and the recorded values must be adjusted to
reflect uniform conditions in order that the frequency analysis will conform to the basic
2-7
EM 1110-2-1415
5 Mar 93
2-8
EM 1110-2-1415
5 Mar 93
.s
6
5
r;:~ "'
3
1 2 I I
•
g., 1 1111111111
0 11111111111111
1 10 15 20 2"'
Arrayed Sequence Number
.
_g 6
...
.s 5
0
r;:
"3
1 2
•
g., 1 I 11111111111111111 II II I II I Ill II Ill I
0
I Ill 11111111111111 IIIII II lllln-lrrr
5 10 15 20 25 30 40 45 50
Arrayed Sequence Number
2-9
EM 1110-2-1415
5 Mar 93
b. Selection and Arrangement of Peak Flow Data. General principles in the selection
of frequency data are discussed in Section 2-3. Data used in the construction of
frequency curves of peak flow consist of the maximum instantaneous flow for each year
of record (for annual-event curve) or all of the independent events that exceed a selected
base value (for partial-duration curve). This base value must be smaller than any flood
flow that is of importance in the analysis, and should also be low enough so that the total
number of floods in excess of the base equals or exceeds the number of years of record.
Table 2-2 is a tabulation of the annual peak flow data with dates of occurrence, the data
arrayed in the order of magnitude, and the corresponding plotting positions.
(2-2a)
The plotting position for the smallest event (P") is the complement (1-P 1) of this value,
and all the other plotting positions are interpolated linearly between these two. The
median plotting positions can be approximated by
For partial-duration curves, particularly where there are more events than years (N),
plotting positions that indicate more than one event per year can also be obtained using
2-10
EM 1110-2-1415
5 Mar 93
Table 2-2. Annual Peaks, Sequential and Arrayed with Plotting Positions.
Equation 2-2b. This is simply an approximate method used in the absence of knowledge
of the total number of events in the complete set of which the partial-duration data
constitute a subset.
d. Plotting Grids. The plotting grid recommended for annual flood flow events is
the logarithmic normal grid developed by Allen Hazen (ref 13) and designed such that a
logarithmic normal frequency distribution will be represented by a straight line, Figure
2-5. The plotting grid used for stage frequencies is often the arithmetic normal grid. The
plotting grids may contain a horizontal scale of exceedance probability, exceedance
frequency, or percent chance exceedance. Percent chance exceedance (or nonexceedance)
is the recommended terminology.
2-11
EM 1110-2-1415
5 Mar 93
v
...u•
...
J
0
-4
1&. ./
v
113 L
v
.:t
••
A.
-4
•c
::J
c
~
2-12
EM 1110-2-1415
5 Mar 93
2-13
EM 1110-2-1415
5 Mar 93
" ~
'r--.
f'
.•u
~
~
D
..;
L
D c
~~
""•• D
~
~I'll
Q.
~~
..•
u ll
ll
..
L
;
D 1e 4
""•
•
IL /
/c
,_ ~~ ~
I""
~
D
... ...
118.88 88.8 88 1 .1 • e1
2-14
EM 1110-2-1415
5 Mar 93
CHAPTER 3
FLOOD FREQUENCY ANALYSIS
3-1. Introduction.
The procedures that federal agencies are to follow when computing a frequency
curve of annual flood peaks have been published in Guidelines for Determining Flood
Flow Frequency, Bulletin 17B (46). As stated in Bulletin 17B, "Flood events ... do not fit
any one specific known statistical distribution." Therefore, it must be recognized that
occasionally, the recommended techniques may not provide a reasonable fit to the data.
When it is necessary to use a procedure that departs from Bulletin 17B, the procedure
should be fundamentally sound and the steps of the procedure documented in the report
along with the frequency curves.
This report contains most aspects of Bulletin 17B, but in an abbreviated form.
Various aspects of the procedures are described in an attempt to clarify the computational
steps. The intent herein is to provide guidance for use with Bulletin 17B. The step by
step procedures to compute a flood peak frequency curve are contained in Appendix 12 of
Bulletin 17B and are not repeated herein.
3-2. Log-Pearson Type III Distribution.
a. General. The analytical frequency procedure recommended for annual maximum
streamflows is the logarithmic Pearson type III distribution. This distribution requires
three parameters for complete mathematical specification. The parameters are: the mean,
or first moment, (estimated by the sample mean, X); the variance, or second moment,
(estimated by the sample variance, S2); the skew, or third moment, (estimated by the
sample skew, G). Since the distribution is a logarithmic distribution, all parameters are
estimated from logarithms of the observations, rather than from the observations
themselves. The Pearson type III distribution is particularly useful for hydrologic
investigations because the third parameter, the skew, permits the fitting of non-normal
samples to the distribution. When the skew is zero the log-Pearson type III distribution
becomes a two-parameter distribution that is identical to the logarithmic normal (often
called log-normal) distribution.
b. Fitting the Distribution.
(1) The log-Pearson type III distribution is fitted to a data set by calculating the
sample mean, variance, and skew from the following equations:
}:X
X • -- (3-1)
N
2 Lx 2 L(X-X) 2
s = = (3-2a)
N-1 N-1
3-1
EM 1110-2-1415
5 Mar 93
. (3-2b)
N-1
N(L (X-X)3 )
G • = (3-3a)
(N-l){N-2)53 (N-I)(N-2)53
in which:
X = mean logarithm
X == logarithm of the magnitude of the annual event
N • number of events in the data set
52 = unbiased estimate of the variance of logarithms
X - X-X, the deviation of the logarithm of asingle event from the mean
logarithm
G • unbiased estimate of the skew coefficient of logarithms
The precision of the computed values is more sensitive to the number of significant digits
when Equations 3-2b and 3-3b are used.
(2) In terms of the frequency curve itself, the mean represents the general
magnitude or average ordinate of the curve, the square root of the variance (the standard
deviation, 5) represents the slope of the curve, and the skew represents the degree of
curvature. Computation of the unadjusted frequency curve is accomplished by computing
values for the logarithms of the streamflow corresponding to selected values of percent
chance exceedance. A reasonable set of values and the results are shown in Table 3-1.
The number of values needed to define the curve depends on the degree of curvature (i.e.,
the skew). For a skew value of zero, only two points would be needed, while for larger
skew values all of the values in the table would ordinarily be needed.
(3) The logarithms of the event magnitudes corresponding to each of the selected
percent chance exceedance values are computed by the following equation:
log Q = X+ K5 (3-4)
where X and 5 are defined as in Equations 3-1 and 3-2 and where
3-2
EM 1110-2-1415
5 Mar 93
log Q = logarithm of the flow (or other variable) corresponding to a specified value
of percent chance exceedance
K = Pearson type III deviate that is a function of the percent chance exceedance
and the skew coefficient.
c. Examole Computation.
( 1) As shown i.!l. the following example, Equation 3-4 is solved by using the
computed values of X and S and obtaining from Appendix V -3 the value of K
corresponding to the adopted skew, G, and...!}le selected percent chance exceedance (P).
An example computation for P-= 1.0, where X, S and G are taken from Table 3-1, is:
Q • 11500 cfs
(2) It has been shown (36) that a frequency curve computed in this manner is biased
in relation to average future expectation because of uncertainty as to the true mean and
standard deviation. The effect of this bias for the normal distribution can be eliminated
by an adjustment termed the expected probability adjustment that accounts for the actual
sample size. This adjustment is discussed in more detail in Section 3-4. Table 3-1 and
Figure 3-1 shows the derived frequency curve along with the expected probability
adjusted curves and the 5 and 95 percent confidence limit curves.
3-3
EM 1110-2-1415
5 Mar 93
--. -.
Computed Frequency Curve
With Expected Probability Adjuatment
.,<+>
~
Confidence Limit Curvea
•u
v
/
~
18 4
v
3
....D
1&.
.:t
••
/
0..
....
•:J A!f!"'
c 18 3
c
q
~-
99.99 99.9 99 98 68 18 1 .1 . 81
d. Broken Record. A broken record results when one or more years of annual peaks
are missing for any reason not related to the flood magnitude. In other words, the missing
events were caused by a random occurrence. The gage may have been temporarily
discontinued for budgetary or other reasons. The different segments of the record are
added together and analyzed as one record, unless the different parts of the record are
considered non-homogeneous. If a portion of the record is missing because the gage was
destroyed by a flood or the flood was too low to record, then the observed events should
be analyzed as an incomplete record.
e. Incomplete Record. An incomplete record can result when some of the peak flow
events were either too high or too low. Different analysis procedures are recommended
for missing high events and for missing low events. Missing high events may result from
the gage being out of operation or the stage exceeding the rating table. In these cases,
every effort should be made to obtain an estimate of the missing events. Missing low
floods usually result when the flood height is below the minimum reporting level or the
bottom of a crest stage gage. In these cases, the record should be analyzed using the
conditional probability adjustment described in Appendix 5 of Bulletin 17B and Section
3-6 of this report.
3-4
EM 1110-2-1415
5 Mar 93
f. Zero-flood years. Some of the gaging stations in add regions record no flow for
the entire year. A zero flood peak precludes the normal statistical analysis because the
logarithm of zero is minus infinity. In this case the record should be analyzed using the
conditional probability adjustment described in Appendix 5 of Bulletin 17B and Section
3-6 of this report.
g. Outliers.
(I) Guidance. The Bulletin 17B (46) defines outliers as "data points which depart
significantly from the trend of the remaining data." The sequence of steps for testing for
high and low outliers is dependent upon the skew coefficient and the treatment of high
outliers differs from that of low outliers. When the computed (station) skew coefficient is
greater than +0.4, the high-outlier test is applied first and the adjustment for any high
outliers and/or historic information is made before testing for low outliers. When the
skew coefficient is less than -0.4, the low-outlier test is applied first and the adjustment
for any low outlier(s) is made before testing for high outliers and adjusting for any
historic information. When the skew coefficient is between -0.4 and +0.4, both the high-
and low-outlier tests are made to the systematic record (minus any zero flood events)
before any adjustments are made.
(3-5)
where:
X = mean logarithm (may have been adjusted for high or low outliers, and/or
historical information depending on skew coefficient)
(3) High Outliers. Flood peaks that are above the upper threshold are treated as high
outliers. The one or more values that are determined to be high outliers are weighted by
the historical adjustment equations. Therefore, for any flood peak(s) to be weighted as
high outlier(s), either historical information must be available or the probable occurrence
of the event(s) estimated based on flood information at nearby sites. If it is not possible
to obtain any information that weights the high outlier(s) over a longer period than that of
the systematic record, then the outlier(s) should be retained as part of the systematic
record.
3-5
EM 1110-2-1415
5 Mar 93
(4) Low Outliers. Flood peaks that are below the low threshold value are treated as
low outliers. Low outliers are deleted from the record and the frequency curve computed
by the conditional probability adjustment (Section 3-6). If there are one or more values
very near, but above the threshold value, it may be desirable to test the sensitivity of the
results by considering the value(s) as low outlier(s).
( 1) Definitions. Historic events are large flood peaks that occurred outside of the
systematic record. Historical information is knowledge that some flood peak, either
systematic or historic, was the largest event over a period longer than that of the
systematic record. It is historical information that allows a high outlier to be weighted
over a longer period than that of the systematic record.
(2) Eauations. The adjustment equations are applied to historic events and high
outliers at the same time. It is important that the lowest historic peak be a fairly large
peak, because every peak in the systematic record that is equal to or larger than the lowest
historic peak must be treated as a high outlier. Also a basic assumption in the adjusting
equations is that no peaks higher than the lowest historic event or high outlier occurred
during the unobserved part of the historical period. Appendix D in this manual is a
reprint of Appendix 6 from Bulletin 17B and contains the equations for adjusting for
historic events and/or historical information.
MSEG(G) + MSEG(G)
G\J = (3-6)
MSEG + MSEG
where:
G\1 = weighted skew coefficient
G = computed (station) skew
G = generalized skew
MSEG = mean-square error of generalized skew
MSEG = mean-square error of computed (station) skew
3-6
EM 1110-2-1415
5 Mar 93
(I) The mean-square error of the computed skew coefficient for log-Pearson type
III random variables has been obtained by sampling experiments. Equation 6 in Bulletin
17B provides an approximate value for the mean-square error of the computed (station)
skew coefficient:
(3-7b)
where:
Appendix F-1 0 provides a table of mean-square error for several record lengths and skew
coefficients based on Equation 3-7a.
(2) The mean-square error (MSE) for the generalized skew will be dependent on the
accuracy of the method used to develop generalized skew relations. For an isoline map,
the MSE would be the average of the squared differences between the computed (station)
skew coefficients and the isoline values. For a prediction equation, the square of the
standard error of estimate would approximate the MSE. And, if an arithmetic mean of
the stations in a region were adopted, the square of the standard deviation (variance)
would approximate the MSE.
3-7
EM 1110-2-1415
5 Mar 93
the confidence limit curves around a frequency curve that has the following assumed
statistics: N= 10, X•O., S= 1.0.
b. If one wished to design a flood protection work that would be exceeded, on the
average, only one time every 100 years (one percent chance exceedance), the usual design
would be based on the normal standard deviate of 2.326. Notice that there is a 0.5 percent
chance that this design level may come from a "true" curve that would average 22
exceedances per 100 years. On the other side of the curve, instead of the expected one
exceedance, there is a 99.5 percent chance that the "true" curve would indicate 0.004
exceedances. Note the large number of exceedances possible on the left side of the curve.
This relationship is highly skewed towards the large exceedances because the bound on the
right side is zero exceedance. A graph of the number of possible "true" exceedances
versus the probability that the true curve exceeds this value, Figure 3-3, provides a curve
with an area equal to the average (expected) number of exceedances.
c. The design of many projects with a target of 1 exceedance per 100 years at each
project and assuming N= I 0 for each project, would actually result in an average of 2.69
exceedances (see Apper. . ix F-8).
4 .50
....
~
u
z
UJ
-
c
1-
< ~
u
.75 :
>
w ~
c
l .90
c !
a:
< .950 :
c .975
z
< .990
1-
en .995
3-8
EM 1110-2-1415
5 Mar 93
Z5
I I I I
TARGET Of 1 EXCEEDA~CE PER 100 EVE~TS
- AND IASEC ON SAMPlE SIZE OF 10 EVENTS
zo
en
cr
...:
w
>-
-
0
0 15
-
\
cr
w
a..
en
w
(.,)
z
c 10
c
w
w
(.,)
><
w
\ AREA UNDER CURVE IS EXPECTED NUMBER
OF EXCEEDANCES PER 100 EVENTS
5
~-....... 2.69 (Appendix V·8l
---
-
- ~
r--..,
0
0.0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1.0
PROBABILITY OF EXCEEDING
Figure 3-3. Cumulative Probability Distribution of Exceedances per 100 Years.
d. There are two methods that can be used to correct (expected probability
adjustment) for this bias. The first method, as described above, entails plotting the curve
at the "expected" number of exceedances rather that at the target value, drawing the new
curve and then reading the adjusted design level. Appendix F-8 provides the percentages
for the expected probability adjustment.
e. The second method is more direct because an adjusted deviate (K value) is used in
Equation 3-4 that makes the expected probability adjustment for a given percent chance
exceedance. Appendix F-7 contains the deviates for the expected probability adjustment.
These values may be derived from the t-distribution by the following equation:
where:
3-9
EM 1110-2-1415
5 Mar 93
N! (3-9)
I!(N-I)!
where:
R1 • risk (probability) of experiencing exactly I flood events
N = number of years (trials)
I .. number of flood events (successes)
P • exceedance probability, percent chance exceedance divided by 100, of the
annual event (probability of success)
(The terms in parentheses are those usually used in statistical texts)
3-10
EM 1110-2-1415
5 Mar 93
and the probability of experiencing one or more floods is easily computed by taking the
complement of the probability of no floods:
D
n(1 or more)
,. ,I -,1 ,I _ 'P, 1\N
c. Aoplication.
= 1-.9930
= 1-.74
.. .26 or 26% chance
(2) Appendix F-12 provides a table for risk as a function of percent chance
exceedance, period length and number of exceedances. This table could also be used to
check the validity of a derived frequency curve. As an example, if a frequency curve is
determined such that 3 observed events have exceeded the derived 1% chance exceedance
level during the 50 years of record, then there would be reason to question the derived
frequency curve. From Appendix F-12, the probability of this occurring is 0.0122 or
about 1%. It is possible for the situation to occur, but the probability of occurring is very
low. This computation just raises questions about the validity of the derived curve and
indicates that other validation checks may be warranted before adopting the derived
curve.
3-6. Conditional Probability Adjustment. The conditional probability adjustment is made
when flood peaks have either been deleted or are not available below a specified
truncation level. This adjustment will be applied when there are zero flood years, an
incomplete record or low outliers. As stated in Appendix 5 of Bulletin 17B, this
procedure is not appropriate when 25 percent or more of the events are truncated. The
computation steps in the conditional probability adjustment are as follows:
1. Compute the estimated probability (P) that an annual peak will exceed the
truncation level:
P = N/n (3-lla)
3-11
EM 1110-21415
5 Mar 93
where N is the number of peaks above the truncation level and n is the total number
of years of record. If the statistics reflect the adjustments for historic information,
then the appropriate equation is
H- WL
p .. (3-11 b)
H
where H is the length of historic period, W is the systematic record weight and L is
the number of peaks truncated.
2. The computed frequency curve is actually a conditional frequency curve. Given
that the flow exceeds the truncation level, the exceedance frequency for that flow
can be estimated. The conditional exceedance frequencies are converted to annual
frequencies by the probability computed in Step 1:.
(3-12)
where P is the annual percent chance exceedance and Pd is the conditional percent
chance exceedance.
3. Interpolate either graphically or mathematically to obtain the discharge values
(QP) for I, 10 and 50 percent chance exceedances.
4. Estimate log-Pearson type III statistics that will fit the upper portion of the
adjusted curve with the following equations:
log (Q 1/Q 10 )
G5 • -2.50 + 3.12 - - - - - (3-13)
log (Q 1/Q 50 )
ss .. (3-14)
(3-15)
where G 5 , S5 and X are the synthetic skew coefficient, standard deviation and
mean, respectively; Q1, Q1D and Q50 the discharges determined in Step 3; and K 1 and
K 50 are the Pearson Type III deviates for percent change exceedances of I and 50
and skew coefficient G 5 •
3-12
EM 1110-2-1415
5 Mar 93
5. Combine the synthetic skew coefficient with the generalized skew by use of
Equation 3-6 to obtain the weighted skew.
6. Develop the computed frequency curve with the synthetic statistics and compare
it with the plotted observed flood peaks.
3-7. Two-Station Comparison.
a. Purpose.
(I) In most cases of frequency studies of runoff or precipitation there are locations
in the region where records have been obtained over a long period. The additional period
of record at such a nearby station is useful for extending the record at a short record
station provided there is reasonable correlation between recorded values at the two
locations.
(3) The procedure for adjusting the statistics at a short-record station involves three
steps: (1) computing the degree of correlation between the two stations, (2) using the
computed degree of correlation and the statistics of the longer record station to compute
an adjusted set of statistics for the shorter-record station, and (3) computing an equivalent
"length of record" that approximately reflects the "worth" of the adjusted statistics of the
short-record station. The longer record station selected for the adjustment procedure
should be in a hydrologically similar area and, if possible, have a drainage area size
similar to that of the short-record station.
(3-16)
where:
3-13
EM 1110-2 1415
5 Mar 93
For most studies involving streamflow values, it is appropriate to use the logarithms of the
values in the equations in this section.
c. Adjustment of Mean. The following equation is used to adjust the mean of a
short-record station on the basis of a nearby longer-record station:
(3-17)
where:
Y • the adjusted mean at the short-record station
Y 1 == the mean for the concurrent record at the short-record station
X 3 == the mean for the complete record at the longer-record station
X1 • the mean for the concurrent record at the longer-record station
R = the correlation coefficient
the standard deviation for the concurrent record at the short-record station
the standard deviation for the concurrent record at the longer-record
station
All of the above parameters may be derived from the logarithms of the data where
appropriate, e.g., for annual flood peaks. The criterion for determining if the variance of
the adjusted mean will likely be less than the variance of the concurrent record is:
where N 1 equals the number of years of concurrent record. If R 2 is less than the
criterion, Equation 3-17 should not be applied. In this case just use the computed mean at
the short-record station or check another nearby long-record station. See Appendix 7 of
Bulletin 17B for procedures to compare the variance of the adjusted mean against the
variance of the entire short-record period.
d. Adjustment of Standard Deviation. The following equation can be used to adjust
the standard deviation:
(approximate)
3-14
EM 1110-2-1415
5 Mar 93
where:
Sy
1
= the standard deviation for the period of concurrent record at the short-record
station
Sx = the standard deviation for the complete record at the base station
Sy
1
= the standard deviation for the period of concurrent record at the base station
All of the above parameters may be derived from the logarithms of the data where
appropriate, e.g., for annual flood peaks. This equation provides approximate results
compared to Equation 3- I 9 in Appendix 7 of Bulletin 17B, but in most cases the
difference in the results does not justify the additional computations.
f. Equivalent Record Length. The final step in adjusting the statistics is the
computation of the "equivalent record length" which is defined as the period of time
which would be required to establish unadjusted statistics that are as reliable (in a
statistical sense) as the adjusted values. Thus, the equivalent length of record is an
indirect indication of the reliability of the adjusted values of Y and Sy. The equivalent
record length for the adjusted mean is computed from the following equation:
Ny
Ny"" _ _ _ _ _ _ _ _ _._ _ _ _ _ _ _ __
(3-20)
2 2
I - [(Nx- Ny )/NxJ [R - (l - R )/(Ny - 3)]
1 1
where:
Figure 3-4 shows the data and computations for a two-station comparison for a short
record station with 21 events and a long record station with 60 systematic events. It can
be seen that the adjustment of the frequency statistics provides an increased reliability in
the mean equivalent to having an additional 17 years of record at the short-record station.
3-15
EM 1110-2-1415
5 Mar 93
Year
River
flow
-1915- -----------
12600
River
--
Chattooge hllulell
-- -- r-~otal-~~~::~~~--=~~~::
Systeutic Events 60 21 39
STATISTICS OF DATA
long ~ecord
--------
Not
;r-- ort
ord
----
21
1917
1918
14000
5900 Historic Period 71 0 0
------0
1919 16000
1920 8200 log lleen 3.8866 3.9310 3.8650 3.4796
1921
1922
4100
6200 Standard Deviation 0.2322 0.2678 0.2140
---
0.2609
1923
1924
5300
9200 CCM~PUted 0.5015 0.1588 0. 7340
--
·0. 0811
1925 3900 Skew Generalized . - . .
1926 6200 Adopted . - . ·0.1000
1927 3600
1928 20100
1929 11400 CQMPUTATIONS FOR TWO-STATION COMPARISON:
1940 29000•
1941 7530 Slope: b • 0.790495
1942 6870 Correlation Coefficient: R •0.811351 CSee Figure 9·01>
1943 6870
1944 3840 Fr~ equations 19, 21, and 22:
1945 2930
1946 6650 '1 • 3.4796 • 0.790495 (3.8866 • 3.9310)
1947 6440 ,, • 3.4445
1948 12400
1949 13900
1950 4740 s,2 • ((0.2322)2 - (0.2678) 2 Jc0.8114) 2 c0.2609/0.2678> 2
1951 5220 .. (0.2609)2
1952 13400
1953 ~020 s, • 0.23~
1954 6230
1955 5820 1 --~~-71"11'1'""
II --.....---oj,---.:2..:. • 38
1956 5820 ,, • - !!.;_1 (0.6583 - - 0.6583,
1957 5820 I I 21 " l
1958 5620
1959 5620
1960 6440 ii • o.o IISEG • 0.302 G • ·0.0811 MSEG • 0.142
1961 7310
1962 9660 G (0.302)(·0.0811) + (0.0142)(0.0)
1963 5420 w• o.3o2. o.142
1964 9880
1965 2noo 7440 Gw. -0.055 -~ -0.1
1966 13400 5140
1967 15400 2800
1968 5620 3100 FREQUEIICY CURVE, TALLULAH RIVER NEAR CLAYTON, GA
1969 14700 2470
1970 3480 2010 •••••• FLOW,CFS •••••••• PERCE liT ..COIIFIDENCE LIMITS ••
1971 3290 976 EXPECTED CHANCE
19n 7440 2160 COMPUTED PIIIOIAIILITY EXCEEDANCE .05 LIMIT .95 LIMIT
1973 19600 8500
1974 6400 4660 12700 14200 .2 18800 9580
1975 6340 2410 10900 11900 .5 15600 8410
1976 18500 6530 9590 10300 1.0 13400 7540
1977 13000 3580 8350 8800 2.0 11300 6680
1978 7850 4090 6760 6990 5.0 8810 5550
1979 14800 6240 5590 5710 10.0 7040 4680
1980 10900 2880 4430 4480 20.0 5370 3780
1981 4120 1600 2810 2810 50.0 3260 2420
1982 5000 1960 1760 1740 80.o 2060 1450
1983 7910 3260 1370 1340 90.0 1640 1080
1984 4810 2000 1110 1070 95.0 1360 847
1985 4740 1010 745 686 99.0 957 524
• Histor-Ic inforution, peak Iargen since 1915.
3-16
EM 1110-2-1415
5 Mar 93
(Figure 9-1 shows the computations for b and R and Figure 9-2 shows the Tallulah River
annual peaks plotted against the Chattooga River peaks.) Figure 3-5 shows the resulting
unadjusted and adjusted frequency curves based on the computed and adjusted statistics
in Figure 3-4. Although N is actually the equivalent years of record for the mean, the
value is used as an estimateyequivalent record length in the computation of confidence
limits and the expected probability adjustment.
(2) Compute Y 1 and Sy for the entire record at the short-record station.
1
(3) Compute X and Sx for the entire record at the longer-record station.
(4) Compute X 1 and Sx for the portion of the longer-record station which is
concurrent with the shott-record station.
(6) Com put£' Y and Sy for the short-record station using Equations 3-17 and 3-18.
.,
.,
...•u 184
., l.,...- ~
3' I" ~
.;t ~
••
Q. ~ ~
... ~
fi""
•:::Jc 18 3 ~
c
G:
,
3-17
EM 1110-2-1415
5 Mar 93
(7) Calculate the equivalent length of record of the mean for the short-record station
using Equation 3-20.
(8) Compute the frequency curve using adjusted values of Y and S in Equation 3-4
and K values from Appendix F-2 corresponding to the adopted skew coefficient.
(9) Compute the expected probability adjustment and the confidence limits.
a. Nature and Purpose. Flood volume frequency studies involve frequency analysis
of maximum runoff within each of a set of specified durations. Flood volume-duration
data normally obtained from the USGS WATSTORE files consists of data for I, 3, 7, 15,
30, 60, 90, 120, and 183 days. These same values are the default values in the HEC
computer program ST ATS (Table 3-2). Runoff volumes are expressed as average flows in
order that peak flows and volumes can be readily compared and coordinated. Whenever it
is necessary to consider flows separately for a portion of the water year such as the rain
season or snowmelt season, the same durations (up to the 30-day or 90-day values) are
selected from flows during that season only. Flood volume-duration curves are used
primarily for reservoir design and operation studies, and should generally be developed in
the design of reservoirs having flood control as a major function.
3-18
EM 1110-2-1415
5 Mar 93
XA > Xa
XA + (-3.719) SA > X 8 + (-3.719) S8
(SB - SA)
> 0.269
(XB- XA)
where:
XA == Value of frequency curve- A at 99.99 percent chance exceedance
X 8 == Value of frequency curve Bat 99.99 percent chance exceedance
XA = Mean of frequency curve A
X8 = Mean of frequency curve B
SA = Standard deviation of frequency curve A
S8 = Standard deviation of frequency curve B
3-19
EM 1110-2-1415
5 Mar 93
0.30
~
> It
<
~
0.26 . - - - - - - - - - 4 -
i ..,~-~.-
"'g-
-l I
0
c:
.2 ,..
s> 0.20 <
t! ~
C)
<')
"'0
jIll
(i.j 0.15
0.10
2.4 2.8 2.8 3.0 3.2 3.4
1.0
•
/
s-
Ill
t:ll
0.5
v
- /.
0
....J
0
v
1 0.0
~
8
0
3:
G)
l}lj / •
~.5
v
-1.0
/
2.4 2.8 2.8 3.0 3.2 3.-4
3-20
EM 1110-2-1415
5 Mar 93
(2) When the skew changes between durations, it is probably easiest to adopt
smoothed relations for the standard deviation and skew and input the statistics into a
computer program that computes the ordinates. The curves can then be inspected for
consistency.
(3) If the statistics for the peak flows have been computed according to the
procedures in Bulletin 178, the smoothing relations should be forced through those points.
The procedure for computing a least-squares line through a given intersection can be
found in texts describing regression analyses.
d. Frequency Curves for Comorehensive Series.
( 1) General Procedure. Frequency curves of flood volumes are computed
analytically using general principles and methods of Chapters 2 and 3. They should also
be shown graphically and compared with the data on which they are based. This is a
general check on the analytic work and will ordinarily reveal any inconsistency in data
and methodology. The computed frequency curves and the observed data should be
plotted on a single sheet for comparison purposes, Figure 3-7.
(f)
~
u PEAK
1-0AY
.z.
-""
0
1-
a:
10~ 3-DAY
7-DAY
:::>
0 15-DAY
a:
0 30-0AY
~
3: 60-0AY
0
-1 10 3
u...
z
""
UJ
::E
99 90 50 10 •1 . 01
PERCENT CHANCE EXCEEDANCE
3-21
EM 1110-2-1415
5 Mar 93
(2) Interpolation Between Fixed Durations. The runoff volume for any specified
frequency can be determined for any duration between 1-day and 365-days by drawing a
curve on logarithmic paper relating mean discharge (or volume) to duration for that
specified frequency (see Figure 3-8a). When runoff volumes for durations shorter than 24
hours are important, special frequency studies should be made. These could be done in
the same manner as for the longer durations, using skew coefficients interpolated in some
reasonable manner between those used for peak and 1-day flows.
e. Applications of Flood Volume-Duration Frequencies.
(I) Volume-duration Curves. The use of flood volume-duration frequencies in
solving reservoir planning, design, and operation problems usually involves the
construction of volume-duration curves for specified frequencies. These are drawn first
on logarithmic paper for interpolation purposes, as illustrated on Figure 3-8a. The mean
discharge values are multiplied by appropriate durations to obtain volumes and are then
replotted on an arithmetic grid as shown on the Figure 3-8b. A straight line on this grid
represents a constant rate of flow. The straight line represents a uniform flow of 1,500
cfs, and the maximum departure from the 2% chance exceedance curve demonstrates that
a reservoir capacity of 16,000 cfs-days (31,700 acre-feet) is required to control the
indicated runoff volumes by a constant release of 1,500 cfs. The curve also indicates that
a duration of about 8 days is critical for this project release rate and associated
flood-control storage space.
(2) Apolication to a Sjngle Reservoir. In the case of a single flood-control reservoir
located immediately upstream of a single damage center. the volume frequency problems
are relatively simple. A series of volume-duration curves, similar to that shown on Figure
3-8, corresponding to selected exceedance frequencies should first be drawn. The project
release rate should be determined, giving due consideration to possible channel
deterioration, encroachment into the flood plain, and operational contingencies. This
procedure can be used not only as an approximate aid in selecting a reservoir capacity, but
also as an aid in drawing filling-frequency curves.
(3) Application to a Reservoir System. In solving complex reservoir problems,
representative hydrographs at all locations can be patterned after one or more past floods.
The ordinates of these hydrographs can be adjusted so that their volumes for the critical
durations will equal corresponding magnitudes at each location for the selected frequency.
A design or operation scheme based on regulation of such a set of hydrographs would be
reasonably well balanced. Some aspects of this problem are described in Section 3-9g.
3-9. Effects of Flood Control Works on Flood Frequencies.
a. Nature of the Problem. Flood control reservoirs are designed to substantially
affect the frequency of flood flows (or flood stages) at various downstream locations.
Many land use changes such as urbanization, forest clearing, etc. can also have significant
effects on downstream flood flows (see Section 3-1 0). Channel improvements (intended to
reduce stages) and levee improvements (intended to confine flows) at specified locations
can substantially affect downstream flows by eliminating some of the natural storage
effects. Levees can also create backwater conditions that affect river stages for a
considerable distance upstream. The degree to which flows and stages are modified by
various flood control works or land use changes can depend on the timing, areal
distribution and magnitude of rainfall (and snowmelt, if pertinent) causing the flood.
Accordingly, the studies should include evaluations of the effects on representative flood
3-22
EM 1110-2-1415
5 Mar 93
___
j_
~ ~+ 1-
-. l -- _j
, --r- r--.!'t-- I i
I i I
:
r- 1-r-- ~--~ I
..."'u I I
~~
II
~ ~ZI
I
,;
. i
0
...
...J sol
t-- --- - t "·~ j--0--
......
z
:E I
I
I I I
!
I I I
I l 10
i
I
ii
DURATION. DAYS
,."' / ~
.c.t.
60000
/ ,.,t./ i-1-'i.~
"'>-...
0 ~
l,4v"
v. .,.
I
...u
"'
'
/
•oooo r:.<:)
..; //
;F
:E
::l •4...";;
...J
0
>
/
. . . ;.,..;t-
//: .. fc,."","' "-"'
STDIAGE IEOUIIEMEMT COMPUTATION
21000 1200 • 16000 CfS·DATS_
1/)-
20000
' ' 01 31700 AC · F T
'
0 '
10 zo ]0 50 60
DURATION. DAYS
3-23
EM 1110-2-1415
5 Mar 93
events. with careful consideration given to the effects of different temporal and areal
distributions.
b. Terminology.
(I) Natural Conditions. Natural conditions in the drainage basin are defined as
hydrologic conditions that would prevail if no regulatory works or other works of man
were constructed. Natural conditions, however, include the effects of natural lakes,
swamp areas, etc.
(2) Present Conditions. Present or base conditions are defined as the conditions that
exist as of the date of the study or some specified time.
(3) Without- Project Conditions. Without-project conditions are defined as the
conditions that would exist without the projects under consideration, but with all existing
projects and may include future projects whose construction is imminent.
(4) With-Project Conditions. With-project conditions are defined as the conditions
that will exist after the projects under consideration are completed.
c. Reservoir-Level Frequency Computation.
(I) Factors to be Considered. Factors affecting the frequency of reservoir levels
include historical inflow rates and anticipated future inflow rates estimated by
volume-frequency studies, the storage-elevation curves, and the plan of reservoir
regulation including location and size of reservoir outlets and spillway. A true frequency
curve of annual maxima or minima can only be computed when the reservoir completely
fills every year. Otherwise, the events would not be independent. If there is dependence
between annual events, the ordinate should be labeled "percent of years exceeded" for
maximum events and "percent of years not exceeded" for minimum events.
(2) Comoutation and Presentation of Results. A frequency curve of annual
maximum reservoir elevations (or stages) is ordinarily constructed graphically. using
procedures outlined in Section 2-4. Observed elevations (or stages) are used to the extent
that these are available, if the reservoir operation will remain the same in the future.
Historical and/or large hypothetical floods may also be routed through the reservoir using
future operating plans. A typical frequency curve is illustrated on Figure 6-4.
Elevation-duration curves are constructed from historical operation data or from routings
of historical runoff in accordance with procedures discussed in Section 2-2, Figure 3-9.
Such curves may be constructed for the entire period of record or for a selected wet
period or dry period. For many purposes, particularly recreation uses, the seasonal
variation of reservoir elevation (stages) is important. In this case a set of frequency or
duration curves for each month of the year may be valuable. One format for presenting
this information is illustrated on Figure 3-10.
3-24
EM 1110-2-1415
5 Mar 93
TOP OF liiDIJCII)
SUIQIA&CE
TOP OF FL.OOD
CONTROl. POOL
------
'
1030
z0
...
....c
....>-1 ' ___ _ j _____ ..
1020 ; --- --- i
.... 1
I
I
TOP OF VATU
SUPPLY POOl.
......... .. 10
PERCENT OF TIME EXCEEDED
10 10
LL ,, "" 1\
- "' "-
334
333
............Vj' \,
....... Percent of
v/ ~ ~
~ .... Daya Exoeedlld
"'
t: ~
e ~~
~ c:--~
~ ..,....., ~ A
a 0.1'1.
-r--
331 + 1.0'4
c(
~ / -...::::: ...........
~
~ rj/
>
,.,.,. --- --
w 330
...........
al ~ ...;
~~
328
v ~ ......._-.......,
""
0 5.0'4
321 ........... 1:>. 10.0'1.
~ ~
~
3'L7 ~ X 20.0'1.
328 ~
"' v 50.0'1.
3-25
EM 1110-2-1415
5 Mar 93
(3) Incidental Control by Water Supoly Space. In constructing fre :-mency curves of
regulated flows, it must be recognized that reservoir operation for purr es other than
flood control will frequently provide incidental regulation of floods. H-wever, the
availability of such space cannot usually be depended upon, and its value is considerably
diminished for this reason. Consequently, the effects of such space on the reduction of
floods should be estimated very conservatively.
(4) Allowance for Operational Contingencies. In constructing frequency curves of
regulated flows, it should be recognized that actual operation is rarely perfect and that
releases will frequently be curtailed or diminished because of unforeseen operation
contingencies. Also, where flood forecasts are involved in the reservoir operation, it must
be recognized that these are subject to considerable uncertainty and that some allowance
for uncertainty will be made during operation. In accounting for these factors, it will be
found that the actual control of floods is somewhat less than could be expected if full
release capacities and downstream channel capacities were utilized efficiently and if all
forecasts were exact.
(5) Runoff from Unregulated Areas. In estimating the frequency of runoff at a
location that is a considerable distance downstream from one or more reservoir projects, it
must be recognized that none of the runoff from the intermediate areas between the
reservoir(s) and the damage center will be regulated. This factor can be accounted for by
constructing a frequency curve of the runoff from the intermediate area, and using this
curve as an indicator of the lower limit for the curve of regulated flows. Streamflow
3-26
EM 1110-2-1415
5 Mar 93
~ •u "
...-"Spillway Duign Flood
116
,-
......
0
Standard Project Flood
/
.X.
•,•
.-[7
a. . ~
.•
...• u6 ""
"
.. "'"
~
&
•
a:
" ~~
" all ~"'""
. ap
" c
"
" ""
u"'
.X. 1.16
•
a.•
...
•c
J
c
~
3-27
EM 1110-2-1415
5 Mar 93
routing and combining of both the flows from the unregulated area and those from the
regulated area is the best procedure for deriving the regulated frequency curve.
e. Effects of Channel. Levee and Floodway Improvements. The effect of channel,
levee and floodway improvements on river stages at the project location and on river
discharges downstream from the project location can generaJiy be evaluated by routing
several typical floods through the reaches of the improvement and the upstream reaches
affected by backwater. The stages or discharges thus derived can be plotted against
corresponding without- project values, and a smooth curve drawn. This curve could be
used in conjunction with a frequency curve of without-project values to construct a
frequency curve of with-project values as discussed in Paragraph 3-09d( 1)b.
Corresponding stages upstream from the selected control point can be estimated from
water-surface profile computations.
(2) Computer Program. It is generaJiy impossible to make all of the flood routings
necessary to evaluate the effect of a reservoir system by hand computations. Computer
programs have been developed to route floods through a reservoir system with complex
operational criteria (55).
a. General Effects. Urbanization has tw9 major effects on the watershed which
influence the runoff characteristics. First, there is a substantial increase in the impervious
area, which results in more water entering the stream system as direct runoff. Second, the
drainage system coJiecting the runoff is generally more efficient and tends to concentrate
the water faster in the downstream portion of the channel system. It is important to keep
these two effects in mind when considering the changes in the flood peak frequency curve
caused by increasing urbanization.
3-28
EM 1110-2-1415
5 Mar 93
~ •u ~~
183
..
3
~/ "'
~
....u..0
~
Ill
• ' 4
/
Q.
II
c
~ """
.... . ~f.·' ..,
II
J
c: 182 , ~·~·
c:
<I
99.99 99.9 99 90 50 10 1 .1 . 01
3-29
EM 1110-2-1415
5 Mar 93
3-30
EM 1110-2-1415
5 Mar 93
(a) Develop peak-discharge frequency curve for specific land use conditions
from available gaged data and/or regional relationships.
(b) Develop balanced storms of various frequencies using data from generalized
criteria, a nearby gage or the equivalent.
(c) Develop rainfall-runoff model for the specific watershed with the adopted
land-use conditions. Calibrate runoff and routing parameters by reproducing
observed hydrographs occurring under natural conditions.
(a) Simulate all major historic events with a relatively simple model to establish
the ranking of events and an approximate peak discharge for each. The
approximate peaks could be developed by using a multiple linear regression
approach, by using a very simple rainfall-runoff model, or by any other
approach that will capture the hydrologic response of the basin.
(d) Use the relationship developed in step c to determine the desired frequency
curve. The same approach can be followed for both existing and future
conditions.
3-31
EM 1110-2-1415
5 Mar 93
CHAPTER 4
LOW- FLOW FREQUENCY ANALYSIS
4-1. ~. Low-flow frequency analyses are used to evaluate the ability of a stream to
meet specified flow requirements at a particular location. The analysis can provide an
indication of the adequacy of the natural flow to meet a given demand with a stated
probability of experiencing a shortage. Additional analyses can indicate the amount of
storage that would be required to meet a given demand, again with a stated probability of
being deficient. The design of hydroelectric power plants, determination of minimum
flow requirements for water quality and/or fish and wildlife, and design of water storage
projects can benefit from low-flow frequency analysis.
4-2. Interpretation.
a. Analytical frequency techniques are usually not applicable to low-flow data
because most theoretical frequency distributions cannot satisfactorily fit the recorded data.
It is recommended that graphical techniques be used and that known geologic and
hydrologic conditions be kept in mind when developing the relationships. As the low
values are the major interest, the data are arranged with the smallest value first. The
probability scale is usually labeled "percent chance nonexceedance."
b. Annual low flows are usually computed for several durations (in days) with the
flow rate expressed as the mean flow for the period. For example, the USGS WATSTORE
output provides the mean flow values for daily durations of I, 3, 7, 14, 30, 60, 90, 120
and 183 days. The default values for the HEC program STATS are the same with the
exception of using a 15-day duration instead of 14 (Table 4-1 ). Often a climatic year
from April 1 to March 31 is specified to provide a definite separation of the seasonal
low-flow periods. Figure 4-1 is a plot of the data in Table 4-1.
4-3. Aoolication Problems.
a. Basin Develooment. The effects of any basin developments on low flows are
usually quite significant. For example, a relatively moderate diversion can be neglected
when evaluating flood flow relations, but it would reduce, or even eliminate, low flows.
Accordingly, one of the most important aspects of low flows concerns the evaluation of
past and future effects of basin developments.
b. Multi- Year Events. In regions of water scarcity and where a high degree of
development has been attained, projects that entail carryover of water for several years
are often planned. In such projects it is desirable to analyze low-flow volume frequencies
for periods ranging from 1-1/2 to 8-1/2 years or more. Because the number of
independent low-flow periods of these lengths, in even the longest historical records, is
very small and because the concept of multi- annual periods is somewhat inconsistent with
the basic concept of an "annual event;" there is no truly satisfactory way for computing
the percent chance nonexceedance for low-flow periods that are more than 1 year in
length. One procedure described in reference (37) has been used with long sequences of
synthetically generated streamflows to derive estimates of drought frequency. Although
4-1
EM 1110-2 1415
5 Mar 93
lfot.a - Data baaed on Climatic Year of April 1 of &ivan year throuah March 31 of next year.
the results obtained through the use of this procedure seem reasonable, it is impossible to
verify the accuracy of the frequency estimates.
c. Regionalization. Regionalization of low-flow events is usually not very
successful. The variations in geologic conditions such as depth to ground water, size of
ground water basin, permeability of the aquifer, etc., are not easily quantifiable to enable
translation into probable low-flow rates. It may be possible to estimate low-flow rates on
a per unit area basis for a given exceedance frequency if the study area is relatively
homogeneous with respect to geology, topography, and climate. If information is needed
at several ungaged sites, the procedures described by Riggs (28) should be reviewed for
applicability.
4-2
EM 1110-2-1415
5 Mar 93
z
0 102~~----~----~--~~~~~~~--~~---------r--------~----4-----~
....c
1-4
a:
::1
0
a:
0
Y..
z 113-DAY
0 6
...J 10~~----~----~--------~--------_,~~~~~~~~~120-DAY--~----~
Y.. 0
z 110-DAY
c 110-DAY
w
% 30-DAY
1!5-DAY
7-DA'I'
3-DA'I'
0
s-ou
11!1.11!1 111.11 II 10 50 tO •t . 01
PERCENT CHANCE NONEXCEEDANCE
4-3
EM 1110-2-1415
5 Mar 93
CHAPTER 5
PRECIPITATION FREQUENCY ANALYSIS
•• ~
.r:.
u
c: /
/
·-I
..
c:
~
0
...,
•-4
...,• 181
•-4
~
~
•-4
u
•
I.
c
/
-4
c
•J
¢
c:
99.99 99.9 99 1 .1 • 91
5-1
EM 1110-2-1415
5 Mar 93
5-2. Available Regional Information. Where practical, use should be made of previous
precipitation-frequency-duration studies that have incorporated regional information. For
durations of 5 to 60 minutes in an area generally east of 105th meridian, see Hydro-35 (9).
For durations of 2 to 24 hours in the same area see Technical Paper 40 (10). Because of
the orographic effect, individual reports have been prepared for each of the 11 western
states (24). These reports have maps for 6- and 24-hour durations with extrapolation
procedures to obtain durations Jess than 6 hours. Longer duration events (2- to 10-days)
are presented in references (21), (22) and (23).
5-3. Derivation of Flood-Freauency Relations from Precipitation.
a. Application. Precipitation-frequency relations are often used to derive
flood-frequency relations where inadequate flow data are available or where existing (or
proposed) watershed changes have modified (or will modify) the rainfall-runoff
relationships. Guidelines for developing runoff frequencies from precipitation
frequencies are presented in references (10) and (44). Flood-frequency curves developed
by rainfall-runoff procedures often have tess variance (lower standard deviation) than
those developed from annual flood peaks. This results because not all the possible loss
rates for a given magnitude of precipitation are modeled. If extensive use wiJI be made of
frequency curves derived by rainfall-runoff modeling, an appropriate ratio adjustment for
the standard deviation should be developed for the region.
b. Calibration. Reference (44) describes the procedures involved in calibrating a
HEC-1 model to a flow-frequency curve based either on gaged data from a portion of the
basin or on regional flood-frequency relations. The coefficients from the calibrated
model must be consistent with those from nearby basins that have also been modeled. It
must be remembered that a frequency curve computed from observed flood peaks is based
on a relatively small sample. It is possible that the flow-frequency curve derived from
precipitation-frequency data is more representative of the population flow-frequency
curve than the one computed from the statistics of the observed flood peaks. But, there
are also errors in calibrating the model and establishing loss rates approximate with the
different frequency events. Therefore, the derivation of frequency relations by
rainfall-runoff modeling requires careful checking for consistency at every step.
c. Partial ouration. The precipitation-frequency relations presented in the National
Weather Service publications represent all the events above a given magnitude; therefore,
these relations are from a partial-duration series. The resulting flood frequency relations
must be adjusted if an annual peak flood frequency relationship is desired. Or, more
typically, the partial-duration series precipitation estimates are adjusted to represent
annual series estimates prior to use.
5-2
EM 1110-2-1415
5 Mar 93
CHAPTER 6
STAGE (ELEVATION) - FREQUENCY ANALYSIS
6-1. l.llli. Maximum stage-frequency relations are often required to evaluate inundation
damage. Inundation can result from a flooding river, storm surges along a lake or ocean,
wind driven waves (runup), a filling reservoir, or combinations of any of these. Minimum
elevation-frequency curves are used to evaluate the recreation benefits at a lake or
reservoir, to locate a water supply intake, to evaluate minimum depths available for
navigation purposes, etc. (Stages are referenced to an arbitrary datum; whereas, elevations
are generally referenced to mean sea level.)
6-2. Stage Data.
a. The USGS WATSTORE Peak Flow File has, in addition to annual peak flows,
maximum annual stages at most sites. Also, some sites located near estuaries have only
stage information because the flow is affected by varying backwater conditions.
b. River stages can be very sensitive to changes in the river channel and floodway.
Therefore, the construction of levees, bridges, or channel modifications can result in stage
data that is non-homogeneous with respect to time. For riverine situations, it is usually
recommended that the flow-frequency curve (Figure 6-1) and a rating curve (stage versus
1···-..---~----P------.--------P-------~------~--~----~
t
i .· .......·/
v
.....
]~
. . . . ·. . . vv
11 6'-ti---~----t------+--------~~----~-,~--+---_,~--ti
6-1
EM 1110-2-1415
5 Mar 93
6-2
EM 1110-2-1415
5 Mar 93
~
~
... i!
•• I""
IL
•
~
Cl
...Ill• i!
/
1
I 0 Obeerved Peak Flow• and Stare•, 19-'9-85
•
I -
seeel
Ratinr Curve for Gace Location, Nov 12, 19-'8 to Pruent
111111
I I I
a••••• 2151111
I
311111
.. · /
3
.. v
...
, ~
...
.
....
: i!
I
IL
•
.
Cl
•
.... 2oe
...
/
Ill
~
...E . ..
•" 2
.....
..
/
~
.
.
.. .
2
.....
... v
Derived Stare-Frequenc:y Curve, Unrerulated Condition•
- - Derived Stare-Frequenc:y Curve, Rerulated Condition•
1
I I I I I I
SUI. SISI SUI. II IISI Sll 151 11 1 .1 .11
P•~•nt Ch.nce Exceea.nce
6-3
EM 1110-2-1415
5 Mar 93
IDS0
-- '
I' I 1 I
---
TOP OF INDUCED
'
I 1047. s 1- SWCHARCE
-
;
!
!
oys 104/o. s 1- TOP OF FLOOD
I CONTllOL POOL
IU 0
--
z"
-- ! ~f. I
I
......c
0
1030
i i
- l
Cl
I
> i
IU
..J I
IU I I
1020
1 '
IDS ~
~ I
1010.0 ~ TOP OF IIATER
I~ SUPPLY POOL
---
1000
..... •••• " 10 so
PERCENT CHANCE EXCEEDANCE
ID I .I .DI
971.0 C....•TOP OF liiACTIV!
6-4
EM 1110-2-1415
S Mar 93
CHAPTER 7
DAMAGE-FREQUENCY RELATIONSHIPS
7-1. Introduction.
a. There are three methods that may be used to compute average annual damage and
are herein termed the historic method, the simulation method and the frequency method.
If 50 years of damage information were available for an area that has remained in
essentially the same land use with a reasonably constant level of economic activity,
historic damage could be scaled to the present to account for price differences (inflation)
and the average simply computed. This approach is termed the historic method and is the
most direct but is seldom used because sufficient data usually do not exist and the land use
and economic activity of an area are usually changing.
c. The most widely used approach within the Corps of Engineers is the frequency
technique. This technique is described in detail in Section 7-2. This technique addresses
the disadvantages of the previous two methods, and yet is fairly easily applied.
Experience in the development and application of damage functions is essential to
computation of reasonable estimates. Care should be taken to assure the rating curve is
not looped so that discharge is a unique function of stage. Otherwise more complex
functions that correctly relate stage and discharge should be developed and applied.
Damage functions in agricultural areas are often a function of the season and the duration
of flooding. Sensitivity analysis may be useful in determining the reliability of the
computed expected annual damage considering the uncertainties involved.
a. Figure 7-1 shows a schematic of the application of the three basic damage
evaluation functions used to compute the expected value of the annual damage. The term
"expected" is used rather than "average" because a frequency curve is used to represent the
distribution of future flood events and the expected value of damage is computed by the
summation of probability weighted estimates of damage.
b. The steps involved in determining the reduction in annual damage due to project
measures are:
7-1
EM 1110-2-1415
5 Mar 93
Rating Cw"ve
Discharge. cfs
Scale
0.118 0.110
Damage,$
Exceedanee Probabllty
Clrve
0 Al'tltiiMtlc Scilla
Exceedance Probablty
7-2
EM 1110-2-1415
5 Mar 93
(2) Combine the stage-damage and the stage-discharge relations into an intermediate
discharge-damage function. Make certain that the stage datum for the stage-damage
and stage-discharge functions is consistent for the index location.
(3) Combine the discharge-exceedance frequency (in events per year) and
discharge-damage function into a damage-exceedance frequency relationship.
(4) Compute the area beneath the damage-exceedance frequency relation (expected
annual damage) for each index location and sum to obtain the total expected annual
flood damage.
(5) Repeat step (I) for each alternative flood plain management plan under
investigation, i.e., revise the three basic evaluation functions as necessary.
(6) Repeat steps (2)-(4).
(7) Subtract results of step (4) (with project) for each plan from results of step (4)
for without-project measures. The differences will be expected annual damage
reduction (raw damage reduction benefits) for each plan.
7-3
EM 1110-2-1415
5 Mar 93
CHAPTER 8
8-1. Objective. One principal advantage of analytical frequency analysis is that there are
means for evaluating the reliability of the parameter estimates. This permits a more
complete understanding of the frequency estimates and provides criteria for
decision-making. For instance, a common statistical index of reliability is the standard
error of estimate, which is defined as the root-mean-square error. In general, it is
considered that the standard error is exceeded on the positive side one time out of six
estimates, and equally frequently on the negative side, for a total of one time in three
estimates. An error twice as large as the standard error of estimate is considered to be
exceeded one time in 40 in either direction, for a total of one time in 20. These
statements are based on an assumed normal distribution of the errors; thus, they are only
approximate for other distributions of errors. Exact statements as to error probability
must be based on examination of the frequency curve of errors or the distribution of the
errors. Both the standard error of estimate and the confidence limits are discussed in this
chapter.
8-2. Reliability of Frequency Statistics. The standard errors of estimate of the mean,
standard deviation, and skew coefficient, which are the principal statistics used in
frequency analysis, are given by the following equations:
Ss • S/(2N)"' (8-2)
SG = {6N(N-1)/[(N-2)(N+l)(N+3)]}Vz (8-3)
where:
8-1
EM 1110-2-1415
5 Mar 93
were calculated. While the expected frequency is that shown by the middle curve, there is
one chance in 20 that the true value for any given frequency is greater than that indicated
by the .05 curve and one chance in 20 that it is smaller than the value indicated by the .95
curve. There are, therefore, nine chances in 10 that the true value lies between the .05
and .95 curves. Appendix E and Example 1 in Appendix 12 of Bulletin 17B (40) provide
additional information and example computations.
/
...
~
Confidence Limit Cu"•
•
~
u 18 4
/
J .·· ~-··
...
D
~
IL. .···
.Jt
••
~
ll ... ···
... .
•
:I
c 18 3
f ......
... ·
....,.. .
c ~-
G:
l.........·
8-2
EM 1110-2-1415
5 Mar 93
CHAPTER 9
REGRESSION ANALYSIS AND APPLICATION TO REGIONAL STUDIES
a. General. Regression analysis is the term applied to the analytical procedure for
deriving prediction equations for a variable (dependent) based on given values of one or
more other variables (independent). The dependent variable is the value sought and is to
be related to various explanatory variables which will be known in advance, and which
will be physically related to the dependent variable. For example, the volume of
spring-season runoff from a river basin (dependent variable) might be correlated with the
depth of snow cover in the watershed (explanatory variable). Recorded values of such
variables over a period of years might be graphed and the apparent relation sketched in by
eye. However, regression analysis will generally permit a more reliable determination of
the relation and has the additional advantage of providing a means for evaluating the
reliability of the relation or of estimates based on the relation.
b. Definitions. The function relating the variables is termed the "regression
equation," and the proportion of the variance of the dependent variable that is explained
by the regression equation is termed the "coefficient of determination," which is the
square of the "correlation coefficient." Correlation is a measure of the association between
two or more variables. Regression equations can be linear or curvilinear, but linear
regression suffices for most applications, and curvilinear regression is therefore not
discussed herein. Often a curvilinear relation can be linearized by using a logarithmic or
other transform of one or more of the variables.
9-2. Calculation of Regression Equations.
a. Simple Regression. In a simple regression (one in which there is only one
independent, or explanatory, variable), the linear regression equation is written:
Y = a+ bX (9-1)
in which Y is the dependent variable, X is the independent variable, "a" is the regression
constant, and "b" is the regression coefficient. The coefficient "b" is evaluated from the
tabulated data by use of the following equations:
(9-2a)
or
(9-2b)
in which y is the deviation of a single value y; from the mean ( Y) of its series, x is
similarly defined, Sy and Sx are the respective standard deviations and R is computed by
9-1
EM III0-2-14J5
5 Mar 93
Equation 9- J J. The regression constant is obtained from the tabulated data by use of the
following equation:
a "" Y- bX (9-3)
All summations required for a simple linear regression can be obtained using Equations 9-
8 and 9-9a.
(9-4)
In the case of two explanatory variables, the regression coefficients are evaluated from the
tabulated data by solution of the following simultaneous equations:
In the case of three explanatory variables, the b coefficients can be evaluated from the
tabulated data by solution of the following simultaneous equations:
For cases of more than three explanatory variables, the appropriate set of simultaneous
equations can be easily constructed after studying the patterns of the above two sets of
equations. In such cases, solution of the equations becomes tedious, and considerable time
can be saved by use of the Crout method outlined in reference (51) or (52). Also,
programs are available for solution of simple or multiple linear regression problems on
practically any type of electronic computer. For multiple regression equations, the
regression constant is determined as follows:
(9-7)
In Equations 9-2, 9-5 and 9-6, the quantities I<x) 2 , L(yx) and I<x 1x2 ) can be determined
by use of the following equations:
9-2
EM 1110-2-1415
5 Mar 93
(9-11)
c. Standard Error. The adjusted standard error (Se) of a set of estimates is the
root-mean-square error of those estimates corrected for the degrees of freedom. On the
9-3
EM I I 10-2-1415
5 Mar 93
average, about one out of three estimates will have errors greater than the standard error
and about one out of 20 will have errors greater than twice the standard error. The
adjusted error variance is the square of the adjusted standard error. The adjusted
standard error or error variance of estimates based on a regression equation is calculated
from the data used to derive the equation by use of one of the following equations:
L(Y) 2
- - b 1 L(yx 1)- b2
df
lJYX 2) ••• - b" L(Yxn)
(9-13a)
(9-13b)
b. Physical Relationship. The values in the table are the annual peak flows for the
water years 1965-1985 (21 values). These two stations are less than 20 miles apart and are
likely to be subject to the same storm events; therefore, the first requirement of a
9-4
EM 1110-2-1415
5 Mar 93
se • 0.15646
9-5
EM 1110-2-1415
5 Mar 93
10s~----------~T~W~O_-~ST~A~T~I~O~N~C~O~MP~A~R~I~S~I~O~N~E~XA~M~P~L~E~-----------.
en
u.
u
Curve 8
::.:::
c
IJJ
Q.
Curve A
...J
c
:::.
z
z
c
a:
UJ
....>a:
~
c(
...J
:::.
...J
...J
c(
~
9-6
EM 1110-2-1415
5 Mar 93
9-7
EM 1110-2-1415
5 Mar 93
d. A verase Slooe. If it is obvious that all of the pertinent variables are included in
the analysis, then the variance of the points about the regression line is due entirely to
measurement errors, and the resulting difference in slope of the regression lines is entirely
artificial. In cases where all pertinent variables are considered and most of the
measurement error is in one variable, that variable should be used as the dependent
variable. Its errors will then not affect the slope of the regression line. In other cases
where all pertinent variables are considered, an average slope should be used. An average
slope can be obtained by use of the following equation:
(9-14)
9-7. Partial Correlation. The value gained by using any single variable (such as April
precipitation) in a regression equation can be measured by making a second correlation
study using all of the variables of the regression equation except that one. The loss in
correlation by omitting that variable is expressed in terms of the partial correlation
coefficient. The square of the partial correlation coefficient is obtained as follows:
9-8
EM 1110-2-1415
5 Mar 93
INPUT DATA
STATISTICS OF DATA
STANDARD
VARIABLE AVERAGE VARIANCE DEVIATION
REGRESSION R~~YLT~
PARTIAL
INDEPENDENT REGRESSION DETERMINATION
VARIABLE COEFFICIENT COEFFICIENT
UNBIASED STANDARD
REGRESSION R R ERROR OF
CONSTANT SQUARE SQUARE ESTIMATE
•• 223698 .9437 .9226 .0375
9-9
EM 1110-2-1415
5 Mar 93
(9-15)
in which the subscript to the left of the decimal indicates the variable whose partial
correlation coefficient is being computed, and the subscripts on the right of the decimal
indicate the independent variables. An approximation of the partial correlation can
sometimes be made by use of beta coefficients. After the regression equation has been
calculated, beta coefficients are very easy to obtain by use of the following equation:
(9-16)
The beta coefficients of the variables are proportional to the influence of each variable on
the result. While the partial correlation coefficient measures the increase in correlation
that is obtained by addition of one more explanatory variable to the correlation study, the
beta coefficient is a measure of the proportional influence of a given explanatory variable
on the dependent variable. These two coefficients are related closely only when there is
no interdependence among the various explanatory variables. However, some explanatory
variables naturally correlate with each other, and when one is removed from the equation,
the other will take over some of its weight in the equation. For this reason, it must be
kept in mind that beta coefficients indicate partial correlation only approximately.
9-8. Verification of Regression Results. Acquisition of basic data after a regression
analysis has been completed will provide an opportunity for making a check of the results.
This is done simply by comparing the values of the dependent variable observed, with
corresponding values calculated from the regression equation. The differences are the
errors of estimate, and their root-mean-square is an estimate of the standard error of the
regression-equation estimates (Paragraph 9-3). This standard error can be compared to
that already established in Equation 9-13. If the difference is not significant, there is no
reason to suspect the regression equation of being invalid, but if the difference is large,
the regression equation and standard error should be recalculated using the additional data
acquired.
9-9. Regression by Graphical Techniaues. Where the relationships among variables used
in a regression analysis are expected to be curvilinear and a simple transformation cannot
be employed to make these relationships linear, graphical regression methods may prove
useful. A satisfactory graphical analysis, however, requires a relatively large number of
observations and tedious computations. The general theory employed is similar to that
discussed above for linear regression. Methods used will not be discussed herein, but can
be found in references 8 and 27.
9-10. Practical Guidelines. The most important thing to remember in making
correlation studies is that accidental correlations occur frequently, particularly when the
number of observations is small. For this reason, variables should be correlated only when
there is reason to believe that there is a physical relationship. It is helpful to make
preliminary examination of relationships between two or more variables by graphical
plotting. This is particularly helpful for determining whether a relationship is linear and
in selecting a transformation for converting curvilinear relationships to linear
relationships. It should also be remembered that the chance of accidentally high
9-10
EM 1110-2-1415
5 Mar 93
correlation increases with the number of correlations tried. If a variable being studied is
tested against a dozen other variables at random, there is a chance that one of these will
produce a good correlation, even though there may be no physical relation between the
two. In general, the results of correlation analyses should be examined to assure that the
derived relationship is reasonable. For example, if streamflow is correlated with
precipitation and drainage area size, and the regression equation relates streamflow to
some power of the drainage area greater than one, a maximum exponent value of one
should be used, because the flow per square mile usually does not increase with drainage
area when other factors remain constant.
- Elevation of station
- General slope of surrounding terrain
- Orientation of that slope
- Elevation of windward barrier
- Exposure of gage
- Distance of leeward controlling ridge
9-11
EM 1110-2-1415
5 Mar 93
Statistics based on runoff measurements might be correlated with the following factors:
Q • aD+ bA + c (9-17)
then it can be seen that one inch change in precipitation would add the same amount of
flow, regardless of the size of drainage area. This is not reasonable, but again a
transformation to logarithms would yield a reasonable relation:
or transformed:
(9-19)
Thus, if logarithms of certain variables are used, doubling one independent quantity will
multiply the dependent variable by a fixed ratio, regardless of what fixed values the other
independent variables have. This particular relationship is reasonable and can be easily
visualized after a little study. There is no simple rule for deciding when to use
9-12
EM 1110-2-1415
5 Mar 93
logarithmic transformation. It is usually appropriate, however, when the variable has a
fixed lower limit of zero. The transformation should provide for near-uniform variance
throughout the range of data.
e. Example of Regional Correlation. An illustrative example of a regional
correlation analysis for the mean log of annual flood peaks (Y) with several basin
characteristics is shown on Figure 9-4. In this example, the dependent variable is
primarily related to the drainage area size, but precipitation and slope added a small
amount to the adjusted determination coefficient. The regression equation selected for the
regional analysis included only drainage area as an independent variable.
f. Selection of Useful Variables. In the regression equations shown on Figure 9-4,
the adjusted determination coefficient increases as variables are deleted according to their
lack of ability to contribute to the determination. This increase is because there is a
significant increase in the degrees of freedom as each variable is deleted for this small
sample of 20 observations. Both the adjusted determination coefficient and standard error
of estimate should be reviewed to determine how many variables are included in the
adopted regression equation. Even in the case of a slight increase in correlation obtained
by adding a variable, consideration of the increased unreliability of R as discussed in
Section 9-3 might indicate that the factor should be eliminated in cases of small samples.
The simplest equation that provides an adequate predictive capability should be selected.
In this example,· there is some loss in determination in only using drainage area, but this
simple equation is adopted to illustrate regional analysis. The adopted equation is:
or
Y • 38.5 AREA · 962 (9-21)
9-13
EM 1110-2-1415
5 Mar 93
JII~T RATA
015110 OIS ID AREA SLOPE LEIIGTN LAICES ELEY FC.EST PIECIP SOILS MEAN
1 5090 292.0 6.3 31.5 1.5 1.230 35.0 40.0 3.5 3.783
2 5140 185.0 14.3 30.3 1.0 1.024 52.0 38.2 3.2 3.783
3 5180 282.0 44.0 29.8 1.0 1.740 64.0 35.0 3.3 4.030
4 5200 298.0 zo. 1 37.3 1.0 1.600 30.0 36.5 3.3 4.044
5 5205 171.0 24.4 44.8 1.0 1.447 33.0 36.0 3.2 4.333
6 5260 114.0 35.8 17.5 1.0 1.3&3 30.0 34.0 3.0 3.751
7 5270 52.2 29.9 17.4 1.0 1.489 21.0 33.0 3.0 2.637
I 5280 66.8 12.6 20.1 1.0 1.305 41.0 33.6 3.0 3.186
9 5305 17.5 45.2 18.9 1.0 1.123 27.0 34.6 3.0 3.348
10 5320 215.0 17.7 27.5 1.1 .966 57.0 35.5 3.0 3.99S
11 5340 383.0 21.3 36.7 3.5 1.370 54.0 42.0 3.3 4.122
12 5375 15.7 291.0 5.6 1.0 1.350 11.0 40.0 3.5 2.722
13 5380 43.8 52.2 22.1 1.0 1.300 15.0 43.0 3.5 3.078
14 5390 274.0 39.6 32.3 1.6 1.200 70.0 43.0 2.8 3.930
15 5445 136.0 37.4 22.7 1.0 1.100 94.0 43.0 4.9 3.590
16 5415 604.0 22.1 44.5 1.0 1.900 13.0 37.0 3.2 4.092
17 5495 37.7 54.8 15.2 1.0 1.350 67.0 37.8 3.2 3.284
11 5500 173.0 45.7 24.2 1.0 1.700 65.0 39.0 3.9 3.116
19 5520 443.0 24.4 56.0 1.3 1.600 89.0 44.0 3.2 4.275
20 5525 23.1 115.7 10.0 1.0 1.100 80.o 47.0 4.3 3.249
STATJ~TI'~ QF ~ATA
STAIIDAIIO
VARIABLE AVEIIAGE YARIAIICE DEVIATIOII
VARIABLE AREA SLOPE LEIIGTN I.MES ELEV FOREST PIECJP SOILS lEAN
AREA 1.0000 ·.6327 .9304 .2749 .0000 .0000 .0000 .0000 .9159
SLOPE ·.6327 1.0000 ·.7144 ·.1318 .1147 .3635 .1867 .2053 ·.4521
LENGTH .9304 ·.7144 1.0000 .2345 .0000 .0000 .0000 •• 1096 .8263
STORAGE .2749 ·.1311 .2345 1.0000 .0000 .0000 .2112 .0000 .2596
ELEV .DODO .1117 .0000 .0000 1.0000 .27'91 .0000 .5297 .0000
FatEST .0000 .3635 .0000 .0000 . .27'91 1.0000 .6117 .4304 .0000
PIECIP • 0000 .1867 .0000 .2112 .0000 .6117 1.0000 .5412 .0000
SOILS .0000 .2053 ·.1096 .0000 .5297 .4304 .5412 1.0000 .0000
NEAll .9159 ·.4521 .8263 .2596 .0000 .0000 .0000 .0000 1.0000
9-14
EM 1110-2-1415
5 Mar 93
R!XHOW. ANALYSIS Win! RmU:SSictl ct1 0RAIN1GE M1.F.A ctn.Y
9-15
EM 1110-2-1415
5 Mar 93
(8) Eliminate variables in turn that contribute the least to the determination
coefficient, recomputing the determination coefficient each time, and select the
regression equation having the highest adjusted determination coefficient, or one
with fewer variables if the adjusted determination coefficient is nearly the same.
(9) Compute the regression errors for each station, plot on a suitable map, and draw
isopleths of the regression errors for the regression equations of X (see Figures 9-5
and 9-6 for an example) and S considering the standard error for each computed, or
adjusted, X and S. Note that an alternate procedure is to add the regression constant
to each error value and develop a map of this combined value. This procedure
eliminates the need to keep the regression constant in the regression equation as the
mapped value now includes the regression constant.
(10) A frequency curve can be computed for any ungaged basin in the area covered
within the mapped region by using the adopted regression equations and appropriate
map values to obtain X and S, and then using the procedures discussed in Section 3-2
to compute several points to define the frequency curve. (It may also be necessary to
develop regional (generalized) values of the skew coefficient if the Pearson type III
distribution is considered appropriate. The next section describes the necessary steps
to compute a generalized skew coefficient.)
i. Generalized Skew Determinations. Skew coefficients for use in hydrologic studies
should be based on regional studies. Values based on individual records are highly
unreliable. Figure 9-7 is a plot of skew coefficients sequentially recomputed after adding
the annual peak for the given year. Note that, after 1950, the skew coefficient was at a
at a minimum of about 0.5 in 1954 and maximum of about 1.9 in 1955, only one year
9-16
EM 1110-2-1415
5 Mar 93
apart. The procedures for developing generalized skew values are generally set forth in
Bulletin 17B (pages I 0-15).
In summary, it is recommended that:
(I) the stations used in the study have 25 or more years of data,
(2) at least 40 stations be used in the analysis, or at least all stations surrounding the
area within 100 miles should be included,
(3) the skew values should be plotted at the centroid of the basins to determine if
any geographic or topographic trends are present,
(6) then select the method that provides the most accurate estimation of the skew
coefficient (smallest mean-square error).
In addition to the above guidelines, care should be taken to select stations without
significant man-made changes such as reservoirs, urbanization, etc.
9-17
EM 1110-2-1415
5 Mar 93
a.h Kl at Shoenw~Wra. PA
Year
IUh Kl at ShoerniiU'a. PA
2.1
2.0
1.1
r'
~~
1.1
1.7
l
., I. I
II 1\
&-
...s 1.11 ~
0 1.4 ~
1:
.2 l.ll
1\ 1\
\ ~
~ 1.2
8 1.1 \ .A
0 v
lit
CD 1.0
II\ \
ll5 0.1 1
0.1
lA
0.7 \...
0.1
v ...t"""\
0.1
0.4
1106 1111 ,... ""
Year
..
, 11171 111116
9-18
EM 1110-2-1415
5 Mar 93
CHAPTER 10
ANALYSIS OF MIXED POPULATIONS
I 0-1. Definition. The term mixed population, in a hydrologic context, is applied to data
that results from two or more different, but independent, causative conditions. For
example, floods originating in a mountainous area or the northern part of the United
States at a given site could be caused by melting snow or by rain storms. Along the Gulf
and Atlantic coasts, floods can be caused by general cyclonic storms or by intense tropical
storms. A frequency curve representing the events caused by one of the climatic
conditions may have a significantly different slope (standard deviation) than for the other
condition. A frequency plot of the annual events, irrespective of cause, may show a
rather sudden change in slope and the computed skew coefficient may be comparatively
high. In these situations, a frequency curve derived by combining the frequency curves
of each population can result in a computed frequency relation more representative of the
observed events.
10-2. Procedure.
a. The largest annual event is selected for each causative condition. As Bulletin 17B
(46) cautions, "If the flood events that are believed to comprise two or more populations
cannot be identified and separated by an objective and hydrologically meaningful
criterion, the record shall be treated as coming from one population." Also, Bulletin 17B
states, "Separation by calendar periods in lieu of separation by events is not considered
hydrologically reasonable unless the events in the separated periods are clearly caused by
different hydrometeorologic conditions."
b. The frequency relations for each separate population can be derived by the
graphical or analytical techniques described in Chapter 2·and then combined to yield the
mixed population frequency curve. The individual annual frequency curves are combined
by "probability of union." For two curves, the equation is:
(10-1)
where:
10-1
EM 1110-2-1415
5 Mar 93
(10-2)
NOTE: The exceedance probability (percent chance exceedance divided by 100) must be
used in the above equations.
10-3. Cautions.
a. If annual flood peaks have been separated by causative factors, a generalized skew
must be derived for each separate series to apply the log-Pearson Type III distribution as
recommended by Bulletin 17B. Plate 1 of Bulletin 17B or any other generalized skew map
based on the maximum annual event, irrespective of cause, will not be applicable to any
of the separated series. Derivation of generalized skew relations for each series can
involve much effort.
b. Some series may not have an event each year. For example, tropical storms do not
occur every year over most drainage areas in the United States, and quite often there are
only a few flood events for the series. Extensive regionalization may be necessary to
reduce the probable error in the frequency relations which results from small sample sizes.
10-2
EM 1110-2-1415
5 Mar 93
.
D Annual Nonhurricane Event1, 19%9-85
~~
""•• ~
..•
II.
:I ~
__,p-
lt. "'
c 114
c
« ~ I
lc
I
I
u:J
88.811 88.8 88 1 .1 .11
•
~
0 11 15
/
/
iD
....
""••
II.
... ~
•
:I
11 4
~
«
E ~a
Ia
10-3
EM III0-2-I4I5
5 Mar 93
CHAPTER II
FREQUENCY OF COINCIDENT FLOWS
b. Selection of pominant Varjable. The variable that has the largest influence on
variable C is designated as variable A; the less influential variable is designated as variable
B. The significance of "influential" will be indicated by means of an example. Figure
11-1 shows water surface profiles along a tributary near the junction with a main river.
Stage on the tributary (variable C) is a function of main river stage and tributary
discharge. In Region I, main river stage, will tend to have the dominant influence on
tributary stage, whereas in Region II, tributary discharge will tend to dominate. The
boundary between Regions I and II cannot be precisely defined and will vary with
exceedance frequency. Stage-frequency determinations will be least accurate in the
vicinity of the boundary where both variables have a substantial impact on the combined
result.
c. Procedure.
(1) Construct a duration curve for variable B. Discretize the duration curve with a
set of "index" values of B. Index values should represent approximately equal ranges
of magnitude of variable B. The area under the resulting discretized duration curve
should equal the area under the original duration curve. The number of index values
of B required for discretization depends on the range of variation of B and the
sensitivity of variable C to B. Therefore, the number of points selected should
adequately define the relationships.
(2) For each of the index values of variable B, develop a relationship between
variable A and the combined result C. In the illustration (Figure 11-1) the
relationship linking variables A, B and C would be obtained with a set of water
surface profile calculations for various combinations of main river stage and
tributary discharge.
11-1
EM II 10-2-1415
5 Mar 93
IUIOI I
IOIIIAL IIAJI
IJYEI STAGE
11-2
EM 1110-2-1415
5 Mar 93
CHAPTER 12
STOCHASTIC HYDROLOGY
12-1. Introduction.
a. A stochastic process is one in which there is a chance component in each
successive event and ordinarily some degree of correlation between successive events.
Modeling of a stochastic process involves the use of the "Monte Carlo" method of adding a
random (chance) component to a correlated component in order to construct each new
event. The correlated component can be related, not only to preceding events of the same
series, but also to concurrent and preceding events of series of related phenomena.
b. Work in stochastic hydrology has related primarily to annual and monthly
streamflows, but the results often apply to other hydrologic quantities such as
precipitation and temperatures. Some work on daily streamflow simulation has been done.
12-2. Apolications.
a. Hydrologic records are usually shorter than 100 years in length, and most of them
are shorter than 25 years. Even in the case of the longest records, the most extreme
drought or flood event can be far different from the next most extreme event. There is
often serious question as to whether the extreme event is representative of the period of
record. The severity of a long drought can be changed drastically by adding or
subtracting I year of its duration. In order that some estimate of the likelihood of more
severe sequences can be made, the stochastic process can be simulated, and long sequences
of events can be generated. If the generation is done correctly, the hypothetical sequence
would have as equal likelihood of occurrence in the future as did the observed record.
b. The design of water resource projects is commonly based on assumed recurrence
of past hydrologic events. By generating a number of hydrologic sequences, each of a
specified desired length, it is possible to create a much broader base for hydrologic design.
While it is not possible to create information that is not already in the record, it is possible
to use the information more systematically and more effectively. In selecting the number
and length of hydrologic sequences to be generated, it is usually considered that 10 to 20
sequences would be adequate and that their length should correspond to the period of
project amortization.
c. It must be recognized that the more hydrologic events that are generated, the
more chance there is that an extreme event or combination of events will be exceeded.
Consequently, it is not logical that a design be based on the most extreme generated event,
but rather on some consideration of the total consequences that would prevail for a given
design if all generated events should occur. The more events that are generated, the less
proportional weight each event is given. If a design is tested on 10 sequences of
hydrologic events, for example, the benefits and costs associated with each sequence
would be divided by 10 and added in order to obtain the "expected" net benefits.
12-3. Basic Procedure. Successful simulation of stochastic processes in hydrology has
been based generally on the concept of multiple linear regression, where the regression
12-1
EM 1110-2-1415
5 Mar 93
equation determines the correlated component, and the standard error of estimate
determines the random component. Figure 12-1 illustrates the general nature of the
process. In this case, a low degree of correlation is illustrated, in order to emphasize
important aspects of the process. It can be seen that, if every estimate of the dependent
variable is determined by the regression line (Figure 12-1 a), the estimated points would be
perfectly correlated with the independent variable and would have a much smaller range
of magnitude than the actual observed values of the dependent variable. In order to avoid
such unreasonable results, it is necessary to add a random component to each estimate
(Figure 12-1 b), and this random component should conform to the scatter of the observed
data about the regression line.
6.0
n
__,...
4.0
0 /p
~--""
0
J,/~ 0
~~
0
1.0
v
0.0
0.0 1.0 2.0 3.0 4.0 e.o
lndlpelldent variable
4.0
FIMdarnEmlr
() . ~
__,...
i ~·.
~
.. 0
3.0
-----
I .. I 0
• vv i--""'
• 0
~
~ 0
1.0
~ •
0.0
0.0 1.0 2.0 3.0 4.0 6.0
kldapelldent variable
12-2
EM 1110-2-1415
5 Mar 93
12-4. Monthly Streamflow Model.
a. In accordance with the above basic procedure, a simulation model for generating
values of a variable which can be defined only partially by a deterministic relation is:
( 12-1)
where:
Y = dependent variable
a = regression constant
b 1,b 2 = regression coefficients
b. This type of simulation model can be used to generate related monthly streamflow
values at one or more stations. Multiple linear regression theory is based on the assumed
distribution of all variables in accordance with the Gaussian normal distribution.
Therefore, mathematical integrity requires that each variable be transformed to a normal
distribution, if it is not already normal. It has been found that the logarithms of
streamflows are approximately normally distributed in most cases. For computational
efficiency it is convenient to work with deviations from the mean which have been
normalized by dividing by the standard deviation. This deviate is sometimes called the
Pearson Type Ill deviate and can be computed as follows:
(12-2)
where:
12-3
EM 1110-2-1415
5 Mar 93
(12-3)
where:
• month number
G = skew coefficient
= Pearson Type III deviate as defined in Equation 12-2
(12-4)
where:
12-4
EM 1110-2-1415
5 Mar 93
(12-5)
d. Note that Equation 12-5 is very similar to Equation 12-1. The differences result
from using normal standard deviates. When this is done, the regression constant, a, equals
zero, the regression coefficients, b, become beta coefficients, /3, and the standard
deviation, S, does not appear in the random component since it equals I. Note also that
one of the independent variables is the flow for the preceding month in order to preserve
the inherent serial correlation. The flow value in the original units is computed by
reversing the transformation process, i.e., from normal standard deviate to Pearson Type
III deviate, to logarithm of flow and finally flow value.
e. A step-by-step procedure for generating monthly streamflows for a number of
interrelated locations having simultaneous records is as follows:
(1) Compute the logarithm of each streamflow quantity. If a value of zero
streamflow is possible, it is necessary to add a small increment, such as 0.1
percent of the mean annual flow, to each monthly quantity before taking the
logarithm.
(2) Compute the mean, standard deviation and skew coefficient of the values
for each location and each month, using equations given in Chapter 2.
(3} For each month and location, subtract the mean from each event and divide
by the standard deviation (Equation 12-2).
(4) Transform these "standardized" quantities to a normal distribution by use of
Equation 12-3.
(5) Arrange the locations in any sequence, and compute a regression equation
for each location in turn for each month. In each case, the independent
variables will consist of concurrent monthly values at preceding stations and
preceding monthly values at the current and subsequent stations.
(6) Generate standardized variates for each location in turn for each month,
starting with the earliest month of generated data. This is accomplished by
computing a regression value and adding a random component. The random
component, according to Equation 12-5, is a random selection from a normal
distribution with zero mean and unit standard deviation, multiplied by the
alienation coefficient which is (I - R 2 )'1r.
(7) Transform each generated value by reversing the transform of Equation 12-
3 with the appropriate skew coefficient, multiplying by the standard deviation
and adding to mean in order to obtain the logarithm of streamflow.
(8) Find the antilogarithm of the value determined in step (7) and subtract the
small increment added in step ( 1). If a negative value results, set it to zero.
12-5
EM 1110-2-1415
5 Mar 93
f. It is obviously not feasible to accomplish the above computations without the use
of an electronic computer. A computer program, HEC-4 Monthly Streamflow Simulation
(51) can be used for this purpose.
12-5. Data Fill ln. Ordinarily, periods of recorded data at different locations do not
cover the same time span, and therefore, it is necessary to estimate missing values in order
to obtain a complete set of data for analysis as described above. In estimating the missing
values, it is important to preserve all statistical characteristics of the data, including
frequency and correlation characteristics. To preserve these characteristics, it is necessary
to estimate each individual value on the basis of multiple correlation with the preceding
value at that location and with the concurrent or preceding values in all other locations. A
random component is also required, as indicated in Equation 12-1.
12-6. Application In Areas of Limited Data. The streamflow generation models discussed
so far have assumed that sufficient records were available to derive the appropriate
statistics. For instance, the monthly streamflow model requires four frequency and
correlation coefficients for each of the 12 months, or 48 values for one station simulation.
A model has been developed (51) that combines the coefficients into a few generalized
coefficients for the purpose of generating monthly streamflow at ungaged locations.
(Procedures for determining generalized statistics for use in generating daily flows have
not yet been developed.) The generalized model considers the following:
- season of maximum runoff
- lag to season of minimum runoff
- average runoff
- variation between maximum and minimum runoff
- standard deviation of flows
- interstation and serial correlations of flows
12-6
EM 1110-2-1415
5 Mar 93
precipitation over the basin, and (2) using a precipitation- runoff model to derive the
resulting streamflow.
12-8. Reliability. While the simulation of stochastic processes can add reliability in
hydrologic design, the techniques have not yet developed to the stage that they are
completely dependable. All mathematical models are simplified representations of the
physical phenomena. In most applications, simplifying assumptions do not cause serious
discrepancies. It is important at this "state of the art," however, to examine carefully the
results of hydrologic simulation to assure that they are reasonable in each case.
12-7
EM 1110-2-1415
5 Mar 93
APPENDIX A
SELECTED BIBLIOGRAPHY
7. Draper, N. R. and H. Smith, Applied Regression Analysis, John Wiley & Sons, Inc.,
1966.
8. Ezekiel, M. and Karl A. Fox, Methods of Correlation and Regression Analysis, 3rd
ed., John Wiley and Sons, 1959.
9. Frederick, R. H., et. al., "Five-To Sixty-Minute Precipitation Frequency for the
Eastern and Central U.S.," NOAA Technical Memorandum NWS Hydro- 35, Silver
Spring, MD, June 1977.
11. Haan, C. T., Statistical Methods jn Hydrology, The Iowa State University Press,
Ames, IA, 1977.
12. Hardison, C. H., "Storage Analysis for Water Supply," Chapter B2 of Hydrologic
Analysis and Interpretation, Techniques of Water-Resources Investigations, U.S.
Geological Survey, GPO, Washington, D.C., 1973.
13. Hazen, A., "Storage to be Provided in Impounding Reservoirs for Municipal Water
Supply", Trans. ASCE 77, p 1539-1559, 1914.
14. Hershfield, D. M., "Rainfall Frequency Atlas of the U.S.," Technical Paper No. 40,
Weather Bureau, U.S. Department of Commerce, Washington, D.C., May 1961.
A-I
EM 1110-2-1415
5 Mar 93
A-2
EM 1110-2-1415
5 Mar 93
31. Searcy, J. K., and C. H. Hardison, "Double-Mass Curves," U.S. Geological Survey
Water-Supply Paper 1541-B, GPO, Washington, D.C., 1960.
32. Searcy J. K., "Flow-Duration Curves," U.S. Geological Survey, Water-Supply Paper
1541-A, GPO, Washington, D.C., 1960.
33. Sokolov, A. A., et al., Flood Flow Computation. Methods Compiled from World
Experience, The UNESCO Press, Paris, France, 1976.
34. U.S. Army Corps of Engineers, "Stream Flow Volume-Duration-Frequency Studies,"
CW-152 Technical Report No. I, Washington District, Washington, D.C., June 1955.
35. U.S. Army Corps of Engineers," Frequency of New England Floods," CW-151
Research Note No. I, Sacramento District, Sacramento, CA, July 1958.
40. U.S. Army Corps of Engineers, "Regional Frequency Study, Upper Delaware and
Hudson River Basins, New York District," Hydrologic Engineering Center, Davis,
CA, November 1974.
A-3
EM 1110-2-1415
5 Mar 93
43. U.S. Army Corps of Engineers, "Generalized Skew for State of New Jersey," Special
Project Memo 480, Hydrologic Engineering Center, Davis, CA, December 1977.
44. U.S. Army Corps of Engineers, "Hydrologic Analysis of Ungaged Watersheds with
HEC-1 ," Hydrologic Engineering Center, Davis, CA, April 1982.
45. U.S. Army Corps of Engineers, Engineering Manual EM 1110-2-1413, "Hydrologic
Analysis of Interior Areas," January 1987.
46. Water Resources Council, "Guidelines for Determining Flood Flow Frequency,"
Bulletin 17B, Hydrology Committee, Washington, D.C .• March 1982.
47. Water Resources Council, "A Uniform Technique for Determining Flood Flow
Frequencies," Bulletin 15, Hydrology Committee, Washington, D.C .• December 1967.
48. World Meteorological Organization, "International Glossary of Hydrology," First
Edition, WMO No. 385, Geneva, Switzerland, 1974.
49. World Meteorological Organization, "Guide to Hydrometeorological Practices," WMO
No. 168, Third Edition, Geneva, Switzerland, 1974.
B. COMPUTER PROGRAMS
50. U.S. Army Corps of Engineers, "Multiple Linear Regression," Computer Program
704-G l-L2020, Hydrologic Engineering Center, Davis, CA, September 1970.
51. U.S. Army Corps of Engineers, "HEC-4, Monthly Streamflow Simulation," Computer
Program 723-X6-L2340, Hydrologic Engineering Center, Davis, CA, February 1971.
52. U.S. Army Corps of Engineers, "Regional Frequency Computation," Computer
Program 723-X6-L2350, Hydrologic Engineering Center, Davis, CA, July 1972.
53. U.S. Army Corps of Engineers, "Flood Flow Frequency Analysis, "Computer Program
723-X6-L 7550, Hydrologic Engineering Center, Davis, CA, February, 1982.
54. U.S. Army Corps of Engineers, "Expected Annual Flood Damage Computation,"
Computer Program 761-X6-L7580, Hydrologic Engineering Center, Davis, CA, June
1977.
55. U.S. Army Corps of Engineers," HEC-5, Simulation of Flood Control and·
Conservation Systems." Computer program 723-X6-L2500, Hydrologic Engineering
Center, Davis, CA, April 1982
56. U.S. Army Corps of Engineers, "HEC-1, Flood Hydrograph Package," Computer
Program 723-X6-L2010, Hydrologic Engineering Center, Davis, CA, September
1981.
A-4
EM 1110-2-1415
5 Mar 93
57. U.S. Army Corps of Engineers, "Storage, Treatment, Overflow, Runoff Model
'STORM'," Computer Program 723-58-L 7520, Hydrologic Engineering Center,
Davis, CA, August 1977.
58. U.S. Army Corps of Engineers, "Statistical Analysis of Time Series Data 'ST ATS' ,"
under development, Hydrologic Engineering Center, Davis, CA, May 1987.
A-5
EM 1110-2-1415
5 Mar 93
APPENDIX B
GLOSSARY
These definitions have been collected from four major sources: (I) Guidelines for
Determining Flood Flow Frequency ( 17B), reference 46; (2) International Glossary of
Hydrology (WMO), reference 48; (3) General Introduction and Hydrologic Definitions
(USGS), reference 18; and (4) Mathematics Dictionary (MD), reference 16.
TERM DEFINITION
Annual Event The most extreme event, either maximum or minimum, in the
year. (17B)
Annual Series A general term for a set of any kind of data in which each item is
the maximum or minimum in a year. (17B)
Array A list of data in order of magnitude; in flood- frequency analysis
it is customary to list the largest value first, in a low-frequency
analysis, the smallest first. ( 17B)
Autocorrelation See "Serial Correlation."
Base Discharge Usually refers to the discharge above which independent
instantaneous peak flows are collected for a partial duration
frequency analysis.
Biased The expected value of a statistic obtained from random sampling
is not equal to the parameter or quantity being estimated. (MD)
Broken Record A systematic record which is divided into separate continuous
segments because of deliberate discontinuation of recording for
significant periods of time. (17B)
Chi-Square The distribution of sample variances drawn from a normal
Distribution distribution. Used to compute confidence intervals for the
population variance estimated from a sample.
Class Interval A convenient sized interval into which data may be grouped. The
upper and lower bounds of the class interval are called "class
limits." (MD)
B-1
EM 1110-2-1415
5 Mar 93
Double Mass Curve Plot of successive accumulated values of one variable against the
contemporaneous accumulated values of another variable. (WMO)
Duration Curve A cumulative frequency curve that shows the percent of time that
specific values are equalled or exceeded. (USGS)
B-2
EM 1110-2-1415
5 Mar 93
Geometric Mean The Nth root of the product of N values or the antilogarithm of
the mean logarithm of a set of values.
Graphical Frequency The development of a frequency curve by drawing a smooth
Analysis curve through plotted points while considering known constraints.
Plotting positions are computed based on the order number and
the total number of values represented, and then plotted on the
appropriate probability paper.
Historic Data Information about significant events before or after the period of
"systematic" data collection. (derived from 17B)
Incomplete Record A streamflow record in which some peak flows are missing
because they were too low or high to record or the gage was out
of operation for a short period because of flooding. ( 17B)
Level of Significance The probability of rejecting a hypothesis when it is in fact true,
At a "10-percent" level of significance the probability is 1/10.
(17B)
B-3
EM 1110-2-1415
5 Mar 93
Log-Pearson Type Ill Application of the Pearson Type III distribution to the logarithms
Distribution of the data.
Mass Curve Curve of an accumulative quantity versus time. (WMO)
Mean The expected value of a random variable, the first moment. The
arithmic mean (or average) of a sample is an estimate of the
population mean.
Mean Daily The mean of daily values in a specified period, or the mean of
values within one day. (derived from WMO)
Median The value at which one-half of ordered observations lie on either
side. If there is no middle value, the median is the average of the
two middle values.
Method of Moments A standard statistical computation for estimating the moment of a
distribution from the data of a sample. (17B)
Mixed Populations A sample whose events have come from two or more different
populations, i.e., data not homogenous.
Mode The most frequent value of a set of numbers. (MD)
Non-Central t A distribution that combines the probable error in the mean and
Distribution the standard deviation for samples from a normal distribution.
Used in the development of confidence limit curves about a
frequency curve computed from sample statistics.
Nonstationary Not stationary with respect to time. See "Stationary Process."
Normal Distribution A probability distribution that is symmetrical about the mean,
median, and mode (bell shaped). It is the most studied
distribution in statistics, even though most data are not exactly
normally distributed, because of its value in theoretical work and
because many other distributions can be transformed into the
normal. It is also known as Gaussian, the Laplacean, the
Gauss-Laplace, or the Laplace-Gauss distribution, or the Second
Law of Laplace. (17B)
Outlier Outliers (extreme events) are data points which depart from the
trend of the rest of the data. (17B)
Parameter A characteristic descriptor of the population, such as mean or
standard deviation. Parameter estimates are called statistics.
Parent Population See "Population".
Partial-Duration Series A list of flood peaks that exceed a chosen base stage or discharge,
regardless of the number of peaks occurring in a year. (USGS)
B-4
EM 1110-2-1415
5 Mar 93
Pearson Type III Family of asymmetrical, theoretical frequency distributions of
Distribution which the normal distribution is a special case.
Percent Chance The probability, expressed as a percentage, with which values
Exceedance exceed a specified magnitude.
Percent Chance The probability, expressed as a percentage, with which values
Non-Exceedance will not exceed a specified magnitude.
B-5
EM 1110-2-1415
5 Mar 93
Sampling Error The difference between a random sampling statistic and the
parameter of the population from which the random sample was
drawn. (MD)
Serial Correlation A measure of the interdependence between an observation at a
given time period and that of a preceeding time period. Also
called autocorrelation.
Skew Coefficient See "Coefficient of Skewness."
Stage The height of a water surface above an established datum plane.
(USGS)
Standard Deviation A measure of the dispersion or precision of a series of statistical
values such as precipitation or stream flow. It is the square root
of the sum of squares of the deviations from the arithmetic mean
divided by the number of values or events in the series. It is now
standard practice in statistics to divide by the number of values
minus one in order to get an unbiased estimate of the variance
from the sample data. (17B)
Standard Error Standard deviation of the sampling distribution of a statistical
parameter. (WMO)
Stationary Process All of the generating moments of the frequency distribution
remain fixed with respect to time.
Statistic An estimate of a population parameter obtained from a sample of
the population.
Stochastic Process Process in which both the probability and the sequence of
occurrence of the variables are taken into account. (WMO)
Student's t- A distribution used in evaluation of variables which involve
Distribution sample standard deviation rather than population standard
deviation. ( 17B)
Systematic Record Information collected by a systematic data collection program.
(derived from 17B)
t-Distribution See "Student's t-Distribution."
Test of Significance A test mode to learn the probability that a result is accidental or
that a result differs from another result. For all the many types
of tests, there are standard formulae and tables. ( 17B)
Transformation The change of numerical values of data to make later
computations easier, to linearize a plot or to normalize a skewed
distribution by making it more nearly a normal distribution. The
most common transformations are those changing ordinary
numerical values into their logarithms, square roots or cube roots;
many others are possible. (17B)
B-6
EM I I 10-2-1415
5 Mar 93
Watershed The divide separating one drainage basin from another and in the
past has been generally used to convey this meaning. However,
over the years, use of the term to signify drainage basin or
catchment area has come to predominate, although drainage basin
is preferred. (USGS)
Water Year Continuous twelve-month period selected in such a way that all
solid and liquid precipitation runs off during this period. Thus,
carryover is reduced to a minimum. (WMO) In U.S. Geological
Survey reports, it is the twelve-month period. October 1 through
September 30. The water year is designated by the calendar year
within which most of the twelve-months occur. (USGS)
B-7
EM 1110-2-1415
5 Mar 93
APPENDIX C
COMPUTATION PROCEDURE
FOR
EXTREME VALUE (GUMBEL) DISTRIBUTION
(5.1)
where X is the mean value, and s,. is the standard deviation of the variable
being studied. Value XTr denotes the magnitude of the event reached or exceeded
on an average once in Tr years. K is the frequency factor. I£ X is not normally
distributed, K depends on frequency and skewness coefficient. A commonly used
distribution of extreme values (annual series) is the double exponential distribu-
tion, which bas been widely applied by Gumbel (see Bibliography), and often
bears his name. In this method
YTr _
K== __ -YD
_
(5.2)
Sn
C-1
EM 1110-2-1415
5 Mar 93
where Y,, the rl'cluccd mf'an, and So, the rt•Juccd standard !leviation, are functions
only of samplt· sizf'; and YT,, the rcdut·l'd variate, is rl'latcd to r·cturn period by
TaLl!' 5.3 a-i''''s \'&lues of K computed by means of Eq. (5.:!) using Gumbel's
values for )·n, sn, and YT.r
Therc arc two basic methods for fitting data to the extrcmc value distribution.
One consists in computation of XTr by means of Eq. (5 .1), after a previous
computation of the values of Xand Sz (Table 5. 4 ). The other consists in plotting
data on suitable graph paper, known as extreme probability paper, and drawing
a line by inspection.
TABLE 5.3
Valu,.s of K based on Eq. (5. 2)
(continued)
C-2
EM l l 10-2-1415
5 Mar 93
TABLE 5. 3 (continued)
Return period (years)
n 2 5 10 25 50 100
C-3
EM 1110-2-1415
5 Mar 93
C-4
EM 1110-2-1415
5 Mar 93
Extreme probability paper has a linear ordinate for the variable being
5tudied. and the abscissa is a linear scale of the reduced variate [Eq. (5.3)]. For
convenience in plotting, the warped scale of Tr is also shown along tlw top of
Fig. 5. 9. Plotting positions are commonly determined by the formulae [16):
Tr = n +1 (5 .4!
m
or
n 0.4
Tr=--~
+ (5.51
m-0.3
where n is the number of years of record (the number of items in the annual
series) and m is the rank of the item on the aeries, m being 1 for the largest.
To illustrate the steps in numerical computation of the rainfall value for a
given return period, hypothetical values of ll serie:; of annual rainfall maxima are
given in the upper part of Table 5. 4. Computations are illustrated in the lower
part of the table for Tr of 10. Rainfall depths for return periods other than
10 years can be computed in a similar manner.
/+
10
~
I
.b:(j,
."69
I
A-66
_!;;_"I
+'BD
I I
20 f--
0
I I I I I I I I I I I
t.O ell 10 M M ft "·'
~IIOIAIIUTY (' 00~+11]
1111111111!1!!1 1 11' d II' 1 It!!' ''11'1111'1' !IIIII' dt "''II' rl l!!!ll"'i.Lw.u.u.1.l.w.J.J
-l.tl II l.tl l.tJ ~0 4.0 ~0 4.11
IIEIHICI. I'AtfiArC
Figure 5.9- Example of extreme probability plot using data of Table 5.4.
C-5
EM 1110-2-1415
5 Mar 93
TABLE ; •. 4
Computation of eztreme values
II+ 1 f'2
"'l"t>ar I' Ill 1'-P (P-I'f
m
1960 37 i 1.72 -9 81 1,369
1961 20 11 1.09 -26 676 400
1962 32 8 1.5 -14 196 1,024
1963 60 3 4.0 + 14 196 3,600
1964 25 9 1.33 -21 441 625
1965 52 4 3.0 +6 36 2,704
1966 46 6 2.0 0 0 2,116
1967 70 2 6.0 + 24 576 4,900
1968 92 1 12.0 + 46 2,116 8,464
196~ 48 5 2.4 +2 4 2,304
1970 24 10 1.2 -22 484 576
Total 506 4,806 28,082
- n 50() ,
P
Sx
= l:P/n = IT =
. .
by square or devaataons:
4u.O
vi (P)Z-
n-1
Pi: p -- ~ -- 21.92
10
C-6
EM 1110-2-1415
5 Mar 93
within which, with given confidence levels, one may expect to find the true
precipitation value xT.Values of t(a) for selected confidence levels are as
follows:
In most cues, values of s., the standard error of estimate, can be computed
by means of the formula
s,.
Se = P-rr·--
'1/n
(5.6)
In particular, for the Gumbel distribution the following relation [19] exists:
where K is the numerical value defined by Eq. (5. 2) and readily obtainable from
Table 5.3. However, for convenience in the use of Eq. (5.6) values of lh-./v'~
can be obtained directly from Table 5. 5. Thus, in determining the 80 per cent
~:onfidence interval for P 10 = 85.7 in Table 5.4, for example,
C-7
EM 1110-2-1415
5 Mar 93
The lower and upper limits of the confidence interval are therefore 85.7- 21.9
and 85.7 + 21.9, respectively, which means that there is an 80 per cent probability
that the true value or PJO lies between 63.8 and 107.6. Similarly, the lower and
upper limits of the same confidence interval for P 1 == 43.0 are 35.1 and 50.9,
respectively, the value of fh.!V~ being 0.2803 (Table 5. 5).
TABLE 5.5
Valuu of fh-./Vn for tue in Eq. (5.6)
Return period (years)
n 2 5 10 25 50 100
10 .2942 .5863 .8285 1.1472 1.3873 1.6273
11 .2803 .5522 .7733 1.0761 1.3007 1.5252
12 .2681 .5232 .7358 1.0161 1.2275 1.4389
13 .2574 .4982 .6992 .9645 1.1646 1.3648
14 .2479 .4763 .6673 .9196 1.1100 1.3005
15 .2393 .4569 .6392 .8801 1.0620 1.2439
16 .2316 .4397 .6142 .8450 1.0193 1.1937
17 .2246 .4242 .5918 .8136 .9811 1.1488
18 .2182 .4102 .5716 .7853 .9467 1.1083
19 .2123 .3974 .5532 .7596 .9155 1.0716
20 .2068 .3857 .5365 .7361 .8871 1.0382
21 .2018 .3750 .5211 .7146 .8610 1.0075
22 .1971 .3651 .5069 .6948 .8370 .9793
23 .1927 .3559 .4937 .6765 .8148 .9532
24 .1886 .3473 .4815 .6595 .7942 .9290
25 .1847 .3394 .4702 .6437 .7750 .9064
26 .1811 .3319 .4595 .6289 .7571 .8854
27 .1777 .3249 .4496 .6150 .7403 .8657
28 .1745 .3183 .4402 .6020 .7245 .8472
29 .1714 .3121 .4314 .5898 .7097 .8298
30 .1685 .3062 .4230 .5782 .6957 .8134
31 .1657 .3007 .4152 .5673 .6825 .7978
32 .1631 .2954 .4077 .5569 .6699 .7831
33 .1606 .2904 .4006 .5471 .6581 .7692
34 .1582 .2856 .3938 .5377 .6468 .7559
35 .1559 .2811 .3874 .5289 .6360 .7433
36 .1537 .2767 .3813 .5204 .6257 .7313
37 .1516 .2726 .3754 .5123 .6159 .7198
38 .1496 .2686 .3698 .5045 .6066 .7088
39 .1476 .2648 .3645 .4971 .5976 .6983
( conlinU«i)
C-8
EM 1110-2-1415
5 Mar 93
C-9
EM 1110-2-1415
5 Mar 93
TABLE 5. 5 (continued)
Return period (years)
•I
5 i•J 25 50 100
"
/;} .lOll:.! .18tit .:.!544 .J45G .4148 .484:.!
j{i .to:>:) .1848 .2526 .3431 .4118 .480i
ii .1041'\ .1835 .2508 .3407 .408!.1 .4ii3
78 .104:.! .1823 .2491 .3383 .4060 .473B
79 .103:> .1810 .2474 .3300 .4032 .4 7011
C-10
EM 1110-2-1415
5 Mar 93
where p = 1/Tr. Assuming, for example, that t = T, = 100 years, then the
probabilities for various values of n are:
n 0 1 2 3 4 5
Pranoo 0.366 0.370 0.185 0.061 0.015 0.003
C-11
EM 1110-2-1415
5 Mar 93
APPENDIX D
HISTORIC DATA 1
Flood information outside that in the systematic record can often be
used to extend the record of the largest events to a historic period much
longer than that of the systematic record. In such a situation, the follow-
ing analytical techniques are used to compute a historically adjusted log-
Pearson Type III frequency curve.
1. Historic knowledge is used to d~fine the historically longer period
of "H" years. The number "Z" of events that are known to be the largest in
the historically longer period "H" are given a weight of 1.0. The remaining
"N" events from the systematic record are given a weight of (H-Z)/(N+L) on the
assumption that their distribution is representative of the (H-Z) remaining
years of the historically longer period.
2. The computations can be done directly by applying the weights to
each individual year's data using equations 6-1, 6-2a, 6-3a, and 6-4a.
Figure 6-1 is an example of this procedure in which there are 44 years of
systematic record and the 1897, 1919 and 1927 floods are known to be the
three largest floods in the 77 year period 1897 to 1973. If statistics have
been previously computed for the current continuous record, they can be
adjusted to give the equivalent historically adjusted values using equations
6-1, 6-2b, 6-3b, and 6-4b, as illustrated in Figure 6-2.
3. The historically adjusted frequency curve is sketched on logarithmic-
probabiiity paper through points established by use of equation 6-5. The
individual flood events should also be plotted for comparison. The histor-
ically adjusted plotting positions for the individual flood events are
computed by use of equation 6-8, in which the historically adjusted order
number of each event 11 rii" is computed from equations 6-6 and 6-7. The com-
putations are illustrated in Figures 6-1 and 6-2, and the completed plotting
1s shown in Figure 6-3.
4. The following example illustrates the steps in application of the
historic peak adjustment only. It does not inc 1ude the fi na 1 step of
weighting with the generalized skew. The historically adjusted skew developed
by this procedure is appropriate to use in developing a generalized skew.
1
Reproduction of Appendix 6 of Bulletin 17B.
D-1
EM 1110-2-1415
5 Mar 93
DEFINITION OF SYMBOLS
E = event number when events are ranked in order from greatest magnitude
to smallest magnitude. The event numbers "E" will range from 1 to
(Z + N).
-G
G = skew coefficient of the X's
= historically adjusted skew coefficient
K = Pearson Type III coordinate expressed in number of standard devia-
tions from the mean for a specified recurrence interval or percent
chance
Q = computed flood flow for a selected recurrence interval or percent
chance
.,
pp = plotting position in percent
..,
p = probability that any peak will exceed the truncation level (used
in step 1, Appendix 5)
z = number of historic peaks including high outliers that have historic
i nfonna ti on
H = number of years in historic period
l = number of low values to be excluded, such as: number of zeros,
number of incomplete record years (below measurable base), and low
outliers which have been identified
a =constant that is characteristic of a given plotting position formula.
For Weibull formula, a= 0; for Beard formula. a = 0.3; and for
Hazen formula, a = 0.5
w = systematic record weight
D-2
EM 1110-2-1415
5 Mar 93
EQUATIONS
w = t!....:..1.
N+ L ( 6-1)
- _ WlX +lXz
M- H-WL (6-2a)
- 2 ~ 2
-2 Wl (X - M) + l (X - M}
s -- {H-WL-1)
z (6-3a)
G = H-WL
(H-WL-1) (H-WL-2)
[Wl(X- M) +l(Xz- M)
3
s
3
J (6-4a)
~ WNM +l X
M - z (6-2b)
- H-WL
2 - 2 - 2
_2 W (N - 1 )S + WN (M - M) + l (Xz - M)
S = (H-WL-1) (6-3b)
Log Q = M + KS (6-5)
rii = E; when: 1 ~ E ~ Z (6-6)
rii =WE- (W- 1) (Z + 0.5); when: (Z +1) ~ E ~ (Z + N+L) (6-7)
pp = rii - a
H + 1 - 2a 100 (6-8)
D-3
EM 1110-2-1415
5 Mar 93
Fii!W'! n-1 HISTORICALLY WEIGHTED LOG PEARSON TYPE !II ·ANNUAL PEAKS
Station: 3-6065, Bic s.u.dy Riwr at Bruceton, TN. D.A. 205 oq...,... mil.,.
Record: 1897, 1919. !927, 1930-1973 <H,......!
Hiftoncal penod: 1897-1973 (77 ,.....,
N • 44; Z • 3; H • 77
....
Yaar
-x -lotr -w •E Nlllllber (Weibu.ll!
pp
s•OC-ill
1897 25,000 4.39794 0.68212 1.00 I 1.00 1.28
1919 21,000 4.32222 0.60640 1.00 2 2.00 2.56
1927 18,600 4.26717 0.65136 1.00 3 3.00 3.85
....
>::::
c;;
19
20
21
22
29.57
31.25
32.94
34.62
37.91
40.06
42.23
44.38
1944 5,340 3.72754 0.01173 23 36.30 46.54
r-'
1951 5,230 3.71850 0.00219 !:: 24 37.98 48.69
1957 5,150 3.71181 ~.00400 25 39.66 50.65
II
1971 5,080 3.70586 ~-00995 26 41.35 53.01
:z; 43.03 55.17
1953 5,000 3.69897 ~.011184 27
1949 4,740 3.675'78 ~.04003
N 28 44.71 57.32
D-4
EM 1110-2-1415
5 Mar 93
figure 6-1. HI5TORICALLY W[IGHT[O LOG PEARSO~·TYPE II! - AN~UAL PEAKS (Continued)
D-5
EM 1110-2-1415
5 Mar 93
Year Y
2
3
(ft ;s) Log v = x lx -
z z I z
M l
I z
(X - M)2
~ ~--- -- ,_ .-.--y--- - -
""5=0.289
- . -- --~-- -t-- -~-- ·: v .I ;
0
~ ~
~ 10000!-----
-
r----- -·-·· r---·· - - - - - - - - -
-·~· --·- - :J;·~-.
,:':
. . I i
I
/v.-
..
APPENDIX E
Mean
t =
X-~.£
(S 2/N)y,
.. 3.8876 - 3. 7782
((0.4681 )2/19)"'
Standard Deviation
Find the 90% confidence interval for the population variance from a sample for which
S2 • 14.4818 and N ""20. Sample values of variance (standard deviation squared) are
distributed like the Chi-square (x2) distribution. From a table (Appendix F-5) for the
Chi-square distribution with 19 degrees of freedom:
Exceedance
Probability Chi-square
0.95 10.117
0.05 30.144
The 90% confidence interval is computed as:
(N - 1) s2 (N-1) 2 s
<O<
X~
2
Xs
19(14.4918) 19(14.14818)
<O<
30.144 10.117
E-1
EM 1110-2-1415
5 Mar 93
APPENDIX F
STATISTICAL TABLES
F-1
EM 1110-2-1415
5 Mar 93
Table F-1
F-2
EM 1110-2-1415
5 Mar 93
F-3
EM 1110-2-1415
5 Mar 93
Table F-2
Positive Skew
Sit. .
Coe!fici~t . . . . . . . . . . . . . . . . . . . . . . ~~t Chance Exeeedance . . . .
<G> 99.0 95.0 90.0 80.0 SO.O 20.0 10.0 5.0 2.0 1.0 0.5 0.2 0.1
0.0 -2.32535 -1.64485 -1.281.55 -o.84162 0.00000 0.114162 1.281.55 1.64485 2.05375 2.32.635 2.57583 2.87816 3.09023
0.1 -2.25258 -1.61594 -1.27037 -o.84611 -o.Ol662 0.83639 1.29178 1.67279 2.10697 2.39961 2.66965 2.99978 3.23322
0.2 -2.17840 -1.58607 -1.25824 -0.84986 -o.03325 0.83044 1.3010.5 1.69971 2.15935 2.47226 2.76321 3.12169 3.37703
0.3 -2.10394 -1.55527 -1.24516 -o.85285 -o.04993 0.82377 1.30936 1.72562 2.21081 2.54421 2.85636 3.24371 3.52139
0.4 -2.02933 -1.5Z357 -1.23114 -0.85S08 -0.06651 0.81638 1.31671 1.7S048 2.26133 2.61539 2.94900 3.36566 3.66608
0.5 -1.95472 -1.49101 -1.21618 -0.85653 -0.08302 0.80829 1.32309 1.77428 2.31084 2.68572 3.04102 3.48737 3.81090
0.6 -1.88029 -1.45762 -1.20028 -0.85718 -o.09945 0.799SO 1.328SO 1.79701 2.35931 2.75514 3.13232 3.60872 3.95567
0.7 -1.80621 -1.42345 -1.18347 -o.85703 -o.11578 0.79002 1.33294 1.81864 2.40670 2.82359 3.22281 3.72957 4.10022
0.8 -1.73271 -1.38855 -1.16574 -o.85607 -o.13199 0.77986 1.33640 1.83916 2.45298 2.89101 3.31243 3.114981 4.24439
0.9 -1.66001 -1.35299 -1.14712 -o.85426 -o.14807 0.7e902 1.33889 1.858S6 2.49811 2.95735 3.40109 3.96932 4.38807
1.0 -1 . .58838 -1.31684 -1.12762 -o.8S161 -o.l6397 0.757~ 1.311039 1.87683 2.54206 3.02256 3.48874 4.08802 4.53112
1.1 -1.51808 -1.28019 -1.10726 -o.84809 -o.17968 0.74537 1.34092 1.89395 2.58480 3.08660 3.57530 4.20582 4.67344
1.2 -1.44942 -1.24313 -1.08608 -0.84369 -o.l9S17 0.73257 1.34047 1.90992 2.62631 3.14944 3.66073 4.322.63 4.81492
1.3 -1.3112.67 -1.20578 -1.06413 -o.83841 -o.21040 0.7191.5 1.33904 1.92472 2.66657 3.21103 3.74497 4.43839 4.95549
1.4 -1.3181.5 -1.16827 -1.04144 -o.83223 -o.22535 0. 70.512 1.33665 1.93836 2. 70.556 3.27134 3.82798 4.55304 5.09S05
1.5 -1.25611 -1.13075 -1.01810 -o.82516 -o.23996 0.690SO 1.33330 1.9S083 2.7432.5 3.33035 3.90973 4.66651 5.23353
1.6 -1.19680 -1.09338 -0.99418 -o.81720 -o.25422 0.67.532 1.32900 1.96213 2.77964 3.38804 3.99016 4.77875 5.37087
1.7 -1.14042 -1.05631 -o.96977 -o.80837 -o.26808 0.659.59 1.32376 1.97227 2.81472 3.44438 4.06926 4.88971 5.50701
1.8 -1.08711 -1.01973 -o.94496 -o.79868 -o.281SO 0.64335 1.31760 1.98124 2.84848 3.49935 4.14700 4.99937 5.64190
1.9 -1.03695 -o.98381 -o.91988 -o.78816 -o.29443 0.62662 1.31054 1.98906 2.88091 3.SS295 4.22336 5.10768 5.77549
2.0 -o.98995 -o.94871 -o.89464 -o.77686 -o.3068S 0.60944 1.30259 1.99573 2.91202 3.60517 4.29832 5.21461 5.90776
2.1 -o.94607 -o.91458 -o.86938 -0.76482 -o.31872 0.59183 1.29377 2.00128 2.94181 3.65600 4.37186 5.32014 6.03865
2.2 -o.90521 -o.88156 -o.84422 -o.7S211 -o.32999 0.57383 1.28412 2.00570 2.97028 3.70543 4.44398 5.42426 6.16816
2.3 -0.86723 -o.84976 -o.81929 -o.73880 -o.34063 o.SSS49 1.27365 2.00903 2.99744 3.7.5347 4.51467 5.52694 6.29626
2.4 -o.83196 -o.81927 -o.79472 -o.7249S -o.3S062 0.53683 1.26240 2.01128 3.02330 3.80013 4.58393 5.62818 6.42292
2.5 -o.79B21 -o.79015 -o.77062 -o.71067 -o.359BZ 0.51789 1.25039 2.01247 3.04787 3.84540 4.65176 5.72796 6.54814
2.6 -o.76878 -o.76242 -o.74709 -o.69602 -o.~ 0.49872 1.23766 2.01263 3.07116 3.88930 4.71815 5.82629 6.67191
2.7 -o.74049 -o.73610 -o.72422 -o.681ll -o.37640 0.47934 1.22422 2.01177 3.09320 3.93183 4.78313 5.92316 6.79421
2.8 -o.71415 -o.71ll6 -o.70209 -o.66603 -o.38353 0.4S980 1.21013 2.00992 3.11399 3.97301 4.84669 6.01858 6.91505
2.9 -o.689S9 -o.68759 -o.6807.5 -o.65086 -o.38991 0.4401.5 1.19539 2.00710 3.13356 4.01286 4.90884 6.112S4 7.03443
3.0 -o.66663 -0.66532 -o.66023 -0.63S69 -o.39SS4 0.42040 1.18006 2.00335 3.15193 4.05138 4.96959 6.20506 7.15235
F-4
EM 1110-2-1415
5 Mar 93
Negative Skew
Sk...
CceWcimt. . . . . . . . . . . . . . . . . . . . . . . Pvcmt Qwx:e ~. . . . . . . . . . . . . .
CGl 99.0 95.0 90.0 80.0 50.0 20.0 10.0 5.0 2.0 1.0 0.5 0.2 0.1
0.0 -2.32535 -1.154485 -1.28155 -0.81ol62 0.00000 0.84162 1.28155 1.154485 2.05375 2.32535 2.57583 2.87816 3.09023
-o.1 -2.391161 -1.67279 -1.29178 -0.83639 0.01662 0.84611 1.27037 1.61.594 1.99973 2.252511 2.48187 2.75706 2.94834
-o.2 -2.47226 -1.69971 -1.30105 -o.83044 o.o3325 0.8491!6 1.2.5824 1.58607 1.94499 2.17840 2.38795 2.63672 2.80786
-0.3 -2.54421 -1.72562 -1.30936 -0.82377 0.04993 0.85285 1.24516 1.55527 1.88959 2.10394 2.29423 2.51741 2.66915
-o.4 -2.61539 -1.75048 -1.31671 -0.81638 0.06651 0.85508 1.23114 1.52357 1.83361 2.02933 2.20092 2.39942 2.53261
-0.5 -2.68572 -1.77428 -1.32309 -0.80829 0.08302 0.85653 1.21618 1.49101 1.77716 1.95472 2.10825 2.28311 2.39867
-0.6 -2.75514 -1.79701 -1.32850 -0.79950 0.09945 0.85718 1.20028 1.45762 1.72033 1.88029 2.011544 2.16884 2.26780
-0.7 -2.82359 -1.818154 -1.33294 -0.79002 0.11578 0.85703 1.18347 1.42345 1.66325 1.80621 1.92580 2.05701 2.14053
-0.8 -2.89101 -1.83916 -1.331540 -o.77986 0.13199 0.85607 1.16574 1.38855 1.60604 1.73271 1.83660 1.94806 2.01739
-o.9 -2.95735 -1.~ -1.33889 -o.76802 0.14807 0.85426 1.14712 1.35299 1.54886 1.66001 1.74919 1.84244 1.89894
-1.0 -3.02256 -1.87683 -1.34039 -o.75752 0.16397 0.85161 1.12762 1.31684 1.49188 1.58838 1.66390 1.74062 1.78572
-1.1 -3.08660 -1.89395 -1.34092 -0.74537 0.17968 0.84809 1.10726 1.28019 1.43529 1.51808 1.58110 1.154305 1.67825
-1.2 -3.14944 -1.90992 -1.34047 -o.73257 0.19517 0.84369 1.08608 1.24313 1.37929 1.44942 1.50114 1.55016 1.57695
-1.3 -3.21103 -1.92472 -1.33904 -0.71915 0.21040 0.83841 1.06413 1.20578 1.32412 1.38267 1.42439 1.46232 1.48216
-1.4 -3.27134 -1.93836 -1.33665 -0.70512 0.22535 0.83223 1.04144 1.16827 1.26999 1.31815 1.35114 1.37981 1.39408
-1.5 -3.33035 -1.95083 -1.33330 -o.69050 0.23996 0.82516 1.01810 1.13075 1.21716 1.25611 1.28167 1.30279 1.31275
-1.6 -3.38804 -1.96Zl3 -1.32900 -0.67532 0.25422 0.81720 0.99418 1.09338 1.16584 1.19680 1.21618 1.23132 1.23805
-1.7 -3.44438 -1.97227 -1.32376 -o.8S959 0.26808 0.80837 0.96977 1.05631 1.11628 1.14042 1.15477 1.16534 1.16974
-1.8 -3.49935 -1.98124 -1.31760 -o.64335 0.28150 0.79868 0.94496 1.01973 1.06864 1.08711 1.09749 1.10465 1.10743
-1.9 -3.5529S -1.98906 -1.31054 -o.62662 0.29443 0.78816 0.91988 0.98381 1.02311 1.03695 1.04427 1.04898 1.05068
-2.0 -3.60517 -1.99573 -1.30259 -o.60944 0.30685 0.77686 0.89464 0.94871 0.97980 0.98995 0.99499 0.99800 0.99900
-2.1 -3.65600 -2.00128 -1.29377 -o.59183 0.31872 0.76482 0.86938 0.91458 0.93878 0.94607 0.94945 0.95131 0.95188
-2.2 -3.70543 -2.00570 -1.28412 -o.57383 0.329119 0.75211 0.84422 0.88156 0.90009 0.90521 0.90742 0.90854 0.90885
-2.3 -3.75347 -2.00903 -1.27365 -o.55549 0.34063 0.73880 0.81929 0.84976 0.86371 0.86723 0.86863 0.86929 0.86945
-2.4 -3.80013 -2.01128 -1.26240 -o.S3683 0.35062 0.72495 0.79472 0.81927 0.82959 0.83196 0.83283 0.83320 0.83328
-2.5 -3.84540 -2.01247 -1.25039 -o.51789 0.35992 0.71067 0.77062 0.79015 0.79765 0.79921 0.79973 0.79994 0.79998
-2.6 -3.88930 -2.01263 -1.23766 -o.49872 0.36852 0.69602 0.74709 0.76242 0.76779 0.76878 0.76909 0.76920 0.76922
-2.7 -3.93183 -2.01177 -1.22422 -0.47934 0.37640 0.68111 0.72422 0.73610 0.73987 0.74049 0.74067 0.74073 0.74074
-2.8 -3.97301 -2.00992 -1.21013 -o.4S980 0.38353 0.66603 0. 70209 o. 71116 0.71377 0. 71415 0. 71425 0. 71428 0. 71428
-2.9 -4.01286 -2.00710 -1.19539 -o.44015 0.38991 0.65086 0.68075 0.68759 0.68935 0.68959 0.68964 0.68965 0.68965
-3.0 -4.05138 -2.00335 -1.18006 -o.42040 0.39554 0.63569 0.66023 0.66S32 0.66649 0.66663 0.66666 0.66667 0.66667
F-5
EM 1110-2-1415
5 Mar 93
Table F-3
NQRMAL QISTRIBU:TIQN
Percent Chance Exceedance for Given Normal Standard Deviate (K)
6.0 .o 7987
Note- Values have been generated by use of computer routines for the normal
distribution.
F-6
EM 1110-2-1415
5 Mar 93
Table F-4
6 .265 .553 .906 1.440 1.943 2.104 2.612 3.143 3.708 4.525 5.208 5.959 7.074 8.025
7 .263 . .549 .896 1.415 1.895 2.046 2.517 2.998 3.500 4.207 4.786 5.408 6.311 7.063
8 .262 . .546 .889 1.397 1.860 2.004 2.449 2.896 3.355 3.991 4.501 5.042 5.812 6.442
9 .261 . .543 .883 1.383 1.833 1.973 2.398 2.821 3.250 3.835 4.297 4.781 5.461 6.011
10 .260 . .542 .879 1.372 1.812 1.948 2.359 2.764 3.169 3.716 4.144 4.587 5.202 5.694
11 .260 . .540 .876 1.363 l. 796 1.928 2.328 2.718 3.106 3.624 4.025 4.437 5.004 5.453
12 .259 .539 .873 1.356 1. 782 1.912 2.303 2.681 3.055 3.550 3.930 4.318 4.847 5.263
13 .259 .538 .870 1.350 1.771 1.899 2.282 2.650 3.012 3.489 3.852 4.221 4.721 5.111
14 .258 .537 .868 1.345 1.761 1.887 2.264 2.624 2.977 3.438 3.787 4.140 4.616 4.985
15 .258 .536 .866 1.341 1.753 1.878 2.249 2.602 2.947 3.395 3.733 4.073 4.528 4.880
16 .258 .535 .865 1.337 1.746 1.869 2.235 2.583 2.921 3.358 3.686 4.015 4. 4.54 4.791
17 .257 .534 .863 1.333 1.740 1.862 2.224 2.567 2.898 3.326 3.646 3.965 4.390 4.714
18 .257 .534 .862 1.330 1.734 1.855 2.214 2.~ 2.878 3.298 3.610 3.922 4.334 4.648
19 .257 .533 .861 1.328 1.729 1.850 2.205 2.539 2.861 3.273 3.579 3.883 4.285 4.590
20 .257 .533 .860 1.325 1.725 1.844 2.197 2.528 2.845 3.251 3.~ 3.850 4.241 4.539
21 .257 .532 .859 1.323 1.721 1.840 2.189 2.518 2.831 3.231 3.527 3.819 4.203 4.493
22 .256 .532 .858 1.321 1.717 1.835 2.183 2.508 2.819 3.214 3.505 3.792 4.168 4.452
23 .256 .532 .858 1.319 1.714 1.832 2.177 2.500 2.807 3.198 3.485 3.768 4.137 4.415
24 .256 .531 .857 1.318 1.711 1.828 2.1n 2.492 2.797 3.183 3.467 3.745 4.109 4.382
25 .256 .531 .856 1.316 1.708 1.825 2.167 2.485 2.787 3.170 3.450 3.725 4.083 4.352
26 .256 .531 .856 1.315 1.706 1.822 2.162 2.479 2.779 3.158 3.435 3.707 4.060 4.324
27 .256 .531 .855 1.314 1.703 1.819 2.158 2.473 2.771 3.147 3.421 3.690 4.038 4.299
28 .256 .530 .855 1.313 1.701 1.817 2.154 2.467 2.763 3.136 3.408 3.674 4.018 4.275
29 .256 .530 .8.54 1.311 1.699 1.814 2.150 2.462 2.756 3.127 3.396 3.659 4.000 4 .2.54
30 .256 .530 .8.54 1.310 1.697 1.812 2.147 2.457 2.750 3.118 3.385 3.646 3.983 4.234
40 .255 .529 .851 1.303 1.684 1.796 2.123 2.423 2.704 3.055 3.307 3.551 3.864 4.094
50 .255 .528 .849 1.299 1.676 1.787 2.109 2.403 2.678 3.018 3.261 3.496 3.795 4.014
60 .2.54 .527 .848 1.296 1.671 1.781 2.099 2.390 2.660 2.994 3.232 3.460 3.750 3.962
70 .2.54 .527 .847 1.294 1.667 1.776 2.093 2.381 2.648 2.977 3.211 3.435 3.719 3.926
80 .2.54 .526 .846 1.292 1.664 1.m 2.088 2.374 2.639 2.964 3.195 3.416 3.696 3.899
90 .2.54 .526 .846 1.291 1.662 1.771 2.084 2.368 2.632 2.9.54 3.183 3.402 3.678 3.878
100 .2.54 .526 .845 1.290 1.660 1.769 2.081 2.364 2.626 2.946 3.174 3.390 3.664 3.862
110 .2.54 .526 .845 1.289 1.659 1.767 2.078 2.361 2.621 2.940 3.166 3.381 3.652 3.848
120 .2.54 .526 .845 1.289 1.658 1.766 2.076 2.358 2.617 2.935 3.160 3.373 3.642 3.837
IliF .253 .524 .842 1.282 1.645 1.751 2.0.54 2.326 2.576 2.878 3.090 3.291 3 . .540 3.719
Note- Values have been generated by use of computer routines for the inverse
t-distribution. A few values for exceedance probabilities of 0.005 and less
may differ plus or minus 0.001 from published tables.
F-7
EM 1110-2-1415
5 Mar 93
Table F-5
Exceecllnce Prcability . ..
OF 0.99 0.95 0,90 0. 75 0 . .50 o,~ 0.10 0.05 0.01 0.005 0.001
1 .o3 ts7 .00393 .0158 .102 .~55 1.323 2.706 3.841 6.635 7.879 10.828
2 .020 .103 .211 .575 1.386 2.m 4.605 5.991 9.210 10.597 13.816
3 .llS .352 .584 1.213 2.366 4.108 6.251 7 .81.5 11.345 12.838 16.256
4 .2S7 . 711 1.064 1.923 3.357 5.385 7.779 9.488 13.277 14.860 18.467
5 .554 1.145 1.610 2.675 4.351 6.626 9.236 11.070 1.5.086 16.7.50 20.51.5
6 .872 1.635 2.204 3.455 5.348 7.841 10.845 12.592 16.812 18.548 22.458
7 1.239 2.167 2.833 4.255 6.346 9.037 12.017 14.067 18.475 20.278 24.322
8 1.846 2.733 3.490 5.071 7.344 10.219 13.362 1.5 . .507 20.090 21.955 26.125
9 2.088 3.325 4.168 5.898 8.343 11.389 14.684 16.919 21.666 23.589 27.877
10 2.558 3.940 4.865 6.737 9.342 12.549 1.5.987 18.307 23.209 25.188 29.588
11 3.053 4.575 5.578 7.584 10.341 13.701 17.275 19.675 24.725 26.757 31.264
12 3.571 5.226 6.304 8.438 11.340 14.845 18.549 21.026 26.217 28.299 32.909
13 4.107 5.892 7.042 9.299 12.340 1.5.984 19.812 22.362 27.688 29.819 34.528
14 4.660 6.571 7.790 10.165 13.339 17.117 21.064 23.685 29.141 31.319 36.123
1S 5.229 7.261 8.547 11.037 14.339 18.245 22.307 24.996 30.578 32.801 37.697
16 5.812 7.962 9.312 11.912 1.5.338 19.369 23.542 26.296 32.000 34.267 39.252
17 6.408 8.672 10.085 12.792 16.338 20.~89 24.769 27.587 33.409 35.719 40.790
18 7.01S 9.390 10.865 13.675 17.338 21.605 25.989 28.869 34.805 37.1S6 42.312
19 7.633 10.117 11.651 14.562 18.338 22.718 27.204 30.144 36.191 38.582 43.820
20 8.260 10.851 12.443 1.5.~52 19.337 23.828 28.412 31.410 37.566 39.997 45.31.5
21 8.897 11.591 16.240 13.344 20.337 24.935 2S.61S 32.671 38.932 41.401 46.797
22 9.542 12.338 14.041 17.240 21.337 26.039 30.813 33.924 40.289 42.796 48.268
23 10.196 13.091 14.848 18.137 22.337 27.141 32.007 35.172 41.638 44.181 49.728
24 10.856 13.848 1.5.659 19.037 23.337 28.241 33.196 36.41.5 42..980 45.559 51.179
25 11.524 14.611 16.473 19.939 24.337 29.339 34.382 37.652 44.314 46.928 52.620
26 12.198 1.5.379 17.282 20.843 25.336 30.435 35.563 38.885 45.842. 48.290 54.052
27 12.879 16.1.51 18.114 21.749 26.336 31.528 36.741 40.113 46.963 49.845 55.476
28 13.!165 16.928 18.939 22.657 27.336 32.621 37.916 41.337 48.278 .50.993 56.892
2S 14.256 17.708 19.758 23.567 28.336 33.711 39.087 42.557 49.588 52.336 58.302.
30 14.953 18.493 20.!189 24.478 29.336 34.800 40.256 43.773 .50.892 53.672 59.703
40 22.184 26 . .509 29.051 33.660 39.335 45.616 51.1105 55.759 63.691 66.766 73.402.
50 29.707 34.764 37.689 42.942 49.335 56.334 63.167 67.505 76.154 79.490 86.661
60 37.485 43.188 lt6.459 52.294 59.335 66.981 74.397 79.082 88.379 91.952 99.607
70 45.442. 51.739 55.32S 61.698 69.334 77.577 85.527 90.531 100.425 104.2.1.5 122.317
80 53.540 60.392 84.2.78 71.1~5 79.334 88.130 96.578 101.879 112.329 116.321 124.839
90 61.754 69.126 73.291 80.625 89.334 98.6.50 107.565 113.145 124.116 128.299 137.208
100 70.065 77.930 82.358 90.133 99.334 109.141 118.498 124.342 135.807 140.169 149.449
F-8
EM 1110-2-1415
5 Mar 93
Table F-6
F-9
EM 1110-2-1415
5 Mar 93
Table F-7
DEVIATES FQR THE &;XP&C:If:Q PRQBAIHLITY ADJUSTMENT
11 .272 •.566 .918 1.433 1.893 2.035 2.464 2.887 3.310 3.881 4.328 4.791 5.434 5.947
12 .270 .562 .911 1.419 1.869 2.007 2.423 2.829 3.233 3.772 4.189 4.618 5.208 5.676
13 .269 .559 .906 1.407 1.850 1.985 2.390 2.782 3.170 3.684 4.078 4.481 5.030 5.462
14 .258 .556 .901 1.398 1.833 1.966 2.362 2.743 3.118 3.611 3.987 4.369 4.886 5.290
15 .267 .554 .897 1.389 1.819 1.949 2.338 2.711 3.074 3.551 3.912 4.276 4.768 5.149
16 .266 .552 .893 1.382 1.B07 1.936 2.318 2.683 3.037 3.499 3.848 4.198 4.668 5.030
17 .265 .551 .890 1.376 1.796 1.923 2.300 2.658 3.005 3.455 3.793 4.131 4.583 4.930
18 .264 .549 .887 1.370 1.787 1.913 2.285 2.637 2.978 3.417 3.746 4.074 4.510 4.844
19 .264 .548 .884 1.365 1.779 1.904 2.271 2.619 2.953 3.384 3.704 4.024 4.446 4.769
20 .263 .546 .882 1.361 Lm 1.895 2.259 2.602 2.932 3.354 3.668 3.979 4.391 4.703
21 .263 .545 .880 1.357 1.765 1.888 2.248 2.587 2.912 3.328 3.635 3.940 4.341 4.645
22 .262 .544 .878 1.353 1.759 1.881 2.239 2.574 2.895 3.304 3.606 3.905 4.297 4.594
23 .262 .544 .877 1.350 1.754 1.875 2.230 2.562 2.879 3.283 3.580 3.874 4.258 4.548
24 .262 .543 .875 1.347 1.749 1.869 2.222 2.551 2.865 3.264 3.557 3.845 4.222 4.506
25 .261 .542 .874 1.344 1.745 1.864 2.215 2.542 2.1152 3.246 3.535 3.820 4.190 4.469
26 .261 .541 .873 1.341 1.741 1.860 2.208 2.532 2.841 3.230 3.516 3.796 4.161 4.435
27 .261 .541 .871 1.339 1.737 1.855 2.202 2.524 2.830 3.216 3.498 3.775 4.134 4.403
28 .260 .540 .870 1.337 1.733 1.851 2.196 2.516 2.820 3.202 3.482 3.755 4.110 4.375
29 .260 .539 .869 1.335 1.730 1.848 2.181 2.509 2.811 3.190 3.466 3.737 4.087 4.349
30 .260 .539 .868 1.333 1.727 1.844 2.186 2.503 2.802 3.178 3.452 3.720 4.066 4.324
40 .258 .535 .862 1.320 1.706 1.820 2.151 2.456 2.742 3.097 3.354 3.602 3.921 4.1.56
50 .257 .533 .857 1.312 1.693 1.B06 2.131 2.429 2.707 3.051 3.298 3.535 3.838 4.061
60 .257 .532 .855 1.307 1.685 1.796 2.118 2.411 2.684 3.020 3.261 3.492 3.785 3.999
70 .256 .531 .853 1.303 1.679 1.nKI 2.108 2.399 2.688 2.999 3.235 3.462 3.748 3.957
ao .256 .530 .851 1.300 1.675 1.785 2.101 2.389 2.656 2.984 3.217 3.439 3.721 3.925
90 .255 .529 .850 1.298 1.671 1.781 2.096 2.382 2.647 2.971 3.202 3.422 3.700 3.901
100 .255 .529 .849 1.297 1.669 1.778 2.092 2.376 2.640 2.962 3.190 3.408 3.683 3.882
110 .255 .528 .849 1.295 1.666 1.775 2.088 2.372 2.634 2.954 3.181 3.397 3.670 3.867
120 .255 .528 .848 1.294 1.665 1.m 2.085 2.368 2.629 2.948 3.173 3.388 3.658 3.854
INF .253 .524 .842 1.282 1.645 1.751 2.054 2.326 2.576 2.878 3.090 3.291 3.540 3.719
Note- Values have been generated by use of computer routines for the inverse
t-distribution. This table is based on samples drawn from a normal distribution.
The above values (after adjustment for skew) may be used as approximate
adjustments to Pearson type-III distributions having small skew coefficients.
F-10
EM 1110-2-1415
5 Mar 93
Table F-8
PER!;:ENT AQ&S FQR THE &XPE!;:TED PRQBABILITY ADJUSTMENT
2 43.5 37.1 30.8 24.3 20.37 19.43 17.12 15.43 14.13 12.81 12.011 11.342 10.602 10.125
3 42.3 34.7 27.1 19.1 14.~ 13.44 10.86 9.08 7.77 6.51 5.793 5.212 4.598 4.219
4 41.8 33.6 25.3 16.7 11.88 10.77 8.18 6.45 5.23 4.11 3.496 3.018 2.531 2.242
5 41.4 32.9 24.3 15.3 10.38 9.26 6.70 5.05 3.92 2.92 2.389 1.990 1.596 1.370
6 41.2 32.4 23.6 14.4 9.41 8.30 .5.78 4.19 3.14 2.23 1.768 1.427 1.101 .919
7 41.0 32.1 23.1 13.8 8.74 7.63 .5.16 3.62 2.63 1.80 1.384 1.086 .809 .658
8 40.9 31.8 22.7 13.3 8.24 7.14 4.70 3.22 2.28 1 . .50 1.127 .863 .623 .495
9 40.8 31.6 22.4 12.9 7.86 6.77 4.36 2.92 2.02 1.29 .947 .710 .498 .388
10 40.7 31.5 22.1 12.6 7.S6 6.47 4.09 2.69 1.82 1.13 .816 .599 .409 .313
11 40.7 31.3 22.0 12.4 7.32 6.23 3.88 2 . .50 1.67 1.01 .716 .516 .345 .259
12 40.6 31.2 21.8 12.2 7.12 6.04 3.71 2.36 1.54 .92 .638 .453 .296 .219
13 40.6 31.1 21.7 12.0 6.9.5 5.87 3.S6 2.23 1.44 .84 .577 .403 .258 .188
14 40.5 31.0 21 . .5 11.9 6.80 5.73 3.44 2.13 1.36 .78 . .525 .363 .228 .164
15 40.5 31.0 21.4 11.8 6.68 5.61 3.33 2.04 1.29 .73 .485 .330 .204 .145
16 40 . .5 30.9 21.3 11.6 6 . .57 5 . .50 3.24 1.97 1.23 .68 .4.50 .303 .184 .129
17 40.4 30.9 21.3 11.5 6.47 5.41 3.16 1.90 1.18 .6.5 .421 .280 .168 .116
18 40.4 30.8 21.2 11.5 6.39 .5.33 3.09 1.85 1.13 .61 .396 .261 .154 .106
19 40.4 30.8 21.1 11.4 6.31 5.26 3.03 1.80 1.09 .59 .374 .244 .143 .097
20 40.4 30.7 21.1 11.3 6.25 5.19 2.98 1. 7.5 1.06 .S6 .355 .230 .133 .089
21 40.4 30.7 21.0 11.2 6.19 5.13 2.93 1. 71 1.02 .54 .339 .217 .124 .083
22 40.3 30.7 21.0 11.2 6.13 5.08 2.88 1.67 1.00 -~ .324 .206 .117 .077
23 40.3 30.6 20.9 11.1 6.08 5.03 2.84 1.64 .97 . .50 .311 .197 .110 .072
24 40.3 30.6 20.9 11.1 6.03 4.99 2.80 1.61 .9.5 .49 .299 .188 .104 .068
25 40.3 30.6 20.9 11.0 .5.99 4.9.5 2.77 1.59 .93 .47 .289 .180 .099 .064
26 40.3 30.6 20.8 11.0 .5.9.5 4.91 2.74 1.S6 .91 .46 .279 .173 .094 .061
27 40.3 30 . .5 20.8 11.0 .5.92 4.87 2.71 1.54 .89 .4.5 .271 .167 .090 .058
28 40.3 30 . .5 20.8 10.9 5.88 4.84 2.68 1.~ .88 .44 .263 .161 .086 .0.55
29 40.3 30.5 20.7 10.9 5.85 4.81 2.66 1. .50 .86 .43 .25.5 .1S6 .083 .0~
30 40.2 30 . .5 20.7 10.9 5.82 4.78 2.63 1.48 .85 .42 .249 .151 .080 .0.50
40 40.2 30.4 20 . .5 10.7 .5.61 4.58 2.47 1.3.5 .75 .3.5 .204 .119 .060 .036
.50 40.1 30.3 20.4 10 •.5 .5.49 4.47 2.37 1.28 .70 .32 .179 .102 .049 .029
60 40.1 30.2 20.4 10.4 5.41 4.39 2.31 1.23 .66 .30 .164 .092 .043 .025
70 40.1 30.2 20.3 10.4 .5.35 4.33 2.26 1.19 .64 .28 .154 .085 .039 .022
80 40.1 30.2 20.3 10.3 .5.30 4.29 2.23 1.17 .62 .27 .146 .080 .036 .020
90 40.1 30.2 20.2 10.3 .5.27 4.26 2.20 1.15 .61 .26 .141 .076 .034 .019
100 40.1 30.1 20.2 10.3 .5.24 4.23 2.18 1.13 .59 .26 .136 .073 .032 .018
110 40.1 30.1 20.2 10.2 5.22 4.21 2.17 1.12 .59 .25 .132 .071 .031 .017
120 40.1 30.1 20.2 10.2 .5.20 4.19 2.15 1.11 .58 .25 .130 .069 .030 .016
INF 40.0 30.0 20.0 10.0 .5.00 4.00 2.00 1.00 . .50 .20 .100 .0.50 .020 .010
Note- Values have been generated by use of computer r.outines for the inverse normal
and inverse t-distributions. This table is based on samples drawn from a normal
distribution. The above values may be used as approximate adjustments to
Pearson type-III distributions having small skew coefficients.
F-11
EM 1110-2-1415
5 Mar 93
Table F-9
Sqll.e
Size 0 Pvcmt. Clun:e Exceedmce o
N 9900 9.500 9000 eooo ~00 2000 10,0 5,0 ~.o 1,0 o5 02 01
Ccllfis!!!!s:• ~! - ,oos
10 -1o232 -0731 -0436 -0027 10028 2o490 3o367 4.121 4.991 So.581 60126 60790 702.58
15 -1.38.5 -0866 -0569 -017.5 0769 lo996 20718 3.336 4.048 40530 40975 5o517 50898
20 -10483 -.951 -08.51 -0259 .640 1.768 2o423 2.982 3.624 4o059 40460 40948 5o2.92
25 -1 •.554 -1.012 -0708 -0317 .SS9 1.632 2.249 2.774 3.378 3078.5 40161 4o618 40940
30 -10608 -1.0.58 -07.51 -0359 .S03 1.539 2.132 2.636 3.214 3o603 30962 40399 40707
40 -1.688 -1.125 -0813 -0419 o428 1.421 1.983 20460 3.00.5 3o372 30710 40122 40412
50 -10 744 -1.172 -08.57 -0460 .379 1.346 1o890 2o350 2.87.5 3.229 3.5.55 3o951 40229
60 -1.787 -1.207 -0890 -0492 .344 1.294 1.825 2.273 2078.5 3.129 30446 30832 40103
70 -1.822 -1.236 -0916 -0.516 .317 10254 1.776 2.216 2o719 30056 30366 3.744 40009
eo -10850 -1.259 -0937 -053.5 .295 1.224 1.738 20172 2o666 2.998 3o304 30675 30936
90 -1.874 -1.279 -09.55 -.552 .277 1o199 10708 2.136 2o624 2.952 30254 3o620 30878
100 -1.894 -1.295 -0970 -osss .263 1.177 1.682 2.106 2o.589 2.914 3o212 30.574 30829
Ccllficltllce lAval. • .~
10 -.5 . .581 -4.121 -3.367 -20490 -1.028 o027 .436 o731 1.038 1.232 1o405 1o61l 10 754
15 -4.530 -3.336 -2.718 -1.996 -0769 .17.5 .569 .81!6 1.182 1.38.5 1.567 1.785 1.936
20 -4.059 -2.982 -2.423 -1.768 -.640 .259 .8.51 .951 1.274 1.483 1.611 10897 2o054
25 -3.78.5 -2.774 -2.21o9 -1.632 -.SS9 o317 0708 1.012 1.340 1. .554 1. 747 1.979 2o140
30 -3.603 -2.636 -20132 -1.539 -.503 .359 o7.51 1.0.58 1.391 1.608 10804 2.040 2.205
loO -3.372 -2.460 -1.983 -10421 -o428 .lo19 o813 1.125 1.1o64 1.688 1o890 2.132 20302
50 -30229 -2.350 -1.890 -1.348 -0379 .460 .8.57 1.172 1 •.518 1.71o4 1.950 2o198 2.371
60 -3.129 -2.273 -1.825 -1.294 -o341o .492 o890 1.207 1.S.58 1.787 1.996 2.248 2o424
70 -3.056 -2.216 -1.776 -1.254 -.317 . .516 o916 1.236 1.590 1.822 2.033 2.288 20466
eo -2.998 -2.172 -1.738 -1.224 -0295 .53.5 o937 1.259 1.615 1.850 20063 2.320 2.SOO
90 -2.952 -2.136 -1.708 -1.199 -.277 .552 .955 1.279 1.637 1.874 2.089 2.348 2.529
100 -2.914 -2.106 -1.682 -1.177 -.263 .S6S .970 1.295 1.8.56 1.894 2.110 2.371 2.554
Ccllficltllce X..V.l • ,010
10 -1.314 -.804 -.508 -.107 .8112 2.243 3.048 3.738 lo.53S .5.074 .5 . .572 6.178 6.605
15 -1.1o.58 -.929 -.829 -.236 .678 1.841 2.521 3.102 3.770 4.222 4.639 .5.147 5 . .504
20 -1.550 -1.008 -.70S -.313 •.568 1.8.51 2.276 2.808 3.419 3.832 4.212 4.67.5 50001
25 -1.616 -1.064 -.757 -.364 .lo98 1.!136 2.129 2.633 3.211 3.601 3.960 4.398 4.706
30 -1.667 -1.107 -.797 -.lo03 .lo!50 1.457 2.030 2.515 3.071 3.1o47 3.792 4.212 4.508
40 -1.741 -1.169 -.854 -.457 .384 1.35.5 1.902 2.364 2.893 3.249 3.577 3.97.5 4.255
50 -1.793 -1.212 -.894 -.496 .340 1.291 1.821 2.2&9 2.781 3.125 3.442 3.826 4.097
60 -1.833 -1.24.5 -.924 -.524 .309 1.21o4 1.764 2.202 2.702 3.038 3.347 3.723 3.987
70 -1.865 -1.272 -.948 - . .545 .285 1.210 1.722 2.153 2.641o 2.974 3.278 3.647 3.906
80 -1.891 -1.293 -.968 -.563 .285 1.183 1.688 2.114 2.599 2.924 3.223 3 . .587 3.842
90 -1.913 -1.311 -.984 -.578 .250 1.160 1.661 2.082 2.561 2.883 3.179 3.538 3.791
100 -1.932 -1.325 -.998 -.391 .236 1.142 1.639 2.056 2.531 2.850 3.143 3.498 3.748
F-12
EM 1110-2-1415
5 Mar 93
F-13
EM 1110-2-1415
5 Mar 93
Cal!idlrlce ~~ • ,900
10 -3.532 -2.S68 -2.066 -1.474 -.437 .429 .828 1.144 1.489 1.715 1.919 2.16.5 2.337
15 -3.212 -2.329 -1.867 -1.320 -.347 .499 .901 1.222 1.576 1.808 2.019 2.273 2.451
20 -3.052 -2.208 -1.765 -1.240 -.297 .541 .91t6 1.271 1.630 1.867 2.082 2.342 2.523
25 -2.952 -2.132 -1.702 -1.190 -.264 .570 .978 1.306 1.689 1.908 2.126 2.390 2 . .574
30 -2.884 -2.080 -1.657 -1.154 -.239 .593 1.002 1.332 1.698 1.940 2.160 2.426 2.612
40 -2.793 -2.010 -1.598 -1.106 -.206 .624 1.036 1.369 1.740 1.986 2.209 2.479 2.668
so -2.73.5 -1.965 -1.559 -1.075 -.184 .645 1.059 1.3116 1.no 2.018 2.244 2.517 2.708
60 -2.694 -1.933 -1.532 -1.052 -.167 .1562 1.on 1.415 1.79Z 2.042 2.269 2.544 2.737
70 -2.1562 -1.909 -1 . .511 -1.03.5 -.155 .674 1.091 1.431 1.810 2.061 2.290 2.567 2.761
80 -2.638 -1.8 -1.495 -1.021 -.144 .684 1.103 1.444 1.824 2.on 2.307 2.585 2.780
90 -2.618 -1.874 -1.481 -1.010 -.136 .693 1.112 1.454 1.836 2.090 2.321 2.600 2.796
100 -2.601 -1.861 -1.470 -1.001 -.129 .701 1.120 1.463 1.847 2.101 2.333 2.613 2.810
Cal!idlrlce ~~ • .~
10 -2.008 -1.382 -1.043 -.625 .222 1.155 1.671 2.104 2.~ 2.927 3.231 3.599 3.858
15 -2.055 -1.422 -1.081 -.661 .179 1.083 1.sn 1.991 2.460 2.775 3.064 3.415 3.661
20 -2.085 -1.448 -1.104 -.683 .154 1.04.5 1 •.528 1.1132 2.390 2.697 2.978 3.320 3.560
25 -2.106 -1.466 -1.121 -.689 .137 1.020 1.497 1.895 2.346 2.648 2.925 3.261 3.497
30 -2.123 -1.479 -1.133 -.710 .125 1.002 1.47.5 1.8611 2.31.5 2.614 2.888 3.220 3. 4.53
40 -2.147 -1.499 -1.1.51 -.726 .108 .978 1.44.5 1.834 2.274 2.S68 2.838 3.165 3.395
50 -2.183 -1 •.513 -1.164 -.738 .096 .962 1.425 1.811 2.247 2 . .538 2.805 3.129 3.3.57
60 -2.176 -1 •.523 -1.173 -.747 .088 .9.50 1.411 1.795 2.227 2•.517 2.783 3.105 3.331
70 -2.186 -1.532 -1.181 -. 7.53 .081 .942 1.401 1.782 2.213 2.S01 2. 76.5 3.085 3.310
80 -2.194 -1 . .538 -1.187 -.759 .076 .93.5 1.392 1.772 2.202 2.489 2.752 3.070 3.294
90 -2.201 -1.544 -1.192 -.763 .071 .929 1.386 1.764 2.192 2.478 2.740 3.058 3.281
100 -2.207 -1.549 -1.196 -.767 .068 .925 1.380 1.758 2.184 2.470 2.731 3.048 3.270
Cal!idge ~1 • ,750
10 -2.927 -2.104 -1.671 -1.155 -.222 .625 1.043 1.382 1.759 2.008 2.235 2.508 2.700
1.5 -2.ns -1.991 -1 •.577 -1.083 -.179 .661 1.081 1.422 1.803 2.0.5.5 2.284 2.562 2.756
20 -2.697 -1.1132 -1 •.528 -1.04.5 -.154 .683 1.104 1.448 1.831 2.08.5 2.317 2.597 2.793
25 -2.648 -1.89.5 -1.497 -1.020 -.137 .699 1.121 1.466 1.6.51 2.106 2.339 2.621 2.819
30 -2.614 -1.8611 -1.475 -1.002 -.125 .710 1.133 1.479 1.867 2.123 2.357 2.641 2.839
40 -2 . .568 -1.834 -1.445 -.978 -.108 .726 1.151 1.499 1.888 2.147 2.383 2.668 2.868
so -2 . .538 -1.811 -1.425 -.962 -.096 .738 1.164 1•.513 1.903 2.183 2.400 2.688 2.889
60 -2 . .517 -1.795 -1.411 -.9.50 -.088 .747 1.173 1•.523 1.916 2.176 2.414 2.702 2.904
70 -2.S01 -1.782 -1.401 -.942 -.081 . 7.53 1.181 1.532 1.925 2.186 2.425 2.714 2.917
80 -2.489 -1.772 -1.392 -.935 -.076 .759 1.187 1..538 1.1132 2.194 2.434 2.724 2.927
90 -2.478 -1.764 -1.386 -.929 -.071 .763 1.192 1.544 1.938 2.201 2.441 2.731 2.935
100 -2.470 -1.758 -1.380 -.925 -.068 .767 1.196 1.549 1.944 2.207 2.447 2.739 2.943
F-14
EM 1110-2-1415
5 Mar 93
F-15
EM 1110-2-1415
5 Mar 93
Tab1eF-IO
F-16
EM 1110-2-1415
5 Mar 93
Table F-11
Note- Table contains one sided 10% significance level deviates for the normal
distribution. Source: Appendix 4, Bulletin 17B, reference (46).
F-17
EM 1110-2-1415
5 Mar 93
Table F-12
BINOMIAL RISK TABLES
Exceedances 10 20 30 40 50 60 80 100
F-18
EM 1110-2-141:
5 Mar 9:
F-19
EM 1110-2-1415
5 Mar 93
F-20