A Model-Based Tool for Assessing the Impact of Land Use Change Scenarios on Flood Risk in Small-Scale River Systems—Part 2: Scenario-Based Flood Characteristics for the Planned State of Land Use

Land use changes can significantly influence the water balance, and thus especially the development of flood-triggering runoff peaks. Hence, it is advisable to assess possible changes already at the level of municipal planning. Since many different actors are usually involved in spatial planning, it is useful to provide a shared platform where stakeholders can access the same information to analyze and evaluate flood hazards. Therefore, a GIS routine for the prediction of soil sealing induced runoff peaks and resulting potential flooding in the watercourse was developed, which is embedded in a GIS based decision support system (GIS-DSS). The so-called storm water routine (SWR) is founded on preprocessed flood characteristics, calculated by means of hydrological/hydraulic models (described in part 1). The potential impact of land use change is assessed purely in GIS as flow difference which is routed through the river system. To validate this simplified method, a process model was set up with an exemplary land use change and its results were compared with the GIS-based results. For 16 of the 18 rainfall scenarios tested, the SWR provided very good to good agreement with the detailed model. For short and highly dynamic rain events the SWR approach is less reliable. Several supplements like the integration of LID are conceivable.


Introduction
Since the beginning of human settlement history, humans have preferred to settle near rivers and streams. Small settlements became large settlements and these became today's cities. As settlements progress, areas formerly close to nature are rezoned as building land and cities are redensified. This inevitably leads to increased soil sealing. Soil sealing strongly modifies the water cycle by reducing infiltration, evapotranspiration, and groundwater recharge, thus intensifying (rapid) surface runoff. As a result, the risk of flooding is increased [1,2] especially where no retention areas exist. For a region in China it has been shown that over the last 20 years urbanisation has led to an increase in the runoff coefficient of 13.4%, resulting in an increase in the maximum flood discharge of 12.9% on average [3]. In [4], a methodology is described which relates urban growth studies based on the analyzation of satellite data to distributed hydrological runoff modeling using the soil conservation service (SCS) method. The results of the corresponding case study indicate that annual runoff depth had increased by 8.10 mm between 1989 and 1997 due to urban sprawl.
The social and economic 'costs' of flooding are expected to rise in the future, not only due to land use changes, but also due to climate change [5][6][7]. In order to minimize flood risk by a proper communal planning, the impact of land use changes on runoff and river surcharge must be considered at an early stage. In this context, the cooperation and exchange of information between all relevant authorities, organisations, associations, and

•
The DSS has to address the actual goals of the decision-makers and stakeholder wishes • The DSS should provide a user-friendly interface and good visualization capabilities for a real participatory use by the stakeholders and decision-makers • The stakeholder and decision-makers should have a clear understanding of the model concept and should ideally be able to edit it by themselves e.g., for scenario analysis purposes There is a wide variety of approaches to address these, depending on the regarded topics, end-users and scientific background of the model developers. They reach from statistical and/or conceptual models processing meta information [16] to high resolute and well-considered 2D/3D flood risk mapping [17][18][19]. Partly, they are collaboratively developed with the end-users [20] or even use virtual reality environment to increase information depth and awareness. In tendency, the participatory character of the model decreases with increasing model accuracy. The concept of [16] is highly participative but does not provide physically exact information. Particularly, 3D maps are mainly prepared for citizens.
Our concept is ranging between these two edges: it is highly cooperative and based on deterministic physical modeling. More precisely, a tool was developed that allows an interactive calculation of the runoff from a particular planned land use and the subsequent propagation of the flow in the river system. The so called 'storm water routine' (SWR) was developed, as part of a larger GIS-DSS for regional and municipal planning activities [21]. The development is primarily based on open-source tools to enable barrier-free use by relevant stakeholders. It combines: • a physical process model of the catchment hydrology and river hydraulics, set up for the current state of land use [22] and • a GIS-routine calculating the additional runoff for land use change scenarios and its routing through the stream system.
The objectives of this tool have been defined in a participatory process with regional water authorities, water associations, and wastewater operators. They identified the lack of information on land use change induced changed runoff and missing information on specific hydraulic capacity at the point of discharge as major deficiencies in the approval of new storm water discharge permits. In distinct cases with expected larger impacts, individual studies are commissioned to consultancies, but a consistent and continuous assessment for a total river basin is missing. While most of the partners do not have their own modeling expertise, the application of GIS is meanwhile daily practice. GIS also provides the required visualization capabilities and interface for editing scenarios in the area of interest. Accordingly, it was agreed to develop a tool to holistically assess the hydrologic and hydraulic impacts of land use changes for the total river system. In contrast to commissioned studies, the tool should be operable and editable by the decision-makers and stakeholders themselves by using the well-known GIS environment. Since the tool is based on the pre-processed results of hydrologic/hydraulic models, but is not directly linked to the models, the end user does not necessarily need to have modeling experience.
While part 1 of this paper has described the set-up and parametrization of the process model and the definition of rainfall scenarios, this article focuses on the GIS-based calculation of additional runoff and the routing. The reliability of the approach has been reviewed by a comparative calculation of river flow with the detailed process model and the simplified GIS-based approach, revealing its strengths and weaknesses.

