Urban Flood Simulation Using MODCEL—An Alternative Quasi-2D Conceptual Model

: Urban ﬂood modelling has been evolving in recent years, due to computational facilities as well as to the possibility of obtaining detailed terrain data. Flood control techniques have also been evolving to integrate both urban ﬂood and urban planning issues. Land use control and ﬂow generation concerns, as well as a set of possible distributed measures favouring storage and inﬁltration over the watershed, also gained importance in ﬂood control projects, reinforcing the need to model the entire basin space. However, the use of 2D equations with highly detailed digital elevation models do not guarantee good results by their own. Urban geometry, including buildings shapes, walls, earth ﬁlls, and other structures may cause signiﬁcant interference on ﬂood paths. In this context, this paper presents an alternative urban ﬂood model, focusing on the system behaviour and its conceptual interpretation. Urban Flood Cell Model-MODCEL is a hydrological-hydrodynamic model proposed to represent a complex ﬂow network, with a set of relatively simple information, using average values to represent urban landscape through the ﬂow-cell concept. In this work, to illustrate model capabilities, MODCEL is benchmarked in a test proposed by the British Environmental Agency. Then, its capability to represent storm drains is veriﬁed using measured data and a comparison with Storm Water Management Model (SWMM). Finally, it is applied in a lowland area of the Venetian continental plains, representing ﬂoods in a complex setup at the city of Noale and in its surroundings.


Introduction
Urban floods are complex phenomena that usually require a mathematical model to support diagnosis and design procedures. The first mathematical flood flow models that were developed tended to focus mainly on the main flow paths or on the drainage network and its hydraulic structures. In this early period, computational limitations imposed the use of one-dimensional (1D) approach. Even nowadays, it may be a possible (and useful) approach, mainly when designing storm drains and channels, once it is not expected that overflows occur in the design process. However, if a drainage system is already failing and flooding a city portion, it is very difficult to work with a 1D branched network to make it representative of the phenomenon. In this case, the use imprecisions, especially if one relies exclusively in the mathematical model complexity as a guarantee of quality results.
Leandro et al. [1] also highlighted this problem, recognising that special key features found in urban areas, such as buildings and roads for example, can obstruct the expected natural flow paths and can cause sudden vertical drops, changes in the flow directions and local head losses. Thus, despite the fact that highly detailed Digital Elevation Models (DEM) may be available, containing detailed information about the morphology of such key features, they cannot be simply inserted into the 2D grid, without any further overland flow interpretations.
Considering this initial picture, this paper aims to present the use of the Urban Flood Cell Model-MODCEL, a mathematical model developed to couple with the urban diversity discussed in the previous lines, as an alternative for urban flood modelling. The early development of MODCEL started in a research about urban flood modelling developed at the Federal University of Rio de Janeiro, in Brazil. This research pointed out the importance of representing hydrologic surface processes and surface flows, jointly with storm drains, drainage channels and rivers. The main challenge was to put these demands together, giving primary importance to the physical interpretation, while maintaining data needs as simple as possible. MODCEL was developed as an interpretative conceptual model that uses a quasi-2D approach to represent physical reality to fulfil the proposed objectives of its construction. This model was presented in a first version in Mascarenhas and Miguez [8] and its evolution is shown here as an alternative tool for aiding in the process of urban water management and flood control design, especially when urban diversity and hydraulic structures play an important role in defining flood patterns.
MODCEL concept was based on the original work from Zanobetti and Lorgerè [9]. The flow cell concept was initially developed to overcome computation limitations when simulating large flood plains. However, due to a modular construction, this concept can be very useful, even nowadays, to simulate urban systems, where natural topography patterns, urban landscape, and hydraulic structures merge in a complex setting.
The modelling process partitions the watershed into a number of cells, which work as ponds (with storage capacity) interconnected by a vast array of hydraulic links, where the dynamic De Saint Venant equation occupies a prominent place.
In this work, MODCEL will be presented and applied, first in a benchmark test, then in a comparison with a storm drain model and, at last, in a complex urbanised flood plain, where the Italian city of Noale lies.

Urban Floods and Mathematical Modelling Importance
Urban floods are usually associated with a complex terrain set-up, where several structures take part. Urban environments may be responsible for a multitude of flow possibilities when storm drains fail. It is quite usual that water spilling out of the drainage system may cause inundation over vast areas and city structures may interact with hydraulic structures, composing an unplanned flow network that includes urban surfaces, mainly through streets acting as channels. At the same time, several undesired reservoirs are addicted to the system, with parks, public squares and buildings temporally retaining waters, without having been designed for this function (and in an undesirable way).
This situation usually demands the aid of a mathematical model to help the design team understand how the resulting system behaviour makes the drainage network interact with the urban landscape.
Another important and current aspect regarding urban drainage discussions refers to the changing paradigm associated with designing flood control alternatives. Burian and Edwards [10] reviewed the evolution of the drainage systems from 3000 BC to present days and pointed out interesting aspects:

•
The urban drainage systems often evolved through trial-and-error modifications after the systems were initially constructed; Water 2017, 9,445 4 of 28 • Changes in perspective of urban drainage in a city were most often caused by disease outbreaks, scientific discoveries, or technical advances in planning, design, and construction.
The traditional drainage design approach works to improve channel conveyance, which often falls on the canalisation or rectifications of the watercourses (focusing on adapting the drainage net to the generated discharges). This approach arouse during the industrial city development, when the drainage systems were designed to face sanitation problems, conveying jointly storm waters and wastewaters. The urban growth of the industrial city occurred with very few controls [11], leading to several urban infrastructure gaps. Although the hygienist concept related with the traditional urban drainage design was important to equate public health problems at that time, this approach showed to be unsustainable, especially when the urbanisation process became more intense. Urbanisation itself limits river canalisation enlargements. Once city growth increases flow generation, this approach tends to transfer problems to downstream reaches of the basin.
On the other hand, Andoh [12] reviewed the urban drainage and wastewater practices, stating the need for a shift from a reactive framework of the traditional approach, which acted on the consequences of the increasing flows, to a proactive-preventative approach, centred on sustainable principles and involving the social dimension. To do this, this author discusses the adoption of distributed facilities to attenuate and/or store and manage urban waters, dealing with the problem at its early stages, near the sources of flow generation.
In this context, in the last decades, several approaches were developed to better adjust flow patterns in space and time. Among these approaches, it is possible to mention (not being exhaustive): low impact development [13,14]; sustainable urban drainage systems [15][16][17][18]; and water sensitive cities [19][20][21]; among others, with similar objectives.
Thus, the new trends in urban drainage points to distributed measures on the watershed (not only in the drainage net) intending to manage flow generation and minimising impacts of urbanisation over natural flow patterns. Once again, mathematical models arise as important tools to help in the design process of flood control alternatives, allowing to act in a distributed way over the watershed and to combine effects of the proposed measures in space and time.

