Estimation of Wastewater Discharges by Means of OpenStreetMap Data

: For the optimization of sewer networks and integration of water management in urban planning, estimations of wastewater discharges at a high spatial resolution are a key boundary condition. In many cases, these data are not available or, for reasons of data protection and company secrecy, the data are not accessible for research purposes. Therefore, procedures are needed to determine the volume of wastewater with high spatial resolution, based on freely accessible data. The approach presented here uses mainly OpenStreetMap (OSM) data, combined with a dataset of the German official topographic–cartographic Information System (ATKIS), to estimate the volume of wastewater on a building level. By comparison with daily values of the dry weather inflow at pumping stations and sewage treatment plants, it is shown that the method can generate realistic results, if target inflows exceed 50 m³/d. Difficulties due to the effect of commuting and the individual use of the buildings have to be considered, as well as data ‐ quality issues in the OSM dataset. As an application example, the generated wastewater discharges are spatially joined with land ‐ use plans. The resulting wastewater yield factors serve as input data for decision ‐ support tools in urban water planning or modeling tasks.


Introduction
Sewage networks for the transport of wastewater from industrial, commercial, and residential areas to sewage treatment plants are essential elements of the water infrastructure. As the construction and maintenance of wastewater networks are cost-intensive, a large number of research projects are dedicated to the optimization of such infrastructures [1].
In recent years, the focus has been on the development of optimization algorithms for layout and component size of sewer networks [1][2][3][4]. Wastewater discharges as a main boundary condition for optimization are often regarded as "given" [1,3] without going into details in terms of data sources. In other works, the amount of wastewater is merely estimated, e.g., based on the drinkingwater consumption of larger supply areas and provided with a peak factor [2]. Here, wrong assumptions about the amount of wastewater may lead to an oversized or undersized system [5].
In the work of Willuweit and O'Sullivan [6] a combination of models is used to simulate water demand, water supply, wastewater, and runoff in urban areas, under changing land-use and climate scenarios. In the water balance approach, it is assumed that the amount of wastewater produced corresponds to the domestic water consumption. To estimate the water demand, flows monitored in district metered areas are evenly distributed to cells of 4 ha size, using GIS software.
To determine the amount of wastewater inflow into existing drainage systems, the worksheet 118 [7] of the German Association for Water, Wastewater and Waste (DWA) recommends measurements of dry weather runoff with long measurement periods, if possible, in different seasons, as well as additional measurements in commercial, industrial, or tourist areas. However, the collection of data with a temporally and spatially high-resolution involves a great deal of effort. In addition, data protection and company secrecy restrict the availability of such datasets. Hence, procedures are needed to determine the volume of wastewater with high spatial resolution, based on freely accessible data.
A common approach to determine the domestic water use is by population figures. For instance, Schiller and Bräuer [8] present a method using population figures from official statistics, disaggregated by building footprints from ATKIS (official German topographic-cartographic dataset). The main part of their method is the classification of three settlement/municipality types and ten building types based on topographic geodata and building footprints from the ATKIS basic digital elevation model (DEM) with GIS tools. Building footprints and information about the use of buildings can also be derived from OpenStreetMap (OSM) datasets. Bakillah et al. [9] use points of interest (POI) from OSM as indicators for high or low population densities, to disaggregate population figures of the city of Hamburg into a grid (500 m² per grid cell). They then distribute the population of each grid cell proportionally to the OSM building footprint size within it. Kunze [10] uses building attributes stored in OSM data to describe the share of non-residential use in existing buildings. Fan et al. [11] derive different building types from semantic information and the shape of the building floor area in OSM.
OSM geodata are maintained and continuously expanded by the OSM community. This way of generating and collecting data is called volunteered geographical information (VGI) [12]. Since the data are collected by nonprofessionals, the preparation and quality assurance of VGI data, in general, and of OSM data, in particular, are much more demanding than of data collected specifically for a particular research question. A large number of articles is devoted to quality assessment and assurance at VGI in general [13][14][15][16][17][18][19]. Particularly for OSM building datasets, the data quality is the subject of various studies regarding completeness [20,21], building density [22], and accuracy of geometries [23]. Heterogeneity and incompleteness of the data is a common challenge. A higher data density and accuracy of OSM buildings in cities compared to rural regions can be explained by the number of local participants in the OSM project, which is generally larger in cities [15,16,24]. The same findings are noted for land-use data [16]. In addition, incorrect data collection by nonprofessionals is a common difficulty, as is the lack of metadata or the ability to verify data [17].
Nevertheless, VGI has proved to be a useful data source in many fields of application and scientific questions, since-as the word "volunteered" suggests-they are voluntarily made available to the public and, thus, the user does not have to spend any time on data collection. OSM data are freely available and may be used under the Open Database License [25]. Application examples for VGI are the mapping of flood hazards [26], disaster management [27], modeling of energy infrastructure [25,28], validation of land-use maps [29], or ecological monitoring [30]. The application of OSM data in urban-planning contexts has so far been investigated primarily in the field of transport infrastructure planning [31,32] and the analysis of settlement structures [8,33]. In the field of urban water management, the OSM road network was used to generate virtual wastewater networks [34,35]. The connecting element between settlement structure and sewage network was not considered in the mentioned work: wastewater from residential buildings, industry and commerce, or other uses that is transported to a treatment plant along the sewage network.
In this article, a method for the estimation of the wastewater volumes based on OSM data is presented. The parameters used for the calculations are optimized by comparing the estimated wastewater volumes in the catchment area of wastewater treatment plants (WWTPs) and sewage pumping stations (SPSs) with measured dry weather inflows. As a sample application, estimated wastewater discharges are spatially intersected with land-use plans, to develop a scenario-capable tool for the integral planning of settlement structural measures.

