Coupling Field Observations and Geographical Information System (GIS)-Based Analysis for Improved Sustainable Urban Drainage Systems (SUDS) Performance

: Urbanization and increased precipitation volumes and intensities due to climate change add pressure to the urban drainage system, resulting in increased ﬂooding frequencies of urban areas and deteriorating water quality in receiving waters. Inﬁltration practices and the use of blue green infrastructure, also called Sustainable Urban Drainage Systems (SUDS), can limit, and, in some cases, reverse the effects of urbanization. However, adequate inﬁltration capacity is an essential parameter for the successful implementation. In this paper, a Geographical Information System (GIS)-based hydrology analysis for SUDS placements is coupled with ﬁeld measurements using Modiﬁed Phillip Dunne inﬁltrometer tests. The case study area is the expansion of the campus at the Norwegian University of Science and Technology (NTNU) over the next decade. Inﬁltration in urban soils can be highly heterogenous over short distances. When comparing measured inﬁltration rates with physical characteristics of the soils showed that the physical characteristics are not a good indication of the inﬁltration potential in urban soils with a large degree of compaction. The results showed that measuring the inﬁltration potential combined with ﬂow path analysis can greatly enhance the beneﬁts of blue green infrastructure, with an up to 70% difference in area required for SUDS solutions for managing 90% of the annual precipitation.


Introduction
The performance of sustainable urban drainage systems (SUDS), such as rain gardens, swales, and bioretention areas are driven mainly by the infiltration rate, the process in which water enters the soil media. The infiltration process may be described by the saturated hydraulic conductivity, K sat , which measures the rate at which water can enter into the soil media at a given viscosity (temperature dependent) and soil permeability [1]. Soil maps can be used to obtain an estimate through soil type, which again can give parameters, such as K sat and permeability. However, in urban areas the heterogeneity of soils is a major issue in using such maps. Soils in urban areas can vary significantly from soils with the same classification in natural non urban areas due to; (1) the degree of compaction over time due to construction and heavy loads; (2) the amount of organic matter; and, (3) contamination from construction debris and unspecified construction fill [2][3][4][5]. This makes it challenging to classify urban soils in the normal soil taxonomy groups, and assume water infiltration properties that are based on values given in tables.

Materials and Methods
In this study, a method that couples field measurements of infiltration protentional with GIS analysis of the flow directions and accumulation into drainage lines is introduced as a case study at the NTNU campus in Trondheim, Norway. First, the method for the GIS analysis is described, followed by the field Modified Phillip Dunne (MPD) measurements and subsequent data analysis. Before the combined new method coupling the analysis of the two data sets is described.
The NTNU campus is a centrally located area at Gløshaugen and Lerkendal in the city of Trondheim, Norway. The campus area is more than 50 ha large in total. The area has a high degree of impervious surfaces, comprising buildings, parking lots, and roads. In total, the sealed surfaces make 58% of the total area, Figure 1. There are also large areas of park and undeveloped green space in the catchment. The surface slope changes from flat areas in the central core to hilly parts with up to 60% surface slope along the edges. The studied area has its highest point at the Gløshaugen plateau, from where the elevation decreases in all directions. Over the next decade, the campus building mass will be expanded to house several parts of the university that today is located across the city at other campuses. This will add up to 92,000 m 2 of building mass in one or more clusters of buildings in parts of what is currently undeveloped areas. The Campus Development Project has a strong focus on development in line with the SDGs and preserving green spaces, which makes it relevant as a reference for green campus initiatives.
Sustainability 2018, 10, x FOR PEER REVIEW 3 of 13 will be expanded to house several parts of the university that today is located across the city at other campuses. This will add up to 92,000 m 2 of building mass in one or more clusters of buildings in parts of what is currently undeveloped areas. The Campus Development Project has a strong focus on development in line with the SDGs and preserving green spaces, which makes it relevant as a reference for green campus initiatives. Figure 1. The Norwegian University of Science and Technology (NTNU) campus area at Gløshaugen and Lerkendal in Trondheim showing the contour lines with elevations in meters above mean sea level, undeveloped areas in green, roads and parking lots in white, and buildings in grey.

