Stormwater Harvested from Permeable Pavements as a Means to Save Potable Water in Buildings

The main objective of this work is to analyse the potential for potable water savings in university buildings by using stormwater collected from permeable pavements. Six buildings located on the campus of the Federal University of Santa Catarina (UFSC) were selected to obtain monthly water consumption patterns and parking lot areas. The same six buildings were then evaluated considering their location in eight different cities in Brazil, with different rainfall patterns. Simulations using the computer programme Netuno were run to obtain the potential for potable water savings in each building and city combined. The structural design of permeable pavements was also assessed using two methods available in the literature, that is, the American Association of State Highway and Transportation Officials (AASHTO) and Brazilian Portland Cement Association (ABCP). The hydrological-hydraulic design of the permeable pavement was also carried out. The designed thicknesses were compared with the thicknesses obtained using the computer programme Permeable Design Pro. The potential for potable water savings between 18.4% and 84.8% was obtained, depending on the city, building and non-potable water demand considered. For the structural design, the thicknesses obtained by using both methods were similar; however, it was observed that the AASHTO method better represents the pavement model. Regarding the hydrological-hydraulic design, the differences obtained show that the simplification performed for the pavement drainage was in favour of safety. In conclusion, the use of permeable pavements in stormwater harvesting systems is promising, aligning the drainage aid, structural capacity and potential for saving potable water.


Introduction
Climate change is evident in the urban environment. The increase in global temperature correlates with the increase in rainfall peak intensity, which facilitates the occurrence of urban flooding [1]. On the other hand, changes in rainfall patterns have made water management more difficult, affecting the amount of water available. As an example, Cape Town in South Africa had the worst drought in its history in 2015, limiting water consumption to 50 L per day per inhabitant [2]. São Paulo also experienced one of the worst droughts in its history in 2014, which resulted in actions on water management and rationing among users [2,3].
According to the Intergovernmental Panel on Climate Change (IPCC) [4], in 2012, Earth had a temperature increase of 0.85 • C compared to the pre-industrial average. According to the IPCC [5], if fast and deep changes are not adopted, the temperature increase may be higher than 1.5 • C by 2100. Amidst the uncertainties of future weather conditions, the use of alternative sources of water will become necessary. Alternative sources, such as rainwater, greywater, seawater and the use of water-efficient systems, are examples of measures that can be applied [6]. Rainwater harvesting in buildings is already a legal

Simulation Variables
The study began with the selection of six buildings, eight cities and the corresponding water end-uses. After obtaining the water consumption data and the built area of the six buildings, daily rainfall and Intensity-Duration-Frequency (IDF) equations for the eight cities were analysed. The potential for potable water savings in each scenario was evaluated using the Netuno 4 computer programme. The results obtained were grouped and analysed in order to assess the influence of the three parameters on the potential for potable water savings.
After the simulation of potable water savings, the hydrologic/hydraulic design was performed according to a method from the literature. Computer routines were performed using the GNU Octave programming language in order to simplify the process. The re-sults obtained were then compared with those obtained from the Permeable Design Pro programme.
A computational routine was also proposed for the structural design of the pavement, based on the American Association of State Highway and Transportation Officials (AASHTO) method [26]. The results obtained with the routine for the AASHTO method were compared with the thicknesses indicated by the Permeable Design Pro computer programme in order to validate the routine. The design was also performed using the Brazilian Portland Cement Association (ABCP) method [27], a method used in Brazilian designs.

Rainwater Harvesting Model
The rainwater harvesting model was carried out based on models available in the literature. Some authors [24,25] have proposed a permeable pavement water harvesting system composed of: a catchment area, the permeable pavement, underground and upper rainwater tanks and water pumping station. Rainwater falls over the pavement, infiltrates through the permeable layers and it is conducted via drains towards the underground tank. When the upper rainwater tank is about to become empty, that is, less than 10% of its capacity, rainwater from the underground tank is pumped to fill it up. If there is no rainwater available in the underground tank, potable water from the building's water tank is conducted towards the upper rainwater tank by gravity. If both upper and underground rainwater tanks are full, rainwater is discarded in the drainage system. Figure S1 in the Supplementary Materials shows a schematic of the system. Water treatment-chlorination only-is applied in the upper rainwater tank in order to avoid waste of treated rainwater.

