Urbanization—Its Hidden Impact on Water Losses: Prądnik River Basin, Lesser Poland

Urban development causes multiple water losses. Some of them may be ignored but some could have a huge influence on the whole catchment, including soil drought. As urban sprawl rises, space for unaffected infiltration and retention is increasingly limited. The objective of this study was to backcast and to estimate water-retention loss due to urbanization during the period of 1990–2018. We used landcover data, meteorological and hydrological data and data on soil water-holding capacity. Water-retention loss was expressed as soil water retention capacity loss, net precipitation loss and total sum of precipitation loss. Historical change in urban extension has led to large impacts on the hydrological cycle of the study area. Progressive urban development caused water-retention losses which range from 3.380 to 14.182 millions of cubic meters—depending on the methodology used. Hydrological analysis showed the lack of a significant trend (decrease trend) of low flow which is caused by the high percentage of natural land use in the upper part of catchment. Our results show that backcasting of water retention change using CLC data (a) brings new and plausible data on retention loss, (b) is possible to replicate and (c) data used are common and easy-to-get.


Introduction
Water resources are characterized by spatial and temporal variability. Changing availability, quantity and quality of water, leads to problems which mankind has faced for centuries. Floods and droughts occur irregularly, but with major impact on countries across Europe [1][2][3]. River deterioration, the legacy of decades of scarce water management, demonstrates itself in the decreased water quality detected in most streams [1]. The mentioned problems are largely the consequence of climate and land use changes (amplified by the absence or non-observation with, environmental legislation; [4]). Consequently, numerous struggles over water resources are increasing.
To reverse or mitigate this trend is the object of many adaptation strategies. In Poland, policy arrangements have been implemented to reduce discharge from both unaffected and artificial, man-made running water since the 1960's. These actions were intended to store runoff in pools and terrain hollows and to support the retention of rainfall in soil environment and aquifers. These arrangements are called "small retention" contrary to water storage in large reservoirs and barrages [5,6]. Technical solutions of traditional water retention make such arrangements cost-ineffective [7], and moreover, cause considerable environmental [8] and health impacts [9].
Land use/landcover is the crucial landscape characteristic transformed by man in the most significant way and therefore modifies runoff and the resulting landscapes' hydric potential (LHP) [10][11][12]. Land use can subvert the hydrologic balance and the dividing of rainfall into evapotranspiration, runoff and groundwater flow [13]. Outflow and stream discharge increase when natural ecosystems (particularly forests) are cleared [14,15]. Landcover alteration can considerably alter maximum flows in river basins [16], especially within the areas where most of the outflow is formed on sites with fast runoff generation [17]. Landcover change is often related to several categories of drought: meteorological, hydrological, agricultural and socioeconomic droughts [18]. Once subsurface water systems are affected by droughts, foremost groundwater recharge, subsequently groundwater levels and groundwater discharge reduce. Groundwater droughts usually observe on a timescale of weeks to years [19]. Droughts impact both surface and underground water resources. Moreover, decreased water supply, worsened water condition, harvest failure, reduced electricity generation, disrupted aquatic ecosystems, as well as postponed recreation or leisure activities are linked with drought [20].
One of the prevalent risks related to land use change and LHP is rapid urban extension [21]. Urbanization comes with a variety of environmental problems for both the local and downstream ecosystems [22]. These environmental challenges are a direct consequence of the physical, chemical and biologic changes to hydrological systems [23,24]. The reduction of permeable surfaces reduces the infiltration into soil [25]; the application of artificial drainage alters or substitutes natural pathways. Impervious surfaces elevate the dangers of augmented discharge and lower infiltration rate. The hydrological response results in the reduction of baseflow and groundwater recharge and more frequently induces the low flow of urban streams. Replacing vegetated areas with impervious surfaces causes more rapid precipitation flows into adjacent streams [26]. As a consequence, the volume and frequency of floods, as well as the peak discharge, increase in nearby streams. Impervious surface cover is positively correlated with flashiness, particularly for small catchment gauges [27]. As the impermeable surface cover rises to 10-20%, runoff increases two-fold; 35-50% impervious surface cover increases runoff three-fold compared to forested catchment [23]. Transformation of natural, semi-natural and suburban landscape to impermeable surfaces, can cause a variety of unfavorable aftermaths [28]. Land use/landcover is the important factor that affects the infiltration and holding of a considerable volume of rainfall. At the same time, there is a high threat of discharge resulting in a faster response [29] and a greater magnitude of stream flow [30], which are associated with deluges from intra-urban waters [21].
The goal of this study is to detect whether a high extent of urban growth upon formerly rural landscape has resulted to serious decreases of water retention. This will be completed through four sequential steps. Foremost, records from hydro-meteorological observing are utilized to characterize the amount of precipitation and evapotranspiration. Next, the retention capabilities of soil systems will be analyzed. Further, the temporal and spatial change in the urban expansion is surveyed for each time lapse. Finally, we will explore how the urbanization and related impervious surfaces have affected the water retention over the surveyed period of urban expansion. The observations and modeling are focusing on the Prądnik river basin, Lesser Poland. The outcomes of this study are discussed to provide information on the effects of urban growth on the water retention (hydric) potential of the river basin.

