The performance of onsite rainwater harvesting (RWH) system in Mediterranean climate was assessed. A stochastic model quantifying the necessary storage, as a function of rainfall (frequency, depth), roof area, residents’ number, specific water use (toilet flushing, laundry) and the required efficiency was developed. Two performance indicators were calculated: water saving efficiency (RSE) – proportion of water used supplied by the RWH system; and rainwater use efficiency (RUE) – proportion of rainwater actually used. The maximum storage capacity and WSE decreased with increasing number of residents for a given roof area, and with an increasing roof area for constant number of residents. For variable storage volume, RUE increased with increasing storage capacity and reached a maximum with an increase in residents’ number and a decrease in the roof area. The model enables to determine WSE and RUE for specific storage volumes or to determine the desired WSE and calculate the necessary storage.
Onsite rainwater harvesting (RWH) is an ancient method which served as an alternative source of water in many places in the Middle East and all around the world. However, with the establishment of central water supply systems, the use of onsite RWH systems has generally stopped. Today due to increased water shortage on one hand, and urban flooding on the other, there is a renewed interest in onsite RWH. Interest in onsite RWH extends from water-scarce regions where the motivation is increasing the amount of available water, to water-ample ones where the motivation is primarily prevention and reduction of urban runoff as well as environmental awareness.
RWH has been acknowledged as a potential source to supply water and to promote significant potable water savings (Ghisi et al. 2007; Gires & de Gouvello 2009). Rainwater, which is a renewable freshwater source, may be used in various non-potable applications at the household level in urban areas. Rainwater, being the main source of freshwater in both natural and human-managed ecosystems, has significant untapped potential for being harvested (Umapathi et al. 2013). Numerous studies investigating the harvested rainwater quality were conducted in Australia, Canada, Denmark, Germany, India, Japan, Spain, New Zealand, Thailand, and the United States (Uba & Aghogho 2000; Evans et al. 2006; Despins et al. 2009; Jones & Hunt 2010; Farreny et al. 2011). However, less information and clear definition on rainwater tank sizing are available (Ghisi 2010; Ward et al. 2010, 2012; Campisano & Modica 2015). The correct tank sizing is important in order to avoid extra costs when the tank is oversized and low efficiency when it is undersized. Several tools were developed for estimating the required tank size and to predict the system performance. For instance, Jenkins et al. (1978) developed two behavioral algorithms to describe the operation of a RWH system during a given time interval. The first algorithm is yield after spillage (YAS), where the amount of water provided by the rainwater collection system, in which the withdrawal occurs after the rainfall has been added to the storage facility and spillage, has been determined. Whereas the second, yield before spillage (YBS) algorithm, assumes that the demand is withdrawn before spillage is determined. Mitchell (2007) investigated the impact of these two computational operating rules and reported that the YAS operational rules underestimated yield and volumetric reliability for a given set of storage, while the YBS operational rule produced an overestimate. Therefore, the author recommended using the YAS operating rule in preference to the YBS, as its estimates are less sensitive to variations in storage and provide conservative values for yield and volumetric reliability. Fewkes (1999) used collected data to verify and refine a rainwater collection sizing model based on the YAS algorithm. The refined model was used to develop a series of dimensionless design curves relating collection area, water demand, rainfall level, system efficiency and storage volume. Fewkes & Butler (2000) evaluated the accuracy of behavioral models, for the sizing of rainwater collection systems using different time intervals and different reservoir operations. Villarreal & Dixon (2005) generated a computer model to quantify the water saving potential of rainwater collection by analyzing the water saving efficiency (WSE). The analysis of several scenarios allowed the authors to suggest suitable sizes of rainwater tanks. Khastagir & Jayasuriya (2010) presented a methodology for optimal sizing of rainwater tanks considering the annual rainfall at the geographic location, the demand for rainwater, the roof area and the desired supply reliability. Ghisi (2010) analyzed the influence of rainfall, roof area, number of residents, potable water demand and rainwater demand on rainwater tank sizing, by using computer simulations. The author indicated that rainwater tank sizing for houses must be performed for each specific situation, i.e., considering local rainfall, roof area, potable water demand, rainwater demand and number of residents. Raimondi & Becciu (2014) who estimated the tank size by using an analytical probabilistic approach showed that the probability of complete rainwater use strongly depends on the period of regulation.
The objective of this study was to develop a novel stochastic model to estimate the optimal rainwater tank size depending on specific demand for rainwater, number of residents and the catchment size (roof area), when the daily rainfall at the location area was considered as the stochastic parameter. Daily rainfall was taken from historical data, and probability functions were derived for each calendar day. The model was developed for Mediterranean climate (Haifa, Israel), characterized by long dry summers (literally) and winters with highly variable rainfall patterns. Nevertheless, a similar methodology may well be implemented to other climatic regions.
MATERIALS AND METHODS
R(t) is the predicted daily rainwater depth (mm/d) for day t (t = 1 → 365).
The Max function was added in order to ascertain that rainwater depth is always non-negative (the result of the polynomial approximation can become negative below a certain probability threshold).
– at, bt, ct, dt, et, ft and gt are the polynomial coefficients for day (t). These were obtained for each day by minimizing the squared error function.
– P(t) is the probability (0 ≤ P ≤ 1) - a number which is randomly chosen by the model (uniform distribution).
Roof area size and type. Five roof area sizes were simulated in the model: 75, 100, 150, 200 and 400 m2.
where - the specific daily rainwater runoff generated (generated runoff divided by the roof area) (L/(m2·d)) for each roof type (i); ai – the slope of the line (L/(mm·m2)), for each roof type; - daily rainwater depth at day t (mm/d) and bi – the intercept with the Y-axis (L/m2·d), for each roof type.
Six possible residents population sizes in a single house were examined in the model 4, 8, 12, 24, 48 and 64 residents.
The model was written in MATLAB and based on YAS algorithm, in which the water supplied from the storage tank after rainfall has been added to the storage facility (Fewkes 1999; Fewkes & Butler 2000). The model is a daily model, i.e. it uses a daily time-step. Since the rainy season starts after a long dry summer (105–220 consecutive dry days), the storage tank at the beginning of winter (the rainy season) was considered as empty.
If the storage tank is completely filled (V(t) > VMax), the excess rainwater generated is released as overflow () and the water volume available for the next day is VMax.
In this case (demand for rainfall is higher than available rainwater), after the above calculation V(t) is set to 0.
The study was divided into two stages. In the first stage the maximum volume in which all the roof rainwater runoff is used (in other words, system utilization efficiency of 100%) was calculated by the simulation model. Figure 2 depicts a schematic flowchart for this simulation. A total of 30 combinations were simulated and analyzed: five roof sizes (75, 100, 150, 200 and 400 m2) and six residents’ population sizes (4, 8, 12, 24, 48 and 64 people). Each option was run for a whole year with a daily time-step. 100 year-long random simulations were performed with stochastic (Monte-Carlo) rainfall input as described above yielding 100 possible maximum storage tank volumes (for each option) that were statistically analyzed.
In the second stage, lower volumes of the storage tank (lower than the maximum values found in the first one) were set prior to the simulation for evaluating the WSE and rainwater use efficiency (RUE) of the system. WSE is defined as the proportion of water demand (of the two relevant uses) provided by the harvested rainwater (Equation (8)). It should be noted that the WSE was calculated only for the rainy season (September–April) The RUE quantifies the proportion of the rainwater harvested that was actually used (Equation (9)).
RESULTS AND DISCUSSION
As aforementioned, the model randomly selected probability and calculated rainwater depth for each day (Equation (1)). The average simulated rainwater depth was 574 mm/y (100 model runs), fell in line with measured values at Haifa Port meteorological station (538 mm/y). The range of the simulated rain depth was 300–900 mm/y very similar to the range of the measured data (292–925 mm/y). Further, the stochastic simulation results were not found to be statistically different from the measured data, indicating satisfactory representation of measured data.
The effect of roof type on the generated runoff, as calculated from the measured data, is presented in Table 1 which presents the parameters of the linear empirical equation (Equation (3)). The correlations between runoff and rainfall were high for all examined roof types (R2 ≥ 0.93 and p < 0.05). ‘a’, the regression line slope, expresses the relationship between rainfall and the generated roof runoff after runoff commenced. Hence, the closer ‘a’ is to 1, the higher the proportion of rainfall that is converted to runoff. Of the three roof materials examined tiles had the highest rain to runoff conversion rate (a = 0.91) while concrete had the lowest (a = 0.78). R(y = 0) is the minimum amount of rainfall needed for runoff to commence (R(y = 0) = −b/a). R(y = 0) actually represents the depression storage of the roof, which is analogous to depression storage in open spaces. For the examined roof types, runoff from the concrete roof started after 2.3 mm of rain as compared with 0.37 and 0.041 mm for the tile and steel-sheets roofs, respectively. The findings were expected as steel-sheets have less and smaller crevices and less water is consumed for wetting the roof material than tile or concrete roofs. To summarize, the runoff from the concrete roof started after the largest rainfall depth as it required the largest amount of rainfall for filling small depressions in the roof before runoff commenced and it generated the lowest volume of runoff for each rainfall event. The tile roof generated the largest volume of runoff for each rainfall event (largest ‘a’), although runoff from the steel-sheets roof started after the lowest rainfall depth (lowest R(y = 0)). The high runoff generated by the tile roof is probably attributed to its high longitudinal slope (30%). It should be noted that apart from depression storage on the roof itself delivery losses from the roof to the storage tank were not considered.
|Roof type||a l/(mm·m2)||b l/(m2·d)||R(y=0)a mm||nb||R2|
|Roof type||a l/(mm·m2)||b l/(m2·d)||R(y=0)a mm||nb||R2|
aR(y=0) – rainfall depth above which runoff commences.
bn – number of rainfall events (observed).
Maximum volume of the storage tank and potential WSE
For brevity in the following sections only the results for the concrete roof (most common roof type in Israel) are presented.
RUE, WSE and storage tank volume
In many cases the maximum storage volume (for which RUE is 100%) leads to a large storage tank volume. However, RWH systems do not always require the maximum volume and in most cases, lower volume tanks generate high (or at least satisfactory) WSE. Therefore, as aforementioned, in the second stage of the research the model was run with varying storage tank volumes and the WSE and RUE were calculated for each of the 30 combinations (roof area × number of residents).
The RUE in a four-person house generally decreased with an increase of the roof area, since the daily water demand was lower than the generated runoff (Figure 6 bottom left). The RUE was found to be more sensitive to tank size than the WSE, with a much larger decrease of the RUE in the small tank volume (0.1·VMax) from 65% to 34% (∼50% decrease) than in the larger tank volume (0.9·VMax, from 94% to 89%, ∼5% decrease). As the number of residents in the house increased the decrease in the RUE diminished, yet here again the decrease in the small tank volume (0.1·VMax) became significantly larger than all other tank volumes. For example, for a 64-person house (Figure 6 bottom right) the RUE of a 0.1·VMax tank decreased from 99% to 68% as the roof area increased from 75 to 500 m2 (32% decrease), while the RUE of a 0.3·VMax tank decreased from 99% to 84% (15% decrease) for the same increase in the roof area. It should be noted that for the combinations examined (roof area × number of residents) the WSE decreased significantly with an increase of number of residents (meaning that a lower proportion of the water consumption was supplied by the rainwater harvested, while the RUE increased (but in a less pronounced manner) with increasing number of residents, meaning that higher proportion of the roof runoff was used. The results demonstrate the importance of using a model for determining the right tank volume (avoiding extra costs due to oversizing of the storage tanks, while keeping a satisfactory efficiency of RWH system).
A stochastic model to quantify the optimal size of rainwater storage tanks for residential homes was developed based on daily rainwater depth, non-potable domestic water demand (toilet flushing and laundry), number of residents and roof area, where rainfall was considered as the stochastic parameter. Daily rainwater depth was calculated from historical data, and probability functions were derived for each calendar day. Using this, the effect of the variable daily rainwater was studied while keeping the seasonal patterns.
Quantifying the storage tank volume based on the WSE, emphasizes the importance of considering the rainwater pattern, roof area, specific water demand and the number of residents. The model output exhibited good correlation between the WSE and storage tank volume, following a saturation curve pattern. This relationship is significant since it can be used for estimating the required storage tank depending on the desired WSE. It was demonstrated that in many cases the maximum storage volume is not really needed and smaller volumes can achieve almost the same efficiencies (WSE and RUE). For example, one can assume a specific storage tank volume and by running the model predict the WSE, or determine the desired WSE and calculate the required tank volume. The model demonstrated that no single optimal storage tank volume exists, since it depends on local weather patterns, roof size, specific demand for the harvested rainwater, number of residents in the house and the desired WSE. This means that for each setting the tank volume may be different and should be determined by the model or a similar simulation tool. The model was developed for Mediterranean climate (Haifa, Israel), but the same methodology may well be implemented for other climatic regions. In addition, the model can be used for examining the effects of extreme weather (or climate change) on WSE and RUE, using the developed probability curves to generate rainfall values for different return periods. Further development of the model would include representation of domestic water uses in a stochastic manner.