Next Article in Journal
Spray Characteristics and Parameter Optimization of Orifice Arrangement for Micro-Sprinkling Hoses
Previous Article in Journal
Zooplankton Compositions in the Danjiangkou Reservoir, a Water Source for the South-to-North Water Diversion Project of China
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Capacity of River Valleys to Retain Nutrients from Surface Runoff in Urban and Rural Areas (Southern Poland)

1
Institute of Nature Conservation, Polish Academy of Sciences, Al. Mickiewicza 33, 31-120 Kraków, Poland
2
Department of Land Reclamation and Environmental Development, Faculty of Environmental Engineering and Land Surveying, University of Agriculture in Krakow, Al. Mickiewicza 24–28, 30-059 Krakow, Poland
*
Author to whom correspondence should be addressed.
Water 2022, 14(20), 3259; https://doi.org/10.3390/w14203259
Submission received: 6 September 2022 / Revised: 7 October 2022 / Accepted: 10 October 2022 / Published: 15 October 2022
(This article belongs to the Section Hydrology)

Abstract

:
Studies on water quality are necessary, as catchments of small watercourses are exposed to anthropogenic influences associated with agricultural activities, settlement, transport and other undertakings, leading to water pollution. There has been insufficient research performed on the valley’s ability to retain nutrients during floods, contributing to water accumulation. The main object of the study was to identify the retention capacity of river valleys under various aspects of human urbanization. To represent soil water retention, the Soil Conservation Service Curve Number (SCS-CN) method was used. Spatiotemporal autoregressive models were exploited to investigate the relationship between pollutants in precipitation and surface water in rivers. In contrast, multivariate analysis was used to identify and reveal patterns of land use for specific chemical compounds in the headwaters. The canonical-correlation analysis (CCA) showed that Mg+2 and Ca+2 cations in rainwater and surface waters play the main roles in the geochemical cycle in urban and rural areas. In the urban catchment area, the strongest relations were found for NO3, K+ and Na+. The average NO3 concentration in urban headwater was 8.3 mg·dm−3, the highest in the study area. The relationship between NO3 concentration in headwater and rainwater was found for all study catchments using spatial autoregression (SAR). High concentrations of SO42− in surface water have been identified in urban areas. Severe water erosion raises the risk of nutrient leaching in soils prone to surface runoff. As a consequence of low soil permeability and urbanization, retention capacity is significantly reduced in areas with low soil permeability. Land development plans should take spatial retention capacity into consideration. To ensure that large reservoirs can retain water in the face of climate change, riparian buffer zones (protective zones in valleys for small water bodies as well as Nature-based Solution) are important.

Graphical Abstract

1. Introduction

Recently, the interest in the retention capacity of catchments has grown significantly [1]. This results from the increasing frequency of floods and their negative consequences, resulting in economic and social losses [2]. Thus, it is important to note the retention capacity of a catchment [3]. The soil and its filtration capacities, i.e., the permeability of the soil or the surface layer of the soil, represent one of the most important factors in this respect. The type of soil may significantly influence the retention potential, as its ability to absorb water plays a particularly important role [4,5]. The role of spatial planning in catchment retention capacity has gained increasing recognition. It is emphasised in various programme documents and legal acts, both at the national and European levels [6]. This issue is gaining in importance, especially in the context of the transformation of arable lands into built-up and urban areas, and this significantly influences the increase in surface runoff and reduces water absorption by the ground [7].
An improvement in the soil structure improves the soil retention capacity. Appropriate agricultural practises, the use of organic fertilisers and other activities increase the soil water resources [8]. Agromelioration practises are also important, as they improve the soil structure and eliminate layers of poor permeability when performed correctly. The agromelioration practises have the following double role: they increase the soil retention capacity and facilitate water penetration into the aquifer [9]. The water parameters of the soil significantly influence the retention parameters of the substrate. Rainwater may move within the soil; however, this depends on many factors. Porosity is influenced by various factors, such as moisture content, precipitation, methods of use and farming practises [10]. Absorption of rainwater means mainly vertical water movement [11]. The initial phase of water absorption involves soil wetting [12].
Researchers all over the world study issues associated with water quality in large and small areas. The management of water resources is decisive for sustainable development [13,14]. One of the methods to retain water in an urban catchment is through activities termed ‘small retention measures’, which are mainly associated with the construction of dams on local watercourses and the creation of water reservoirs [15]. The aim of measures supporting small retention is to increase the water volume in the soil through technical (agricultural, agromelioration and phytomelioration) activities and water retention by plants [16]. Apart from storage of water for direct use, small retention aims at regulating and controlling water circulation in the environment [17]. It facilitates better protection and renewal of water resources. It also enables reliable water management without disrupting the balance in nature and develops urban resilience [18]. Reservoirs of dams play a very important role in water management because they help to reduce flood waves in urbanised areas [19], while they can be used to supply water to a river during a drought [20]. Furthermore, small water reservoirs are very popular in bare, arid areas where water for drinking, household use and irrigation in agriculture is particularly scarce [21,22]. In many countries, small retention tanks are very useful in rural areas because they gather water that is later used in agriculture, mainly to irrigate crop areas [23].
Surface waters in micro and small catchments are studied by research centres, which determine the effect of human activities on the quality of water resources in previously unmonitored catchments [24]. On a global scale, changes that took place in land cover and use are named as the cause of poor quality of water [25]. An increase in the area of arable land and pastures, as well as an increasing number of paved non-permeable surfaces, with a simultaneous decrease in the area of green lands and forests, resulting in a global reduction in evapotranspiration [26], with a simultaneous increase in the surface runoff [27,28]. This had a significant influence on the quality of surface waters, as their salinity increased and parameters deteriorated [29,30] due to fertilisers used in agriculture and contaminations flowing from urbanised areas [31].
One of the methods to solve complicated problems of water management may be various environmental, technical and organisational actions supporting water accumulation, delaying runoff and increasing its availability to the economy, including agricultural production and land development, such as the fashionable ‘Sponge City’ concept [32]. Recently, nature-based solutions have been particularly popular [33,34]. In this respect, geographic information systems and remote sensing procedures are helpful [35], enabling the description of catchment characteristics and integrating data within one system [36]. Various programmes are proposed to improve the quality of water through correct development and management of river catchments [37]. Therefore, the paper verified the physical and chemical parameters of a surface watercourse in urban and rural areas, and of rainwater. It is very important to ensure that the supply of a sufficient quantity of water is combined with its appropriate quality, as polluted water may cause the degradation of a water reservoir within a short time [38]. An analysis of the spatial retention capacity of three small catchments used in different ways is very important in terms of current policies related to climate change [39]. Cost-effective engineering solutions can reduce flood risk by optimising drainage and detention [40]. In response to climate change, it has become necessary to improve the water retention capacity of catchments and protect valleys [41].
Areas with varying degrees of anthropogenization were calculated based on catchment characteristics and soil type. Nutrient leaching into small catchments needs to be quantified to ensure the forest environment’s productivity. For drained forest catchments, it can be implemented, for instance, through improved water runoff measures. A rural area’s social interests should be managed and protected in the most effective way. Insufficient research has been conducted on the valley’s ability to retain water during floods, contributing to water accumulation overall. A vegetative buffer zone should be established around local water supply networks to prevent valley retention. Consequently, the following research objectives were established: (i) measuring the effect of anthropopressure on water quality in a catchment area; (ii) identifying soil retention rates in urban, suburban and rural areas; (iii) calculating the optimal balance between retention capacity and drainage of river valleys; (iv) conducting spatial analysis of rainfall and surface water quality.

2. Materials and Methods

2.1. Research Area

The studied catchments are located in three different regions of the Małopolskie Voivodeship, the Western Carpathians province, and the Outer Western Carpathians subprovince [42]. A catchment of the urban watercourse (6.69 km2) is located in the direct vicinity of the northern border of the city of Kraków in the Wyżyny Polskie province. A suburban catchment (of 7.21 km2) is located to the southeast of Tarnów. A rural area catchment (4.92 km2) is in the southwest part of Małopolska, away from large urban complexes (Figure 1).
Sub-urban catchment is located in the Central Beskidian Piedmont (Pogórze Środkowobeskidzkie) macroregion, the Ciężkowice Piedmont (Pogórze Ciężkowickie) mesoregion, which stretches between valleys of two rivers. The Cretaceous-Ciężkowice sandstones and conglomerates are distinguished by their resistance amongst rock series of the Silesian nappe, forming west-east folds [42].
The urban catchment is located in the south-western part of the Nida Basin (Niecka Nidziańska) macroregion and the Proszowicki Plateau (Płaskowyż Proszowicki) mesoregion. Its catchment is a low-raised upland (up to 250 m a.s.l.) of Miocene rocks (clays and sands) covered with a relatively thick layer of loess. This area, with numerous flat hills, is divided by lower sections of Nidzica, Szreniawa and Dłubnia valleys. Appropriate black soils and advantageous morphological conditions make the Proszowicki Plateau one of the most fertile regions of Poland, with very good conditions for farming and agricultural production [42]. The geological substrate in the study area of the rural catchment is formed by sedimentary rocks. The highest parts of this catchment are formed by flysch rocks of the Magura nappe, of which the southern part, similar to this entire region, is made of undulating, alternating and resistant to denudation layers of sandstones and conglomerates, with interbeddings of less resistant argillaceous schists and, occasionally, marly slates. Hollows in the land are scoured in siltstones and argillaceous schists, more susceptible to erosion [43].

2.2. Examined Indices of Surface Water and Rainwater Quality