Study Area
To show the temporal and spatial change of water retention we chose the Prądnik river basin, close to Kraków, Poland. The city of Kraków and its vicinity have long-lasting problems with water quality, repeated deluges or droughts [31]. The Prądnik river discharge has changed considerably. In the 1960s the mean annual discharge was 0.490 m 3 ·s −1 and towards the beginning of the 21st century it fell down to 0.330 m 3 ·s −1 . Such a change results from both natural and human influences. Situated in the headspring part of the Prądnik river basin, significant yield water sources have been turned into water intakes planned for economic use that directly reduced the discharge by more than 0.1 m 3 ·s −1 [32]. The Prądnik river, is important for relevant for Kraków. It has a direct influence on the water related issues of the built-up area such as deluges, sediment loads and pollution. Significant flood threat is distinctive particularly in the urbanized part of the study area. This part of the study area has a dense storm sewage system and is typical of the impervious cover [33].
The Prądnik river basin ( Figure 1) includes rural, suburban, heavily urbanized and natural landscape. The study area is drained primarily by the Prądnik river. It is a left-bank tributary of the Vistula River. The spring area of the river is found near the Sułoszowa village at an altitude of 450 m ASL; the river is 33.4 km long. The Prądnik river basin is situated in the south-eastern part of the Kraków-Częstochowska Upland and covers an area of 193.1 km 2 . The basin is typical of agricultural landscape with 55% of its area. The rest is covered by urbanized areas (20.1%) and forests (13.7%) and other types of land use/landcover (LU/LC).
Water 2020, 12, x FOR PEER REVIEW  3 of 18 into water intakes planned for economic use that directly reduced the discharge by more than 0.1 m 3 •s −1 [32]. The Prądnik river, is important for relevant for Kraków. It has a direct influence on the water related issues of the built-up area such as deluges, sediment loads and pollution. Significant flood threat is distinctive particularly in the urbanized part of the study area. This part of the study area has a dense storm sewage system and is typical of the impervious cover [33]. The Prądnik river basin ( Figure 1) includes rural, suburban, heavily urbanized and natural landscape. The study area is drained primarily by the Prądnik river. It is a left-bank tributary of the Vistula River. The spring area of the river is found near the Sułoszowa village at an altitude of 450 m ASL; the river is 33.4 km long. The Prądnik river basin is situated in the south-eastern part of the Kraków-Częstochowska Upland and covers an area of 193.1 km 2 . The basin is typical of agricultural landscape with 55% of its area. The rest is covered by urbanized areas (20.1%) and forests (13.7%) and other types of land use/landcover (LU/LC). The Prądnik river basin is formed typically by loess, limestone and river deposits with favorable water transmissivity (0.001-0.01 m 2 •s −1 ). On these substrates are formed cambisols (61.0%), podzols (24.8%) and fluvisols (7.3%). Sandy loams (90.1%) dominate together with silt (4.7%) and loamy sands (4.7%). The Prądnik river basin is situated at an altitude of 220-481 m ASL. with an average annual air temperature of +7.5 °C. The average annual precipitation is 580 mm in the eastern parts and 760 mm in the north-west. The average annual potential evapotranspiration ranges from 300-320 mm in the valleys to more than 400 mm in the plateau region of Ojców. Further data on landscape characteristics are shown below ( Figure 2). The Prądnik river basin is formed typically by loess, limestone and river deposits with favorable water transmissivity (0.001-0.01 m 2 ·s −1 ). On these substrates are formed cambisols (61.0%), podzols (24.8%) and fluvisols (7.3%). Sandy loams (90.1%) dominate together with silt (4.7%) and loamy sands (4.7%). The Prądnik river basin is situated at an altitude of 220-481 m ASL. with an average annual air temperature of +7.5 • C. The average annual precipitation is 580 mm in the eastern parts and 760 mm in the north-west. The average annual potential evapotranspiration ranges from 300-320 mm in the valleys to more than 400 mm in the plateau region of Ojców. Further data on landscape characteristics are shown below ( Figure 2).