Permeable Pavement Model
The permeable pavement model used in the study was based on previous research [28], which made models of a permeable pavement using permeable concrete pavers. The main reason for using data from actual physical models is to be able to guarantee that the simulated pavers exist on the market. The infiltration rate obtained by Ghisi et al. for concrete pavers was 710 cm/h [28], which is approximately twice the lower limit of the Brazilian standards [29], that is, 360 cm/h for permeable coating. Figure 1 shows the pavement layers used in the study with the following specifications: •  [28], in order to comply with Brazilian standards [29].

Buildings
Six university buildings were chosen for the study. They were chosen to represent different built areas and water consumption, and also for having a parking lot area nearby. All buildings were part of the Federal University of Santa Catarina (UFSC). Table 1 shows the six buildings and their parking lot areas and water consumption. The decision to choose university buildings takes into account two points. The first is that the non-potable water end-uses in this type of building are usually higher than in residential buildings. As examples, the studies of Botelho [32] and Rainmap [33] found that the non-potable water end-uses of the UFSC Technology Centre were equal to 69% and 70%, respectively. The second point is the current Brazilian context, in which there are budgetary difficulties in the academic environment and the need for the reduction of expenses. In 2018, for example, UFSC spent approximately 7.3 million reais (R$) on water bills [30].

Cities
The eight Brazilian cities were chosen taking into account that they have a public university. Table 2 shows the different characteristics of the cities and Figure 2 shows the location of the cities on the map of Brazil.

Potable Water Savings
The simulations to evaluate the potential for potable water savings were carried out using the Netuno 4 computer programme [35]. Eight variables were considered, that is, daily rainfall, harvesting area, total water demand, rainwater demand, infiltration coefficient, upper rainwater tank capacity and underground rainwater tank capacity.
Netuno 4 is a computer programme, created by Ghisi and Cordova [35], used for sizing rainwater tanks and estimating the potential for potable water savings when there is rainwater harvesting. The simulation uses defined variables, that is, it does not rely on a stochastic model of variables and statistical outputs. The main use and definitions are for rooftop rainwater harvesting, but it has all the necessary inputs for a water balance simulation, which can be applied similarly in the permeable pavement system design.

Daily Rainfall
Rainfall data from 2003 to 2017 were obtained through the Meteorological Database for Education and Research (BDMEP) [36] of Brazil, in CSV format. The first flush was set to zero since the permeable pavement is meant to be used to filter the stormwater and there are no diverters for the first flush.

Harvesting Area
The harvesting area refers to all areas that contribute to the volume of rainwater harvested. For all six buildings, the areas of the parking lot and impermeable surrounding surfaces were considered, that is, the total harvesting area was equal to the sum of permeable and impermeable surfaces. The water on the impermeable surface flows towards the permeable pavement and then infiltrates it.

Total Water Demand
The total water demand in the building refers to all water consumption during a given period-daily, monthly or yearly. Thus, for each of the six buildings, the average monthly consumption was obtained through a historical series of water consumption data. Data were obtained via compiled files from the UFSC's monitoring of water supply and sanitation [30].
Monthly water consumption was converted to daily data in order to be used in Netuno 4. This procedure was carried out by dividing the monthly consumption by the number of working days in each month plus one day. This extra day was to take into account the weekend; it was considered that all the weekend days in a month represent one working day. This consideration was made based on the fact that the university has a higher occupation over the week. Thus, consumption on a weekend day represents between 10.0% and 12.5% of the consumption on a working day. Equation (1) shows the consumption on a weekday, and Equation (2) shows the consumption on a weekend day.
where: Dcwd i is the daily consumption on working days in the month i (L/day); Mc i is the monthly consumption in the month i (L/month); Dcwk i is the daily consumption on weekends in the month i (L/day); wd i is the number of working days in the month i (days); wk i is the number of weekend days in the month i (days).

