A return level answers the inverse question to a return period: which threshold has the specified rarity? Let be the common cumulative distribution function of annual sample maxima in a model with stationarity. For a continuous, strictly increasing and , the -year return level isIt has the same units as , whereas is a count of annual blocks. Independence between years is unnecessary for this one-year quantile definition; it is needed for the geometric waiting and multi-year formulas under return period. For a discrete probability distribution, use the generalized quantile function . Then , with equality not guaranteed at jumps; an exact requested exceedance probability may be unattainable.
For a concrete standard model, choose the maxima convention of the Gumbel distribution, with location and scale :The second formula follows by taking logarithms twice and solving . The maxima cumulative distribution function and its inverse are given in NIST: extreme value type I distribution, maximum case. With the deliberately chosen toy values and in arbitrary units, , , and , rounded to three decimals. In particular, is approximately the annual quantile at cumulative probability ; it is not a prediction of the largest observation in the next 100 years. Under independence, the probability that all those annual sample maxima stay below the exact is .
The time block matters. If a toy year contains independent daily blocks, each with exceedance probability at one fixed threshold, its annual exceedance probability is , giving annual return period . The reciprocal counts daily blocks. Dependence among daily exceedances requires their joint law instead; a threshold crossed on several successive days still produces only one exceedance indicator for that year's sample maximum. The peaks-over-threshold method and declustering of extremes address a different choice of extreme-event sampling.
New to topics? Read the docs here!