Data and Methods
The historical change in urban extent and impervious cover was analyzed from CORINE landcover (CLC) inventory in the years 1990, 2000, 2006, 2012 and 2018 [34]. Artificial surfaces were considered as impermeable. Impervious cover included the following CLC classes: 111-continuous urban fabric; 112-discontinuous urban fabric; 121-industrial or commercial units; 122 roads and rail networks and associated land; 123-port areas; 124-airports; 131-mineral extraction sites; 133-construction areas and 142-sport and leisure facilities. The area of artificial man-made surfaces within listed CLC classes cannot be considered as 100% impervious. To estimate exact data on impervious cover we determined the correction factor (cf ) based on the weighted average of imperviousness degree for each CLC class.
The imperviousness degree was obtained from IMD status layer with 100 m pixel size for the reference year 2006 [35] and the correction factor was determined as follows: CLC 112-cf 0.56, CLC 121-cf 0.74, CLC 133-cf 0.52 and CLC 142-cf 0.46.

Data and Methods
The historical change in urban extent and impervious cover was analyzed from CORINE landcover (CLC) inventory in the years 1990, 2000, 2006, 2012 and 2018 [34]. Artificial surfaces were considered as impermeable. Impervious cover included the following CLC classes: 111-continuous urban fabric; 112-discontinuous urban fabric; 121-industrial or commercial units; 122 roads and rail networks and associated land; 123-port areas; 124-airports; 131-mineral extraction sites; 133-construction areas and 142-sport and leisure facilities. The area of artificial man-made surfaces within listed CLC classes cannot be considered as 100% impervious. To estimate exact data on impervious cover we determined the correction factor (cf) based on the weighted average of imperviousness degree for each CLC class. The imperviousness degree was obtained from IMD status layer with 100 m pixel size for the reference year 2006 [35] and the correction factor was determined as follows: CLC 112-cf 0.56, CLC 121-cf 0.74, CLC 133-cf 0.52 and CLC 142-cf 0.46.
The soil-retention potential of the study area was estimated according to the total available (holding) capacity. Water-holding capacity is primarily controlled by soil texture. A review of waterholding capacity of the main soil texture groups is shown in Table 1. Data are averages and vary with structure and organic matter differences. We used data on water-holding capacity according to Kutílek and Nielsen [36].  The soil-retention potential of the study area was estimated according to the total available (holding) capacity. Water-holding capacity is primarily controlled by soil texture. A review of water-holding capacity of the main soil texture groups is shown in Table 1. Data are averages and vary with structure and organic matter differences. We used data on water-holding capacity according to Kutílek and Nielsen [36].
Water-retention loss due to urbanization was described in three ways; (i) soil water-holding capacity (SWHC) loss, (ii) volume of average annual climatic water balance-an expression of net precipitation which could be retained under natural conditions (potential water loss-PWL) and (iii) volume of mean annual precipitation falling on the total impervious cover of the study area (total water loss-TWL).
Slope inclination was derived from a digital elevation model (DEM). The DEM was created after analyzing photogrammetric stereoscopic air images at a scale of 1:26,000. It was prepared on the basis of a height points network with a 20 m resolution. A soil vector map was used to determine soil type and texture. The DEM and the soil-agronomic vector map were obtained from the Provincial Center for Geodesic and Cartographic Documentation. CORINE Landcover 1990, 2000 data in 100 m resolution [34], were used to acquire data about the landscape structure of the study area. Data on imperviousness degree in 100 m resolution were obtained from the IMD status layer for the reference year 2006 [35]. The spatial analysis, data processing, synthesis and evaluation of individual sub-basins were realized by ArcGIS 10.3 software. Data from 32 meteorological stations from the period of 1971-2000 were used to estimate the spatial distribution of mean annual precipitation, mean annual potential evapotranspiration and climatic water balance. The hydrological data from the period of 2002-2017 and meteorological data from the period of 1990-2017 from the station in Ojców were used for drought phenomenon detection. The hydrometeorological data (daily precipitation, daily temperature and daily discharges) were obtained from the Institute of Meteorology and Water Management of the National Research Institute in Warsaw.
Climatological analyses of precipitation time series were performed. The relative precipitation index (RPI) was utilized as one of the widely used descriptors of meteorological drought. The RPI index is calculated by dividing actual precipitation by normal precipitation for the time being considered and multiplied by 100 according to the equation: where: Analyses using the RPI index are very effective when used for a single region or a single season. It can be calculated for a variety of time scales, including monthly, seasonal or annual. The pluvial characteristics were based on the Kaczorowska [41] classification the RPI index, which propose two comparative criteria: one related to monthly sums of precipitation and the other related to the sum of precipitation in longer periods (Table 2).
Additional, hydrological analysis trends of low flow time series were performed. The low flow was used as one of the descriptors of hydrological drought. This hydrological characteristic is significant for human impact as LU/LC change or water take. The hydrological analysis covers the determination of descriptive statistics and verification of the significance of the trend for low-flows time series for the Prądnik basin. The trend of low flow data were analyzed for Ojców gauge station (catchment area equal 67.5 km 2 ) for the 2002-2017 period. The low flows characteristic was assumed due to the fact that it is one of the drought hydrological indicators. Descriptive statistics included the following measures: location (minimum, maximum, average, median), dispersion (standard deviation, coefficient of variation) and measures of distribution shapes (skewness coefficient, kurtosis).
Additionally, the trend signification was verified; three tests were used: Mann-Kendall, Pettitt and Kruskal-Wallis. The calculations were carried out for the significance level α = 0.05. The significance of the trend of the low flow time-series, was carried out using the Mann-Kendall (MK) test. The null hypothesis H 0 of the test assumed no monotonic data trend, while an alternative stated that such a trend exists. The calculations were carried out for the significance level α = 0.05. The S Mann-Kendall statistics were determined based on the equation [42,43]: where: n-number of elements of the time series.
Normalized statistics Z were calculated according to the equation: where: Var(S)-variance S, determined from the equation: The main premise of the MK test used is the lack of autocorrelation in a series of data. In the case of the analysis time series such relations could arise. When autocorrelation is found, this phenomenon leads to an underestimation of the Var(S) variance. Therefore, an adjustment for correction of variance has been included, calculated only for data with significant partial autocorrelation [44]: where: n n * s -effective number of observations calculated as: where: k-shift; ρ k -value of the next significant autocorrelation coefficient.
The Pettitt test allows the detection of a single shift at an unknown time t. The null hypothesis assumed no change in the distribution of a sequence of random variables. The test is based on formulas [45]: where: The test statistics U t,T depends on D ij and is calculated as: The test statistics is assessed for all random variables from 1 to T. The most significant change point is selected where the value of |U t,T | is the largest: A change point occurs at time t when the statistic K T is significantly different from zero at a given level.
The Kruskal-Wallis test allows the detection of statistical homogeneity. The inference involved assigning rank to the ordered elements of all samples and determination of rank sum for each sample. The test statistics H is calculated as [46]: where: N-total number, R i -the total sum of ranks in the i-th group, n i -the number in the i-th group.
When differences between the analyzed sums are small, the null hypothesis (H0) assuming that all samples originate from the same general population (are homogeneous) is true. A critical region was defined by Pearson's statistic χ 2 with n − 1 degrees of freedom [47].

