Assessment of the Impact of Forestry and Settlement-Forest Use of the Catchments on the Parameters of Surface Water Quality: Case Studies for Chechło Reservoir Catchment, Southern Poland

The aim of the study was to determine the impact of natural and anthropogenic factors on the values of 22 quality indicators of surface waters flowing out of two small catchments differing in physiographic parameters and land use, in particular forest cover and urbanization of the area. The research was carried out in the years 2012–2014 at four measurement-control points located on the Chechło river and the Młoszówka stream (Poland), which are the main tributaries of the retention reservoir. Basic descriptive statistics, statistical tests, as well as cluster analysis and factor analysis were used to interpret the research results. The water that outflowed from the forestry-settlement catchment of the Młoszówka stream contained higher concentrations of total phosphorus, phosphates, nitrite, and nitrate nitrogen and salinity indicators than outflow from the Chechło river. Water from the Młoszówka stream was characterized by more favourable oxygen conditions. Higher oxygen concentration in the catchment influenced a large slope of the watercourse and thus higher water velocity, which is promoted by the mixed process. In the case of the forest catchment of the Chechło river, the water quality was generally better than in the Młoszówka stream, mainly in cases of total suspended solids TSS, total phosphorus TP, phosphates PO43−, total nitrogen TN, nitrite N–NO2−, nitrate N–NO3−, and salinity parameters. Despite it being a short section of the river taken into the study, favourable self-purification processes like mixed, nitrification, and denitrification were observed in its water. The research shows that forest areas have a positive effect on the balance of most substances dissolved in water, and natural factors in many cases shape the quality and utility values of surface waters on an equal footing with anthropogenic factors. In the case of a large number of examined parameters and complex processes occurring in water, the interpretation of the results makes it much easier by applying multivariate statistical methods.


Introduction
Anthropogenic activities most often have a negative impact on the environment, including the aquatic environment. As a result of these activities, changes in the relief and soil cover take place, land is drained or flooded, soil, forests, natural ecosystems, and surface and groundwater are degraded. All these changes, resulting from direct or indirect human influence on the environment, are called anthropogenic changes, and the actions causing them are called anthropogenic factors [1,2].