Study Area
The study area, approximately 530 km 2 in size, is located in the northeastern German lowlands and includes the Hanseatic City of Rostock and its 18 surrounding municipalities ( Figure 1). The latter are part of the Rostock rural district. While the core of the city is dominated by uses such as housing, industry, and commerce as well as transport areas, the surrounding rural district is mainly characterised by arable land and grassland. The study area overlaps 19 largely independent surface catchments, with the headwaters of the streams mostly in the rural district and the mouths in the urban area of Rostock or the Baltic Sea. of new storm water discharge permits. In distinct cases with expected larger impacts, individual studies are commissioned to consultancies, but a consistent and continuous assessment for a total river basin is missing. While most of the partners do not have their own modeling expertise, the application of GIS is meanwhile daily practice. GIS also provides the required visualization capabilities and interface for editing scenarios in the area of interest. Accordingly, it was agreed to develop a tool to holistically assess the hydrologic and hydraulic impacts of land use changes for the total river system. In contrast to commissioned studies, the tool should be operable and editable by the decision-makers and stakeholders themselves by using the well-known GIS environment. Since the tool is based on the pre-processed results of hydrologic/hydraulic models, but is not directly linked to the models, the end user does not necessarily need to have modeling experience. While part 1 of this paper has described the set-up and parametrization of the process model and the definition of rainfall scenarios, this article focuses on the GIS-based calculation of additional runoff and the routing. The reliability of the approach has been reviewed by a comparative calculation of river flow with the detailed process model and the simplified GIS-based approach, revealing its strengths and weaknesses.

Study Area
The study area, approximately 530 km 2 in size, is located in the northeastern German lowlands and includes the Hanseatic City of Rostock and its 18 surrounding municipalities ( Figure 1). The latter are part of the Rostock rural district. While the core of the city is dominated by uses such as housing, industry, and commerce as well as transport areas, the surrounding rural district is mainly characterised by arable land and grassland. The study area overlaps 19 largely independent surface catchments, with the headwaters of the streams mostly in the rural district and the mouths in the urban area of Rostock or the Baltic Sea.  The history of the city of Rostock goes back to the 13th century. Since then, the city has grown not only in its territorial area but also in its population, reaching nearly 210.000 inhabitants today [23]. For Rostock, an increase of nearly 15.000 inhabitants is expected by 2035 [24]. In contrast, projections indicate that the rural district will experience a population decline of 1.8% by 2030 [25]. However, as in many developed countries, population and sealing do not develop in parallel, but the demand for building ground is increasing, even in areas with shrinking population. Main drivers are the demand for residential units, the development of logistic centers and industrial parks.

Basic Data
In the course of the study, spatial data from different sources were used and integrated into a homogeneous and consistent data set. The data relevant for the application of the storm water routine are presented in the following.

Land Use Map
The Hanseatic City of Rostock provides a spatially high-resolution geodata set that was created on the basis of aerial photographs and supplemented by existing geodata. Since such a dataset did not exist for the rural district, it was generated from the data of the Official Real Estate Cadastre Information System (German abbreviation: ALKIS) in combination with satellite data. The result is a polygon shape file that depicts the land use in its actual state on the basis of 48 classes. In the GIS-DSS, this land use map represents the reference state for all future changes of land use. In addition, different levels of protection (risk classes) were defined for each land use class according to their vulnerability (see [22]). The level of protection is expressed by the return period which describes the probability of recurrence of a natural event, in this case a statistical model rainfall of a certain duration. For hydrologic/hydraulic modeling [22], the 48 classes were generalized to 13 land use classes ( Figure 2a). The history of the city of Rostock goes back to the 13th century. Since then, the city has grown not only in its territorial area but also in its population, reaching nearly 210.000 inhabitants today [23]. For Rostock, an increase of nearly 15.000 inhabitants is expected by 2035 [24]. In contrast, projections indicate that the rural district will experience a population decline of 1.8 % by 2030 [25]. However, as in many developed countries, population and sealing do not develop in parallel, but the demand for building ground is increasing, even in areas with shrinking population. Main drivers are the demand for residential units, the development of logistic centers and industrial parks.

Basic Data
In the course of the study, spatial data from different sources were used and integrated into a homogeneous and consistent data set. The data relevant for the application of the storm water routine are presented in the following.

Land Use Map
The Hanseatic City of Rostock provides a spatially high-resolution geodata set that was created on the basis of aerial photographs and supplemented by existing geodata. Since such a dataset did not exist for the rural district, it was generated from the data of the Official Real Estate Cadastre Information System (German abbreviation: ALKIS) in combination with satellite data. The result is a polygon shape file that depicts the land use in its actual state on the basis of 48 classes. In the GIS-DSS, this land use map represents the reference state for all future changes of land use. In addition, different levels of protection (risk classes) were defined for each land use class according to their vulnerability (see [22]). The level of protection is expressed by the return period which describes the probability of recurrence of a natural event, in this case a statistical model rainfall of a certain duration. For hydrologic/hydraulic modeling [22], the 48 classes were generalized to 13 land use classes (Figure 2a).