In the selected catchments, regular monthly analyses of water quality indicators were conducted over a four-year period. In the studied brooks, in the terrain directly in the cross-section of dams and in planned reservoirs, water temperature, specific conductivity, pH, dissolved oxygen and oxygen saturation were determined on a randomly selected day of the month. The majority of indices were determined once a month, and only heavy metals (CrTot., Zn2+, Cd2+, Cu2+, Ni2+, Pb2+) once a quarter, according to the Polish departmental regulation. Rainwater samples were taken once a month on the same day as surface water samples. A trap container near the planned water reservoir was used to collect this water. Rainwater collection traps use a tray with a hole in the bottom were joined by a hopper.

2.3. Apparatus and Methods for Determination of Studied Indices

Water temperature, dissolved oxygen and water saturation with oxygen were measured directly in field using the apparatus of CO-411 type. Its pH was measured with a pH-meter CP-104, and specific conductivity with an electrical conductivity meter CC-102.
Parameters determined in the laboratory included total suspended solids by the oven test; dissolved substances by evaporation; SO42− level by precipitation; BOD5 by the Winkler method; COD by the permanganate method [44] and levels of Ca2+, Na+, K+, Mg+, Mn2+, CrTot., Zn2+, Cd2+, Cu2+, Ni2+, Pb2+ and FeTot by the atomic absorption spectrometry method using the spectrometer UNICAM SOLAR 969. Biogenic indicators such as ammoniacal (N–NH4+), nitrite (N–NO2) and nitrate (N–NO3) nitrogen, PO43− and Cl were determined by the flow colourimetric assay using FIAstar 5000. The results of N–NH4+, N–NO2 and N–NO3 determinations were used to calculate levels of NH4+, NO2 and NO3. Coliform and faecal coliform counts were determined on media containing lactose, following incubation at 37 and 44 °C, to within ±0.5 °C.

2.4. Land Use and Development

Maps of land cover were drawn on a basis of 1:2000 scale cadastral maps and 1:13,000 scale ortho-photos purchased at the Provincial Centre for Geodetic and Cartographic Documentation, which were then verified during onsite visits. Hypsometric and slope maps were drawn on the basis of a digital terrain model (DTM). Obtained DTM data, in form of triangulated irregular networks (TIN) created with a spatial precision corresponding to a 1:10,000 scale map, were transformed to a raster file in the form of a regular square GRID, of a resolution of 20 m/pixel. The study was divided into stages, of which the first includes the literature review, describing soil retention capacity, modelling of surface runoff and a rainfall-runoff model. Then the following three studied areas: Rygliczanka brook, Osielczyk and Sudół Dominikański catchments were characterised. The characteristics of the study area include descriptions of the catchment’s administrative and geographical locations, geology and geomorphology, climate and weather conditions, catchment hydrography, soil map and land cover and use. The final stage presents test results. The study used the SCS-CN method [45], for which map studies were prepared using QGIS 3.8. During development of features at the study sites, spatial analyses and creation of maps, the GIS software: ArcGIS 10, MapInfo Professional 8.0, Surfer 8, Erdas Image 8.7 and Corel Draw Graphics Suite X3 were used.

2.5. Hydromorphological Evaluation of Fresh Waters for Nature-Based Solution

A hydromorphological evaluation of running waters was performed according to the River Habitat Survey (RHS) method adapted from the Frame Water Directive (FWD) [46]. The illustrative picture shows the hydrological processes studied, including catchment retention, detention and drainage of soil water in the valley and nutrient flux to surface headwater and river (Figure 2).
The studies performed according to the RHS system allowed for the collection of ca. 400 parameters determining hydromorphological riverbed conditions were collected in 10 control profiles located every 50 m, and in a synthetic description of elements not registered in profiles. The evaluation took place at the location of the planned dam, and a five-hundred-meter-long section was examined upstream along the brook bed. The qualitative parameters describing morphological features were analysed and then used as a basis for the calculation of the following two synthetic watercourse quality indices for all studied sections:
(1)
The Habitat Quality Score (HQS), based on a presence and diversity of natural elements of the course and the river valley;
(2)
The Habitat Modification Score (HMS) determining the range of modifications in the watercourse morphology [47].
Values of the Habitat Quality Score, HQS, for a given section are calculated on a basis of a sum of partial values for the following categories listed below:
  • Flow type (waterfall, spill, boil, torrential, chaotic, rapid, rising, smooth, invisible, or dry bed);
  • Bed bottom material (rock outcrop, boulders, stones, small stones/gravel, sand, mud, clay/loess, peat/muck, concrete, mesh and stone gabions, claddings and pavings, riprap, or synthetic covering);
  • Natural morphological elements of the bed (rock outcrop, exposed boulders, rock outcrops/boulders covered with vegetation, a central bar not stabilised by vegetation, a central bar stabilised by vegetation, an island, a natural dam);
  • Natural morphological elements of banks (lateral erosion, stable bank erosion, a meander fluvial bar not stabilised by vegetation, a meander fluvial bar stabilised by vegetation, a point bar not stabilised by vegetation, a point bar stabilised by vegetation and natural embankment);
  • Structure of bank vegetation (none, uniform, simple, or complex);
  • Meander fluvial bars (stabilised and not stabilised by vegetation);
  • Groups of water plants (liverworts and mosses, emergent broad-leaf plants, emergent narrow-leaf plants, submerged plants with floating leaves, free-floating plants and plants rooted on the bank with stems floating on the water, submerged broad-leaf plants and submerged plants with narrow and with strongly lobed leaves),
  • The use of the land within a belt of 50 m from the bank top (deciduous/mixed forests, coniferous forests, wetlands);
  • Tree stands and elements morphologically accompanying them (isolated/dispersed, regularly distributed, continuous and semi-continuous);
  • Valuable elements of the river environment (waterfall, side channel, leaf heaps, natural reservoir, reed bed, peat bog, or bog).
  • The value of the Habitat Modification Score, HMS, is based on a sum of partial values for the following categories listed below:
  • Modifications observed in control points (bank stabilisation; bed stabilisation; bank or bed profiling; braided channel; an embankment on a river bank slope; a culvert; a damming structure; a crossing; a bank trampled by livestock);
  • Structures not observed in control profiles (a pedestrian crossing; a road or railway bridge; a groyne; a damming structure; a crossing; a culvert);
  • Modifications observed during a synthetic evaluation not registered in the control profiles (bed bottom material of anthropogenic origin; a stabilised entire bank profile; a profiled bank; a sectional bank profile; an embankment on a river bank slope; an embankment outside a river bank slope; plants removed from the bed; bank mowing) [48].

2.6. Statistical Analyses

Statistical analyses were performed in STATISTICA 9 [49]. An arithmetic mean and a median, as well as minimum and maximum values were calculated for all indices studied in surface and rain waters, and the sample sizes were specified. Three tests: Student’s t-test, Cochran-Cox test and non-parametric U Mann–Whitney (Wilcoxon) test were used for analyses of the seasonal variability in structures in surface and rain waters.
When at least one data set (sample) of a given indicator from the summer or the winter half of the year was not characterised by a normal distribution, the U Mann–Whitney test was used for analyses, as it is considered a non-parametric equivalent of the Student’s t-test. When both samples were characterised by the normal distribution, the hypothesis about homogeneity of variance of relevant populations was verified. As in all cases the size of the studied sample was always ≤50, the z-test was not used. When variances were equal, the Student’s t-test was used; otherwise, the Cochran-Cox test was selected. The level of significance α = 0.05 was assumed for all statistical analyses.
In all cases, also in the case of inter-object analyses, the normality of distribution of studied samples was evaluated through analyses of histograms and the Kolmogorov–Smirnov test, which was the least powerful test verifying normality, versus the Shapiro–Wilk and Lilliefors tests. The Kolmogorov–Smirnov test was selected as a specific compromise between a will to examine the significance of differences with parametric tests (more powerful than non-parametric ones), and requirements of these tests assuming a normal distribution of the analysed parameter.
In all cases, the homogeneity of variance was verified with the Brown–Forsythe test, which conducts an analysis of variance (ANOVA) for each variable on a basis of the absolute deviation of a value from its individual median. The Brown–Forsythe test determined the homogeneity of variance. In multiple comparisons, a one-way ANOVA was used. When it established that there was a significant difference between the analysed objects, Tukey’s post-hoc test was performed, which precisely indicated the objects between which that significant difference existed. When at least one of conditions for normality of the distribution in a given sample or for the homogeneity of variance was not met, then the Kruskal–Wallis Nonparametric ANOVA on Ranks was performed, which indicated significance of differences and led to multiple comparison of mean ranks for all samples. The PCA (Principal component analysis) was also used to demonstrate the most important physical and chemical variables in surface waters using software PQ stat ver. 1.6.8. The canonical correspondence analysis (CCA) was applied to compare physical and chemical parameters in rainwater and surface waters for each catchment, separately and independently in Cancoco for Windows ver. 4.51. Spatial relations between studied catchments and physical and chemical parameters of surface and rain waters were demonstrated on a basis of the simultaneous autoregressive model (SAR), using the Spatial Analysis in Macroecology (SAM; version 4.0) software [50].

3. Results

3.1. Physical and Chemical Properties of Surface Waters

