Dynamic Distributed Storage of Stormwater in Sponge-Like Porous Bodies: Modelling Water Uptake

: An innovative concept of dynamic stormwater storage in sponge-like porous bodies (SPBs) is presented and modelled using ﬁrst principles, for down-ﬂow and up-ﬂow variants of SPBs. The rate of inﬂow driven by absorption and / or capillary action into various porous material structures was computed as a function of time and found to be critically dependent on the type of structure and the porous material used. In a case study, the rates of inﬂow and storage ﬁlling were modelled for various conditions and found to match, or exceed, the rates of rainwater inﬂow and volume accumulation associated with two types of Swedish rainfalls, of 60-min duration and a return period of 10 years. Hence, the mathematical models indicated that the SPB devices studied could capture relevant amounts of water. The theoretical study also showed that the SPB concepts could be further optimized. Such ﬁndings conﬁrmed the potential of dynamic SPB storage to control stormwater runo ﬀ and serve as one of numerous elements contributing to restoration of pre-urban hydrology in urban catchments. Finally, the issues to be considered in bringing this theoretical concept to a higher Technological Readiness Level were discussed brieﬂy, including operational challenges. However, it should be noted that a proper analysis of such issues requires a separate study building on the current presentation of theoretical concepts.


Introduction
Urbanization dramatically alters the hydrological cycle of developing areas by reducing hydrological abstractions and accelerating runoff, which leads to increased runoff volumes and flow peaks, and ultimately the risk of water ponding or flooding [1]. Such changes of the urban environment, accompanied by deteriorating runoff quality and geomorphology of urban streams, lead to reduced biodiversity and contribute to unsustainability in urban areas [2]. Therefore, restoration of the pre-development catchment by measures enhancing water abstraction (e.g., infiltration in soakaways and evapotranspiration on green roofs) and slowing the speed of runoff (e.g., cross-berms in runoff swales) have been pursued in modern drainage design [3,4]. In practice, the above measures, or their functional principles, are applied via broadly varying scales, typically described as lot-level [5], neighborhood and catchment scales [6]. The discussion here focuses on lot-scale measures (LSMs), which in common terminology are also called distributed control measures, and when such LSMs serve to dissipate runoff volume they are referred to as runoff source controls. The importance of these measures follows from their key features: The objectives of the paper are to: (a) describe the theoretical sponge-like porous bodies (SPB) water storage concept on the basis of first principles, (b) derive and verify semi-analytical models for entry of stormwater into SPBs, (c) use such models to demonstrate the theoretical capacity of SPBs to fully capture selected Swedish design rainfalls, and (d) discuss practical challenges that may be encountered in field applications in stormwater control, including steps in advancing technology readiness level (TRL) of SPBs.

SPB Storage Concept Description
In the first SPB storage variant, Down-flow SPB storage, a relatively large area, such as sections of a roof, parking lot, playground, or football field, would be covered with a material that absorbs the rainwater directly upon contact. The material then swells in a vertical direction, keeping the intercepted water in place (Figure 1a). Note that the SPB storage sketches presented in the figure are intended just to elucidate the theoretical concept of such storage, without any aesthetic, practical or placement considerations.
One possible type of material to consider in this application are hydrogels, recognizing that such materials can satisfy several demands of Down-flow SPB storage. For maximum effectiveness, the influx of water into a hydrogel layer should match the influx of rainwater, but partial interceptions One possible type of material to consider in this application are hydrogels, recognizing that such materials can satisfy several demands of Down-flow SPB storage. For maximum effectiveness, the influx of water into a hydrogel layer should match the influx of rainwater, but partial interceptions with some water bypassing are fully acceptable, as demonstrated by similar widely used control measures, e.g., green roofs [11]. In fact, there may be opportunities for installing Down-flow SPB storage units in conjunction with green or gray roofs.
In the second variant, water ponded on the ground moves upwards and is stored in two types of pre-installed Up-flow SPB storage structures: (i) intermediate height (1-2 m) supported vertical structures that expand only in the horizontal direction when absorbing water (Figure 1b), and (ii) low-height (<0.2 m) unsupported structures that grow from the ground as they absorb water (Figure 1c). The former structures can be wrapped around trees (Figure 1b) or lamp posts, or attached to walls and concrete bridge abutments, to mention a few examples. The supported storage structures may swell in a horizontal direction when water first moves horizontally in and then upwards through the structures. For the unsupported structures that grow from the ground as they absorb water (Figure 1c), the material not only swells, but also stiffens. This process is analogous to the rising of flowers, when irrigated after a dry period. The water impoundments formed by these structures may also store water from adjacent areas. Examples of these storage structures include hydrogel fibers and porous rods made of natural fibers [31].
The design of unsupported storage structures may require the use of a cluster of compliant rods, which stiffen as they become saturated by capillary action, or by applying a pressure gradient that makes them grow vertically as a function of the inflow rate. However, there must be spaces for storing the expanding material during dry weather periods, which can be a challenge with respect to space limitations and safety. Depending on site conditions, individual rods may grow in isolation, one by one, or in clusters forming various patterns serving as barriers around water impoundments holding significant volumes of stormwater ( Figure 1c). Such impoundments could be designed to feed soakaway pits, or similar infiltration structures.

Governing Equations
The governing equations are presented below for the Down-flow and Up-flow SPB storage variants.

Down-Flow SPB storage: Governing Equations
The flow into absorbing materials, like hydrogels, may be described by the diffusion equation: where θis the concentration of water within the swelling material at location rand time t and D ij is the diffusion coefficient. Assuming that D is constant in space and independent of θ, Equation (1) reduces to: where the product Dθ ,j is the flux of water per unit area. From these equations the volumetric uptake V of water per unit time ∆t, Q ab , was derived by Sweijen et al. [32] as or, in a different form, after expressing the volumetric water uptake as a function of time ∆t, the following equation is obtained where A ab is the area of the interface between the water and the absorbing material, n is the vector normal to the interface, and V ab is the volume of water that penetrated through the surface A ab into the absorbing material. From this equation it is obvious that the flow rate will increase with the size (area) of the interface, the magnitude of D, and the gradient of θ(r,t), which should be as large as possible. The governing equation for the up-flow SPB storage is assumed to be the Darcy law, together with the condition of water incompressibility, according to: where u i is the velocity, φ the porosity, K ij the permeability and p the pressure. The constant g is the acceleration due to gravity that is set to act in the negative x 1 direction and ρ is the density of water. This equation is valid for Newtonian flow through a stationary porous medium, up to particle Reynolds numbers Re p = Ud/ν~40 [33], where U is the average velocity in the porous medium, d is a characteristic length of the solids in the porous medium, and ν is the kinematic viscosity of the fluid. The porous medium is not allowed to deform in this first model of Up-flow SPB storage. In our model, the flow is driven by capillary action described by where ∆p is the pressure jump over the curved water surface (the capillary pressure), γ is the surface tension, Θ is the contact angle, and R is the pore radius of the porous medium. By replacing ∆p on the left-hand side of Equation (7) with ρ gh max , the maximum height, h max , within a capillary tube can be derived, and this expression is also known as Jurin's law. Important for Up-flow SPB storage is that the main driving mechanism for the upward flow, capillary action, is promoted by small scale capillaries, while the water moves much more easily through large-scale channels. The reason for this is that the capillary pressure increases with 1/R, while the permeability K ij and thus the flow rate is proportional to R 2 , where R is the typical length dimension of the problem, e.g., a pore radius in a porous medium or the size of particles of the porous medium. In the following, we have chosen to use the size of the particles (fibers), while it is also possible to include the porosity in the expression [27,28].