Watercourse Cadastre
The watercourse cadastre contains stream networks ( Figure 2b) and their catchments ( Figure 2c). The watercourse lines were brought together from several different sources and corrected in order to create a tree structure. Thus, each tributary flows to a receiving stream and is connected to it at a certain point. These hierarchical routes were divided into 50 m sections and equipped with an individual and hierarchical identification number (12-digit watercourse identification number: WIN) and chainage.
Each watercourse was assigned a surface catchment. These were calculated on the basis of a digital elevation model (DEM) with a cell size of 5 m. The stream catchments were also further subdivided so that each 50 m section is assigned its own subcatchment. Due to the strict hierarchy of both watercourses and catchments, which is maintained by means of the WIN, sections or subcatchments located upstream or downstream can be quickly identified and aggregated if necessary [26].

Flood Characteristics
For the watercourses and catchments of the study area, spatially high-resolution flood characteristics for the current state were preprocessed using combined rainfall-runoff and hydrodynamic stream models [22]. The flood characteristics are available for 18 precipitation events of different durations and return periods (see Table 1).

Watercourse Cadastre
The watercourse cadastre contains stream networks ( Figure 2b) and their catchments ( Figure 2c). The watercourse lines were brought together from several different sources and corrected in order to create a tree structure. Thus, each tributary flows to a receiving stream and is connected to it at a certain point. These hierarchical routes were divided into 50 m sections and equipped with an individual and hierarchical identification number (12-digit watercourse identification number: WIN) and chainage.
Each watercourse was assigned a surface catchment. These were calculated on the basis of a digital elevation model (DEM) with a cell size of 5 m. The stream catchments were also further subdivided so that each 50 m section is assigned its own subcatchment. Due to the strict hierarchy of both watercourses and catchments, which is maintained by means of the WIN, sections or subcatchments located upstream or downstream can be quickly identified and aggregated if necessary [26].

Flood Characteristics
For the watercourses and catchments of the study area, spatially high-resolution flood characteristics for the current state were preprocessed using combined rainfall-runoff and hydrodynamic stream models [22]. The flood characteristics are available for 18 precipitation events of different durations and return periods (see Table 1). The return period expresses the level of protection that has been defined for a particular land use class [27]. For example, if a return period of 100 a has been defined for a residential area, the area should not be flooded by a rainfall event that statistically occurs every 100 a. However, for reasons of economic efficiency, the area is allowed to be flooded by an even rarer and thus larger event. In order to detect the worst case, different durations were taken into account, as long rainfall durations are usually critical for large areas and shorter durations are relevant for smaller areas with shorter flow paths. Table 2 shows the preprocessed flood characteristics used in the storm water routine of the DSS.

Maximum Rainfall Intensities
Within the framework of the heavy rainfall regionalisation, the KOSTRA Atlas of the German Weather Service provides precipitation amounts as a function of duration D and annuality T (return period) [28]. While a block rain was assumed for the duration stages greater than or equal to 3 h, the intensity course for the duration stage 1 h was determined statistically based on DWA-A 118 [29] and includes a clear rain peak after about half of the time. The resulting maximum rainfall intensities (ri max ) were listed in tabular form and linked to the geodata of the hydrologic catchments.

Pre-Processing of Functions to Calculate Peak Runoff Coefficients
For the development of the storm water routine, functions were derived beforehand on the basis of the existing model results (see Table 2), which enable the calculation of the peak runoff coefficient for new land use polygons on the basis of a specified degree of sealing. In the first step, the peak runoff coefficients are determined for each subcatchment of a reference model site (Schmarler Bach site, see [22]) for the 18 simulated model rainfall events. The coefficient describes the ratio of the maximum runoff rate per area to the corresponding maximum rainfall rate per area (Equation (1)) and depends primarily on the proportion of paved areas (degree of sealing), the slope of the terrain and the rainfall intensity [29]. As the study area is located in the lowlands, the gradient plays a minor role here. This was supported by the model results, as the dependence of the coefficient on the gradient was tested and evaluated as not significant. Subsequently, a comparison of the peak runoff coefficients and the degrees of sealing was performed in a correlation diagram. Thereby polynomial functions were derived, which describe the correlation between the two parameters. The functions are used to calculate the peak runoff coefficient for a new land use polygon on the basis of the defined degree of sealing and the selected rainfall scenario.

The Storm Water Routine
The storm water routine aims to provide potential future runoff and flood risk associated with planned land use changes. For this purpose, a method was developed that enables the calculation of the direct runoff and resulting potential flows for planned land uses on the basis of pre-processed flood characteristics. The procedure takes place in the Web GIS-DSS and is divided into eight basic steps ( Figure 3). These were automated using Python programming language.
Subsequently, a comparison of the peak runoff coefficients and the degrees of sealing was performed in a correlation diagram. Thereby polynomial functions were derived, which describe the correlation between the two parameters. The functions are used to calculate the peak runoff coefficient for a new land use polygon on the basis of the defined degree of sealing and the selected rainfall scenario.