Rainwater Demand
Rainwater demand represents the part of the total water demand that requires nonpotable water. It is one of the parameters required by Netuno 4 to indicate the maximum potential for potable water savings that the system can obtain. It is related to the percentage of water end-uses that can be supplied with non-potable water, such as toilets, urinals and garden irrigation.
The Brazilian standard for non-potable water uses in buildings is NBR 16783 [37]. The possible non-potable uses included in the standard are car washing, toilet flushing, urinal flushing, cleaning of outside areas, irrigation, ornamental use and cooling. For these specific uses, the standard generally addresses seven quality parameters (pH, Escherichia coli, turbidity, biochemical oxygen demand, residual chlorine, electrical conductivity or total dissolved solids, and total organic carbon); for other uses it requires a specific study of quality parameters depending on the water source.
In this study, three rainwater demands were considered based on the studies of Botelho [32], Rainmap [33] and also on a survey carried out at UFSC library in 2018, which was not published. Based on such studies, three rainwater demand scenarios were analysed, that is, 69%, 77% and 85% of the total water demand.

Infiltration Coefficient
The infiltration coefficient refers to the percentage of stormwater that is harvested in relation to the amount of stormwater that falls onto the harvesting area. Since this study deals with permeable pavements, the infiltration coefficient refers to the relationship between the amount of stormwater that falls onto the pavement and the amount that reaches the lower rainwater tank. The value used in the study was the same as that obtained by Ghisi et al. [28] in their study on permeable pavements with pervious concrete pavers, that is, 88.1%.

Upper Rainwater Tank
It was defined that the upper rainwater tank in each building has a capacity equal to the average daily rainwater demand for that building. For the pumping pattern, it was defined that the pump will be turned on every time the amount of water in the upper rainwater tank reaches 10% of its total capacity.

Underground Rainwater Tank
The pavement will serve as a filter for stormwater and also as temporary storage, but the bottom drains in the reservoir course direct the stormwater to the underground rainwater tank. The capacity of the underground rainwater tank was sized using Netuno 4.
As the programme iterates in the simulations, it calculates the potential for potable water savings and, from the results, indicates the optimal rainwater tank capacity.
The simulations were performed for underground rainwater tank capacities ranging from 0 to 50,000 L at intervals of 1000 L. For each lower tank capacity, a different potential for potable water savings was obtained. An optimal rainwater tank capacity was defined as a function of the potential for potable water savings and the underground rainwater tank capacity, that is, an added potential lower than 1%/m 3 .

Hydrologic/hydraulic Design
The hydrologic/hydraulic design consists of defining the thickness of the reservoir course layer. It was performed for each of the eight Brazilian cities in order to size for the city's design rainfall, which was obtained via the Intensity-Duration-Frequency (IDF) equation or the design rainfall height. The ratio between the permeable area and the impermeable area of the pavement was assessed for each building. Another relevant point is the definition of the pavement bottom drain, which serves for the specific outlet flow and influences the thickness required for the reservoir course.

IDF Equations
The IDF equations are necessary for assessing design rainfall, which is used in the design of any type of hydraulic structure. The knowledge about design rainfall is necessary for making the structures compatible with future rainfall demands. As far as design rainfall is concerned, according to Back [38], some methods can be used to define the rainfall curve mathematically. In Brazil, the IDF curve is usually used.
The IDF equations were used in the potential format, as shown in Equation (3). Table 3 shows the IDF curve parameters of each of the eight cities.
where: K, m, n and d are coefficients that vary according to the place; T is the return period (years); i is the average maximum intensity of the design rainfall (mm/h); and t is the duration of the design rainfall (minutes). To calculate the average maximum intensity using the IDF equations, it is also necessary to define the duration of the rainfall and the return period. The return period and duration of rainfall have been chosen according to Brazilian standards [29]. The minimum return period was taken to be ten years, with a minimum duration of rainfall equal to 60 min.