The highest values for levels of the studied indices were noted in rainwater from the urban area. The inter-object analysis of the significance of differences in surface water quality indices showed that the highest differences were found between the waters of the mountains and the urban areas. No significant differences between objects were found for water temperature, total suspended solids and heavy metal levels (Table 1). The temperature of the water running off in the winter half of the year was significantly lower than that in the summer half of the year, and this results from climate conditions. This led to significantly higher levels of dissolved oxygen during the winter. Nitrite, phosphate, dissolved substances, magnesium and potassium levels and specific conductivity values were significantly higher in the summer, while nitrate levels were significantly higher in the winter. Seasonal differences in the group of metals were statistically significant only in a very few cases. In most cases, the object and seasonal variability were consistent with the results of inter-object comparisons. Significant differences were noted for pH, dissolved oxygen, COD-Mn, PO43−, Pb2+ and microbiological indices. The inter-object analysis of differences in values of indices in rainwater shows that the majority of them are statistically significantly higher in the urban area versus the suburban and rural areas (Table 2). PCA demonstrated the importance of physical and chemical parameters in each studied catchment. On the basis of PCA, the most important physical and chemical parameters of surface waters were selected for the spatial model. For selected parameters, the SAR model demonstrated that NO3 played a large role in the studied catchments, while in the urban catchment, the SO42− level was also an important factor (Table 3).
The hydrochemical analysis of running waters—RHS—adapted to the Polish conditions demonstrated that the suburban and the rural catchments can be classified as unmodified watercourses, while the urban catchment, with the HMS score of 28 points, can be called a partly modified course.

3.2. Catchment Retention Capacity

In the urban catchment area (Figure 3), the curve number (CN) ranged from 55 (forests on soils of permeability above average) to 98 (non-permeable areas). In the suburban catchment area, the calculated CN (Figure 4) ranged from 55 (forests on soils of permeability above average) to 98 (non-permeable areas). In the rural catchment area, the CN (Figure 5) assumed values from 77 (forests on soils of low permeability) to 98 (non-permeable areas).
The highest potential retention in the urban area was noted in the study points in forests (207.82 mm) and meadows (183.93 mm) growing on soils of permeability above average (Figure 6). The lowest potential retention was found in areas with non-permeable cover (5.18 mm). The areas of the lowest possible retention are located mainly in the western part of the catchment. The eastern and central parts of the catchment, except for small clusters, are characterised by much higher potential retention.
The highest potential retention in the suburban area was characterised by areas with forests (207.82 mm) and meadows (183.93 mm) growing on soils with a permeability above average (Figure 7). The lowest potential retention was found in areas with non-permeable cover (5.18 mm). The highest potential retention in the rural area (Figure 8) was observed in areas with forests (75.87 mm) and meadows (71.64 mm). The lowest potential retention was found in areas with non-permeable cover (5.18 mm).
In the urban area, surfaces with a moderate retention capacity (Figure 9) predominate (69.33%). The areas with very high, low, very low and high retention capacity have a much lower share (11.83%, 9.19%, 5.65% and 4%, respectively). In the suburban catchment, areas with very high retention capacity (Figure 10) have the highest share (55.24%). The areas qualified as characterised by high and very low retention capacity have the lowest share of the total surface area (1.43% and 1.55%, respectively). In the rural area, surfaces with a moderate retention capacity (Figure 11) strongly predominate (88.30%). The areas with low and very low retention capacity were also noted (8.51% and 3.19%, respectively). The areas of very high and high retention capacity were not confirmed.

4. Discussion

4.1. Relationship between a Land Use and Physical and Chemical Quality of Running Waters

As has already been mentioned, the excessive content of physical and chemical components in surface waters of small catchments also depends on the soil type, its pH, plant cultivation methods [51] and the intensity of erosion processes. This last factor is determined not only by the soil type and the way of land use but also by morphological characteristics of the catchment, such as land slopes, exposure, or slope length [52]. It also depends on the climate, especially the rainfall amount and intensity, and the temperature, with its low values causing freezing of the soil [53]. Water analyses at characteristic points along the entire length of the watercourse in the urban catchment between July and February were conducted by Kanownik and Rajda [54]. When analysing the indices (temperature, specific conductivity, dissolved substances, dissolved oxygen, NH4+, NO3, NO2, PO43−, SO42−, FeTot, Mn2+, Ca2+, Mg2+ and Cl), they noted an increase in twelve out of sixteen analysed parameters in the section where a wastewater discharge from a local treatment plant was located. The poor quality of water discharged from the sewage treatment plant was confirmed by own tests, as during their performance the sewage treatment plant in Węgrzce was closed because the local sewage system was connected to the Kraków-Płaszów treatment plant, so the values of the indices before and after closing of the treatment plant could be compared. It was found that the levels of the majority of the indices improved significantly, but levels of oxygen, NO3, Mg2+, Zn2+ and Ni2+ deteriorated (Table 1).
The worst quality of rainwater in the town located in the urban catchment can be associated within a close vicinity of the city of Kraków, and in particular, of its industrial district, Nowa Huta, with the operating large steelworks. The studies conducted in areas of northern France, also subjected to the influence of a similar plant manufacturing steel, conducted before and after its closure, demonstrated that the quality of rainwater significantly improved after the closing of the steelworks [55]. When compared to rainwater studied in large global agglomerations, the rainwater from our facilities was not acidic as, e.g., in China [56], Japan [57] or the United States [58].
Studies on the quality of water in a lowland river on the Szczecin sea coast, whose catchment was used mainly for agricultural purposes (no forests and a very small share of grasslands), showed that the nitrate (NO3) levels ranged from 0.9 to 42.0, with a mean level of 6.8 mg∙dm−3, ammonia (NH4+) levels ranged from 0.06 to 1.80 with a mean level of 0.32 mg∙dm−3, and phosphates ranged from 0.01 to 1.60 with a mean level of 0.35 mg∙dm−3. The mean potassium level was 13.4 mg∙dm−3 (ranging from 6.5 to 28.6) [59]. Comparing our results to those shown above, it can be said that only K+ levels in lowland rivers were higher than in all three sites; in all remaining cases, the values of biogenic indices (PO43−, NH4+, NO3) in a lowland river on the Szczecin sea coast were similar to levels noted in Węgrzce, yet they were always lower by at least 5%.
When our results are compared to studies of Olszewska and Krzemińska conducted in an agricultural catchment of the Jeziorka river, a right-bank tributary of the river Odra in the Malczyce region, it can be said that mean values of PO43−, NH4+, NO3, Ca2+ and Mg2+ levels were between mean values calculated for the catchments in the rural and urban areas, while the mean values of Mn2+, FeTot, SO42− and Cl levels in the Jeziorka river were much higher than levels noted in our monitoring (Table 1) [60]. The mean NO3 level was the highest in Rumia; in the remaining cases, it ranged from 0.9 to 3.2. The specific conductivity in all sites was not high and ranged from 19 to 52 µS∙cm−1; similarly, the following remaining indices: SO42−, Cl, Ca2+, Mg2+, Na+ and K+ were characterised by very low levels.
PCA demonstrated that TDS had the greatest influence on the studied variables (Figure 12). The total suspended solids (TSS) could enter the surface watercourse together with Na+ (Figure 13). In the catchment of the suburban watercourse, Na+ also had a high influence, but it did not necessarily enter the water together with suspended solids (TSS) as in the urban catchment described above. In the rural area, TSS played the main role in the surface watercourse. The remaining studied indices, especially Ca2+ and SO4−2, were characterised by significant correlations (Figure 14).
A hydrochemical evaluation of running waters, RHS adapted to the Polish conditions, was used as an assay supporting evaluation of water quality. Szoszkiewicz and Gebler (2011) conducted a statistical analysis on 950 sites all over Poland, analysed using the said method [61]. On the basis of the HMS indicator, they classified the rivers as modified (HMS > 45) and unmodified (HMS—from 0 to 2). The studies concerning the quality of water in small catchments using the same measurement methods and determination of physical and chemical indices, as in our study, were also conducted by [62]. They conducted their observations on watercourses flowing through variously used areas for two hydrological years. In those studies, pH, dissolved oxygen, NH4+, NO3-, NO2-, PO43-, SO42-, Cl-, Ca2+, Mg2+ and specific conductivity were determined in the surface water every month. The studies were conducted in microcatchments of the following: forest (75% forest, 13% arable land, 7% grasslands, 3% orchards, a few houses); agricultural (81% arable land, 7% grasslands, 5% orchards, 4% forests and tree stands, scarce dispersed houses), and settlement and agricultural (59% arable land, 19% forest, 8% grassland, 4% orchards, numerous dispersed houses) character. The CCA shows that Mg+2 and Ca+2 cations in rainwater and surface watercourses play the main roles both in urban and suburban areas (Figure 15). No large dependencies were noted in the rural areas. Additionally, NO3 influences rural areas. In the urban catchment area, the strongest relations were found for NO3, K+ and Na+. It has been shown that rainfall can have an effect on the biogeochemical cycle of NO3 (Figure 16). Overall, higher amounts of NO3 flow into the watercourse statistically significantly (Table 3). Watercourses in rural areas retain significant concentrations of NO3 as well as K+ (Figure 15).
Spatial autoregressive parameter (rho): 0.43
The land use and cover are of fundamental importance for the quality of water resources in the catchment. Long-term studies of water temperature, pH, BOD5, COD, dissolved substances, total phosphorus and total nitrogen, determined in water samples collected from 118 points in the Han River catchment (South Korea), showed that there is a correlation between the urban development of the land and water pollution [63]. In the catchment of the Little Miami River (USA), it was also found, using geographic information systems (GIS), that agricultural areas and paved non-permeable surfaces in urbanised areas cause much higher water pollution when compared to other ways of land use and cover [64]. Physical and chemical indices of surface water are correlated with farming activities [65]. This group consists of phosphates related to phosphate fertilisers, chlorides associated with potassium salt fertilisers, as well as nitrogen in its various forms of NH4+, NO2 and NO3 caused by the use of nitrogen fertilisers [66,67]. The partial influence of agriculture is emphasised by indices such as BOD5, COD and sulphates [68,69]. Studies in Slovakia have shown that with a decrease in field fertilisation with nitrogen compounds, its levels in brooks flowing through the studied areas decreased [70]. Swedish researchers studying 35 small (2–35 km2) agricultural catchments for five years noted a strong correlation between inorganic nitrogen levels and the way of land use [71].