Results and Discussion
The presented study provides a novel perspective on the nature of the urban development-water relationship. This section is organized to provide the outcomes from: (i) pluvial conditions and trends of low flow, (ii) determining the water retention potential of the Prądnik river basin, (iii) mapping any historical change of urban extent and associated impervious cover within the Prądnik river basin and (iv) backcasting of water-retention loss using historical levels of impervious areas. Results are discussed in relation to urban extent and how urban growth has changed the retention potential of the study area.

Pluvial Conditions
The annual precipitation in the period 1990-2017 for the Prądnik basin was 743 mm and was higher than the normal precipitation for Poland. According to Czarnecka and Lencewicz [48] the annual normal precipitation in Poland for the period 1951-2000 was 594 mm. Higher precipitation in the Prądnik basin results from the location of the catchment area on the higher part of the Krakowsko-Częstochowska Upland. The longtime course of the precipitation is characterized by a significant fluctuation of the values. However, the statistical trend has not been significant (Figure 3).  The instability of precipitation in the Prądnik basin is evidenced by considerable monthly and annual precipitation. The calculated values of the RPI index (Table 3) demonstrate the high variability of the precipitation. The indicator shows the atmospheric precipitation deficiencies and excesses and may be useful for risk assessment of drought or higher river flow.  The instability of precipitation in the Prądnik basin is evidenced by considerable monthly and annual precipitation. The calculated values of the RPI index (Table 3) demonstrate the high variability of the precipitation. The indicator shows the atmospheric precipitation deficiencies and excesses and may be useful for risk assessment of drought or higher river flow.
Extremely dry years did not occur in the analysis period (1990-2017); only one very dry year was observed (2011) in which the precipitation amounted to only 64% of the normal precipitation. There were also eight dry years (1990-1993, 2003, 2004, 2006 and 2015). In the whole analysis period, extremely dry months occurred 14 times, very dry 49 times and dry 77 times and were most common in autumn and winter. The extremely dry months, very dry and dry represented 42% of the analysis period of 1990-2017. Similar results were observed in other parts of Poland [49,50].