Descprition of the Study Catchments
The study included surface waters outflows out of two catchments-the Młoszówka stream and the Chechło river, both of which discharge to Chechło retention reservoir ( Figure 1). The reservoir was established in 1944-1945 mainly for water retention for industrial plants and fire-fighting purposes. In addition, it is used to equalize outflow in the Chechło river bed, reduces the effects of flooding in the valley below impounding level, and serves as a recreational area. With a normal level of impounding, the floodplain has an area of 38.50 ha and the amount of retained water reaches 587.5 thousand m 3 .
The catchment area of the Chechło reservoir is located (Latitude 50 • 8.29 N, Longitude 19 • 30.24 E) in the north-western part of the Małopolskie Voivodeship (Poland)-in the districts of Chrzanów and Krakow ( Figure 1). Several important communication routes run parallel through the catchment area, such as the A4 motorway, national road No. 79 on the Krakow-Katowice route, and the railway line connecting the Krakow and Silesian agglomerations. According to the Kondracki's [30] division of Poland into physio-geographical regions, the catchment area of the Chechło reservoir is located in province 34 (Polish Upland), in the south-eastern part of the Silesian-Krakow Upland-subprovince 341. In geological terms, the catchment of the Chechło reservoir belongs to the Silesian-Krakow monocline, which is a large plate gently inclined towards the northeast. In its geological structure there are Tertiary, Mesozoic, and Paleozoic formations. The southern part of the plate forming the ground is intersected with leaps running from NW to SE, between which the lowered rock strata form the Krzeszowicki rift fault and Dulowska Basin.
The Chechło river originates from ponds, wetlands, and peat bogs of the Dulowska Forest. It is a right-bank tributary of the Vistula River. About halfway along the Chechło river, i.e., at the 16 + 000 km, there is a front dam of the Chechło water reservoir. In abiotic terms, it is a type 6 watercourse, i.e., it belongs to the group of upland carbonate streams with fine-grained substrate on loess and loess-like soils [31]. The route of the upper running of Chechło river, from its sources to the reservoir, is not very winding and runs from the south-east to the north-west ( Figure 1). From the reservoir to the intersection with the highway, the few-meters wide riverbed of the Chechło river was strongly transformed by man by straightening its route and securing the base of the slopes with a fascine fence and as a result of the construction of a damming structure which enables the supply of fish ponds with water. Clearly more watercourses drain the right-hand side of the catchment ( Figure 1). The average levelled slope of the Chechło river is 0.5% and the density of the river network is 3.11 km·km −2 ( Table 1). From the hydrological point of view and in terms of protection of quantitative and qualitative resources of the Chechło river, attention should be paid to the presence of marshy areas and small water reservoirs in its catchment area. The Młoszówka stream currently discharges to the water reservoir, but originally it was one of the largest right-bank tributaries of the Chechło river. The water in it flows in a bed that shows little winding; in the upper part it has been paved or protected with prefabricated openwork slabs and concrete walls, while in the lower part it is turfed and locally protected with fascine fences. In the catchment of the Młoszówka stream, the river network is much less developed than in the case of the Chechło river, as its density is only 1.97 km·km −2 . Clearly more watercourses are located in the lower part of the catchment, while in the upper part there is only the Młoszówka stream with a rather long dry valley. The average levelled slope of the bottom of Młoszówka is 2.5 times greater than the Chechło river and amounts to 1.2% (Table 1, Figure 1).   The Chechło reservoir catchment at the site of the dam is 42.46 km 2 , of which 76% of the area belongs to the catchment of the Chechło river, 20% is the basin of the Młoszówka stream, and the remaining 4% of the area is the littoral (Figure 1). The Chechło river catchment covers an area of 32.32 km 2 and has an average and maximum length of 8.67 and 9.04 km. The shape indicator values allow considering the catchment as firm, medium-elongated, and strongly close to the wheel. The catchment area of the Młoszówka stream equals 8.35 km 2 and is almost four times smaller than the catchment area of the Chechło river. For these reasons, other geometrical parameters are also much smaller (Table 1). The value of the elongation index (0.41) indicates the lanceolate shape of the catchment, which may theoretically influence the extension of the runoff path and thus the smoother formation of freshets.
The catchments of both of the analyzed watercourses are upland in character, with the area of the Chechło river catchment being less sculpted, and the northern part of the Młoszówka stream catchment being the most sculpted ( Figure 1). The highest point in the Chechło river catchment is located on the watershed in the northern part of the catchment-Kowalikowa Mountain of height 418.4 m above sea level. The height difference is 143.4 m above sea level, the weighted average altitude is 300.9 m above sea level and the weighted average slope is 3.5%. The highest and at the same time the northernmost point on the watershed of the Młoszówka stream catchment is a hill without a name (419.2 m above sea level), situated at the foot of Ostra Mountain. The height difference of this catchment is similar, while the average weighted height is 46.5 m above sea level higher, and the weighted average slope is more than twice as high, compared to the Chechło river catchment (Table 1, Figure 1).
The soils of the Chechło reservoir catchment area are characterized by great diversity, both typologically and in terms of species. There are mineral soils of temperate (brown, lessive, and limestone soils) and dry (e.g., podzolic soils) habitats; there are marshy soils as well as organic soils at different stages of moorshification. On more than half of the area of the Chechło river catchment, mainly under the Dulowska Forest (Figure 2), poorly drained clayey sands developed into podzolic soils. In the valleys of watercourses, there are constantly or periodically over-wetted soils. These are mineral formations with grain size of medium clays, from which gley soils and fen soils were formed, organic-mineral formations in the form of muck-mineral, mucky, and peat-muddy soils, as well as organic formations in the form of shallow low peats. These soils, irrespective of their origin, belong to the medium grassland complex. In the central-southern part of the catchment, there are also lessive soils, brown leached soils and brown specific soils, black degraded soils, grey soils, and limestone soils. Depending on the moisture content and terrain, these soils belong to the wheat faulty or rye very good and good arable land complexes or form weak and very weak grassland complexes. The catchment area of the Młoszówka stream is dominated by brown specific and brown leached soils, which were formed mainly from strong and light clayey sands or loess formations. In the northern and southern part of the catchment, under large and dense forest complexes (Figure 2), podzolic soils were formed from poorly drained clayey and light clayey sands. In the lower part of the Młoszówka stream catchment there are shallow low peats and muck-mineral and mucky soils. Organic formations form marshy forest habitats or middle class grassland.
Land use and development have a significant influence on the balance and hydrologic and biogeochemical cycling in the catchment area. In this respect, it can be stated that the Chechło river catchment is forest dominated, as forests and areas of afforestation occurring mostly in a concentrated form in the central and eastern part (Dulowska Forest) occupy almost 67% of the catchment area. The second largest form of land use is grassy ecosystem occupying about 15% of the catchment area. These are both natural and anthropogenic habitats, which in some places are periodically too dry or over-wetted. Arable land occupies only 8.5% of the catchment area, with the largest area thereof in the northern, south-western, and source south-eastern part of the catchment. The built-up areas occupy 6.3% of the catchment area and are located mainly in the vicinity of asphalt roads (Table 1, Figure 2).
The catchment of the Młoszówka stream is also dominated by forests, but the forest cover of over 40% is much smaller than that of the Chechło river catchment. The forest areas are partially scattered, with the largest share in the northern and southern part of the catchment area. Grass ecosystems and arable land covering 24.3 and 14.8% of the total area respectively, are situated mainly in the central and northern part of the catchment. The built-up areas, mostly in the middle part, cover 15.1% of the area and together with the areas under the paved roads form 17.6% of the sealed catchment area (Table 1, Figure 2). In terms of land use, the catchment area of the Młoszówka stream can be considered as a mixture of forestry and settlement areas.
Compared to the Młoszówka stream, the maximum flows likely in all variants in the Chechło river are about 3.8 times higher (Table 1). Moreover, on the basis of the calculated maximum flows with a probability of 1%, it was found that the corresponding unit outflows calculated for the Chechło river catchment and the Młoszówka stream are at a very similar level and amount to about 1.02 m 3 ·s −1 ·km −2 .