Urban Drainage Modelling
Urban floods, as previously mentioned, involve a great variety of flow patterns over a very complex surface, passing through different hydraulic structures. Waters spilled out of the system may flow on the streets as though they were channels. Depending on water levels reached on the streets, open spaces and even buildings may store and retain waters, acting as reservoirs. Walls may work as weirs, for example, dividing a certain portion of the urban area, with different flow depths and discontinuous flow paths at each side. In such cases, flows may occur quite differently from the original (confined) expected behaviour, picturing their own patterns, in various directions and deeply influenced by interactions with local structures.
One of the most important aspects to consider in flood modelling of urban environments is the way which the resistance caused by buildings or other structures is represented in the model [22].
In this context, a model capable to fulfil the requisites for urban floods simulation should be able to represent hydrologic and hydraulic processes in a distributed way, integrating different possibilities of superficial flows and their interactions with urban landscape and with the formal drainage system.
In this situation, if the superficial flow is mainly governed by topography and urban occupation patterns, it is difficult to adequately represent the flooding behaviour with strict 2D models. High imperviousness and short time of concentration are common elements of urban hydrology. Under such conditions, sewer drainage networks play a primary role in modern cities when transporting runoff during storm events [23]. However, 1D models focused on the drainage network will also miss the representation of important superficial psychical processes.
Djordjevic et al. [3] presented the concept of dual drainage modelling, providing a more realistic alternative for representing urban flooding. In this concept, the network of storm drains is able to Water 2017, 9, 445 5 of 28 interact with superficial flows over open spaces, through the streets or between houses, composing a high-varied model in space and time. Besides that, underground pipes, in this case, might present free or surcharged flows. Note that, in 1999, when this work was published, this was an evolution from the simple 1D modelling, bringing 2D needs into discussion.
Following this tendency, Nasello and Tucciarelli [24] proposed a dual multilevel urban drainage model, representing the system as a double network, formed by an upper network of open channels (the street gutters) and a lower network of closed conduits (the sewer pipes), mainly focusing on this interaction. This model considered that a significant part of the flow moves along streets, until reaching an available inlet. Two overlapping cells occupy the same space location, in different levels, representing the street and the sewer network at each inlet basin, connected by a vertical link. Mass and momentum conservation laws are applied to each pipe/channel section, and the diffusive approach is used. When the sewer is under pressure, the hydrodynamic behaviour of the lower layer connects with the upper one. Leon et al. [25], for example, proposed a model to be developed in four modules: (1) hydrology; (2) street flows; (3) flow interception at inlets; and (4) storm-sewer flows in urban areas. More recently, Chang et al. [23] developed a "novel approach" to represent flow interactions between storm sewer system and overland surface flows. Noh et al. [26] compared urban flood modelling approaches with results obtained in laboratory scale experiments for the interaction of manholes, sewer pipes and surface flows. These examples show that this matter is still a current modelling challenge.
Maksimovic and Prodanovic [27] stressed that modelling improvement of surface-subsurface interactions requires not only change in the modelling concepts, but also detailed spatial data resolution. This data refers to: land use, terrain elevation and other features affecting surface runoff, available in digital form at a desirable horizontal and vertical resolution. At the time of this publication, the authors stressed that this could be a limiting fact in evolving models, once the required information could be not always easy or cheap to obtain. In fact these limitations are being surpassed, although detailed information on terrain modelling is still not cheap. Leandro et al. [28] discussed 1D/1D and 1D/2D models, respectively, for sewer and surface representations. One of the observations made is that 2D models are computationally much more expensive than 1D models, resulting in a greater computation time. In general, Leandro et al. [28] say that the choice between using a 1D surface network model or a 2D surface system model depends on the physical behaviour mapped for the case study. This choice determines the reliability of the results and the computational time required to obtain them. Considering the simplicity of 1D models, it is easy to adjust their parameters and the result may be even better than that obtained by 2D models, where the uncertainties in the required data are greater and the difficulties for calibration are higher. In the presence of major overbank flows and storm drains failure, however, the flows occurring over the urbanised flood plain may determine the use of 2D models as the best option.
The results obtained by Leandro et al. [28] showed evidence that it is possible to set up an accurate 1D/1D model, depending on the problem characteristics. To accomplish this task it is necessary to define a detailed 1D surface network of pathways and possible ponds based on the physical reality (expressed, in this case, by a digital elevation model). Manholes should be linked to pathways, pathways to ponds and ponds to ponds. The connection between sewer and surface used a multiple linking element, depending on the characteristics of the flow. This conception is similar to that of the flow cell.
If the overland flow paths are well known, 1D surface models may be built in a sound and effective way, being an economical alternative to 2D models [29].
Simões et al. [2] presented a proposal that joined the 1D/1D and the 1D/2D approaches, in a hybrid configuration, aiming to take advantage of the benefits of each one and overcoming the drawbacks. A case study was developed for the Canbrook catchment in London (United Kingdom). Several tests were performed. The 1D/1D modelling was used to represent higher areas of the basin, while 1D/2D concept was applied in the areas where larger flooded areas were expected. Results showed that the hybrid model was almost as fast as the 1D/1D model, with results similar to the 1D/2D model written to the entire area.
Eleutério and Mosé [30] compared 1D, 1D/2D and 2D modelling for mapping riverine flooding in the town of Fislis, eastern France. Complementarily, 28 scenarios were also built varying the choices made by the modeller when representing the topography of the river-floodplain system. Considering the modelled results for the flooding areas and water depths, the authors highlighted that special attention should be paid both on modelling and modeller choices for representing physical situation. Although defining the model dimensions (1D, 1D/2D and 2D) appeared to be the most important choice, the fact is that physical modelling interpretation to accommodate different information scales and detailing levels might also significantly influence the precision of the modelled results.
Abily et al. [31] raised a pertinent and important discussion about possibilities, performances and limits of the use of standard modelling tools for high-resolution runoff simulations over an industrial site as an example of a dense environment with complex aboveground structures that affect (and alter) drainage paths. These are common characteristics also found on an urban environment. This work highlights (like previously cited works) that strict 2D models are being used to represent urban areas and high-resolution topography data are becoming available to support these models. However, these authors introduce new concerns arguing that the modelling of complex sites deserves special attention on topics such as rapid changes in flow regime, small water depths and high gradient properties, and vertical effects in runoff hydrodynamics introduced by some of the aboveground components. We can add to this discussion some other potential problems such as the wetting and drying of model elements, and discontinuities in 2D flood surface moving over the urban watershed (generated by physical obstacles). Walls or earthworks may produce independent surface flows connected by weirs or orifices, for example.
In particular, Abily et al. [31] tested two 2D models and mapped some of the difficulties mentioned. The authors also stress that field measurements would have been desirable to validate models. This final remark seems to be valid in most urban cases. Although 2D modelling is becoming more common and efforts are being made to benchmark models [32,33], the lack of observed data seems to be a common problem when comparing model simulations to real floods in spatially distributed real world events (e.g., [34,35]).

