Assessment of Causes and E ﬀ ects of Groundwater Level Change in an Urban Area (Warsaw, Poland)

: Monitoring the data of groundwater level in long-term measurement series has allowed for assessment of the impact of natural and anthropogenic factors on groundwater recharge. It allows for assessing the actual groundwater quantity, which constitutes the basis for balanced and sustainable groundwater planning and management in an urban area. Groundwater levels in three aquifers were studied: the shallow and deeper Quaternary aquifers and the Oligocene aquifer in Warsaw (Poland). Statistical analysis was performed on a 27-year (1993–2019) cycle of daily measurements of groundwater levels. The studies focused on determining the range and causes of groundwater level changes in urban-area aquifers. The groundwater table position in the Quaternary aquifer pointed to variable long-term recharge and allowed for the identiﬁcation of homogenous intervals with identiﬁcation of water table ﬂuctuation trends. A decrease in the water table was observed within the Quaternary aquifers. The Oligocene aquifer displayed an opposite trend.


Location of the Study Area
The groundwater monitoring system is located in the center of Warsaw (Poland) in the Research Station of the Faculty of Geology at the University of Warsaw. The research station is located on the crossroads of two busy streets,Żwirki i Wigury and Banacha, in the Ochota district. The district is characterized by the prevalence of multifamily buildings, similarly to most of left-bank Warsaw. In turn, single-family houses dominate right-bank Warsaw. Larger areas of undeveloped zones (park and forest complexes) are found in the urban area peripheries (Figure 1). Warsaw is located on two sides of the Vistula River and is slightly elongated along its banks (extending ca. 30 km along the N-S direction and at ca. 28 km along the W-E direction). The urban area covers an area of 517.24 km 2 . The population is 1.8 million at a density of 3500 people/km 2 . The city is located on two sides of the Vistula River, the largest river in Poland. Two types of terraces can be distinguished along the river: a modern floodplain covering a larger surface area in the southern part of the city and older terraces with the largest surface area in the northern part of right-bank Warsaw. Moraine plateaus extend from the river to the west as the Warsaw Plain and to the east as the Wołomin Plain. The research station is located on the Warsaw Plain ( Figure 1).
Groundwater recharge and drainage are variable within the urban area. In the centre of Warsaw, infiltration and its drainage are restricted by spatial management and a sewage system (stormwater drainage) typical for urban infrastructure. Precipitation from surfaces covered by non-permeable materials (roads, pavements, and roofs) is distributed by the sewage system, and thus effective infiltration is completely reduced. Groundwater recharge takes place through lateral inflow from the neighboring areas and to a large degree through water loss from the water-sewage system of the urban area [21,22]. Infiltration has a significant contribution to groundwater recharge and resource formation in other parts of the urban area.

Geology and Hydrogeological Conditions
The research station is located within a denuded post-glacial plateau located at an elevation of 108-115 m a.s.l. and incised by a glacial valley. The oldest series drilled in the Warsaw region are Upper Cretaceous strata, developed as grey marls interbedded with limestones and covered by Oligocene strata developed as sands, silts, and clays with glauconite and phosphorites. Miocene deposits include sands, silts, and clays with interbeds of lignite. The Pliocene comprises variably colored compact clays, the so-called variegated clays, and clay and sandy silts. The Pleistocene includes glacial, fluvioglacial, and ice-dammed deposits of the Podlasian, Sanian, Odranian, and Vistulian glaciations, and the Cromerian, Mazovian, and Eemian interglacials. In the vicinity of the research station, these deposits are developed as glacial tills located on silts, sands, and ice-dammed clays ( Figure 2). The geological succession of the subsurface zone comprises (from the top) anthropogenic embankment deposits, sandy silts, clayey sands, glacial tills, and fine-grained and silty sands (Figure 2). At a depth of about 7.6 m b.g.l. within fine-grained and silty sands is the first Quaternary aquifer (first aquifer) with a water table monitored by piezometers. The deeper Quaternary aquifer (second aquifer) with a confined water table is located within fine-grained and poorly sorted sands at a depth interval from 26 to 46 m b.g.l. The water table drilled during the well construction at the depth of 18 m b.g.l. became stabilized at 9.95 m b.g.l. (Figure 2).
Within the Warsaw agglomeration, waters of the Quaternary aquifer remain in hydraulic contact and are combined in one aquifer in the northern part of the urban area. Generally, groundwater in Quaternary deposits with hydraulic conductivity in the range of 25-30 m/day, whose thickness increases northwards and flows from the east with a velocity of 30-100 m/year [23]. The city is located on two sides of the Vistula River, the largest river in Poland. Two types of terraces can be distinguished along the river: a modern floodplain covering a larger surface area in the southern part of the city and older terraces with the largest surface area in the northern part of right-bank Warsaw. Moraine plateaus extend from the river to the west as the Warsaw Plain and to the east as the Wołomin Plain. The research station is located on the Warsaw Plain ( Figure 1).
Groundwater recharge and drainage are variable within the urban area. In the centre of Warsaw, infiltration and its drainage are restricted by spatial management and a sewage system (stormwater drainage) typical for urban infrastructure. Precipitation from surfaces covered by non-permeable materials (roads, pavements, and roofs) is distributed by the sewage system, and thus effective infiltration is completely reduced. Groundwater recharge takes place through lateral inflow from the neighboring areas and to a large degree through water loss from the water-sewage system of the urban area [21,22]. Infiltration has a significant contribution to groundwater recharge and resource formation in other parts of the urban area.

