Dead Run watershed at Franklintown. Small enough for storm-core placement and drainage geometry to matter sharply.
Radar rainfall resolution, bias corrected using a network of 54 rain gauges around Baltimore.
Three-hour extreme-storm catalog drawn from a 16-year radar record and resampled through SST.
Paper's headline average increase in flood peaks when realistic spatial rainfall heterogeneity is retained instead of forcing uniform rainfall.
real space-time structure
resample + transpose
distributed runoff + routing
SST → distributed hydrology
1. Build the storm catalog
Select the 200 largest 3-hour “Dead-Run-shaped” storms from 2000–2015 radar fields over a ~7000 km² transposition domain.
2. Generate synthetic years
Draw a stochastic number of storms per year from a Poisson model, then translate storm fields within the regional domain using a nonuniform occurrence distribution.
3. Keep annual maxima
For each synthetic year, retain the largest basin-average 3-hour rainfall. Repeat to create frequency curves.
4. Preserve the fields
The actual transposed space-time rainfall field, not merely its basin average, becomes hydrologic model input.
5. Route through GSSHA
Run overland flow, channels, storm sewers, infiltration, and detention through the distributed watershed model.
6. Compare tails
Ask how flood distributions change with return period, basin, rainfall structure, and a uniform-rainfall counterfactual.
Why SST can estimate rarer events than the local record
SST substitutes regional space for local time. A storm observed elsewhere in a meteorologically coherent domain becomes a plausible candidate over the target basin after transposition. This is powerful but assumption-heavy: transposition must not violate regional storm climatology.
Dead Run subwatersheds
| Basin | Area km² | Impervious % | Detention-controlled % |
|---|---|---|---|
| DR1 | 1.32 | 73.6 | 41.9 |
| DR2 | 1.92 | 55.5 | 18.5 |
| DR3 | 4.95 | 62.2 | 24.4 |
| DR4 | 6.29 | 51.5 | 12.2 |
| DR5 | 2.05 | 47.9 | 3.2 |
| Franklintown | 14.3 | 52.3 | 25.1 |
Rainfall feature explorer
M(t)
Basin-average rainfall rate at time t. Conventional magnitude descriptor.
Mmax
Maximum basin-average rain rate during the event.
Rsum
Storm-total basin-average rainfall depth.
Z(t)
Fraction of basin covered by a storm core above 25 mm h⁻¹.
RWD
Rainfall-weighted flow distance. Rain far upstream along long drainage paths differs from rain near the outlet.
S(t)
Dispersion of rainfall-weighted flow distance; a compact descriptor of spatial organization.
Uniform-rainfall counterfactual
Interpretation
Uniform rainfall erases storm-core placement and some tributary synchronization. The lost geometry changes peak discharge even when basin-scale rainfall magnitude is comparable.
Rainfall return period ≠ flood return period
Spearman rank correlation between rainfall and modeled flood return periods. Related, but far from one-to-one.
Outlet validation
Median Nash–Sutcliffe efficiency at Franklintown for 21 warm-season validation events.
Spatial heterogeneity
Largest tabled uniform-versus-distributed median peak difference: DR5 at the 50-year rainfall return period.
The paper as an argument under review
Interactive discussion begins.
Referee 1 issues a strongly skeptical review: novelty, short records, model errors, generalizability, and paper organization.
Referee 2 requests major revisions focused on statistical interpretation, physical explanation, return-period scale, drainage representation, equations, and validation.
Detailed author replies uploaded.
Editor: publish subject to revisions and further review.
Revised manuscript and tracked changes submitted.
Referee 2: comments satisfactorily addressed; only the missing Δt in Eq. 3 remains.
Referee 1 recommends publication; editor requests technical corrections.
Final paper published.
Challenge matrix
Critical audit
Strongest evidence
The distributed-versus-uniform rainfall experiment directly tests a major design-storm assumption and finds large peak-flow effects across several basins and return periods.
Most assumption-heavy step
SST extrapolates a 16-year radar archive into long synthetic frequency estimates through regional space-for-time substitution.
Hydrologic weakness
Validation is uneven. Franklintown is respectable, while DR1 has poor NSE and roughly +57% median peak error, probably related partly to incomplete sewer mapping.
Residual paper issue
Equation 3 remains dimensionally incomplete as printed without Δt, despite the second reviewer explicitly flagging it again after revision.
Statistical caution
Random-forest feature-importance shifts are physically interesting but moderate, and correlated predictors can complicate feature-importance interpretation.
Transferable lesson
For quick-response urban basins, precipitation validation should preserve or test storm geometry, not merely basin totals.