Permeable to Impermeable Area Ration (R)
Another important parameter for performing the hydrologic/hydraulic design is the ratio between the total area and the permeable area of the pavement. Equation (4) shows how the ratio parameter (R) was calculated. R values between 1 and 2 can be used, that is, the impermeable area can range from zero to a figure equal to the permeable pavement area, not exceeding it [47].
where: R is the ratio of the total area to the permeable area (dimensionless); A per is the area of the permeable surface (m 2 ); and A imp is the area of the impermeable surface (m 2 ). For each building, parking lot sketches were made in order to estimate R. The pavement slope turns the runoff towards the edges of the track, where there is drainage through the permeable pavement. Figure S2 in the Supplementary Materials shows the sketch of the Technology Centre parking lot, as an example.

Bottom Drains
The specific output flow is the hydrologic/hydraulic design parameter that defines the drainage capacity after infiltration into the pavement. Hammes et al. (2018) used a simplified method to evaluate the specific output flow, defining it as a constant figure, which corresponds to the total amount of water falling on the floor divided by the permeable area over 24 h. The time necessary to empty the reservoir course was defined by Hammes et al. [20], based on Tomaz [48], as at least 24 h. Equation (5) shows the simplified method for calculating the specific output flow, which was used herein.
where: q s is the specific output flow (m/h); i is the average maximum rainfall intensity (mm/h); t is the duration of the design rainfall, considered equal to 1 h (h); A tot is the total pavement area (m 2 ); A per is the permeable pavement area (m 2 ); and TE is the time necessary to empty the reservoir course (h).

Hydrologic/hydraulic Design Method
The method of Hammes et al. [20] was used herein for hydrologic/hydraulic design. The method is very similar to that in the Tennessee Permanent Stormwater Management and Design Guidance Manual [47] and the Brazilian standard [29]. The result of the hydrologic design is the thickness of the reservoir course required for temporary rainwater storage. After the storage, the water is drained by gravity through underground drains to a lower rainwater tank; if the rainwater tank is full, the water flows into the drainage system. Equation (6) shows the method of Hammes et al. [20].
where: H rc is the required thickness of the reservoir course (mm); t is the duration of the design rainfall (min); R is the ratio of the total area to the permeable area (dimensionless); q s is the constant specific outlet flow of the pavement (m/h); i is the average maximum intensity of the design rainfall (mm/h); and η is the porosity of the reservoir course (dimensionless).

