1. Introduction
Changes in the natural environment are determined, with varying degrees of intensity, by the diverse human use of space. Each form of human activity has a specific impact on the quality of the environment. The most significant effects are related to the development of industry and the human settlement network. However, agricultural land use has a specific impact on the environment too. On the other hand, the intensity of impacts is determined by the resilience manifested in semi-natural areas, among other factors. Understanding the determinants of the impacts of diverse land use is a key element of environmental management [
1,
2,
3].
The quality of surface water varies in space and time, and human impacts can improve or worsen negative changes [
4,
5,
6]. Land cover and land use influence the individual environmental components, including their chemical properties. In the case of river water, this impact is multi-faceted and applies to the volume of flux as well as the physical and chemical properties of water [
7,
8,
9,
10,
11]. The influence of particular types of land use—woodland, agriculture, the built environment and industrial areas—on non-point source pollution of river water is varied and depends largely on their mutual proportions within a catchment [
12,
13]. What is highlighted is the particular role of urbanisation and suburbanisation processes leading to the expansion of the built environment [
14,
15]. In general, the greater the proportion of anthropogenic areas associated with human settlement and economic activity, the higher the level of river water contamination [
16,
17].
At present, considerable land cover changes can be observed in Europe, which can influence changes in the quality of river water in the future [
18,
19,
20]. A distinctive feature of these changes is the increased area of forests at the expense of areas used for agriculture as a result of land abandonment processes [
21,
22,
23,
24].
The study objective was to determine the impact of land use and land cover in a first-order catchment on the selected parameters of river water within a large area. While research on this scale is rarely conducted, it enables the development of a water resources management strategy at the province level; such a strategy can be based on controlling the type of land cover, for example. Furthermore, eastern Poland is a unique region in the European context because of the preserved mosaic of land use and land cover which is not encountered elsewhere [
25,
26].
2. Study Area
Lubelskie Province is located in eastern Poland (
Figure 1). It covers 25,123 km
2, which accounts for 8% of Poland’s territory. It is inhabited by 2.1 million people. The size of the population successively decreases as the depopulation processes and their consequences are increasingly visible, particularly in rural areas. The negative migration balance (general balance: −6495 people, including −5383 people of working age) is a very serious problem. A high negative natural increase rate occurs (−1.77‰; the national rate for Poland being −0.99‰). The population density, 84 people/km
2, is distinctly lower than the national average [
27]. About 46% of the inhabitants live in towns and cities (which total 48). Lubelskie Province is among the least developed regions not only in Poland but also in the entire EU. Its peripheral location is not conducive to economic development. This is directly manifested in the unemployment rate of 7.5% (compared to Poland’s national rate of 5.2%), which exceeds 15% in counties in the eastern part of the region. The value of gross domestic product per capita in Lubelskie Province is 32.2% lower than the national average, whereas the gross value added per one employed person in Lubelskie Province is 28.7% lower than the national average. The number of entities of national economy per 10,000 inhabitants (879) is lower than the national number (1175). The region’s economic development is also hampered by relatively poor transport accessibility, as evidenced by indices that are below Poland’s national average: public roads with hard paving, 91 km per 100 km
2 (Poland: 98.2 km per 100 km
2), among which national roads account for 4.6%, provincial roads for 10.2%, county roads for 40.7%, and district roads for 44.5%. The situation is even worse in the case of railway lines (3.9 versus the national average of 6.2 km per 100 km
2) [
28].
The land use structure is dominated by agricultural land that accounts for 70.5% of the province’s total area (60% is the national average). Forests and woodlots occupy 23.7% of the province’s territory, compared to the national average of 30.5%. The proportion of built-up and urbanised land (3.9%) is also lower than the average for all of Poland (5.5%). There is little underwater land (0.8% versus 2.1% for Poland) and wasteland (0.9% versus 1.5% for Poland). Capital expenditure on fixed assets for environmental protection amount to PLN 323 per inhabitant, which is the same as the national expenditure. Expenditure on water management is lower in Lubelskie Province (PLN 58) than in Poland as a whole (PLN 84). The share of capital expenditure on environmental protection in the total capital expenditure of the province is 5.3% [
28].
The diverse natural environmental conditions and the not very intensive economic use have enabled the fragmentary preservation of valuable natural assets. All legal forms of nature conservation occur in the Lublin region. Areas under legal protection account for a total of 22.69% (excluding Natura 2000 areas), national parks account for 0.72%, nature reserves for 0.47%, landscape parks for 9.56%, areas of protected landscape for 12.03%, while documentation sites, ecological lands, nature and landscape complexes account for 0.32% of the area of the province. Natura 2000 areas comprise the following: special bird protection areas (SPA) account for 13.4%, while special habitat protection areas (SAC) account for 6.6% [
27].
In national-level strategic documents, Lubelskie Province is recognized as an agricultural region and its development is programmed accordingly. Thanks to the occurrence of fertile soils (Luvisols, Cambisols), the province is a major centre of agricultural production in Poland. Lubelskie Province ranks second in Poland in cereal production (11.3%), second in rape and agrimony production (12.6%), third in sugar beet production (14.5%), and second in the production of fruits from trees (15.3%). Livestock breeding is of much lower significance as the province ranks seventh in Poland in cow milk production (5.7%) and ninth in slaughter livestock (4.3%). The fragmentation of agriculture is characteristic for most of the study area. Small family-owned farms are inefficient and incapable of commercial production. The average size of an agricultural holding is 9.05 ha (including 7.94 ha of farmland), while among the private (individually-owned) farms alone, it is 8.76 ha (including 7.68 ha of farmland) [
27].
Lubelskie has never been an industrial region. The political and economic changes in the late 1980s/early 1990s and the economic crisis led to the collapse of production and closing down or restructuring of many industrial plants. Small- and medium-sized towns were particularly affected by these problems. The collapse of the agri-food industry, including fruit and vegetable processing, led to changes in the structure of agricultural land and crops. The monitoring of rural area development in Poland [
29] shows that the indices in the individual components of socio-economic development—spatial, economic, social and those determining the quality of life—are low in Lubelskie Province in comparison with regions in western Poland in particular. Only in the case of districts in the vicinity of large cities (mainly Lublin), the level of development is higher.
The province has a varied natural environment. The northern part is a lowland composed of glacial sediments (sand, till and silt). Loess and limestone uplands occur in the central and SE part. The south-western part of the province is occupied by an old glacial plain (sub-Carpathian basins). The annual precipitation is about 550 mm, the mean annual air temperatures range from 7.0 to 7.6 °C, and the mean specific runoff rate is about 4.0 l/s·km
2 [
30].
Lubelskie Province lies in the basin of the Vistula river and is regarded as an area with small resources of surface water. The main rivers are: Wisła (Vistula), Bug, Wieprz, Krzna, Bystrzyca, and Huczwa. The river network is particularly sparse in the south-western part of the Lublin Upland and in Roztocze. The region’s water consumption needs are met by groundwater. Water consumption for the needs of the national economy and the population amounts to 130.8 m
3 per resident per year [
27]. The total water consumption for industrial needs is 99,195 dam
3, of which 18,827 dam
3 comes from underground water and 74,210 dam
3 from surface water. The following pollution loads in wastewater are discharged to the water or the soil: BOD
5 79,729 kg/year, COD 932,541 kg/year, total suspended sediment 241,024 kg/year, sum of chloride and sulphate ions 10,164,442 kg/year, total nitrogen 1,905,580 kg/year, total phosphorus 4477 kg/year. Total discharged wastewater amounts to 86,170 dam
3, including 5171 dam
3 of wastewater discharged to the sewage network and 80,999 dam
3 of wastewater discharged directly to the water or soil. Total treated wastewater amounts to 23,764 dam
3, 62,080 dam
3 are cooling waters (not requiring treatment), while untreated wastewater amounts to 326 dam
3 [
27].
Lubelskie Province has extremely varied conditions respecting the occurrence and circulation of groundwater. For most of the province’s territory—Lublin Upland, Roztocze, and southern Polesie (
Figure 1)—the main aquifer is in Cretaceous rocks and, locally, in Tertiary and Quaternary formations [
30]. Water in Upper Cretaceous carbonate rocks has the greatest use value. Some of it is highly mineralised. The groundwater resources are large despite their intensive economic exploitation. The region is rich in water springs (about 1600 inventoried springs, among which 1280 function continuously), most of which occur in its upland parts. The river valleys in the Lublin Upland and Roztocze dissect the main aquifer, and springs flow from the bases of the valley slopes. Highly efficient springs fed by infiltration water from the plateau tops also occur in the valley floors. The flow from the largest spring, in Zaporze, is 300 l/s. Springs play a significant role in river runoff and account for 19% of the runoff in the Vistula and Bug interfluvial zone, and for 38% in the western part of the Lublin Upland and Roztocze. The role of springs is particularly important during periods without rain, when springs supply more than 50% of the water and, in the case of very small catchments, nearly 100%. Lakes occur in the north-western part of the province—in the Łęczna-Włodawa Lake District. They are small (a few hectares), shallow, mesotrophic, and eutrophic [
30].
At present, the quality of surface water is considerably influenced by the following: excessive water uptake, discharges of urban and industrial wastewater, area source pollution (mainly from agricultural areas), and pollution related to the development of tourism and recreation [
31]. On the other hand, recent years have seen a decreasing trend in pollution introduced to water or soil. Biogenic substance content in treated wastewater decreases systematically. The predominant area source pollutants in Lubelskie Province are those generated by agricultural production [
28,
30].
The sewage system in the province is poorly developed, particularly in rural areas where a considerable part of the population uses septic tanks that do not always meet the applicable technical requirements. In general, wastewater treatment plants are used by 58.1% of the province’s inhabitants (94% in towns and cities, 27% in rural areas, according to [
28]). The difficulties with connecting to the collective sewage system result from the dispersed nature of settlements in rural areas, which hinders the complete provision of technical infrastructure (this also applies to water supply and mains gas). Another significant factor is the demographic situation in these areas: the outflow of young people and a high proportion of elderly people; social exclusion and poverty are not conducive to improving the quality of life. Lubelskie Province is a region in which there are serious deficiencies in household sanitary facilities, which, in extreme cases, leads to a peculiar paradox: the lack of access to running water from the mains means lower water consumption in a household, which means a smaller discharge of pollutants to the environment. Owing to high prices, particularly for collecting sewage, and low incomes, households connected to the water supply and sewage system are forced to be frugal. The number of pollutants discharged to the environment has definitely decreased due to the decreasing livestock population, also.
3. Materials and Methods
Data collected by the Chief Inspectorate of Environmental Protection as part of the national monitoring of surface water in Poland were used in the study. The forms and manner of conducting monitoring of surface water bodies is regulated by the 9 October 2019 Ordinance of the Minister of Maritime Economy and Inland Navigation [
32], and includes:
(1) The kinds of sampling points and criteria for their designation;
(2) The scope and frequency of assaying the individual aspects of the classification of ecological states and chemical conditions;
(3) The scope of assaying the individual aspects of the classification of ecological potential and chemical conditions;
(4) The reference methodology and conditions for ensuring the quality of measurements and assays;
(5) The number of sampling points, the criteria for their designation, and the scope and frequency of testing for substances particularly harmful to the water environment.
Approximately 34 catchments (surface water bodies) located in Lubelskie Province were encompassed by the study (E Poland) (
Figure 1). Overall, water quality monitoring covers 95 catchments. However, for most of them, only the most basic parameters are determined. Additionally, only first-order catchments were studied in order to obtain clear correlations between water quality indices and a clearly defined area of water supply. This method of catchment selection eliminates the influence of point sources of pollution located at greater distances from the sampling locations.
At the same time, the compared catchments have less diverse areas. Fifteen indices describing the state of water were considered: (i) biological parameters (2); (ii) physicochemical parameters: physical condition (2), oxygen conditions (5), salinity (3), acidity (1), and biogenic conditions (2). The methods of sampling and sample analysis are presented in
Table 1. Only indices that are available for a larger number of catchments were investigated. This allows us to recognize the relationships as representative.
The present study has used water quality data for 2018 from the Report prepared by the General Directorate for Environmental Protection “Assessment of the condition of water bodies of rivers and dam reservoirs in 2014–2019 based on monitoring” [
33]. The data presents averaged water quality results based on assays of samples collected six times a year at locations at the mouth of the selected river catchments.
Data related to land cover were obtained from the Topographic Objects Database (BDOT10k) [
34]. Six types of land cover and land use were distinguished: anthropogenic areas; orchards and plantations; arable land; meadows and pastures (grassland); water and wetland areas; and forests and woodlots. We also calculated the indicator showing the diversity of terrain functions (land use and land cover) proposed by Śleszyński (2013) [
35]. The share of specific LULC forms within catchments has been calculated using GIS software (ArcGIS).
Statistical analyses were then carried out to detect and identify correlations between surface water quality and land use/land cover in a catchment (correlation coefficients, clustering method, redundancy analysis). Pearson’s correlation coefficients were calculated for the percentage share of the individual types of land use/land cover and water quality parameters. The level of correlation significance was also determined. The spatial similarity of the studied catchments in terms of land cover was identified using Ward’s method [
36], i.e., agglomerative hierarchical clustering. In the calculation, Euclidean distances were used as measures of similarity between objects. A dendrogram showing a graphical representation of the similarities was thereby obtained. The calculations were carried out using the XLSTAT add-in in Microsoft Excel. Then, the mean values of water quality parameters were calculated for clusters of catchments with similar land cover characteristics. The correlation between land cover categories and river water quality indices was examined by means of redundancy analysis (RDA) [
37], which enables the investigation of correlations between two datasets. RDA was carried out using the XLSTAT add-in to Microsoft Excel. Owing to their different orders of magnitude, the variables were standardised for the purposes of the analysis. In the graphic image, type II scaling focusing on the correlations between the explanatory variables was used.
5. Discussion
The water in the catchments under study are of low quality in terms of biological indices, of quite a good quality with regard to oxygen indices, and of good quality in terms of salinity and biogenic substance content. The analyses indicate the occurrence of various correlations (including the lack of correlation) between the quality of river water and land use in catchments. It should also be noted that not all correlations were confirmed by more than one statistical method. The most important relationships found between LULC and river water quality are as follows:
- (1)
A larger number of macrophytes is fostered by the presence of meadows and pastures, water and wetland areas, as well as diversity of land use in catchments. This indicates that such forms of land use are conducive to the development of organisms that are indicative of better water quality.
- (2)
The mean water temperature is the lowest in catchments with varied land use.
- (3)
Total organic carbon content is greater in catchments with a high proportion of meadows and pastures, as well as water and wetland areas; it is lower in catchments used as arable land. This pattern results from the greater supply of organic matter to river channels in peat bog and grassland areas. Organic carbon is generally treated as a contaminant influencing the oxygen conditions of water.
- (4)
A similar situation occurs in the case of COD indices: they are higher in catchments with a high proportion of forests, meadows and pastures, as well as water and wetland areas. An increase in dissolved substance content is favoured by the use of catchments as arable land, while a clear decrease occurs in forest-covered catchments.
- (5)
The same pattern occurs in the case of chlorides and sulphates. The agricultural use of catchments is conducive to the supply of contaminants in dissolved form, while forests effectively mitigate this phenomenon.
- (6)
The presence of arable land, orchards, and plantations is conducive to increased nitrogen content in water. Nitrogen content is clearly lower in the water of catchments with a greater number of forests. The negative influence of agricultural land use and the positive impact of forests can also be observed in this case.
- (7)
No statistical correlations were found between types of land use and the following water quality properties: phytobenthos, total suspended sediment, total phosphorus, dissolved oxygen, BOD5.
Land cover is one of the factors influencing water quality, and the lack of unequivocal correlations found in some cases may result from the impact of other factors, such as topography, hydrogeological conditions, or the presence of point sources of pollution in a catchment [
10,
14,
39,
40,
41]. The poor hydrological connectivity between slope systems and river channels may have a particularly important influence. Water and the substances in it—originating from agricultural areas—do not always reach rivers. This phenomenon is particularly intensive in the case of a land use mosaic, which is linked to a varied microtopography conducive to local accumulation of contaminants [
41,
42]. What is quite significant in the case of small catchments in the Lublin Upland and Roztocze is the high proportion of water supply from the springs, which may exceed 50% in periods with less precipitation. This significantly reduces the impact of land cover on the quality of river water. It should be emphasized that the functional diversity index used in the study only partially reflects the land use mosaic (based on the main CLC classes) because it does not take into account the diversity resulting from the occurrence of a field mosaic.
Eastern Poland is extremely varied in terms of land use. Historical factors and changes in agriculture following changes in the political and economic system determine the present-day characteristics of the landscape. A mosaic of highly diverse forms of agricultural use has developed: meadows and pastures, orchards and plantations, and wasteland (
Table 7,
Figure 5). In some areas, individual fields are less than 10 m wide. The more diverse the land relief, the smaller the area covered by the individual forms of land use. Forests and woodlots add to this diversity. Field boundaries are barriers to the formation of surface runoff. Thus, the patchwork of different forms of land use—including an increased number of small areas of wasteland, woodland, woodlots, and shrubbery in a mosaic-like pattern—reduces the delivery of pollutants of agricultural origin to river water (
Figure 5). Therefore, this kind of land use should be recommended in policies concerning the sustainable development of rural areas both in Poland and other countries. Since a mosaic of land use reduces the threat posed not only by soil erosion but also the pollution of rivers, it should be preserved in the selected areas [
26].
Investigations conducted so far indicate that urban areas have an important influence on the deterioration of water quality [
17]. In the catchments under study, however, no correlation was found between the size of anthropogenic areas and the characteristics of surface water. It should be emphasised, though, that there are few typical urban areas in the region, hence the impact of this category of land cover is smaller than in the case of highly urbanised areas. The lack of clear correlations between LULC and water quality in the areas under study can generally result from the relatively weak anthropogenic pressure.
When searching for the determinants of water quality, the correlation between the type of catchment (type of sediments occurring in the substrate) and characteristics of river water was also analysed. For some indices, poorer oxygen indices were found for catchments where peat bog areas occur (
Table 8). This is related to the large supply of organic matter to river channels, which leads to the deterioration of water quality. On the other hand, the lowest indices (the best water quality) occur in rivers that are direct tributaries of the Vistula (
Table 8). At the same time, no statistically significant correlations were found between water quality parameters and catchment characteristics, such as area, density of the river network, or density of the road network (
Table 9).
The results of studies conducted in other areas are not unequivocal either. For example, Ding et al. (2015) [
11] found that in a subtropical river basin, the proportion of forests was positively associated with dissolved oxygen concentration but negatively associated with water temperature, total phosphorus, and total nitrogen. The proportion of urban land was strongly positively associated with total nitrogen. According to the study, agricultural land use did not have a significant impact on water quality. De Carli et al. (2018) [
44] found strong but not significant correlations between short-term agriculture and total phosphorus, long-term agriculture and turbidity, pasture and total nitrogen, pasture and total phosphorus forest and total phosphorus, and urban buildings and total nitrogen. Based on the analysis of publications concerning research into the impact of land cover on water quality, Camara et al. (2019) [
17] found that while a positive correlation was found between the area of arable land and quality of river waters, the correlation was not significant. Urban areas showed the most significant correlations with hydrological parameters (decline in water quality). A significant positive correlation was demonstrated between the size of forested areas and changes in the physicochemical parameters of water quality. Studies conducted in the USA indicated a significant correlation between land use and river water quality, particularly in the case of nitrogen and phosphorus [
9]. Research conducted in other areas of Poland indicates a positive impact of the share of meadows and pastures and a negative impact of arable lands on the quality of river water [
45,
46,
47]. All researchers emphasise that this influence is determined by several regional or local factors related to human activity and the characteristics of the natural environment.
Lee et al. (2009) [
10] point out that the structure of LULC patches, and not just the proportion of individual forms of land use, is an important influence on water parameters. During investigations conducted in Korea, they found that a greater deterioration of water in reservoirs occurs in areas where there is greater mosaicism of land use types. This particularly concerns the built environment: its negative impact on water quality is weaker when it is more compact. Large forest complexes, on the other hand, are conducive to an improvement of water quality. In the case of Lubelskie Province, forest areas are usually small, and this can be the reason for the lack of their positive impact on water quality. It should be stressed that the results obtained by [
10] contradict, to a certain extent, the previously suggested positive influence of land use mosaicism. However, one should take into consideration the differences in the structure of land use in Poland (agricultural country) and Korea (predominance of forest areas) and in the entirely different forms of agricultural land use.
Small farms and the associated mosaic of land use are not economically efficient. The process of land consolidation and intensification of agriculture may lead to their disappearance. Institutions responsible for land use policy should therefore pursue a rational policy of spatial management and creating a specific structure of LULC patches, and thus influence the quality of water. Spatial planning can be a tool for maintaining a mosaic of land use, and thus for water quality protection. The structure of land use and land cover in the Lublin region is the result of historical processes operating from the Middle Ages to the present day (including planned village layouts, farm management, parcelling, inheritance divisions, suburbanization, and deagrarisation). Nowadays, this mosaic of use shows high landscape values and, due to the stability of the ecological systems, it is conducive to maintaining the natural balance [
26]. At the same time, excessive fragmentation of the landscape may lead to the breakdown of the ecological network, and small, isolated units may not always fulfill their functions [
48]. In the case of many rural areas, there is no comprehensive shaping of their functional and spatial structure. The basis of a correct spatial policy is local spatial development planning, with different plans treating environmental issues in different ways [
49,
50]. In many territorial units, investment processes are carried out on the basis of decisions about land development conditions which do not always fully take into account the entire spectrum of natural conditions [
51,
52,
53].
6. Conclusions
A negative influence of land use on water quality (i.e., higher numbers of pollutants) has not been detected in the area under study. This influence is visible only in the case of some indices, and the statistical correlations found are usually not strong. No statistical correlations with the type of land use/land cover were found in the case of some indices.
In agricultural catchments, a greater number of substances dissolved in water were found (including sulphates and chlorides) and a greater amount of total nitrogen. Orchards and plantations and anthropogenic areas (the built environment) did not result in the decline of the indices analysed. Natural catchments (with forests, meadows and pastures, water and wetland areas) had better biological parameters (macrophytes) and smaller quantities of dissolved substances (forests). At the same time, water in catchments with a high proportion of these types of land use is characterised by poorer oxygen conditions (greater content of organic substances in water).
The quality of river water is a resultant of several natural factors and factors related to human activity. The lack of correlations with land cover may result from the diversity of the environmental characteristics of catchments, and this diversity (e.g., point sources of pollution) was not taken into account in the present study. A synergistic effect of multiple factors on water quality parameters is also possible.
One of the factors that can mask the impact of land use and land cover on river water quality is its unique structure. The mosaic of fields with different types of use results in a peculiar system of water circulation, i.e., increased micro-retention and infiltration within slopes and dry valleys. Therefore, the direct supply of water and substances suspended and dissolved in water is more limited.
The assessment of the impact of the land use mosaic on the physical and chemical parameters of river water requires further detailed investigation. The results of this study may suggest its positive influence in reducing the adverse effects of human activity in the environment. This problem requires further study. However, it seems that land cover mosaics can be regarded as a model form of land use. Spatial planning can be a tool for maintaining the mosaic of land use and, as a result, fertile soils and water quality.