4.2. Hydrochemical Evaluation of Rainwater

CCA has shown that there are relationships between the chemical composition of rainwater and surface waters. In the mountain catchment, the content of Mg2+, Ca2+ and NO3 should be noted (Figure 16). In the rural catchment, additionally, Na and Cl levels are of importance, while in the urban catchment, NO3, K+ and Na+ should also be noted. Many factors contribute to the quality of surface waters [72,73], and their influence on excessive levels of physical and chemical components in those waters is different and varies over time [74]. In small catchments, the anthropogenic influence on the quality of surface waters is partly combined with the effect of natural components such as rocks, soil humus substances [75], and, to some extent, precipitation [76]. The main source of pollution is of anthropogenic origin [77]. A deterioration in the quality of water in the urban catchment is associated with human existence and economic activities [78], including incorrect water and wastewater management, excessive use of artificial and organic fertilisers, and air pollution [79].
When compared to our study, a two-year study by Rajda et al. (2001) on selected physical and chemical components of rainwaters in a region of a small forest catchment of the Jagódka brook located in an eastern part of the Little Beskids, as well as the study by Ostrowski et al. (2000) in three microcatchments located near the city of Andrychów, also in the Little Beskid, demonstrated that only in two cases mean values of the indices were higher in the urban catchment than mean levels noted in the Little Beskid, while the lowest mean values predominated in the urban catchment waters, and this confirms that rainwaters of the best quality were found in the area located the furthest from large agglomerations [80,81].
The SAR model has shown that NO3 in the urban and suburban areas and SO42− in the urban area may be of significance when determining the influence of rain on the physical and chemical quality of water in the catchment (Table 3). Even partial settlement development of land in small catchments causes greater water pollution than is the case in typical agricultural catchments [82]. Factors associated with settlements such as the number of inhabitants, the build-up area share, the arable land share, the number of livestock and quantities of mineral and organic fertilisers used within the catchment area have a significant influence on levels of certain indices of surface water quality [75,83].
It is difficult to compare the results of physical and chemical determinations conducted by different authors in other sites, as methods used for measurements and analyses of the results vary. Furthermore, complete information on studied sites is not available, and it is difficult to establish precisely which factors significantly influence many simultaneous processes decisive for water quality [84].

4.3. Modelling of Surface Runoff

Surface runoff is one of the processes associated with water circulation in the environment. It is that part of rainwater that flows on the surface of the ground. It occurs after intensive rainfalls or is the result of snow melting in spring [85]. The surface runoff speed depends on the land slope and coverage [86]. The integration models are of importance for the evaluation of thermal cooling in river revitalisation processes [87]. During intense rainfalls and a small resistance to water flow, the runoff intensifies and floods develop [88]. The results showed that issues associated with hydrological modelling mainly concern the transformation of the effective precipitation in the surface runoff. Land development must be planned to take into account soil permeability so that any excess rainwater can be removed. In areas covered by soils of low permeability, the share of non-permeable surfaces should be limited [89]. This can be achieved, for example, by using open slabs, using the greatest possible areas such as green belts, parks and green squares, and maintaining the existing natural plant coverage to the greatest extent possible. The CN parameter showed that, in terms of the rainwater-runoff relationship, urban soils (Figure 3) are haracterized by low permeability at the same low level as suburban (Figure 4) or rural (Figure 5) areas.
Soil permeability is a significant element influencing surface runoff. Soil permeability influences the absorption of rainwater and depends on the physical properties of a substrate [90]. On permeable soils, rainwater is absorbed by the soil, reducing surface runoff. In areas where non-permeable soils are found, water evaporates [91].
A very important problem that needs to be solved is how to prevent the pollution of running water. One of the ways to protect water is the natural structure of bank vegetation, which helps to maintain biological diversity and improve the condition of the natural, especially aquatic, environment [92]. The creation of buffer zones with forest tree stands has a positive effect, as they provide conditions for the development of flora and fauna, improving the ecological state of a watercourse [93]. Grasses also have a positive influence on the reduction of non-point pollution of agricultural origin. It is recommended to plan grasslands in the vicinity of watercourses and reservoirs to form a natural buffer to protect against an excessive supply of nutrients leading to water eutrophication as a nature-based solution [94]. Moreover, waters in courses with forests, grasslands and organic soils located along them were characterised by good quality, as they captured a significant part of non-point pollution running off from fields [95,96]. Rainwaters, reaching the soil surface, partly infiltrate it depends on the ground permeability, supplying groundwaters [97]. The maximum retention capacity of a catchment can depend on the way it is used and the distribution of trees and shrubs in the valley (Figure 6). The results of the study showed that areas with large groups of trees located on the borders of cities increase the retention capacity of the catchment. Rainwater accumulation and runoff were similar to agricultural catchment areas in rural areas (Figure 8). The slope is also an important factor, and the steeper the slope, the greater the water runoff [98]. When rain falls on nonurbanised areas covered with vegetation, the water soaks into the soil and slowly infiltrates into the surface and groundwater. In urbanised areas, where paved surfaces predominate, much of the rainwater ends up in canals and watercourses [99].

4.4. Catchment Retention Capacity

The majority of urban catchment areas are located on soils with above-average permeability. The soils in the suburban brook catchment were characterised by above-average permeability. Forests grew on low-permeable soils in the rural brook catchment, resulting in a moderately retained catchment area that was not affected by the drought. Catchment retention capacity is significantly influenced by the type of soil that is used. As non-permeable surfaces (buildings, roads, etc.) increase, catchment retention capacity declines.
The conducted studies showed that the way of the catchment area’s development significantly influences its retention capacity. The highest retention capacity was observed in forests and meadows, i.e., non-urbanised areas covered with natural vegetation (Figure 9). The more urbanised the catchment area is, the lower its maximum retention capacity. Additionally, the type of soil found in a given area is very important, as the low soil permeability significantly reduces the retention capacity (Figure 10).
Built-up, urbanised areas found on soils of low permeability have a very low retention capacity, and this intensifies the water runoff [100]. Urbanisation resulted in a reduction of the soil infiltration capacity in the catchment, as described by the CN parameter (Figure 3). This study contains observations indicating that the retention capacity may be a very important factor, which should be considered in the spatial planning process. The soil permeability study enables correct decisions concerning land cover, so the risk of disruptions in water runoff leading to floods in urban [101] and rural [102] areas is reduced.
The small retention (natural water retention measures), particularly in rural areas, aim at increasing water resources. One of the actions included in natural water retention measures may include procedures aiming at catching water running off from slopes of significant incline, as was the case in the Spanish mountains Sierra de Gador, where water flowing from slopes was harvested and directed to cisterns through a system of canals [103,104]. In the presented maps, areas where values of the CN parameters changed significantly due to water erosion and a reduction in water retention, as well as those related to urbanisation and the development of road infrastructure in recent years (Figure 4) and in rural areas (Figure 5) were identified.
Ponds, i.e., small hollows in the land where water accumulates, play an important role in natural water retention, as they support its circulation [105]. Apart from their landscape and environmental values contributing to areas where they can be found, they also increase water resources in their vicinity through a network of blue and green infrastructure [106]. In cities, urban wetland parks and specific natural areas should be created by using the ecohydrological approach or development of wetlands [107]. Variability in the catchment retention capacity plays an important role in the shaping of the surface runoff (Figure 6). The maximum retention capacity can be used to estimate the direct runoff level in catchments characterised by soils of low permeability (Figure 3). Due to deviations in the runoff size calculated according to NRCS, CN requires an adjustment to be taken into account, e.g., a correction for initial losses or water content conditions. In urban conditions (Figure 3), the variability in the formation of the surface runoff should take into account the local characteristics of the catchment.
Construction of reservoirs for limited retention should be preceded by analyses of water quality, as actions to improve water parameters should be initiated in catchments with water of poor or unsuitable quality. Taking this aspect of water quality into account when making decisions about the construction of water reservoirs will contribute to their multifunctional and optimal use. Maintaining good quality water in planned reservoirs for low retention will require proper operation, control over water and wastewater management, and the elimination of unlicensed waste landfills in the catchment areas. The analysis of water quality and functional parameters in the context of its long-term retention in reservoirs will enable us to specify how to use it for local purposes. Enhancing the catchment retention capacity of urbanised watercourses could be achieved by limiting the construction of entirely sealed surfaces and by increasing biologically active surfaces. The results of this study may also contribute to the development of guidelines for establishing new nature conservation sites in urban riverside parks. supporting a process of planning blue and green infrastructure in the cities and designing water reservoirs for small retentions in rural areas.

5. Conclusions

Soil drainage affects nutrient leaching in all studied catchments. A variety of land-use practises have also contributed to the retention of nutrients in the river valley. Urbanised areas and valley drainage were associated with NO3 concentrations in surface water based on the spatial analysis of the autoregressive SAR model. Results of CCA indicated that rainwater and surface watercourses accumulating Mg2+ and Ca2+ tend to be found in urban, suburban and rural areas. A rural area’s land use may also influence NO3 concentrations in surface water. The strongest correlations were revealed in the urban catchment for NO3, K+ and Na+ concentrations in surface waters. The quality of nutrients in the catchment was affected by varying levels of anthropogenic pressure due to land use and development. It was found that there was no statistically significant difference in most of the rainwater indicators tested across all catchments. At each catchment, significant differences were displayed only for pH, Na+ and FeTot. The pH values were lower, and Na2+ and FeTot values were higher in the summer half of the year. Further, periodic soil and valley susceptibility research is needed to address climatic change’s effects on catchments’ capacity to retain nutrients.

Author Contributions