GIS Analysis
The GIS analysis is performed by the open source program QGIS ® v. 2.18.13 with implementations from the open source program SAGA© GIS v. 2.3.2. The output of the GIS analysis are several maps, which show urban planners potential and high potential SUDS areas, as well as main flow paths for surface runoff. Complete reference guidelines, including files of the developed tools to perform the described GIS analysis, can be downloaded from www.klima2050.no.
The GIS process is illustrated in Figure 2. The input data are a digital elevation map (DEM) (

GIS Analysis
The GIS analysis is performed by the open source program QGIS ® v. 2.18.13 with implementations from the open source program SAGA© GIS v. 2.3.2. The output of the GIS analysis are several maps, which show urban planners potential and high potential SUDS areas, as well as main flow paths for surface runoff. Complete reference guidelines, including files of the developed tools to perform the described GIS analysis, can be downloaded from www.klima2050.no.
The GIS process is illustrated in Figure 2. The input data are a digital elevation map (DEM) (  Illustrated geographical Information Systems (GIS) process, (a) input digital elevation model (DEM) (white = high to black = low), (b) input data orthophotos, (c) calculated flow accumulation (blue = low to red = high), (d) calculated slope (green = low to red = high), (e) determined undeveloped land (white), (f) selected main flow paths (red), (g) resulting potential sustainable urban drainage systems (SUDS) areas (green), (h) resulting high potential SUDS areas (green) next to main flow paths (red).
To ensure uninterrupted flow lines when calculating the main flow paths, local surface depressions in the DEM are filled, using a sink filling algorithm [11]. Flow accumulation shows for each cell in a DEM raster the number of cells that flow lines pass through that grid cell [12]. This means that the number of flowlines from the entire map that have segment endpoints within that cell will be shown for each DEM raster cell. The flow direction is calculated using the Rho8 algorithm [13] and the flow direction calculated for the steepest downslope neighbor. The Rho8 algorithm is well suited for steep topography [13], typically found in many Norwegian catchments, which is also the case for the case study presented here. To ensure uninterrupted flow lines when calculating the main flow paths, local surface depressions in the DEM are filled, using a sink filling algorithm [11]. Flow accumulation shows for each cell in a DEM raster the number of cells that flow lines pass through that grid cell [12]. This means that the number of flowlines from the entire map that have segment endpoints within that cell will be shown for each DEM raster cell. The flow direction is calculated using the Rho8 algorithm [13] and the flow direction calculated for the steepest downslope neighbor. The Rho8 algorithm is well suited for steep topography [13], typically found in many Norwegian catchments, which is also the case for the case study presented here.
To identify the highest flow concentrations from the flow accumulation results, a threshold for the main flow paths was calculated. The root mean square flow accumulation is calculated according to Equation (1), while considering all flow accumulation values: where x n is the flow accumulation value in grid cell n and x RMS is the root mean square used as a flow accumulation threshold. Higher values are given increased weight by the RMS, due to the squaring of each value. This is necessary since there are many low values that are slightly larger than 1 (which is not of interest for the calculation) and only a few large values of interest for the calculation. All cells with a larger flow accumulation than the root mean square are defined as the main flow paths. These paths represent the connection of many flow paths, meaning that the majority surface runoff will follow these paths. Furthermore, the potential SUDS areas were defined as undeveloped areas having a slope less than 10%, whereas high potential areas in addition should be located close to main flow paths and have an elevation that is equal to or lower than the nearby runoff areas. More details about the criteria for the potential SUDS areas and high potential SUDS areas are given in Section 2.4.

Field Measurements
Field infiltration measurements were conducted using the MPD Infiltrometer method [14,15]. This single ring infiltrometer test was performed using a 50 cm tall, 110 mm inner diameter acrylic tube. A measuring tape was attached to the tube in order to read off water levels as the water infiltrated into the under lying soils. The tube was placed five centimeters into the soil surface and filled with 40 cm of water (3.8 L). The water level was read off at 2 min intervals the first 30 min and then on a 5 min interval until all the water had infiltrated. The reading was done by reading off the height of water using the measuring tape glued on to the tube and a stop watch. Prior to the infiltration test, a soil sample of a minimum 1 L of soil was collected close to the location of the tube, this was used to measure initial soil moisture in the soil. After the infiltration test, a soil sample of minimum 1 L of soil was collected where the tube had been placed. This latter sample was used to measure the soil moisture after the infiltration test. In addition to the manual method described above, an automated MPD Infiltrometer (Upstream Technologies, New Brighton, USA) was used for some of the samples. The theory behind and the procedure is the same as for the fully manual method, however for the automated method the change in head is automatically recorded based in a water depth/pressure measurement. The infiltration rate is further automatically calculated for every 200 s using the same method, as described below.
All field samples were done with three replicates at each site location. The final K sat was calculated as an average of the three measured values. The hydraulic conductivity is given by: where: K = hydraulic conductivity (m/s), k = permeability (m 2 ), µ = dynamic viscosity (kg/m·s), ρ = density of water (kg/m 3 ), and g = gravitational acceleration (m/s 2 ). The manual field measurements were post processed using a MATLAB scrip to compensate for the horizontal flow direction of the water as it enters the soils [16].

Laboratory Analysis
The collected soil samples were used to determine the soil moisture content before and after the test and for particle size sieving. The samples were first weighed followed by oven drying at 105 • C for 24 h. Following the oven drying, the samples were weighed again. The difference in wet and dry weight gives the water weight of the sample, which can be used to calculated moisture content by mass using Equation (3): where: θ is the soil moisture (%); W wet and W dry are the wet and dry weight of the soil sample, respectively, in grams. The dried soil samples were then sieved in a AS200 sieving machine with the following fractions; 4 mm, 2 mm, 1 mm, 500 µm, 250 µm, 125 µm, 63 µm, and 45 µm. This produced soil classifications curves for each sample. A minimum of 350 g dried soil was used for the sieving.

Potential SUDS Area Detection
Potential SUDS areas were identified based on the following criteria; first, by locating SUDS to areas where the flow path naturally accumulates water it minimizes the regrading needed to collect the surface runoff. This provides the opportunity to treat stormwater effectively using infiltration based SUDS with minimal regrading of the terrain. Secondly, priority was given to areas with an elevation of the potential SUDS areas equal or lower than the nearby runoff areas, enabling natural drainage and avoiding costly modification of nearby terrain. Areas fulfilling these two criteria were identified as high potential SUDS areas, as natural areas where water concentrated in the flow path and has available space for SUDS implementation at a favorable elevation, resulting in areas highly recommended for SUDS implementation.
Two main criteria for SUDS placement were identified; first, the location in connection to flow path of area to maximize the potential efficiency of the system; second, the infiltration capacity in the native soils, which determines the sites' suitability for infiltration based SUDS and to the extent possible using the on-site soil. The infiltration potential from the field measurement was made into an infiltration potential map using the point measurements from the field measurements to make a Thiessen polygons [17], where the areas closest in proximity to each point measurement were classified with the infiltration of the point measurement.
The GIS analysis of flow accumulation creating runoff paths was coupled with landuse data and field measurements to determine potential SUDS placements using the following algorithm. Overlaying the surface runoff flow path, the point field observation map identifies the surface runoff paths with the highest infiltration rates, which gives the largest potential for SUDS implementation. The second step overlays the landuse map to identify which of the highest potential infiltration areas are available for use. This can either be based on property ownership or by undeveloped area. Undeveloped area in this case could be unused green areas, inactive courtyards, or storage areas. Also, areas like parking lots, where a reduction of the pervious area is planned, could be "undeveloped" in the broad sense of this method.

SUDS Selection and Design Storms
Norway has adopted a three-step approach to stormwater management, where SUDS are classified into three categories.
Step 1 is infiltration based SUDS intended to infiltrate all small events, typically events with less than 1-2 year return period.
Step 2 is based on the detention of all medium events and subsequent infiltration and slow release downstream. While, the third step is focused on safe floodways for all large and extreme events [10]. The definition of the size of the storm events is left to local authorities to decide. This case study used a requirement that SUDS in step 1 and 2 should handle 90% of the annual volume, i.e., by using daily observed rainfall data the daily depth of rainfall that will only be exceeded 10% of the time. Using data from 2013-2017 from the nearby Risvollan station (data available for download at www.eklima.no) it was determined that on a daily basis the SUDS would have to handle 19.3 mm of precipitation in order to manage 90% of the annual volume. In other words, a design specification of 19.3 mm for the specific SUDS, like a bioretention area or a combination of SUDS, like green roofs coupled with bioretention areas, would infiltrate 90% of the annual precipitation, and only 10% would become overflow from the system. The calculation uses a daily timestep, meaning that if the 19.3 mm came in a shorter duration event, the infiltration capacity would be reached faster. Trondheim is a wet coast area with frequent but typically low intensity events, which makes this assumption reasonable.
The design of the bioretention cells was sized using Darcy's law [4]. The design recommendation for raingardens in Norway is given as 10 cm/h [18,19], which is significantly higher than the commonly recommended 3 cm/h in more temperate climates. This is due to the need to maintain a sufficient infiltration capacity also in the winter months with partly frozen soils.
where: A (m 2 ) equals the needed surface area of the bioretention area, A drain equals the total drainage area draining into the bioretention area (m 2 ), C is the rational formula runoff coefficient, P is the precipitation (m), t f is the time to drain (h), h f is the average height of water above the filter bed (m), d f is the depth of the filter media (m), and k is the saturated hydraulic conductivity (m/h). The second group of SUDS chosen in this study was green roofs. Green roofs have an array of positive benefits; including a positive effect on stormwater detention and retention metrics, serves to protect roof membranes for UV exposure, reduce temperature fluctuations in the roof, and reduce energy consumptions, especially in warmer climates [20][21][22][23][24][25]. The performance of green roofs can be divided into detention and retention performance. Detention performance is an event based metrics, while the retention performance is a longer term performance metric that is dependent on the evapotranspiration function of the roof [19,25,26]. In this study, green roofs were used to show case the potential performance enhancing effect of the overall SUDS system performance on an annual basis, when green roofs and bioretention areas are coupled. This leads to a focus on the retention performance of the green roofs. The retention performance was based on the retention performance from studies by Johannessen et al. [19], where Trondheim was one of the field measurement locations.

Drainage Lines
The estimated drainage lines on campus are shown in Figure 3 and clearly show the hill shaped area. The main drainage lines are away from the central high point down relatively steep slopes to the lower laying areas around the outer edges of the case studies. The landuse in the area (cf. Figure 1) has most of the impervious surface areas located on the higher plateau of the case study areas. would be reached faster. Trondheim is a wet coast area with frequent but typically low intensity events, which makes this assumption reasonable. The design of the bioretention cells was sized using Darcy's law [4]. The design recommendation for raingardens in Norway is given as 10 cm/h [18,19], which is significantly higher than the commonly recommended 3 cm/h in more temperate climates. This is due to the need to maintain a sufficient infiltration capacity also in the winter months with partly frozen soils.
where: A (m 2 ) equals the needed surface area of the bioretention area, Adrain equals the total drainage area draining into the bioretention area (m 2 ), C is the rational formula runoff coefficient, P is the precipitation (m), tf is the time to drain (h), hf is the average height of water above the filter bed (m), df is the depth of the filter media (m), and k is the saturated hydraulic conductivity (m/h). The second group of SUDS chosen in this study was green roofs. Green roofs have an array of positive benefits; including a positive effect on stormwater detention and retention metrics, serves to protect roof membranes for UV exposure, reduce temperature fluctuations in the roof, and reduce energy consumptions, especially in warmer climates [20][21][22][23][24][25]. The performance of green roofs can be divided into detention and retention performance. Detention performance is an event based metrics, while the retention performance is a longer term performance metric that is dependent on the evapotranspiration function of the roof [19,25,26]. In this study, green roofs were used to show case the potential performance enhancing effect of the overall SUDS system performance on an annual basis, when green roofs and bioretention areas are coupled. This leads to a focus on the retention performance of the green roofs. The retention performance was based on the retention performance from studies by Johannessen et al. [19], where Trondheim was one of the field measurement locations.

Drainage Lines
The estimated drainage lines on campus are shown in Figure 3 and clearly show the hill shaped area. The main drainage lines are away from the central high point down relatively steep slopes to the lower laying areas around the outer edges of the case studies. The landuse in the area (cf. Figure  1) has most of the impervious surface areas located on the higher plateau of the case study areas.

Field Infiltration Measurements
The infiltration capacity has been measured at a total of 20 sample sites, with three replicas at each site. The case study area is an area where there have been several construction stages, making it a typical urban area. From the sieving analysis of the soil samples, all samples were found to have a high sand content; in the range of 46-80% of the particles were in the sand particle range (0.05-2 mm), the complete sieving analysis results are attached in the supplementary information (Table S1). The predominant soil types were found to be sandy loams and loamy sands according to the soil texture triangle. K sat values ranged from less than 1 cm/h to over 20 cm/h (Figure 4). When comparing the simple mean to the geometric mean, it can be seen than while most of the replica were similar, some sampling points, like sampling point number 4, shows a larger variation between the replicas, resulting in a larger difference between the simple mean and the geometric mean. This indicates that additional measurement should be undertaken prior to SUDS placement decisions. The spatial distribution of the K sat values indicated generally higher infiltration capacities in the central areas than in the north and south ends of the case study area ( Figure 5).

Field Infiltration Measurements
The infiltration capacity has been measured at a total of 20 sample sites, with three replicas at each site. The case study area is an area where there have been several construction stages, making it a typical urban area. From the sieving analysis of the soil samples, all samples were found to have a high sand content; in the range of 46-80% of the particles were in the sand particle range (0.05-2 mm), the complete sieving analysis results are attached in the supplementary information (Table S1). The predominant soil types were found to be sandy loams and loamy sands according to the soil texture triangle. Ksat values ranged from less than 1 cm/h to over 20 cm/h (Figure 4). When comparing the simple mean to the geometric mean, it can be seen than while most of the replica were similar, some sampling points, like sampling point number 4, shows a larger variation between the replicas, resulting in a larger difference between the simple mean and the geometric mean. This indicates that additional measurement should be undertaken prior to SUDS placement decisions. The spatial distribution of the Ksat values indicated generally higher infiltration capacities in the central areas than in the north and south ends of the case study area ( Figure 5).

Field Infiltration Measurements
The infiltration capacity has been measured at a total of 20 sample sites, with three replicas at each site. The case study area is an area where there have been several construction stages, making it a typical urban area. From the sieving analysis of the soil samples, all samples were found to have a high sand content; in the range of 46-80% of the particles were in the sand particle range (0.05-2 mm), the complete sieving analysis results are attached in the supplementary information (Table S1). The predominant soil types were found to be sandy loams and loamy sands according to the soil texture triangle. Ksat values ranged from less than 1 cm/h to over 20 cm/h (Figure 4). When comparing the simple mean to the geometric mean, it can be seen than while most of the replica were similar, some sampling points, like sampling point number 4, shows a larger variation between the replicas, resulting in a larger difference between the simple mean and the geometric mean. This indicates that additional measurement should be undertaken prior to SUDS placement decisions. The spatial distribution of the Ksat values indicated generally higher infiltration capacities in the central areas than in the north and south ends of the case study area ( Figure 5).

SUDS Implementation Potential
Looking at the existing development of the campus case study area, the pervious area makes up 44% of the total area. This is area that is currently available for SUDS implementation. Choosing three different strategies for the implementation of SUDS in this case demonstrates the importance of coupling field measurements with GIS based analysis.

GIS Data Driven Approach
Using GIS-based analysis, it was found that the green area of the campus represents 44% of the total study area. Further, eliminating areas with slopes exceeding 10% reduced the potential infiltration area to 36% of the green area (or 16% of the total case study area) (cf. Figure 6a). The flow accumulation was also considered in the analysis, identifying the potential infiltration areas that are located at the end of the drainage lines for the case study areas. This left 20% of the potential infiltration area (or 3.2% of the total case study area) (cf. Figure 6b). The location of these 3.2% of the case study area are found in drainage areas 6,8,9,10,11,12,13,20, and the lower parts of area 14 and 18 (cf. Figure 5).

SUDS Implementation Potential
Looking at the existing development of the campus case study area, the pervious area makes up 44% of the total area. This is area that is currently available for SUDS implementation. Choosing three different strategies for the implementation of SUDS in this case demonstrates the importance of coupling field measurements with GIS based analysis.

GIS Data Driven Approach
Using GIS-based analysis, it was found that the green area of the campus represents 44% of the total study area. Further, eliminating areas with slopes exceeding 10% reduced the potential infiltration area to 36% of the green area (or 16% of the total case study area) (cf. Figure 6a). The flow accumulation was also considered in the analysis, identifying the potential infiltration areas that are located at the end of the drainage lines for the case study areas. This left 20% of the potential infiltration area (or 3.2% of the total case study area) (cf. Figure 6b). The location of these 3.2% of the case study area are found in drainage areas 6,8,9,10,11,12,13,20, and the lower parts of area 14 and 18 (cf. Figure 5). Figure 6. GIS analysis of the NTNU campus case study with highlighted the potential SUDS areas (a) and highlighted the high potential SUDS areas (b) next to main flow paths.

Field Measurement Based Approach
Using the field measurement results and the recommendation of 10 cm/h infiltration rates in bioretention areas in cold climates given by Paus et al. [18], the best suited areas that are based on infiltration rates can be identified. This yields five drainage areas; 6, 7, 10, 12, and 16 (cf. Figure 5). The total area draining into these areas include 176,327 m 2 from the areas around.

Field Measurement Based Approach
Using the field measurement results and the recommendation of 10 cm/h infiltration rates in bioretention areas in cold climates given by Paus et al. [18], the best suited areas that are based on infiltration rates can be identified. This yields five drainage areas; 6, 7, 10, 12, and 16 (cf. Figure 5). The total area draining into these areas include 176,327 m 2 from the areas around.

Combined GIS Analysis and Field Measurement Approach
If we couple the results from the GIS analysis with the field measurement approach, the best fitted area for SUDS given both drainage lines and infiltration potential can be identified. This can further be used to estimate the volume that can be managed in the SUDS based on a requirement of, for example, 90% of the annual volume, or a specific return period volume. The best fitted areas by overlaying the two approaches are the areas identified in drainage area 10, 12, and 6. The south end of the case study area, drainage area 8, 9 though identified through the GIS analysis as well suited for SUDS, show poor infiltration capacity, well below the recommended minimum 3 cm/h rate. The same is true for the north west end of the case study area, area 20. Area 10 and 12 are part of the green area in front of the existing main building on campus. The buildings around these areas have high historic value limiting the options for exterior retrofitting of downspouts. However, there are also 22 544 m 2 of roads and parking spaces. The runoff from these areas can with minor retrofitting be drained into SUDS, like bioretention cells in area 12 and 10. In area 6, there are several existing buildings with flat roofs; however, though the buildings are not of historic value the potential for retrofitting by disconnecting the roof drains is limited. The cold climate in the case study areas means that all flat roofed buildings have interior roof drains, and leading water over the parapet is generally avoided due to freezing problems (Norwegian National Building code, Tek17, www.dibk.no). However, retrofitting with green roofs could be a good option here. This will also reduce the area needed for bioretention to 13.3 mm given the 6 mm retention in the green roofs.
Looking at the required areas for SUDS for the case study area using bioretention areas as an example and using the design equation that is based on Darcy's Equation (4), it is possible to evaluate the total area needed for bioretention if the whole case study area was retrofitted with bioretention areas to manage all runoff from the impervious surfaces (Table 1). The importance of these results is the large variation in the required area as a result of the K sat values from the various field measurements. This shows the potential in improving SUDS performance by localizing the best placements early in the planning and design process. Significant reduction in required area and cost savings could be achieved by this method.

Discussion
The two first methods applied to analyze the placement of SUDS; the GIS analysis, and the field measurements approach gave different results, which is not necessarily surprising. Applying a combination of both methods, a partial overlapping area can be identified as the third method. Demonstrated through the case study, it can be seen that the third method yields a result that is superior to either of the two by identifying the areas where criteria from both the GIS analysis and the field measurements are met. Following strictly the GIS analysis of drainage lines and slopes less than 10% yields predominantly areas at the bottom of the steep hills that drop down from the center plateau of the study areas. Though this would serve to collect the maximum amount of runoff in an efficient manner, the high runoff volumes that are generated from the impervious areas on the plateau that would be transported down these steep hills would yield high flow velocities, which could lead to erosion problems. In addition, when the field measurement results are added, it can be seen that the areas at the bottom of the hills; areas 19, 20, 8, and 9 are also the areas with the poorest infiltration rates, where infiltration based SUDS are not recommended without replacing the in-situ soils. Infiltrating more of the runoff directly on the plateau where the infiltration rates are higher would decrease the runoff volumes down the steep hills.
The variation in infiltration rates are typical for urban areas, and they are largely due to soil compaction, as shown in several previous studies, among others [27]. This variation demonstrated the need for a combined GIS and field measurement approach for SUDS placement. The NTNU campus will add a significant number of buildings to the case study area over the next decade. Applying the demonstrated combined method of an in-depth GIS analyses coupled with the field measurements approach can optimize the locations of building with respect to stormwater management early in the planning process. This will ensure a sustainable stormwater management with SUDS placed at locations with favorable performance.

Conclusions
This case study investigated how coupling GIS analysis that is based on standard landuse maps and DEMs with targeted field measurements of infiltration potentials can significantly increase the capacity of SUDS to manage stormwater. The method that is proposed by coupling drainage line analysis with field infiltration measurements demonstrated how the area that is required for SUDS could be significantly reduced by considering the infiltration potential in addition to the GIS analysis.
The field measurement from the case study area of 50 ha can vary greatly from practically zero infiltration to relatively high rates. When comparing the K sat values with the soil physical analysis, it was shown that physical characteristics is not a good indication of infiltration potential in urban soils with a large degree of compactions.
The site specific placement was shown to have a large impact on the area required for SUDS and the potential for using local in situ soils. The area requirement for handling 90% of the annual rainfall volume (19.3 mm) was up to 70% between two neighboring drainage areas (10 and 12). This is important, as it has cost implications, as for development cases there is always a pressure to minimize area required for SUDS implementation. There is also a cost incentive to use in-situ soil to the greatest extent possible, where the local infiltration potential become an important parameter. However amending the in-situ soil with small amounts of sand of compost or digging it up and disaggregate it to make it less compacted will also go a long way, especially in the case study that is presented here where the sand content was high.
Supplementary Materials: The following are available online at http://www.mdpi.com/2071-1050/10/12/4683/ s1, Table S1: Soil sieving analysis of the sampling sites. Acknowledgments: This study was initiated and supported by Klima 2050 and DRENSSTEN with the aim to reduce the societal risk caused by stormwater floods associated with climate change.

Conflicts of Interest:
The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript, or in the decision to publish the results.