SPB Process and Geometric Parameters
As briefly mentioned in the previous section, there are several parameters of SPBs that can be tuned up for optimal performance. Relevant values of the most significant among these parameters are discussed here in more detail.

Down-Flow SPB Storage: Parameters
The magnitude of the diffusion coefficient of water entering the Down-flow SPB storage, D, is determined by the properties of the absorbing material, the composition of water, and conditions at the site. It is, therefore, desirable to find a combination of such characteristics that maximizes D. One group of efficient absorbers of water are hydrogels. For example, Bajpai [34] derived values of D between 6.9 × 10 −9 and 1.3 × 10 −8 m 2 /s for a hydrogel based on Acrylamide doped with different concentrations of maleic acid. Doll et al. [35] measured values between 1.3-1.4 × 10 −9 m 2 /s for two types of bio-based hydrogels, and El-Hamshary [36] tested a number of variants of poly(acrylamide-co-itaconic acid) hydrogels and obtained values between 3.0 × 10 −10 -1.2 × 10 −8 m 2 /s. Hence, there is a broad range of D values in the literature, depending on the hydrogel composition. Following Bajpai [34], D is here approximated as 1.0 × 10 −8 m 2 /s, which falls into the upper end zone of the results reviewed but is short of the maximum.
Regardless of D, the rate of water uptake is enhanced by a large area of contact between the absorbing material and the water, A ab , as shown in Equation (3). So instead of using a flat surface Water 2020, 12, 2080 6 of 21 area of the absorbing material the interface with the water should be corrugated in some way. As an example, let us study a prismatic block of absorbing material, with a plan view area, A flat , height (thickness) H, and a grid arrangement of relatively small cavities, with a square cross-section b × b in plan, and depth h, shown in Figure 2. In such a simplified porous medium body, repeating square cells can be described by a plan dimension (side B), and height H, in Figure 2. Now also let b = b/B and h = h/B and observe that 0 < b < 1 and 0 < h < H/B. The total contact area between the water and the absorbing material may now be derived as: and for a maximum value of h , Equation (8) becomes The absorbing area increases with b but, as stated above, the theoretical maximum value of b equals one. In practice, the actual value should be smaller, around 0.5, to maintain a sufficient volume of the absorbing material. More interesting is the ratio H/B. Realistically, the upper limit H is constrained, while B is fairly arbitrary and sets the scale of the cavities. Finally, Equation (9) may be generalized to: where C is a constant equal to 4 for the geometry in Figure 2 and π for cylindrical cavities, for instance. As shown by Equations (3) and (4) the flow rate into the absorbing material is directly related to the gradient of θ(r,t), which will decrease as the absorbing material becomes saturated.
Water 2020, 12, x FOR PEER REVIEW 6 of 22 Regardless of D, the rate of water uptake is enhanced by a large area of contact between the absorbing material and the water, Aab, as shown in Equation (3). So instead of using a flat surface area of the absorbing material the interface with the water should be corrugated in some way. As an example, let us study a prismatic block of absorbing material, with a plan view area, Aflat, height (thickness) H, and a grid arrangement of relatively small cavities, with a square cross-section b × b in plan, and depth h, shown in Figure 2. In such a simplified porous medium body, repeating square cells can be described by a plan dimension (side B), and height H, in Figure 2. Now also let b' = b/B and h' = h/B and observe that 0 < b' < 1 and 0 < h'< H/B. The total contact area between the water and the absorbing material may now be derived as: and for a maximum value of h', Equation (8) becomes The absorbing area increases with b' but, as stated above, the theoretical maximum value of b' equals one. In practice, the actual value should be smaller, around 0.5, to maintain a sufficient volume of the absorbing material. More interesting is the ratio H/B. Realistically, the upper limit H is constrained, while B is fairly arbitrary and sets the scale of the cavities. Finally, Equation (9) may be generalized to: where C is a constant equal to 4 for the geometry in Figure 2 and π for cylindrical cavities, for instance. As shown by Equations (3) and (4) the flow rate into the absorbing material is directly related to the gradient of θ(r,t), which will decrease as the absorbing material becomes saturated. In Section 2.6, the sensitivity of the above parameters will be discussed with the combinations of values listed in Table 1 as a base and with h = H. In Table 1 it can be seen that, as the unit cell length scale B decreases, the number of square cavities increases for the same plan-view area. Hence, the area exposed to water, Aab, increases substantially. In Section 2.6, the sensitivity of the above parameters will be discussed with the combinations of values listed in Table 1 as a base and with h = H. In Table 1 it can be seen that, as the unit cell length scale B decreases, the number of square cavities increases for the same plan-view area. Hence, the area exposed to water, A ab , increases substantially. Up-flow SPB storage is schematized as a model cylinder consisting of a solid cylinder with radius δ and height H s , a porous inner annulus with a larger radius a and height H i , and a porous outer annulus with a larger radius b and a height H o (see Figure 3). The solid cylinder mimics a tree or lamp post, and the inner and outer annuli confine the porous media, cf. Figures 1b and 3. The porosities are φ i and φ o , respectively, the geometrical scales of the porous media are R i and R o , respectively, and the origin of a cylindrical co-ordinate system (x, r, θ) is located at the center of the solid cylinder base ( Figure 3). The actual dimensions of the cylinders are arbitrary, and so are the values of R and φ, but to make the cylinder an effective absorber (R i R o ), we set R i < R o . By applying this model, the water uptake, as a function of time, Q in (t), can be derived for various geometries and conditions, and the Up-flow SPBs can be evaluated. Now, as an additional assumption, let the porous media consist of solid vertical rods with radii R i and R o , respectively, let the porosity φ be 0.7 in both annuli and let the temperature be 20 • C implying that the surface tension (water-air) is 72.8 mN/m. To simplify the model further, the radii are set to remain constant during the filling, and the porous media do not deform; hence, there is no swelling and perfect wetting is assumed, θ = 0. The sensitivity of the parameters of Up-flow SPB storage will be discussed in Section 2.6, based on the combinations of the values listed in Table 2. Up-flow SPB storage is schematized as a model cylinder consisting of a solid cylinder with radius δ and height Hs, a porous inner annulus with a larger radius a and height Hi, and a porous outer annulus with a larger radius b and a height Ho (see Figure 3). The solid cylinder mimics a tree or lamp post, and the inner and outer annuli confine the porous media, cf. Figures. 1b and Figure 3. The porosities are ϕi and ϕo, respectively, the geometrical scales of the porous media are Ri and Ro, respectively, and the origin of a cylindrical co-ordinate system (x, r, θ) is located at the center of the solid cylinder base ( Figure 3). The actual dimensions of the cylinders are arbitrary, and so are the values of R and ϕ, but to make the cylinder an effective absorber (Ri ≠ Ro), we set Ri < Ro. By applying this model, the water uptake, as a function of time, Qin(t), can be derived for various geometries and conditions, and the Up-flow SPBs can be evaluated. Now, as an additional assumption, let the porous media consist of solid vertical rods with radii Ri and Ro, respectively, let the porosity ϕ be 0.7 in both annuli and let the temperature be 20 °C implying that the surface tension (water-air) is 72.8 mN/m. To simplify the model further, the radii are set to remain constant during the filling, and the porous media do not deform; hence, there is no swelling and perfect wetting is assumed, θ = 0. The sensitivity of the parameters of Up-flow SPB storage will be discussed in Section 2.6, based on the combinations of the values listed in Table 2.