Geology and Hydrogeological Conditions
The research station is located within a denuded post-glacial plateau located at an elevation of 108-115 m a.s.l. and incised by a glacial valley. The oldest series drilled in the Warsaw region are Upper Cretaceous strata, developed as grey marls interbedded with limestones and covered by Oligocene strata developed as sands, silts, and clays with glauconite and phosphorites. Miocene deposits include sands, silts, and clays with interbeds of lignite. The Pliocene comprises variably colored compact clays, the so-called variegated clays, and clay and sandy silts. The Pleistocene includes glacial, fluvioglacial, and ice-dammed deposits of the Podlasian, Sanian, Odranian, and Vistulian glaciations, and the Cromerian, Mazovian, and Eemian interglacials. In the vicinity of the research station, these deposits are developed as glacial tills located on silts, sands, and ice-dammed clays ( Figure 2). The geological succession of the subsurface zone comprises (from the top) anthropogenic embankment deposits, sandy silts, clayey sands, glacial tills, and fine-grained and silty sands (Figure 2). At a depth of about 7.6 m b.g.l. within fine-grained and silty sands is the first Quaternary aquifer (first aquifer) with a water table monitored by piezometers. The deeper Quaternary aquifer (second aquifer) with a confined water table is located within fine-grained and poorly sorted sands at a depth interval from 26 to 46 m b.g.l. The water table drilled during the well construction at the depth of 18 m b.g.l. became stabilized at 9.95 m b.g.l. (Figure 2). Below the Quaternary aquifers is a 100-m thick series of Pliocene clays and Miocene silty sands, silts, and clays with lignite interbeds, with variable thicknesses. The deepest Oligocene aquifer (third aquifer) recognized in Warsaw in the monitoring system is found in fine-and medium-grained sands with glauconite at a depth interval from 225 to 263 m b.g.l. The groundwater table drilled at a depth of 221 m b.g.l. became stabilized at the depth of 35.75 m b.g.l. during the drilling.

Water Supply in Warsaw
The Warsaw agglomeration is supplied about 96% of its water from surface water sources, and the remaining 4% is from groundwater sources. Surface water comes from the Vistula River and the Zegrze Reservoir (located 7 km north of the Warsaw border). The water system, with a total length of 4215.7 km, covers almost the entire area of the Warsaw urban agglomeration. The population of Warsaw uses about 340,000 m 3 of water per day [24], and the daily use of water is about 135 L/person [25]. Groundwater is exploited mainly by industrial plants. About 500 intakes from the Quaternary aquifers and about 100 intakes from the Oligocene aquifer, exploiting almost 25,000 m 3 /day, are located in the Warsaw agglomeration. Due to its quantity and quality, groundwater represents the strategic groundwater resources of the urban area. In the entire Mazovian voivodeship, of which Within the Warsaw agglomeration, waters of the Quaternary aquifer remain in hydraulic contact and are combined in one aquifer in the northern part of the urban area. Generally, groundwater in Quaternary deposits with hydraulic conductivity in the range of 25-30 m/day, whose thickness increases northwards and flows from the east with a velocity of 30-100 m/year [23].
Below the Quaternary aquifers is a 100-m thick series of Pliocene clays and Miocene silty sands, silts, and clays with lignite interbeds, with variable thicknesses. The deepest Oligocene aquifer (third aquifer) recognized in Warsaw in the monitoring system is found in fine-and medium-grained sands with glauconite at a depth interval from 225 to 263 m b.g.l. The groundwater table drilled at a depth of 221 m b.g.l. became stabilized at the depth of 35.75 m b.g.l. during the drilling.