Urban Flood Cell Model-MODCEL
MODCEL is a quasi-2D model [36], as defined in the classical literature. It represents the two-dimensional characteristics of the watershed, but uses only 1D equations. MODCEL is able to describe natural or artificial watercourses and elements of the urban fabric (streets, squares, roofs, etc.), the flow in the subterranean storm drains, and the mutual "upper layer-lower layer" connections amongst such elements, including overflows. Because it integrates an upper and a lower layer of flows (superficial flows and storm drains), MODCEL can be seen (in a particular interpretation) as a quasi-3D model. However, by doing the flow representation through 1D equations written for pre-defined possible flow paths, the model preserves simplicity and spares computational time. On the other hand, this is also one of the weaknesses of this representation, once it fails to accurately model real 2D flow surfaces by not considering the cross influences of flow velocities in x,y-Cartesian axis. MODCEL also conjugates a simple hydrological distributed representation, performing rainfall-runoff transformations in each cell.
Topographic and hydraulic representations of the physical reality are two core elements in MODCEL. It is important to point out, however, that MODCEL does not need a Digital Terrain Model (DTM) in the strict sense-it is the responsibility of the modeller to interpret terrain features to provide synthetic geometrical information characterizing each individual cell and each connection between cells to the model. This way, it is important to note that, MODCEL cannot function by a blind, automatic application of geomatics tools and the objective description of boundary conditions. It rather obliges the modeller to investigate and understand how the real system works and how can each component be represented to simulate system functioning. At a first glance, this could Water 2017, 9, 445 7 of 28 be seen as another weakness (because it is neither automatic nor easy), but, through this process, indeed, inevitably, the particularities of the drainage system come to light and an actual, very useful understanding of the system behaviour is gained.

Background
In the 1960s, the first relevant mathematical model capable to describe two-dimensional flow patterns was proposed and implemented. This model was constructed for the Mekong delta river area by Societé Grenobloise d'Etudes et Applications Hydrauliques (SOGREAH), in a work for UNESCO [9,37]. This model represented the basin using storage cells, which were able to connect river reaches and floodplains to simulate flow patterns in the inundated Mekong delta areas, considering both natural floods and floods modified by a regulation dam, which would act in flood damping, flow regularization for navigation purposes and irrigation. The modelled area had about 50,000 km 2 and was inhabited by approximately 10 million people.
In 1990, a different kind of cell model, considering a mesh of rectangular cells, was developed for heavy floods simulation, including dam break [42]. In this latter case, zones near the dam were treated in a particular way, considering the shock wave and privileging inertia effects. Distant dam zones were considered to be characterized by great flood plains, terrain topography, land use, presence of dykes, among other factors, so the shock waves could be considered to have suffered significant dissipation. In this way, cells could represent distant flood plains.
The first version of a cell model developed in Brazil was built to represent Pantanal Matogrossense, a large marsh in west-centre region of Brazil [43]. Later, this model was adapted to an urban environment, gaining a series of new facilities, being presented by Mascarenhas and Miguez [8], in an initial version, that is in continuous improvement and which is presented in this paper on its current stage of development. This model was called MODCEL. Some applications of this model can be found in the literature (see, for instance, [44][45][46][47]).

