How Many Small Agglomerations Do Exist in the European Union, and How Should We Treat Their Wastewater?

The European Union (EU)’s legislation on urban wastewater requires all agglomerations with a population equivalent (PE) above 2000 people to undergo a secondary (mechanical/physical and biological) wastewater treatment. Agglomerations below 2000 PE, though, fall outside the scope of the current EU’s legislation. As such, their regulation is heterogeneous across the various EU member states, and there is no systematic collection nor reporting of data enabling an estimation of their actual significance as a source of pollution for the receiving water bodies. Here we present a spatial model to delineate agglomerations in a GIS, based on population distribution and land cover. From the model results, in the EU, we identify 364,650 agglomerations with 2000 PE or less, housing a cumulative population of about 75 million inhabitants. We then calculate the organic matter and nutrient loads these agglomerations can discharge, assuming they presently undergo primary wastewater treatment, and the reduction of loads that can be expected under different treatment scenarios, together with the corresponding treatment costs based on a simple cost model. Using a conventional shadow price for the organic matter and nutrients removed, we show that all treatment scenarios show a benefit-to-cost ratio (B/C) above (or close to) 1. However, only a scenario of secondary treatment applied to all agglomerations above 1000 PE provides sufficient safety margins on the B/C. This suggests the opportunity to expand the scope of the current legislation down to agglomerations of this size, while addressing smaller agglomerations depending on their actual impacts on the receiving water bodies, through “appropriate treatments” defined by the local authorities.

with 2000 population equivalents (PE) or more. Smaller agglomerations (less than 2000 PE) are not regulated, although Art. 7 of the UWWTD requires member states to implement an appropriate wastewater treatment for agglomerations below 2000 PE served by a sewer network. In all other cases, member states are not due to act nor report on the matter in the absence of well-identified impacts on the receiving water bodies. Because of this, we have limited evidence in order to quantify the pollution coming from small agglomerations at European scale. In this contribution, we present an estimation of the population living in agglomerations smaller than 2000 PE, the associated loads of pollutants and the costs of treatment under scenarios where we extend the obligations of the UWWTD to agglomerations of less than 2000 PE. Based on the results, we propose recommendations for the regulation of small agglomerations in the EU.

Delineation of Small Agglomerations
As a first step in the analysis, we generate a distribution of agglomerations statistically consistent with the agglomerations delineated by the EU member states compliant with the UWWTD. To this end, we examine the spatial distribution of population in Europe and we group it into agglomerations through a mathematical morphology algorithm. We refer to the 100-m population density map of Europe of Freire et al., (2016). We denote this map as Pop, and we process it through the following steps: 1. We define a threshold population density 0 , above which we assume collection of wastewater is justified, and generate a Boolean map 2. We expand the zones where C 0 = 1 by a number of cells n. We call the resulting map C1. 3. We apply a clumping (or region-grouping) algorithm to map C1 (Pistocchi, 2014). We consider the regions grouped with this operation as agglomerations, and we assign each of them a univocal identifier. 4. For each of these simulated agglomerations, we can compute the total population from map Pop, using a zonal statistics algorithm. All the above calculations are performed using ESRI ArcGIS 10.7 Spatial Analyst software.
The steps of the procedure are illustrated in Fig. 1. The procedure yields different results depending on the parameters n and 0 . These were calibrated in order to reproduce as close as possible the reported number of agglomerations and cumulative population by agglomeration size at EU level.
The results of the calibration are provided in Annex 4.. The values finally chosen were n = 2, 0 = 10 persons/ha. This implies that the smallest possible agglomeration is composed of 10 persons and has an extent of 1 ha.