The Storm Water Routine
The storm water routine aims to provide potential future runoff and flood risk associated with planned land use changes. For this purpose, a method was developed that enables the calculation of the direct runoff and resulting potential flows for planned land uses on the basis of pre-processed flood characteristics. The procedure takes place in the Web GIS-DSS and is divided into eight basic steps ( Figure 3). These were automated using Python programming language. In the first step, a new land use polygon is either directly drawn into the map of the graphical user interface or uploaded as a shape file. In this process, the user specifies the planned degree of sealing and selects the precipitation scenario with the corresponding maximum rainfall intensity per area. For the determination of the peak runoff coefficient, preprocessed polynomial functions were used that allow a calculation based on the degree  In the first step, a new land use polygon is either directly drawn into the map of the graphical user interface or uploaded as a shape file. In this process, the user specifies the planned degree of sealing and selects the precipitation scenario with the corresponding maximum rainfall intensity per area. For the determination of the peak runoff coefficient, preprocessed polynomial functions were used that allow a calculation based on the degree of sealing. The peak runoff R max,plan for the new land use polygon can thus be determined with the equation The change in the peak runoff of the affected areas is computed by intersecting the maximum runoff of the actual land use (results from scenario simulations, see Table 2) with the maximum runoff of the new polygon both in raster format. For this purpose, the runoff rates are converted into runoff rates per raster cell (with a size of 100 m 2 ) .
in which: After intersecting the new polygon with the 50 m subcatchments, the change of peak runoff can finally be aggregated by segment ID and thus assigned to each receiving stream segment. There are two possibilities for the allocation of the runoff: (a) If no receiving segment is specified, the discharge is automatically distributed to the segments according to its area share in the subcatchment (Figure 4a,b) It is also possible to assign the discharge change to a defined segment and to throttle it via a maximum permissible discharge (Figure 4b). In the case of Figure 4b, the corresponding stream segment receives the peak runoff difference and the neighbouring upstream segment loses runoff compared to the actual state (Q red ) due to the reduced drained area.
The change in the peak runoff of the affected areas is computed by intersecting the maximum runoff of the actual land use (results from scenario simulations, see Table 2) with the maximum runoff of the new polygon both in raster format. For this purpose, the runoff rates are converted into runoff rates per raster cell (with a size of 100 m 2 ).
in which: After intersecting the new polygon with the 50 m subcatchments, the change of peak runoff can finally be aggregated by segment ID and thus assigned to each receiving stream segment. There are two possibilities for the allocation of the runoff: a) If no receiving segment is specified, the discharge is automatically distributed to the segments according to its area share in the subcatchment (Figure 4 a,b) It is also possible to assign the discharge change to a defined segment and to throttle it via a maximum permissible discharge (Figure 4b). In the case of Figure 4b, the corresponding stream segment receives the peak runoff difference and the neighbouring upstream segment loses runoff compared to the actual state (Qred) due to the reduced drained area.  The new inflow (Q diff ) is added to the maximum stream flow of the actual state (Q max,act ) and routed downstream according to the WIN From the full flow (Q full , Table 2), which represents the maximum possible flow for a specific cross profile, and the new flow of the planned state (Q max,plan ), the flow that would still fit into the profile (Q free,plan ) is determined. If the value is positive, the profile still offers enough capacity for the additional flow. If it is negative, the segment is overloaded and flooding occurs.
Q f ree,plan = f low, that would still f it in pro f ile in plan state l s −1 Q max,act = maximum f low in actual state l s −1 Both the stream flows for the 50 m segments as well as the capacities of the planned state are mapped in the GIS-DSS. For the visualisation, the line width is adjusted according to the total flow. As well, six colour levels have been assigned for the capacity to indicate how much space the profile still offers or whether it is already overloaded.

Validation of the Storm Water Routine
The validation of the storm water routine was performed by comparing its results for an exemplary land use change with the corresponding results from the combined hydrological/hydraulic model in SWMM. SWMM was considered to be a suitable benchmark, as it takes water body shapes at individual points into account and is able to calculate precisely dynamics of the flow, including effects of throttling and retention [22].
The exemplary land use change is the conversion of almost 70.000 m 2 of an orchard into a virtually planned residential area ('planned area') with a degree of sealing of 70%. The site is located in the catchment of the Schmarler Bach in the north-west of the investigation area (Figure 5a). Within the catchment area, the site is located in the southern-central part, just beyond the city boundary in the administrative area of the rural district (Figure 5b).
The discharge point into the river is set to a user-defined stream segment (red arrow in Figure 5c) as described in Figure 4b. Hence, the new residential area will drain in northern direction towards the open ditch just below the culvert (Figure 5b,c). On the flow path towards the mouth into the Warnow River, numerous structures such as culverts and pipes act as throttling elements in particular at high flows.
For the validation, representative precipitation scenarios were tested (cf. Table 1). The comparison between SWR and hydrological/hydraulic model was conducted for peak runoff from the development site, peak flows in the affected river sections.

Derived Functions for the Determination of Peak Runoff Coefficients
Using the simulation results for the maximum surface runoff and the maximum rainfall intensity, the peak runoff coefficients for the current state were calculated according to Equation (1). Figure 6 shows the correlation of the peak runoff coefficient and the degree of sealing. Functions for the determination of the peak runoff coefficients on the basis of the degree of sealing are displayed in Table 3. These are used in the storm water routine to determine peak discharges for planned land use changes.

Derived Functions for the Determination of Peak Runoff Coefficients
Using the simulation results for the maximum surface runoff and the maximum rainfall intensity, the peak runoff coefficients for the current state were calculated according to Equation (1). Figure 6 shows the correlation of the peak runoff coefficient and the degree of sealing. Functions for the determination of the peak runoff coefficients on the basis of the degree of sealing are displayed in Table 3. These are used in the storm water routine to determine peak discharges for planned land use changes. Table 3. Functions relevant for the study area for the determination of the peak runoff coefficient Ψ p (y) on the basis of the degree of sealing (x).