Sampling and Analytical Methods
Monthly water samplings were conducted in years 2012-2014 at four measurement and control points, for which the following coding was applied (Figures 1 and 2): P1-A point located on the Młoszówka stream, 150 m above the mouth of the Chechło reservoir (km 0 + 150); P2-A point located on the Chechło river, 10 m above the mouth of the reservoir (km 17 + 710); P3-A point located on the Chechło river, 820 m above the mouth of the Chechło reservoir and 200 m below a small post-mining reservoir (km 18 + 520); P4-A point located on the Chechło river, 1750 m above the mouth of the reservoir and 50 m below a bridge on the A4 highway (km 19 + 450).
During the sampling for analysis, the following parameters were measured in situ: pH with the CP-104 pH meter and water temperature, dissolved oxygen concentration (DO) and water saturation (DO sat ) using a CO-411 oxygen meter. Also, in situ measurements of velocity of flow were conducted using current meter type HeGa. The following parameters were determined in the laboratory: with the gravimetric method for total suspended solids (TSS) and dissolved substances (DS); ion concentration of calcium (Ca 2+ ), magnesium (Mg 2+ ), sodium (Na + ), potassium (K + ), and chlorophyll "a" (Chl a) were determined by absorption atomic spectrometry on UNICAM SOLAR 969 spectrometer; concentration of nitrogen in ammonium (N-NH 4 + ), nitrite (N-NO 2 − ), and nitrate form (N-NO 3 − ) as well as total phosphorus, phosphates (PO 4 3− ), and chlorides (Cl − ) were determined using colorimetric method using the FIAstar 5000 apparatus; concentration of sulphates (SO 4 2− ) was determined using the precipitation gravimetric method; the five-day biochemical oxygen demand (BOD 5 ) was determined using the Winkler titration method, and the chemical oxygen demand (COD Cr ) and organic nitrogen (Org-N) using the titration method. The total nitrogen (TN) concentration was calculated by summing up the results of N-NH 4 + , N-NO 2 − , N-NO 3 − , and Org-N concentrations. All analytical tests were carried out in three replicates using a certified equipment.

Basic Data Analysis
For all the analyzed water quality indicators, measured separately at each gauging location, basic descriptive statistics were determined: minimum and maximum values (ranges), arithmetic mean values, and standard deviations. The coefficient of variation (CV) values were calculated from the ratio of standard deviation and mean value. On the basis of the CV, the assessment of variation of data in the data sets was made, with four classes of variation used: small variation (0-20%), medium variation (21-40%), large variation (41-60%), and very large variation (>60%). In order to test whether the values of surface water quality parameters examined at individual locations differ significantly, a statistical test was carried out. Since most of the compared empirical data groups were found to possess no normal distributions (W Shapiro-Wilk test) or homogeneity of variance (Brown-Forsythe test), non-parametric Mann-Whitney U test designed for independent samples was used assuming a significance level α = 0.05. Using the Vollenweider and Organisation for Economic Co-operation and Development OECD classification [32], the trophic condition of the examined waters was assessed.
A detailed analysis of physiographic characteristics of the researched catchments were performed with the use of GIS geoinformatics software tools and probable flows were calculated using regional formulas.