Structural Design
The structural design of the permeable pavement is necessary in order to guarantee that the thickness of the layers is adequate to dissipate the stress generated by the traffic of vehicles. Two methods used for design were assessed: the ABCP method [27] and the AASHTO method [26]. The ABCP method is a simple design method that has as its input parameters the Equivalent Single Axle Load (ESAL) data, estimated for pavement life, and the California Bearing Ratio (CBR) of the subgrade. Only these two parameters are necessary to design the reservoir course. As a restriction for the subgrade material, the method indicates a minimum CBR equal to 20% for less than 500,000 ESALs, and 30% for more than 500,000 ESALs.
The AASHTO method requires more input parameters than the ABCP. The AASHTO parameters are: the number of ESALs, also estimated for the pavement's lifespan; subgrade CBR; thickness and material of all layers; pavement's drainage condition; Present Serviceability Index (PSI) and statistical parameters. It is a more comprehensive method, with application to flexible pavements in general, while the ABCP method is restricted to pavements with interlocking concrete pavement. The thickness and materials of the pavers, the bedding layer and choker course were chosen following the permeable pavement construction manuals [47,49] and the Brazilian standard [29]. Drainage coefficient, PSI and statistical parameters were chosen according to the AASHTO manual [26] in order to agree with the type of pavement and its use. Equations (7) and (8) show the AASHTO method for structural design.
where: W 18 is the number of ESALs over the pavement lifespan (dimensionless, 18 stands for the reference axle load of 18,000 lb); and Z r is the number of standard deviations needed to obtain, in a standard normal distribution, the probability of structural failure. For example, at 95% reliability and 5% chance of failure, Z is −1.645 (dimensionless); S o is the standard deviation of the supported value of traffic over the service life, taking into account the variability of construction, assumed parameters of materials and the assumed number of future vehicles the pavement will support (dimensionless); SN is the Structural Number of the pavement, which takes into account the thickness and materials of layers (inches); ΔPSI is the difference of Present Serviceability Index in the beginning and at the end of the pavement's life cycle (dimensionless); M r is the resilient modulus of the subgrade layer (psi); a (i) is the structural coefficient of the pavement layer i (dimensionless); d (i) is the thickness of the pavement layer i (inches); m (i) is the drainage coefficient of the pavement layer i (dimensionless); n is the number of layers in the pavement. As the AASHTO structural design method uses the resilient modulus of the layers, it was necessary to make a correlation between the resilient modulus and CBR to be able to compare with the ABCP method. Thus, the correlation shown in Equation (9) between CBR and the resilient modulus was used [50,51].
where M r is the Resilient Modulus (psi) and CBR is the California Bearing Ratio (dimensionless).

Model Validation
The permeable pavement was also designed using Interlocking Concrete Pavement Institute (ICPI) programme Permeable Design Pro [51]. The purpose was to compare the results obtained through the hydrologic/hydraulic design method found in the literature with those obtained from the programme. The programme uses an hourly analysis of mass balance in the hydrologic/hydraulic design.

Potential for Potable Water Savings
The potential for potable water savings was estimated for each combination of building, city and rainwater demand, totalling 144 simulations. The resulting analysis was performed in order to understand how these three parameters influence the potential for potable water savings.

Rainfall Data
Rainfall data were obtained for each of the eight cities and were treated and organised for better visualisation. Three main characteristics that differ among the cities were observed. The first is the average annual rainfall, which varied between 1454 and 2287 mm/year. Porto Alegre has the lowest annual rainfall and Manaus has the highest. The second characteristic is the distribution of monthly rainfall during the year. Some of the cities had a fairly constant monthly rainfall, while others had well defined dry and wet periods. As an example, Brasília has a dry period between May and September. The last characteristic observed is the maximum rainfall in one day, which was compared with the results obtained using the IDF curves, in order to understand whether the curves match the local rainfall profile. Table 4 shows the characteristics observed for each of the eight cities.

Water Consumption
Water consumption was assessed in order to understand two points: differences in the consumption of the buildings over time and distribution of consumption during the year. The analysis of the variation over the years serves to check divergences due to leaks or retrofits in plumbing systems. The analysis of distribution during the year serves to create the daily water consumption model, used in the computer simulations. Figure S9 in the Supplementary Materials shows the water consumption for the six buildings. Figure 5 shows the monthly average water consumption during the year for the six buildings. Based on the monthly average water consumption in each building, it was possible to estimate the daily water consumption. Figure 6 shows, as an example, the daily water consumption of the Technology Centre building used as input data in the computer simulations. Figures S10-S14 in the Supplementary Materials show the daily water consumption for the remaining buildings. Daily water consumption in Technology Centre building considered in the computer simulations. Figure 7 shows, as an example, the potential for potable water savings for the six buildings in the city of Florianópolis considering the three different rainwater demands. The ideal underground tank capacity for the Technology Centre was 24,000, 21,000 and 23,000 L for rainwater demands equal to 69%, 77% and 85%, respectively. Figures S15-S21 in the Supplementary Materials show the results for the other cities.  Figure 8 summarises the average potential for potable water savings for the six buildings, eight cities and three rainwater demands. Such savings ranged from 18.48% to 82.81%. The lowest potential occurred in the architecture and urbanism building when the rainwater demand was equal to 69% of the total water consumption in Belo Horizonte. The highest potential for potable water savings occurred in the Department of Architecture and Engineering Projects building when the rainwater demand was equal to 85% of the total water consumption in Florianópolis. Table 5 shows, as an example, the optimal underground tank size for each building in Florianópolis.  Three factors were observed in the results:

Potable Water Savings
The first factor is the water consumption in each building. The Department of Civil Engineering and the Department of Architecture and Engineering Projects buildings presented a much lower consumption compared to the other buildings. Both buildings also obtained the highest potential for potable water savings. All other buildings presented water consumption of the same order of magnitude with different distributions throughout the year, which resulted in different potentials for potable water savings for each city. It should be noted that the harvesting area in these buildings was not proportional to water consumption, and therefore buildings with smaller water consumption showed greater availability of rainwater.
The second factor is the rainfall pattern. The cities of Belo Horizonte and Brasília had a distinction between the rainfall volume in the dry period, which occurred in winter, and in the rainy period, which occurred in summer. This characteristic can be seen by analysing the standard deviation of the monthly average rainfall, which is, respectively, 96.12 and 113.15 mm/month (as shown in Table 4). For both cities, the monthly average rainfall in June, July and August did not exceed 20 mm/month, while in November, December and January it was over 200 mm/month. The cities of São Paulo, Manaus and Recife also did not show consistency in the rainfall pattern, with standard deviations equal to 84.71, 100.31 and 124.55, respectively. Finally, the cities in the south, that is, Florianópolis, Curitiba and Porto Alegre, showed low variation in the monthly average rainfall; the standard deviations were much lower than in the other cities, that is, 39.79, 33.49 and 16.73, respectively.
The third factor is the annual rainfall. The cities of Manaus and Recife have the highest annual rainfall, that is, 2287 and 2230 m 3 /year, respectively. Both cities have considerably higher rainfall than other cities, which varied between 1454 and 1789 m 3 /year. Thus, even if they have a variation in the monthly rainfall pattern, Recife and Manaus showed the greatest potential for potable water savings in some of the combinations between building and rainwater demand.
By combining these three factors, one can explain the results obtained. Low consumption buildings, such as the Department of Architecture and Engineering Projects, obtained the best results in all combinations of city and rainwater demand. Within the cities, the best results were obtained for Florianópolis, Porto Alegre and Curitiba due to low variation in rainfall. High consumption buildings, such as the Technology Centre, showed the best results in the cities with higher rainfall, that is, Manaus, Recife and Florianópolis. It can also be observed that Florianópolis presented the best result for low water consumption buildings and the third-best for high consumption buildings, as it has high annual average rainfall and low variation between months.
Differences in the water consumption pattern between buildings also influenced the results. As an example, the Technology Centre presented the best results in Recife and the University Library presented the best results in Manaus. The difference in the results occurred as rainfall and water consumption distribution are different in both buildings and cities.
One also assessed how different rainwater demands influenced the results obtained for each building. In order not to extend the analysis of the results too much, the average potential for potable water savings for all cities was used. Table S1 in the Supplementary Materials shows the potential for potable water savings for each building and the rainwater demand. A linear regression was performed and the slope was calculated to assess the impact of the rainwater demand on the potable water savings potential. The linear regression slope was higher for the low water consumption buildings; therefore the rainwater demand has a greater impact on the potable water savings potential of low water consumption buildings. Although this impact is lower in high water consumption buildings, rainwater harvesting in such buildings can lead to significant monetary savings.