Return Period 2a
10a 25a 50a 100a  Basically, a clear polynomial dependence of the peak runoff coefficients on the degree of sealing is visible. This is also confirmed by the coefficients of determination in Table 4. In the case of block rainfall (3 h to 12 h), the strength of the dependency increases with duration, as conditions become more and more static. Within a duration category, the dependency tends to decrease from small to large return periods (respectively small to large rain intensities). For the duration 1 h with a clear intensity gradient, the scatter of the peak runoff coefficients around the function curve is significantly greater than for the block rainfall, which suggests that other factors, such as the area size or the flow length or the soil conductivity, superimpose the influence of the degree of sealing. In the case of the most severe event (1h-100a), it is noticeable that even the unsealed areas react as if they were sealed, since the conductivity of the soil is increasingly exceeded.  Basically, a clear polynomial dependence of the peak runoff coefficients on the degree of sealing is visible. This is also confirmed by the coefficients of determination in Table 4. In the case of block rainfall (3 h to 12 h), the strength of the dependency increases with duration, as conditions become more and more static. Within a duration category, the dependency tends to decrease from small to large return periods (respectively small to large rain intensities). For the duration 1 h with a clear intensity gradient, the scatter of the peak runoff coefficients around the function curve is significantly greater than for the block rainfall, which suggests that other factors, such as the area size or the flow length or the soil conductivity, superimpose the influence of the degree of sealing. In the case of the most severe event (1 h-100 a), it is noticeable that even the unsealed areas react as if they were sealed, since the conductivity of the soil is increasingly exceeded.

Comparison of Model Results for Actual and Plan State in SWMM
In order to better interpret and evaluate the results of the GIS-DSS storm water routine, this section first compares the model results for the current and the planned state. To start with a less extreme event, Figure 7 shows the difference in flows (Q max,plan − Q max , act ) in the planned and actual state for the rainfall event 6h-10a simulated with SWMM. The turquoise and yellow sections show on which stretches the discharge change has an effect.

Comparison of Model Results for Actual and Plan State in SWMM
In order to better interpret and evaluate the results of the GIS-DSS storm water routine, this section first compares the model results for the current and the planned state. To start with a less extreme event, Figure 7 shows the difference in flows (Qmax,plan -Qmax,act) in the planned and actual state for the rainfall event 6h-10a simulated with SWMM. The turquoise and yellow sections show on which stretches the discharge change has an effect. As expected, the largest differences in stream flow occur immediately downstream of the discharge point and are of the same magnitude as the peak runoff from the planned area (Figure 8). At the same time, the new runoff also affects an approx. 130 m long stretch upstream of the discharge point but with a less prominent amplitude. As expected, the largest differences in stream flow occur immediately downstream of the discharge point and are of the same magnitude as the peak runoff from the planned area (Figure 8). At the same time, the new runoff also affects an approx. 130 m long stretch upstream of the discharge point but with a less prominent amplitude.
With downstream distance from the discharge point, the flow difference flattens from 0.071 m 3 s −1 to 0.062 m 3 s −1 , which can be explained by the natural retention of an open channel. A pipeline section follows, which is about 95% full in the actual state, but already overloaded in the planned state. As a result, the flows are noticeably throttled here, which is why the additional discharge of the new residential area from this section onwards has relatively little effect.
A much more intense event is the 1h-100a rainfall (Figures 9 and 10). Here, 0.82 m 3 s −1 of runoff is generated on the new residential site. Since the profile downstream can only capture 0.42 m 3 s −1 , the conduits capacity limit is already reached. However, due to the rise in water level in the area of the discharge point and the resulting inversion of the upstream gradient, water flows in both directions: 0.42 m 3 s −1 flows downstream and approx. 0.26 m 3 s −1 flows upstream ( Figure 11).
The remaining discharge (approximately 0.14 m 3 s −1 ) is excess water and is temporarily stored in a conceptual reservoir above the banks. In SWMM, this is called the ponded volume. The speed with which the water flows out of the reservoir back into the crosssection depends on the size of the base area of the reservoir. If the surface area is small, the hydraulic pressure in the reservoir is relatively high and the water flows back again correspondingly quickly. This in return influences the shape and the maximum value of the stream flow. Since the study area is located in the lowlands, the footprint of the reservoir was chosen to be relatively large (distance to the next cross section below multiplied by 40 m). However, this conceptual approach has uncertainties from the moment flooding occurs. With downstream distance from the discharge point, the flow difference flattens from 0.071 m 3 s −1 to 0.062 m 3 s −1 , which can be explained by the natural retention of an open channel. A pipeline section follows, which is about 95% full in the actual state, but already overloaded in the planned state. As a result, the flows are noticeably throttled here, which is why the additional discharge of the new residential area from this section onwards has relatively little effect.
A much more intense event is the 1h-100a rainfall (Figures 9 and 10). Here, 0.82 m 3 s −1 of runoff is generated on the new residential site. Since the profile downstream can only capture 0.42 m 3 s −1 , the conduits capacity limit is already reached. However, due to the rise in water level in the area of the discharge point and the resulting inversion of the upstream gradient, water flows in both directions: 0.42 m 3 s −1 flows downstream and approx. 0.26 m 3 s −1 flows upstream ( Figure 11).   With downstream distance from the discharge point, the flow difference flattens from 0.071 m 3 s −1 to 0.062 m 3 s −1 , which can be explained by the natural retention of an open channel. A pipeline section follows, which is about 95% full in the actual state, but already overloaded in the planned state. As a result, the flows are noticeably throttled here, which is why the additional discharge of the new residential area from this section onwards has relatively little effect.
A much more intense event is the 1h-100a rainfall (Figures 9 and 10). Here, 0.82 m 3 s −1 of runoff is generated on the new residential site. Since the profile downstream can only capture 0.42 m 3 s −1 , the conduits capacity limit is already reached. However, due to the rise in water level in the area of the discharge point and the resulting inversion of the upstream gradient, water flows in both directions: 0.42 m 3 s −1 flows downstream and approx. 0.26 m 3 s −1 flows upstream ( Figure 11).   The remaining discharge (approximately 0.14 m 3 s −1 ) is excess water and is temporarily stored in a conceptual reservoir above the banks. In SWMM, this is called the ponded volume. The speed with which the water flows out of the reservoir back into the crosssection depends on the size of the base area of the reservoir. If the surface area is small, the hydraulic pressure in the reservoir is relatively high and the water flows back again correspondingly quickly. This in return influences the shape and the maximum value of the stream flow. Since the study area is located in the lowlands, the footprint of the reservoir was chosen to be relatively large (distance to the next cross section below multiplied by 40m). However, this conceptual approach has uncertainties from the moment flooding occurs.  The remaining discharge (approximately 0.14 m 3 s −1 ) is excess water and is temporarily stored in a conceptual reservoir above the banks. In SWMM, this is called the ponded volume. The speed with which the water flows out of the reservoir back into the crosssection depends on the size of the base area of the reservoir. If the surface area is small, the hydraulic pressure in the reservoir is relatively high and the water flows back again correspondingly quickly. This in return influences the shape and the maximum value of the stream flow. Since the study area is located in the lowlands, the footprint of the reservoir was chosen to be relatively large (distance to the next cross section below multiplied by 40m). However, this conceptual approach has uncertainties from the moment flooding occurs.