Solutions of the Governing Equations
The governing equations are solved numerically for the two geometries considered. For the model of the Down-flow SPB, this is rather straightforward, while for the model of the Up-flow SPB some further analysis is required. Most of the derivations for the Up-flow SPB are presented in Supplementary Materials.

Down-Flow SPB Storage: Numerical Set-Up
A numerical scheme for Equations (2)-(4) and (10) is set-up in COMSOL Multiphysics ® for the unit cell defined by dashed lines in Figure 4. Equation (2) is then solved with the boundary conditions specified in Equations (11)- (18). Notice that it is assumed that the height of the absorbing material is constant throughout the absorption process, hence the material will not swell in the model. Instead it is assumed that only 50% of the cell volume is occupied by the absorbing material making room for the water. These restrictions may be relaxed in future studies.
Water 2020, 12, x FOR PEER REVIEW 8 of 22

Solutions of the Governing Equations
The governing equations are solved numerically for the two geometries considered. For the model of the Down-flow SPB, this is rather straightforward, while for the model of the Up-flow SPB some further analysis is required. Most of the derivations for the Up-flow SPB are presented in Appendix A.

Down-Flow SPB Storage: Numerical Set-Up
A numerical scheme for Equations (2)(3)(4) and (10) is set-up in COMSOL Multiphysics ® for the unit cell defined by dashed lines in Figure 4. Equation (2) is then solved with the boundary conditions specified in Equations (11)- (18). Notice that it is assumed that the height of the absorbing material is constant throughout the absorption process, hence the material will not swell in the model. Instead it is assumed that only 50% of the cell volume is occupied by the absorbing material making room for the water. These restrictions may be relaxed in future studies.  The problem is set up with periodic boundary conditions according to and the condition of no flux at the bottom can be expressed as Dirichlet boundary conditions for a concentration θ(r,t) = 0.5 at the boundaries specified by Water 2020, 12, 2080 9 of 21 As an initial condition, it is assumed that θ = 0 in the volume V ab specified by the boundaries given by Equations (13)- (18) and z = 0.

Up-Flow SPB Storage: Numerical Set-Up
The capillary pressure, Equation (7), created between the rods within the annular clusters induces both radial and axial components of the flow. These components need to be calculated for each time step, depending on the positions of the flow front in both the inner and outer annular porous media, as illustrated in Figure 3. The pressures in the annular porous media are p i (x, r, t) and p o (x, r, t), respectively, where the flow and pressure are assumed to be independent of the azimuthal co-ordinate and the fluid is set to be incompressible. The perpendicular, i.e., in the radial direction, and parallel permeabilities, i.e., in the axial direction, in the annular porous media are denoted as K i,⊥ , K i, , K o,⊥ , and K o, , respectively. Darcy's law, Equation (5), may now be expressed as: where the modified pressures are ϕ i = p i + ρgx and ϕ o = p o + ρgx, respectively, and g is the gravitational acceleration. The permeabilities in Equation (19) are chosen to be given by where Π max = π/ 2 √ 3 in the case of hexagonal packing of the fibers [37] and Π b is the fraction of fibers. This is a simplification, but more advanced expressions, as e.g., derived in [38][39][40][41], can be used in future studies.
Equation (19) are set-up for the geometry and boundary conditions of interest according to Supplementary Materials, in which the derivations in [26] were followed, with two main differences: in the present model the finite widths of the channels are taken into account, and the gravity is considered. From this derivation the problem was solved in a semi-analytical fashion using MATLAB.

Verification of Numerical Solutions
By definition, the numerical derivations are approximations and need to be verified, as presented below.