Cluster Analysis
With a large number of physicochemical indicators sampled in many locations, an accurate interpretation of water quality assessment results is quite difficult. In such a case, multivariate statistics can be used. One of the most frequently used is the hierarchical method of cluster analysis, in which single clusters are gradually combined by adding up smaller clusters into larger ones. The formation of clusters is a process in which the vectors of parameters are combined into cluster groups so that in each cluster the vectors are as similar as possible, and-in relation to other groups-as different as possible [33]. Cluster analysis was performed on empirical data standardized according to the following formula [34,35]: where: x ij -value for the i-th object and j-th feature; x i -arithmetic mean relative to the i-th feature; s i -standard deviation relative to the i-th feature. A square Euclidean measure (metric) of the distance between the objects described by their respective features was used as a measure of similarity for calculations. This metric is often used to classify objects due to its convenient graphical interpretation and simple mathematical properties [36]. It is a geometric distance in a multivariate space, calculated on the basis of the formula: Euclidean distances were calculated on the basis of raw data and not standardized data. To estimate the distance between clusters, the Ward method [37] was used, which is based on the analysis of variance. This method seeks to minimize the sum of the squares of deviations of any two clusters that can be formed at each stage. This method is deemed as very effective, although it aims at creating small clusters. It is a commonly used technique for regionalization in hydrology, meteorology, and hydrochemistry [37,38].
Cluster analysis was performed for six variants and presented graphically in the form of dendrograms, for the following groups of indicators: 1.
As a result of the cluster analysis, a cluster hierarchy was obtained, which was presented in the form of dendrograms. The interpretation of the results obtained was based on diagrams of the trend of agglomeration of data, which, due to the limitation of volume, were not included in this study.

Factor Analysis
The purpose of factor analysis was, among others, to detect the structure and general regularities in relationships between the examined physicochemical parameters of water, to reduce the number of variables, and to describe and classify the examined variables in the separated spaces described by the new factors [35]. This analysis uses mathematical models in the form of a system of linear equations. In general, in factor analysis, variables are transformed in order to replace them by a smaller number of factors on which the tested variables are linearly dependent and can explain the correlations between them in the best possible way. Factor analysis attempts to explain the correlations between the observations in terms of the underlying factors, which are not directly observable [39]. The analysis was performed for the following water pollution parameters: water temperature, pH, dissolved oxygen, BOD 5 , COD Cr , PO 4 3− , TP, Org-N, N-NH 4 + , N-NO 2 − , N-NO 3 − , SO 4 2− , Cl -, Ca 2+ , Mg 2+ , Na + , and K + . Interdependent parameters, such as TSS, O sat , TN, and DS, were omitted. In the group of pollution parameters, physical, aerobic, biogenic, and mineral parameters influencing water salinity may be included. The question arises whether there are any connections and regularities between these parameters in the sampled water and whether it is possible to replace these parameters with new factors that will represent the quality of the sampled surface water.
In the first stage of the analysis, a correlation matrix between all the variables was performed, which allowed to draw preliminary conclusions as to the regularities occurring in the sample. Then, to assess the significance of the correlation matrix, Bartlett's sphericity test [40] was used to verify the null hypothesis that the correlation matrix is a unitary matrix (H 0 : R = I), all correlation coefficients of which are equal to zero. When assuming a null hypothesis, there is no point in performing factor analysis. However, the hypothesis can be tested using chi-square distribution statistics with p(p − 1)/2 degrees of freedom: where: p-number of variables; n-number of cases; λ i -i-th eigenvalue.
In the present case, the value of U statistics was −233.113, which corresponded to a test probability below 0.0000001, meaning the assumption of a null hypothesis was rejected. Thus, it was confirmed that the analyzed correlation matrix is not a unitary matrix. In the next step, the eigenvalues informing about the hierarchy of significance of the respective factors in explaining the overall variance were determined. Based on the recommendations given by Pekey et al. [40], only those factors whose eigenvalues exceeded 1 were selected for further analysis. For a selected number of new factors (variables), factor loadings, interpreted as correlation coefficients of the individual parameters, and their corresponding new factors (variables) were calculated. Also a communality expressed as part of the variance of the variable explained by common factors was calculated. The results of the calculations are presented in graphical form, taking into account the process of rotation of the coordinate system made using the Varimax method [41].