Comparison of Storm Water Routine Results with Model Results for the Plan State
Regarding peak runoff, the deviations of the results achieved by the SWR from those of the model tend to be larger the shorter the event is. As validation criterion, SWR results for events with a relative deviation smaller than 10% were classified as 'good'. As shown in Table 5, this applies to all SWR peak runoffs with a duration of 3 h and longer which have a coefficient of determination higher than 0.9.
The two 1 h events, whose triggering rains show a varying intensity course, deviate considerably with 30% and 32%, respectively. It should be noted here that more accurate results for those events would have been expected if the degree of sealing for the planned area had been <50% (cf. correlation for 1 h events in Figure 6). However, the peak runoffs for the 1 h events calculated with SWR can be taken as the mean area response if other parameters (gradient, flow length, roughness, etc.) are still unknown in the process of planning.  Figure 12 illustrates the effect of the user-defined discharge point. The sealed residential area overlaps four subcatchments where in this precipitation scenario (6 h-10 a) hardly any runoff was generated in the initial state. In the planning scenario, the surface runoff increases by 0.071 m 3 s −1 in the selected subcatchment (highlighted in yellow), while in the other three subcatchments the runoff slightly decreases or remains unchanged. During more extreme events such as 1h-100a (Figure 13), runoff also occurs on the unsealed subcatchments in the initial state. As the planned area drains only into the central user-defined segment (+1.169 m 3 s −1 ), the resulting difference in runoff of adjacent and intersected subcatchments is negative in the planning scenario (blue and green areas in Figure 13). During more extreme events such as 1h-100a (Figure 13), runoff also occurs on the unsealed subcatchments in the initial state. As the planned area drains only into the central user-defined segment (+1.169 m 3 s −1 ), the resulting difference in runoff of adjacent and intersected subcatchments is negative in the planning scenario (blue and green areas in Figure 13). During more extreme events such as 1h-100a (Figure 13), runoff also occurs on the unsealed subcatchments in the initial state. As the planned area drains only into the central user-defined segment (+1.169 m 3 s −1 ), the resulting difference in runoff of adjacent and intersected subcatchments is negative in the planning scenario (blue and green areas in Figure 13).  When comparing SWR-calculated peak flows with those simulated by SWMM the results directly downstream of the discharge point (300 m reach) are nearly identical (cf. Figure 12). With increasing distance, however, the differences become larger. This can be explained by different hydraulic effects such as wave flattening and retention in the hydraulic model, which are not implemented in the simplified routing algorithm of the SWR. The same applies for the deviations between SWR and model directly upstream of the discharge point. Due to backwater, a part of the runoff flows upstream for a short time in the model (cf. Figure 11), which the SWR does not take into account. After heavy rainfall events, the deviations of the peak flows can be even larger than the indicated runoff difference from Table 5, due to the flow redistribution (flow direction downstream and upstream), which the SWR cannot consider.
Additionally, the different spatial resolutions of both approaches have an impact here: In order to be consistent with other data sets of the DSS, the SWR routes the peak runoff from the discharge point along river sections of 50 m length. By contrast, the hydraulic model in SWMM uses river segments of varying length to address any relevant changes (profile, infrastructure).
Despite the mentioned simplifying approaches, the SWR enables a rapid evaluation of the impact of land use change scenarios on flood risk. Figure 14 shows an exemplary result for flow capacities (Q max, SWR /Q full * 100 %) calculated by the SWR for the 6 h-10 a event. The colors of the stream segments indicate how critical its water level could rise in the planning scenario. Pink sections mark the areas that could potentially be affected by flooding, especially if the capacity is greater than 100%. changes (profile, infrastructure).
Despite the mentioned simplifying approaches, the SWR enables a rapid evaluation of the impact of land use change scenarios on flood risk. Figure 14 shows an exemplary result for flow capacities (Qmax, SWR/Qfull * 100 %) calculated by the SWR for the 6h-10a event. The colors of the stream segments indicate how critical its water level could rise in the planning scenario. Pink sections mark the areas that could potentially be affected by flooding, especially if the capacity is greater than 100 %.