Down-Flow SPB Storage: Verification
To verify the robustness of the numerical results for the Down-flow SPB storage an approximate analytical solution was found. This approximate solution is valid as long as the boundary layers of thickness ∆, that are initially developing on all the boundaries given by Equations (13)- (18), are much smaller than the characteristic dimensions of the Down-flow SPB storage, i.e., ∆ << b, B, H and h. Considering the different boundaries in Equations (13)-(18) separately, a similarity solution for each boundary is found. To exemplify, in the region from x = b/2 to x = B/2 the approximate solution is given by Similar solutions can be obtained for the other four boundaries. The total volume of water absorbed into the volume V(t) as a function of time is then found by integrating the water concentration over the total volume V ab , according to Equation (4), yielding an approximate result given by As an additional verification, a simple asymptotic solution valid for short times can also be found by expanding Equation (22) in the limit t → 0 according to In addition to these expressions, the numerical results can be compared to the maximum possible uptake of water.
In Figure 5 the absorbed water volume diffusing into the volume V ab is plotted as a function of time for the three cases in Table 1 and for the three solutions. The solid line is the numerical result, the dotted line is the approximate analytical solution, Equation (22), and the dashed line is an asymptotic solution valid for short times, Equation (23). For cases Down1 (blue lines) and Down2 (red lines), the numerical and approximate analytical solutions overall agree quite well, but for short times there is a slight disagreement since the numerical solution cannot capture the rapid t 1/2 initial development of the analytical solution. For case Down3 (black line), the agreement for short times is good, while for long times only the numerical solution fulfills the correct asymptotic volume of 75 L. The analytical solution here becomes invalid, because the boundary layers are no longer small compared to the dimensions of case Down3. It should be noted that Down1 and Down2 also approach the correct asymptotic volume, which takes a long time. The agreement with the analytical solutions and the final volume of water for all cases indicates that the numerical solution is correct.
Similar solutions can be obtained for the other four boundaries. The total volume of water absorbed into the volume V(t) as a function of time is then found by integrating the water concentration over the total volume Vab, according to Equation (4), yielding an approximate result given by As an additional verification, a simple asymptotic solution valid for short times can also be found by expanding Equation (22) in the limit → 0 according to In addition to these expressions, the numerical results can be compared to the maximum possible uptake of water.
In Figure 5 the absorbed water volume diffusing into the volume Vab is plotted as a function of time for the three cases in Table 1 and for the three solutions. The solid line is the numerical result, the dotted line is the approximate analytical solution, Equation (22), and the dashed line is an asymptotic solution valid for short times, Equation (23). For cases Down1 (blue lines) and Down2 (red lines), the numerical and approximate analytical solutions overall agree quite well, but for short times there is a slight disagreement since the numerical solution cannot capture the rapid t 1/2 initial development of the analytical solution. For case Down3 (black line), the agreement for short times is good, while for long times only the numerical solution fulfills the correct asymptotic volume of 75 L. The analytical solution here becomes invalid, because the boundary layers are no longer small compared to the dimensions of case Down3. It should be noted that Down1 and Down2 also approach the correct asymptotic volume, which takes a long time. The agreement with the analytical solutions and the final volume of water for all cases indicates that the numerical solution is correct.

Up-Flow SPB Storage: Verification
The results of the model can be verified and validated for the case of flow in a single channel with no porous medium, but including the effects of gravity. This is the famous experiment of capillary rise in a tube, with the well-known approximate analytical solution by Washburn [23]. More recently, Fries and Dreyer [42] found an exact analytical solution for this case and also tested it experimentally, with good agreement. For the case of flow in two annular regions, the maximum heights reached in the two channels (see Equation (25), in Supplementary Materials) can be verified by energy arguments. In the initial phase of the development of the fronts in each channel, when the interaction between the channels is zero, the expressions for the fronts Equation (24) in Supplementary Materials also agree with the results of Fries and Dreyer [42]. In Figure 6 a comparison is shown between the numerical results for the inner and outer channels given by Equations (22) and (23), in Supplementary Materials and the result obtained by using the analytical solution of Fries and Dreyer without any interaction. It is seen that, during this initial phase, in which there is a small interaction between the channels, the agreement is good. In the work by Zarandi et al. [27] the analytical solution by Fries and Dreyer [42] was also validated for porous media. In their case the porous medium consisted of several different kinds of glass-fiber wicks. The agreement between theory and experiments is good for several cases while for other cases there is a clear difference. The conclusion is that the disagreement in some cases can be derived from inhomogeneities, caused by 'kinking' fibers in the porous medium. Hence, the theory of Fries and Dreyer [42] provides an upper bound of the experimental data. In this context it is also important to mention that the surface energies (contact angle) may play an important role and that there are measures to reduce this angle as shown by Caglar et al. [29]. Including additional mechanisms, Caupin et al. [43] derived an upper limit for the capillary rise, which on the nano-scale was very large and exceeded Jurin's law. Experimental evidence is still required. The results of the model can be verified and validated for the case of flow in a single channel with no porous medium, but including the effects of gravity. This is the famous experiment of capillary rise in a tube, with the well-known approximate analytical solution by Washburn [23]. More recently, Fries and Dreyer [42] found an exact analytical solution for this case and also tested it experimentally, with good agreement. For the case of flow in two annular regions, the maximum heights reached in the two channels (see Equation (25), in Appendix A) can be verified by energy arguments. In the initial phase of the development of the fronts in each channel, when the interaction between the channels is zero, the expressions for the fronts Equation (24) in Appendix A also agree with the results of Fries and Dreyer [42]. In Figure 6 a comparison is shown between the numerical results for the inner and outer channels given by Equations (22) and (23), in Appendix A and the result obtained by using the analytical solution of Fries and Dreyer without any interaction. It is seen that, during this initial phase, in which there is a small interaction between the channels, the agreement is good. In the work by Zarandi et al. [27] the analytical solution by Fries and Dreyer [42] was also validated for porous media. In their case the porous medium consisted of several different kinds of glass-fiber wicks. The agreement between theory and experiments is good for several cases while for other cases there is a clear difference. The conclusion is that the disagreement in some cases can be derived from inhomogeneities, caused by 'kinking' fibers in the porous medium. Hence, the theory of Fries and Dreyer [42] provides an upper bound of the experimental data. In this context it is also important to mention that the surface energies (contact angle) may play an important role and that there are measures to reduce this angle as shown by Caglar et al. [29]. Including additional mechanisms, Caupin et al. [43] derived an upper limit for the capillary rise, which on the nano-scale was very large and exceeded Jurin's law. Experimental evidence is still required.  (22) and (23) in Appendix A against the analytical solution of Fries and Dreyer [42] within the limits of short times (t ≤ 4 s).

Parameter Sensitivity
Before comparing the modelled SPB capacities to rainfall data, a brief parameter sensitivity analysis is performed.

Down-Flow SPB Storage: Parameter Sensitivity
The value of B, which is the scale of water/SPB interface corrugations shown in Table 1 and Figure 2, influences the volumetric flow rate to a large extent, as evident from comparing the solid lines of different colors in Figure 5. Hence, the smaller the scale studied, the faster is the overall absorption process. For the smallest scale, Down3, where the sides of the open squares are 10 mm,

Parameter Sensitivity
Before comparing the modelled SPB capacities to rainfall data, a brief parameter sensitivity analysis is performed.