Basic Concepts and Mathematical Modelling Structure
MODCEL is based in the basic principles of mass conservation, energy conservation and momentum conservation laws. The entire basin and its different elements are represented by cells, including the rivers, channels, storm drains, flood plains, hilly areas, urban areas, reservoirs, and other structures. Therefore, several urban structures and flow patterns can be simulated by the combination of a pre-defined set of cell types and cell links. The cells act as the storage elements of the model, representing terrain elevation and land occupation characteristics, while the links, on the other side, activate the flows between each pair of cells. Thus, cells and links are the computational basic elements of MODCEL.
The mass balance is applied to all cells. Thus, at each time step, the stored amount of water in one cell depends on the discharges exchanged between this cell and its neighbours, as well as on the rainfall contribution added through a rainfall run-off transformation. Two simple hydrologic models are available in MODCEL: the use of the Rational Method per cell, applying a runoff coefficient, evaluated according to the land use characteristics of each cell; or the use of a simple hydrological model to represent infiltration, vegetal interception and depression retentions, being the two latter parcels considered in a combined way as a limited reservoir. On the other hand, infiltration can occur as long as there is water accumulated over the surface of the modelled cell, even though the rainfall has finished. Abstraction (vegetal interception + retention) and infiltration are considered to occur simultaneously in the model, as part of the abstraction may occur in impervious surfaces.
The mesh of cells composes a hydrodynamic looped network, in a spatial representation that links surface flows, channel flows and storm drains, allowing discharges to occur in different directions on the modelled watershed, depending on the calculated water levels. All mathematical relations written are one-dimensional. Figure 1 illustrates the basic concepts of a cell. There are some important characteristics associated to every cell to adequately represent land surface properties: the total plain area, where the rainfall occurs; the storage area, where the mass balance is applied; and the land use and occupation characteristics, which is essential to estimate run-off generation. In urban cells, the different levels established by streets, sidewalks and buildings affect the storage capacity. Different patterns may be pre-defined in MODCEL, accounting for different characteristics of neighbourhoods in the urban tissue. These patterns influence storage availability (and consequent flooding levels) by establishing superficial areas associated with different terrain elevations. Figure 2 illustrates the different levels defined by the urban patterns inside an urbanised cell. The street level refers to the lowest area available for storage purposes. The sidewalk level refers to sidewalks themselves, but also to gardens, parks, and parking lots. The building level indicates the threshold related with initiating the flooding of houses and their contents. total plain area, where the rainfall occurs; the storage area, where the mass balance is applied; and the land use and occupation characteristics, which is essential to estimate run-off generation. In urban cells, the different levels established by streets, sidewalks and buildings affect the storage capacity. Different patterns may be pre-defined in MODCEL, accounting for different characteristics of neighbourhoods in the urban tissue. These patterns influence storage availability (and consequent flooding levels) by establishing superficial areas associated with different terrain elevations. Figure 2 illustrates the different levels defined by the urban patterns inside an urbanised cell. The street level refers to the lowest area available for storage purposes. The sidewalk level refers to sidewalks themselves, but also to gardens, parks, and parking lots. The building level indicates the threshold related with initiating the flooding of houses and their contents.   total plain area, where the rainfall occurs; the storage area, where the mass balance is applied; and the land use and occupation characteristics, which is essential to estimate run-off generation. In urban cells, the different levels established by streets, sidewalks and buildings affect the storage capacity. Different patterns may be pre-defined in MODCEL, accounting for different characteristics of neighbourhoods in the urban tissue. These patterns influence storage availability (and consequent flooding levels) by establishing superficial areas associated with different terrain elevations. Figure 2 illustrates the different levels defined by the urban patterns inside an urbanised cell. The street level refers to the lowest area available for storage purposes. The sidewalk level refers to sidewalks themselves, but also to gardens, parks, and parking lots. The building level indicates the threshold related with initiating the flooding of houses and their contents.   The set of pre-defined cell types used in MODCEL is listed below: • River/channel cells: This type of cell is used to model the main free open channel flows, in which the cross section is taken as a rectangular equivalent shape. • Storm drain cells: This type of cell represents enclosed sections in the drainage network, with a top limit to the flow depths. It allows superficial or surcharged flows. • Urban surface cells: These cells are used to represent urban flood plain surfaces, where run-off occurs, as well as flood flows in inundated areas. They are larger areas, when compared with river cells, and they act as storage areas linked to each other by a hydraulic link (a street functioning as a channel, for example). They also may represent slope areas, if a little storage area is assigned, leaving the most part of the total area just to receive rainfall contribution. • Natural surface cells: These cells are similar to the preceding ones, but they do not consider any kind of urbanisation pattern. The connection between two of these cells may be done by natural channels or thalwegs. When there are local elevations inside this kind of cell, the storage area is reduced to calculate adequately water levels inside the cell. • Reservoir cells: These cells are used to simulate a temporary pond or reservoir, represented by the relation of the terrain elevation with the surface area. Departing from this curve, it is possible to evaluate the stored volume variation from the water depth variation. The reservoir cell type may play the role of damping inflow discharges when representing flood control measures.
The water level variation in a cell i, at a time interval t, is given by the continuity equation applied for that cell as stated in Equation (1), in discrete terms.
where Q i,k is discharge between two neighbouring cells i and k; Z i is the water surface level at the centre of the cell i; A Si is the water surface area for the cell i; and P i is the discharge related to the rainfall over the cell. The discrete time interval (n + 1)·∆t, represented by the index t + 1, is taken as the calculation time, when the variables are unknown. On the other hand, at the discrete time n·∆t, represented by the index t, all variables are known, due to their previous calculation or because they were prescribed as initial conditions. The water surface area of the cell i, A t Si , is taken as a function of the known water level Z t i , implying that a first order approach, (∆A Si /A Si ) 1, is adopted. P i is a known value, since the rainfall is considered a known entry for all time intervals. Z t+1 i and Q t+1 i,k are unknowns and the numerical scheme is implicit. To solve this equation, the unknown discharges may be written as a function of the water levels in the cell i and all its k neighbour cells. Thus, to avoid the presence of a non-linear term, it is possible to develop Q t+1 i,k in Taylor series, taking only the first order terms, as shown in Equation (2).
Taking into account the previous discussion, Equation (1) may be re-written as Equation (3).
This way, analysing Equation (3), discharges are written in the time interval n·∆t, when all the variables are already known. The variables ∆Z t+1 i and ∆Z t+1 k are related respectively to the water levels Z t+1 i and Z t+1 k , written for each cell, and they are the only unknowns. They refer to the considered cell i and to its immediate adjacent k neighbours. Therefore, this system can be solved in a relatively simple way, once the mathematical relations for the discharges between cells are written for values of Z t i and Z t k . Cells are arranged in a topological scheme, which is numerically solved by a double sweep method [36,37].
The discharge links between cells can be expressed through known hydraulic laws. The possibility of introducing different mathematical relations is one of the features that inspire the use of the cell model to represent the urban floods diversity. The types of links considered in the current version of the model are briefly presented and discussed in the following.
Types of Links: • River/Channel link: This link is related to river and channel free flows corresponding to the De Saint Venant dynamic equation. Equation (4) results from the consideration of rectangular cross sections.
where A i,k is the wetted cross-section area between cells i and k; S f is the energy line slope; and the other variables were previously defined.
The parameter A i,k , and any other that appears in the flow section between cells i and k, may be evaluated through a weighting procedure departing from the values of water levels Z i and Z k .
Note that, when working with natural river stretches, the hypothesis of a rectangular cross section may be strong and should be used carefully. However, when modelling channelized urban rivers and considering artificial drainage systems, this is a reasonable approach. Besides that, the cell concept allows the modeller to compose a cross section using multiple cells with different terrain levels, representing the main channel, the secondary channels and the flood plains.
Multiplying Equation (5) Remembering that the discharges relations were written in the time interval n· ∆t, the expression to evaluate this discharge explicitly provides the entrance required in the modified continuity Equation (3). Writing Equation (5) in discrete terms leads to Equation (6), considering cell i in an upstream position: S f , by its turn, may be approximated by Expression (7): where: n is Manning's roughness coefficient; R i,k is the hydraulic radius of the flow cross-section between cells i and k.
Combining Equations (6) and (7), the term Q t i,k is the focus of interest. To explicate this term and obtain a direct solution for this equation, the quadratic terms were "factored", using a numerical simplification as it can be seen in Expression (8).
From this discussion, Equation (9) for the river flow link is obtained and may be used in the mass conservation balance: • Surface flow link: This link corresponds to the free surface flow without inertia terms, as presented in Zanobetti et al. [37]. MODCEL uses this link frequently to represent flow between surface cells (natural and/or urban). There are two possibilities to use this link. The first one, which refers to the original representation, considers the main drainage line entering in a closed cross section. In this case, MODCEL considers the closed storm drain as cells linked between each other. They also communicate with superficial cells, by manholes to a street above them. In this situation, flow over the street, represented by an urban surface cell, occurs with free surface and the water level associated to this flow can be considered equal to the piezometric line level for the drowned storm drain flow. In this case, cross sections are also rectangular. It is important to note that mass balance is calculated accordingly with real geometry, and the model keeps track of the cross section water levels and possible surcharge.
In the second possibility, which is used to represent minor drainage network, the manholes are represented by cells, but the storm drains are simple links, without the need to introduce a formal cell. Cross sections are taken as circular ones.
Thus, departing from Bernoulli Equation (10), Equation (11) is developed, where cell i is in the upstream position. Figure 3 shows the storm drain link representation. It is important to stress that, with these considerations, depending on the water levels involved in the calculation, flow might be forced backwards (what can really happen in physical reality). In MODCEL, junctions of river reaches are treated as special cells with "Y" shape, where mass balance equates inflows and outflows and, therefore, this kind of confluence does not need a special link.