Trends of Low Flow
As part of the initial data analysis, basic descriptive statistics were determined for the low flow time-series. Figure 4 shows the course of low flows in the analyzed period: 2002-2017. The results of the calculations are summarized in Table 4.
Based on the results summarized in Table 4 it was found that the difference between minimum and maximum observed low flows values was 29%. The dynamics of changes in the low flows remained at an average level. This is evidenced by the coefficient of variation whose value was almost 30%. This variation coefficient can be explained by the stationary conditions of factors affecting the course of low flows in the analyzed period (meteorological factors, human activities). Analyzing the values of the skewness coefficient, it was shown that it is greater than 0. Therefore, right-sided asymmetry of the empirical distribution of random variables was found. This is due to the fact that in the analyzed time-series most observations are smaller than their average value. Therefore, the average value is greater than the median. In turn, the kurtosis values for time-series indicated the leptokurtic (positive value) empirical distribution of the studied random variables. The analyses of the significance of low flows trend were conducted using the Mann-Kendall, Pettitt and Kruskal-Wallis tests. The results of this analysis are summarized in Table 5. Extremely dry years did not occur in the analysis period (1990-2017); only one very dry year was observed (2011) in which the precipitation amounted to only 64% of the normal precipitation. There were also eight dry years (1990-1993, 2003, 2004, 2006 and 2015). In the whole analysis period, extremely dry months occurred 14 times, very dry 49 times and dry 77 times and were most common in autumn and winter. The extremely dry months, very dry and dry represented 42% of the analysis period of 1990-2017. Similar results were observed in other parts of Poland [49,50].

Trends of Low Flow
As part of the initial data analysis, basic descriptive statistics were determined for the low flow time-series. Figure 4 shows the course of low flows in the analyzed period: 2002-2017. The results of the calculations are summarized in Table 4.     Based on the results summarized in Table 5 it was found that, for the analyzed period, there has not been a significant statistical trend of low flows in the Prądnik basin, but a decreasing of low flow is visible over the 2002-2017 period. This is confirmed by the values of the probability p for obtained test statistics. In each case the probability p was lower than significant level (α = 0.05). The lack of statistically significant trends means that the investigated random variables derive from the same general population. Thus, in the analyzed multiyear period, no factor appeared to be significantly affecting the processes that shape low flows. The similar results of trend analysis for low flows, in other river basins of the upper Vistula basin were obtained by Młyński et al. [51]. The results obtained from the research can be related to the outcomes regarding the trend analysis of the factors significantly affecting the flow of waters in the river basins of the Upper Vistula water region. Results brought by Niedźwiedź et al. [52] presented no significant trends concerning the indicators describing daily precipitation. Analyses carried out by Falarz [53] confirmed the lack of significant trends for the duration of snow cover remaining in mountainous catchments. The supply of rivers in the Upper Vistula water region comes primarily from snowmelt runoff and from precipitation-hence it is assumed that the rhythm of effluents is repeated by the low flows. The lack of significant trend of low flow in the Prądnik river basin is caused by the high percentage of natural land use in the upper part of the catchment and the positive value of climatic water balance. In result, the upper part of the Prądnik river basin has a higher water storage potential that moderately affects the impact of climatic conditions on flow regime. Higher water storage potential in the upper part of the study area is linked with the higher value of LHP compared to rest of the river basin [11].

