Education and communication 14 May 2020
Education and communication  14 May 2020
Wetropolis extreme rainfall and flood demonstrator: from mathematical design to outreach
 ^{1}School of Mathematics and Leeds Institute for Fluid Dynamics, University of Leeds, LS2 9JT, Leeds, UK
 ^{2}Wowlab, Enschede, the Netherlands
 ^{1}School of Mathematics and Leeds Institute for Fluid Dynamics, University of Leeds, LS2 9JT, Leeds, UK
 ^{2}Wowlab, Enschede, the Netherlands
Correspondence: Onno Bokhove (o.bokhove@leeds.ac.uk)
Hide author detailsCorrespondence: Onno Bokhove (o.bokhove@leeds.ac.uk)
Government and consulting experts on flood mitigation generally face difficulties when trying to explain the science of extreme flooding to the general public, in particular the concept of a return period. Too often, for example, people perceive they are safe for the next 100 years after a 1:100year returnperiod flood has hit their town. UK flood practitioners therefore gave us the challenge to design an outreach tool that conceptualises the science of flooding in a way that is accessible to and directly engages the public, and in particular demonstrates what a return period is. Furthermore, we were tasked with designing a live 3D physical model rather than a graphical or animated 2D game on a screen. We show here how we tackled that challenge by designing, constructing, and showcasing the Wetropolis Flood Demonstrator. Wetropolis is a transportable and conceptual physical model with random rainfall, river flow, a flood plain, an upland reservoir, a porous moor, representing the upper catchment and visualising groundwater flow, and a city which can flood following extreme and random rainfall. A key novelty is the supply of rainfall every Wetropolis day. Several aspects of Wetropolis are considered.
 i.
We present the modular mathematical and numerical design on which Wetropolis is based. It guided the choice of parameter values of Wetropolis, which was loosely inspired by the Leeds Boxing Day floods of the River Aire in 2015. The design model further serves as the building block and inspiration for adaptations suited to particular local demands. Moreover, the model is purposely lean and therefore quick to compute, serving flexibility in the outreachtool design, but is less suitable for any detailed scientific validation.
 ii.
The constructed Wetropolis is described here in broad terms, but we include a link to a GitHub site with details to inspire other bespoke designs. The goal, again, is to facilitate new adaptations of Wetropolis for particular catchments different to the Leeds River Aire case.
 iii.
Our experience in showcasing Wetropolis is summarised and discussed, with the purpose of giving an overview as well as inspiring improved and bespoke adaptations. While Wetropolis should be experienced live, with videos found on the GitHub site, here we provide a photographic overview. To date, Wetropolis has been showcased to 500 to 1000 people at public workshops and exhibitions on recent UK floods, as well as to flood practitioners and scientists at various research and stakeholder workshops.
 iv.
We conclude with some ongoing design changes, including how people can experience natural flood management in a revised Wetropolis design. Finally, we also discuss how Wetropolis, although originally focussed solely on outreach, led to a new costeffectiveness analysis and protocol for assessing floodmitigation plans and inspired other physical models for use in education and water management.
The Boxing Day flood of 2015, and more recently Storm Ciara in 2020, caused widespread damage in Yorkshire, UK, due to extreme flooding of the River Aire in and around Leeds and the River Calder in and around Todmorden, Mytholmroyd, and Hebden Bridge. Thankfully no fatalities occurred, but the economic damage was severe and estimated to be around GBP 500 million (West Yorkshire combined Authority, 2016). November 2015 was the third wettest month on record in terms of precipitation and December 2015 was the wettest on record (Environment Agency, 2016). As such, the soil was already saturated when 48 h of extreme rainfall fell in Yorkshire, leading to the Boxing Day floods of 2015. After this and other recent floods, many flood victims asked why extreme flood events, seemingly occurring more and more often, were causing such havoc in their communities.
To provide some mathematical background on modelling, mitigation, and statistics of extreme flood events, we were asked to disseminate the scientific background on “risk in the age of extremes” at the “Science of floods” citizens' conference in Hebden Bridge (Science of Floods, 2016). In the first 2 decades of the 2000s, including the Boxing Day and Storm Ciara floods, Hebden Bridge was hit by both summer flash floods and winter floods, leading to concerns amongst flood victims that their lives and properties were insufficiently protected. It led to further and intense discussions with environmental agencies on the need for more and different types of flood defences. Important questions in this discussion are the following.

Is it going to rain more in the future in the UK?

Can we define extreme precipitation and flooding events?

How (well) can we predict heavy precipitation and floods?

Finally, how can we elucidate these questions, their answers, and uncertainties, in an interactive tabletop demonstration? The latter question was posed to us by experts on flood mitigation, in their wish to communicate better with the general public. In particular, experts faced troubles explaining what a return period is in flooding, alongside general difficulties explaining statistical notions involved in flood mitigation.
We will discuss some answers to these questions in turn given their relevance to Wetropolis' design.
1.1 Is it going to rain more in the future in the UK?
Both the IPCC report (Stocker et al., 2013) and Sanderson's UKCP09 report of the Met Office (Sanderson, 2010) show that there is no increase in significance of average annual rainfall foreseen in climate projections, not across the globe on average and also not in the UK. However, there are geographical and seasonal variations foreseen: winter rainfall will generally increase and summer rainfall will thus decrease, but with more intense downpours.
1.2 Can we define extreme precipitation and flooding events?
Extreme events tend to be expressed as the chance that an extreme event occurs on one day in a year. If that chance is 1 %, for example, then we say that this extreme event has a return period of 1:100 years. Flooding events can be classified in terms of such return periods as e.g. 1:10, 1:20, 1:50 and 1:100 or 1:200year events, with the latter two considered to be extreme. The uncertainty in an event with a 1:100year return period will be larger than one with a 1:20year return period. An extreme event occurs in the tail of a probability distribution and may have never been observed. Events with return periods longer than the data record cannot be classified directly from that data and must be determined using theoretical probability distributions. However, the low number of extreme events in a finitetime data set means it is difficult to establish accurately the tail of the distribution. Accordingly, there is a great deal of uncertainty associated with extreme events, which should ideally be quantified and communicated effectively. By assuming a suitable (parametric) probability distribution function (pdf), one can use the data to fit the parameters of the pdf and subsequently generate data (“sample the pdf”) in the tail for events with return periods beyond the length of the data record. Classical pdfs, such as (half a) Gaussian distribution or a Gamma distribution, are typically used to model rainfall intensities, but are not suitable for extreme events. Extremevalue theory offers a family of distributions to overcome the limitations of the classical pdfs. Given a sufficient amount of (rainfall or riverlevel) data over a certain threshold, the generalised Pareto distribution (GPD) attempts to model the probabilities of extreme events beyond the data record. It is also possible to capture both the bulk and the extreme tail of the distribution in a socalled mixture model by combining a Gamma and generalised Pareto distribution (GPD) (cf. Wong et al., 2014). The above is relevant because within Wetropolis we will use a “discrete” distribution with a “rare” tail.
1.3 How (well) can we predict heavy precipitation and floods?
The Boxing Day floods of 2015 were caused by largescale winter rainfall. The 48 h of consecutive rainfall in the days leading up to the flood were the wettest on record: in Bradford 69.4 mm and in Bingley 93.6 mm of rainfall were measured over 48 h. We will use this direct response of floods driven by 1 or 2 d of heavy rainfall in our design. It resulted in the flooding of the River Aire, with riverlevel records reached in Leeds and elsewhere along the river. In Armley, Leeds, the gauge station measured a maximum river level of 5.22 m, while the previous electronic record was 4.03 m in the autumn of 2000 (Environment Agency, 2016). The river level during the 1866 flood was roughly around 4.5 m; cf. Bokhove et al. (2020) (their Fig. 3). Both rainfall and river levels were by and large wellpredicted by the UK Met Office via numerical weather prediction and by the Environment Agency (– EA; cf. a presentation by an EA – Yorkshire leader in Leeds). Predictions are generally quite good for largescale winter rainfall and the resultant changes in river levels. Downpours, e.g. in the summer, tend to be more localised and are therefore much more difficult to predict in terms of location, intensity, and duration. The same holds for resulting flash floods and downpourinduced surfacewater flooding. Hence, simply put, fluvial or river flooding in winter tends to be easier to predict than pluvial or surfacewater flooding events in the summer.
1.4 Finally, how can we elucidate these questions, their answers, and uncertainties in an interactive tabletop demonstration?
The decision to design a tabletop demonstrator was triggered by the desire of flooding experts to have a 3D physical setup as opposed to 2D animations or computer graphics on a screen. It follows a recent trend to use transportable physical setups, such as the coastal wave tank of JBA Trust (an online video has attracted to date over 5.7 million views)^{1}. The expression of extreme events in terms of return periods is difficult to grasp and often misunderstood, especially by the public. The Boxing Day flood of 2015 was classified as an event with a $\mathrm{1}:{\mathrm{200}}^{+}$year return period – including the unclear meaning of the plus sign in 200^{+} (Environment Agency, 2016)^{2}. That does not mean that it has to take another 200^{+} years, so until after 2215, before the next Boxing Daytype flood might occur in Leeds. It does, however, mean that the average time between events of similar magnitude will be 200^{+} years, given a sufficiently long record of “stationary” statistical data. To let people experience such an extreme event in a tabletop setup, they can of course not be asked to wait for 200 years on average, so our design for a flood demonstrator with rainfall must be scaled down in both size and duration. Miniature river flooding has been demonstrated in smallscale experiments (e.g. as in the Lego model of Pampaloni et al., 2018), but these all tend to involve deterministically imposed extreme water input – with water inflow supplied and adjusted deterministically and/or manually. The key novelty in our design lies in the way rainfall is supplied randomly to our tabletop hydrodynamic setup for both river and groundwater flow. We have modified the classical symmetric Galton board, inspired by such a setup used at Leeds' School of Mathematics open days. A typical Galton board has a tilted surface in which a (steel) ball falls down under gravity and encounters a series of symmetric pins or channel corners, each determining with a p and (1−p) chance whether the ball continues or falls to the left or to the right. The design is usually such that $p\approx \mathrm{1}/\mathrm{2}$ but small variations can occur in practice, and after a series of n rows of splittings a binomial distribution arises, given a sufficiently large number of trials. Moreover, for n→∞, a Gaussian distribution emerges. The Galton board is often used to visualise and demonstrate statistical distributions in real time during e.g. outreach events. To obtain an asymmetric discrete distribution with a discrete tail representing relatively extreme events, the standard symmetric Galton board described above was modified as follows. For $p=\mathrm{1}/\mathrm{2}$ and n=1, the first and only split leads to a (1, 1)∕2 distribution. The first split of the second row for n=2 is now eliminated, while the second split is not, leading to a (3, 1)∕4 distribution. Continuing to the third row of splittings as usual, for n=3, we obtain a (3, 4, 1)∕8 distribution. The last and fourth row for n=4 yields the final (3, 7, 5, 1)∕16 distribution. An image of such a asymmetric Galton board is given in Fig. 1. Two of these Galton boards will be used to supply rain to our tabletop river and groundwater flow model, one concerning the duration and amount (0.1, 0.2, 0.4, 0.9)r_{0} of rainfall during a Wetropolis day, with its unit wd, and another one concerning the location of the rainfall. Rain duration will be either (10 %, 20 %, 40 %, 90 %) of the amount r_{0} of rainfall per Wetropolis day and rain location will be either rainfall (i) in a reservoir with generally instant runoff into the river; (ii) in both a reservoir and on a moor; (iii) on the moor with groundwater flow and its nonlinear, delayed release of water into the river; and (iv) no rain, in the Wetropolis catchment. See Fig. 2 for a plan view of Wetropolis. Both duration and location are determined by the outcome of one trial through two Galton boards per Wetropolis day, together yielding a 4by4 matrix of joint probabilities, with the norain case having a rare chance (rare for the “UK”) of 1∕16 comprised of 4 of those 16 outcomes. By design, an extreme event occurs when it rains 90 % in both locations (i.e. in the reservoir and moor) with a chance of $\mathrm{7}/\mathrm{256}\approx \mathrm{0.0273}=\mathrm{2.73}$ %, which in our construction will by design lead to flooding of a city further downstream along a (winding) river in the setup. A Boxing Daytype event with 2 consecutive days of extreme 90 % rainfall then has a chance of $\mathrm{49}/\left({\mathrm{256}}^{\mathrm{2}}\right)\approx \mathrm{0.000748}=\mathrm{0.075}$ %. The next and crucial step in the design is to identify and determine the various unknowns in order to ascertain whether a feasible design for a working physical setup is possible at all.
Given a (winding) river of length L and curvilinear coordinate s along this river, these remaining key unknowns are as follows:

the influx discharge Q_{0} at the upstream boundary at s=0;

the locations s_{res} and s_{m} where the reservoir and moor enter into the river (with the distance along the river given by a curvilinear coordinate s), with a section further downstream along the river comprising a city plain that is prone to (extreme) flooding;

the rainfall amount r_{0} (dimensionally a speed, as we will explain in the next section), determining the strength of the pumps required and whether pumping rates can be realistic at all; and

the length of a Wetropolis day, in relation to the extreme rainfall and corresponding extreme flooding event, such that the viewer experiences some irritation in having to wait for a randomly induced extreme event but on average will experience such an extreme event within a reasonable time, i.e. on average within several minutes.
We have chosen s_{res}, s_{m}, and Q_{0} a priori and will determine wd and r_{0} via visual optimisation of a series of simulations of a simplified mathematical and numerical model (as explained in Sect. 2.3). Note that the latter model is a lean design model exclusively geared towards obtaining quick estimates of the design parameters. We emphasise that it is not intended as a predictive model for validation of measured data, but to facilitate and assess design changes efficiently. While the individual components of the mathematical and numerical model are not new in separation, their holistic combination and coupling with our random rainfall delivered by a Galton board are indeed novel. A plan view of a sketch of Wetropolis is given in Fig. 2. The river–canal combination established is inspired by the River Aire and Leeds–Liverpool canal sharing a large part of the same river valley with the canal allowing some minor flood alleviation via (manual) flood control. Furthermore, implicit motivation may stem from the fact that Onno Bokhove and Wout Zweers are Dutch citizens and from the context of the Dutch Deltaworks, of which smallscale test versions were built in the “NoordOost” (NorthEast) polder after the 1953 North Sea flood; conceptual modelling of river flood components like in Wetropolis is perhaps natural for Dutch engineers and designers. See also (online) literature regarding the “Waterloopkundig laboratorium” in the “NoordOost” polder.
This introduction has given both an anecdotal and scientific background to Wetropolis' inception. Our intention here is to document its journey from design to outreach, and in doing so to inspire and enable readers to redesign Wetropolis bespoke to their own local catchment characteristics; the remainder of the article has the following outline. The above unknowns were determined via visual optimisation of simulations of an idealised mathematical and numerical model of Wetropolis before any design and construction of the tabletop setup were undertaken. This mathematical and numerical modelling is therefore explained in detail and used to determine (some of) the design unknowns in Sect. 2. The resulting tabletop design of the Wetropolis Flood Demonstrator is disseminated in Sect. 3. Our experience in demonstrating Wetropolis to the general public and to flood practitioners is summarised in Sect. 4, including the a priori surprising outcome that professionals in flood prediction and mitigation have also been inspired by Wetropolis, despite our primary aims having been public engagement. A discussion is found in Sect. 5, in which we outline some new designs and future directions, as gathered from such publicengagement sessions.
The mathematical model of Wetropolis comprises random rainfall and space–time continuous hydraulic modelling of interconnected river flow, reservoir and canallevel changes, as well as groundwater flow in the moor. While the individual modelling elements in separation are known or straightforward, their holistic combination with the statistical rain modelling as well as the subtle massconserving coupling between the elements is nontrivial and new. In addition, dissemination of the model is also required to facilitate adaptations by the readers. One other reason to be quite pedagogical is to reach a wider readership of enlightened and interested members of the public, especially educators. Subsequently, we will establish a numerical discretisation of this space–time continuous model and use numerical simulations to determine the a priori unknown parameters of rainfall amount r_{0} and wd. Other parameters will be determined heuristically in order to obtain desirable and practical dimensions of the experimental setup. We present the model here completely; full details of the numerical discretisation can be found in Appendix A.
2.1 Statistical modelling of randomised rainfall
As discussed in the introduction, rainfall is modelled stochastically via the outcome of draws from two Galton boards. In the mathematical model these outcomes are simulated, while in the physical flood demonstrator we have either used two actual Galton boards with two steel balls or one Galton board with one steel ball running through two consecutive Galton board channels. We have discretised rainfall into two categories: location and rain amount, per Wetropolis day. Rain location has four outcomes: rain in reservoir, moor and reservoir, moor, or no rain in the catchment with a discrete distribution of (3, 7, 5, 1)∕16. Independently, rain amount has per location four outcomes (1, 2, 4, 9)r_{0}∕wd, with again the discrete distribution (3, 7, 5, 1)∕16; hence, there are $\mathrm{4}\times \mathrm{4}=\mathrm{16}$ outcomes determined as a direct product of these two independent distributions, given in Table 1. The possible rain amounts per Wetropolis day are therefore (0, 1, 2, 4, 8, 9, 18)r_{0} and the value r_{0} will be determined in the subsequent modelling such that there is no flooding for (1, 2, 4)r_{0}∕wd rainfall, with potentially limited flooding for (8, 9)r_{0}∕wd and generally major flooding in the city plain for 18r_{0}∕wd. The resulting distribution of the rainfall per Wetropolis day in the “Wetropolis catchment” (to be read with Table 1) therefore becomes
A pdf of this discrete distribution for a computer trial over 500 wd's is shown in Fig. 3 (blue bars) with the expected values overlaid (crosses). Suitable values for r_{0} and wd are established by further modelling, described next.
2.2 Mathematical modelling in space–time
The key components of Wetropolis are a river channel with a onesided flood plain, a groundwater moor, a reservoir, and canals with three segments separated by lock weirs. The canal flows into the river in the city plain which lies at the downstream end of the setup. We refer to Fig. 2 for a plan view locating these elements in the actual tabletop experiments. In the original design model, the locations of the reservoir and moor have been swapped, and we have used a shorter river channel. Simplified mathematical submodels of these different elements are derived next in isolation before being coupled into one complete and novel mathematical model of rainfall and flooding via suitable boundary and interface conditions.
2.2.1 River dynamics
River flow is often modelled as 1D flow in a channel with a cross section $A=A(s,t)$ as a function of space, with a horizontal curvilinear coordinate s∈[0, L] along a winding river channel of length L, and time t. Both $A(s,t)=A(h;b(s\left)\right)$ and the in situ water depth $h=h(s,t)$, above a fixed river bottom b=b(s), and the mean flow velocities u(s,t) are all averaged over the cross section of the river. Hence, the river bottom lies at z=b(s), with vertical coordinate z and the river surface at $z=b\left(s\right)+h(s,t)$. For example, for a rectangular river channel we have A=w_{r}h with fixed width w_{r}, but in general A(h;s) can depend in a complicated fashion on the river depth h and directly on s – the latter via b=b(s). The governing equations are the SaintVenant equations (e.g. Bates et al., 2010) consisting of continuity and momentum equations, s_{m}, s_{res} and canal1 location L_{1c} along the river where the water flows from the moor and reservoir into the river, as well as a parameterisation of the channel friction, i.e.
with the phenomenological Manning friction coefficient C_{m}∈[0.01, 0.15] m${}^{\mathrm{1}/\mathrm{3}}$ s, cf. Munson et al. (2005), the hydraulic radius $R\left(h\right)=\frac{\mathrm{wet}\phantom{\rule{0.25em}{0ex}}\mathrm{area}}{\mathrm{wetted}\phantom{\rule{0.25em}{0ex}}\mathrm{perimeter}}$ (in metres), acceleration due to gravity g, and the discharge rates Q_{m}, Q_{res}, and Q_{1c} of the water flows from the moor, reservoir and canal into the river at locations s_{m}, s_{res}, and L_{1c}. Partial derivatives have been denoted by ${\partial}_{s}=\partial /\partial s$ and ${\partial}_{t}=\partial /\partial t$. These discharge rates are modelled as point sources using delta functions δ(s−s_{m}), δ(s−s_{res}) and δ(s−L_{1c}). In reality the same inflow rates will occur along finitelength, short strips along the river, centred around s_{m}, s_{res} and L_{1c}. The two unknown fields are A and u, with $h=h(A;s)$ an (often) implicit relation at every location s. For the above example with a river channel of a rectangular cross section, we find $R\left(h\right)=\frac{{w}_{\mathrm{r}}h}{\mathrm{2}h+{w}_{\mathrm{r}}}$ and $h=A/{w}_{\mathrm{r}}$. Initial conditions A(s,0) and u(s,0) have to be imposed at t=0 as well as boundary conditions A(0,t), u(0,t) and A(L,t), u(L,t) at s=0 and s=L. These latter boundary conditions are used (partially) according to the way the characteristics at $s=\mathrm{0},L$ of the hyperbolic equations (Eq. 2) determine whether the boundary data are flowing into the domain or not, cf. Toro (2001).
Even though the actual winding river channel with its onesided flood plain and city plain has a varying cross section, for the design calculations we made the simplification to model only a rectangular river channel. In addition, a zerothorder kinematic model approximation to the SaintVenant equations (Eq. 2) has been used for the limiting case with positive velocity u>0 and a constant slope ${\partial}_{s}b>\mathrm{0}$ of the river channel. To zeroth order, we assume that the bed slope and friction are locally in balance, i.e. the underlined terms in Eq. (2) are negligible, such that we obtain
This is a classical approximation used in hydraulics (Munson et al., 2005): the river flow is thus modelled by a kinematic or scalar hyperbolic equation in A, as follows:
with an upwind information speed $\mathrm{d}{Q}_{\mathrm{f}}\left(A\right)/\mathrm{d}A>\mathrm{0}$ for flux Q(A)=Au and inflow $A\left(h\right(\mathrm{0},t);s=\mathrm{0})$. Note that the flux Q_{f}=Q_{f}(A) is an implicit function of A since h=h(A). For u>0 with A=w_{r}h, it is a kinematic model or nonlinear, scalar conservation law in the water depth h given by
with the flux Q_{f}=Q_{f}(h) rewritten in terms of h, initial condition h(s,0) and upstream influx ${Q}_{\mathrm{0}}(\mathrm{0},t)={w}_{\mathrm{r}}h(\mathrm{0},t)u(\mathrm{0},t)$ defining h(0,t) since u is expressed in terms of h through Eq. (3). The SaintVenant equations are more advanced than the above kinematic model and allow both subcritical and supercritical flows. An interim and better model arises when we instead of Eq. (3) use the following balance:
which, after substitution into the continuity equation, Eq. (2a), yields an advection–diffusion equation; cf. Bates et al. (2010). Both these more advanced models are not required for the design estimates, but when one wishes to perform (more) accurate predictions of the hydrodynamics in Wetropolis, then such advanced models may be required.
2.2.2 Groundwater flow
Groundwater levels after rainfall are made visible in a transparent and elongated rectangular box filled to a high level with porous small lava rocks. The box is open at the top and one side and has walls at the remaining three sides and the bottom. Rain falls uniformly along this box via a copper pipe with a series of equidistant holes. The groundwater dynamics are modelled to zeroth order by assuming that the surface of the grains is flat, the rainfall uniform per surface area, that there is no surface runoff and that the fallen rainwater infiltrates sufficiently fast to contribute instantly to the groundwater level h_{m}(y,t) with coordinate y in a different direction, locally orthogonal to s at the location s_{m} where the groundwater flows into the river; cf. the delta function in the continuity equation (Eq. 2a) and kinematic equation (Eq. 5). A depthaveraged groundwater model with level h_{m}(y,t) from Barenblatt (1996) is used in a cell of width w_{v} and length L_{y}, e.g. w_{v}=0.095 m and coordinate $y\in [\mathrm{0},{L}_{y}=\mathrm{0.925}$ m]. The nonlinear diffusion equation for the groundwater level h_{m}(y,t), taken to be uniform in the lateral direction, is
with moor rainfall ${R}_{\mathrm{m}}(y,t)={R}_{\mathrm{m}}\left(t\right)$, porosity ${m}_{\mathrm{por}}\in [\mathrm{0.1},\mathrm{0.3}]$, the fraction ${\mathit{\sigma}}_{e}\in [\mathrm{0.5},\mathrm{1}]$ of pores filled with water, $\mathit{\alpha}=k/\left(\mathit{\nu}{m}_{\mathrm{por}}{\mathit{\sigma}}_{e}\right)$ with permeability $k={\mathrm{10}}^{\mathrm{8}}$ m^{2}, and viscosity $\mathit{\nu}={\mathrm{10}}^{\mathrm{6}}$ m^{2} s^{−1}. The boundary conditions are no flux through the wall at y=L_{y} such that ${\partial}_{y}h{}_{y={L}_{y}}=\mathrm{0}$, while at y=0 the moor is held at the level h_{3c}(t) of canal3, the upstream branch of the canal running in parallel to the river; i.e. this is a timedependent Dirichlet condition ${h}_{\mathrm{m}}(\mathrm{0},t)={h}_{\mathrm{3}\mathrm{c}}\left(t\right)$. The mass flux of moor water running in the river (at $s={s}_{\mathrm{m}}=\mathrm{2.038}$ m) is
where Q_{tm} is the mass flux following from integration of Eq. (7) over the domain $y\in [\mathrm{0},{L}_{y}]$ and $\mathrm{0}<\mathit{\gamma}<\mathrm{1}$ the fraction of moor water entering into the river. The reason to multiply by m_{por}σ_{e} is that the water volume in the matrix of particles in the moor changes suddenly from a space filled with pores into free space.
2.2.3 Reservoir
The reservoir is a rectangular box of dimensions ${h}_{\mathrm{res}}\times {w}_{\mathrm{res}}\times {L}_{\mathrm{res}}$ with timedependent water level h_{res}(t), e.g. L_{res}=0.293 m and w_{res}=0.123 m. In the physical model, the random rainfall enters either via a pipe or a long pipe with numerous holes visualising the rainfall, and it leaves the reservoir via an overflow pipe, here modelled simply as a straight weir; cf. Munson et al. (2005). Overflow of the reservoir once it is overfilled is not modelled. The reservoirlevel dynamics are governed by
in which L_{res}w_{res} is the area of the reservoir such that L_{res}w_{res}h_{res} is its timedependent volume, P_{wr} the overflow height of the weir, R_{res}(t) the reservoir rainfall, and Q_{res} the flux down into the river. Note that the coefficient C_{f} is dimensionless. The weir is located at $s={s}_{\mathrm{res}}=\mathrm{0.932}$ m, where water flows into the river; cf. the delta function in the continuity equation (Eq. 2a) and kinematic equation (Eq. 5).
2.2.4 Canal sections
One canal of uniform width w_{c} runs alongside the river, cf. the Leeds–Liverpool canal and the River Aire in central Leeds, in which the three canal sections are separated by lock weirs and have mean, timedependent water depths h_{1c}(t), h_{2c}(t), and h_{3c}(t). We are thus ignoring currents and height changes along the separate canal sections. Each canal section has a certain depth and is separated from the river by a berm. Canal3 is the highest and is blocked off on one end, at s=0, and has a weir located at $s={L}_{\mathrm{3}\mathrm{c}}=\mathrm{1.724}$ m. Its level is modelled as the variation of its volume due to partial inflow from the moor and outflow of water via a weir in canal2:
with weir height P_{3w}. Canal2 resides from $s\in [{L}_{\mathrm{3}\mathrm{c}},{L}_{\mathrm{2}\mathrm{c}}]$ with $s={L}_{\mathrm{2}\mathrm{c}}=\mathrm{3.608}$ m and is modelled likewise, but with inflow Q_{3c} from canal3 and outflow Q_{2c} into canal1:
The section of canal1 runs from $s\in [{L}_{\mathrm{2}\mathrm{c}},{L}_{\mathrm{1}\mathrm{c}}]$ with L_{1c}=3.858 m, width w_{c} and depth h_{1c}(t). It is modelled in the same manner with inflow from canal2 and outflow Q_{1c} into the river, as follows:
The weir at s=L_{1c} where water flows into the river is assumed to be subcritical; i.e. we assume there is a sufficient drop from canal1 to the river level. In terms of height levels, canal3 has a berm at z=0.06 m and its bottom resides at z=0.04 m; canal2 has a berm at z=0.04 m and its bottom resides at z=0.02 m, and canal1 has a berm at z=0.021 m and its bottom resides at z=0.001 m. To wit, the outflow at the two weirs into canal2 and canal1 is based on Bernoulli's relation and flow criticality; cf. Munson et al. (2005). At s=L_{1c}, e.g. for subcritical flow with flow depth h_{2c} and flow speeds V_{2c}≈0 upstream as well as critical flow V_{c} of height h_{c} over the weir, we therefore derive the following:
with ${C}_{\mathrm{f}}=(\mathrm{2}/\mathrm{3}{)}^{\mathrm{3}/\mathrm{2}}$. Similar derivations with suitable adaptations of the quantities involved determine the fluxes Q_{1c}, Q_{3c}, and Q_{res} over the other weirs.
2.2.5 Fully coupled system
When all of the above models for the individual components are combined, we obtain the entire coupled model, including its initial and boundary conditions, for the unknowns h(s,t), h_{m}(y,t), h_{res}(t), h_{1c}(t), h_{2c}(t), and h_{3c}(t).
with
as well as timedependent rainfall functions R_{res}(t),R_{m}(t), and upstream inflow Q_{0}(t). The remaining parameters are constants, with units and typical values listed in Table 2. The rainfall functions are defined such that, in the absence of other effects, e.g. unit porosity in the moor, they directly lead to a linear increase in the moor's groundwater level and the reservoir depth. Note that the mass fluxes involving the different Q^{′}s between the various systems are consistent. A space–time numerical discretisation of Eq. (14) is given in Appendix A. It involves a secondorder finitedifference approximation of the groundwater equation (Eq. 14b) in y, a firstorder finitevolume discretisation of the river equation (Eq. 14a) in s, and straightforward firstorder forwardEuler time discretisations of the time derivatives involved.
The rainfall functions are constant during every Wetropolis day and generally vary from Wetropolis day to day. On a given Wetropolis day,
in which ${n}_{r}=\mathrm{1},\mathrm{2},\mathrm{4},\mathrm{9}$ is drawn daily with probability $(\mathrm{3},\mathrm{7},\mathrm{5},\mathrm{1})/\mathrm{16}$ via one Galton board, while one of the combinations,
is drawn daily with probabilities $(\mathrm{3},\mathrm{7},\mathrm{5},\mathrm{1})/\mathrm{16}$ via the other Galton board, as explained in Sect. 2.1. The rainfall speed r_{0} will be determined by trialanderror and has the units of ∂_{t}h, i.e. m s^{−1}. Hence, the volumetric rate of rainfall per Wetropolis day on the moor for unit n_{r}=1 can then be calculated, yielding
2.3 Numerical results
Given the choice of parameter values with (or near) the values given in Table 2, the goal is to determine a suitable rainfall speed r_{0} and duration wd via trialanderror through numerical simulation. As initial conditions we take $h(s,\mathrm{0})=\mathrm{0.0135}$ m, ${h}_{\mathrm{m}}(y,\mathrm{0})=\mathrm{0}$, zero canal, and reservoir levels ${h}_{\mathrm{10}}={h}_{\mathrm{20}}={h}_{\mathrm{30}}={h}_{\mathrm{res}\mathrm{0}}=\mathrm{0}$ and an upstream influx of river water corresponding to the mass flux ${Q}_{\mathrm{0}}={Q}_{\mathrm{f}}\left(h\right(\mathrm{0},\mathrm{0}\left)\right)$ associated with h(0,0). In reality, rainfall will be varied daily by changing the action of two pumps, which require a fraction of a second to change gear. For someone viewing Wetropolis, a length of day between 5 and 20 s seems reasonable, so we have chosen wd=10 s as a first guess. In the design phase, values of r_{0} have been chosen and tuned in simulations of 100 to 500 wd's, i.e. 1000 to 5000 s, which can be simulated in about 10 % of real time. Hence, with the choice wd=10 s, the return period of extreme rainfall in Wetropolis is therefore (256∕7)10 s = 6:06 min, and the chance of 2 consecutive days of extreme rainfall has a return period of (256∕7)^{2}10 s = 223 min ≈ 3:43 h.
To monitor whether r_{0} has the (approximately) desired value during a simulation, major flooding is defined as occurring when the river level significantly, i.e. by circa 0.01 m or more, exceeds the canal1 berm along the strip of river bordering the city plain. This is monitored visually in daily snapshots, one of which is given in Fig. 4b. It contains a compound of levels to enable this flood monitoring, which requires some explanation. While the canal water enters the river in the city, for simplicity flood waters of the river are not modelled numerically to enter the canal or city, which suffices for our design purposes. The information displayed in Fig. 4b is as follows. The vertical axis has units in metres, so the range across the length L=4.21 of the setup is about 0.06 m, with the horizontal s axis lying along the river and canal. The zero level of the canals and rivers in the vertical is put at the river exit $(s,z)=(L,\mathrm{0})$ with vertical coordinate z. In reality, the lengths of river and canal differ slightly owing to the curved channels; for design and our numerical model, this difference is considered negligible, and it is sufficient to assume they are of the same length. The lower, solid, and thick red line displays the fixed sloped bottom of the river (with 1 % downhill gradient, i.e. ${\partial}_{s}b=\mathrm{0.01}$). The thinner solidblack line displays the river level with the upstream input depth of $h(\mathrm{0},t)={h}_{\mathrm{0}}\left(t\right)$. This water profile is generally adjusting dynamically to become uniform except at the reservoir influx (s_{res}=0.932 m), the moor influx (s_{m}=2.038 m), and the canal1 influx ($s={L}_{\mathrm{1}\mathrm{c}}=\mathrm{3.858}$ m), which cause sudden increases in the river levels. These larger and smaller jumps are indeed visible and identifiable in Fig. 4b. The flux into the river from canal1 is comparatively small, so the rise in the river level here is much smaller than the timevarying influx of water from the reservoir and moor. The bottom of the three canal sections is displayed by the (stepped) thickblue dashed line, with the upper canal3 level at $z={b}_{\mathrm{3}}=\mathrm{0.04}$ m, the middle canal2 level at $z={b}_{\mathrm{2}}=\mathrm{0.02}$ m, and the lower canal1 level at $z={b}_{\mathrm{1}}=\mathrm{0.0}$ m. The three berm (or dike) heights are 0.02 m higher at $\mathit{\{}{d}_{\mathrm{3}},{d}_{\mathrm{2}},{d}_{\mathrm{1}}\mathit{\}}=\mathit{\{}\mathrm{0.06},\mathrm{0.04},\mathrm{0.02}\mathit{\}}$ m. Canal berm or dike levels are displayed with a (stepped) thick solidred line, while the three varying canal levels are displayed as the (stepped) solidblue line. Steps in the berms occur where the weirs are placed and the jumps in the varying canal levels are determined by the hydraulic weir relations at these weirs, cf. Eq. (13). Some river flooding can occur when the river level, the (stepped) solidblack line, exceeds the canal2 berm downstream of the second weir, as is visible in Fig. 4b. Major flooding is defined when the river levels exceeds the canal1 berm in the city section, i.e. at $s={L}_{\mathrm{1}\mathrm{c}}=\mathrm{3.858}$ m the water depth h(L_{1c},t) significantly exceeds 0.02 m, visible as a snapshot in Fig. 4b, where the solidblack line of the river level is seen to exceed the solidred canal1 berm level downstream of the last weir around s=3.7 m. Via visual optimisation, i.e. monitoring when major flooding occurred in the city for the extreme or rare events of 90 % rainfall in both the reservoir and moor, a suitable value of the rainfall speed is found to lie around the value
The corresponding water volumes for the various Galton board outputs required in the moor and reservoir per Wetropolis day are then
Consequently, the pump supplying the rainfall on the moor should have a maximum discharge of about 324 mL s^{−1}, which is a manageable amount from a design perspective. Such a discharge is feasible by using inexpensive offtheshelf aquarium pumps, both for the supply of the upstream river influx Q_{0} and the varying rainfall amounts in reservoir and moor.
An example of a simulation over 500 wd's is summarised in Fig. 5. Since reservoir and canals are empty at t=0, we observe that it takes about 25 wd's before they are filled. During this time major flooding is lessened, or prevented completely, because the reservoir and canal in essence act as floodattenuation storage sites, supplying passive flood control. Extreme rainfall in this startup period tends to be buffered such that city flooding is prevented. Reservoir and canal levels are displayed in Fig. 5a versus time. The (constant) upstream river level h(0,t) and city river level h(L_{1c},t) are displayed as a function of time t in Fig. 5b, in which extreme events with n_{r}=18 are clearly identifiable as flood peaks at time t when $h({L}_{\mathrm{1}\mathrm{c}},t)>\mathrm{0.02}$ m. A snapshot of the groundwater level h_{m}(y,t) in the moor is displayed in Fig. 5c; it shows the noflux upstream boundary condition at y=L_{y} and the gradual decrease in the groundwater level towards its outflow location at y=0. The rain units for the moor (being 1, 2, 4, or 9), reservoir (being 1, 2, 4, or 9), and their summation n_{r}, with the discrete values of 1, 2, 4, 8, 9, or 18, are displayed in Fig. 5d. The peaks of extreme rainfall are, of course, seen to match the peaks in extreme flooding in Fig. 5b and d except, possibly, during the first circa 25 wd's, when the canals and reservoir tend to act as floodattenuation buffers.
After the design calculations commenced on 29 May 2016 and were completed on 8 June 2016^{3}, the Wetropolis Flood Demonstrator was constructed and finalised between 4 June and 31 August by Onno Bokhove and Wout Zweers^{4}. The final design was limited by the demand to transport it in the back of a car. We note that the reservoir and moor have been swapped in the actual setup, relative to the mathematical design, and that the riverchannel length has been increased to 5.2 m.
The physical Wetropolis model consists of several elements which we describe next.