•
Inlet link: This link promotes the interface between the street gutters of the surface cells and the underground storm drain cells. When not drowned, it acts as an equivalent weir conveying flow from streets to storm drains. This weir has the length of the perimeter of a single inlet multiplied by the number of inlets along the street modelled by the considered cell. When drowned, this link considers flow occurring through a certain number of orifices associated to the inlet grates in the street. It is useful to simulate the damping effect of a reservoir, in the design condition, and to verify reservoir operation in more severe conditions (those in which the weir can start to be used).  In MODCEL, junctions of river reaches are treated as special cells with "Y" shape, where mass balance equates inflows and outflows and, therefore, this kind of confluence does not need a special link.

•
Inlet link: This link promotes the interface between the street gutters of the surface cells and the underground storm drain cells. When not drowned, it acts as an equivalent weir conveying flow from streets to storm drains. This weir has the length of the perimeter of a single inlet multiplied by the number of inlets along the street modelled by the considered cell. When drowned, this link considers flow occurring through a certain number of orifices associated to the inlet grates in the street. • Broad crested weir link: This link represents the flow over broad-crested weirs. It is used, mainly, to represent the flow between a river and its riverine areas. The classic formula of flow over broad-crested weirs is used here. Flow over a weir may be free or drowned, depending on water levels of the cells connected by this link. • Orifice link: This link represents the classic formula for flow through orifices. • Reservoir link: This link combines orifices (at different possible heights), as the outlet discharge of a reservoir, with a weir, which can enter or not in charge, depending on reservoir operation. It is useful to simulate the damping effect of a reservoir, in the design condition, and to verify reservoir operation in more severe conditions (those in which the weir can start to be used).

This section intends to show MODCEL performance in a benchmark test proposed in a research commissioned by the British Environment Agency in the context of the Flood and Coastal Erosion
Risk Management Research and Development Programme. The resulting report of this research [33] describes "the results from a benchmarking exercise assessing the latest generation of 2D hydraulic modelling tools for a variety of purposes in Flood and Coastal Risk Management to support Environment Agency decision making". This benchmarking exercise involves 10 test cases and one of the objectives of this research is to provide a data set against which a model can be evaluated by its developer.
Test 8, in particular, was designed to benchmark models for urban flood modelling and it is divided in two parts, 8A and 8B.
Test 8A assumes that the flood arises from two sources: a uniformly distributed rainfall event (peaking at 400 mm/h over a time base of 3 min), applied to the modelled area; and a point discharge source occurring over a time base of approximately 15 min, reaching a peak at 5 m 3 /s, 35 min after the rainfall event. This point source occurs at the location represented by the "INFLOW" indication in Figure 4. In Test 8B, a culverted watercourse of circular section is assumed to run through the site, with a single manhole at the location indicated as "MANHOLE" in Figure 4. An inflow boundary condition is applied at the upstream end of the pipe, forcing a surcharge at the manhole, aiming at assessing the model capability to simulate shallow inundation resulting from a surcharged underground pipe.
The modelled area is a rectangle of approximately 0.4 km by 0.96 km as it is shown in Figure 4. Model grid resolution is 2 m, resulting in~97,000 nodes (0.388 km 2 ). The test considers initial dry bed, Manning's coefficient of 0.02 for roads and pavements and 0.05 elsewhere Water 2017, 9,445 13 of 28

This section intends to show MODCEL performance in a benchmark test proposed in a research commissioned by the British Environment Agency in the context of the Flood and Coastal Erosion
Risk Management Research and Development Programme. The resulting report of this research [33] describes "the results from a benchmarking exercise assessing the latest generation of 2D hydraulic modelling tools for a variety of purposes in Flood and Coastal Risk Management to support Environment Agency decision making". This benchmarking exercise involves 10 test cases and one of the objectives of this research is to provide a data set against which a model can be evaluated by its developer.
Test 8, in particular, was designed to benchmark models for urban flood modelling and it is divided in two parts, 8A and 8B.
Test 8A assumes that the flood arises from two sources: a uniformly distributed rainfall event (peaking at 400 mm/h over a time base of 3 min), applied to the modelled area; and a point discharge source occurring over a time base of approximately 15 min, reaching a peak at 5 m 3 /s, 35 min after the rainfall event. This point source occurs at the location represented by the "INFLOW" indication in Figure 4. In Test 8B, a culverted watercourse of circular section is assumed to run through the site, with a single manhole at the location indicated as "MANHOLE" in Figure 4. An inflow boundary condition is applied at the upstream end of the pipe, forcing a surcharge at the manhole, aiming at assessing the model capability to simulate shallow inundation resulting from a surcharged underground pipe.
The modelled area is a rectangle of approximately 0.4 km by 0.96 km as it is shown in Figure 4. Model grid resolution is 2 m, resulting in ~97,000 nodes (0.388 km 2 ). The test considers initial dry bed, Manning's coefficient of 0.02 for roads and pavements and 0.05 elsewhere  Figure 5 shows the test area divided into cells, resulting in a mesh of 163 cells (contrasting with the expected 97,000 simulation nodes). Roads were individualised to act as channels, while blocks (with buildings and open green areas) were mainly taken as storage areas, although they are also allowed to represent surface flows. Figures 6 and 7 show the results obtained with MODCEL, respectively, for Test 8A and Test 8B, compared with the published results [33]. Time of simulation using MODCEL to obtain these results was 3 min and 10 s. The simulation time step gave answers at each 1 s.  Figure 5 shows the test area divided into cells, resulting in a mesh of 163 cells (contrasting with the expected 97,000 simulation nodes). Roads were individualised to act as channels, while blocks (with buildings and open green areas) were mainly taken as storage areas, although they are also allowed to represent surface flows. Figures 6 and 7 show the results obtained with MODCEL, respectively, for Test 8A and Test 8B, compared with the published results [33]. Time of simulation using MODCEL to obtain these results was 3 min and 10 s. The simulation time step gave answers at each 1 s.     Observing the results obtained with MODCEL, when compared with the set of 2D models that participated in the test and that had their results published in Néelz and Pender [33], it is possible to note that there are no significant differences. In Test 8A, the major differences refer to the velocity measures-and this is, probably, the major weakness of MODCEL, due to the hydrodynamic 1D relations, that are not able to capture transversal interferences. Test 8B showed equivalent results among the different models and MODCEL, both for discharge and water levels observations.