Water Retention Potential of the Prądnik River Basin
The capability of soil environment to control water resources is an essential ecosystem service. Infiltrated precipitation through soil is filtered, accumulated for vegetal utilization and distributed through flow pathways to groundwater and surface water bodies. Hence, the sustainability of water resources (both amount and quality) is inherently affected by soil [54]. Based on the variability of spatial soil texture data ( Figure 2) and SWHC data [36], we estimated a mean soil water-holding capacity of the Prądnik river basin 29.221 mil·m 3 (maximum estimated soil water-holding capacity is 67.831 mil·m 3 ).

Historical Change in Urban Extent and Impervious Cover
Mapping the change in urban extent across the Prądnik river basin from CORINE Landcover inventory over the period from the 1990's to 2018 confirmed that the study area has undergone a significant increase of the impervious areas ( Figure 5). The historical change in urban development and related impervious cover illustrates the quantitative influence of the urbanization across the study river basin over the period 1990-2018. Within the Prądnik river basin, the impervious cover rose considerably from 6.7% to 11.5%, with the greatest development seen in 2006 and 2012. For further data on the change of impervious cover see Table 6; the presented area of impervious surfaces is after correction factor (cf ) was used. With all estimates of historical urban extension and related impervious cover, there is some associated uncertainty. As seen from Table 6 there is uncertainty of CLC 2000 data. Compared to CLC 1990, there is a deficiency of 234 ha within the CLC class 112-discontinuous urban fabric. The total reliability of CLC 2000 is 87.0% ± 0.8%; within the class 112, there is the highest class reliability >95.0% [55].

Backcasting the Water Retention Change of Progressive Urban Development
The mapping of urban development provided the basis for estimating historical values of impervious cover. The backcasting of the water-retention loss using a mean annual precipitation, climatic water balance (CWB), soil water-holding capacity (SWHC) and the extent of impervious surfaces (IS) change enables a comparison for every year (1990, 2000, 2006, 2012 and 2018). Results are summarized in Table 7. An increase in total impervious cover within the Prądnik river basin from 1297 ha in the 1990´s to 2230 ha in 2018, has led to a substantial/consequential water loss from 8.122