The topographic landscape with a winding river channel, onesided flood plain, canals, and the city plain has been routed out of two standard polystyrene foam plates each of dimension $\mathrm{5}\times \mathrm{60}\times \mathrm{120}$ cm^{3} (plus a small extra foam plate) with an overlay to fit two plates together. A smaller third piece was added to extend the river length after the city which enhanced flooding in the city plain. An overview is given in Fig. 2 and a photograph is found in Fig. 6. Drawings have been made in computeraided design (CAD) programs Rhino and Grasshopper and used to steer the router. Routing precision is circa 0.8 mm. After the routing, the river channel and its flood plain have been roughened by varnishing fine sand to the base.

A framework of wooden support slabs has been made that fits on four Aframes. This framework is put together with a bolt–nut system such that it can be disassembled for transport. Wooden wedges are used to squeeze and level the foam pieces within the slab framework. The three foam pieces are squeezed together to limit leakage. Aluminium onesidesticky tape is used locally to seal two sections of the river channel together and thus bridge two adjacent channels. Rather than sealing off all leakage, which in practice becomes impossible, an “aquifer” system of two interconnected gutters underneath the seams of the foam pieces leads leaked water back to the holding reservoir with the three aquarium pumps. Hence, the water budget is closed in the absence of evaporation, the latter of which is negligible on the timescale of operation (typically a few hours) but not on the timescale of a day or more.

