The impact of climate change on rainwater harvesting in households in Poland

In water management, climate changes require adaptation, protection of existing resources and the search for alternative water sources. Rainwater harvesting (RWH) is increasingly becoming an alternative water source, applicable for many directions of its use. The aim of the research was to analyse the influence of long-term climate change on the potential for rainwater harvesting in households in Central Europe in Poland. The analysis of long-term climate changes impact on the household rainwater harvesting potential was conducted for the 50-year period for 19 cities in Poland. The water balance model was operated in a day-step mode, and the research was carried out in four “time scales approach”. For the purpose of all analyses, a standard weekly water demand profile was developed. It includes the daily sub-profiles for every working day and weekend, covering the washing and toilet flushing needs for a 4-person household. In order to evaluate the long-term changes in residential harvesting in Poland, a hypothetical residential RWH was investigated in 19 locations and rainfall conditions. To illustrate the time-spatial variation of annual rainfall amounts in the analyzed period, the heat map was prepared as the pre-simulation stage. Results show that the design of RWH systems should be based on archival data and take into account the many years of rainfall changes. This will improve performance and secure benefits for the users of this system.


Introduction
In recent years, more and more attention is paid to climate change, especially in the context of global warming and of the increasing occurrence of severe weather phenomena (Dai 2011, Kundzewicz 2012, Schiermeier 2011. An increase in the average annual temperature on the globe causes an increased water circulation in the hydrologic cycle, as the maximum water vapour mass in the atmosphere increases, as per the Clausius-Clapeyron relation, by approximately 7% along with a 1 °C increase in temperature (Held 2006, Lenderink 2011, Pall 2006). Therefore, one should expect that an increase in the availability of moisture in the atmosphere will cause a change in the rainfall regime, and this also includes an increase in the intensity of extreme rainfall (Allan 2008, Fleig 2015, Saboia 2017, Trenberth 2011. According to the European Environment, the annual precipitation trends have changed over the last century. Additionally, it depended on location. Northern Europe has become wetter (10-40% growth), and southern Europe has allow for maintaining a balance between water supply and demands (Czapczuk 2015, Adham 2016). However, due to the limited capacity of the tank, some quantities of rainwater would not be collected and this remaining water will be lost through the overflow. On the other hand, when a shortage of rainwater will occur, the remaining water demand will be covered by tap water. The simulation of reservoir operation with respect to time is the so-called behavioural analysis (Butler 2005).
Poland is one of the European countries threatened with a water deficit. The area of Central and Eastern Europe, including Poland and Spain, has on average less than 200 mm of fresh water available each year, while demand for water is 3-10 times higher. It is important that 97% of water in Poland comes from precipitation. Therefore, water resources are characterized by high variability in the annual cycle (in the summer season, the precipitation can be several times higher than in winter), between individual years and areas of the country. The southern areas of Poland (upland and mountainous areas) are characterized by a greater number of them, while the central and northern part of Poland is suffering from their deficit. In recent years, the problem of water deficit has been intensified both by increasing anthropopressure and climate changes.
Rainfall variability in Poland has been described, among others, in the work (Czarnecka 2012), in which there was an analysis of the monthly precipitation value from 38 weather stations (excluding mountain stations) from the years 1951-2010. Slightly increasing tendencies in the spring and autumn seasons and decreasing tendencies in the summer and winter seasons have been indicated in the majority of the country (however, these were statistically insignificant). A decrease in the total precipitation values in the summer months has also been confirmed in the work (Skowera 2014), on the basis of research from 16 stations from the years 1971-2010 and in the work (Żmudzka 2009) on the basis of research from 50 stations from the years 1951-2000. In monograph (Pińskwar 2010), there has been an analysis of the value of precipitation from 22 stations with daily data from the years 1951-2006 and from 6 stations with monthly data from the years 1956-2009. On the basis of the analysis of the precipitation, one has noted a decrease in the total amount of summer precipitation in comparison with the winter precipitation, as well as of the total values of precipitation in the warm half-year in relation to the sum of precipitation in the preceding cooler half-year (especially in the western parts of Poland). In the case of the total value of the annual precipitation, one can observe trends of both signs.
The goal of the work is an analysis of the influence of long-term climate change (the observed changes in the precipitation values) on the potential rainwater harvesting for washing and toilet flushing needs for a 4-person household in Poland. For the purpose of all analyses, a standard weekly water demand profile was developed. The analysis of longterm climate changes impact on the household rainwater harvesting potential was conducted for the 50-year period for 19 cities in Poland.

Materials and methods
A wide range of archival rainfall data obtained from the Institute of Meteorology and Water Management-National Research Institute (IMGW-PIB) was used as input data for the model. The foundation for the analysis was the implementation of the water balance model based on the behavioural algorithm to simulate the flows of water through the harvesting system. For all analyzed locations, the simulation input database was prepared based on historical precipitation data from 1970-2019. The water balance model was processing the hypothetical RWH system, having regard to the 50-year daily rainfall data, runoff capture, rainwater storage, and residential water demand. Benchmarking metrics for assessing the long-term performance and behaviour of all 19 RWH systems were represented by the simulation output data.
The water balance model was operated in a day-step mode, as described in the diagram in Fig. 1. The developed algorithm was carried out for all 19 locations (marked as i = 1…19) and 50 years (marked as j = 1970…2019). In practice, two algorithms are often used to describe the behaviour of the collecting system during a given time interval: the yield after spillage (YAS) and the yield before spillage (YBS). In this investigation, each simulation was tested by YAS algorithm, in which the withdrawal occurs before the rainfall: where Q o (m 3 ) is the volume discharged as overflow from the storage tank at time step t , V t (m 3 ) is the volume stored at time step t , R t (m) is the rainfall at time t , A (m 2 ) is impervious surface area (m 2 ), Ψ-surface runoff coefficient, Y t (m 3 ) is the yield of rainwater from the storage tank at time step t , D t (m 3 ) is the water demand at time step t , and S (m 3 ) is the tank storage capacity.
The usual indicator used to evaluate the performance of a rainwater collecting system is the water saving efficiency ( E ). It is equal to the percentage of the overall demand supplied by the tank and it is given by the following equation (Gires 2009): where T is the total time period under consideration. Equation provides a measure of how much water has been conserved in comparison with the overall demand and is also referred to as water saving efficiency (Fewkes 2000).
In the analyzed RWH model, the rainwater will be collected from the rooftop area A = 100 m 2 , with surface runoff coefficient Ψ = 0.95, stored in a typical rainwater tank with S = 1000 dm 3 and locally used. The developed water balance model provides the RWH metrics in monthly and annual The diagram of the water balance model for rainwater harvesting system scale. The basic calculation step is 24 h. The water balance model uses a daily sum of precipitation from historical data and a daily water demand profile. On this basis monthly and annual metrics of rainfall and operation of RWH system are developed. The statisticaelaboration allows to determine and assess significant trends of changes 1970-2019 according to climate changes in Poland and to draw appropriate conclusions on a multiannual, annual, monthly and daily basis.
For the purpose of all analyses, a standard weekly water demand profile was developed. It includes the daily subprofiles for every working day and weekend, covering the washing and toilet flushing needs for a 4-person household (Fig. 2). Each member of the family flushes a 3 or 6 dm 3 (on the average 4.5 dm 3 ) 4 times a day on a working day and 6 times on weekends. The laundry is washed four times a week, with a water demand of 40 dm 3 per cycle (Struk-Sokołowska 2020, Willis 2010, Freni 2019. The sum of water needs characterizes the water demand profile in litres per day (dm 3 /d) used in the long-term simulations. The annual water demand profile was imported into the water balance model as input data.
In order to evaluate the long-term changes in residential harvesting in Poland, a hypothetical residential RWH was investigated in 19 locations and rainfall conditions (Fig. 3). Precipitation is highly dependent on the shape of the surface. The average annual precipitation (from 1970-2019) in the analysed cities is 597.3 mm, with the lowest precipitation generally occurring in the central part of Poland and

POLAND
Koszalin Suwałki the highest on the coast and in the Tatras (over 1000 mm). Among the towns analysed in the study, the lowest precipitation occurs in Gdańsk (512.7 mm) and Poznań (518.1 mm), and the highest in Koszalin (731.6 mm) and Katowice (724.6 mm). To illustrate the time-spatial variation of annual rainfall amounts in the analyzed period (1970-2019) for all 19 cities, the heat map was prepared as the pre-simulation stage (Fig. 4). The preparation of the archive data was necessary as input data for the model used in the simulations. The analysis of the colour distribution of annual rainfall amounts shows the history of rainfall in Polish cities. Both dry and most humid years are visible.
The colour scheme indicates the dry year 1982 and the most humid year 2010, with average annual precipitation amounted to 424 mm and 791 mm, respectively. The longterm values of average annual precipitation ranged from 513 mm/year in Gdańsk to 732 mm/year in Koszalin, with a median of all 19 stations 598 mm/year. Absolutely the lowest precipitation amounts 275 mm were recorded in Poznan in 1982, whereas absolutely the highest precipitation amounts 1021 mm were recorded in Krakow in 2010. Many extreme meteorological and hydrological events took place in Poland at that time, and they were taken into account in daily water balance modelling. Meteorological drought and hydrological drought in Poland occur cyclically in periods of 5-10 years. All these factors describe the working conditions of RWH affecting its behaviour and performance. They are included in the prepared simulation.
The next step of the investigation was the preparation and analysis of the distribution of monthly precipitation levels  recorded in the Polish cities under consideration. An average annual sum of precipitation was also determined for all cities. The results are presented in Table 1.
In the case of monthly sums of precipitation, the lowest values were recorded in the winter months (an average of 40.9 mm in December, 35.6 mm in January and 28.2 mm in February) and the highest in the summer months (69.3 mm in June, 85.1 mm in July and 65.7 mm in August). This distribution is characteristic of Poland's climatic zone and occurs in all cities. It should be noted that the difference in the mid-month total precipitation between the towns with the largest and the smallest sums of precipitation is about double, i.e. in January from 25.6 mm in Gdańsk to 48.5 mm in Koszalin; in February-May from 18.0 mm, 23.1 mm, 28.6 mm and 48.6 mm in Gdańsk to 37.0 mm, 43.8 mm, 49.5 mm and 76.7 mm in Katowice; in June from 57.2 mm in Zielona Góra to 85.1 mm in Kraków; in July from 69.8 mm in Szczecin to 101.7 mm in Katowice; in August from 53.2 mm in Gorzów Wielkopolski to 80.1 mm in Koszalin; in September-November from 41.3 mm, 35.7 mm and 34.8 mm in Poznań to 76.7 mm, 67.6 mm and 63.2 mm in Koszalin, respectively; and in December from 31.9 mm in Wrocław to 58.7 mm in Koszalin.
The differential and distribution of monthly and annual sums of precipitation directly affect the performance RWH, differentiating them in particular localities and years. It should be noted that total precipitation is not the only factor determining the precipitation-related water saving efficiency.
The high amount of precipitation, which is characterised by a large amount of precipitation in a short period of time, will usually be used only to a small extent-the greater part of it will be transferred out of the system. Therefore, a proportionally smaller water saving efficiency in relation to the sum of precipitation is to be expected in the months in which the most frequent precipitation occurs in Poland, i.e. June-August (VI-VIII) (Kaźmierczak 2015, Pińskwar 2019. Due to the hourly simulation step, this phenomenon is correctly taken into account by the proposed water balance model. Simulations covering the period of 50 years in 19 Polish cities, taking into account the mentioned precipitation conditions, provided results determining the performance of hypothetical RWH. In order to determine trends of changes in water saving efficiency, linear regression was used (Wdowikowski 2016). Changes (increases or decreases) at a significance level above 95% were considered as statistically relevant. Moreover, changes at the significance level from 90 to 95% were assumed to be close to statistical significance, while changes in the significance level from 75 to 90% were assumed to be a tendency to change. Changes at the level of significance below 75% were considered statistically insignificant and consequently without a specific direction of change (Pińskwar 2010).

Water saving efficiency
The simulation results were graphically presented in the form of box plots, in which the height of the box corresponds to the value of the quarter range, i.e. the difference between the first and third quartiles (50% of all observations). The second quartile (median) was presented as the horizontal line inside the box. The whiskers are limited to the minimum and maximum values. Figure 5a shows the annual water saving efficiency as a coverage of analysed annual water demand for flushing toilets and washing expressed-for all 19 locations in the investigated multi-year period. The average degree of demand coverage at the national level amounted to favourably as much as 80.7% (minimum 74.4%, maximum 85.1%), whereby the differences between cities were usually minor. While analysing individual cities, the highest average coverage of water demands was recorded in Koszalin-average equal to 85.1% (minimum value 65.2% in 1982, maximum value 95.2% in 1998), while the lowest in Gdańsk-average 74.4% (minimum value 53.3% in 1982, maximum value 94.7% in 1970). The highest annual coverage of needs equal to 97.9% (with annual rainfall of 857.4 mm) was recorded in Kielce in 1970, with the lowest one in already mentioned Gdańsk (53.3%-with annual rainfall of 345.5 mm). The greatest diversity of annual performances was recorded in Gdańsk, where the minimum value is almost half the maximum one. This implies large discrepancies in tap water savings in individual years of RWH operation. The highest median and short box oscillating around high values is favourable for achieving the highest water saving efficiency of the RWH system. The whiskers are not decisive, as reflecting the recorded outliers-maximum and minimum values. Figure 5b illustrates the dependence of the annual RWH performance on the total annual precipitation in all simulated years and cities. The relationship is very clear, but the height and width of the cloud indicate that, as expected, the height of annual precipitation is not the only factor that affects water saving efficiency. For example, in 1993 and 1996, Warsaw had an identical annual rainfall of 453.6 mm (green and red dot in the figure), and the water saving efficiency was as high as 83.7% and only 66.9%, respectively (with a multi-year average of 77.5% at 538 mm rainfall). The course of the 1993 cumulative precipitation curve in Warsaw (Fig. 6a) is much more linear than in 1996 (Fig. 6b), when many days with precipitation exceeding 10 mm were recorded. This translated into an average water saving efficiency per 100 mm of precipitation −18.5% and 14.7%, respectively.
The above analysis proves that the best performance from the point of view of RWH is years with not only a large amount of precipitation but also its successive growth during the year. Extreme cases in this respect were recorded in Toruń in 1989 and in Katowice in 1974, where in the first case, 310.4 mm of precipitation gave 70.5% water saving efficiency (on average 100 mm of precipitation covered as much as 22.7% of annual demand) and in the second case, 1011.2 mm of precipitation gave 83.6% water saving efficiency (on average 100 mm of precipitation covered only 8.3% of demand).
The results of calculations of the average water saving efficiency on an annual and monthly basis for 19 analysed cities in the period 1970-2019 are presented in Table 2.
Same as the Fig. 5b, Fig. 7 presents the water saving efficiency in individual months of the year depending on the sum of precipitation recorded in the years 1970-2019 in 19 analysed localities. The graphical presentation of the results shows both the varied nature of precipitation in individual months described by the Polish climate and the difference in the water saving efficiency obtained for the same monthly sums of precipitation. For the analyzed hypothetical RWH system, the same in all locations, for the same average monthly precipitation, differentiation of harvesting performance was obtained. This means that, similarly to the annual sum of precipitation, the sum of precipitation in given months is not the only factor determining the water saving efficiency (although this dependence is also clear at this point).
Analysing the above results, it must be concluded that in order to determine the impact of climate changes on water saving efficiency, it is not possible to examine only trends in annual or monthly amounts of precipitation. The dependence of the water saving efficiency on the amount of precipitation, although clear, cannot be determined unequivocally, as an equally important factor as the sum of precipitation is its course over time (e.g. Fig. 6). Therefore, in order to analyse the impact of climate change on the water saving efficiency, it is not necessary to examine the variability of precipitation totals, but the variability of the water saving efficiency of the modelled system (with assumptions).

The trends of change
Based on the simulated efficiency of the RWH system, its long-term variability in the context of climate change was determined. Calculated on the basis of archival data of RHW efficiency from 1970 to 2019 are shown in Table 2. The variation of water saving efficiency was calculated by means of linear regression and compiled in Table 3 as ΔE in % per decade. Statistically significant and close to statistically significant trends with a significance level above 90% were distinguished by colours. It should be noted that positive (blue) and negative (yellow) trends appeared in the table. The results were compared in two "time scales" as multiyear trends in monthly and annual precipitation changes in all 19 analysed Polish cities.
In the case of only one of the analysed locations, the annual water saving efficiency is characterised by a significant trend. In Opole, a decrease in −1.44% per decade was recorded. For this city in Fig. 8 presents the values and trend of annual water saving efficiency in the years 1970-2019 against the background of annual rainfall.
For the other cities, the changes in annual performance are statistically insignificant. This means that in the face of At the level of analysis of monthly water saving efficiency averages the situation turns out to be more dynamic. In all analysed cities, there are positive and negative trends in monthly sums of precipitation, including as many as 24 statistically significant ones in the analysed period. Statistically significant changes in monthly efficiency appear in cities where no statistically significant change in annual efficiency was recorded. This interesting phenomenon means significant trends in the annual precipitation distribution.
The distribution of statistically significant trends (positive blue, negative yellow) in Table 3 shows no statistically significant annual changes (except for Opole, where a statistically significant decreasing trend was noted) and a statistically significant monthly increase in winter and decrease in summer. This confirms the observation of other authors from On the basis of the data from Table 3, graphs for statistically significant change trends have been developed. Figure 9 presents 24 statistically significant trends of changes in monthly efficiency compared to monthly precipitation totals (also with a trend line-blue). Figure 9 shows evident trends in monthly water saving efficiency over the 50 years analysed from 1970 to 2019. Most of the trends in water saving efficiency correlate with changes in precipitation, although there are exceptions to this rule. Statistically significant trends of changes in water saving efficiency in Białystok, Kielce and Toruń in September are negative, although trends of changes in monthly total precipitation are positive in these locations in September. This clearly confirms that total precipitation is not the only factor determining water saving efficiency.
The use of valuable long-term data from 1970 to 2019 on precipitation provided full information on long-term changes in the potential of the RWH in 19 large Polish cities. Time scale approach in water balance modelling results investigation provided interesting and unexpected conclusions about the impact of perennial climate change on rainwater harvesting and trends in RWH operation and performance.

Conclusions
On the basis of the results of the research carried out the following conclusions were drawn: • In the context of undisputed climate changes, statistically significant trends of changes in the average annual efficiency of the RWH system in 18 out of 19 analysed cities in Poland in the years 1970-2019 were not found. A statistically significant impact of climate changes on the annual water saving efficiency was found only in Opole. A decrease in efficiency by 1.44% per decade means clearly decreasing economic and ecological benefits for users of such a system during its operation. This is a negative trend. • The analysis on a monthly basis provided more detailed information. In the analysed half of the century, statistically significant changes in the monthly efficiency of the RWH system were shown. A total of 24 significant monthly trends show an increase in efficiency in the winter months and a decrease in the summer months. For the residential RWH systems, this is a positive change reducing the seasonal difference between summer and winter precipitation. For garden watering RWH systems, it is an unfavourable change due to the reduced availability of water for watering at the beginning of the growing season. Research on water collected from roofs should be extended to include quantitative and qualitative analysis of its composition, including identification of heavy metals. Their results will confirm the suitability or disqualification of water collected in particular regions (especially highly industrialised ones) for domestic use and garden watering. • Differences in statistical significance of annual and monthly trends of RWH efficiency in the analysed cities indicate that the sum of precipitation is not the only indicator influencing its efficiency. The course of precipitation is also important. Therefore, water saving efficiency trends should be analysed, not precipitation trends. • The design of RWH systems should be based on archival data and take into account the many years of rainfall changes. This will improve performance and secure benefits for the users of this system.

Conflict of interest
The authors declare that they have no conflict of interest.
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/.