Backcasting the Water Retention Change of Progressive Urban Development
The mapping of urban development provided the basis for estimating historical values of impervious cover. The backcasting of the water-retention loss using a mean annual precipitation, climatic water balance (CWB), soil water-holding capacity (SWHC) and the extent of impervious surfaces (IS) change enables a comparison for every year (1990, 2000, 2006, 2012 and 2018). Results are summarized in Table 7. An increase in total impervious cover within the Prądnik river basin from 1297 ha in the 1990´s to 2230 ha in 2018, has led to a substantial/consequential water loss from 8.122 mil·m 3 to 14.182 mil·m 3 . An increase of almost 75% was reached only in 28 years. A comparison of total annual water retention potential of the Prądnik river basin (29.221 mil·m 3 resp. 67.831 mil·m 3 ) and total water-retention loss due to urbanization (14.182 mil·m 3 ) describes how impervious areas affect the hydrologic cycle. It is noticeable that impervious areas covering only 11.5% (2230 ha) of the study area are responsible for water loss which equals 48% of mean river basin retention capacity (or 21% of maximum retention capacity). Table 7. Change of water retention potential during 1990-2018, the Prądnik river basin.
Year IS (ha) IS-impervious surfaces; PWL-potential water loss; SWHCL-soil water-holding capacity loss; TWL-total water loss.
Progressive urban development is responsible for soil sealing which averts precipitation to infiltration. SWHCL expresses the loss of soil´s space to retain water. Every single square meter of impervious cover decreases soil water-holding capacity whose volume ranges from 100 to 400 mm·m −3 [36]- Table 1. Soil sealing led to a decrease of mean soil water-holding capacity by 1.375 mil·m 3 (or by 3.360 mil·m 3 of maximum SWHC) during the period of 1990-2018 ( Table 7). Loss of SWHC corresponds to LHP loss and, hence, to the loss of one of the key ecosystem services.
The influence of climate must be involved in all estimations of water availability in the Prądnik river basin. Anthropogenic activity aside, surface water bodies and groundwater quantity is determined by the balance between atmospheric inputs and losses related to evapotranspiration. The relationship between precipitation and evapotranspiration is a well-established indicator of water availability [56]. Potential water loss describes the depth of net precipitation which cannot be retained in the soil environment due to urbanization. PWL states the depth of CWB in the natural environment, where infiltration and retention of water is unblocked. The annual climate water balance varies from 196 mm in the east and north of the study area, to 496 mm in the valleys of the Prądnik river ( Figure 6). As the area of impervious surfaces increased the volume of water loss rose as well. In the 1990's the area of IS was 1297 ha, the PWL equaled to 3.095 mil·m 3 ; nowadays, with a sealed area of 2230 ha, the PWL raised to 5.232 mil·m 3 .
However, the total water loss is far greater. In urbanized areas, precipitation falling on artificial surfaces either evaporates or runs off. The runoff water drains quickly, and, in the long term, less surface water remains available for evapotranspiration [57]. Urban evapotranspiration can be neglected because it is considerably less than from neighboring rural areas, because of the contrasts in the hydrological properties of building materials and vegetation-covered soils [58]. Hence, it would be better to use the mean annual precipitation as the description of total water loss. It is obvious that TWL exceeds the water retention potential of the soil environment. TWL is almost 3 times greater than PWL and clearly illustrates the actual water loss of the urban environment. It is evident ( Table 7) that the historical impacts from urbanization are massive and complex. Sealed surfaces covering 2230 ha (reference year 2018) are concerned with the runoff and/or evaporation of 14.182 mil·m 3 per year. Such progressive urban development and consequent TWL is often linked with the prevention of groundwater recharge [59][60][61] and decreased low flows between rainfall events [62]. As urban areas develop, the resulting LU/LC changes affect hydrologic processes through alterations in recharge, runoff and evapotranspiration which, in turn, affect the quantity of water bodies [63]. However, the total water loss is far greater. In urbanized areas, precipitation falling on artificial surfaces either evaporates or runs off. The runoff water drains quickly, and, in the long term, less surface water remains available for evapotranspiration [57]. Urban evapotranspiration can be neglected because it is considerably less than from neighboring rural areas, because of the contrasts in the hydrological properties of building materials and vegetation-covered soils [58]. Hence, it would be better to use the mean annual precipitation as the description of total water loss. It is obvious that TWL exceeds the water retention potential of the soil environment. TWL is almost 3 times greater than PWL and clearly illustrates the actual water loss of the urban environment. It is evident ( Table 7) that the historical impacts from urbanization are massive and complex. Sealed surfaces covering 2230 ha (reference year 2018) are concerned with the runoff and/or evaporation of 14.182 mil. m 3 per year. Such progressive urban development and consequent TWL is often linked with the prevention of groundwater recharge [59][60][61] and decreased low flows between rainfall events [62]. As urban areas develop, the resulting LU/LC changes affect hydrologic processes through alterations in recharge, runoff and evapotranspiration which, in turn, affect the quantity of water bodies [63].
The urban area extent almost doubled during the observed time period of 1990-2018. Nowadays, covering 2230 ha (11.5%) of the total, impervious surfaces may have a significant effect on the hydrology of the Prądnik basin. Booth and Jackson [64] state that changes in the hydrologic responses of river basins begin at very low levels of urban development; with measured effects when the area of impervious cover reaches about 10%. In addition, Oudin et al. [65] stated that high, low and mean flows were impacted at a threshold of a 10% impermeable extent. Yang et al. [66] indicated that 3%-5% impervious cover is the threshold value-beyond which urban development effects start to have a statistically significant impact on streamflow regime. Model simulations at 6% of the total impervious area confirmed increases in maximum flows in urbanized river basins when compared to rural river basins in comparable physiographic locations. For shorter recurrence periods (1-2 years), the increase in maximum flow due to urban development is more evident than for longer recurrence interval peak flows (10 years or greater) [67]. Walega et al. [68] showed that the increase of the impervious area in Suduł Dominikanski (right tributary of Prądnik river) increased peak flow and flood volume and also decreased the response of catchment of heavy storm. Bhaduri et al. [69] found a linear relationship between the extent of the impervious cover and average annual runoff, later confirmed by Assouline and Mualem [70]. A research on the city of Leipzig, Germany, demonstrated that runoff had more than doubled in period of 1940-2003 due to the progressive urban The urban area extent almost doubled during the observed time period of 1990-2018. Nowadays, covering 2230 ha (11.5%) of the total, impervious surfaces may have a significant effect on the hydrology of the Prądnik basin. Booth and Jackson [64] state that changes in the hydrologic responses of river basins begin at very low levels of urban development; with measured effects when the area of impervious cover reaches about 10%. In addition, Oudin et al. [65] stated that high, low and mean flows were impacted at a threshold of a 10% impermeable extent. Yang et al. [66] indicated that 3-5% impervious cover is the threshold value-beyond which urban development effects start to have a statistically significant impact on streamflow regime. Model simulations at 6% of the total impervious area confirmed increases in maximum flows in urbanized river basins when compared to rural river basins in comparable physiographic locations. For shorter recurrence periods (1-2 years), the increase in maximum flow due to urban development is more evident than for longer recurrence interval peak flows (10 years or greater) [67]. Walega et al. [68] showed that the increase of the impervious area in Suduł Dominikanski (right tributary of Prądnik river) increased peak flow and flood volume and also decreased the response of catchment of heavy storm. Bhaduri et al. [69] found a linear relationship between the extent of the impervious cover and average annual runoff, later confirmed by Assouline and Mualem [70]. A research on the city of Leipzig, Germany, demonstrated that runoff had more than doubled in period of 1940-2003 due to the progressive urban development [71]. Similarly, Perry and Nawaz [72] demonstrated comparable results in the urbanized area of the city of Leeds, Great Britain. A 75% expansion of the urban area throughout the Prądnik river basin has led to a 75% increase in total water loss. This finding is consistent with studies that reported urbanization promoted quick runoff [73][74][75][76].
The LU/LC change has a direct impact on the high quantity of water loss-not only theoretical, but with visible impacts on ecosystems, groundwater table, crop yields and climate as well. It is not only due to urbanization, but also agriculture and forestry affect hydrologic balance significantly. Many draining facilities have been established across the landscape. In Poland, about 2 million hectares of grasslands are drained with a system of open ditches. Most of them are located in floodplains [76]. Networks of open ditches contribute to the decrease of the groundwater table level. In the study area, there are 53.3 km of ditches. When the river basin is created by pervious geological structures, this network of open ditches decreases the groundwater level in the adjacent areas as well as in the case of the Prądnik river basin. Quick runoff from the drainage facilities is observed typically during the spring. Consequently, for a period of the growing season, there is not sufficient water for the vegetation.
Soil environment under natural forests is usually porous, typical with favorable infiltration and subsequently low rates of runoff. Such conditions are not automatically the case for managed forest stands, especially where no natural understory of vegetation is maintained. Forest management activities concerning site preparation, drainage, road network development and logging may have detrimental effects on water resources [77,78] and the transport of sediments into streams [79].