Three aquarium pumps with a maximum pumping capacity of 0.375 L s^{−1} are placed in a holding reservoir, a rectangular bucket with dimensions $\sim \mathrm{0.3}\times \mathrm{0.40}\times \mathrm{0.22}$ m^{3}, which is hung underneath the wooden framework in a rectangular area adjacent to the upstream inflow point of river water. Plastic tubing with inner and outer diameters of circa (1.8, 2.2) mm leads the water from this reservoir to the upstream point and, depending on whether it rains or not, to the reservoir and the moor. Rainfall on the moor is spread out and visualised using a copper pipe with numerous downwardfacing holes over the lava grit.

The moor unit is made of acrylic and on the open side face of the box a gauze prevents the lava grit from avalanching into the river. The acrylic reservoir is open from the top and water can enter through a hole near the top edge. Outflow of water in the river is regulated via an internal pipe whose outflow level can be manually adjusted. Hence, active flood control can be demonstrated by manually adjusting this outflow level. Outflow into the canal can be arranged separately via an adjustable valve. Note that this is slightly different from the setup in the mathematical and numerical design model, where the outflow of moor water was partitioned between the river and canal.

The setup for the Galton boards, including the Galton boards themselves, the accompanying Arduino control units, power sockets, and plugs to operate the three aquarium pumps. The two draws from the discrete probability distribution are either computer generated or determined from the random paths of (a) steel ball(s) through the asymmetric Galton boards. In the latter case, the steel ball triggers a signal by interrupting optical sensors in one of the four channels on each Galton board, cf. Fig. 1. The signal subsequently steers either the reservoir pump, moor pump, both, or none as arranged via the Arduino technology.
Further specifications have been provided in Appendix B.
Illustrative images of Wetropolis in action are shown in the photographs of Figs. 7 and 8. It includes closeups of excessively flooded river bends and the city plain, the reservoir, and its outlets as well as the moor under heavy rainfall in Fig. 7. During extremely heavy rainfall (90 % per Wetropolis day) after a relatively wet period, the moor becomes supersaturated and the groundwater level can rise through the lava grit and trigger fast surface runoff. In other situations the groundwater level is below the surface of the lava grit. Under varying rainfall the rising and falling groundwater level can be observed through the transparent acrylic walls. This visualisation was inspired by cartoons in hydrological textbooks, in which we often find similar cross sections of the earth and its groundwater levels.
To date Wetropolis Flood Demonstrator has been showcased to flood victims at public events, attendees of science fairs, and scientists and flood professionals at bespoke workshops on natural hazards. Wetropolis has been showcased^{5}:

first at the general assembly of EPSRC UK network “Maths Foresees” in Edinburgh in September 2016 to an audience of scientists with expertise in environmental fluid dynamics and representatives from stakeholders such as the Met Office, JBA Trust, and the Environment Agency; Wetropolis has been created as an outreach project within this Maths Foresees network;

at the Churchtown Flood Action Group conference in January 2017 in Lancashire to circa 140 flood victims and several experts on flooding, after a public lecture by Onno Bokhove on the statistics of extremes; see Fig. 8;

two afternoons as part of the public exhibition on the Boxing Day 2015 floods in Leeds' Armley Industrial Museum in December and March 2017 to an audience of flood victims, family, and friends, including children;

at a bespoke Canal and River Trust workshop in Liverpool in March 2017 by students of the Leeds Centre for Doctoral Training in Fluid Dynamics;

at the “Be Curious” science festival in March 2017 at the University of Leeds for a general audience;

to scientists at the Maths Foresees “Environmental Modelling in Industry” study group at the Newton Gateway to Mathematics, Cambridge, in April 2017, which particular event focussed on solving mathematical challenges related to flooding; see https://gateway.newton.ac.uk/event/tgmw41 (last access: 4 May 2020); and

at the Yorkshire iCASP confluence (integrated catchment program) in June 2018 for a range of scientists, flood professionals, stakeholders and politicians; see Fig. 1b.
Based on interactions during the above events with audience members and organisers, as well as on formal discussion during workshops regarding public engagement, we consider some of Wetropolis' strengths and weaknesses. Most of these considerations are anecdotal, except for the discussions at Maths Foresees meetings and the Newton Gateway to Mathematics, which were based on formal notes of the indepth round table and workshop discussions as well as the studygroup host. We stress, therefore, that the outreach component of this study is lacking a proper scientific method (from a socialscience perspective). Indeed, none of these discussions concerns formal questionnaires, wellbalanced questions and subsequent statistical analysis, which would constitute a formal investigation of the feedback and Wetropolis' impact. There are two reasons why such a formal analysis is lacking: the authors do not have the expertise to undertake such an analysis and, more importantly, during the showcasings for atthetime recent flood victims we did not want to further bother these victims by conducting intrusive questionnaires in a potentially unhelpful manner. However, this is something to be considered for future demonstrations, in collaboration with the necessary experts. It also implies that the conclusions suggested below are preliminary in nature.
The strength of Wetropolis is that it is a physical visualisation of the probability of extreme rainfall and flooding events including actual and visual river hydraulics, groundwaterlevel changes, and interactive flow control. We recall that the reservoir has valves such that the audience can store and release water interactively into the canal and river in order to control and possibly prevent flooding in the city. Wetropolis is, however, a conceptual model of flooding rather than a literal scale model of a specific catchment. It has, however, been inspired by the Boxing Day floods of 2015 of the River Aire, in and upstream of Leeds, UK. This conceptuality is both a strength and weakness because one needs to explain the translation of a 1:200year return period for a realistic extreme flooding and rainfall event such as the Boxing Day 2015 flood of the River Aire into one in Wetropolis with its one in a 6:06 min return period, and one also needs to explain that the moor and reservoir are conceptual valleys where all the rain falls, since rain cannot fall everywhere in the Wetropolis catchment, in contrast to rainfall in real catchments. This scaling and translation step is part of the conceptualisation, which the audience, whether public or scientific, needs to grasp. The visualisations of flooding in the city and the groundwater level also involve learning steps. Hitherto, this conceptualisation step was either explained by the Wetropolis wardens in attendance at a demonstration, via our bespoke poster, or both. Alternatively, we aim to arrange bespoke audiovisual material.
Due to this learning curve, the most receptive public audiences have been flood victims or people with friends or family who went through the unpleasant and potentially devastating experience of being flooded (cf. two showcasings at Armley Musuem flood exhibition and the Churchtown Flood Action Group workshop). We have perceived such audiences to be the most receptive, inquisitive and interactive because they have an intrinsic interest in flooding phenomena and wish particular questions to be addressed in order to gain more understanding as to what causes flood hazards, how these hazards can be predicted and how such floods can possibly be tackled through flood mitigation and/or management. Recall that these were exactly the questions with which we started off in the introduction. In particular, Wetropolis aids in raising awareness of the probabilistic character and randomness of rainfall and flooding events, also in connection with the difficulties in predicting some of these extreme events. Combining showcasing Wetropolis with a general public lecture on the science of flooding has proven to be particularly successful, cf. Science of Floods (2016) and Potter (2016), owing to such a presentation whetting the appetite to view a scale model with rainfall and river flooding. While Wetropolis was designed as a public outreach project, the reception from flood practitioners and scientists working in environmental fluid dynamics has been surprisingly positive.
Finally, we recently showcased Wetropolis II (a new and improved demonstrator based on the same design principles presented here) as part of the (biannual) Mathematics of Planet Earth exhibition^{6} organised by the corresponding Centre for Doctoral Training at Imperial College in London for about 500 to 1000 visitors over 9 d, from 15 to 23 February 2020. In line with the directives of the organisers, the audience was encouraged to volunteer bespoke feedback on two postit boards, one for the larger exhibition and one specifically for Wetropolis. We received 10 feedback posts with positive feedback as well as suggestions, the latter ranging from (a) build more of these, (b) make an exhibition setup for tsunamis to (c) please add a fullfledged rain and river flow predictive model of Wetropolis and compare the two^{7}. While this feedback is still lacking an indepth statistical basis, given that the organisers felt that a formal questionnaire would be too intrusive, it does provide additional insights into the audience's perception of Wetropolis.
We have demonstrated how the Wetropolis demonstrator of extreme rainfall and floods was constructed after an efficient mathematical and numerical design enabled us to estimate the characteristic components of the envisioned setup. This efficient mathematical model was first presented as a coupled system of ordinary and partial differential equations which we subsequently solved numerically to define a nearoptimal design. While that mathematical model is close to a prediction model for river and groundwater flows in Wetropolis, due to its relatively minimalist nature and purpose to facilitate only the design, it is likely not quite sophisticated enough to make bona fide predictions. In a final modelling step, we determined the reasonable rainfall and flow rates through numerical simulation, on which rates we based the actual design and construction of Wetropolis. We highlight that we have provided all these design steps and made all software and design drawings fully available in order to facilitate design adaptations by the reader, to conceptualise their particular catchments of interest in their own Wetropolis. Several improvements and extensions of Wetropolis are under exploration, as follows:

to accentuate a flooding event in the city more prominently, e.g. by measuring the flood waters of each flood event via a separate drain into a measuring cup or by triggering flashlights in the city to light up by closing an electric circuit by the flood water and to go off when the circuit is broken; in Wetropolis II, we have added a drain in one half of the city, with its land lying somewhat lower and protected by a dike, such that the volume of each flood is visually captured in a beaker;

to visualise key principles of natural flood management (NFM) or more broadly naturebased solutions (NBS) (e.g. Hankin et al., 2017; Lane, 2017; Potter, 2016; Harrabin, 2016; Cabaneros et al., 2018; Bokhove et al., 2019, 2020) by visualising the effects of different riverbed roughness to slow down the flow in a cutout riverbed segment via removable riverchannel inserts and by including a porous upland catchment with various smallscale river channels and flowattenuation features to enhance water storage, the features of which can be manipulated by the audience (operational in Wetropolis II);

to include droughts by modification of the dry days visualised by drying up of a water supply pipeline from the moor to the city, e.g. given the extremely dry summers of 1976 and 2018 in Europe; and

to include climate change by making certain flood events more extreme (i.e. by varying the discrete probability distribution); climate change is included in Wetropolis II by the addition of an extra upstream reservoir and corresponding pump, which is synchronised with the moor operation through the Galton boards, while the primary reservoir has been moved downstream near the city plains; via a switch this extra water influx due to “climate change” can be switched on or off.
Droughts can be modelled by extending 1:4 of the dry days/periods, i.e. with their 1∕16 chance, to 4 or 5 dry days, while the other three outcomes are kept at 1 dry day, with suitable changes in the overall statistics. Since the Galton board outcome on rainfall amount is not used on dry days, its outcome can be recycled to randomly assign for the drought period as 1:4. Climate change can be modelled by extending extreme rainfall to 2 or 3 d for on average 1:3 extreme rainfall outcomes, or per another ratio, in a similar fashion. These drought or climate change adaptations can be included in the same physical setup via changes in the Arduino programming.
5.1 Games
One of the shortcomings of the current Wetropolis setup is that it lacks bespoke educational material and games. The game suggestions for the current Wetropolis setup, which have arisen in evaluation sessions during some of the various workshops mentioned above, are as follows.

The notion of (theoretical and sampled) pdfs can be developed, based on recording histograms of actual Galton board outcomes and comparison of these outcomes against the theoretical outcomes. Games can be created to determine whether the outcomes have been tampered with by human intervention, e.g. a game in which one team is allowed to trigger extreme flooding by tricking the electronic eyes by a finger or by purposely misaligning the Galton boards without telling the other team whether or not such unnatural interventions took place. While the audience generally likes massive flooding to occur more often in Wetropolis, by secretly triggering daily 90 % rainfall in both moor and reservoir, recording the Galton board outcomes would immediately reveal that such tampering is unrealistic in that it makes no rain and low to intermediate rainfall into rare rather than common events.

Building a game on flood prevention in the city by controlling the valves on the reservoir, say over 30 to 100 wd's, with the winning team having the least or zero amount of flooding in the city. This can also include a discussion on the possibility that the winning team can win by chance over a limited set of trials rather than by virtue of optimal flood control.

The audience can play with the setup. The setup can namely be modified by interchanging the two out of three locations where rainfall is random, changing the locations of the reservoir and moor, for example by bringing one unit in close proximity to the city, including an investigation as to what consequences these changes entail in observed spatiotemporal rainfall and flooding patterns.