Costs of Wastewater Collection and Treatment
For the sake of a quantification at European scale, we assume a bottom line scenario where all households in agglomerations below 2000 PE are served by a primary treatment, most often consisting of a single household septic tank, with no collection network. Under this assumption, treating wastewater of small agglomerations to the level required by the UWWTD for agglomerations above 2000 PE entails a cost, which will include the collection network and the treatment plant. The cost will generally depend on location-specific factors such as the extent of the network and the design of the plant. As a first approximation, we estimate the additional cost of collecting wastewater in a small agglomeration and treating it to a secondary level (compliant with art. 4 of the UWWTD) as the difference of the cost of the network combined with the cost of a secondary plant and the cost of primary treatment. For the cost of the network and the primary treatment (assumed to coincide with a septic tank), we refer to OECD's FEASIBLE model expenditure functions (COWI, 2010;OECD, 2004) while for the secondary treatment of small agglomerations, we refer to the average of cost functions for various types of treatment wetlands (TWs) described in Pistocchi et al. (2020). TWs are common solutions for small plants, together with other technologies including, e.g. sequential batch reactors (SBR) and other small technical aerobic systems (see, e.g. Langergraber et al., 2018). TWs can be designed to achieve performances at least comparable to biological treatment processes for larger agglomerations (Dotro et al., 2017;Langergraber et al., 2019).
Eventually, the additional combined costs of collection and secondary treatment are given by the average of costs functions of TWs, plus the cost function of the collection network, minus the cost function of the septic tank. The combination is described by the following cost function: where C Cremoval is the annualised cost (Euro/PE) and Pop the agglomeration's population. All details of the derivation of the above equation are provided in Annex 5..
In addition to secondary treatment, we consider the additional cost of TWs in order to provide denitrification. Assuming that a hybrid TW, providing up to 80% N removal, corresponds to the highest cost C Cremoval = 218.36Pop −0.37 function among those described in Pistocchi et al. (2020), such additional cost is estimated as: where C Nremoval is the annualised cost (Euro/PE) in addition to C Cremoval , in order to have denitrification. All details of the derivation of the above equation are provided in Annex 5.. Finally, for P removal, we refer to the cost function: C Premoval being the additional cost of P removal. The above equation derives from the OECD FEASIBLE model expenditure functions as explained in Annex 5..

Cost-effectiveness Analysis of Treating Small Agglomerations
We assume the load of organic matter (as 5-day biochemical oxygen demand, BOD) associated with untreated wastewater equal to 60 g/PE/day, while those of total nitrogen (N) and total phosphorus (P) are estimated through country-specific emission factors as in Malagò and Bouraoui (2021). In order to compare scenarios, we assume a constant removal efficiency for BOD, N and P for primary, secondary and tertiary (nutrient removal) treatment processes.