Hydrologic/hydraulic Design
The hydrologic/hydraulic design was performed using the method proposed by Hammes et al. [20], with the author's simplification for the specific output flow. Rainfall data were assessed for the last 15 years using data from BDMEP, based on Geraldi and Ghisi [52]. IDF curve data were obtained from published works that show local coefficients based on a locally measured rainfall database, as shown in Table 3. Table 6 shows the rainfall intensity, the number of rainfall events above the design rainfall and the extreme rainfall observed in the BDMEP data. R values equal to 1.50, 1.75 and 2.00 were used. They were obtained from sketches of the parking lots for each building. Table 7 shows the areas and the R values obtained from the sketches. Figure 9 shows the results obtained for the design of the reservoir course using the method of Hammes et al. The reservoir course thickness varied between 170 and 493 mm, depending on the city and building.

Structural Design
Structural design using the AASHTO method [26] was performed in a routine written in GNU Octave, while the ABCP method [27] was assessed using an abacus. Parameters used in the methods are shown in Table 8. The number of ESALs was chosen based on light traffic, defined according to PMSP [27], which has design traffic of 100,000 ESALs. CBR of the subgrade was chosen to be 5%, according to the soil profile at the site, which presented very weak structural support capacity with low Standard Penetration Test (SPT) values. Thus, a choice was made to replace the material of the subgrade with soil with a bearing capacity equal to the minimum recommended for permeable pavements [51]. It is also worth mentioning that the AASHTO method uses the resilient modulus and not the CBR value, and thus a correlation was used to calculate the required thickness, according to Equation (8). In an executive project, it is necessary to find out the real resilient module of the layers materials. Figure 10 shows the variable results according to the CBR of the sub-bed and drainage coefficient. The ABCP method resulted in a good estimation of the pavement thickness, presenting results similar to those of AASHTO for good and excellent drainage. However, in cases where the accumulation of rainwater for longer periods in the reservoir course is desired, the method becomes inefficient as it does not consider the drainage period. For CBR between 6% and 8%, the ABCP method led to results lower than those from AASHTO with excellent drainage and may represent insufficient thickness for the layer. It is also interesting to analyse the fact that, for CBR above 16%, regardless of the quality of drainage, all structural design is minimal.
The system was defined to be of fast storage type, that is, it does not maintain the rainwater for long periods. Thus, it was considered that the pavement should present good drainage, which indicates a thickness of 20.3 cm according to the AASHTO method and 17.5 cm according to the ABCP method. We chose a thickness equal to 21.0 cm in order to address the design of the AASHTO method.

Model Validation
The pavement was then modelled using the Permeable Design Pro computer programme. Two criteria are analysed through the design method of the programme, that is, structural and hydrologic/hydraulic requirements. The structural method is the same as that used in our study, that is, the AASHTO method. The drainage considered in the programme is good drainage, which corresponds to the outflow of water from the reservoir course in up to one day. Thus, the thickness obtained using the AASHTO method, 20.3 cm, did not differ from the result obtained from the Permeable Design Pro, that is, 20.5 cm. The final design thickness remained equal to 21.0 cm.
The hydrologic/hydraulic design method used in the programme is more precise than the method of Hammes et al. [20], requesting additional data. The programme creates a water balance on an hourly basis. In order to obtain the water balance, it uses the permeability of the pavement layers, drainage system parameters, pavement geometry and daily rainfall distribution.
Layer permeabilities were based on the programme standard figures for the choker course and reservoir course materials [51] and on Ghisi et al. for the paver layer [28]. Drainage data were defined according to the recommendations of permeable pavement manuals [47,49]. The pavement geometry was obtained through the parking lot sketches obtained through Google Earth Pro photos [31]; the streets were considered an impermeable area and the parking lots as a permeable area. Finally, the daily rainfall distribution was defined as the NRCS type III rainfall, because, among the options available in the programme, it presents the highest rainfall concentration, that is, it is the critical model. For each city, we adopted the rainfall curve with the intensities obtained from the IDF equations. Table 9 shows the parameters used in the Permeable Design Pro computer programme (ICPI, 2019). Table 9. Hydrologic/hydraulic design parameters.