A particular case can be chosen by changing a switch between normal, drought, and climate settings; and a game can be created for teams to determine which probability distribution is used from the observed outcomes. Given the inevitable bias in the analogue Galton boards and the nature of statistics, such a game will include some uncertainty. Finally, when the Galton boards fail, which happens occasionally if a steel ball balances exactly on a split, then an automatic routine takes over. Calculating the changes in the return periods for the droughts, extreme events, and superextreme climatechange events is an interesting yet straightforward (classroom) exercise.

Recently, it has been suggested to make a full predictive numerical model and include a few measurements with data assimilation in order to facilitate a live scientific display of the weather and river flow predictive capabilities.

Finally, a formal statistical evaluation of the response of the public to Wetropolis has not yet been undertaken.
Each of these suggestions requires further development.
5.2 New approaches in science and water management
Flood practitioners from various stakeholders have been quite positive about Wetropolis' novel way of visualising the probability of extreme rainfall and flooding events via the asymmetric Galton boards and how the outcomes from these boards directly lead to observable rainfall and river dynamics in the setup. Stakeholders such as JBA Trust and the Environment Agency see Wetropolis as a potentially useful tool to trigger discussions about innovations in flood mitigation and water management, as part of workshops and brainstorm sessions. To date, Wetropolis has triggered two innovations: one on the use of the revisited concept of floodexcess volume (FEV) in devising a novel and graphical costeffectiveness analysis to flood mitigation, in particular meant for decision makers, and one on education in water engineering and management.
Floodexcess volume (FEV) concerns the volume of flood waters that caused flood damage. It is the flood volume of the river flow beyond a certain, chosen, and relevant threshold water level at a certain target, the critical river location. This FEV, expressed in cubic metres (m^{3}), or expressed more visually and comprehensively as a square lake of 2 m depth with a certain side length, is a useful measure to devise floodmitigation strategies. It allows us to quantify what fraction of the FEV is mitigated by a certain strategy. Our costeffectiveness analysis (Bokhove et al., 2019, 2020) results in a series of squarelake graphs, one for each floodmitigation scenario envisioned, which expresses the flood volume mitigated by a particular floodmitigation measure, its cost, its cost per percent mitigated, and the overall costs. When an accumulation of floodmitigation measures captures the entire FEV, the FEV is essentially reduced to zero. Building higher flooddefence walls in a city at or just above the maximum river level to be mitigated does, for example, reduce the FEV to zero in one fell swoop. However, building high walls around a river in a city is only one type of floodmitigation scenario, one that may be undesirable, so in general floodmitigation scenarios, expressed visually in squarelake graphs, will consist of an accumulation of measures such as riverbed widening, i.e. givingroom to the river (GRR), active floodstorage plains, higher flood walls, and NBS. Each measure cuts a certain fraction as a rectangular or quadrilateral strip off the square flood lake, with accompanying costs displayed. The graphical costeffectiveness analysis has been developed in a series of papers by Bokhove et al. (2018a, b, 2019, 2020) for several extreme river floods in the UK and France in order to facilitate and improve evidencebased decisionmaking by city councils and citizens' groups.
We further note that Wetropolis has inspired a project on tangible models for education and water management in Enschede, the Netherlands, involving a consortium of Dutch SMEs, schools and universities – see https://www.wetropolis.nl (last access: 4 May 2020). It consists of parallel activities with an overarching Wetropolis theme, including

development of educational content and tools to raise awareness in primary schools and extend the Dutch GRR programme at secondary school levels; cf. Sect. 5.1;

buildup of “citizensensing” experiments, experiments to measure climaterelated indicators such as groundwater, drought, and temperature in urban areas; it is a relatively recent form of communitybased participatory environmental monitoring with success in the Netherlands, Spain, and Kosovo; cf. Woods et al. (2018); and

development of outreach models on the water cycle, drought, and heat stress phenomena in the local environment suitable for handson exploration in public settings such as museum exhibits.
Allowing for irregular time steps $\mathrm{\Delta}{t}_{n}={t}^{n+\mathrm{1}}{t}^{n}$ with t^{0}=0, the entire system (Eq. 14) has the following space–time discretisation, using regular finite differences for the groundwater equation, a firstorder finitevolume method with upwinding for the river equation, and a firstorder forwardEuler time discretisation for all differential equations involved, as follows (cf. Morton and Mayers, 2005; Leveque, 1990).
with
in which N_{y} is the number of regular grid points in the moor across L_{y}, $\mathrm{\Delta}y={L}_{y}/{N}_{y}$, and ${h}_{j}^{n}={h}_{m}(j\mathrm{\Delta}y$, t^{n}) and ${h}_{j+\mathrm{1}/\mathrm{2}}=({h}_{j+\mathrm{1}}+{h}_{j})/\mathrm{2}$; N_{x} for k=1, …, N_{x} is the number of finitevolume cells $\mathrm{\Delta}{x}_{k}={s}_{k+\mathrm{1}/\mathrm{2}}{s}_{k\mathrm{1}/\mathrm{2}}$ for the river with cell faces ${s}_{k\pm \mathrm{1}/\mathrm{2}}$ and cell average
and in which the Kronecker delta symbol is δ_{km}=1 for the cell k in which s_{m} resides and is zero elsewhere and likewise δ_{kr}=1 in the cell in which s_{res} resides, etc. A stable, explicit time step is determined using suitable CFL conditions based on the information speed and the nonlinear diffusion. )
A GitHub site contains information on the materials used, building instructions, as well as a historical timeline of its construction^{8}. Some design tools and materials are briefly outlined as follows.

Matlab programs of the numerical model are available on GitHub, concerning three versions. The third version, “tabletopt3v2019.m”, was used here, which equals “tabletopt2v2016.m”, except for some relabelling of figure axes.

Blue polystyrene plates were used, “Isolatieplaat polystyreen XPS” of dimensions $\mathrm{120}\times \mathrm{60}\times \mathrm{5}$ cm^{3} to route the terrain. Yacht varnish was mixed with fine calcinated and sieved sand and shells (with holes of 0.9 m and a wire thickness of 0.1 m). The following CAD programs were used: Solidworks for the designs, saved as a Step file for import in Rhino (V5); plugin in Rhino for routing: Rhinocam 5, which generates routing/NC files; and foams were routed on a BZT 1400 PF router/frees with a Winpcnc driver.

Aquarium pumps are used: Syncra 1.5, 234–240 V, 50 Hz, 23 W, $Q\mathrm{max}=\mathrm{1350}$ L h^{−1}, $H\mathrm{max}=\mathrm{1.8}$ m.