Fig. 2
Comparison of modelled and reported population equivalents (PE) and number of agglomerations for different size classes above 2000 PE. Reported data are gathered from the 8th UWWTD Implementation Report. Only for Poland and Italy, and due to their unavailability, the information is based on data gathered from the 7th Implementation Report (https:// ec. europa. eu/ envir onment/ water/ water-urban waste/ imple menta tion/ imple menta tionr eports_ en. htm) Vol.: (0123456789) These values correspond to the assumptions made in the GREEN model used for EU scale assessment of nutrients (Grizzetti et al., 2021;Pistocchi et al., 2019).
While septic tanks may remove virtually no N, we assume their effluents are discharged to small ditches or ponds, or subsurface drainage before reaching the receiving water bodies, and this is sufficient to achieve some natural attenuation of nutrient concentrations.
The removal of N and P assumed for secondary treatment may require already the design of hybrid, multistage TW systems including a P trap, while the removal of N at tertiary level may require a more expensive TW of similar type.
We compute the loads of pollutants discharged by settlements below 2000 PE under scenarios of increasing stringency of treatment, and we compare them with a scenario with all settlements undergoing primary treatment.

Number and Population of Small
Agglomerations in the EU The calibrated model allows delineating agglomerations throughout Europe and an estimation of their population. Figure 1D shows the cumulative frequency distribution of agglomerations by size and the cumulative population as a function of agglomeration size. The model allows reproducing quite well the number and resident population of agglomerations larger than 2000 PE in Europe (Fig. 2). For agglomerations smaller than 2000 PE, we estimate a total population in the EU of about 75 million persons, distributed by country as shown in Table 1, in a total of 364,650 agglomerations distributed by country as shown in Table 2. The number of agglomerations by size is reasonably consistent with a benchmark of independently estimated data (Wood plc, pers.comm., 2022), as shown in Fig. 3. This indicates the simulation model may represent the number of small agglomerations in the EU within a factor 2 of accuracy in more than 75% of the cases and within a factor 10 in all cases, with the exception of Croatia. In this case, modelled agglomerations below 500 PE are a factor > 60 more than reported. We find larger errors on the number of smaller agglomerations and a clear tendency to overestimation (Fig. 3). The modelled distribution of agglomerations below 2000 PE indicates these house around 15 to 20% of the population in the EU, with some variability among member states (Fig. 4). Population in small agglomerations is a particularly small percentage of the national total in Malta, and below average in Belgium, Spain, Finland and Ireland, while it is higher than average in the Czech Republic, Slovakia, Croatia, Slovenia, Romania and Poland. Usually agglomerations between 1000 and 2000 PE account for 20 to 30% of the population in small agglomerations, while very small agglomerations (below 100 or 50 PE) usually account for 10 to 20%.

Discharges of Pollutants with Wastewater from Small Agglomerations
We could estimate the discharges of BOD, N and P through wastewater from small agglomerations, based on the emission factors described in Section 2.3 and whether they fall in areas defined as "sensitive" according to the UWWTD (Table 3). For the identification of sensitive areas, we refer to the available map discussed in Bouraoui et al. (2022), and Pistocchi et al. (submitted). In the most extreme scenarios, we assume that agglomerations above a size threshold undergo tertiary treatment for N and P in the whole territory of the EU. With these assumptions, we calculate the discharges of BOD, N and P under each scenario, as shown in Fig. 5. Obviously, discharges decrease as we lower the size threshold above which treatment is required. BOD discharges do not change significantly after adding tertiary treatment, while the reduction of N and P is substantial. For BOD, the maximum reduction of discharges is in the order of 80% under the most extreme scenarios (requirement of a secondary treatment or higher for all agglomerations above 50 PE). Limiting the requirement of secondary treatment to agglomerations above 1000 and 500 PE would reduce baseline BOD discharges of about a quarter and a half, respectively. For N and P, Fig. 6 presents the marginally avoidable discharges attainable under increasingly demanding scenarios. Requiring secondary treatment for agglomerations above 1000 PE would   reduce discharges of about 12% for N and 13% for P. Under the most extreme scenario of tertiary treatment in the whole territory for all agglomerations above 50 PE, discharges could be reduced of about 70% for N and 80% for P.

Costs and Benefits of Treatment for Smaller Agglomerations
Using the cost model described in Section 2.2, we can estimate the total costs of the various scenarios. These are presented together with the discharges of N, P and BOD in Fig. 5. While the costs increase quite mildly under scenarios entailing secondary treatment only, tertiary treatment causes a much sharper rise as shown by the graphs. We can appraise the cost-effectiveness of the various scenarios by referring to a conventional benefit-to-cost ratio (B/C) by valuing the removal of 1 kg of BOD, N and P through a shadow price of € 0.05, € 20 and € 30 respectively, in line with UNEP, 2015. A conventional benefit is calculated for each scenario by multiplying the removed quantities of BOD, N and P (Figs. 5 and 6) by the respective shadow prices. The benefit value of each scenario is then divided by the corresponding cost shown in Fig. 5. The calculated B/C for all scenarios is shown in Table 3. While B/C for a scenario of secondary treatment for all agglomerations above 1000 PE exceeds 2, this ratio decreases with more demanding scenarios and comes close to 1, or even below 1 when more advanced treatment is required for the smallest agglomerations. Fig. 6 Incremental reduction of discharges of total P and total N from small settlements, depending on the level of treatment enforced. The colour coding reflects increasingly demanding scenarios, and the grey sector indicates the share of baseline discharges that would remain under the most demanding scenario

Discussion and Conclusions
Based on the model discussed in Section 2.1, we have identified 364,650 agglomerations below 2000 PE in the EU and the corresponding population of about 75 million people. This estimate is a factor 3 higher than a previous pan-European estimate by Vigiak et al. (2020), who estimated agglomerations not covered by the UWWTD to account for about 23 million people. Our estimates are rather in line with information recently collected and extrapolated from EU member states (Wood plc, pers. comm., 2022) as shown in Fig. 3. However, data on small agglomerations are not systematically collected and reported, limiting at present the accuracy of estimations within a factor of about 2. The GIS model used for the delineation of agglomerations is likely to overestimate the number and population of small agglomerations, because it identifies as small agglomerations any cluster of population meeting the continuity and separation criteria discussed in Section 2.1. In this way, we count as small agglomerations also many small clusters at borderline distance from larger agglomerations that in reality are likely to be connected to the wastewater treatment plant of the latter. For example, data for Austria (where all agglomerations above 50 PE can be assumed to always have a WWTP) indicate a number of 1040 agglomerations between 50 and 500 PE, 135 between 500 and 1000 PE and 120 between 1000 and 2000 PE (Lenz et al., 2021), much smaller than our estimates (5111, 667 and 447 agglomerations respectively In principle, the model could be refined through a better calibration with known and reported data on the number and served population of small agglomerations. However, it is unlikely that we can achieve better results unless we use statistics for the majority of countries. At the same time, if these were effectively available, the model itself would not be needed anyway.
In spite of these limitations, the model provides a means to estimate the consistency of small agglomerations in the EU in the absence of reported data. Based on our calculations, we can conclude that small agglomerations in the EU generate significant loads of BOD and nutrients, and may represent important pressures on water bodies.
The scope of the UWWTD covers agglomerations above 2000 PE, requiring secondary treatment unless they discharge in sensitive areas, in which case nutrient removal is also required. While various member states have already regulated smaller agglomerations, in many cases, quality standards correspond to primary treatment, if not even untreated wastewater. In principle, by enforcing an appropriate level of treatment, we could reduce pollution loads from small agglomerations to a significant extent.
However, in order to decide on appropriate treatment standards, we need to compare the costs and effectiveness of investments in treatment of small agglomerations with those in the control of other sources of pollution. The benefit-to-cost ratio (B/C) of various scenarios of treatment for small agglomerations suffers from the variability of the value of BOD removal and, most importantly, N and P removal depending on the conditions of the receiving water body. Moreover, the costs of wastewater treatment for small agglomerations depend on the initial conditions of the infrastructure and local factors such as the availability of space and the possibility to automate the management of processes. Still, the simple conventional calculation presented above provides, as a first approximation, some support for decisions in this matter. Particularly, the B/C shows that investments in secondary and even more in tertiary treatment for small agglomerations may not be the most cost-effective options to improve water quality. Nevertheless, requiring secondary treatment for agglomerations above 1000 PE shows a B/C around 2; hence, lowering the current threshold of the UWWTD above which secondary treatment is required, from 2000 to 1000 PE, may be justified by a reasonable safety margin.
An alternative to tertiary treatment for small agglomerations, yielding comparable results in terms of avoiding N and P discharges to the receiving water bodies, could be a treatment configuration enabling the reuse for water and nutrients for agricultural fertilisation via irrigation ("fertigation"). This could help avoid use of mineral fertilisers to some extent, while reducing the investment and operation costs of plants. At the same time, fertigation requires an appropriate control of pathogens and micropollutants released with wastewater, which may prove expensive and undermine the feasibility of this solution.
In addition to the costs of action and the benefits from BOD, N and P removal, in principle, it is relevant to consider the potential change of greenhouse gas (GHG) emissions under the various scenarios. To this end, a GHG emission assessment was performed including different small plant typologies (see Annex 6.). The evidence available (Fig. 11) does not suggest a clear advantage of one treatment system over the other for small agglomerations: while replacing septic tanks and similar primary systems with secondary systems significantly contributes to reducing direct emissions of methane, increased direct emissions of nitrous oxide may offset the improvement. Moreover, more advanced treatments may entail larger energy consumption, which usually corresponds to higher GHG emissions. However, we can observe that some TW typologies combined with a primary settler and sludge stabilisation in a reed bed offer lower emission factors among the considered options. Thus, hybrid multistage systems need to be applied to improve N removal, triggering higher methane emissions. The off-site treatment and disposal of sludge may add to the GHG emissions and must be properly considered. Last but not least, the infrastructure for the collection of sewage entails use of concrete and other resources and is associated to relatively high GHG emissions (Morera et al., 2016). Changes in GHG emissions must be evaluated on a case-by-case basis in order to account for the specific factors that may make a solution preferable over the others. Furthermore, it must be stressed that GHG emissions are affected by a high uncertainty related to the estimation of N 2 O and CH 4 generation in passive systems such as treatment wetlands. There is also evidence that poorly managed or overloaded systems can lead to much higher CH 4 emissions than assumed in theoretical calculations.
All things considered, more stringent treatment of smaller agglomerations is not likely to change the GHG balance significantly.
Measures for smaller agglomerations, and generally a requirement of more stringent treatment for agglomerations below 2000 PE, may not be justified in all cases, although usually our calculated B/C remains above 1. In these cases, it may be most appropriate to make a decision case by case also depending on the conditions of the receiving water bodies, as currently required by the Water Framework Directive 2000/60/EC. Data Availability A GIS layer of agglomerations delineated with the calibrated model of Section 2.1 is made available as Supplementary material to this article.

Conflict of Interest The authors declare no competing interests.
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http:// creat iveco mmons. org/ licen ses/ by/4. 0/.

Annex 1. Model calibration
For different combinations of n and 0 , we check how the corresponding simulation of agglomerations compares with the reported numbers and populations in agglomerations above 2000 PE (Table 4), assumed to be the reference.
If x is the vector of reported PEs (first row in Table 4), y is the vector of reported number of agglomerations (second row in Table 4) and (x k , y k ) is the tensor of the corresponding data coming from statistics of the kth simulation of agglomerations, we compute distance metrics as: Ideally, a simulation should have D x = D y = 0 . We choose the simulation with minimum distance to the ideal point D = √ D x 2 + D y 2 . We consider simulations with n = 1 and 0 = 1, 10, 20, 30, 40 and 50, in addition to simulations with 0 = 10 and n = 0, 2 and 5. Their performance in the D x , D y space is shown in Fig. 7. While the simulations are reasonably realistic, they obviously present discrepancies compared to the reported data (Fig. 8).

Annex 2. Details of the cost model
We estimate the cost entailed by providing a collection system and a wastewater treatment system for all small agglomerations below 2000 PE under the following assumptions.
1. The cost of a collection system (Euro/PE) depends on the population in the agglomeration, P, according to the FEASIBLE model (COWI, 2010;OECD, 2004) expenditure function: We assume an operation and maintenance (O&M) cost of 1% of the investment per year. The lifetime of the network is supposed to be 100 years, and we use a discount rate of 2.5% for the investment. 2. The cost of a septic tank, according to OECD, 2004 andupdate in COWI, 2010, is: We assume an O&M cost of 3%, a lifetime of 30 years and a discount rate of 2.5%. 3. The cost of an improved wastewater treatment system is given by the average of costs of various typologies of treatment wetlands according to Pistocchi et al. (2020). The cost of investment (Euro/PE) and operation (Euro/PE/year) is given in the form cost = a P b , with a (coefficient) and b (exponent) given in Table 5. For the plants, we also assume a lifetime of 30 years and a discount rate of 2.5%. With these parameters, we compute an annualised cost function for each typology of treatment wetland (see Fig. 9), and we consider the average of the ensemble.

Fig. 9
Annualised costs for the typologies of treatment wetlands in Table 5 (FBA, forced bed aeration; FRB, French reed beds; HF, horizontal flow; VF, vertical flow; Hybrid denotes a combination of VF and HF; "reuse" or "no reuse" denote typologies where effluents are/are not required to meet standards for water reuse) 4. The additional cost of collection and treatment, compared to the cost of septic tanks without collection, is given by the sum of collection network costs and the above average, minus the cost of a septic tank system. This can be approximated by the function: This is equal to an estimated cost of a treatment wetland system, plus the estimated cost of a small collection network servicing the agglomeration, minus the cost of a septic tank sized for the whole agglomeration.
The cost of N removal in excess of simple carbon removal (secondary treatment) is computed as the cost of a vertical flow treatment wetland as per Table 5, minus the average of costs assumed for secondary treatment. This is well approximated by the function: For P removal, we assume the difference of costs between secondary plants and tertiary plants with P removal only according to OECD, 2004, andCOWI, 2010. This is well approximated by the function: The cost functions of the above additional costs are plotted in Fig. 10.

Annex 3. Greenhouse gas emissions
We conducted an analysis of greenhouse gas emissions (GHG) using a life cycle assessment approach in line with the methods outlined in Parravicini et al. (2022). This allowed defining emission factors for various typologies of wastewater treatment plants for small agglomerations. For technical plants down to a capacity of 500 PE, we refer to typologies discussed in Parravicini et al. (2022) (typologies 4.x in Fig. 1 therein). For plants below 500 PE, we consider the additional typologies listed in Table 6.
For each typology, we compute a PE-specific GHG emission load in kg CO 2 e per PE per year, including direct and indirect GHG emission sources. The balance boundaries of the assessment comprise the building and the operation phases of the wastewater treatment system (assumed life time of 30 years). Nitrous oxide and methane emitted in the water bodies receiving the treated effluents were also included. Sludge stabilisation options at the technical plants were aerobic simultaneous sludge stabilisation (SASS), post-aerobic sludge stabilisation (ASS) and reed beds (RB). Sludge transport and further treatment and/or disposal were out of the assessment boundaries. In the option with primary sludge being transferred to a WWTP for stabilisation, CO 2 e emitted by transportation and treatment in the WWTP were accounted for. For transportation, we applied an emission factor of 0.134 kgCO 2 e per ton per kilometre (Ecoinvent 3.7, market group for transport, freight, lorry, unspecified, Global, 2020). For the sludge stabilisation in the WWTP, we apply the average total CO 2 e emission load estimated for WWTP typologies 5000-30,000 PE with SASS or ASS or RB in Parravicini et al. (2022).
GHG emissions from the construction phase of sewer network down to 50 PE were also included in the assessment and estimated to about 38 g CO 2 e/PE/year, as described in Parravicini et al. (2022). The estimation of methane emissions from sewer networks is challenging due to their geographical extension and variability of conditions ( Ye et al., 2022). Of the methane produced and emitted in the sewer network, only the load remaining dissolved in the sewage and then stripped in aerated treatment systems was considered, applying the approach described in IPCC, 2019. The emitted methane load was estimated at 4 to 6 g CO 2 e/PE/year. Tailored COD, N and P mass-balances for each plant typology, capacity and treatment performance provided the inventory data for the GHG assessment, Electricity demand for the operation of the technical aerobic systems was derived from Ganora et al. (2019) (see also Parravicini et al., 2022). Electricity demand for non-aerated TW (e.g. for wastewater pumping) was also estimated, although being comparatively small. In all options featuring P removal, iron chloride was used per default as flocculant/precipitation agent.
Emission factors for the impact assessment were mainly derived from the IPCC Guideline Refinement 2019 (IPCC, 2019) or other literature sources presented in details in Parravicini et al. (2022).
In addition to the methodological approach described in in Parravicini et al. (2022), we made the following additional assumptions: • Small technical aerobic systems  For typology 5.x.y (x=a,b; y=1 to 6), we set the same assumptions as for the larger wastewater treatment plants 500-5000 PE in Parravicini et al. (2022). An exception is the electricity demand, which was assumed higher for the smaller plants. We assumed, based on literature values, that plants striving N removal exhibit lower direct N 2 O emissions than plants without targeted N removal (detailed information in Parravicini et al., 2022).
• Treatment wetlands (< 2.000 PE) For typologies 6.x.y (x=a,b; y=1,2) were modelled as treatment wetlands, in the options vertical flow (VF) and horizontal flow (HF) wetlands. The removal performance of the two systems was derived from the IWA books on treatment wetlands Dotro et al. (2017) and Langergraber et al. (2019): The removal of COD, NH4+, N and P for VF wetlands are 90%, 90%, 20% and 15% respectively; for HF wetland they are 80%, 25%, 40%, 15% respectively. Applying hybrid multistage systems combining different wetland typologies would allow for higher treatment performance, especially in term of N removal. This option was not evaluated separately, and the results can be extrapolated from the ones for VF and HF wetlands. Enhanced P removal was achieved through flocculation/precipitation with iron chloride. We assumed that the wastewater is treated first mechanically in a primary settler (50% COD removal) and the primary sludge is stabilised and dewatered in a RB. In an additional option, a septic tank was used for primary treatment. In this case, sludge is transported to a WWTP in order to undergo stabilization. We assume a transport distance of 40 km both ways. The additional CO 2 e emissions released at the WWTP were added to the CO 2 e emissions directly related to the TW. Electricity demand for wastewater pumping (calculated assuming a flow of 150 L/PE/ day, and a pumping head of 2 metres) is small compared to the consumption for aeration in technical aerobic systems. Emission embedded in the infrastructure were estimated considering resource use for excavation works (4-5 m 2 /PE, 1-1.5-m deep according to Mander et al., 2014 andDotro et al., 2017). The emission factor for direct N 2 O and CH 4 emissions were derived from the 2013 Supplement on Wetlands to the 2006 IPCC Guidelines for Fig. 11 Population-specific CO 2 e emission loads of the wastewater treatment typologies considered National Greenhouse Gas Inventories (IPCC, 2014, Chapter 6, table 6.4/equation6.2 and table 6.7): for VSSF 0.0025 kg CH 4 /kg COD influent and 0.00023 kg N 2 O-N/kg N influent and for HSSF 0.025 kg CH 4 / kg COD influent and 0.0079 kg N 2 O-N/kg N influent .
Interestingly, rather aerobic systems as VSSF, providing ammonium oxidation but lower N removal, exhibit so far lower N 2 O emission than systems with pronounced anoxic/anaerobic milieu (HSSF). This is not in line with experiences at centralised WWTP, where denitrification can act as a sink for the generated N 2 O (Valkova et al., 2021).
To predict N 2 O and CH 4 emissions in passive systems as treatment wetlands is even more challenging than in centralised WWTPs, where operating conditions can be described and controlled. The IPCC (2014) gives, e.g. for VSSF an uncertainty range of ± 56% for CH 4 and of ±70% for N 2 O. It is pointed out that emission factors measured so far in TW are lower than for centralised WWTPs, so far the plants are properly designed and managed. This can be explained perhaps by the low loading conditions typically applied in these systems.
• Septic tanks with overflow and ground/filter infiltration For typologies 7.x.y (x=a,b; y=1,2), removal performance was set to 50%, 33% and 33% for COD and N and P removal respectively. The sludge production was estimated to 60 g COD/PE.d, 42 g VS/PE.d (COD/VS = 1.42, DWA, 2016), 48 g TS/ PE.d (assuming organic content being ~90% of TS), and 1600 g sludge/PE.d (3%TS). Two options with sludge stabilisation in RB or transport and treatment at a WWTP (assumption: 40 km both ways) were computed. The emission factor for direct methane emissions from septic tank was derived from IPCC, 2019 (table 6.3). GHG emissions embedded in the infrastructure were evaluated on the basis of concrete and stainless steel demand estimated for fictive septic tank (0.25 m 3 /PE, Dotro et al., 2017).
Based on the above assumptions, we obtain the population-specific CO 2 e emission loads shown in Fig. 11.