Conceptualization, T.S. and W.H.; methodology, T.S.; software, T.S.; validation, W.H., T.S. and W.F.; formal analysis, W.H.; investigation, W.F.; resources, T.S.; data curation, W.H.; writing—original draft preparation, T.S.; writing—review and editing, W.F.; visualization, T.S.; supervision, W.H.; project administration, W.H. All authors have read and agreed to the published version of the manuscript.

Funding

The research was supported by the Ministry of Education and Science with funds granted to the University of Agriculture for 2022.

Data Availability Statement

Data is contained within the article.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. La Rosa, D.; Pappalardo, V. Planning for spatial equity—A performance based approach for sustainable urban drainage systems. Sustain. Cities Soc. 2020, 53, 101885. [Google Scholar] [CrossRef]
  2. Goh, K. Urban Waterscapes: The Hydro-Politics of Flooding in a Sinking City. Int. J. Urban Reg. 2019, 43, 250–272. [Google Scholar] [CrossRef]
  3. Kelleher, C.; Golden, H.E.; Burkholder, S.; Shuster, W. Urban vacant lands impart hydrological benefits across city landscapes. Nat. Commun. 2020, 11, 1563. [Google Scholar] [CrossRef] [Green Version]
  4. Halecki, W.; Stachura, T. Evaluation of soil hydrophysical parameters along a semiurban small river: Soil ecosystem services for enhancing water retention in urban and suburban green areas. Catena 2021, 196, 104910. [Google Scholar] [CrossRef]
  5. Walton, C.R.; Zak, D.; Audet, J.; Petersen, R.J.; Lange, J.; Oehmke, C.; Wichtmann, W.; Kreyling, J.; Grygoruk, M.; Jablonska, E.; et al. Wetland buffer zones for nitrogen and phosphorus retention: Impacts of soil type, hydrology and vegetation. Sci. Total Environ. 2020, 727, 138709. [Google Scholar] [CrossRef] [PubMed]
  6. Sylla, M.; Solecka, I. Highly valued agricultural landscapes and their ecosystem services in the urban-rural fringe—An integrated approach. J Environ. Plann. Manag. 2020, 63, 883–911. [Google Scholar] [CrossRef]
  7. Sylla, M.; Hagemann, N.; Szewranski, S. Mapping trade-offs and synergies among peri-urban ecosystem services to address spatial policy. Environ. Sci. Policy 2020, 112, 79–90. [Google Scholar] [CrossRef]
  8. Zheng, W.; Shen, C.; Wang, L.P.; Jin, Y. An empirical soil water retention model based on probability laws for pore-size distribution. Vadose Zone J. 2020, 19, e20065. [Google Scholar] [CrossRef]
  9. Wu, Y.Y.; He, G.J.; Ouyang, W.; Huang, L. Differences in soil water content and movement drivers of runoff under climate variations in a high-altitude catchment. J. Hydrol. 2020, 587, 125024. [Google Scholar] [CrossRef]
  10. McFarland, A.R.; Larsen, L.; Yeshitela, K.; Engida, A.N.; Love, N.G. Guide for using green infrastructure in urban environments for stormwater management. Environ. Sci. Water Res. Technol. 2019, 5, 643–659. [Google Scholar] [CrossRef]
  11. Uygun, B.S.; Albek, M. Determination of Climate Change Effects of Impervious Areas in Urban Watershed. Water Air Soil. Pollut. 2020, 231, 1–12. [Google Scholar] [CrossRef]
  12. Wang, P.J.; Zheng, H.F.; Ren, Z.B.; Zhang, D.; Zhai, C.; Mao, Z.X.; Tang, Z.; He, X.Y. Effects of Urbanization, Soil Property and Vegetation Configuration on Soil Infiltration of Urban Forest in Changchun, Northeast China. Chin. Geogr. Sci 2018, 28, 482–494. [Google Scholar] [CrossRef] [Green Version]
  13. Yu, Q.; Wu, C.; Qiao, S.; Wu, B.; Wang, J. The Application of Watershed System Strategy on Urban Water Environment Treatment. J. Water Resour. Res. 2019, 08, 191–198. [Google Scholar] [CrossRef]
  14. Mgelwa, A.S.; Hu, Y.L.; Ngaba, M.Y. Patterns of nitrogen concentrations and their controls in two southern China urban river ecosystems. Glob. Ecol. Conserv. 2020, 23, e01112. [Google Scholar] [CrossRef]
  15. Teurlincx, S.; Kuiper, J.J.; Hoevenaar, E.C.M.; Lurling, M.; Brederveld, R.J.; Veraart, A.J.; Janssen, A.B.G.; Mooij, W.M.; Domis, L.N.D. Towards restoring urban waters: Understanding the main pressures. Curr. Opin. Environ. Sustain. 2019, 36, 49–58. [Google Scholar] [CrossRef]
  16. Shifflett, S.D.; Newcomer-Johnson, T.; Yess, T.; Jacobs, S. Interdisciplinary Collaboration on Green Infrastructure for Urban Watershed Management: An Ohio Case Study. Water 2019, 11, 738. [Google Scholar] [CrossRef] [Green Version]
  17. Wise, J.L.; Van Horn, D.J.; Diefendorf, A.F.; Regier, P.J.; Lowell, T.V.; Dahm, C.N. Dissolved organic matter dynamics in storm water runoff in a dryland urban region. J. Arid Environ. 2019, 165, 55–63. [Google Scholar] [CrossRef]
  18. McClymont, K.; Fernandes Cunha, D.G.; Maidment, C.; Ashagre, B.; Vasconcelos, A.F.; De Macedo, M.B.; Nobrega dos Santos, M.F.; Gomes Junior, M.N.; Mendiondo, E.M.; Barbassa, A.P.; et al. Towards urban resilience through Sustainable Drainage Systems: A multi-objective optimisation problem. J. Environ. Manag. 2020, 275, 111173. [Google Scholar] [CrossRef]
  19. Verol, A.P.; Battemarco, B.P.; Merlo, M.L.; Machado, A.C.M.; Haddad, A.N.; Miguez, M.G. The urban river restoration index (URRIX)—A supportive tool to assess fluvial environment improvement in urban flood control projects. J. Clean. Prod. 2019, 239, 118058. [Google Scholar] [CrossRef]
  20. Xiong, W.; Li, Y.; Pfister, S.; Zhang, W.L.; Wang, C.; Wang, P.F. Improving water ecosystem sustainability of urban water system by management strategies optimization. J. Environ. Manag. 2020, 254, 109766. [Google Scholar] [CrossRef]
  21. Reyes, B.; Hogue, T.S.; Maxwell, R.M. Urban irrigation in the modeling of a semi-arid urban environment: Ballona Creek watershed, Los Angeles, California. Hydrol. Sci. J. 2020, 65, 1344–1357. [Google Scholar] [CrossRef]
  22. Tang, J.F.; Wang, W.D.; Yang, L.; Qiu, Q.L.L.; Lin, M.X.; Cao, C.L.; Li, X.H. Seasonal variation and ecological risk assessment of dissolved organic matter in a peri-urban critical zone observatory watershed. Sci. Total Environ. 2020, 707, 136093. [Google Scholar] [CrossRef]
  23. Xia, F.; Zhang, C.; Qu, L.Y.; Song, Q.J.; Ji, X.L.; Mei, K.; Dahlgren, R.A.; Zhang, M.H. A comprehensive analysis and source apportionment of metals in riverine sediments of a rural-urban watershed. J. Hazard. Mater. 2020, 381, 121230. [Google Scholar] [CrossRef]
  24. Ruan, X.F.; Huang, J.Y.; Williams, D.A.R.; Harker, K.J.; Gergel, S.E. High spatial resolution landscape indicators show promise in explaining water quality in urban streams. Ecol. Indic. 2019, 103, 321–330. [Google Scholar] [CrossRef]
  25. Sinha, J.; Jha, S.; Goyal, M.K. Influences of watershed characteristics on long-term annual and intra-annual water balances over India. J. Hydrol. 2019, 577, 123970. [Google Scholar] [CrossRef]
  26. Shao, Z.F.; Fu, H.Y.; Li, D.R.; Altan, O.; Cheng, T. Remote sensing monitoring of multi-scale watersheds impermeability for urban hydrological evaluation. Remote Sens. Environ. 2019, 232, 111338. [Google Scholar] [CrossRef]
  27. Favreau, G.; Cappelaere, B.; Massuel, S.; Leblanc, M.; Boucher, M.; Boulain, N.; Leduc, C. Land clearing, climate variability, and water resources increase in semiarid southwest Niger: A review. Water Resour. Res. 2009, 45, L16402. [Google Scholar] [CrossRef] [Green Version]
  28. Rost, S.; Gerten, D.; Hoff, H.; Lucht, W.; Falkenmark, M.; Rockstrom, J. Global potential to increase crop production through water management in rainfed agriculture. Environ. Res. Lett. 2009, 4, 044002. [Google Scholar] [CrossRef]
  29. Peter, K.T.; Hou, F.; Tian, Z.Y.; Wu, C.; Goehring, M.; Liu, F.M.; Kolodziej, E.P. More Than a First Flush: Urban Creek Storm Hydrographs Demonstrate Broad Contaminant Pollutographs. Environ. Sci. Technol. 2020, 54, 6152–6165. [Google Scholar] [CrossRef]
  30. Sarkar, K.; Majumder, M. Application of AHP-based water quality index for quality monitoring of peri-urban watershed. Environ. Dev. Sustain. 2021, 23, 1780–1798. [Google Scholar] [CrossRef]
  31. Luo, L.L.; Mei, K.; Qua, L.Y.; Zhang, C.; Chen, H.; Wang, S.Y.; Di, D.; Huang, H.; Wang, Z.F.; Xia, F.; et al. Assessment of the Geographical Detector Method for investigating heavy metal source apportionment in an urban watershed of Eastern China. Sci. Total Environ. 2019, 653, 714–722. [Google Scholar] [CrossRef] [Green Version]
  32. Nguyen, T.T.; Ngo, H.H.; Guo, W.S.; Wang, X.C.C.; Ren, N.Q.; Li, G.B.; Ding, J.; Liang, H. Implementation of a specific urban water management—Sponge City. Sci. Total Environ. 2019, 652, 147–162. [Google Scholar] [CrossRef]
  33. Oral, H.V.; Carvalho, P.; Gajewska, M.; Ursino, N.; Masi, F.; Hullebusch, E.D.V.; Kazak, J.K.; Exposito, A.; Cipolletta, G.; Andersen, T.R.; et al. A review of nature-based solutions for urban water management in European circular cities: A critical assessment based on case studies and literature. Blue-Green Syst. 2020, 2, 112–136. [Google Scholar] [CrossRef] [Green Version]
  34. Rigotti, J.A.; Pasqualini, J.P.; Rodrigues, L.R. Nature-based solutions for managing the urban surface runoff: An application of a constructed floating wetland. Limnetica 2020, 39, 441–454. [Google Scholar] [CrossRef]
  35. Oliver, S.; Corburn, J.; Ribeiro, H. Challenges Regarding Water Quality of Eutrophic Reservoirs in Urban Landscapes: A Mapping Literature Review. Int. J. Environ. Res. 2019, 16, 40. [Google Scholar] [CrossRef] [Green Version]
  36. Liu, J.H.; Shao, W.W.; Xiang, C.Y.; Mei, C.; Li, Z.J. Uncertainties of urban flood modeling: Influence of parameters for different underlying surfaces. Environ. Res. 2020, 182, 108929. [Google Scholar] [CrossRef]
  37. Parajulee, A.; Wania, F.; Mitchell, C.P.J. Hydrological transit times in nested urban and agricultural watersheds in the Greater Toronto Area, Canada. Hydrol. Process. 2019, 33, 350–360. [Google Scholar] [CrossRef]
  38. Park, H.; Beom, J.; Jeung, M.; Choi, W.; Her, Y.; Shirmohammadi, A.; Yoon, K.S. Identifying feasible nonpoint source pollutant sampling intervals for watersheds with paddy field and urban land uses. Water Supply 2020, 21, 780–790. [Google Scholar] [CrossRef]
  39. Porse, E.; Pincetl, S. Effects of Stormwater Capture and Use on Urban Streamflows. Water Resour. Manag. 2019, 33, 713–723. [Google Scholar] [CrossRef]
  40. Paredes, I.; Otero, N.; Soler, A.; Green, A.J.; Soto, D.X. Agricultural and urban delivered nitrate pollution input to Mediterranean temporary freshwaters. Agric. Ecosyst. Environ. 2020, 294, 106859. [Google Scholar] [CrossRef]
  41. Pistocchi, A. A preliminary pan-European assessment of pollution loads from urban runoff. Environ. Res. 2020, 182, 109129. [Google Scholar] [CrossRef]
  42. Kondracki, J. Geografia Regionalna Polski; Wydawn Naukowe PWN: Warsaw, Poland, 2009. [Google Scholar]
  43. Bednarek, R.; Prusinkiewicz, Z. Geografia Gleb; PWN: Warsaw, Poland, 1999. [Google Scholar]
  44. Hermanowicz, W.; Dojlido, J.R.; Dożańska, W.; Koziorowski, B.; Zerbe, J. Fizyczno-Chemiczne Badanie Wody i Ścieków; Wydawnictwo Arkady: Warsaw, Poland, 1999. [Google Scholar]
  45. Baiamonte, G. SCS Curve Number and Green-Ampt Infiltration Models. J. Hydrol. Eng. 2019, 24, 04019034. [Google Scholar] [CrossRef]
  46. Jusik, S.; Szoszkiewicz, K. Możliwości wykorzystania systemu River Habitat Survey (RHS) na potrzeby wdrażania Ramowej Dyrektywy Wodnej w Polsce. Possibilities of use of the River Habitat Survey to the needs of Water Framework Directive in Poland. Acta Sci. Polonorum. Form. Circumiectus 2010, 9, 23–34. [Google Scholar]
  47. Raven, P.; Holmes; Dawson, F.; Fox, P.J.A.; Everard, M.; Fozzard, K.J.R. River Habitat Quality: The Physical Character of Rivers and Stream in the UK and Isle of Man; River Habitat Survey Report No. 2; Environment Agency: Bristol, UK, 1998. [Google Scholar]
  48. Szoszkiewicz, K.; Jusik, S.; Lawniczak, A.E.; Zgola, T. Macrophyte development in unimpacted lowland rivers in Poland. Hydrobiologia 2010, 656, 117–131. [Google Scholar] [CrossRef] [Green Version]
  49. StatSoft. Elektroniczny Podręcznik Statystyki PL, Kraków. 2006. Available online: http://www.statsoft.pl/textbook/stathome.html (accessed on 30 September 2022).
  50. Rangel, T.F.; Diniz, J.A.F.; Bini, L.M. SAM: A comprehensive application for Spatial Analysis in Macroecology. Ecography 2010, 33, 46–50. [Google Scholar] [CrossRef]
  51. Kincaid, D.W.; Seybold, E.C.; Adair, E.C.; Bowden, W.B.; Perdrial, J.N.; Vaughan, M.C.H.; Schroth, A.W. Land Use and Season Influence Event-Scale Nitrate and Soluble Reactive Phosphorus Exports and Export Stoichiometry from Headwater Catchments. Water Resour. Res. 2020, 56, e2020WR027361. [Google Scholar] [CrossRef]
  52. Lepeska, T.; Wojkowski, J.; Walega, A.; Mlynski, D.; Radecki-Pawlik, A.; Olah, B. Urbanization-Its Hidden Impact on Water Losses: Pradnik River Basin, Lesser Poland. Water 2020, 12, 1958. [Google Scholar] [CrossRef]
  53. Aboelnour, M.; Gitau, M.W.; Engel, B.A. Hydrologic Response in an Urban Watershed as Affected by Climate and Land-Use Change. Water 2019, 11, 1603. [Google Scholar] [CrossRef] [Green Version]
  54. Kanownik, W.; Rajda, W. Źródła zanieczyszczenia wód powierzchniowych w zlewni potoku Sudół Dominikański. Sources of surface water pollution in Sudol Dominikanski stream catchment. Acta Sci. Pol. Form. Circumiectus 2008, 7, 3–14. [Google Scholar]
  55. Miler, A.; Liberacki, D.; Plewiński, D. Ilościowa i jakościowa ocena odpływu z dwu kontrastowych pod względem zalesienia mikrozlewni nizinnych. Zesz. Nauk. AR Kraków 2000, 365, 435–444. [Google Scholar]
  56. Yi, L.; Xiaolan, Y.; Hongbing, C.; Weili, L.; Jie, T.; Shufeng, W. Chemical characteristics of precipitation at three Chinese regional background stations from 2006 to 2007. Atmos. Res. 2010, 96, 173–183. [Google Scholar] [CrossRef]
  57. Ozeki, T.; Ihara, T.; Ogawa, N. Study of pollutants in precipitation (rain and snow) transported long distance to west coasts of Japan Islands using oblique rotational factor analysis with partially non-negative constraint. Chemom. Intell. Lab. Syst. 2006, 82, 15–23. [Google Scholar] [CrossRef]
  58. Butler, T.J.; Likens, G.E.; Vermeylen, F.M.; Stunder, B.J.B. The relation between NOx emissions and precipitation NO3 in the eastern USA. Atmos. Environ. 2003, 37, 2093–2104. [Google Scholar] [CrossRef]
  59. Durkowski, T.; Wesolowski, P. Kształtowanie się odpływu wody i zanieczyszczeń z małych zlewni rolniczych. Zesz. Probl. Postępów Nauk Rol. 2008, 528, 41–47. [Google Scholar]
  60. Olszewska, B.; Krzeminska, A. Jakość wód rzeki Jeziorki w latach 1995–2003. Jeziorka stream water quality in the years 1995–2003. Zesz. Probl. Postępów. Nauk. Rol. 2008, 528, 105–114. [Google Scholar]
  61. Szoszkiewicz, K.; Gebler, D. Ocena warunków hydromorfologicznych rzek w Polsce metodą River Habitat Survey. Ochr. Sr. Zasobów Nat. 2011, 47, 70–81. [Google Scholar]
  62. Misztal, A.; Kuczera, M. An Impact Assessment of Built-Up Residential Areas on Selected Water Quality Indexes. Pol. J. Environ. Stud. 2008, 17, 985–991. [Google Scholar]
  63. Chang, H. Spatial analysis of water quality trends in the Han River basin, South Korea. Water Res. 2008, 42, 3285–3304. [Google Scholar] [CrossRef]
  64. Liu, M.; Han, G.L.; Li, X.Q. Comparative analysis of soil nutrients under different land-use types in the Mun River basin of Northeast Thailand. J. Soils Sediments 2021, 21, 1136–1150. [Google Scholar] [CrossRef]
  65. Blaszczak, J.R.; Delesantro, J.M.; Urban, D.L.; Doyle, M.W.; Bernhardt, E.S. Scoured or suffocated: Urban stream ecosystems oscillate between hydrologic and dissolved oxygen extremes. Limnol. Ocean. 2019, 64, 877–894. [Google Scholar] [CrossRef] [Green Version]
  66. Puckett, L.J. Identifying the major sources of nutrient water-pollution. Environ. Sci. Technol. 1995, 29, 408A–414A. [Google Scholar] [CrossRef]
  67. Parry, R. Agricultural phosphorus and water quality: A US Environmental Protection Agency perspective. J. Environ. Qual. 1998, 27, 258–261. [Google Scholar] [CrossRef]
  68. Banasik, K.; Rudzka, E.; Smietanka, A. Eutroficzne wskaźniki zanieczyszczeń w odpływie rzecznym z małej zlewni rolniczej. Przegląd Nauk. Wydziału Melior. Inżynierii Sr. 1996, 11, 23–28. [Google Scholar]
  69. Kopeć, S.; Krzanowski, S. Ocena udziału zanieczyszczeń rolniczych Beskidu Wyspowego i Żywieckiego w zanieczyszczeniach rzek tego regionu. In Rolnicza Przestrzeń Produkcyjna Beskidu Wyspowego i Żywieckiego, a Jakość Wody; Wydawnictwo AR Kraków: Krakow, Poland, 1994; pp. 143–156. [Google Scholar]
  70. Pekarova, P.; Pekar, J. The impact of land use on stream water quality in Slovakia. J. Hydrol. 1996, 180, 333–350. [Google Scholar] [CrossRef]
  71. Arheimer, B.; Liden, R. Nitrogen and phosphorus concentrations from agricultural catchments—Influence of spatial and temporal variables. J. Hydrol. 2000, 227, 140–159. [Google Scholar] [CrossRef]
  72. Choo, Y.M.; Jo, D.J.; Yun, G.S.; Lee, E.H. A Study on the Improvement of Flood Forecasting Techniques in Urban Areas by Considering Rainfall Intensity and Duration. Water 2019, 11, 1883. [Google Scholar] [CrossRef] [Green Version]
  73. Balthazard-Accou, K.; Emmanuel, E.; Agnamey, P.; Raccurt, C. Pollution of Water Resources and Environmental Impacts in Urban Areas of Developing Countries: Case of the City of Les Cayes (Haiti). In Environmental Health—Management and Prevention Practices; Makan, A., Ed.; IntechOpen: London, UK, 2020. [Google Scholar] [CrossRef] [Green Version]
  74. Faridatul, M.I.; Wu, B.; Zhu, X.L. Assessing long-term urban surface water changes using multi-year satellite images: A tale of two cities, Dhaka and Hong Kong. J. Environ. Manag. 2019, 243, 287–298. [Google Scholar] [CrossRef]
  75. De Figueiredo, H.P.; De Figueiredo, C.R.P.; Barros, J.H.D.; Constantino, M.; Magalhaes, F.J.C.; De Moraes, P.M.; Da Costa, R.B. Water quality in an urban environmental protection area in the Cerrado Biome, Brazil. Environ. Monit. Assess. 2019, 191, 117. [Google Scholar] [CrossRef]
  76. Alamdari, N.; Sample, D.J.; Ross, A.C.; Easton, Z.M. Evaluating the Impact of Climate Change on Water Quality and Quantity in an Urban Watershed Using an Ensemble Approach. Estuaries Coasts 2020, 43, 56–72. [Google Scholar] [CrossRef]
  77. Bao, L.L.; Li, X.Y.; Su, J.J. Alteration in the potential of sediment phosphorus release along series of rubber dams in a typical urban landscape river. Sci. Rep. 2020, 10, 2714. [Google Scholar] [CrossRef] [Green Version]
  78. Blaszczak, J.R.; Delesantro, J.M.; Zhong, Y.; Urban, D.L.; Bernhardt, E.S. Watershed urban development controls on urban streamwater chemistry variability. Biogeochemistry 2019, 144, 61–84. [Google Scholar] [CrossRef]
  79. Alamdari, N.; Sample, D.J. A multiobjective simulation-optimization tool for assisting in urban watershed restoration planning. J. Clean. Prod. 2019, 213, 251–261. [Google Scholar] [CrossRef]
  80. Rajda, W.; Natkaniec, J. Jakość wody odpływającej z mikrozlewni podgórskiej o zróżnicowanym użytkowaniu. Quality of water drain from differently managed submontane watershed. Zesz. Nauk. Akad. Rol. Krakowie. Inżynieria Sr. 2001, 21, 33–40. [Google Scholar]
  81. Ostrowski, K.; Bogdal, A.; Natkaniec, J. Wybrane cechy fizykochemiczne wód opadowych. Selected physicochemical properties of rainwaters. Zesz. Nauk. Akad. Rol. Krakowie. Inżynieria Sr. 2000, 20, 15–27. [Google Scholar]
  82. Fernandes, A.C.P.; Fernandes, L.F.S.; Cortes, R.M.V.; Pacheco, F.A.L. The Role of Landscape Configuration, Season, and Distance from Contaminant Sources on the Degradation of Stream Water Quality in Urban Catchments. Water 2019, 11, 2025. [Google Scholar] [CrossRef] [Green Version]
  83. Baek, S.S.; Ligaray, M.; Pyo, J.; Park, J.P.; Kang, J.H.; Pachepsky, Y.; Chun, J.A.; Cho, K.H. A novel water quality module of the SWMM model for assessing low impact development (LID) in urban watersheds. J. Hydrol. 2020, 586, 124886. [Google Scholar] [CrossRef]
  84. Gao, W.; Duan, Z.Z.; Yan, C.A.; Liu, C.G. Influence of nutrient mitigation measures on the fractional export of watershed inputs in an urban watershed. Environ. Sci. Pollut. Res. 2020, 27, 18521–18529. [Google Scholar] [CrossRef]
  85. Beal, L.; Senison, J.; Banner, J.; Musgrove, M.; Yazbek, L.; Bendik, N.; Herrington, C.; Reyes, D. Stream and Spring Water Evolution in a Rapidly Urbanizing Watershed, Austin, TX. Water Resour. Res. 2020, 56, e2019WR025623. [Google Scholar] [CrossRef]
  86. Bradshaw, J.L.; Osorio, M.; Schmitt, T.G.; Luthy, R.G. System Modeling, Optimization, and Analysis of Recycled Water and Dynamic Storm Water Deliveries to Spreading Basins for Urban Groundwater Recharge. Water Resour. Res. 2019, 55, 2446–2463. [Google Scholar] [CrossRef]
  87. Abdi, R.; Endreny, T.; Nowak, D. A model to integrate urban river thermal cooling in river restoration. J. Environ. Manag. 2020, 258, 110023. [Google Scholar] [CrossRef]
  88. Forrest, N.; Stein, Z.; Wiek, A. Transferability and scalability of sustainable urban water solutions—A case study from the Colorado River Basin. Resour. Conserv. Recycl. 2020, 157, 104790. [Google Scholar] [CrossRef]
  89. Jarosińska, E.; Gołda, K. Increasing natural retention—Remedy for current climate change in urban area. Urban Clim. 2020, 34, 100695. [Google Scholar] [CrossRef]
  90. Fashae, O.A.; Ayorinde, H.A.; Olusola, A.O.; Obateru, R.O. Landuse and surface water quality in an emerging urban city. Appl. Water Sci. 2019, 9, 25. [Google Scholar] [CrossRef] [Green Version]
  91. He, J.; Wang, S.; Liu, H.; Nguyen, V.; Han, W. The critical curve for shallow saturated zone in soil slope under rainfall and its prediction for landslide characteristics. Bull. Eng. Geol. Environ. 2021, 80, 1927–1945. [Google Scholar] [CrossRef]
  92. Frazar, S.; Gold, A.J.; Addy, K.; Moatar, F.; Birgand, F.; Schroth, A.W.; Kellogg, D.Q.; Pradhanang, S.M. Contrasting behavior of nitrate and phosphate flux from high flow events on small agricultural and urban watersheds. Biogeochemistry 2019, 145, 141–160. [Google Scholar] [CrossRef]
  93. De Jesus-Crespo, R.; Ramirez, A. The use of a Stream Visual Assessment Protocol to determine ecosystem integrity in an urban watershed in Puerto Rico. Phys. Chem. Earth 2011, 36, 560–566. [Google Scholar] [CrossRef]
  94. Huang, Y.J.; Tian, Z.; Ke, Q.; Liu, J.G.; Irannezhad, M.; Fan, D.L.; Hou, M.F.; Sun, L.X. Nature-based solutions for urban pluvial flood risk management. WIREs Water 2020, 7, e1421. [Google Scholar] [CrossRef]
  95. Cipoletti, N.; Jorgenson, Z.G.; Banda, J.A.; Hummel, S.L.; Kohno, S.; Schoenfuss, H.L. Land Use Contributions to Adverse Biological Effects in a Complex Agricultural and Urban Watershed: A Case Study of the Maumee River. Environ. Toxicol. Chem. 2019, 38, 1035–1051. [Google Scholar] [CrossRef]
  96. Kang, J.H.; Park, M.H.; Ha, S.J.; Stenstrom, M.K. An empirical modeling approach to predicting pollutant loads and developing cost-effective stormwater treatment strategies for a large urban watershed. Sci. Total Environ. 2020, 760, 143388. [Google Scholar] [CrossRef]
  97. Chen, Z.; Wu, G.; Wu, Y.; Wu, Q.; Shi, Q.; Ngo, H.H.; Vargas Saucedo, O.A.; Hu, H.Y. Water Eco-Nexus Cycle System (WaterEcoNet) as a key solution for water shortage and water environment problems in urban areas. Water Cycle 2020, 1, 71–77. [Google Scholar] [CrossRef]
  98. Yu, W.Y.; Zhao, L.S.; Fang, Q.; Hou, R. Contributions of runoff from paved farm roads to soil erosion in karst uplands under simulated rainfall conditions. Catena 2021, 196, 104887. [Google Scholar] [CrossRef]
  99. Jebamalar, A.; Sudharsanan, R.; Ravikumar, G.; Eslamian, S. Rainwater Harvesting Impact on Urban Groundwater. In Handbook of Water Harvesting and Conservation: Basic Concepts and Fundamentals; Eslamian, S., Eslamian, F., Eds.; John Wiley and Sons: Hoboken, NJ, USA, 2021; pp. 207–224. [Google Scholar]
  100. Choi, J.; Kim, K.; Sim, I.; Lee, O.; Kim, S. Cost-Effectiveness Analysis of Low-Impact Development Facilities to Improve Hydrologic Cycle and Water Quality in Urban Watershed. J. Korean Soc. Water Environ. 2020, 36, 206–219. [Google Scholar] [CrossRef]
  101. Hosseiny, H.; Crimmins, M.; Smith, V.B.; Kremer, P. A Generalized Automated Framework for Urban Runoff Modeling and Its Application at a Citywide Landscape. Water 2020, 12, 357. [Google Scholar] [CrossRef] [Green Version]
  102. Grimley, L.E.; Quintero, F.; Krajewski, W.F. Streamflow Predictions in a Small Urban-Rural Watershed: The Effects of Radar Rainfall Resolution and Urban Rainfall-Runoff Dynamics. Atmosphere 2020, 11, 774. [Google Scholar] [CrossRef]
  103. Frot, E.; van Wesemael, B.; Benet, A.S.; House, M.A. Water harvesting potential in function of hillslope characteristics: A case study from the Sierra de Gador (Almeria province, south-east Spain). J. Arid Environ. 2008, 72, 1213–1231. [Google Scholar] [CrossRef]
  104. Frot, E.; Van Wesemael, B. Predicting runoff from semi-arid hillslopes as source areas for water harvesting in the Sierra de Gador, southeast Spain. Catena 2009, 79, 83–92. [Google Scholar] [CrossRef]
  105. Jimenez, M.; Perez-Belmont, P.; Schewenius, M.; Lerner, A.M.; Mazari-Hiriart, M. Assessing the historical adaptive cycles of an urban social-ecological system and its potential future resilience: The case of Xochimilco, Mexico City. Reg. Environ. Chang. 2020, 20, 7. [Google Scholar] [CrossRef]
  106. Kozak, D.; Henderson, H.; Mazarro, A.D.; Rotbart, D.; Aradas, R. Blue-Green Infrastructure (BGI) in Dense Urban Watersheds. The Case of the Medrano Stream Basin (MSB) in Buenos Aires. Sustainability 2020, 12, 2163. [Google Scholar] [CrossRef] [Green Version]
  107. Song, S.; Albert, C.; Prominski, M. Exploring integrated design guidelines for urban wetland parks in China. Urban For. Urban Green. 2020, 53, 26712. [Google Scholar] [CrossRef]