MODCEL and Minor Drainage Representation at La Riereta Watershed, Saint Boi de Lobregat/Spain
This test was introduced to compare MODCEL with a well-known software for urban storm drains simulation-the Storm Water Management Model (SWMM), developed by the United States Environmental Protection Agency (US-EPA), was used as reference. La Riereta watershed is located in the old town of Sant Boi de Llobregat, near Barcelona. This urban basin has a superficial area of nearly 180,000 m 2 , with a high level of impervious surfaces and medium slope values. Runoff generated by rooftops is dropped directly to the streets where a large set of grate inlets are located on sidewalks to assure the capture of the superficial flows. The sewer of "La Riereta" is a combined system, consisting mostly of circular conduits with variable diameters. The majority of the pipes are made of concrete. Figure 8 shows a map of this watershed. According to topographic characteristics, This Figure also shows 17 sub-catchments in a topological scheme defined for SWMM modelling purposes. The same division was assumed in MODCEL, to test its capabilities in the same conditions and without increasing the level of details. In general terms, the manholes were taken as cells, and storm drain links ware written between each pair of manholes, connecting them. Then, the 17 subcatchments were defined as urban surface cells and they were connected with manholes by inlet links, and they were also connected with each other by surface flow links. Observing the results obtained with MODCEL, when compared with the set of 2D models that participated in the test and that had their results published in Néelz and Pender [33], it is possible to note that there are no significant differences. In Test 8A, the major differences refer to the velocity measures-and this is, probably, the major weakness of MODCEL, due to the hydrodynamic 1D relations, that are not able to capture transversal interferences. Test 8B showed equivalent results among the different models and MODCEL, both for discharge and water levels observations.

MODCEL and Minor Drainage Representation at La Riereta Watershed, Saint Boi de Lobregat/Spain
This test was introduced to compare MODCEL with a well-known software for urban storm drains simulation-the Storm Water Management Model (SWMM), developed by the United States Environmental Protection Agency (US-EPA), was used as reference. La Riereta watershed is located in the old town of Sant Boi de Llobregat, near Barcelona. This urban basin has a superficial area of nearly 180,000 m 2 , with a high level of impervious surfaces and medium slope values. Runoff generated by rooftops is dropped directly to the streets where a large set of grate inlets are located on sidewalks to assure the capture of the superficial flows. The sewer of "La Riereta" is a combined system, consisting mostly of circular conduits with variable diameters. The majority of the pipes are made of concrete. Figure 8 shows a map of this watershed. According to topographic characteristics, This Figure also shows 17 sub-catchments in a topological scheme defined for SWMM modelling purposes. The same division was assumed in MODCEL, to test its capabilities in the same conditions and without increasing the level of details. In general terms, the manholes were taken as cells, and storm drain links ware written between each pair of manholes, connecting them. Then, the 17 sub-catchments were defined as urban surface cells and they were connected with manholes by inlet links, and they were also connected with each other by surface flow links.      It is possible to note that both models performed reasonably well, representing adequately flood peaks, peak phase and hydrograph form. Simulated volumes, however, are a little bit greater than the observed ones.

Case Study of Noale-MODCEL and the Simulation of Complex Flood Plains
The main purpose of bringing here this last case is to illustrate model capabilities when representing flood plains with a complex topology, in which artificial structures heavily changes the natural patterns and creates several local disconnected flows.
Particularly here, the proposed example refers to Noale, a city in the Venetian continental plains, in Italy, where urban flood occurs frequently and this phenomenon is driven by anthropic actions that modified the watershed with artificial hydraulic structures. This study was conducted in the context of The European Project SERELAREFA-SEmillas REd LAtina Recuperación Ecosistemas Fluviales y Acuáticos (Seeds of a Latin American Network for the Restoration of Fluvial and Aquatic Ecosystems), funded by the UE programme FP7 IRSES-PEOPLE 2009. The main aim of this project was to improve the way watercourses are managed by achieving benefits for both the environment and socio-economic activities. SERELAREFA fostered the adoption of River Restoration concept and developed several case studies that were supported by mathematical modelling. The Project started in September 2010 and lasted four years. Partners from Italy, Spain, Mexico, Chile and Brazil were involved in several researches regarding the use of river restoration approach to restore riparian areas, recover hydrological and morphological processes, and control floods, improving river ecosystems in general, as well as city environment, in the urban cases. The general coordination was issued to the Centro Italiano per la Riqualificazione Fluviale-CIRF (Italian Centre for River Restoration). It is possible to note that both models performed reasonably well, representing adequately flood peaks, peak phase and hydrograph form. Simulated volumes, however, are a little bit greater than the observed ones.

Case Study of Noale-MODCEL and the Simulation of Complex Flood Plains
The main purpose of bringing here this last case is to illustrate model capabilities when representing flood plains with a complex topology, in which artificial structures heavily changes the natural patterns and creates several local disconnected flows.
Particularly here, the proposed example refers to Noale, a city in the Venetian continental plains, in Italy, where urban flood occurs frequently and this phenomenon is driven by anthropic actions that modified the watershed with artificial hydraulic structures. This study was conducted in the context of The European Project SERELAREFA-SEmillas REd LAtina Recuperación Ecosistemas Fluviales y Acuáticos (Seeds of a Latin American Network for the Restoration of Fluvial and Aquatic Ecosystems), funded by the UE programme FP7 IRSES-PEOPLE 2009. The main aim of this project was to improve the way watercourses are managed by achieving benefits for both the environment and socio-economic activities. SERELAREFA fostered the adoption of River Restoration concept and developed several case studies that were supported by mathematical modelling. The Project started in September 2010 and lasted four years. Partners from Italy, Spain, Mexico, Chile and Brazil were involved in several researches regarding the use of river restoration approach to restore riparian areas, recover hydrological and morphological processes, and control floods, improving river ecosystems in general, as well as city environment, in the urban cases. The general coordination was issued to the Centro Italiano per la Riqualificazione Fluviale-CIRF (Italian Centre for River Restoration).