Down-Flow SPB Storage: Parameter Sensitivity
The value of B, which is the scale of water/SPB interface corrugations shown in Table 1 and Figure 2, influences the volumetric flow rate to a large extent, as evident from comparing the solid lines of different colors in Figure 5. Hence, the smaller the scale studied, the faster is the overall absorption process. For the smallest scale, Down3, where the sides of the open squares are 10 mm, the 1 m 2 flat surface can absorb around 75 L in 15 min. The sensitivity of the results due to diffusion, D, is given by the absorbed volume dependence on square root of the diffusion coefficient, according to the asymptotic expression (23). It only affects the growth of the volume in time, not the total absorbed volume: see the difference between the solid lines D = 1.0 × 10 −8 m 2 /s, the dashed lines, D = 0.5 × 10 −8 m 2 /s, and the dotted lines D = 1.0 × 10 −9 m 2 /s in Figure 7. In this context, the variation of the parameter b (=b/B, see Figure 2) is also interesting. Although this parameter should be set large for rapid growth, a large value of b yields a smaller total volume of absorption, as best seen by scrutinizing the solid lines of different colors in Figure 7. Here blue lines denote b = 0.4, red b = 0.5, and green b = 0.6. The effect of varying H is close to linear and increases both the initial rate of growth as well as the total absorption.  Figure 7. In this context, the variation of the parameter (=b/B, see Figure 2) is also interesting. Although this parameter should be set large for rapid growth, a large value of yields a smaller total volume of absorption, as best seen by scrutinizing the solid lines of different colors in Figure 7. Here blue lines denote = 0.4, red = 0.5, and green = 0.6. The effect of varying H is close to linear and increases both the initial rate of growth as well as the total absorption.

Up-Flow SPB Storage: Parameter Sensitivity
In all three cases in Table 2, one can anticipate initially fast inflows into both porous media, followed by a slowdown in the inflow (see Figure 8a,b). The graphs, however, reveal that, in the three cases studied, the water uptake differs. Initially the water uptake is fastest in case Up1, during the period 0 < t < 2 s, as shown in Figure 8a. Then, during a period up to slightly more than 100 s, the water uptake is greatest in case Up2, (Figure 8b). After that, the water uptake in the system with almost equal radii (Up3) dominates. Hence, the storage set-up with two porosities can be designed for both a fast water uptake, as well as the maximum volume of water uptake. The maximum water uptake is about 8, 20 and 46 L, for durations of 10, 100 and 1000 s, respectively.
The results may also be interpreted in terms of the water column height in the respective annuli. The maximum height is about 1.4 m for the cases studied, as exemplified by cases Up1-3 in Figure 9. The figure also shows that there is a large difference in the final height between the inner and outer annuli for Up1-2, which is unsatisfactory. Hence, having several annuli with different radii may make the dynamic storage structure more efficient, as a complement to using the porous media with less spread sizes of the fibers (Case Up3).

Up-Flow SPB Storage: Parameter Sensitivity
In all three cases in Table 2, one can anticipate initially fast inflows into both porous media, followed by a slowdown in the inflow (see Figure 8a,b). The graphs, however, reveal that, in the three cases studied, the water uptake differs. Initially the water uptake is fastest in case Up1, during the period 0 < t < 2 s, as shown in Figure 8a. Then, during a period up to slightly more than 100 s, the water uptake is greatest in case Up2, (Figure 8b). After that, the water uptake in the system with almost equal radii (Up3) dominates. Hence, the storage set-up with two porosities can be designed for both a fast water uptake, as well as the maximum volume of water uptake. The maximum water uptake is about 8, 20 and 46 L, for durations of 10, 100 and 1000 s, respectively.  Table 2). Graphs in the two panels correspond to two time intervals: (a) 0-1 min and (b) 0-1 h, respectively. Graph colors: the blue, red and green curves denote cases Up1, Up2 and Up3, respectively.  Table 2. The graph panels represent two time intervals, 0-1 min and 0-1 h, respectively. The blue, red and green curves represent cases Up1, Up2 and Up3, respectively. Dashed and dotted lines represent the height of the water column after the merging of fronts.
Concerning the other parameters in Table 2, it is obvious that the amount of water uptake will increase with an overall increase in δ, a and b (see also Figure 3). The effect of changing porosity (set in Table 2 as 0.7) is, however, less obvious and, therefore, such an effect is illustrated in Figure 10 where the color blue denotes ϕ = 0.6, red ϕ = 0.7 and green ϕ = 0.8 and the different line types represent the cases Up1-3. Hence, using sparser porous media will both speed-up the uptake and increase the maximum uptake of water and makes it evident that the porosity should be as high as possible, while keeping the capillary pressure sufficiently high.  Table 2). Graphs in the two panels correspond to two time intervals: (a) 0-1 min and (b) 0-1 h, respectively. Graph colors: the blue, red and green curves denote cases Up1, Up2 and Up3, respectively.
The results may also be interpreted in terms of the water column height in the respective annuli. The maximum height is about 1.4 m for the cases studied, as exemplified by cases Up1-3 in Figure 9. The figure also shows that there is a large difference in the final height between the inner and outer annuli for Up1-2, which is unsatisfactory. Hence, having several annuli with different radii may make the dynamic storage structure more efficient, as a complement to using the porous media with less spread sizes of the fibers (Case Up3).  (Table 2). Graphs in the two panels correspond to two time intervals: (a) 0-1 min and (b) 0-1 h, respectively. Graph colors: the blue, red and green curves denote cases Up1, Up2 and Up3, respectively.  Table 2. The graph panels represent two time intervals, 0-1 min and 0-1 h, respectively. The blue, red and green curves represent cases Up1, Up2 and Up3, respectively. Dashed and dotted lines represent the height of the water column after the merging of fronts.
Concerning the other parameters in Table 2, it is obvious that the amount of water uptake will increase with an overall increase in δ, a and b (see also Figure 3). The effect of changing porosity (set in Table 2 as 0.7) is, however, less obvious and, therefore, such an effect is illustrated in Figure 10 where the color blue denotes ϕ = 0.6, red ϕ = 0.7 and green ϕ = 0.8 and the different line types represent the cases Up1-3. Hence, using sparser porous media will both speed-up the uptake and increase the maximum uptake of water and makes it evident that the porosity should be as high as possible, while keeping the capillary pressure sufficiently high.  Table 2. The graph panels represent two time intervals, 0-1 min and 0-1 h, respectively. The blue, red and green curves represent cases Up1, Up2 and Up3, respectively. Dashed and dotted lines represent the height of the water column after the merging of fronts.
Concerning the other parameters in Table 2, it is obvious that the amount of water uptake will increase with an overall increase in δ, a and b (see also Figure 3). The effect of changing porosity (set in Table 2 as 0.7) is, however, less obvious and, therefore, such an effect is illustrated in Figure 10 where the color blue denotes φ = 0.6, red φ = 0.7 and green φ = 0.8 and the different line types represent the cases Up1-3. Hence, using sparser porous media will both speed-up the uptake and increase the maximum uptake of water and makes it evident that the porosity should be as high as possible, while keeping the capillary pressure sufficiently high.