Discussion
Before discussing the pros and cons of total approach, the background and overall objectives shall be shortly remembered: The developed approach shall enable regional planners, water authorities, and other stakeholders to assess the potential impact of planned land use changes on runoff, peak flow, and related flood risk in the receiving river system. Focus water bodies are small lowland rivers, often without any monitoring

Discussion
Before discussing the pros and cons of total approach, the background and overall objectives shall be shortly remembered: The developed approach shall enable regional planners, water authorities, and other stakeholders to assess the potential impact of planned land use changes on runoff, peak flow, and related flood risk in the receiving river system. Focus water bodies are small lowland rivers, often without any monitoring data for model calibration. The often missing modeling expertise of the stakeholders is addressed replacing detailed modeling by GIS based analysis.
Setup and parametrization of a detailed hydrologic/hydrodynamic model 2.
Forecasting runoff change induced by land use changes and downstream flood risks applying newly developed GIS routines For consistency, here, we shortly repeat the main findings of step 1. For the combined representation of rainfall-runoff and hydrodynamic streamflow processes, we selected the free-ware model SWMM-UrbanEVA. Since the model is physically based, it can be strictly setup and parametrized using geodata. We proposed an according workflow, including the definition of important model parameters, for lowland river basins. The approach was shown and tested for two river basins representing two extremes of approximately 20 superior catchments in the northeast German lowlands. Even though all are located in the same climatic and geo-morphologic zone, they differ in some characteristics (e.g., average terrain slope, depth to groundwater table, dominant land use, etc.). For the here interesting peak runoff, the most hydrologically significant attribute is land use, specifically the degree of urbanity or the percentage of sealed surfaces. It was found that the model performed better at stronger urbanised river basin. As discussed there, the lower model performance is only partly due to the parametrization but rather to uncertain rain data. Though, it can be postulated that the object data of urbanised areas are in most cases better documented and the rainfall-runoff process is less complex, leading to an expectable better model performance compared to near natural river basins. Despite the proposed strict and highly automated workflow based on geodata, model setup still requires an experienced modeler.
For the flood risk analysis, model based or with the GIS tool, distinct precipitation scenarios have to be chosen. Here, we propose synthetic rains, characterised by reoccurrence interval and duration. The reoccurrence intervals have to be chosen according to the desired protection level of the different land uses. The advantages of this approach are short computation time and the availability of regionalised rain data. However, the approach is very sensitive to the initial condition of the model, when precipitation starts. Here, we preprocessed the model using rain times series to produce average flow. More detailed analysis of rain statistics may lead to more differenced settings for the initial condition according to the precipitation scenario.
Once the model exists, the analysis of land use changes can be performed by any person with basic GIS skills. As the above discussed results show, the simplified calculation of runoff, based on the rain intensity dependent peak runoff functions matches very well with the results of the physical model. Since only a static peak runoff is applied, the calculation accuracy decreases with increasing impact of dynamics, as seen here for the 1 h precipitation scenario. The same applies for the routed stream runoff in river. For rain duration stages between 6 h to 12 h, the results of the benchmark model are met very well (R 2 > 096), while dynamic effects for short duration stages (wave flattening, backwater effects) cannot be addressed by the simplified routing algorithm. Generally, the routing tends to overestimate downstream flow effects and is, regarding the impact assessment, on the safe side.
The method presented here was specifically tailored to the study area and the available data, but if similar data sets are available, it is also very well transferable to other sites. It should be noted that the proposed peak runoff functions have been derived for the characteristics of rather flat morphology and short distances between regarded area and discharge point. In other regions or spatial conditions, the terrain gradient or the flow length may lead to deviating peak runoff coefficients. Case specific functions can be derived as described above.
A great advantage of the method is the evaluation of flood risk for planners without the setup of a complex model, as proposed by [12].
Compared to earlier approaches, assessing runoff induced by land us change in the past ( [3,4]), the SWR can be applied to predict runoff from future settlement areas. Additionally, not only runoff, but the actual flood risk is assessed.
For this, the method goes in two dimensions one step further: i) small scale land use change can be directly allocated and assessed; ii) the spread of runoff in the watercourse is calculated by incorporating high resolution DEM-derived stream geometry and characteristics. Considering that the method is to be used primarily in early planning phases, when it is perhaps just known that an industrial area of a certain size is to be built in a district, for example, it is a great progress if the stream hydraulic impact can be calculated immediately without prior knowledge. The method quickly shows tendencies where the water is overflowing and where additional measures may be necessary.
Following the generic framework of [15], the tool can be regarded partly as collaborative and partly as participatory modeling approach. Arguments for real collaborative modeling are:

•
The modeling concept has been developed jointly • The data for the model are accessible for or provided by the decision-makers and stakeholders • The tool is designed interactively and embedded in a familiar GIS environment • Results are processed and visualized for direct use (interactive planning, decision on storm water discharge applications, etc.) However, the underlying process models set up with SWMM-UrbanEVA will require regular updates (changes of assets, land-use). Still, this requires an expert modeler. Though the highly automated model setup strictly based on geodata (see part 1) simplifies model updates significantly, this task should remain in one hand to avoid inconsistent model versions at the various partners.
Currently, the tool provides one-dimensional information of flood risk. For experienced users, this is already highly informative. However, these data are less suitable for dissemination to citizens, since interpretation of an actual precise flood risk at a certain point would require well-considered 2D flood maps [17] and in best case a clear instruction, what a certain flood level means for the property and which response measures should be taken [30]. This applies even more, since information on flood risk have clear financial consequences for property owners [31]. Therefore, the main purpose of this tool is not to warn citizens, but to avoid flood risks by interactive sustainable land use planning and well-reasoned definition of tolerable storm water discharges.
In validating the results produced using the SWR tool, citizen science can be an important component (cf. [32,33]). For example, the tool shows us stream segments that potentially spill during heavy rain events. Interviewing eyewitnesses to past flooding events can help confirm potential problem areas. More automatable is the evaluation of social media comments on flooding events (e.g., Twitter) [34]. Therefore, the latter is particularly suitable for large-scale observations.

Conclusions and Outlook
This study presents a procedure that uses a simple method to calculate the impact of planned land use changes on river flow and potential flooding. It combines a detailed hydrologic and hydrodynamic model set up for the status quo with a simplified forecast of direct runoff in the planning scenario, which is routed through the downstream watercourse network.
The SWR forecast tool is easy and intuitive to use and does not necessarily demand modeling experience. It requires only the drawing or uploading of one or more new land use polygons, the planned degree of sealing and the precipitation scenario of interest. The calculations in the GIS-DSS take only few seconds, so that several scenarios can be tested and compared efficiently-e.g., size and location of area, degree of sealing or discharge point. The advanced visualization of the results, for example the display of overloaded river sections, improves communication with stakeholders enormously [33].
The underlying physically based process model can be set-up highly automated based on available geodata. For comparable catchment characteristics a parameter transfer from gauged to ungauged catchments is possible with tolerable loss of model accuracy.
Nevertheless, the described method is based on simplifications that may deviate more or less from reality. Three main sources of uncertainty can be mentioned here: -Uncertainty of the input data (object data, time series of rainfall and flow) for set-up and calibration/parametrization of the process model for the current state (in the investigated river basins the resulting peak runoff error was 8-26 %) -Inaccuracies in the simplified calculation of the peak runoff of a newly planned site: 0-32 % (depending on rain scenario) -Inaccuracies in the propagation of additional peak runoff in the watercourse: 0-1.07 m 3 s −1 downstream of discharge point and 0.06-1.05 m 3 s −1 upstream (depending on rain scenario) A flood event is a transient process, whereas the applied calculation method is stationary. Therefore, the presented method performs better for block rains of larger durations. Generally, the deviations between SWR-calculated and simulated peak flows increase with distance from the discharge point downstream due to wave attenuation and flow retention. However, this is not to be seen as a shortcoming-rather, the calculated maximum flows represent the potential changes in planned state and thus have a warning effect on the user. In the exemplary case, downstream of the imaginary discharge point, an enlargement of the pipeline was discussed, which would reduce the retention effect and might produce flooding in the lower reach.
Although the results of the SWR already provide the user with a great deal of information gain, several suggestions for improvement can be derived.
The implementation of a retention factor per stream segment would be conceivable for the routing of peak flows. These could be calculated directly from the model results as a percentage of the respective upstream flow. It would also be possible and simple to implement to stop the routing as soon as the capacity is exceeded at a point in the watercourse (typically a culvert). However, this 'improvement' has to be taken with care since it cannot be ruled out that in reality the water continues to flow on the surface of the terrain and enters the stream again at another point.
As a next step, the intersection of the flooded stream segments with the surrounding land uses (and their required level of protection/return period) would also be conceivable, e.g., by means of a 50 m buffer around the segments. The land use polygons where the required level of protection cannot be achieved could subsequently be marked with a warning color.
What has not been addressed, yet, is the possibility of using the SWR tool for low impact design measures (LID) [35]. The required conceptual approach to link LID measures with direct runoff calculations could be derived from scenario modeling using SWMM-UrbanEVA. This upgrade of SWR would allow a fast preliminary design and allocation of LID to mitigate existing or potential flooding situations.