Study Area
The "Regiopolregion Rostock" is located in the northeast of Germany (see Figure 1, left). A regiopol region consists of a smaller large city with the surrounding rural region. Urban and rural areas differ in building types, infrastructure, land use, and population density [36], leading to a large heterogeneity concerning amount and variability of mass flows, such as wastewater [37]. The socioeconomic center of the area is the city of Rostock, with approximately 200,000 inhabitants. The communities in the surrounding area, by contrast, have a more village-like or small-town character. The area of responsibility of the Warnow Water and Wastewater Association (WWAV), purple boundary line in Figure 1) covers the entire urban area of Rostock with a sewer connection rate of 99.7%, as well as 29 municipalities in the south and east of the city. Wastewater treatment is mainly performed at the central wastewater treatment plant in Rostock (point 1 in Figure 1  For the method presented here, building polygons and land-use polygons are obtained from the OSM provider Geofabrik (for Mecklenburg-Vorpommern: [39]). The dataset clipped to the study area contains 79,053 building polygons. Basic elements in OSM are points/"nodes" (points of interest or centroids of surfaces), open or closed lines/"ways", and "relations" (relationships between objects, such as groups of objects). All OSM elements get their properties from attributes/"tags", each consisting of a key and a value [40]. Polygons are constructed by closed lines with the attribute "area=yes" and with certain tags like "building=*" or "landuse=*" ("*" being a particular attribute), which are assigned exclusively to areal elements. In OSM, buildings, in particular, are assigned by the tag "building=*". For example, a building is already clearly identified by "building=yes", but it can also be more precisely labeled by other values (e.g., "building=hotel"). However, an element declared as a building could also get the property "hotel" by "tourism=hotel". For the assignment of OSM tags, the OSM Wiki website provides a multitude of hints and examples. These can be regarded as a guideline more than a fixed set of rules. New tags can be proposed, changed, and discarded by the OSM community [40].
Additional information on land use, as well as population figures of municipalities and municipal boundaries, is taken from official data (e.g., ALKIS, provided by regional administrations [41]). A digital map (shapefile) of the sewage network and measured values of the daily inflow to WWTPs and SPSs from the regional wastewater association WWAV are used for the optimization process.

Generation of Wastewater Volume per Building
GIS software is used to link OSM building polygons with information on land use (OSM: land use, ALKIS) and location (municipality, WWTP catchment area) by spatial intersection. The further calculation of the wastewater volume is performed in four steps using R scripts [42]. An overview about the method is illustrated in Figure 2. According to the wastewater discharge components described by the German DWA worksheet 118 [7], buildings are classified into "residential buildings", "industrial buildings", and "commercial buildings" in step 1. Buildings that combine different uses (e.g., retail trade with living space on the floors above), are assigned to the type "mixed-use buildings". Buildings that (presumably) do not discharge wastewater into the network, are assigned to the ʺNULLʺ category. The classification process is carried out by various attributes: 1. OSM tags of the buildings, e.g., "building", "leisure", and "amenity" (see [11]). Some examples are shown in Table 1. 3. If none of the former is stated: land-use information derived from ALKIS.
In step 2, the number of floors is-if available-directly extracted from OSM tags ("levels"), derived from the building height ("building:height") or estimated on the basis of the usage class of the building. The usable area of a building is calculated by multiplying the floor area with the number of floors, less 5% as unusable area. The number of inhabitants (NI) of every municipality is then disaggregated proportionally to the usable area of each of its residential or mixed-use buildings (step 3).
Step 4 of the algorithm is the calculation of wastewater discharges per building by Equation (1) (residential buildings), Equation (2)
Discharge rates for residential (qR) (in L/(person*d)), commercial (qC), and industrial buildings (qI) (each in L/(m²*d), related to the building footprint area) are determined by the optimization algorithm described in the following section. It is assumed that a large proportion of the inhabitants of the surrounding municipalities does not work at their place of residence. As a result, a part of the domestic wastewater is not produced in these municipalities, but is already contained in the wastewater from commercial buildings. This fact is taken into account in the calculation by a "rural factor", fR (a fixed value for every municipality; between 0.6 and 1), depending on the population density of the municipality.

Optimization of Parameters
The wastewater volume, Qcatchm,n which is generated in the catchment area, An, of a sewage treatment plant or a pumping station with the number, n, consists of domestic, industrial, and commercial wastewater: Qcatchm,n = ∑QR,n + ∑QI,n + ∑QC,n, Catchment areas of pumping stations and sewage treatment plants are defined along the sewage network in the area of responsibility of the WWAV (see Figure 1). For this purpose, buildings are linked to the nearest sewer by GIS queries. Flow paths are determined by using routing scripts from the QGIS plugin "WaterNetAnalyzer", which was developed by the authors within this study [43].
The amount of wastewater, Qcatchm, generated in the catchment area of 15 WWTPs (Points 1 and 7-20 in Figure 1) and 5 main SPS (Points 2-6 in Figure 1) are compared with corresponding measured values of the daily inflow. The first quartile of dry weather inflows (1 January 2017 to 31 December 2018) is selected as the target value for each of the five main pumping stations and the WWTPs in the study area. Dry weather inflow is assumed, if the total precipitation measured on the day itself and on the previous day is below 0.1 mm/d. Where daily values of the inflow are not available, the average annual inflow of wastewater (2013-2017) to small WWTPs converted to a daily value is used as a target. The parameters qR, qI, and qC are optimized by minimizing the deviation of the generated wastewater quantities from the target value for (Equations (5) and (6)): Replacing Qcatchm,n in Equation (5) with Equations (1)-(4) equals the following: min ∑ (Qtarget,n − (∑NIcatchm,n × qR × fR + ∑AI,usable,catchm,n × qI + ∑AC,usable,catchm,n × qC)) 2 (6) Optimum parameter combinations are determined using the R function "optim()" with the Nelder-Mead method [44], and the R function "lm()" for linear regression. The values used for optimization are presented in detail in Table A (Appendix A).

Application Example: Aggregation at LUP Level
Land-use plans (LUP) represent the intended urban development of a municipality. They are an essential tool for the integral planning of structural settlement measures in regiopol regions. A catchment yield factor, qLUP,i (in L/(s*ha)), is derived for each subarea, i, of the LUP by spatially joining the estimated wastewater volume in GIS with the LUP of the City of Rostock (see Figure C, Appendix C). For this purpose, wastewater discharges are aggregated in every LUP subarea and divided by its area ALUP (Equation (6)).
The range of the catchment yield factors sorted by types of land use (e.g., residential area, mixeduse area, etc.) can be used in planning scenarios. For example, the effect of a planned industrial park on a sewage network can be investigated in models by multiplying its intended area with the median catchment yield factor of industrial areas.

Classification of Buildings
The classification of the buildings into the mentioned categories is based on OSM tags and information on real land use from ALKIS data. Figure 3 shows the result of the classification in the inner-city area of the city of Rostock. In the catchment of WWTP Rostock, the share of residential space is 67%. Moreover, 25% of the building areas are identified as commercial, and 8% are industrial areas. However, the distribution of the area categories also varies within the urban area (cf. Figure 4). In the catchment area of the SPS at point 6, the share of residential areas is only 35%. In contrast, 75% of the building areas are categorized as residential in the catchment of the SPS at point 4, where larger multistory blocks of flats are located. In the rather rural part of the study area, "residential" accounts for an average of 82% of the building areas, 11% are commercial and 7% are industrial building areas. Absolute values for each area type are provided in Table A (Appendix A). An essential element of uncertainty is the unknown completeness and correctness of the OSM building polygons. Götz and Zipf [45] estimated that, at the time of their investigation only approximately 30% of all building geometries in Germany are recorded in OSM. They expected that continuous contributions of the OSM community would lead to an increasing percentage here, up to 90% within a few years. Even without a more precise evaluation, it is assumed that this figure will not yet be reached in the study area, as in certain villages, hardly any buildings have yet been marked (cf. Figure 5). A large number of missing building geometries can be expected in the rural municipalities [20]. A comparison of OSM buildings with official cadastral data can help to identify missing geometries here [21]. Irrespective of this, the OSM dataset for the city of Rostock is very extensive, as the Cadastral Office of Rostock donated a building dataset to the OSM project in 2009, including information on building height and use [20]. Moreover, the number of buildings in OSM is steadily increasing. Compared to an OSM dataset of the same research area obtained two months before, the dataset used here contains already 7695 more (which corresponds to +4.7%) building polygons. Not only completeness (in terms of the number of buildings) but also the number of tags assigned in rural areas in the study area is considerably lower. Properties such as the building height or the number of floors are mainly available in the urban part of the study area. About 82% of the buildings are only marked by the tag "building = yes" [47]. To deal with this lack of information about building attributes, the layer can be intersected with land-use data [25]. In the current version of the method presented here, the attribute "land use" from the ALKIS dataset is also utilized to categorize buildings. In cases where even this information is not available, the category "residential buildings" is chosen, which applies to 0.6% of the buildings, mostly in rural areas. In this way, the algorithm classifies all buildings reliably into wastewater categories. At the same time, however, agricultural buildings in rural areas-stables, halls, etc.-are wrongly identified as residential buildings, when OSM tags are missing. As a result, the disaggregation of the population per municipality distributes the inhabitants among too many houses. The classification process could be improved by similarity measurements of building footprints, e.g., as described by Fan et al. [11], or by interpreting buildings and related "points of interest" [29].

Determination of the Inflow Target Value
A critical point of the method is to find suitable target values for the calibration of the discharge rates qR, qI, and qC. In the approach presented here, daily measurements of wastewater inflow to treatment plants and pumping stations are used. Since the volume of domestic and commercial wastewater is subject to daily fluctuations, the absolute minimum cannot be the target value. In this study, dry weather inflow is assumed if rainfall on the day itself and on the previous day is below 0.1 mm/d. Figure 6 shows the inflow to a pumping station at such conditions. In the periods from July 2017 to September 2017 and October 2017 to February 2018, high precipitation is recorded, which is the cause of higher water inflows to sewage treatment plants and pumping stations (see Figure 6). In general, the year 2017 can be classified as high in precipitation with 739 mm compared to the longterm average of 616 mm (1981-2010, calculated by precipitation data from the German national meteorological service (DWD) [48]). Presumably, due to the finely branched Rostock sewage network, which extends from the urban area to the suburbs and surrounding villages, rainfall events that occurred more than two days before are wrongly detected as dry weather inflow. Additionally, unavoidable infiltration of stormwater into separate drainage systems leads to higher inflows. Groundwater infiltration is likely to occur, as groundwater levels are generally high in many parts of the study area. Another issue that has to be considered here is extraneous water and snowmelt infiltration into the sewer system [49]. This is the case in March 2018 to April 2018. The mean value of the dry weather runoff is increased by such events and is therefore unsuitable as a target inflow value. To take the mentioned conditions into account, the first quartile of the dry weather inflow at sewage treatment plants and pumping stations is chosen as the target value. Other target values and approaches to determine them may be more suitable for other study areas. For example, wastewater volumes can be compared to drinking water intake in a first step. This can be used to exclude inflow values influenced by rainfall infiltration. Alternatively, hydrograph separation could be used to distinguish between wastewater flows and stormwater flows [50]. Where seasonal changes in wastewater inflow are high, the selection of separate target values for optimization for each season might be beneficial. To avoid misinterpretation of seasonal changes caused by single extreme events, the evaluation of time series longer than two years is recommended.

Parameter Optimizing Process
The two optimization methods result in similar parameter combinations, both for the optimization on the basis of the wastewater treatment plant inflows and the pumping station inflows. Since the catchment areas of the pumping stations are mostly located in the urban area and those of the sewage treatment plants (with the exception of the central sewage treatment plant Rostock) cover the rather rural area, differences in the optimized parameters qR, qI, and qC (cf. Table 2) can be attributed to settlement structural conditions. The estimated amount of wastewater, qR, per inhabitant for parameter optimization is between 61.5 and 93.6 L/d, depending on the optimization method and the target values considered. The lower values are obtained for rural areas. As residential use is the main source of wastewater in the study area and the per capita volume of wastewater of 93 L/d matches the average drinking water consumption of 93.8 L/(person*d) indicated by the WWAV in 2017 [38], the parameter set optimized for the entire study area is chosen for further calculations:  qI = 0.6 L/(m²*d).
The optimized parameter, qC, for commercial buildings varies between 2.4 and 3.1 L/(m²*d). For industrial buildings, the optimization result of the parameter qI is 0.6 to 2.9 L/(m²*d). Both qC and qI are higher in the rural part of the study area. These findings are explained by the classification process (cf. Figure 2, step 1). The algorithm identifies large production halls with presumably low wastewater output as industrial buildings in the urban area of Rostock, while agricultural buildings with presumably higher wastewater output (e.g., stables) are classified as "industrial" in rural areas. Accordingly, the actual amount of wastewater per m² in rural industrial buildings would also be higher than in urban industrial buildings (e.g., production halls), which is reflected in a higher optimized qI. Although a part of the drinking water in rural areas is not returned to the sewage network (e.g., for the irrigation of plants in gardens and dwellings) 61.5 L/(person*d) is too low, compared to the average consumption. There are several possible explanations for the low qR value in rural areas. The effect of missing or incorrectly assigned buildings in the OSM dataset, as well as the effect of commuting, has to be considered. A large proportion of the rural population's wastewater may be produced at the workplace in urban areas, so an adjustment of fR is needed. In addition, with this factor, the settlement structure of different localities within a municipality is not taken into account so far. A small village with residential development next to a larger village with business parks (and accordingly jobs) in the same municipally is therefore evaluated equally. Here, further work is needed to establish parameters which can better represent this heterogeneity.

Calculated Wastewater Volume
The comparison of the calculated wastewater inflow with the daily values of dry weather inflow at SPSs and WWTPs shows that the method can generate realistic values mainly in urban areas (cf. Figure 7a, numbers 1-6). Here, the maximum relative error is 31%. The estimated wastewater volume of the central WWTP of Rostock (number 1) is 35,992 m³/d, which differs by 0.1% from the target value of 35,947 m³/d. Moreover, WWTPs in rural areas show relative errors less than 30% if target inflows exceed 50 m³/d. The estimated volume at WWTP number 8 is 341 m³/d, with 391 m³/d set as the expected value. WWTPs with a target inflow smaller than 50 m³/d can have large relative errors. For example, the estimated quantity of 14 m³/d at sewage treatment plant number 17 is more than twice the estimated value of 6 m³/d. An overview about all values is provided in Table A2 (see  Appendix B).
It must be borne in mind that the target value (1st quartile of dry weather inflow) is a parameter defined ad hoc which represents an assumption about the true average inflow in dry weather. In addition, the calculation with a single set of mean discharge rates qR, qC, and qI cannot reflect the individual use of a specific building. For example, individual behavior of water consumption will lead to overestimations or underestimations if catchments are small. Furthermore, in sewage treatment plants with a smaller catchment area in rural areas, only a few missing or incorrectly allocated buildings have a negative effect on the calculated result. Hence, the relative error related to the target value is smaller, the larger the catchment area of the sewage treatment plant or pumping station (cf. Figure 7b Using disaggregation procedures as an example, Schiller and Bräuer [8] also argue that the individual error is greater the higher the level of detail. In the case of values generated on a building level, the error range is largest. The calculated result should therefore be aggregated on the largest possible reference level according to the planning task. For the drainage planning, the mentioned study proposes the level of districts or "solitary settlement units". In the application example described here, the spatial resolution of LUP is chosen, which corresponds the proposed aggregation level by Schiller and Bräuer [8].

Application Examples LUP and Urban Wastewater Management
As described in Section 2.4, yield factors for every type of land use in LUPs are derived by a spatial join and aggregation process, using QGIS. The combination of land-use types and wastewater discharges enables illustrating the effect of land-use change in modeling future scenarios.
The largest wastewater yield factors ( Figure 8) are related to the LUP types "areas for public facilities" and "core development areas" and (medians: 0.14 to 0.22 L/(s*ha)). Industrial parks, mixed building areas, residential areas, and special building areas are in the middle range of yield factors (medians: 0.05 to 0.07 L/(s*ha)). In the work of Willuweit and O'Sullivan [6] , three land-use types are distinguished in the urban area of the Dublin region: For "medium dense to continuous dense urban fabric/commercial", a value of 0.48 to 0.50 L/(s*ha) is estimated. For the category "sparse and discontinuous urban fabric/residential", 0.14 to 0.16 L/(s*ha) is calculated, and 0.10 to 0.14 L/(s*ha) for cells of the land-use type "industrial". Compared to the region of Rostock, the higher values cited here presumably result from a higher per capita water use. The water demand of the population of approximately 1.5 million in the Dublin Region Water Supply Area is 550*10 6 L/d, which corresponds to 367 L/(person*d).
The DWA worksheet A 118 [7] recommends using 0.2 to 0.5 L/(s*ha) for choosing the component size of sewer networks in commercial and industrial areas for companies with low water consumption and up to 1.0 L/(s*ha) for companies with medium to high water consumption, referring to the catchment area connected to the sewage system. These values are higher than the calculated ones (see Figure 8), as they represent an (hourly) maximum instead of daily mean values. However, deriving the variability of wastewater discharge from OSM data or other VGI could be a subject of future research.
In Figure 8, wastewater is also "produced" on road areas, agricultural areas, park and forest areas, and areas for waste management. This can be due to incorrectly classified buildings or inaccurately drawn LUP areas. On the other hand, e.g., road areas actually contribute to the volume of wastewater, as tram depots or ticket counters are located in these areas. The GIS-based approach to accumulate generated wastewater along a sewer network allows the estimation of the actual flows (see Figure 9). With the same approach, wastewater flows can be represented in future urban and regional planning variants. Possible applications are decision support systems for the construction or reconstruction of drainage systems on a conceptual level.
In the rural structured part of the study area, phosphate inputs from wastewater treatment plants into rivers are investigated [51]. When merging decentralized small wastewater systems into a central system, it is worth considering a fourth treatment stage to reduce the input via the treated wastewater. In order to determine the technological optimization potential, the amount of wastewater to be treated is the central data basis, in addition to the sewer system itself and equipment of the wastewater treatment plants.
On the other hand, especially in developing countries, the construction and maintenance of large sewerage networks and cost-effective wastewater treatment are a challenge. With the help of spatially differentiated wastewater data, decentralized wastewater concepts can be created in addition to existing systems, which enable the wastewater treatment and reuse close to the source [52].

Conclusions
The method introduced here represents an instrument to estimate discharges of wastewater from domestic, commercial, and industrial use on a high spatial resolution. The approach adds a new application for VGI to the existing methods. Using OSM as the main data source makes the method transferable to other regions. Where quality issues of OSM building data generate uncertainty, especially in rural regions, additional data from different sources, e.g., authoritative data, may be used to identify missing buildings. To achieve higher accuracy in the classification process, OSM buildings data could be enhanced with POI from OSM, or an analysis of shape and size of building footprints could be applied.
For the calibration of the discharge rates qR, qI, and qC, finding suitable target values is an essential step. Extraneous water, such as groundwater infiltration through leakages and rainwater or snowmelt infiltration into the sewage system, leads to temporarily high discharge rates, which increase the difficulty of finding a target value for the inflow. As this is a crucial step of the method, different approaches for the determination should be tested in future research, e.g., hydrograph separation or calibration with consumption rates of drinking water. Where long-term series of inflow to WWTPs and SPSs are available, the variability of wastewater discharges can be estimated by minimum and maximum values for different spaces of time. Seasonally varying target values derived in this way could improve the accuracy of the results, especially for small catchment areas of WWTPs. The method shown here, which uses one target value per WWTP or SPS, can be applied when inflow data sources are limited.
Wastewater discharges aggregated along the sewer networks in the study area generate realistic inflows to WWTPs and SPSs, which differ by a maximum of 31% from the target inflow when this value exceeds 50 m³/d. In small catchments, individual behavior of water consumption may lead to overestimates or underestimations. Here, the integration of peak factors to represent variability of wastewater discharges could improve the estimation. The estimated values cannot replace timeconsuming measurements. However, the method can fill the gap of input data for model-building and urban water management, especially for research purposes, when measured data do not exist or their availability is legally restricted. By intersection of discharges with land-use plans, wastewater yield factors are derived as a scenario-capable tool, which enables us to simulate the effect of landuse change in modeling future scenarios. Funding: This study was conducted within the framework of the Project PROSPER-RO, funded by BMBF, grant number 033L212. We acknowledge financial support by Deutsche Forschungsgemeinschaft and Universität Rostock within the funding program Open Access Publishing.