Modelling the SPB Storage Interception of Short-Duration Design Rainfalls
From the hydrological point of view, two essential properties of SPB storage are the rates of inflow and the total storage volume, which should be selected according to the design rainfall characteristics. Ultimately, the SPB structures fully intercepting and storing the incoming rainwater would provide the highest attainable runoff control, but such a condition should not be viewed as a prerequisite in assessing the feasibility of SPB storage applications in stormwater control. One can envisage integrated storage designs, in which SPB storage would complement specific features of conventional storage (e.g., when applied as one of the layers of a green roof). Also note that LSMs, and even larger scale measures, are sometimes designed to partly bypass high flows, as noted, e.g., for green roofs by Shafique et al. [11]. With this in mind, capacities of SPB storage inflow rates and volumes are theoretically assessed in this section for specific Swedish short-duration extreme rainfall data adopted from Olsson et al. [44].

Inflow
Two types of inflow into stormwater storage can be distinguished: direct rainfall over the footprint of the storage facility, and indirect inflow diverted from adjacent drainage contributing areas. The indirect inflow generally exceeds the direct by a significant factor given by the ratio of contributing areas (Aindirect/Adirect). Direct inflow in wet weather can be determined from the local rainfall regime [44], and for the purpose of this study the following rainfall data were adopted from Swedish precipitation records (see Table 3): (a) 60-min block rainfalls of return period of 10-years, for southwestern and northern regions of Sweden (i.e., the regions with the highest and lowest annual precipitation depths, respectively), and (b) 60-min block rainfalls, with a 5-min high-intensity rainfall burst (starting at 27.5 min since rain onset), with the return period of 10 years, for the southwestern and northern regions.

Modelling the SPB Storage Interception of Short-Duration Design Rainfalls
From the hydrological point of view, two essential properties of SPB storage are the rates of inflow and the total storage volume, which should be selected according to the design rainfall characteristics. Ultimately, the SPB structures fully intercepting and storing the incoming rainwater would provide the highest attainable runoff control, but such a condition should not be viewed as a prerequisite in assessing the feasibility of SPB storage applications in stormwater control. One can envisage integrated storage designs, in which SPB storage would complement specific features of conventional storage (e.g., when applied as one of the layers of a green roof). Also note that LSMs, and even larger scale measures, are sometimes designed to partly bypass high flows, as noted, e.g., for green roofs by Shafique et al. [11]. With this in mind, capacities of SPB storage inflow rates and volumes are theoretically assessed in this section for specific Swedish short-duration extreme rainfall data adopted from Olsson et al. [44].

Inflow
Two types of inflow into stormwater storage can be distinguished: direct rainfall over the footprint of the storage facility, and indirect inflow diverted from adjacent drainage contributing areas. The indirect inflow generally exceeds the direct by a significant factor given by the ratio of contributing areas (A indirect /A direct ). Direct inflow in wet weather can be determined from the local rainfall regime [44], and for the purpose of this study the following rainfall data were adopted from Swedish precipitation records (see Table 3): (a) 60-min block rainfalls of return period of 10-years, for southwestern and northern regions of Sweden (i.e., the regions with the highest and lowest annual precipitation depths, respectively), and (b) 60-min block rainfalls, with a 5-min high-intensity rainfall burst (starting at 27.5 min since rain onset), with the return period of 10 years, for the southwestern and northern regions.

Comparison of Rainfall Volumes to Storage Capacities: Down-Flow SPB Storage
The Down-flow SPB storage provides a direct capture of rainwater, and the comparison between the numerical solution (black solid line) in Figure 5 (i.e., Down3) and the rainwater inflows from Table 3 reveals that the rainwater volumes for the two events in the two regions are fully captured in the case of Down3 (see Figure 11). It should be noted that the Down-flow SPB can capture more rainwater than that supplied; thus, in the Down3 case, there is some safety available should the storage capacity deteriorate for unknown reasons. This point will be further discussed towards the end of this section. It can be inferred from data in Figure 11 that the SPB storage concept proposed here has a great potential to capture and store all the rainwater associated with the short-duration extreme rainfall events studied.

Comparison of Rainfall Volumes to Storage Capacities: Down-Flow SPB Storage
The Down-flow SPB storage provides a direct capture of rainwater, and the comparison between the numerical solution (black solid line) in Figure 5 (i.e., Down3) and the rainwater inflows from Table 3 reveals that the rainwater volumes for the two events in the two regions are fully captured in the case of Down3 (see Figure 11). It should be noted that the Down-flow SPB can capture more rainwater than that supplied; thus, in the Down3 case, there is some safety available should the storage capacity deteriorate for unknown reasons. This point will be further discussed towards the end of this section. It can be inferred from data in Figure 11 that the SPB storage concept proposed here has a great potential to capture and store all the rainwater associated with the short-duration extreme rainfall events studied.

Comparison of Cumulative Rainfall Volumes and Theoretical Storage Capacities: Up-Flow SPB Storage
The Up-flow SPB storage, that has a minor footprint, is not intended to directly capture rainwater, but to store stormwater from a certain contributing area, draining into individual storage facilities. Hence, various contributing areas A should be considered for individual storage arrangements, formed by sets of SPB cylinders. Assuming that stormwater can easily enter the confined storage area without any losses along the flow path, the Up-flow SPB (green curve in Figure  8, case Up3) is compared to eight design rainfall scenarios in Figure 12. The scenarios are listed and were defined as follows: two regional rainfalls (SW & N), two events (60 min rainfalls with or without a high intensity burst), and two contributing areas, 1 and 5 m 2 . The comparison indicates that one Up-flow SPB storage unit per unit area, Up3/m 2 , has a greater storage capacity than needed to capture all runoff, but one storage unit/(5 m 2 ) is not enough to capture all runoff, regardless of the design event or the geographical region. Nevertheless, the Up-flow SPB storage can collect a great amount of water, and with a high density of storage units, it would be feasible to capture all runoff. It should also be noted that, in terms of capacity, Up3 is merely the "best" case among the three cases studied, rather than an optimized solution.