Results and Discussion
The catchment areas of the Chechło river and the Młoszówka stream are characterized by high retention capacities. This is mainly due to the relatively small land gradient, the dominant northern exposure of the slopes which caused lower evaporation losses and, above all, the dominated forest land cover of the area. Taking into account the last factor, it can be concluded that both catchments have high retention capacity, as they are dominated by forests with areas of afforestation with 66.7% and 40.2% of the respective catchment area, respectively. Such a large share of forest area and their mostly concentrated distribution in the Chechło river catchment (Dulowska Forest) will, by retaining a significant amount of water, influence the levelling of outflow [42,43] and the circulation of mineral and organic components. However, it should be taken into account that in the case of long-term rain-free periods and high air temperatures, such a high forest cover may reduce groundwater resources through an increase in evapotranspiration [44]. This, in turn, may cause the occurrence of deepened and long-lasting low groundwater levels which, by transforming over time into hydrological droughts, may threaten the proper functioning of the Chechło water reservoir. From the point of view of quantitative and qualitative resources, attention should be paid to a few percent of wetlands, which may mitigate the effects of climate change [45] and reduce the amount of nutrients and sediments washed out of the catchment [46]. On the other hand, arable land located mainly in the upper part of the Młoszówka stream catchment, with relatively large slopes of the land, may contribute to the formation of erosion processes and influence the increased leaching of substances from the soils.
Assessing the impact of anthropogenic activities on hydrological regime and habitat condition, the Chechło river catchment can be considered as not at risk, as the share of urbanized areas does not exceed 10% of the catchment area [47]. Thus, the hydrological regime and geomorphological and habitat conditions may be similar to natural catchments. The Młoszówka catchment, on the other hand, can be considered at risk (with 17.6% of the impervious areas), with clear signs of degradation caused by anthropogenic activities, such as larger and more frequent runoffs changing the geometry of watercourses, unstable banks of watercourses, and deterioration of habitat quality. The catchment is also much more likely to undergo rapid freshets and deepened low waters. Unfortunately, the presence of forest areas in it is not very beneficial from the point of view of the retention role, as these areas are scattered and their largest share occurs in the lower part of the catchment area. This situation may potentially result in the largest runoffs, which occur in the headwaters of the catchment, not being sufficiently slowed down by vegetation [48]. This, in turn, can lead to larger and more violent freshets in the lower part of the catchment. The distribution of grasslands in the central part of the catchment should slightly alleviate the unfavourable conditions for the formation of freshets. In the case of the Chechło river catchment, the greatest potential for the formation of significant surface runoffs is in the northern part of the catchment. The shape of the catchment also plays a significant role in the formation of the outflow. The value of the elongation index at the level of 0.71 indicates that the shape of the Chechło river catchment is similar to a circle, which may theoretically influence the more rapid formation of freshets. In the case of the catchment of the Młoszówka stream, on the other hand, the value of the elongation factor of 0.41 proves its lanceolate shape. Theoretically, this may have an effect on the extension of the water run-off route and thus on the alleviation of freshets. However, with this catchment, a somewhat faster reaction to rainfall compared to the catchment of the Chechło river is possible. This is due to the fact that the average slope of the Młoszówka catchment area is more than twice as high, and its area and forest cover are smaller. Undoubtedly, this is also influenced by the greater share of impervious areas in the catchment of the Młoszówka stream, which results in a faster surface runoff and a more rapid outflow of water from the catchment.
Understanding the sources of supply, concentration, and circulation of matter components is essential for the assessment of pollution loads that can lead to deterioration of water quality [49]. Therefore, a detailed analysis of changes in the values of physicochemical indicators of waters outflow of the two catchments differing in physiographic parameters, including land use, was undertaken. Table 2 presents basic descriptive statistics of the examined water quality parameters, and Table 3 presents the results of statistical test carried out with the use of Mann-Whitney U test.
The sampled water in the Chechło river was warmer than in the Młoszówka stream, as shown in Table 2. The largest amount of TSS flowed out of the catchment of the Młoszówka stream, and the smallest amount thereof in the waters of the Chechło river ( Table 2). The reason for this are likely to be higher slopes, lower forest cover, and higher share of soils more susceptible to erosion in the catchment of the Młoszówka stream. Similar observations were noted by Halecki et al. [18], based on research conducted in small Carpathian catchments.
Due to the enrichment of water with alkaline substances, the pH values increased systematically along the length of the examined section of Chechło river. Throughout the research period, the reaction of river waters varied from slightly acidic to slightly alkaline. A similar reaction was also found in the water of the Młoszówka stream; however, in this case, the pH values were significantly higher on α = 0.05 (Tables 2 and 3), at some dates indicating it as strongly alkaline. This state of affairs was definitely influenced by the humus substances contained in the water that are washed away from poor and acidic forest soils located in the bed of the Chechło river catchment. Similar relationships were described in the studies of Muwamba et al. [50] in the catchment WS79, at Santee Experimental Forest in coastal South Carolina in the US.
The most favourable oxygen conditions prevailed in the Młoszówka stream; DO and DO sat (Tables 2 and 3), whose average slope of the riverbed bottom is over twice as high as in the Chechło river, which results in a higher speed of flowing water and intensity of aeration. The biggest problem with the quality of the Chechło river waters are the organic substances contained in them, which were determined on the basis of the values of indicators indirectly informing about the level of organic pollution-BOD 5 and COD Cr . The basic source of organic substances is the Dulowska Forest with its organic soils, swamps, and a large amount of decomposing plant matter. It can, therefore, be concluded that in the case of the waters of the river investigated here we are dealing with a natural source of pollution by organic substances, which contain mainly humus substances, as indicated by the brown color of the water. Analyzing the values of the four oxygen parameters and organic pollution at the individual measuring-control points, it was observed that the higher the values of BOD 5 and COD Cr , the lower the concentrations and saturation of water with dissolved oxygen-and vice versa. Such relationships between these parameters are related to oxygen consumption in the processes of biochemical and chemical mineralization of organic substances contained in water, as confirmed by studies by Singh et al. [51] and Sotomayor et al. [52].
Forest vegetation can function as a natural filter against pollutants dumped in watersheds; hence, vegetation cover, diversity, and structure are key drivers of water quality [53,54]. Taking into account the above, the concentrations of all analyzed nutrients were generally at a medium level, while at some study dates only trace values of phosphates and nitrite nitrogen were recorded in the water of the main tributaries of the Chechło reservoir. The highest concentrations of phosphates, which usually occurred at a relatively low level either way, were found in the outflow from the catchment of the Młoszówka stream. Also in the waters of the same watercourse the highest concentrations of TP were recorded, with the difference, however, that their concentrations were quite visible, but did not differ significantly between the examined watercourses for α = 0.05 (Tables 2 and 3). The reason for the higher content of these two biogenic parameters of pollution in the water of the Młoszówka stream should be sought in the high urbanization of the central part of its catchment and the presence of arable land. Zhao et al. [55] and Wan et al. [56] state that the settlement and agricultural land use is positively correlated with the content of nutrients in waters, in particular with phosphorus compounds. Undoubtedly, the water quality of the Chechło river is influenced by the buffer role of woodland source areas of the catchment in minimizing phosphorus pollution [57].  (Tables 2 and 3).
The higher intensity of these favourable processes taking place in Młoszówka is related to the large amount of oxygen dissolved in water, which was also confirmed by the studies of Kowalik et al. [29], and to the reservoir located in the lower part of the catchment, which is heavily overgrown with wetland vegetation and acts as a wetland sewage treatment plant. The most saline water flowed out of the Młoszówka stream catchment, as the values of all salinity indicators (DS, SO 4 2− , Cl − , Ca 2+ , Mg 2+ , Na + and K + ) were statistically significantly for α = 0.05, higher than in the Chechło river (Tables 2 and 3). For these reasons, the water flowing through Młoszówka was mostly medium-hard, while the water in the Chechło river was soft. This fact mainly stems from land use and build-up of the area, i.e., greater urbanization of the catchment of the Młoszówka stream, much greater forest cover of the Chechło river catchment and quite significant differences in soil types occurring in both catchments. This is consistent with the studies described by Zhang et al. [58], in which it was found that the greatest impact on the deterioration of water quality in watercourses are in urban areas, followed by arable land, while grassland ecosystems and forest areas usually have a beneficial effect on water chemistry. Similar observations were presented in the works of Kosmowska et al. [59], De Mello et al. [60,61], Namugize et al. [62], and Wijesiri et al. [63]. For these reasons, we can speculate on greater anthropogenic activities in the Młoszówka catchment, the central part of which is heavily built-up and the concrete and paved bed of the stream runs in many places parallel to asphalt roads. In addition, above the main road between Trzebinia and Krakow there are several interceptors that drain rainwater directly running off asphalt surfaces into the Młoszówka stream, which undoubtedly results in higher salinity and, consequently, in a deterioration of water quality. Differentiation of the measured pH values during the study period was small at all water sampling locations, which is confirmed by the coefficient of variation values within the range of 4-5%, as shown in Table 2. Low or average random variation occurred for DO, DO sat , and TP, as well as for most salinity indicators. Large differences in empirical data were recorded for temperature, COD Cr and N-NH 4 + , and very large for the other water quality parameters. In general, the physicochemical parameters of water flowing out of the Młoszówka stream were characterized by lower random variation compared to the outflow from the Chechło river catchment (Table 2). On the basis of hydrological measurements mean velocity between the profiles P4 and P2 was equal to 0.39 m·s −1 and distance between profiles was equal 1740 m. Thus the approximated travel time is about 1.23 h. For the reason that water samples were collected from the Chechło river in a few-minute interval, in the order of point P4, P3, and P2, it can be concluded that the quality parameters were measured in the same water. The description of the test results presented in Table 2 are based on average values, but the relationships related to the decrease in the length of the Chechło river were also recorded in each field measurement dates. In a distance of 1740 m of the Chechło river the values of organic pollution indicators decreased: in the case of BOD 5 by 11%, and COD Cr by 33% of the initial value recorded in point P4. The concentration of TN and Org-N slightly increased on the tested section of the river, but this state was probably affected by biodegradation of organic compounds occurring in the self-purification process, as a result of which nitrogen-containing compounds are also formed.
There are only a few small tributaries on the examined part of the Chechło river, in which water flows periodically. Therefore, it can be concluded slightly diluting effect on the decline in the value of some pollution indicators (dilution is also an element of the water self-purification process). On the other hand, the rapid flow of water in the river causing its oxygenation increases and oxygen is used in the process of biodegradation of organic substances, and thus affects the water self-purification.
According to the trophicity criteria proposed by Vollenweider [32], average TP concentrations at all measurement gauging locations allowed the water to be classified as hypertrophic, and-due to mineral nitrogen content-as mesotrophic waters. The assessment based on TP concentrations carried out in accordance with OECD criteria [32] confirmed that the sampled waters are hypertrophic. Chlorophyll "a" concentrations indicate mesotrophy of waters flowing out of the catchments of Chechło river and the Młoszówka stream.
On the basis of the cluster analysis carried out with the Ward method for a total of 18 physicochemical parameters and separately for physical parameters and acidification, aerobic, phosphorus, nitrogen, and salinity parameters, six dendrograms were obtained showing similarity or differences in the quality of the water tested at individual measurement and gauging locations ( Figure 3).
Taking into account all the examined indicators, the similarity of water quality between the locations was less noticeable than in the other analyzed variants, because the whole data aggregation process ended at an intra-cluster distance of 89, and at a distance of 12 two clusters developed (Figure 3a): the first consisting of three points located on the Chechło river (P2-P4) and the second comprising only point P1 located at the mouth of the Młoszówka stream into the water reservoir. This proves that there are large differences in the quality of water flowing in both of the sampled watercourses. In order to analyze the data in detail, a cluster analysis was carried out separately for each group of the examined physicochemical indicators.
The arrangement of clusters in the dendrogram for physical parameters and acidification (temperature, TSS and pH) indicates a very high similarity in the water quality of the Chechło river in points 2 and 3. Smaller values of the three analyzed parameters of water flowing out of the Dulowska Forest resulted in the formation of a single cluster in point P4. One-element cluster also occurred in the case of point P1 located on the Młoszówka stream, which was influenced by the highest values of all analyzed physical parameters and acidification (Figure 3b).
In the case of oxygen parameters (O 2 , BOD 5 , COD Cr ) at the level 2 cut-off point, three clusters were formed, including two single clusters at points P1 and P4, and one double cluster consisting of points P2 and P3 (Figure 3c). This arrangement of clusters was influenced by the highest concentrations of dissolved oxygen and the lowest COD Cr values in the water of the Młoszówka stream. On the other hand, the formation of a separate cluster in point P4 is related to the self-purification processes in the Chechło river water, which led to the decomposition of organic substances and contributed to the reduction of BOD 5 and COD Cr values at the measurement and control points along the river-P3 and P2.  Two clusters (cut-off level 2), including one three-element cluster consisting of P2-P4 points located on the Chechło river and one-element cluster on the Młoszówka stream (Figure 3d)  ions shows that the water flowing out of the Młoszówka stream catchment is much more saline than the water flowing down the Chechło river (Figure 3f). Such a cluster system is characteristic for river catchments with significantly different anthropogenic pressure, which is the case with the two tributaries of the Chechło reservoir under examination. From the presented correlation matrix (Table 4) it is evident that there is a group of interrelated variables. It can be seen that the pH is strongly positively correlated with all salinity ions and PO 4 3 Based on the analysis of eigenvalues, four factors were identified which together explain almost 71% of the total variation of the analyzed water quality indicators. The highest variation of parameters observed in the analyzed sample is explained by the first factor −45% of variation. Another three factors explained over 32% of the variation. In subsequent analyses, the variables were further reduced and eventually three factors were adopted, for which the minimum value is 1.31 and which explain 71% of the total variation (Table 5).
On the basis of a detailed analysis of the results presented in Table 5, it can be concluded that concentrations of SO 4 2− , Cl − , Ca 2+ , Mg 2+ , Na + , and K + have high factor loadings with the first factor, water temperature and DO with the second factor, and BOD 5 and Org-N with the third factor. Three new variables can therefore be created, which will represent water salinity (factor 1), oxygen conditions (factor 2) and water pollution by organic substances (factor 3). Factors selected in this way determine to the greatest extent the quality of water flowing into the Chechło reservoir. The conducted factor analysis largely coincides with the results obtained from the cluster analysis in which it was shown that well oxygenated waters with high salinity (mineralization) flow through the Młoszówka stream, and the waters of the Chechło river contain significant amounts of organic substances. A graphical summary of factor loadings for the analyzed water quality indicators is presented in Figure 4. Figure 4a shows that the highest factor loadings for Factor 1 occur in the case of mineral compounds. For Factor 2 a large share was found in case of water temperature, DO N-NH 4 + , N-NO 2 − , and COD Cr .      All these parameters represent aerobic conditions or are associated with the conversion of pollutants under aerobic conditions (N-NH 4 + and N-NO 2 − ). The opposite position of temperature and DO indicates their inversely proportional relationship. In the case of factor 3 we can observe a large share of Org-N and BOD 5 factor loadings (Figure 4a,b). This can be explained by the fact that BOD 5 is a parameter specifying the amount of oxygen needed for the decomposition of readily decomposable organic compounds contained in water, which is closely related to the mineralization process. Since organic nitrogen compounds decompose relatively slowly, they generally do not cause an increase in BOD 5 in the initial period.

