Seven Minutes to the Border
A catastrophic flood tore down the Nepal — China frontier last August. The official story was a glacial lake bursting its banks. The timing says something else entirely happened.
At 08:37 in the morning on 26 August 2026, seismometers across the Himalaya registered a shock near the Nepal — China border. Within minutes, a wall of water was moving down the Bhote Koshi — Trishuli river system. It destroyed Timure and swept through Rasuwagadhi, Syabrubesi and Betrabati. Damage was reported as far downstream as Galchhi and Malekhu, more than a hundred kilometres away. Hundreds of people were killed or went missing on both sides of the border.
Several of the river gauging stations that would normally record such an event were themselves destroyed by it. That absence is why this reconstruction leans on satellites, seismographs and physics rather than instrument records.
The early explanation, carried by most outlets in the days afterward, was a glacial lake outburst flood: a meltwater lake perched on the Purepu Glacier had drained catastrophically. It was a reasonable first guess. That exact lake had done something similar, on a smaller scale, on 8 July 2025.
But a single number in the record does not fit that story, and it turns out to be the number that matters most.
Nepal time
reached the border
in those minutes
185 km/h
A surveillance camera at the Gyirong Port border crossing recorded floodwaters arriving roughly seven minutes after the seismic signal. The distance from the probable source area to that crossing is about 21.6 kilometres. That gives an average front speed of just over 51 metres per second.
Rivers do not do this. A flood wave moving down a channel is governed by friction against the bed and banks, and that physics imposes a ceiling. To find out where the ceiling sits for this valley, we built a hydraulic model of it — and the model, calibrated against the flood's own arrival time at a town 112 km downstream, moves its flood wave at about 12.6 metres per second.
The observed front was roughly four times faster than the physics of river flooding allows.
Whatever reached the border in seven minutes was not a flood wave travelling down a river. It was something else, moving under different rules.
The investigationWhat we actually did
Answering this required assembling the event from scratch: locating the source, sizing the mass that failed, working out what it built in the valley, and modelling what came downstream. Each step used a different, independent line of evidence, which matters — because when independent methods converge on the same answer, that convergence is itself the result.
Data and tools
Everything here is free and publicly accessible. Nothing proprietary, no commercial imagery, no licensed modelling software.
- Sentinel-2 L2A optical imagery (Element84 Earth Search, AWS) — source-area survey, damage mapping, flood-width measurement.
- Sentinel-1 SAR — attempted cloud-penetrating confirmation of the source.
- Landsat Collection 2 L2 via Microsoft Planetary Computer — gap-filling where Sentinel-2 was cloud-blocked. (USGS's own archive is requester-pays and was unavailable.)
- Copernicus GLO-30 30 m global DEM — all terrain geometry: river centerline, valley cross-sections, dam-height solve, inundation mapping.
- Public seismic catalogues (USGS, GFZ) — event timing and magnitude.
- Python throughout:
rasterio,pysheds,pyproj,scikit-image,numpy. The hydraulic solver was written for this problem rather than taken off the shelf.
Six experiments, in sequence
The reconstruction ran as a chain, each step feeding the next. Where a step failed to confirm something, that failure is reported here as a result rather than quietly dropped — several of them constrain the answer as tightly as the successes do.
Map the damage from orbit
Difference vegetation index (NDVI) between pre- and post-event imagery along the entire 124 km river corridor, to find where flood damage begins and ends.
Hunt for the source lake
Water index (NDWI) scanning across the basin's glacial lakes; snow index (NDSI) change detection at the candidate failure slope, with a control site to test whether the signal was real.
Weigh the landslide with a seismograph
Convert the recorded seismic magnitude to a mass using the Ekstr • m & Stark (2013) force-scaling relation, validated against their own published worked example.
Solve the dam from the terrain
Use real DEM cross-sections to find where that mass, dropped into this valley, produces a physically self-consistent blockage — and how much water would pond behind it.
Model the breach and the flood
Published breach-parameter regressions to set the outflow hydrograph, then a 1D Saint-Venant solver to route it 106 km downstream, calibrated against the real reported arrival time.
Test the model against the satellite
Compare modelled inundation width to Sentinel-2-measured damage width at 36 cross-sections — the independent check that the whole chain has to survive.
Experiment 01 & 02Where did it start?
The received account named a specific glacial lake. Testing that meant looking at the lake itself, and at the corridor of destruction leading away from it.
The lake is real. NDWI scanning confirmed the Purepu supraglacial lake system at roughly 28.404 • N, 85.59 • E — persistent, mapped, and genuinely the source of the July 2025 event. What no available satellite pass could establish was that it drained on 26 August 2026. Cloud cover over the source region in the only usable post-event optical pass ran to 75 — 79%, and the only Sentinel-1 radar pass in range predates the event by two days.
So we approached from the other direction: instead of asking whether a particular lake emptied, we asked how far up the valley the destruction actually reaches. Flood damage strips vegetation, and that leaves a clear signature in NDVI.
A negative result worth reporting
At the candidate failure slope on Langtang Lirung's north face, snow-index change detection found 1.86 km • of apparent snow and ice loss between an April baseline and the post-event pass — which looks like confirmation until you run a control.
We tested an unrelated, undisturbed slope nearby using the identical method. It showed a larger change fraction: 33.8% against 19.4% at the candidate site. Ordinary April-to-August seasonal melt cannot be separated from avalanche disturbance this way. The test did not confirm the source, and reporting it as though it had would have been wrong.
The honest conclusion from remote sensing alone: no specific origin point was visually confirmed. The convergent evidence brackets the origin to roughly the 15 — 25 km reach — consistent with, but not proof of, the Langtang Lirung north-slope location that glaciologists had named within days of the event.
That bracket was enough to proceed. The next question was not where but how much.
Experiment 03Weighing a mountain with a seismograph
A large landslide announces itself seismically. As an enormous mass accelerates downslope and then decelerates against the valley floor, it applies a long-period force to the Earth that seismometers record — and which, crucially, looks different from the sharp rupture of a tectonic earthquake.
Ekstr • m and Stark established the quantitative link between that signal and the mass responsible in a 2013 paper in Science. Their relations connect the peak force to a landslide magnitude, and the force to the mass. Applied here, they turn a seismic reading into an estimate of how much mountain came down.
Before trusting it, we checked it against the paper's own worked example, the 2010 Hunza-Attabad landslide in Pakistan: the formula returns 3.24 • 1011 kg where field and inversion estimates give 1.1 — 1.4 • 1011 kg. That is inside the factor-of-two uncertainty the authors state themselves — so the method is usable, provided nobody pretends the output is precise.
The magnitude will not sit still
Here the record becomes genuinely difficult. The event was first catalogued as a magnitude 4.4 earthquake. It was later reclassified by USGS as a magnitude 5.2 landslide, once analysts determined the shaking came from mass movement rather than fault slip.
A revision of 0.8 magnitude units sounds minor. Run through the scaling relation, it changes the implied mass roughly sixty-fold.
The likely explanation is not a record-breaking landslide but a mismatch of measurement scales. The Ekstr • m — Stark relation is calibrated against their own catalogue of events analysed by a specific long-period single-force inversion. A routine catalogue reclassification from "earthquake" to "landslide" comes from a different, more generic waveform assessment. There is no guarantee the two numbers mean the same thing, and treating them as equivalent would be exactly the kind of error this reconstruction tried to avoid elsewhere.
This work therefore proceeds on the original M4.4 — 4.5 reading, giving a slide volume near 1.3 — 2.2 million cubic metres, and flags the M5.2 figure as a real and unresolved open question rather than silently adopting or discarding it.
Experiment 04What that mass built in the valley
Knowing roughly how much material fell is not the same as knowing what it did on landing. A landslide dam's height depends entirely on the shape of the valley it lands in — the same volume that barely dents a wide floodplain can plug a narrow gorge completely.
So rather than assume a shape, we solved for one. At each candidate location along the bracketed origin reach, the DEM supplies a real cross-section; from that we build an actual area-versus-height relationship for that specific spot, and solve for the dam height at which the local cross-sectional area, filled to that height over its own width, accounts for the estimated slide volume. No geometry is imported from anywhere else in the valley.
Why the constraint matters
An early version of this calculation used a river width measured from satellite imagery in the wide, low-gradient reach 80 km downstream, where the valley floor sits between 450 and 1,200 m elevation. Applied to a dam site in a narrow gorge near 3,850 m, that width is simply the wrong geometry — it describes a different landscape.
Forcing the solve to use only each site's own local cross-section changed the answer materially, and is the reason the result below is trustworthy at all. Tested across candidate sites from 16 to 24 km, only one produced a physically sensible, non-degenerate solution.
| Quantity | Solved value |
|---|---|
| Dam location | Chainage 18 km |
| Dam height | 46.9 m |
| Pool surface elevation | 3,895.5 m |
| Backwater extent upstream | 1.0 km |
| Impounded volume | 1.11 million m • |
An independent check falls into place
That dam site was derived purely from seismology and terrain geometry. It knows nothing about photographs of the mountain.
Separately, satellite-imagery review published in a hazard assessment put the avalanche detachment at roughly 5,200 m elevation, dropping about 1,200 vertical metres to the valley floor. Commentary circulating on social media described the fall as around 1.5 km.
The bed elevation at our independently solved dam site is 3,848.6 m. The drop from 5,200 m to that point is 1,351 metres — sitting almost exactly between the two reported figures. Two completely separate chains of evidence, one seismic and topographic, one photographic, landed in the same place.
The fall height also passes a physics check. For a drop of 1,200 — 1,500 m, frictionless free fall would cap velocity around 150 — 170 m/s. The independently reported descent speed of roughly 54 m/s is about a third of that ceiling — squarely inside the 20 — 60% efficiency that real rock and ice avalanches achieve once friction and entrainment are accounted for.
Experiment 05Modelling the flood that followed
With a dam height and an impounded volume, the next question is what happens when it fails. This needs two things: a description of the breach, and a way to route the resulting flood downstream.
The breach
Rather than prescribe a plausible-looking hydrograph shape, we used published regression equations from Froehlich (1995), obtained through a US Bureau of Reclamation review by Wahl (2001). That source was chosen deliberately: it does not merely cite competing breach formulas, it tests ten of them against a database of 108 real historical dam failures and reports which perform best. Froehlich's peak-outflow relation came out with the smallest demonstrated uncertainty of any method examined.
| Breach parameter | Relation | Result |
|---|---|---|
| Peak outflow | 0.607 • Vw0.295 • hw1.24 | 4,353 m • /s |
| Average breach width | 0.1803 • Ko • Vw0.32 • hb0.19 | 45.1 m |
| Formation time | 0.00254 • Vw0.53 • hb — 0.9 | 7.64 min |
The 45 m breach width sits in the same range as the surveyed 27 • 9 m breach from the 2016 Zhangzangbo outburst in this same basin — larger here, consistent with a taller dam. That analog event, modelled in detail by other researchers, is the single most useful calibration reference available for this valley.
Routing it downstream
The hydraulic model
- Governing equations: full 1D Saint-Venant (dynamic wave) — solving momentum, not a simplified hydrologic routing approximation.
- Numerical scheme: explicit Rusanov (local Lax-Friedrichs) finite volume with a surface-gradient source term, chosen for shock-capturing stability on a dam-break front. An earlier MacCormack implementation oscillated badly at the shock and was abandoned.
- Geometry: centerline traced by D8 flow accumulation, which is monotonically downhill by construction — zero uphill steps across 124 km. A first attempt using least-cost pathfinding required extensive despiking and still wandered out of the valley in places.
- Cross-sections: 249 sections sampled from the DEM at 500 m spacing, idealised to trapezoids in the upper gorge where 30 m posting cannot resolve a channel maybe 10 — 50 m wide.
- Calibration: Manning's roughness tuned against the one hard timing constraint available — the flood was reported at Galchhi, 112 km downstream, within 134 minutes.
Experiment 06The test the model nearly failed
A flood model that gets the timing right can still be wrong about everything else. The independent check is width: how far across the valley did the water actually spread? Sentinel-2 imagery answers that directly, and we measured it at 36 cross-sections through the Betrabati — Malekhu reach.
Run as clear water, the model predicted a median inundation width of 52 m. The satellite says 385 m. The model was seven to eight times too narrow.
That gap is not a rounding error, and no adjustment to the peak discharge fixes it. What it points to is that this was never clear water. Every account describes an ice-and-rock avalanche; the flow that resulted carried an enormous sediment load, and a debris-laden flow behaves differently from water in two specific ways — it is bulkier, and it is far rougher.
A literature-bounded debris correction
Two adjustments, both taken from published ranges rather than fitted to close the gap:
- Volumetric bulking at a solids concentration Cv = 0.65, near the top of Pierson's (2005) debris-flow classification (0.4 — 0.8) — a 2.86 • discharge multiplier.
- Roughness raised to Manning's n = 0.15, the ceiling of the 0.05 — 0.15 range reported for debris flows by Rickenmann (1999) and Costa (1988).
Roughness turns out to be the more effective lever on width. On a steep, confined channel, extra discharge at fixed roughness mostly becomes extra velocity; it is extra roughness that forces the flow into extra depth and lateral spread — which is also what the grain-collision physics of a real debris flow does.
Both parameters were pushed to the top of their published ranges, and the result checked afterward. Median mismatch improved from 7.3 • to 3.0 • — a real gain, but still short.
The remaining factor of three had a more interesting cause than "the physics is wrong". Looking at individual cross-sections rather than the median, the mismatch ranged from 0.17 • to 18 • . At some points the model was already too wide. The worst failures clustered precisely where the observed damage was widest — which is exactly where the measurement method itself breaks down.
Our width extraction walks two straight rays outward from the channel centerline until the terrain rises above the modelled water surface. Where floodwater spread laterally around a low rise into a connected floodplain, those two rays stop at the rise and never find it.
So we re-measured the identical model output a second way: an eight-connected flood fill, taking any terrain topologically connected below the water surface. That method has the opposite failing — it has no concept of momentum or travel time and will happily fill any low ground it can reach.
The answerAn avalanche pushing water, not a dam filling up
Return to the seven minutes.
The instinct that a landslide dam must fill before it bursts is sound physics, and we tested it directly. Filling the 1.11 million m • impoundment at plausible monsoon inflow rates takes somewhere between two and thirty-one hours — around eight hours at a mid-range estimate. That is the right description of a blockage that forms and then slowly ponds up behind itself. It cannot describe water arriving 21.6 km downstream seven minutes after the slope failed.
Three independent lines of evidence resolve it.
Speed. At 51 m/s, the front was moving four times faster than friction-governed channel flow permits. An impulse wave — water displaced by a mass entering the channel — is not bound by that limit, because it is not being routed, it is being shoved.
Appearance. The border footage shows what looks like water, not a visible wall of rock and ice. That fits a leading edge made largely of the river's own pre-existing water, pushed ahead of the slower, heavier debris front, which drags against the bed in a way the water surge does not.
No filling time was needed. This closes it. At the M4.4-derived mass, an ice-rock avalanche that is even 20 — 30% glacier ice by mass carries 0.65 — 1 million m • of water-equivalent; at 50 — 100% ice it reaches 1.6 — 3.2 million m • . Those figures are comparable to, or larger than, the entire impoundment volume solved from the dam geometry. The water did not have to accumulate from the river. Most of it was already sitting frozen on the mountain, and it arrived as ice.
This was not a lake that filled slowly and let go. It was an avalanche that hit a river like a piston — driving the water already in the channel downstream ahead of it, melting as it went — and only afterwards, once the debris settled into a genuine blockage, produced anything resembling a conventional dam-break flood.
The two processes are separated in time, and the confusion between them is what made the original account seem plausible.
The slope fails
An estimated 1.3 — 2.2 million m • of rock and glacier ice detaches around 5,200 m and falls roughly 1,350 vertical metres toward the valley floor.
Impact surge
The mass strikes the channel and displaces the river's water ahead of it. Governed by impulse-wave dynamics rather than channel friction, this front moves near 51 m/s. Order-of-magnitude equivalent discharge: 3,800 — 38,000 m • /s, lasting seconds.
The border camera
Floodwater reaches Gyirong Port, 21.6 km downstream. It looks like water because it largely is water — displaced river flow and meltwater, running ahead of the debris.
Debris front, then impoundment
The heavier ice and rock arrives behind the surge, decelerating against the bed, and settles into a blockage roughly 47 m high that ponds about 1.11 million m • .
The breach flood
The blockage fails: peak outflow near 4,353 m • /s, breach forming over 7.6 minutes with a very sudden recession. This is the classical dam-break flood — and the process the hydraulic model represents.
Galchhi, 112 km downstream
The sustained flood arrives at a far more ordinary 12.6 m/s, 125 minutes after the trigger — matching the reported constraint of 134 minutes.
A separate barrier lake
Chinese authorities identify a newly formed blockage lake, which bursts a full day later. This one is the slow accumulate-then-breach process — a distinct event, and a useful contrast.
That last entry matters. Nobody is claiming landslide-dammed lakes do not form and slowly breach in this basin — one did, forty-eight hours later, in the textbook manner. Having a clean example of the slow process in the same valley days afterward is what makes the case that 26 August was something different, rather than a quarrel about definitions.
LimitsWhat this reconstruction does not know
Ordered by how much each would change the picture if resolved.
- The source is not visually confirmed. The 18 km dam site is the best-supported candidate from seismic and terrain evidence, corroborated by reported fall height — but no clear post-event image of the failure scar was available.
- The magnitude question is open. M4.4 versus M5.2 is a sixty-fold difference in implied mass. This work uses the lower figure and explains why, but that is an argument, not a resolution.
- A 30 m DEM cannot resolve the upper gorge. Where the real channel may be 10 — 50 m wide, idealised trapezoids stand in for unresolvable cross-sections.
- No bathymetry exists for any lake or impoundment here. Every volume is inferred from remote geometry, never surveyed.
- The debris correction is an approximation. Bulking plus elevated roughness is a planning-level substitute for genuine two-phase rheology. A dedicated debris-flow model — RAMMS, FLO-2D or similar, as was applied to the 2016 analog — is the rigorous next step.
- This is a one-dimensional model. Inundation width is draped onto terrain around a 1D water surface. Standard practice, but not a true 2D hydraulic solution, and Figure 6 is where that limitation becomes visible.
- Every discharge figure inherits its method's uncertainty. Froehlich's relations carry 0.3 to 1 log-cycle of scatter even as the best-performing option available.