Comparison of Cumulative Rainfall Volumes and Theoretical Storage Capacities: Up-Flow SPB Storage
The Up-flow SPB storage, that has a minor footprint, is not intended to directly capture rainwater, but to store stormwater from a certain contributing area, draining into individual storage facilities. Hence, various contributing areas A should be considered for individual storage arrangements, formed by sets of SPB cylinders. Assuming that stormwater can easily enter the confined storage area without any losses along the flow path, the Up-flow SPB (green curve in Figure 8, case Up3) is compared to eight design rainfall scenarios in Figure 12. The scenarios are listed and were defined as follows: two regional rainfalls (SW & N), two events (60 min rainfalls with or without a high intensity burst), and two contributing areas, 1 and 5 m 2 . The comparison indicates that one Up-flow SPB storage unit per unit area, Up3/m 2 , has a greater storage capacity than needed to capture all runoff, but one storage unit/(5 m 2 ) is not enough to capture all runoff, regardless of the design event or the geographical region. Nevertheless, the Up-flow SPB storage can collect a great amount of water, and with a high density of storage units, it would be feasible to capture all runoff. It should also be noted that, in terms of capacity, Up3 is merely the "best" case among the three cases studied, rather than an optimized solution.
In both comparisons, the SPB storage units modelled (i.e., Down-flow and Up-flow) could capture and store more water than available early during the rainfall event, which is physically impossible. For a more realistic comparison between the cumulative rainfall volume and storage filling, the numerical model for the Up-flow case was modified to handle events with a constant rate of inflow. Using this new model and focusing on Case Up3 with the largest uptake of water during the first 60 min, it appears that at the beginning of the rainfall event, the inflows into the Up-flow SPB storage follow the SW region runoff inflows for various contributing areas: 1, 2.5, 5 and 10 m 2 (see Figure 13). After some time, depending on the contributing area, the storage can no longer fully capture runoff from respective areas. For the area of 2.5 m 2 , this tipping point occurs just before the 60-min rainfall event ends, and hence, the SPB unit can almost capture runoff from a contributing area of 2.5 m 2, exposed to the SW region design event with a uniform intensity (see the black dashed and dashed lines in Figure 13). Therefore, this finding provides evidence that the presented Up-flow SPB concept has great potential for controlling runoff. A similar comparison for the Down-flow SPBs would yield the same conclusion. Notice that the comparison will hold also for multiple closely time-spaced rain showers. In both comparisons, the SPB storage units modelled (i.e., Down-flow and Up-flow) could capture and store more water than available early during the rainfall event, which is physically impossible. For a more realistic comparison between the cumulative rainfall volume and storage filling, the numerical model for the Up-flow case was modified to handle events with a constant rate of inflow. Using this new model and focusing on Case Up3 with the largest uptake of water during the first 60 min, it appears that at the beginning of the rainfall event, the inflows into the Up-flow SPB storage follow the SW region runoff inflows for various contributing areas: 1, 2.5, 5 and 10 m 2 (see Figure 13). After some time, depending on the contributing area, the storage can no longer fully capture runoff from respective areas. For the area of 2.5 m 2 , this tipping point occurs just before the 60-min rainfall event ends, and hence, the SPB unit can almost capture runoff from a contributing area of 2.5 m 2, exposed to the SW region design event with a uniform intensity (see the black dashed and dashed lines in Figure 13). Therefore, this finding provides evidence that the presented Up-flow SPB concept has great potential for controlling runoff. A similar comparison for the Down-flow SPBs would yield the same conclusion. Notice that the comparison will hold also for multiple closely timespaced rain showers.  Comparison between the cumulative rainfall volume and Up-flow SPB (Up3) storage filling (black solid line) for various inflow scenarios: two climatic regions (N and SW), two shapes of rainfall hyetographs (block rainfalls without and with a high intensity burst; blue and red lines, respectively), and two runoff contributing areas (1 and 5 m 2 , respectively). The figure also demonstrates the storage capacity of the SPB.
In both comparisons, the SPB storage units modelled (i.e., Down-flow and Up-flow) could capture and store more water than available early during the rainfall event, which is physically impossible. For a more realistic comparison between the cumulative rainfall volume and storage filling, the numerical model for the Up-flow case was modified to handle events with a constant rate of inflow. Using this new model and focusing on Case Up3 with the largest uptake of water during the first 60 min, it appears that at the beginning of the rainfall event, the inflows into the Up-flow SPB storage follow the SW region runoff inflows for various contributing areas: 1, 2.5, 5 and 10 m 2 (see Figure 13). After some time, depending on the contributing area, the storage can no longer fully capture runoff from respective areas. For the area of 2.5 m 2 , this tipping point occurs just before the 60-min rainfall event ends, and hence, the SPB unit can almost capture runoff from a contributing area of 2.5 m 2, exposed to the SW region design event with a uniform intensity (see the black dashed and dashed lines in Figure 13). Therefore, this finding provides evidence that the presented Up-flow SPB concept has great potential for controlling runoff. A similar comparison for the Down-flow SPBs would yield the same conclusion. Notice that the comparison will hold also for multiple closely timespaced rain showers.

Discussion
The discussion focuses on the following aspects of SPB storage: potential role of SPB storage in stormwater management, water uptake (speed and volume) and its modelling, filtration of stormwater by SPBs, and other practical considerations.
Analysis of the SPB storage concept presented here places this type of storage into the category of 'best practice', serving to manage stormwater at the smallest scale, i.e., the lot or site level (defined earlier as LSMs). While the quantified effect of a single measure on the catchment runoff is not measurable, the overall benefit of many such measures throughout the catchment is significant as argued, e.g., in the case of green roofs [45]. Furthermore, these localized distributed measures are not required to intercept the incoming stormwater fully, because of their multiplicity and hierarchical placement in the upper reaches of the catchment. For example, another similar LSM, green roofs, on average retain 57-76% of incoming rainwater [45] and bypass the rest, to be intercepted by other downstream measures or be conveyed to the receiving waters. Functionally, Down-flow SPB storage compares favorably to gray roofs (i.e., roof storage without vegetation), but the literature on gray roofs is rather limited. In terms of storage filling, the Down-flow SPB storage was shown here to be capable of fully intercepting design rainfall from 10-year cool temperate climate events. Furthermore, one can also envisage the SPBs as part of an integrated system containing various control measures, in which SPBs complement other measures. Compared to green roofs, gray roofs or SPB storage measures would reduce the load of the storage structure on the roof, eliminate the weight of green roof vegetation substrate (typically 200 mm [45]) and avoid water quality issues connected with elution of green roof media. However, these benefits have to be measured against the ecological benefits of green roofs, even though Francis and Jensen's review [17] points out that reports of green roof benefits in the form of ecosystem services are sometimes insufficiently documented. Up-flow SPBs have a limited footprint and could play a different role in local runoff management than down-flow SPBs, by forming vertical storage units, or, in the case of unsupported rods growing from the ground up, by temporarily retaining stormwater over infiltration soakaways, or ground depressions, and thereby increasing their effectiveness in reducing runoff. Table 3. Direct inflow of rainwater into storage.