Rationale of Case Study
The lowland area of Mestre/Venice is naturally subject to harsh flooding as a significant part of it lies in areas of low absolute elevations, and even below average sea level. Consequently, it is artificially drained and it has been heavily modified over time. Dramatic events in 2006 and 2007 led to a mandate to a special Commissioner who issued a specific law that introduced the "hydraulic invariance" concept. The idea is that if a new building/infrastructure would increase the peak discharge of a reference flood event (return period of 50 years)-owing to the loss of natural storage volume and infiltration capacity-then building permission could be granted only if provided that suitable compensatory measures are put in place to maintain the original peak value under control. This concept has been fully developed by the Piano Territoriale di Coordinamento Provinciale (Territorial Coordination Provincial Plan) [48], approved by Venetian province administration (December 2010), which stated the need to elaborate specific Water Plans.
More recently, guidelines for the ecological restoration of the extensive network of natural and artificial irrigation and drainage canals [49] were issued, and several concepts of sustainable water management and river restoration were included, much in line with SERELAREFA philosophy.
Within the framework of SERELAREFA project, a first study trip happened in September 2010, during which the project partners could understand the problems faced in that region. After this visit, the group of the Federal University of Rio de Janeiro took part in a relevant Italian conference (Acqua e Città held in Venezia, 2011) and then in a specific workshop in Mestre. This workshop, hosted by Consiglio di Bacino Laguna di Venezia (Venetian Lagoon Watershed Council), counted with the participation of several key actors, amongst which the technical staff of hydro-geological defence of the Venetian Province. In this opportunity, the modelling framework of MODCEL was presented and this event opened a door for setting up an Italian case study around the Venetian area, conceived as an academic opportunity of applied research and know-how exchange.

Noale
Noale is a city that has its origins in ancient times. From the middle of the 1800s until today, great changes happened and its population was approximately multiplied by four. The centre of the city is flooded by Marzenego River, where one of its branches overflows the dike of the historic city castle's moat. The discharges that surpass the dyke go superficially to the south, towards Roviego River.
Urban grow and consequent land use modifications generated an increase in the hydraulic risks reaching important central areas and affecting a greater exposed population. Today, flood risk management is of prime concern for the protection of properties and economic assets. The first step in this path refers to understanding and mapping the interrelations between the city and the river. Therefore, this first step is the main objective of this case study, which also intends to illustrate MODCEL capability to map the current flooding situation in this complex flood plain.

General Hydraulic Description
The Marzenego River passes through Noale's centre and defines the main river course of the city. However, this watershed is not as usual and very few natural characteristic (if any) are preserved. The Marzenego River runs from northwest to southeast. In the north region, Dragonziolo River is its main tributary to the left bank. South to Marzenego runs Roviego River. Figure 11 shows a general plan view of this region and Figure 12 highlights some points located and numbered from 1 to 4 in Figure 11.
All Rivers running west to east in this region are contained by dikes on the right bank, so that the riverine areas that should be connected to a river by its right bank are artificially drained to the left riverbank of the next river to the south. The dyke in the centre of the city, referred to the castle moat, has a low height (approximately 50 cm) and frequently overflows.
River running north to south may experience "river crossings", both by siphons or by elevated channels (similar to aqueducts over earth works), as shown in Figure 13. The main urban storm drain passes below Marzenego River, reaching a channel that outflows in the Roviego River.
River running north to south may experience "river crossings", both by siphons or by elevated channels (similar to aqueducts over earth works), as shown in Figure 13. The main urban storm drain passes below Marzenego River, reaching a channel that outflows in the Roviego River. Marzenego River also presents a series of watermills in the main river course, which limit discharges flowing through these reaches. There are also some pumping stations and gates that artificially control (and changes) the flows in this complex network. Marzenego River also presents a series of watermills in the main river course, which limit discharges flowing through these reaches. There are also some pumping stations and gates that artificially control (and changes) the flows in this complex network. Water 2017, 9, 445 20 of 28

Modelling Noale's Case
Marzenego River was modelled from Molino Cosma to Robegano City. The upstream reach of Marzenego River, until reaching the city of Castelfranco, was transformed into a boundary condition. The upper reach of Dragonziolo River and Scolo Fossalta were also taken as boundary conditions. Roviego River was modelled as an outfall. The whole area between Marzenego and Roviego may drain to this river and all the drainage network-conducting waters into this river may discharge there, but the Roviego hydrodynamics was not represented. Figure 13. Example of a complex hydraulic site in the watershed, in a reach upstream the city centre: (a) Marzenego River running between dykes-on the right bank stands a deviation structure working for flood alleviation; (b) Deviation channel-looking to the deviation structure from downstream to upstream-note that on the left bank it is possible to see the outfall of a channel that crosses below Marzenego River; (c) Schematic view of the site, showing the river, its dykes, the deviation weir, the deviation channel and a former tributary passing below Marzenego River, also going to the deviation channel and (d) physical interpretation of the system working and how it is modelled.