Figure 1. Location of the research area and boundaries of the examined catchments ((A)—urban; (B)—suburban; (C)—rural).
Figure 1. Location of the research area and boundaries of the examined catchments ((A)—urban; (B)—suburban; (C)—rural).
Water 14 03259 g001
Figure 2. Illustration representing hydrological processes examined in urban and rural catchments.
Figure 2. Illustration representing hydrological processes examined in urban and rural catchments.
Water 14 03259 g002
Figure 3. Curve number for urban area.
Figure 3. Curve number for urban area.
Water 14 03259 g003
Figure 4. Curve number for suburban area.
Figure 4. Curve number for suburban area.
Water 14 03259 g004
Figure 5. Curve number for rural area.
Figure 5. Curve number for rural area.
Water 14 03259 g005
Figure 6. Maximum potential retention for urban area.
Figure 6. Maximum potential retention for urban area.
Water 14 03259 g006
Figure 7. Maximum potential retention for suburban area.
Figure 7. Maximum potential retention for suburban area.
Water 14 03259 g007
Figure 8. Maximum potential retention for rural area.
Figure 8. Maximum potential retention for rural area.
Water 14 03259 g008
Figure 9. Retention potential of headwater catchment for urban area.
Figure 9. Retention potential of headwater catchment for urban area.
Water 14 03259 g009
Figure 10. Retention potential of headwater catchment for suburban area.
Figure 10. Retention potential of headwater catchment for suburban area.
Water 14 03259 g010
Figure 11. Retention potential of headwater catchment for rural area.
Figure 11. Retention potential of headwater catchment for rural area.
Water 14 03259 g011
Figure 12. Projected data for physicochemical properties of stream in studied urban area using PCA.
Figure 12. Projected data for physicochemical properties of stream in studied urban area using PCA.
Water 14 03259 g012
Figure 13. Projected data for physicochemical properties of stream in studied suburban area using PCA.
Figure 13. Projected data for physicochemical properties of stream in studied suburban area using PCA.
Water 14 03259 g013
Figure 14. Projected data for physicochemical properties of stream in studied rural area using PCA.
Figure 14. Projected data for physicochemical properties of stream in studied rural area using PCA.
Water 14 03259 g014
Figure 15. CCA plot with eigenvectors designated factors affecting physicochemical parameters in rainwater. Open squares show urban area. Red squares represent rural area.
Figure 15. CCA plot with eigenvectors designated factors affecting physicochemical parameters in rainwater. Open squares show urban area. Red squares represent rural area.
Water 14 03259 g015
Figure 16. CCA plot with eigenvectors designated factors affecting physicochemical parameters in river. Open squares show urban area. Red squares represent suburban area.
Figure 16. CCA plot with eigenvectors designated factors affecting physicochemical parameters in river. Open squares show urban area. Red squares represent suburban area.
Water 14 03259 g016
Table 1. Average concentrations of physical and chemical indices for surface water in the entire studied area.
Table 1. Average concentrations of physical and chemical indices for surface water in the entire studied area.
Parameter Urban AreaSuburban AreaRural Area
pH-7.908.108.17
ECµS∙cm−1811454236
DOmg∙dm−38.3810.4711.40
PO43−mg∙dm−30.840.030.04
NH4+mg∙dm−36.850.01<0.001
NO2mg∙dm−30.200.040.01
NO3mg∙dm−38.31.91.8
SO42−mg∙dm−370.440.917.9
Clmg∙dm−354.29.52.1
Ca2+mg∙dm−3107.065.044.0
Mg2+mg∙dm−320.818.25.42
Na+mg∙dm−328.39.85.2
K+mg∙dm−39.92.11.5
FeTotmg∙dm−30.430.550.13
Mn2+mg∙dm−30.140.090.02
TDSmg∙dm−3605324166
TSSmg∙dm−318.26.24.5
Table 2. Average levels of physical and chemical indices for rainwater in studied catchments.
Table 2. Average levels of physical and chemical indices for rainwater in studied catchments.
Parameter Urban AreaSuburban AreaRural Area
pH-5.55.04.9
ECµS∙cm−133.019.019.0
DOmg∙dm−3---
PO43−mg∙dm−30.430.190.15
NH4+mg∙dm−31.771.040.95
NO2mg∙dm−30.100.040.03
NO3mg∙dm−33.202.102.60
SO42−mg∙dm−33.305.002.60
Clmg∙dm−31.700.800.70
Ca2+mg∙dm−32.100.801.00
Mg2+mg∙dm−30.300.100.10
Na+mg∙dm−30.300.200.20
K+mg∙dm−30.500.600.20
FeTotmg∙dm−30.400.200.10
Mn2+mg∙dm−30.050.030.02
TDSmg∙dm−326.021.021.0
TSSmg∙dm−32.903.103.40
Table 3. Results for spatial autoregressive models. A level of the studied element in rainwater was chosen as a dependent variable, while levels of elements in the surface watercourse were used as independent variables. Only significant results were presented.
Table 3. Results for spatial autoregressive models. A level of the studied element in rainwater was chosen as a dependent variable, while levels of elements in the surface watercourse were used as independent variables. Only significant results were presented.
VariableOLS
Coefficient
SAR
Coefficient
Standard
Coefficient
Standard
Error
tp
Urban area
Constant12.532−23.0822.4520.2410.9230.834
NO323.7927.8621.3130.2716.7220.013
SO42−4.55216.9720.8831.3412.5620.045
PO411.4325.2520.1231.8931.0030.039
Cl0.3420.7620.0210.6310.2320.032
Suburban area
Constant53.382−85.9810.831--0.643
NO31.7231.2521.3110.2736.721<0.001
Mg2+24.933−02.973−0.741--0.041
Ca2+−32.55322.8910.883--0.021
PO40.0730.0253<0.0015.1620.0050.009
Cl<0.001<0.001−0.042<0.001−0.240.079
Rural area
Constant12.94256.3420.1320.0421.0920.028
NO311.0522.7520.0920.2726.6420.002
Mg2+14.683−1.083−0.9720.0333.9730.007
Ca2+−2.5421.893−0.7820.3632.6730.043
Note: Result for predictor variables using OLS is R2 = 0.98; Result for predictor + space is R2 = 0.89.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Halecki, W.; Stachura, T.; Fudała, W. Capacity of River Valleys to Retain Nutrients from Surface Runoff in Urban and Rural Areas (Southern Poland). Water 2022, 14, 3259. https://doi.org/10.3390/w14203259

AMA Style

Halecki W, Stachura T, Fudała W. Capacity of River Valleys to Retain Nutrients from Surface Runoff in Urban and Rural Areas (Southern Poland). Water. 2022; 14(20):3259. https://doi.org/10.3390/w14203259

Chicago/Turabian Style

Halecki, Wiktor, Tomasz Stachura, and Wioletta Fudała. 2022. "Capacity of River Valleys to Retain Nutrients from Surface Runoff in Urban and Rural Areas (Southern Poland)" Water 14, no. 20: 3259. https://doi.org/10.3390/w14203259

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop