Return period explained

A return period, also known as a recurrence interval or repeat interval, is an average time or an estimated average time between events such as earthquakes, floods,[1] landslides,[2] or river discharge flows to occur.

It is a statistical measurement typically based on historic data over an extended period, and is used usually for risk analysis. Examples include deciding whether a project should be allowed to go forward in a zone of a certain risk or designing structures to withstand events with a certain return period. The following analysis assumes that the probability of the event occurring does not vary over time and is independent of past events.

Estimating a return period

Recurrence interval

={n+1\overm}

n number of years on record;

m is the rank of observed occurrences when arranged in descending order[3] For floods, the event may be measured in terms of m3/s or height; for storm surges, in terms of the height of the surge, and similarly for other events. This is Weibull's Formula.[4] [5]

Return period as the reciprocal of expected frequency

The theoretical return period between occurrences is the inverse of the average frequency of occurrence. For example, a 10-year flood has a 1/10 = 0.1 or 10% chance of being exceeded in any one year and a 50-year flood has a 0.02 or 2% chance of being exceeded in any one year.

This does not mean that a 100-year flood will happen regularly every 100 years, or only once in 100 years. Despite the connotations of the name "return period". In any given 100-year period, a 100-year event may occur once, twice, more, or not at all, and each outcome has a probability that can be computed as below.

Also, the estimated return period below is a statistic: it is computed from a set of data (the observations), as distinct from the theoretical value in an idealized distribution. One does not actually know that a certain or greater magnitude happens with 1% probability, only that it has been observed exactly once in 100 years.

That distinction is significant because there are few observations of rare events: for instance, if observations go back 400 years, the most extreme event (a 400-year event by the statistical definition) may later be classed, on longer observation, as a 200-year event (if a comparable event immediately occurs) or a 500-year event (if no comparable event occurs for a further 100 years).

Further, one cannot determine the size of a 1000-year event based on such records alone but instead must use a statistical model to predict the magnitude of such an (unobserved) event. Even if the historic return interval is a lot less than 1000 years, if there are a number of less-severe events of a similar nature recorded, the use of such a model is likely to provide useful information to help estimate the future return interval.

Probability distributions

One would like to be able to interpret the return period in probabilistic models. The most logical interpretation for this is to take the return period as the counting rate in a Poisson distribution since it is the expectation value of the rate of occurrences. An alternative interpretation is to take it as the probability for a yearly Bernoulli trial in the binomial distribution. That is disfavoured because each year does not represent an independent Bernoulli trial but is an arbitrary measure of time. This question is mainly academic as the results obtained will be similar under both the Poisson and binomial interpretations.

Poisson

The probability mass function of the Poisson distribution is

P(r;t)={(\mut)r\overr!}e-\mu={(t/T)r\overr!}e-t/T

where

r

is the number of occurrences the probability is calculated for,

t

the time period of interest,

T

is the return period and

\mu=1/T

is the counting rate.

The probability of no-occurrence can be obtained simply considering the case for

r=0

. The formula is

P(r=0;t)=e-\mu=e-t/T

Consequently, the probability of exceedance (i.e. the probability of an event "stronger" than the event with return period

T

to occur at least once within the time period of interest) is

P(t>0;t)=1-P(t=0;t)=1-e-\mu=1-e-t/T

Note that for any event with return period

T

, the probability of exceedance within an interval equal to the return period (i.e.

t=T

) is independent from the return period and it is equal to

1-\exp(-1)63.2\%

. This means, for example, that there is a 63.2% probability of a flood larger than the 50-year return flood to occur within any period of 50 year.

Example

If the return period of occurrence T is 243 years (\mu = 0.0041) then the probability of exactly one occurrence in ten years is

\begin{align} P(r;t)&=

(\mut)r
r!

e-\mu\\[6pt] P(r=1;t=10)&=

(10/243)1
1!

e-10/2433.95\% \end{align}

Binomial

In a given period of

n x \tau

for a unit time

\tau

(e.g.

\tau=1year

), the probability of a given number r of events of a return period

\mu

is given by the binomial distribution as follows.

P(X=r)={n\chooser}\mur(1-\mu)n-r.

This is valid only if the probability of more than one occurrence per unit time

\tau

is zero. Often that is a close approximation, in which case the probabilities yielded by this formula hold approximately.

If

ninfty,\mu0

in such a way that

n\muλ

then
n!
(n-r)!r!

\mur(1-\mu)n-re

λr
r!

.

Take

\mu=

1
T

={m\overn+1}

where

T is return interval

n is number of years on record.

m is the number of recorded occurrences of the event being considered

Example

Given that the return period of an event is 100 years,

p={1\over100}=0.01.

So the probability that such an event occurs exactly once in 10 successive years is:

\begin{align} P(X=1)&=\binom{10}{1} x 0.011 x 0.999\\[4pt] &10 x 0.01 x 0.914\\[4pt] &0.0914 \end{align}

Risk analysis

Return period is useful for risk analysis (such as natural, inherent, or hydrologic risk of failure).[6] When dealing with structure design expectations, the return period is useful in calculating the riskiness of the structure.

The probability of at least one event that exceeds design limits during the expected life of the structure is the complement of the probability that no events occur which exceed design limits.

The equation for assessing this parameter is

\overlineR=1-\left(1-{1\overT}\right)n=1-(1-P(X\ge

n
x
T))
where

{1\overT}=P(X\gexT)

is the expression for the probability of the occurrence of the event in question in a year;

n is the expected life of the structure.

See also

Notes and References

  1. Book: ASCE, Task Committee on Hydrology Handbook of Management Group D of. Hydrology Handbook Books. en. 10.1061/9780784401385. 1996. 978-0-7844-0138-5.
  2. Peres. D. J.. Cancelliere. A.. 2016-10-01. Estimating return period of landslide triggering by Monte Carlo simulation. Journal of Hydrology. Flash floods, hydro-geomorphic response and risk management. 541. 256–271. 10.1016/j.jhydrol.2016.03.036. 2016JHyd..541..256P .
  3. Kumar. Rajneesh. Bhardwaj. Anil. 2015. Probability analysis of return period of daily maximum rainfall in annual data set of Ludhiana, Punjab. Indian Journal of Agricultural Research. en. 49. 2. 160. 10.5958/0976-058X.2015.00023.2. 0367-8245.
  4. Book: https://directives.sc.egov.usda.gov/OpenNonWebContent.aspx?content=17781.wba. 7 February 2023. Chapter 5: Stream Hydrology. National Engineering Handbook, Part 654: Stream Restoration Design. National Resources Conservation Service. U.S. Department of Agriculture. August 2007. Washington, D.C..
  5. Web site: Flood Estimation Handbook. Anonymous. 2014-11-07. UK Centre for Ecology & Hydrology. en. 2019-12-21.
  6. Water Resources Engineering, 2005 Edition, John Wiley & Sons, Inc, 2005.