Design drawings for the topographic foam plates are given in Fig. B1 and photos of the wooden support system in Fig. B2.
Data sets are available by contacting the correspondence author and running the codes provided at the GitHub site given.
OB and WZ made the first calculations, designed the setup, and completed the construction. Further developments, including outreach and writing, were undertaken evenly by all the authors.
The authors declare that they have no conflict of interest.
The Wetropolis Flood Demonstrator originally started as an outreach project in EPSRC UK Living with Environmental Change (LWEC) network “Maths Foresees”, under grant EP/M008525/1 for Thomas Kent, Onno Bokhove, Tiffany Hicks, and the construction materials. The project is also partially stimulated by funding for Onno Bokhove and Thomas Kent from DARE project (Data Assimilation for the REsilient city) EP/P002331/1 and local EU project “Wetropolis decisionmaking, participation, education and outreach with physical and augmented hydrological models” for Wout Zweers. The evaluation of Wetropolis in Sects. 4 and 5 is based on formal discussions within two workshops of the Maths Foresees network, in 2016 in Edinburgh and 2018 in Leeds, as well as feedback from the Churchtown Flood Action Group, Leeds' Armley Museum, the Newton Gateway to Mathematics and the Mathematics of Planet Earth exhibition, which (anecdotal) feedback can be found on the Wetropolis' GitHub site. Finally, we thank Chris Skinner and two anonymous referees for their reviews, which helped improve our original manuscript.
This research has been supported by the EPSRC (grant no. EP/M008525/1) and the EPSRC (grant no. EP/P002331/1).
This paper was edited by Matjaz Mikos and reviewed by Christopher Skinner and two anonymous referees.
Barenblatt, G. I.: Scaling, selfsimilarity, and intermediate asymptotics, Cambridge University Press, Cambridge, 1996. a
Bates, P. D., Horritt, M. S., and Fewtrell, T. J.: A simple inertial formulation of the shallow water equations for efficient twodimensional flood inundation modelling, J. Hydrol., 387, 33–45, 2010. a, b
Bokhove, O., Kelmanson, M. A., and Kent, T.: On using floodexcess volume to assess natural flood management, exemplified for extreme 2007 and 2015 floods in Yorkshire, available at: https://eartharxiv.org/87z6w/ (last access: 4 May 2020), 2018a. a
Bokhove, O., Kelmanson, M. A., and Kent, T.: Using floodexcess volume to show that upscaling beaver dams for protection against extreme floods proves unrealistic, available at: https://eartharxiv.org/w9evx/ (last access: 4 May 2020), 2018b. a
Bokhove, O., Kelmanson, M. A., Kent, T., Piton, G., and Tacnet, J.M.: Communicating (naturebased) solutions using floodexcess volume for three UK and French river floods, River Res. Appl., 35, 1402–1414, 2019. a, b, c
Bokhove, O., Kelmanson, M. A., Kent, T., Piton, G., and Tacnet, J.M.: A costeffectiveness protocol for floodmitigation plans based on Leeds' Boxing Day 2015 floods, Water, 12, 652, 2020, original submnisison available at: https://eartharxiv.org/stc7r/, last access: 4 May 2020. a, b, c, d
Cabaneros, S., Danieli, F., Formetta, G., Gonzalez, R., Grinfeld, M., Hankin, B., Hewitt, I., Johnstone, T., Kamilova, A., Kovacs, A., Kretzschmar, A., Kiradjiev, K., Pegler, S., Sander, G., and Wong, C.: JBA Trust challenge: a riskbased analysis of small sale distributed “naturebased” flood risk management measures deployed on river networks, Report of Environmental Modelling in Industry Study Group, April 2017, of the UK EPSRC Living with Environmental Change Network “Maths Foresees”, Turing Gateway to Mathematics, Cambridge, available at: https://gateway.newton.ac.uk/event/tgmw41 (last access: 4 May 2020), 2018. a
Environment Agency: Hydrology of the December 2015 Flood in Yorkshire, in: April 2016 report of Environment Agency, available at: https://www.ceh.ac.uk/sites/default/files/20152016 Winter Floods report Low Res.pdf (last access: 4 May 2020), 2016. a, b, c
Hankin, B, Arnott, S., Whiteman, M., BurgessGamble, L., and Rose, S.: Environment Agency October, Working with Natural Processes – Using the evidence base to make the case for Natural Flood Management, available at: https://www.gov.uk/government/uploads/system/uploads/attachment_data/file/654435/Working_with_natural_processes_using_the_evidence_base.pdf (last access: 4 May 2020), 2017. a
Harrabin, R.: Pickering leaky dams flood prevention scheme `a success', BBC News Online, available at: http://www.bbc.co.uk/news/ukenglandyorknorthyorkshire36029197 (last access: 4 May 2020), 2016. a
Lane, S. N.: Natural flood management, WIREs Water, 4, 1211, https://doi.org/10.1002/wat2.1211, 2017. a
Leveque, R. L.: Numerical Methods for Conservation Laws, in: Lectures in Mathematics, Birkhäuser, Basel, 1990. a
Morton, K. W. and Mayers, D. F.: Numerical solution of partial differential equations, Cambridge University Press, Cambridge, 2005. a
Munson, B. R., Young, D. F., and Okiishi, T. H.: Fundamentals of Fluid Mechanics, Wiley, Hoboken, NJ, 2005. a, b, c, d
Pampaloni, M., Pileggi, T., and Caporali, E.: Awareness to the mitigation of hydraulic risk through LEGO models, Poster, in: European General Assembly 2018, 8–13 April 2018, Vienna, 2018. a
Potter, M.: Fighting the floods: community knowledge and engagement, presentation at the “Science of Floods workshop”, May 2016, Hebden Bridge, available at: https://www.pennineprospects.co.uk/file_download/752/2016+MIKE+POTTER+PICKERING+.pdf (last access: 4 May 2020), 2016. a, b
Sanderson, M.: Changes in the frequency of extreme rainfall events for selected towns and cities, Met Office “Ofwat” report, available at: https://www.ofwat.gov.uk/wpcontent/uploads/2015/11/rpt_com_met_rainfall.pdf (last access: 4 May 2020), 2010. a
Science of Floods: workshop May 2016, Hebden Bridge, available at: https://www.pennineprospects.co.uk/news/thescienceoffloods (last access: 4 May 2020), 2016. a, b
Stocker, T. F., Qin, D., Plattner, G.K., Tignor, M., Allen, S. K., Boschung, J., Nauels, A., Xia, Y., Bex, V., and Midgley, P. M. (Eds.): IPCC 2013: Summary for Policymakers, in: Climate Change 2013: The Physical Science Basis. Contribution of Working Group I to the Fifth Assessment Report of the Intergovernmental Panel on Climate Change, Cambridge University Press, Cambridge, UK, New York, NY, USA, 2013. a
Toro, E. F.: ShockCapturing Methods for FreeSurface Shallow Flows, John Wiley and Sons, New York, 2001. a
West Yorkshire combined Authority: Leeds city region flood review report, December 2016, available at: https://westyorkshire.moderngov.co.uk/documents/s9687/Item 8iv  Appendix 4  Leeds City Region Flood Review Report (last access: 4 May 2020), 2016. a
Wong, G., Maraun, D., Vrac, M., Widmann, M., Eden, J. M., and Kent, T.: Stochastic model output statistics for bias correcting and downscaling precipitation including extremes, J. Climate, 27, 6940–6959, 2014. a
Woods, M., Balestrini, M., Bejtullahu, S., Bocconi, S., Boerwinkel, G., Boonstra, M., Boschman, D.S., Camprodon, G., Coulson, S., Diez, T., Fazey, I., Hemment, D., van den Horn, C., Ilazi, T., JansenDings, I., Kresin, F., McQuillan, D., Nascimento, S., Pareschi, E., Polvora, A., Salaj, R., Scott, M., and Seiz, G.: Citizen Sensing: A Toolkit. Making Sense, https://doi.org/10.20933/100001112, 2018. a
https://www.youtube.com/watch?v=3yNoy4H2Zo (last access: 4 May 2020) and https://www.jbatrust.org/howwehelp/physicalmodels/wavetank/ (last access: 4 May 2020); this tank was designed by a team, including Onno Bokhove, of the Centre for Doctoral Training in Fluid Dynamics in Leeds, upon a request by JBA Trust.
In a recent personal communication with a Yorkshire flood expert, a 1:300year return period was mentioned, with no formal confirmation yet to date.
The first design and complete model calculations were presented during a seminar at Imperial College London on 1 June 2016. A week later an error in the use of the Manning coefficient C_{m} was fixed, leading to an increase in the riverchannel length L by a factor of 4. Hence, the winding channel.
See public postings in that period around 8 June and 31 August 2016 on https://www.facebook.com/resurging.flows (last access: 4 May 2020) and the design history on https://github.com/obokhove/wetropolis20162020 (last access: 4 May 2020).
For movie footage, see the posts dated 31 August 2016 (with an extremely rare Boxing Day 2015type flood after 2 consecutive days of extreme rainfall), 6 September 2016, 16 January 2017, 8 December 2016, and 7 April 2017 at https://www.facebook.com/resurging.flows (last access: 4 May 2020) or at https://youtu.be/1FIHFOn6IPQ (last access: 4 May 2020) and the Boxing Dayequivalent flood at https://youtu.be/N4Sp5gHXcz4 (last access: 4 May 2020).
See http://www.imperial.ac.uk/news/195539/planetearthaliveusingmathematicsunderstand/ (last access: 4 May 2020).
This feedback and the feedback of all exhibits are found at https://github.com/obokhove/wetropolis20162020/tree/master/feedback (last access: 4 May 2020).
https://github.com/obokhove/wetropolis20162020 (last access: 4 May 2020).
 Abstract
 Introduction
 Mathematical design
 Tabletop design
 Wetropolis illustrated and demonstrated
 Summary and discussion
 Appendix A: Numerical discretisation of the entire system
 Appendix B: Wetropolis' design details
 Data availability
 Author contributions
 Competing interests
 Acknowledgements
 Financial support
 Review statement
 References
 Abstract
 Introduction
 Mathematical design
 Tabletop design
 Wetropolis illustrated and demonstrated
 Summary and discussion
 Appendix A: Numerical discretisation of the entire system
 Appendix B: Wetropolis' design details
 Data availability
 Author contributions
 Competing interests
 Acknowledgements
 Financial support
 Review statement
 References
tabletopdemonstration model with extreme rainfall and flooding, including random rainfall, river flow, flood plains, an upland reservoir, a porous moor, and a city which can flood. It lets the viewer experience extreme rainfall and flood events in a physical model on reduced spatial and temporal scales with an event return period of 6.06 min rather than, say, 200 years. We disseminate its mathematical design and how it has been shown most prominently to over 500 flood victims.
tabletopdemonstration model with extreme rainfall and flooding, including...