Conclusions
The interactions between different forms of land use and water quality are difficult to estimate due to the huge number of parameters and complexity of the processes occurring in the catchment and streamwater [64,65]. In addition to these factors, there are physiographic parameters of the catchment that can positively or negatively influence the properties of surface water. These hypotheses are confirmed by the results of our research. Water outflow from the forest and settlement mixed catchment, which is characterized by higher land slopes and a high share of mineral soils, contained higher concentrations of phosphorus, phosphates, nitrite, and nitrate nitrogen as well as salinity indicators compared to the nearby forest catchment. In spite of a larger share of urbanized areas and roads, the water was characterized by more favourable oxygen conditions, which was a result of a large gradient of the bottom of the investigated watercourse and a smaller amount of organic substances subject to the process of decomposition. In terms of land cover, generally better water quality was found in the forest catchment, which confirmed the assumed research hypothesis about the beneficial effect of this type of land use on surface water quality. It should be noted, however, that water outflow from this catchment is also polluted, but, in this case, with organic substances of natural origin, i.e., from organic soils and from the decomposition of plant matter occurring in the Dulowska Forest. Despite the short section of the river under study, beneficial self-purification processes were observed in the water. The research shows that forest areas have a positive effect on the balance of most substances dissolved in water, and natural factors, in many cases, shape the quality and utility values of surface waters on an equal footing with anthropogenic factors.
The results of this study also allowed to state that in case of a large number of examined physicochemical parameters of water, the interpretation of the results for their relationships and interactions is rather difficult, justifying the use of multivariate statistical methods. Cluster analysis allows, among others, for objective determination of similarities of the examined groups of water parameters between the individual measurement and control points, while factor analysis allows to reduce the number of variables to those actually having greatest impact on the quality of surface waters [66,67].