Water Supply in Warsaw
The Warsaw agglomeration is supplied about 96% of its water from surface water sources, and the remaining 4% is from groundwater sources. Surface water comes from the Vistula River and the Zegrze Reservoir (located 7 km north of the Warsaw border). The water system, with a total length of 4215.7 km, Water 2020, 12, 3107 5 of 17 covers almost the entire area of the Warsaw urban agglomeration. The population of Warsaw uses about 340,000 m 3 of water per day [24], and the daily use of water is about 135 L/person [25]. Groundwater is exploited mainly by industrial plants. About 500 intakes from the Quaternary aquifers and about 100 intakes from the Oligocene aquifer, exploiting almost 25,000 m 3 /day, are located in the Warsaw agglomeration. Due to its quantity and quality, groundwater represents the strategic groundwater resources of the urban area. In the entire Mazovian voivodeship, of which Warsaw is the largest urban area, the exploitable groundwater resources from Quaternary deposits are in the range of about 205,000 m 3 /h, and from the Palaeogene-Neogene deposits, they are in the range of 17,500 m 3 /h [26].

Materials and Methods
Changes in groundwater levels in Warsaw have been studied by means of the following (

Materials and Methods
Changes in groundwater levels in Warsaw have been studied by means of the following (  Measurements of the water table depth are made automatically with continuous control. The research station monitors three aquifers: two Quaternary aquifers (shallow and deeper) and the Oligocene aquifer (Figure 2), using 1 piezometer and 2 wells. The analysis was performed on a 27-year measurement series in hydrogeological years 1993-2019; therefore, the total number of daily observations was 9861 for each aquifer. In the frame of the groundwater monitoring system in Poland, the observation site of water levels is located in Warsaw (II-22-1, the western part of the city; Figure 1) in the deeper Quaternary aquifer. Measurements of groundwater levels have been Measurements of the water table depth are made automatically with continuous control. The research station monitors three aquifers: two Quaternary aquifers (shallow and deeper) and the Oligocene aquifer ( Figure 2), using 1 piezometer and 2 wells. The analysis was performed on a 27-year measurement series in hydrogeological years 1993-2019; therefore, the total number of daily observations was 9861 for each aquifer. In the frame of the groundwater monitoring system in Poland, the observation site of water levels is located in Warsaw (II-22-1, the western part of the city; Figure 1) in the deeper Quaternary aquifer. Measurements of groundwater levels have been performed at this site since 1994 [28].
Analysis of the changes of groundwater levels was performed using statistical methods most commonly applied in assessment studies of monitoring data [29][30][31][32][33][34]. The performed analysis allows the determination and characterization of statistical data and prediction of groundwater level changes. Seasonal variability of the water level results from cyclic climate changes. It was characterized by annual amplitude, median value and analysis of charts with water table fluctuations. The annual amplitude was calculated as the difference between the highest and lowest multiannual mean water level, or as the difference between the absolute highest and the absolute lowest registered water level. The coefficient of variation, i.e., ratio of the standard deviation to the mean value, was also calculated. The measure of multiannual change of groundwater levels is the value of standard deviation of mean annual groundwater levels [35].
The magnitudes of estimated changes in the trend of meteorological variables in this study were estimated by the application of the Sen slope method. This technique calculates the gradient as a change in the measurements correlated with units of temporal change. The advantages of this method include an allowance for missing data, avoidance of assumptions about the distribution of tested data, and averting the effects of gross data errors and outliers [36]. Therefore, the method reduces the consequences of missing data and/or anomalous trends therein using the median values of the time series of various slopes that were detected as an evaluation tool [36,37]: In this equation, Xj and Xi denote values dated at times j and i, and time j is after time i. The estimator β is the median overall combination of the recorded pairs for the entire dataset where trend analysis is performed, and a positive β indicates an increasing trend, while a negative value indicates a decreasing trend [37].
Assessment of the trend of changes in the water table level in the studied aquifers using the Sen slope method was based on the average monthly positions of the water table.
The magnitude of infiltration was determined based on measurements of groundwater levels. The methodology of assessing the magnitude of recharge at a regional scale is the focus of many hydrogeological investigations, and their synthesis has been presented e.g., by de Vries, Simmers [38], Scanlon et al. [39], and Assouline [40].
Observations of the groundwater levels of the shallow Quaternary aquifer (first aquifer) allowed the determination of the recharge using the water table fluctuation (WTF) method (e.g., [41,42]). The method is based on an assumption that the increase in groundwater level is caused by infiltration. In intervals with no recharge, the water level decreases. Recharge is calculated as follows: The method is mainly used for short observation periods in areas with a shallow groundwater level, which causes quite large changes in its position over time [39]. The mean annual sum of precipitation from the multiannual interval 1993-2019 in the Warsaw area at 570 mm was used in the calculations of the infiltration coefficient [43]. The meteorological data came from the Warszawa-Okęcie Station, located in the south of Warsaw. A similar average amount of rainfall (571 mm) was recorded in the years 2001-2010 in Warsaw [44]

Results
The depth of the water table of the shallow Quaternary aquifer was at intervals from 3.62 to 7.76 m and a mean depth of 5.91 m over the period of 1993-2019. The deeper Quaternary aquifer had a depth of 6.40 to 10.13 m and a mean depth of 8.21 m. In the Oligocene aquifer, the extreme values of the depth to the water table were from 21.75 to 32.28 m, at a mean depth of 25.44 m ( Table 1). The characteristic groundwater levels in three aquifers were determined in order to recognize the dynamics of changes indicating similarities and differences in the hydrodynamic conditions. The value of the 1st quartile (5.56) and 3rd quartile (6.28) was exceeded over 2440 times in the case of depth to groundwater of shallow Quaternary aquifer ( Table 1) The Oligocene aquifer is different compared to Quaternary aquifers in the case of the 1st and 3rd quartiles. The 1st quartile was exceeded from October 2013 to October 2019, while the 3rd quartile was exceeded from January 1993 to July 1999.
Taking into consideration the coefficient of variation, it can be seen that variability of depth of the groundwater table is statistically significant in the case of shallow Quaternary and Oligocene aquifers ( Table 1).
A similarity can be observed between the Quaternary aquifers with regard to low groundwater depths. The amplitudes (WNG-NNG; highest low groundwater depth to lowermost low groundwater depth) were similar in the first and second aquifers and reached a maximum of 3 m (Table 2). In both aquifers, there were similar differences between the mean depth to the water table in the multiannual interval and the mean depth at low levels (SNG; average low groundwater depth). In the case of the Oligocene aquifer, the amplitude (WNG-NNG) exceeded 10.5 m, whereas the difference between the mean depth to the water table in the multiannual interval and at low levels was 1 m. This indicator points to the similarity of hydrodynamic conditions in the Quaternary aquifers and the distinctiveness of the Oligocene aquifer. The ratio of annual mean to multiannual mean is 0.17 in the first aquifer, 0.20 in the second aquifer and 0.13 in the third aquifer, which indicates significant multiannual changes at small seasonal amplitudes.
A comparison of daily measurements from both Quaternary aquifers shows that larger amplitudes occur in the shallow aquifer (4.14 m) than in the deeper aquifer (3.73 m). Analysis of the water table position on an annual basis did not reveal cyclic changes related to seasonal fluctuations of the water table, with typical high levels in spring-summer and low levels in autumn-winter. The water table position was not affected by seasonal infiltration changes. There was no direct relationship between precipitation and groundwater levels of the Quaternary aquifers ( Figure 4). The analysis of trends of long-term changes carried out using the Sen's slope method indicates a decreasing trend of the water table in a shallow Quaternary aquifer (Table 3).  The ratio of annual mean to multiannual mean is 0.17 in the first aquifer, 0.20 in the second aquifer and 0.13 in the third aquifer, which indicates significant multiannual changes at small seasonal amplitudes.
A comparison of daily measurements from both Quaternary aquifers shows that larger amplitudes occur in the shallow aquifer (4.14 m) than in the deeper aquifer (3.73 m). Analysis of the water table position on an annual basis did not reveal cyclic changes related to seasonal fluctuations of the water table, with typical high levels in spring-summer and low levels in autumn-winter. The water table position was not affected by seasonal infiltration changes. There was no direct relationship between precipitation and groundwater levels of the Quaternary aquifers ( Figure 4). The analysis of trends of long-term changes carried out using the Sen's slope method indicates a decreasing trend of the water table in a shallow Quaternary aquifer (Table 3).    In the deeper Quaternary aquifer, the trend of changes of the water table position was opposite in the multiannual interval. However, detailed analysis of the groundwater levels in both aquifers shows that only in the first part of the multiannual interval (when the Quaternary aquifer was intensively exploited for industrial requirements), the water table of the deeper aquifer was at a lower level, indicating larger differences between both aquifers. After reducing the exploitation of groundwater due to the economic transformation in Poland in the early 1990s, since 1997 the water table of the deeper Quaternary aquifer began to increase systematically, and since 2001, its position has been analogous to the fluctuations in the shallow Quaternary aquifer. Similar dynamics of fluctuations in the multiannual interval 2001-2019 may indicate that after reducing exploitation in the deeper Quaternary aquifer, they are influenced by the same factors (recharge, drainage) that decide on the analogous water table position. The similarity between the water table positions of both Quaternary aquifers suggests the existence of a hydraulic connection between them in zones beyond the area of the research station ( Figure 5).
Water 2020, 12, x FOR PEER REVIEW 9 of 17 In the deeper Quaternary aquifer, the trend of changes of the water table position was opposite in the multiannual interval. However, detailed analysis of the groundwater levels in both aquifers shows that only in the first part of the multiannual interval (when the Quaternary aquifer was intensively exploited for industrial requirements), the water table of the deeper aquifer was at a lower level, indicating larger differences between both aquifers. After reducing the exploitation of groundwater due to the economic transformation in Poland in the early 1990s, since 1997 the water table of the deeper Quaternary aquifer began to increase systematically, and since 2001, its position has been analogous to the fluctuations in the shallow Quaternary aquifer. Similar dynamics of fluctuations in the multiannual interval 2001-2019 may indicate that after reducing exploitation in the deeper Quaternary aquifer, they are influenced by the same factors (recharge, drainage) that decide on the analogous water table position. The similarity between the water table positions of both Quaternary aquifers suggests the existence of a hydraulic connection between them in zones beyond the area of the research station ( Figure 5). Assessment of the reliability of accumulated data allowing for spatial interpretations in the Warsaw city was achieved by analysing changes of water table positions of the deeper Quaternary aquifer in well no. II-22-1 [28], located in the westernmost part of the city ( Figure 5). Comparison of the water table fluctuations in both Quaternary wells indicates their clear similarity, which allows the local observations in the research station to be transferred to spatial assessments (Figures 4 and 5, Table 3).
The potentiometric surface of the Oligocene aquifer displays an opposite trend compared to both Quaternary aquifers. The trend of changes in the potentiometric surface of the Oligocene aquifer in the multiannual interval 1993-2019 points to the gradual restoration of hydrostatic pressure. In 1993, the potentiometric surface became stabilized at 80.83 m a.s.l., whereas it presently stabilizes at 90.72 m a.s.l., which shows a rise of the potentiometric surface in the multiannual interval by about 10 m ( Figure 6). Assessment of the reliability of accumulated data allowing for spatial interpretations in the Warsaw city was achieved by analysing changes of water table positions of the deeper Quaternary aquifer in well no. II-22-1 [28], located in the westernmost part of the city ( Figure 5). Comparison of the water table fluctuations in both Quaternary wells indicates their clear similarity, which allows the local observations in the research station to be transferred to spatial assessments (Figures 4 and 5, Table 3).
The potentiometric surface of the Oligocene aquifer displays an opposite trend compared to both Quaternary aquifers. The trend of changes in the potentiometric surface of the Oligocene aquifer in the multiannual interval 1993-2019 points to the gradual restoration of hydrostatic pressure. In 1993, the potentiometric surface became stabilized at 80.83 m a.s.l., whereas it presently stabilizes at 90.72 m a.s.l., which shows a rise of the potentiometric surface in the multiannual interval by about 10 m ( Figure 6). Trend analysis of the water table levels using Sen's slope method confirms the rising trend of the Oligocene potentiometric surface ( Table 3).
The gradual restoration of pressure in the Oligocene aquifer results from reduced exploitation. The largest increase in the potentiometric surface was observed in 1997-2003 (Figures 6 and 7). After 2003, the potentiometric surface position changed in a much narrower range. There was excessive exploitation of groundwater from the Oligocene aquifer in Warsaw in the 1970s, reaching 50,000-60,000 m 3 /d, resulting in the formation of a vast depression cone. Following the protection of these waters, in 1977, the construction of new water intakes for industrial purposes was banned, and the intake by existing industrial plants was strongly limited [45,46]. At present, 53 public sources rely on groundwater from Oligocene aquifer in the Warsaw city (as of 13.12.2019) [47]. Water from this aquifer is also used by hospitals and industries requiring water of the highest quality. The introduced restrictions have brought expected effects. A similar case took place in Barcelona, in which the intensive extraction of water from deltaic aquifers caused large drawdowns, leading to seawater intrusion [48]. Since the 1970s, many industries have migrated from the city, which reduced water demand, resulting in a progressive recovery of heads.  Trend analysis of the water table levels using Sen's slope method confirms the rising trend of the Oligocene potentiometric surface ( Table 3).
The gradual restoration of pressure in the Oligocene aquifer results from reduced exploitation. The largest increase in the potentiometric surface was observed in 1997-2003 (Figures 6 and 7). After 2003, the potentiometric surface position changed in a much narrower range. There was excessive exploitation of groundwater from the Oligocene aquifer in Warsaw in the 1970s, reaching 50,000-60,000 m 3 /d, resulting in the formation of a vast depression cone. Following the protection of these waters, in 1977, the construction of new water intakes for industrial purposes was banned, and the intake by existing industrial plants was strongly limited [45,46]. At present, 53 public sources rely on groundwater from Oligocene aquifer in the Warsaw city (as of 13.12.2019) [47]. Water from this aquifer is also used by hospitals and industries requiring water of the highest quality. The introduced restrictions have brought expected effects. A similar case took place in Barcelona, in which the intensive extraction of water from deltaic aquifers caused large drawdowns, leading to seawater intrusion [48]. Since the 1970s, many industries have migrated from the city, which reduced water demand, resulting in a progressive recovery of heads.
Accepting the urban area as the watershed-scale balance area and undertaking the attempt to identify the causes of groundwater level changes on a short-term and long-term basis, the basin water budget was evaluated for the shallow Quaternary aquifer. Recharge was assessed based on the analysis of changes of the water table position monitored in the first aquifer using the water table fluctuation (WTF) method. The specific yield taken for the calculations was at 0.10 [41,[49][50][51]. This value is typical for sediments with low hydraulic conductivity (the aeration zone comprises sands with the addition of the silt fraction and clay).
Due to the lack of seasonal water table fluctuations, annual recharge was not discussed; instead, the method was modified by assessing recharge in intervals similar with regard to the dynamics of the water table. Based on the analysis of the water table position, four intervals were distinguished, differing in water table fluctuations as reflected in the magnitude of recharge (Figure 8).
the intake by existing industrial plants was strongly limited [45,46]. At present, 53 public sources rely on groundwater from Oligocene aquifer in the Warsaw city (as of 13.12.2019) [47]. Water from this aquifer is also used by hospitals and industries requiring water of the highest quality. The introduced restrictions have brought expected effects. A similar case took place in Barcelona, in which the intensive extraction of water from deltaic aquifers caused large drawdowns, leading to seawater intrusion [48]. Since the 1970s, many industries have migrated from the city, which reduced water demand, resulting in a progressive recovery of heads.   Accepting the urban area as the watershed-scale balance area and undertaking the attempt to identify the causes of groundwater level changes on a short-term and long-term basis, the basin water budget was evaluated for the shallow Quaternary aquifer. Recharge was assessed based on the analysis of changes of the water table position monitored in the first aquifer using the water table fluctuation (WTF) method. The specific yield taken for the calculations was at 0.10 [41,[49][50][51]. This value is typical for sediments with low hydraulic conductivity (the aeration zone comprises sands with the addition of the silt fraction and clay).
Due to the lack of seasonal water table fluctuations, annual recharge was not discussed; instead, the method was modified by assessing recharge in intervals similar with regard to the dynamics of the water table. Based on the analysis of the water table position, four intervals were distinguished, differing in water table fluctuations as reflected in the magnitude of recharge (Figure 8). In the first interval, in the hydrological years 1993-1994, and in the third interval, in 2010-2011, infiltration was a long-term process (from 2 to 9 months), which caused a significant rise of the water table (1.2-2.9 m). The recharge coefficient was calculated as the percentage contribution of recharge (R) in the mean precipitation (P), assuming 570 mm/year of the precipitation for Warsaw. In these intervals, the annual recharge was from 130 to 160 mm, at a recharge coefficient from 23 to 28% (Table 4). In the second interval, in 1995-2009, the infiltration time was much shorter (1 to 2 months) and caused smaller seasonal fluctuations of the water table. Annual recharge was about 80 mm, and the recharge coefficient was 14%. For 16 years, the water table fell by almost 2 m. The fourth interval, in 2012-2019, was characterized by slightly lower values of parameters; infiltration was short-term (5-20 days), and despite a few such episodes during the year, systematic fall of the water table was observed. In 2019, the water table was at 7.71 m, i.e., 1.8 m below the mean value for the entire multiannual interval 1993-2019. In the first aquifer, only long-term recharge lasting several months can result in changes in the water table position. The observed differences in the water table position and its gradual decrease from 2012 are not caused entirely by climate change as reflected in precipitation recharge. The distribution of annual totals of precipitation did not differ significantly from the remaining intervals during this period (Figure 4). This factor is not the only reason for long-term and non-cyclic changes. The lack of typical seasonality is the cumulative effect of many elements, which are interdependent and cannot be analyzed separately. In the first interval, in the hydrological years 1993-1994, and in the third interval, in 2010-2011, infiltration was a long-term process (from 2 to 9 months), which caused a significant rise of the water table (1.2-2.9 m). The recharge coefficient was calculated as the percentage contribution of recharge (R) in the mean precipitation (P), assuming 570 mm/year of the precipitation for Warsaw. In these intervals, the annual recharge was from 130 to 160 mm, at a recharge coefficient from 23 to 28% (Table 4). In the second interval, in 1995-2009, the infiltration time was much shorter (1 to 2 months) and caused smaller seasonal fluctuations of the water table. Annual recharge was about 80 mm, and the recharge coefficient was 14%. For 16 years, the water table fell by almost 2 m. The fourth interval, in 2012-2019, was characterized by slightly lower values of parameters; infiltration was short-term (5-20 days), and despite a few such episodes during the year, systematic fall of the water table was observed. In 2019, the water table was at 7.71 m, i.e., 1.8 m below the mean value for the entire multiannual interval 1993-2019. In the first aquifer, only long-term recharge lasting several months can result in changes in the water table position. The observed differences in the water table position and its gradual decrease from 2012 are not caused entirely by climate change as reflected in precipitation recharge. The distribution of annual totals of precipitation did not differ significantly from the remaining intervals during this period ( Figure 4). This factor is not the only reason for long-term and non-cyclic changes. The lack of typical seasonality is the cumulative effect of many elements, which are interdependent and cannot be analyzed separately.

Discussion
Recharge in urban areas depends on a number of factors, including climate conditions and the dynamics of their changes [52,53], lithology of subsurface sediments [54,55], topography [56], land use and land cover [52,57], and thickness of the aeration zone [55,58].
The change of climate conditions in Warsaw was recorded as anomalously high precipitation in 2017 (821 mm) and anomalously low precipitation in 2015 (297 mm) (Figure 4). Measurements of groundwater levels in the research station have indicated that this change has no direct impact on the change of groundwater levels, even in the case of the shallowest Quaternary aquifer. In all analyzed aquifers, the influence of precipitation on water table fluctuations was insignificant due to the rather unfavourable filtration parameters in the aeration zone and the impact of urban infrastructure, especially in the subsurface zone. The reaction of the water table, even after intense precipitation, was greatly delayed. Multiannual trends of changes in the position of the water table in Quaternary aquifers and a stable trend of increasing potentiometric surface in the Oligocene aquifer are notable.
In urban areas, the magnitude of recharge depends on land use. The research station is located in Warsaw's centre and was developed at a time when there was hardly any infrastructure around it. However, in the 27-year interval, numerous changes in its surroundings have taken place, particularly related to the development of the University of Warsaw, Ochota Campus infrastructure, which may have had a significant impact on groundwater recharge and flow in the shallowest zone. The recharge zone also includes more distant regions, covered by typical urban housing. New multi-family buildings, parking lots and streets were constructed, modifying direct recharge from the surface, which is typical of urban agglomerations. The construction of multi-story buildings requiring deep foundations permanently modified the structure of groundwater flow [59,60]. Some buildings required drainage during construction, which periodically changed the hydrodynamic system. This factor is difficult to estimate because of construction activities in the urban area, but has a significant influence on the local position of the water table.
In urban areas with high water demands, the hydrodynamic conditions are modified due to significant, focused groundwater sources. There are only a few groundwater sources near the research station and in all of Warsaw; therefore, restoration of piezometric surfaces was documented in the Oligocene aquifer. Leakage of the water supply system and infiltration of water from the leaking network to the aquifer, documented even in the vicinity of the research station (Figure 9), can be observed in Warsaw.
Water 2020, 12, x FOR PEER REVIEW 12 of 17

Discussion
Recharge in urban areas depends on a number of factors, including climate conditions and the dynamics of their changes [52,53], lithology of subsurface sediments [54,55], topography [56], land use and land cover [52,57], and thickness of the aeration zone [55,58].
The change of climate conditions in Warsaw was recorded as anomalously high precipitation in 2017 (821 mm) and anomalously low precipitation in 2015 (297 mm) (Figure 4). Measurements of groundwater levels in the research station have indicated that this change has no direct impact on the change of groundwater levels, even in the case of the shallowest Quaternary aquifer. In all analyzed aquifers, the influence of precipitation on water table fluctuations was insignificant due to the rather unfavourable filtration parameters in the aeration zone and the impact of urban infrastructure, especially in the subsurface zone. The reaction of the water table, even after intense precipitation, was greatly delayed. Multiannual trends of changes in the position of the water table in Quaternary aquifers and a stable trend of increasing potentiometric surface in the Oligocene aquifer are notable.
In urban areas, the magnitude of recharge depends on land use. The research station is located in Warsaw's centre and was developed at a time when there was hardly any infrastructure around it. However, in the 27-year interval, numerous changes in its surroundings have taken place, particularly related to the development of the University of Warsaw, Ochota Campus infrastructure, which may have had a significant impact on groundwater recharge and flow in the shallowest zone. The recharge zone also includes more distant regions, covered by typical urban housing. New multi-family buildings, parking lots and streets were constructed, modifying direct recharge from the surface, which is typical of urban agglomerations. The construction of multi-story buildings requiring deep foundations permanently modified the structure of groundwater flow [59,60]. Some buildings required drainage during construction, which periodically changed the hydrodynamic system. This factor is difficult to estimate because of construction activities in the urban area, but has a significant influence on the local position of the water table.
In urban areas with high water demands, the hydrodynamic conditions are modified due to significant, focused groundwater sources. There are only a few groundwater sources near the research station and in all of Warsaw; therefore, restoration of piezometric surfaces was documented in the Oligocene aquifer. Leakage of the water supply system and infiltration of water from the leaking network to the aquifer, documented even in the vicinity of the research station (Figure 9), can be observed in Warsaw. The content of water leaking from the water supply system due to leakage or failure is roughly estimated in Warsaw at about 18% (mean value of water loss for Poland according to the Central Statistical Office for 2006), which gives the absolute value of 61,200 m 3 water per day. This value should be increased by another 18% resulting from leakage from the sewage system, which is 4160.7 km long in the Warsaw agglomeration [24]. Taking into account that these are estimates in many cases are understated, it can be concluded that a total of about 125,000 m 3 of water per day recharges The content of water leaking from the water supply system due to leakage or failure is roughly estimated in Warsaw at about 18% (mean value of water loss for Poland according to the Central Statistical Office for 2006), which gives the absolute value of 61,200 m 3 water per day. This value should be increased by another 18% resulting from leakage from the sewage system, which is 4160.7 km long in the Warsaw agglomeration [24]. Taking into account that these are estimates in many cases are understated, it can be concluded that a total of about 125,000 m 3 of water per day recharges the shallow Quaternary aquifer within this urban area. In Barcelona, the water supply network losses provide 22% of the total recharge of aquifers [48].
Water use does not increase proportionally with the increasing population of Warsaw. This may result from more economical water use due to rising costs and ecological lifestyles. This lower water source results in smaller amounts of water potentially recharging the subsurface zone due to leakage of the water system. The cumulative impact of many elements causing periodical rise or fall of recharge in urban areas allows for an assessment of their impact on the groundwater environment. In the case of the three analyzed aquifers, water level fluctuations are slow and occur over longer periods. Quaternary aquifers are characterized by multiannual intervals with higher recharge and much longer episodes of water table fall, which, compared to the entire 27-year interval, results in a falling trend towards decreasing water table levels. Such a case is very distinct in the shallow Quaternary aquifer since 2012. The Oligocene aquifer displays a trend of gradual but uneven restoration of hydrostatic pressure over the analyzed time span. The highest increase in the potentiometric surface was found in 1997-2003, after which the increase was much slower. The restoration of pressure in the Oligocene aquifer in the Warsaw area is confirmed by monitoring observations performed by the Polish Hydrogeological Survey [61].

Conclusions
Analysis of the changes in groundwater levels in several aquifers-two Quaternary ones and an Oligocene one-points to diverse factors controlling their location. In Quaternary aquifers, recharge conditions and dynamics of changes depend on local factors: infiltration conditions in the region of the monitoring system and development of the area, including elements characteristic of urbanized regions. In the case of the deep Oligocene aquifer, the position of the potentiometric surface depends on water circulation in the regional hydrogeological system, which covers a much larger area than the Warsaw agglomeration. Statistical analysis is a useful tool for interpreting multiannual measurements, although numerous natural and anthropogenic aspects need to be considered for its correct application. Groundwater monitoring systems usually exist in urban areas for several years, whereas statistical analysis requires long-term measurement series to analyse the trend of changes. The position of the water tables of the Quaternary aquifers across a long-term interval showed their fall, whereas the deep Oligocene aquifer displayed an opposite trend, with the restoration of hydrostatic pressure resulting in a rise in the potentiometric pressure by about 10 m.
Analysis of the changes in groundwater levels in time and factors influencing them in Warsaw does not show the impact of climate change (e.g., through extreme precipitation) on groundwater. Seasonal trends and long-term changes in levels have not been observed in the studied aquifers in the urban area.
The studies have shown that groundwater recharge and drainage in urban areas are components of natural processes related to hydrogeological and climate conditions and anthropogenic factors resulting from the groundwater source, water distribution, and changes in spatial management. When analysing the cumulative impact of these factors, often with opposite effects with regard to groundwater recharge, a relatively rapid decrease in the water table position may be observed in the shallowest groundwater aquifer. Therefore, it is critical to correctly manage water resources and increase water retention conditions through increased area covered by green zones, increased application of precipitation for urban green management and application of innovative solutions to enhance the infiltration of precipitation in high-density urban zones. Analysis of data sets containing long-term measurement series allows for a reliable assessment of the existing conditions (water efficiency), aids in planning sustainable water management and, to some degree, allows for the prediction of the effects of irrational water management or inappropriate development of urban infrastructure.