The stormwater management model SWMM of the US EPA is widely used to analyse, design or optimise urban drainage systems. To perform advanced analysis and visualisations of model data this technical note introduces the R package swmmr. It contains functions to read and write SWMM files, initiate simulations from the R console and to convert SWMM model files to and from GIS data. Additionally, model data can be transformed to produce high quality visualisations. In accordance with SWMM’s open source policy the package can be obtained through github.com or the Comprehensive R Archive Network (CRAN).

Stormwater runoff quality was measured with online turbidity sensors at four common types of small urban subcatchments: (i) a flat roof; (ii) a parking lot; (iii) a residential catchment; and (iv) a high-traffic street. Samples were taken to estimate site-specific correlations between total suspended solids (TSS) and turbidity. Continuous TSS time series were derived from online turbidity measurements and were used to estimate event loads and event mean concentrations. Rainfall runoff event characteristics were subjected to correlation analysis to TSS loads. Significant correlations were found for rainfall intensities at sites with high imperviousness and decrease with increasing catchment size. Antecedent dry weather periods are only correlated at the parking lot site. Intra-event TSS load distributions were studied with M (V)-curves. M (V)-curves are grouped at runoff quantiles and statistically described with boxplots. All sites show, in general, a more pronounced first-flush effect. While wash-off of the flat roof tends to be source-limited, the parking lot and high-traffic street sites show a more transport-limited behavior. Wash-off process of the residential catchment appears to be influenced by a composition of different subcatchments.

Results from a long-term stormwater quality monitoring program were used to derive total suspended solids (TSS) event load distributions at four small urban environments (flat roof, parking lot, residential catchment, high traffic street). Theoretical distribution functions were fitted to the empirical distribution functions obtained. Parameters of the theoretical distribution functions were optimized with respect to a likelihood function to get both optimized parameters and standard errors. Kolmogorov-Smirnov and Anderson-Darling test statistics were applied to assess the goodness-of-fit between empirical and theoretical distribution. The lognormal distribution function was found to be most expressive to approximate empirical TSS event load distributions at all sites. However, the goodness-of-fit of the statistical model strongly depends on the number of events available. Based on the results of a Monte-Carlo-based resampling strategy, around 40 events should be considered.

Stormwater quality models are usually calibrated using observed pollutographs. As current models still rely on simplified model concepts for pollutant accumulation and wash-off, calibration results for continuous pollutant concentrations are highly uncertain. In this paper, we introduce an innovative calibration approach based on total suspended solids (TSS) event load distribution. The approach is applied on stormwater quality models for a flat roof and a parking lot for which reliable distributions are available. Exponential functions are employed for both TSS buildup and wash-off. Model parameters are calibrated by means of an evolutionary algorithm to minimize the distance between a parameterized lognormal distribution function and the cumulated distribution of simulated TSS event loads. Since TSS event load characteristics are probabilistically considered, the approach especially respects the stochasticity of TSS buildup and wash-off and, therefore, improves conventional stormwater quality calibration concepts. The results show that both experimental models were calibrated with high goodness-of-fit (Kolmogorov–Smirnov test statistic: 0.05). However, it is shown that events with high TSS event loads (>0.8 percentile) are generally underestimated. While this leads to a relative deviation of −28% of total TSS loads for the parking lot, the error is compensated for the flat roof (+5%). Calibrated model parameters generally tend to generate wash-off proportional to runoff, which is indicated by mass-volume curves. The approach itself is, in general, applicable and creates a new opportunity to calibrate stormwater quality models especially when calibration data is limited.