The Science Behind SAROS
Search decisions carry real consequences, so a planning tool must be able to answer two questions plainly: where do its numbers come from, and would it give the same answer twice? This page sets out the models SAROS implements, the published sources they derive from, how the computation is kept reproducible and auditable, and how the results have been tested against real incidents.
Every behavioural distribution, leeway coefficient and detection model traces to a named published source — Koester, the Allen leeway field studies, UK Mountain Rescue research, Grampian Police case data, the IAMSAR Manual. The citations live in the source code beside the numbers they justify.
Monte Carlo does not mean arbitrary. Drift and Aero ensembles run from a recorded random seed: the same inputs and seed reproduce the same particles, the same contours, the same plan — reviewable months later in a debrief or inquiry.
Live drift trials with a police marine unit, 21 real incidents from the published Koester casebook re-run retrospectively, IAMSAR sweep-width tables verified cell-by-cell, and ~3,800 automated tests executed before every release.
SAROS is decision support for qualified operators, not an oracle. Every output declares its assumptions, every simplification is stated, and residual probability outside the searched area is always shown. Absence of probability is never treated as evidence of absence.
How a drifting object moves is one of the most-studied problems in SAR science. SAROS implements the leeway approach established by the US Coast Guard field-study programme and codified in the IAMSAR Manual.
A drifting object moves with the water it sits in, plus a wind-driven component — leeway — that depends on how much of the object faces the wind above and below the waterline. Decades of field experiments — principally Allen & Plourde (1999) and Allen (2005), refined by Breivik et al. (2011) — measured this relationship for real SAR objects by instrumenting them and recording drift against measured wind. SAROS decomposes leeway into downwind and crosswind components relative to the wind, on top of tidal and ocean-current advection.
SAROS carries a versioned object-class library — 15 classes in four groups: persons in water in different postures, ballasted and unballasted liferafts, small craft and kayaks, surfboards and debris — each with its empirically derived downwind slope and offset, leeway uncertainty, crosswind divergence angle δ and jibing rate. One class, the sit-on-top kayak, is additionally calibrated directly against a physical drift trial SAROS ran with NZ Coastguard in the Hauraki Gulf (15.1 km drift over 5.8 hours).
Particles deployed across the Last Known Position uncertainty, from a recorded random seed
Each particle advected by measured wind, tide and current, with field-study leeway coefficients for the object class
Probability contours computed from where the ensemble actually went — priority regions, not a single guess
Two live targets — a kayak and a person in water — drifted 15 km and 7.5 km over 8 hours of Auckland tide and offshore wind. Both were recovered inside the SAROS-planned search area, 0.7 and 0.27 miles from the predicted positions.
Where a lost person goes is not random. It follows patterns measured across thousands of resolved incidents — and that measured behaviour, not intuition, is what SAROS simulates.
The published percentile tables define, for each subject category, how far subjects are found from the point last seen. SAROS fits a log-normal distribution to those percentiles and runs thousands of simulated journeys against the real terrain.
Twenty-one real search incidents documented in Koester's published casebook — dementia, autism, children, hikers, hunters, despondent subjects — were re-run through SAROS from the original initial planning point. The find location fell within the statistically expected range, and the P95 contour encompassed the find, at the rates shown below.
A hiker in the Cairngorms does not behave like a hiker in the Rockies. SAROS is one of the few tools that carries UK-calibrated tables — Perkins & Roberts and the Grampian dataset — alongside the US ISRID-derived values, with a documented fallback chain when a category lacks regional data.
For an overdue aircraft the physics is unforgiving and therefore tractable: a gliding aircraft trades height for distance at a known ratio, drifting with the wind as it descends.
The Monte Carlo layer then perturbs what is genuinely uncertain — position, altitude, heading, glide performance, wind — across ensembles of up to 20,000 particles, producing probability heatmaps and nested containment contours. The probability accounting is exact by construction: cell masses plus the residual outside the grid sum to one; nothing is renormalised away.
The glide engine is verified against a benchmark set of reference cases — including the specification's worked examples — and must reproduce them within 1% on every test run. The Monte Carlo layer is seeded: if no seed is supplied, one is generated and recorded with the results, so any ensemble can be re-run identically.
The International Aeronautical and Maritime SAR Manual is the operational standard. SAROS implements it directly rather than approximating it.
The IAMSAR Vol III uncorrected sweep width tables for merchant vessels, helicopters and fixed-wing aircraft (Tables 3-19 to 3-21) are carried as data with their source recorded — and an automated test validates every cell against the published values, including all weather-correction combinations. Wc = Wu × fw (IAMSAR corrected sweep width)
Probability of detection uses the IAMSAR exponential (random-search) detection model from a single shared implementation across all three domains — the same POD mathematics governs a vessel sweep, a hillside segment and an air search area. POD = 1 − e−C, C = Wc / S (IAMSAR coverage)
Parallel track, expanding square and sector search generated per IAMSAR Vol III Chapter 4 geometry — sized to the particle spread at commence-search time, with waypoints, bearings, leg times and estimated duration on an exportable search card. Sector legs R, R, 2R, R, 2R, R, R (9R total) · Expanding square S, S, 2S, 2S…
"Would it give the same answer twice?" — the property that turns a model into defensible decision support.
Drift and Aero Monte Carlo ensembles run from an explicit random seed carried with the results. Identical inputs and seed produce identical particles, contours and search areas — asserted by automated tests on every release. Land search validation is statistical, measured by percentile containment across the published case benchmark.
Every run is stamped with the model version and the object-class library version that produced it. Case assumptions are versioned append-only — author, time, reason and a full snapshot — so a plan can be reconstructed exactly as it stood at any decision point.
Evidence records (radar, ADS-B, witness, debris…) carry position, time, uncertainty and confidence. Excluding evidence requires a recorded reason — it is flagged out, never deleted. Environmental data provenance records which provider supplied each forcing field and any fallback taken.
SAROS is decision support and is engineered to stay that way: acknowledgement gates on ground-frame patterns in set, nearshore model limits and residual risk; assumptions printed on every output; a full audit trail on all operations; and ~3,800 automated tests run before every release.
Four layers, from physics to the real world.
Physics engines reproduce published worked examples and golden benchmark cases within stated tolerance (glide: 1%) on every automated test run.
IAMSAR sweep-width and weather-correction tables validated cell-by-cell against the manual. Behavioural tables carry their bibliographic source in code.
21 published Koester casebook incidents re-run from the original planning point: 75% of finds within P50, 81% within P95.
Instrumented drift targets with the NZ Police Marine Unit: both recovered inside the SAROS search area. Object-class coefficients calibrated against the measured tracks.
SAROS's own validation reports — the NZ Police live trial results, the 21-case Koester retrospective benchmark and the sweep-width verification suite — are available to prospective customers and professional reviewers on request.
A model you can trust is one that tells you where it stops. These limitations are documented in the SAROS source code and on the relevant outputs; the roadmap items among them are tracked openly.
| Decision support | SAROS assists qualified human operators. It does not replace SAR doctrine, training or command judgement, and its terms of use make acknowledging this a condition of access. |
| Conditional probability | All displayed probability is conditional on the modelled scenario set. The subject may lie outside every modelled scenario — residual probability is always displayed, never hidden. |
| Aero performance data | Aircraft profiles are generic engineering values, not type-certified data — stated on every output. Retrospective validation against published accident investigations is a declared roadmap requirement. |
| Behavioural coverage | Distance tables for all categories are published-source; behaviour-state weights for some categories still use general Koester heuristics pending full per-category literature review — tracked as an open backlog item. |
| Weather uncertainty | Forecast wind and current are applied as provided; forecast-error propagation through the ensemble is a roadmap item. Provider and fallback provenance are recorded on every run. |
We welcome scrutiny from SAR professionals, scientific advisors and procurement teams — including live walkthroughs of the models and validation evidence.