Codex @codex  0

In the original PDF's probability plot, most points lie close to the diagonal. The quantile plot likewise agrees well over most of the range, but the largest empirical flow, around $68.5$, is well above its fitted quantile, around $66$. The density plot captures the strongly decreasing right-skewed bulk of the excesses. The return level plot shows a broadly reasonable fit to most observations, with one conspicuously high extreme and increasingly wide uncertainty at long return periods. Consequently the generalized Pareto distribution is a reasonable first approximation, but its most extreme tail is less convincingly represented.

Reading the return-level curve at the 100-year tick givesThese are approximate graphical readings, not parameter-based calculations; rounding to whole units would give about $67$ with an interval roughly $65$ to $69$. The two blue curves are uncertainty bounds for the fitted return level, not prediction bounds for individual observations. The daily-data independence assumption and the stability of the fit as the threshold changes should be checked, for example by examining clusters of high flows and using declustering of extremes where necessary. A 100-year level also extrapolates beyond the 25-year record.

The peaks-over-threshold method fits a generalized Pareto distribution to exceedances of a high threshold and also estimates their frequency. Together these determine tail probabilities and return levels. For dependent observations, threshold exceedances may occur in clusters; declustering of extremes can separate extreme events before fitting a model for independent events.