Conclusions
The study presents the backcasting of the water retention change due to progressive urban development. Urban development in the Prądnik river basin largely affects the hydrologic balance. Impervious surfaces limit infiltration and space for retention of water while promoting runoff and evaporation. We have provided three ways of insight into urban related water-retention loss: • water-retention loss could be perceived as damage to the soil´s ability to retain precipitation due sealing by impervious surfaces; • water-retention loss due to urbanization could be expressed by volume of net precipitation. It provides the comparison of water volume which could be retained in the natural environment, where infiltration and retention is unblocked; • water-retention loss can be stated as the total sum of runoff and evaporation from impervious areas, as well as retention of sealed soil-the sum of which equals the precipitation depth per square meter.
Urban areas contribute to decreased water capacity in catchment, not only by the decreasing of soil water storage-as is traditionally declared. Total water loss is far greater and comprises not only the loss of soil-retention potential, but runoff and evaporation must be included. This water loss is rapidly changing over space and time. In consequence in this catchment a strong variability of flow regime is visible (decreased low flow and increased high flow). In the case of the Prądnik basin, the lack of significant trend (decrease trend) of low flow is caused by the high percentage of natural land use in the upper part of the catchment.
The European Union's "no net land take" initiative is a systematic way to eliminate the adverse effects of further urbanization. According to the initiative, each new area of land would have to be compensated for by the revitalization of the previously urbanized area while brownfields should be used as a priority for any new construction. An alternative to standard technical solutions is the concept of green and blue infrastructure based on the strengthening of natural processes in the landscape. The main advantage of green and blue infrastructure compared to gray infrastructure is its versatility, multifunctionality and ability to self-regulate. Even though academic ecological approaches highlight the significance of the functional aspects of ecosystems, applied management guidelines have still not been able to introduce them into practice. The restoration of ecosystem functions is frequently the consequence of the revitalization actions performed and is rarely intended as a primary purpose.
Our results showed how progressive urban development has affected the hydrologic cycle in the Prądnik river basin. One of the greatest advantages of this presented research is the methodology used and data availability. Altogether, this study provides a better understanding of the impact of urbanization on the temporal and spatial distribution of water resources. It could serve as a base for future smarter urban development, especially in the time of general drought predicted for Europe due to LU/LC change and climate change.