Two Steps to a Sharper Forecast of Cloudbursts ✎

Original: "A two-step approach to model precipitation extremes in California based on max-stable and marginal point processes"
· Hongwei Shang, Jun Yan, Xuebin Zhang
arXiv:1204.0286 · 2012-04-02 · arXiv non-exclusive · independent analysis · 2 min
How daily weather-station records become a risk map for rare floods.

On an old film strip of rainfall, only the brightest frames are visible, and between them lies noise. A meteorologist assessing flood risk long looked only at such frames, that is, at annual rainfall maxima. Everything else, thousands of daily observations, remained off-screen. The new method develops the entire film, and details that were previously invisible come into view.

The first step is like tuning each microphone in a studio before recording an orchestra. For each of California's 35 weather stations, a curve is fitted individually to describe the tail of the rainfall distribution. In the tail live rare but destructive downpours. What matters here is not just the average rainfall, but how quickly the probability of extreme values decays. The fitting uses daily threshold exceedances rather than annual peaks, so the curve is supported by hundreds of points instead of 55.

Key idea: first clean each signal individually, and only then tune the connections between them. That way, noise from one sensor does not corrupt the calibration of its neighbor.

The second step brings the connections between stations into focus. Rain in Napa and rain in Winters often arrive together because they are born from the same atmospheric front. To measure this connection, the method no longer uses daily records but annual maxima at pairs of stations. Pairwise comparison is like listening to duets: if two instruments always sound loud at the same time, there is a common cause between them. This is how the synchrony of extremes at different points on the map is assessed.

The result is visible on the risk map. The traditional approach, relying only on annual maxima, found a significant El Niño influence at 14 stations. The new method detected this effect at 22 stations. Confidence intervals for joint 50-year return levels narrowed by 27–51%. This means that an insurer or a city planner gets not a blurred patch but a sharp line marking where the zone of acceptable risk ends.

The efficiency of the traditional method was estimated at 23% to 95%, depending on the parameter. For the tail shape, the most sensitive parameter, efficiency dropped to 45%.

The old film is finally fully developed. Daily records, once dismissed as noise, turned out to be signal. The method does not require complex multivariate models and fits easily into the familiar toolkit of hydrologists. The next frame can already be shot with a more sensitive film, and even the rarest downpours will no longer come as a surprise.

The joint 50-year return level for the station pair Napa and Winters at a Southern Oscillation Index value of -1 was about 10–15 cm of precipitation per day, depending on the estimation method.

G(z) = \exp\left\{-\left[1 + \xi\left(\frac{z-\mu}{\sigma}\right)\right]^{-1/\xi}\right\}
Distribution of block maxima, where μ is the location parameter, σ the scale, and ξ the shape.
\Pr\{Z(x_1) \le z, \ldots, Z(x_p) \le z\} = \exp(-\zeta/z)
Measure of extremal dependence: ζ=1 for complete dependence, ζ=p for independence.
Concepts
Multivariate ENSO Index NOAA CPC maximum likelihood Oceanic Niño index ERSST Niño 1+2 probabilistic forecast ENSO precipitation anomaly relative ONI ensemble forecast bayesian evidence zebiak cane model likelihood ratio test spring predictability barrier probabilistic model nested sampling Monte Carlo methods Analogue forecast parameter estimation