Return level ($z_T=F_{\mathrm{annual}}^{-1}(1-1/T)$)

A return level answers the inverse question to a return period: which threshold has the specified rarity? Let $F(z)=F_{\mathrm{annual}}(z)=\mathbb{P}(M_j\le z)$ be the common cumulative distribution function of annual sample maxima in a model with stationarity. For a continuous, strictly increasing $F$ and $T>1$, the $T$-year return level isIt has the same units as $M_j$, whereas $T$ 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 $F^{-1}(q)=\inf \{z:F(z)\ge q\}$. Then $\mathbb{P}(M_j>z_T)\le 1/T$, 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 $\mu$ and scale $b>0$:The second formula follows by taking logarithms twice and solving $F(z_T)=1-1/T$. 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 $\mu =10$ and $b=2$ in arbitrary units, $z_{10}=14.501$, $z_{100}=19.200$, and $z_{1000}=23.815$, rounded to three decimals. In particular, $19.200$ is approximately the annual quantile at cumulative probability $0.99$; 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 $z_{100}$ is $0.99^{100}\approx 0.366$.

The time block matters. If a toy year contains $m$ independent daily blocks, each with exceedance probability $q$ at one fixed threshold, its annual exceedance probability is $p=1-(1-q{)}^m$, giving annual return period $1/[1-(1-q{)}^m]$. The reciprocal $1/q$ 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.

A return level plot graphs $T$ horizontally and $z_T$ vertically, often using a logarithmic scale for $T$. To read it, choose a return period, move up to the curve, and read the corresponding return level. Conversely, a threshold's horizontal coordinate gives the reciprocal of its annual exceedance probability. The usual plot of fitted quantiles against return periods, and optional pointwise confidence intervals, are documented in Mathieu Ribatet, POT package: return level plot.

The near-straight shape has a mathematical explanation: since $-\log (1-1/T)\sim 1/T$ as $T\to \infty$, the Gumbel distribution gives $z_T=\mu +b\log T+o(1)$. Multiplying a large return period by ten therefore increases the return level by approximately $b\log 10$, about $4.605$ units for this toy curve. A curve indexed by rarity puts several thresholds on one common scale while retaining their measurement units.

For a fitted curve, a pointwise confidence interval concerns uncertainty in the estimated quantile at a specified $T$; it is distinct from a prediction interval for a future annual sample maximum, and pointwise intervals need not form a simultaneous confidence band. A far-right curve value can also be an extrapolation beyond the observation record: extending the line adds no observations to support the chosen tail model. The diagram above uses fixed toy parameters, so it shows neither estimation uncertainty nor an observational record. Its interpretation throughout assumes the stated annual block and stationarity.