Modelling Noale's Case
Marzenego River was modelled from Molino Cosma to Robegano City. The upstream reach of Marzenego River, until reaching the city of Castelfranco, was transformed into a boundary condition. The upper reach of Dragonziolo River and Scolo Fossalta were also taken as boundary conditions. Roviego River was modelled as an outfall. The whole area between Marzenego and Roviego may drain to this river and all the drainage network-conducting waters into this river may discharge there, but the Roviego hydrodynamics was not represented.
The topographic information was obtained from maps in the scale of 1:5000, from Segreteria Regionale All'Ambiente e Territorio (Territory and Environment Regional Secretariat). The modelled area covered an area of approximately 42 km 2 and 20 km of Marzenego River, with a medium slope of 0.0006. The time of concentration for the whole watershed was estimated as 12 h, using the Kirpich formula.
The model implemented in this study has 63 channel cells representing Marzenego River and other 54 cells representing the other water courses. There are also 219 cells representing the flood plains and the city of Noale, while 11 cells represent the storm drains, totalising 347 cells. It is important to stress that this model includes only the main storm drains. Figure 14 shows the watershed, including the boundary conditions and the cell division, while Figure 15 shows the cell division over an aerial view.
A compilation of available information on Venetian floods was made, resulting in the map shown in Figure 16. The flooding event of 2006 was used to calibrate the model, using the rainfall measured. The envelope of important flooding events, corresponding to their superimposed effects, was used to validate the model. In this latter case, a rainfall of 50 years of return period was introduced uniformly in the model intending to reproduce this larger flood map.
Two rainfall gauges were taken as representative of the modelled area: Istrana (14 km far from Noale's city centre) and Stra (14 km far from Noale's city centre) [50].
Considering a 24 h event, but introducing a time variation component, using the information of the rainfall values occurred in 1, 3, 6 and 12 h, a design rainfall of 50 years of return period [51] was built, using the alternate blocks method from the US Soil Conservation Service. The resulting rainfall was defined for time intervals of 0.5 h, intending to represent a critical situation to the whole basin and also to the sub-basins with different times of concentration. That is, if we are working with a storm drain, in a street in the middle of the city, the higher rainfall intensity related with a single time interval (of 30 min) will be critical in the local scale, while the composed rainfall intensity for the entire modelled basin scale (12 h) will also be critical. Figure 17 shows the design rainfall built for the return period of 50 years.
The topographic information was obtained from maps in the scale of 1:5000, from Segreteria Regionale All'Ambiente e Territorio (Territory and Environment Regional Secretariat). The modelled area covered an area of approximately 42 km 2 and 20 km of Marzenego River, with a medium slope of 0.0006. The time of concentration for the whole watershed was estimated as 12 h, using the Kirpich formula.
The model implemented in this study has 63 channel cells representing Marzenego River and other 54 cells representing the other water courses. There are also 219 cells representing the flood plains and the city of Noale, while 11 cells represent the storm drains, totalising 347 cells. It is important to stress that this model includes only the main storm drains. Figure 14 shows the watershed, including the boundary conditions and the cell division, while Figure 15 shows the cell division over an aerial view.
A compilation of available information on Venetian floods was made, resulting in the map shown in Figure 16. The flooding event of 2006 was used to calibrate the model, using the rainfall measured. The envelope of important flooding events, corresponding to their superimposed effects, was used to validate the model. In this latter case, a rainfall of 50 years of return period was introduced uniformly in the model intending to reproduce this larger flood map.
Two rainfall gauges were taken as representative of the modelled area: Istrana (14 km far from Noale's city centre) and Stra (14 km far from Noale's city centre) [50].
Considering a 24 h event, but introducing a time variation component, using the information of the rainfall values occurred in 1, 3, 6 and 12 h, a design rainfall of 50 years of return period [51] was built, using the alternate blocks method from the US Soil Conservation Service. The resulting rainfall was defined for time intervals of 0.5 h, intending to represent a critical situation to the whole basin and also to the sub-basins with different times of concentration. That is, if we are working with a storm drain, in a street in the middle of the city, the higher rainfall intensity related with a single time interval (of 30 min) will be critical in the local scale, while the composed rainfall intensity for the entire modelled basin scale (12 h) will also be critical. Figure 17 shows the design rainfall built for the return period of 50 years.          Figure 19 shows the inundation map for a rainfall event with 50-year return period. It is possible to see that the modelled results are adherent in both cases.    Figure 19 shows the inundation map for a rainfall event with 50-year return period. It is possible to see that the modelled results are adherent in both cases.   Figure 19 shows the inundation map for a rainfall event with 50-year return period. It is possible to see that the modelled results are adherent in both cases.   When considering a rainfall with 50-year return period distributed over the entire basin, the basic intention was to reproduce (for this representative rainfall intensity) the enveloped results of all previous important observed floods. This aim was reached and the results shown in Figure 19 are quite similar to those shown in Figure 16. An important thing to be added is that this result was submitted to the appreciation of the Venetian Lagoon Watershed Council and the local staff considered that the general flooding representation and paths of flows were quite accurate.

Concluding Remarks
Mathematical models are useful tools for managing drainage systems, especially when flooding problems produces large inundation areas. In these cases, the drainage system tends to interact with urban landscape and urban structures may act hydraulically composing a new flow pattern distinct from that predicted in the original project. Most times, the superficial processes may play a major role in the flood problem, being affected by urban structures and local obstacles in a complex arrangement. This is where the main potential of MODCEL arises, focusing on the interactions between surface flows and the drainage net flows, also allowing the introduction of typical urban structures in this representation. This is an important feature to consider in this kind of situation, providing the possibility to explore different integrated flood control alternatives by placing different types of interventions over the basin landscape, even outside the drainage network.
In this paper, MODCEL was benchmarked in a 2D test proposed by the British Environmental Agency, showing water level results compatible with 2D Models, but using only 163 cells, against 97,000 elements. Then MODCEL was compared with a storm drain model to show its capacity to represent the 1D branched network. Finally, MODCEL was applied to a complex watershed in the lowland area of Mestre/Venice, where the city of Noale lies in the Marzenego River basin. The results were compared with known flooding events. However, the case of Noale could only allow a spatial calibration (flooded area) and not its temporal evolution. The simulation of this watershed, in particular, was important to show MODCEL capability to deal with a mixed and complex setup, where a plain territory interacts with artificial structures and non-natural flow paths, including the When considering a rainfall with 50-year return period distributed over the entire basin, the basic intention was to reproduce (for this representative rainfall intensity) the enveloped results of all previous important observed floods. This aim was reached and the results shown in Figure 19 are quite similar to those shown in Figure 16. An important thing to be added is that this result was submitted to the appreciation of the Venetian Lagoon Watershed Council and the local staff considered that the general flooding representation and paths of flows were quite accurate.

Concluding Remarks
Mathematical models are useful tools for managing drainage systems, especially when flooding problems produces large inundation areas. In these cases, the drainage system tends to interact with urban landscape and urban structures may act hydraulically composing a new flow pattern distinct from that predicted in the original project. Most times, the superficial processes may play a major role in the flood problem, being affected by urban structures and local obstacles in a complex arrangement. This is where the main potential of MODCEL arises, focusing on the interactions between surface flows and the drainage net flows, also allowing the introduction of typical urban structures in this representation. This is an important feature to consider in this kind of situation, providing the possibility to explore different integrated flood control alternatives by placing different types of interventions over the basin landscape, even outside the drainage network.
In this paper, MODCEL was benchmarked in a 2D test proposed by the British Environmental Agency, showing water level results compatible with 2D Models, but using only 163 cells, against 97,000 elements. Then MODCEL was compared with a storm drain model to show its capacity to represent the 1D branched network. Finally, MODCEL was applied to a complex watershed in the lowland area of Mestre/Venice, where the city of Noale lies in the Marzenego River basin. The results were compared with known flooding events. However, the case of Noale could only allow a spatial calibration (flooded area) and not its temporal evolution. The simulation of this watershed, in particular, was important to show MODCEL capability to deal with a mixed and complex setup, where a plain territory interacts with artificial structures and non-natural flow paths, including the operation of gates, pumping stations and channel crossings (where overflows cause discharges transpositions from one watershed to another), MODCEL showed to be able to reproduce this great variety of hydraulic patterns in the urban landscape, as well as the set of hydraulic structures working jointly with natural and artificial drainage network. Although using a relative simple squasi-2D scheme, the physical interpretation phase that is needed to use this type of model allowed obtaining adequate and representative results, with low computational costs and less information need. The main weakness of this proposal refers to great flooding areas, where significant high flood levels occur (producing a continuous and extensive flooding surface). In this case, real 2D models are expected to produce better results.