Choose how to read

The full view retains the technical content, figures and sources. The plain-language explanation introduces the ideas without changing the findings.

Terms explained

Adrian Dunkley · Physics-informed numerical verification

Can the solver follow the equation?

Actual Jamaican rainfall inputs, an assumed storage model and independently checked numerical results.

What these results mean

This is a retrospective equation-solver study, not a flood forecast. The rainfall inputs are satellite-plus-gauge estimates. Storage and outflow are simulations: there are zero observed runoff targets. Small numerical errors do not validate the assumed runoff fraction, residence time or flood impacts.

What was compared

Twelve configurations were evaluated across three seeds. Selection chose 16 fixed hidden features, boundary weight 1 and ridge row amplitude 0.000001 (the squared-penalty coefficient is 0.000000000001).

  • 25 June–2 July 2024: warm-up history.
  • 3–6 July: configuration selection against an analytic storage solution.
  • 7–10 July: subsequent numerical evaluation.

Each day's output weights are fitted using that day's known rainfall, including evaluation days. This split evaluates configuration choice, not a causal future forecast. The assumed runoff fraction is 0.3 and residence time is one day.

Selected neural storage follows the analytic assumed solution. Wrong decay creates large differences, while the zero-boundary diagnostic remains zero. The second panel shows errors against the analytic scenario on a logarithmic scale, with selection and evaluation periods shaded.
Seed 23 illustration. Open the figure to zoom. All seeds are reported below. These are simulated storage errors, not measured flood errors.

All-seed evaluation results

Root mean square error measures disagreement with the analytic scenario. The physics residual measures violation of the intended storage equation. Water-budget errors include jumps from the propagated initial state.

DiagnosticSeedStorage error (mm)Equation residual (mm/day)Max daily budget error (mm)Max continuity jump (mm)
selected physics234.16e-089.84e-077.45e-081.57e-11
selected physics376.25e-081.24e-066.7e-088.48e-11
selected physics531.46e-083.68e-072.47e-081.57e-11
wrong decay2318.717.920.71.59e-11
wrong decay3718.717.920.77.17e-10
wrong decay5318.717.920.71.49e-11
zero boundary diagnostic234.13.717.280
zero boundary diagnostic374.13.717.280
zero boundary diagnostic534.13.717.280

Download all-seed results (CSV) · Download configurations and metrics (JSON)

The wrong-decay diagnostic deliberately uses an incorrect parameter. The zero-boundary diagnostic is identically zero, not a data-trained competing model. No runtime or predictive superiority is established.

What would establish flood-prediction skill?

Independent streamflow or stage records, rating curves, subdaily rainfall, channel and drainage geometry, and checked inundation observations are needed. Daily rainfall cannot establish a flash flood's within-day peak or arrival time. A future catchment experiment must evaluate unseen events with equal issue-time information and report peak timing, magnitude, inundation, uncertainty and computational cost.

This short, retrospective search uses one rainfall sequence. Boundary weights are nearly tied; the selected configuration is not a universally optimal network. The protocol is generated by the executable, not a preregistered experiment.

Data provenance: CHIRPS daily rainfall estimates, 25 June–10 July 2024; terrain-derived input selection uses Copernicus surface elevation. See the accompanying paper and source records for acquisition details and attribution. Independent checks verify numerical reconstruction and integrated balance, not field performance.

Terms explained

Select a term to read its meaning in this study.

Physics-informed neural network

A neural network fitted partly by checking how well its answer follows a chosen physical equation. A small equation error does not prove the equation describes the real situation.

Manufactured solution

A deliberately chosen mathematical answer used to build a test equation. It allows exact error checks but is not an observed event.

Forcing

An input that drives a model. Rainfall drives the storage calculation; the first mathematical test uses a constructed input instead.

CHIRPS

Climate Hazards Group InfraRed Precipitation with Station data: rainfall estimates combining satellite information and rain-gauge records. They are estimates, not perfect measurements at every point.

Reservoir

Here, a mathematical container that gains water from rainfall and loses stored water at an assumed rate. Each map cell is separate, so water is not routed through a river network.

Residual

The amount by which a calculated answer fails to satisfy the chosen equation. Zero means agreement with that equation, not necessarily agreement with nature.

Root mean square error

A measure of how far answers are from their comparison values, giving more weight to large errors. Here the comparison is a mathematical solution, not measured floodwater.

Seed

A recorded number that makes the network's random starting choices repeatable. Different seeds help check whether a result depends on one lucky starting choice.

Catchment

An area whose rainfall drains toward a shared river or outlet. The selected map cells in this experiment are not a mapped catchment.

Discharge

The volume of water passing a river location per unit of time. The paper's simulated outflow has not been checked against measured river discharge.

Continuity jump

A mismatch between the water storage passed from the previous day and the storage predicted at the start of the next day.

Water budget

An accounting check that water entering minus water leaving matches the change in stored water.

Analytic solution

An answer calculated directly from a mathematical expression. It provides a comparison for the neural solver without proving the model's assumptions are correct.

Surface elevation

The mapped height of the surface, which can include trees and buildings. It is not the same as a survey of the riverbed.

Upwind method

A numerical method that updates a location using information from the side the modelled flow comes from. Here it provides a conventional solver comparison.