Parameter Figure Considered
Pavers Permeability (mm/s) Maximum water height in the reservoir course 100.0% In conclusion, 48 output reports were obtained, one for each combination of city and building. For the 48 combinations, the hydrologic/hydraulic design resulted in a reservoir course thickness smaller than that obtained in the structural design. Figures 11-14 show, as an example, the results obtained for the critical case, which occurred at the Technology Centre parking lot in Florianópolis.    The simplification used by Hammes et al. [20] was much in favour of safety, with a specific output flow much lower than the actual drain capacities. The divergence happens because the simplified model only uses one drain to collect all the branches of the drainage, without distinction among the branches that reach the rainwater tank.

Conclusions
In this work, an analysis of three aspects of permeable pavements was performed to improve knowledge about the technique. The first aspect was the evaluation of the potential for potable water savings as a result of stormwater harvesting. For this purpose, analyses were performed on the rainfall pattern of eight Brazilian cities and the consumption of six university buildings. The potential for potable water savings ranged from 18.48% to 82.81%. Such disparity in results is due to both large differences in water consumption patterns and differences in rainfall patterns. The city of Florianópolis presented the best results in terms of potable water savings potential overall, as it has a fairly constant rainfall distribution pattern and a high annual rainfall.
The second aspect verified was the structural design. For this, a GNU Octave code was written for the AASHTO method. The thicknesses according to the ABCP structural design method were also assessed. The results were validated through those obtained from the Permeable Design Pro programme, which also uses the AASHTO methodology for the structural design of permeable pavements.
There was little variation among the structural results of the ABCP method, the AASHTO code and the Permeable Design Pro, especially for the CBR of 5%. It should be noted that the ABCP method does not indicate the use of a choker course and does not consider drainage conditions, which may modify the result in other system settings. The absence of more parameters makes it less exact for the evaluation of permeable pavement systems with stormwater harvesting purposes.
The results of the hydrologic/hydraulic design using Permeable Design Pro diverged greatly from the values found using the method of Hammes et al. The main difference is that the programme considers an hourly water balance, requesting drainage specifications. The programme specific output flow was higher than the output flow obtained via the simplification that Hammes et al. considered. The difference obtained in drainage resulted in smaller reservoir course thicknesses.
It can be concluded that the results shown herein are examples of how a permeable pavement system can save potable water in buildings. It is believed that, for public universities, the system has a great potential for reducing the water bill costs and relieving the urban drainage systems.
Supplementary Materials: The following are available online at https://www.mdpi.com/article/ 10.3390/w13141896/s1, Figure S1: Rainwater harvesting schematic with water flows, Figure S2: Permeable and impermeable areas of the parking lot in the Technology Centre, Table S1: Potential for potable water savings and rainwater demand influence in the potential for potable water savings, Figure S3: Average, minimum and maximum monthly rainfall in the city of Florianópolis over 2003-2017, Figure S4 Figure S10: Daily water consumption in the Department of Architecture and Engineering Projects considered in the simulation, Figure S11: Daily water consumption in the Civil Engineering Building considered in the simulation, Figure S12: Daily water consumption in the University Library considered in the simulation, Figure S13: Daily water consumption in the University Administration Building considered in the simulation, Figure S14: Daily water consumption in the architecture and urbanism building considered in the simulation, Figure S15: Potential for potable water savings for the six buildings located in São Paulo, Figure S16: Potential for potable water savings for the six buildings located in Curitiba, Figure S17: Potential for potable water savings for the six buildings located in Brasília, Figure S18: Potential for potable water savings for the six buildings located in Belo Horizonte, Figure S19: Potential for potable water savings for the six buildings located in Manaus, Figure S20: Potential for potable water savings for the six buildings located in Recife, Figure S21: Potential for potable water savings for the six buildings located in Porto Alegre.