Skip to contents

Given an Open-Meteo marine forecast tibble, computes per-station summaries: the peak forecast hour, peak H_s, peak P(H_max > 20 m), peak P(H_max > 25 m). Uses p_hmax_exceedance() applied independently per forecast hour, then takes the maximum across the forecast horizon.

This is a forecast-derived risk surrogate — it should always be presented alongside the deterministic-NWP caveat (lead-time skill drops sharply after day 2-3, no ensemble spread).

Usage

summarise_forecast_rogue_risk(
  marine_forecasts,
  thresholds = c(10, 15, 20, 25),
  duration_s = 3600
)

Arguments

marine_forecasts

Tibble from fetch_all_marine_forecasts().

thresholds

Numeric vector of H_max thresholds in metres (default c(20, 25)).

duration_s

Window duration in seconds for each forecast hour (default 3600).

Value

Tibble with one row per station: station_id, peak_time, peak_hs_m, peak_period_s, p_hmax_gt_10, p_hmax_gt_15, p_hmax_gt_20, p_hmax_gt_25, n_forecast_hours. Empty tibble if marine_forecasts is empty.