Direct Rainwater Inflow Velocity (=Intensity) [m/s] per Unit Area [m 2 ] Southwest (SW) North (N)
Uniform intensity (also called the block rainfall) 60-min duration event, with return period 1:10 years 0.68 × 10 −5 no preceding rainfall 0.53 × 10 −5 no preceding rainfall 60-min duration event, with a high-intensity burst of 5 min, return period 1:10 years (the hyetograph was patterned after Berggren [46] The computational models presented in this study focused on the inflow to the storage facility, but the outflow of water during storage drawdown, at rates much lower than those of inflow, is also important and subject to ongoing investigation. There is a significant difference between the time scales of the inflow and outflow rates associated with design events: the former can be as short as 1 h [46] and the latter can be as long as 48 h, depending on the local rainfall regime [9]. Hence, the average outflow rate would be about 50 times smaller than the inflow rate of design events. Theoretically, water may be released, e.g., by changes in capillary action, by mechanisms applying pressure to stored water, and by evaporation of water into the atmosphere [47][48][49][50][51]. All three measures have different time-scales, and their design would need to meet the time limit on the full drawdown of the storage facility, typically specified as 48 h (such limits depend on the type of storage facility and the local rainfall regime).
For Down-flow SPBs, only fine particles will be present in the rainwater entering the SPB storage unit, and filtration may be a minor issue, but airborne materials, dust, leaves, twigs, and bird droppings may interfere with the transport of water into the absorbing material. One practical option would be to apply some pre-treatment of rainwater, in the form of diversion or screening [52]. Even with this simple pre-treatment, it is likely that impurities collected on the SPB surface would form a biological filtration layer (usually called "a cake" or "schmutzdecke"), as described, e.g., by Li and Davis [53] and Tien [54], but the SPB material itself should not be affected by the "particles". However, the functioning of the storage will deteriorate when "particles" start to fill the cavities, serving to increase the interface area and cover the surface of the absorbing material as the cake is formed. This implies that the flow rates into and out of the storage would be impaired, while the amount of water that can be stored would remain about the same. This type of filtration can be modelled on various spatial scales and represents a possible topic for future research, e.g., building on the work of Tien [54], Li and Davis [55], Frishfelds et al. [56], Lundström et al. [57], and Zhou et al. [58]. From the operational point of view, "the cake" needs to be removed at certain time intervals, which is commonly achieved in practice by periodic flushing.
The Up-flow SPBs are also exposed to the entry of water carrying various materials (e.g., solids), chemicals (e.g., nutrients, dissolved organic matter), and fecal microorganisms attributed mostly to pets and wildlife droppings [59]. However, these SPBs should be less susceptible to clogging by larger particles, because of gravity forces acting on such particles in the direction opposite to that of flow (away from the filter membrane). In general, the storage structure needs to be designed to sustain its operation and functioning taking into account such influent quality. The first step for enhancing the storage structure operation would be, again, to include some kind of pre-treatment of the incoming stormwater (mostly settling, or screening). Even with pre-treatment, the SPB storage system needs to be able to handle solids of various sizes. Conceivably, fine particles could be allowed to enter into the vertical structures, in a process mimicking deep bed filtration, while larger particles would be filtered away, in the form of cake filtration, or deposition at the SPB structure foot. The advantage of this approach would be that the cake will have a relatively high permeability since it will consist of larger size particles/objects. When allowing smaller particles to enter the storage, the filtration process should not impair the flow rates to the same extent as in the case of cake filtration of such particles [57]. In addition, in this case, in-depth studies may be carried out for the dynamic storage itself, including porous media with various sizes of pore structures [55,60,61]. However, concerns about the porous material surface fouling will need to be addressed. While the treatment of stormwater by filtration in SPBs is a challenge for flow entry and storage filling, it offers a great opportunity for enhancing the stored stormwater quality. This dual functionality, providing both quantity and quality control benefits, would improve cost-benefit considerations of SPB storage facilities and their economic attractiveness.
In addition to the storage outflow, several practical issues of operating SPB storage may need to be resolved before starting prototype field testing: rainwater or stormwater filtration upon storage entry through the water/SPB interface, formation of biofilms impeding inflow (fouling), environmental concerns, resistance to damage or deterioration of SPB due to freezing, and safety of designs incorporating SPB storage. For filtration, the conditions at the point of inflow into SPBs are of high importance, i.e., the SPB performance in receiving, storing and treating rainwater or stormwater carrying suspended solids and chemicals. Consequently, when looking for practical applications of SPB storage, one should start with relatively "clean" waters (e.g., rainwater on rooftops).
The simulated test cases show that significant volumes of water can be stored relatively quickly in the SPB storage units analyzed and that the geometry of such units could be further optimized. The modelling experiments with relevant inflow rates resulting from short-duration design rainfalls furthermore indicate that the concepts presented could be used to contribute to the restoration of catchment hydrology and reduce runoff even during heavy rainfalls. Although the design rainfall events are generally developed for return periods ranging from 2 to 100 years [44], it is conceivable that the SPBs would be best applied within a lower range of return period, perhaps 2-25 years, and in an integrated manner with other types of storage, existing or planned, in the catchment. For example, the overflows from SPBs could be diverted to conventional storage, and vice versa.

Conclusions
Potential roles of dynamic water storage in sponge-like porous bodies (SPBs) in lot-level stormwater management were proposed and theoretically examined. To this end, mathematical analysis and numerical modelling, based on first principles, demonstrated that SPBs could be designed to fully capture Swedish design rainfalls of 1-h duration and the average return period of 10 years. Furthermore, such analyses demonstrated that the use of absorbing and/or porous media on several scales in self-driven storage designs can be optimized with respect to the theoretical maximum storage capacity and inflow rates. Hence, the potential of the concepts discussed is much greater than what could be presented in this paper. At the same time, it is evident that additional theoretical and experimental work is needed to advance the Technology Readiness Level of the proposed theoretical SPB concept beyond the current level.