1. Introduction
The demand for water, both in quantitative and qualitative terms, is increasingly significant having in consideration the population growth and economic development. Therefore, it is essential to manage water resources and other natural resources to guarantee access to water under suitable usage conditions [
1]. Given this importance, the United Nations (UN) within the framework of Agenda 2030 prioritizes drinking water and sanitation as one of the 17 Sustainable Development Goals (SDGs) for a global sustainable development agenda by 2030. SDG 6 aims to ensure the availability and sustainable management of water and sanitation for all. Its main objectives involve improving water quality, integrated water resource management, universal access to basic sanitation, water use efficiency, and the protection and restoration of aquatic ecosystems such as forests, rivers, lakes, aquifers, and wetlands [
2]. The assessment of water quality, and mitigation of major sources of pollution, is one of the pillars adopted by Brazil to fulfill SDG 6. This assessment is one of the steps in water resources management and is carried out through monitoring physical, chemical, and microbiological parameters via measurements of dissolved substance concentrations, physical properties, and analysis of pathogenic microorganisms [
3].
The spatial distribution of land use is a key factor that affects water quality, where industrial, urban, and agricultural land uses generally pose relatively high pollution risks compared to other types of land use such as forests and wetlands [
4]. In addition to land use, water quality can also be influenced by quantity, as it directly depends on the amount of water available to dissolve, dilute, and transport substances, whether beneficial or harmful to living organisms [
5]. Hence, water quality may be perceived as the outcome of both natural processes and human activities within the watershed, with land use and cover therein acting as determining factors. Understanding its relationship with land use and cover can be useful for identifying potential risks to water quality and coordinating water resource management with other uses of environmental resources [
3].
The Miranda River Watershed (MRW) is one of the main hydrographic regions within the Pantanal Biome, located in the state of Mato Grosso do Sul, Brazil, and entirely falls within the Upper Paraguay Basin (UPB). It is formed by the Miranda River as its main course and the Aquidauana River, and is situated in the transition zone between the Cerrado and Pantanal biomes, as well as the transition from the Plateau to the Pantanal Plain. Natural wetland areas are environments with high biodiversity associated with permanently or seasonally flooded areas. They are responsible for groundwater recharge due to their slow flow and long water retention time above ground, serving as "buffer zones" upstream and downstream of their territory, slowing down water flow and nourishing other ecosystems during dry periods. These areas retain large amounts of sediment and chemicals, improving water quality and the entire biological dynamics downstream. A considerable amount of pollutants that enter wetland areas are captured by vegetation and microorganisms or adhere to suspended particles. Consequently, sedimentation in wetlands contributes positively to water quality [
6].
Understanding the relationship between land use and water quality would help assess water quality in unmonitored basins, as monitoring is often costly and time-consuming. Additionally, this knowledge would provide guidelines for watershed managers and policymakers to prioritize future land use development [
7]. Thus, conducting an analysis of the relation between water quality and different types of land use and cover is an important step in the management of natural resources, developing conversion strategies, and guiding the sustainable use of these resources. In this regard, studies have analyzed the influence of land use and cover on water quality through multiple approaches; utilizing correlation coefficients between quality and land use [
8], comparative analysis between input and output [
9], through the use of regression equations [
10], and multivariate analysis [
11]. Although many approaches are being evaluated to analyze the influence of land use on water quality, most studies rely on limited time series data, often from only a single year or a single monitoring point, or few parameters. The MRW has a monitoring network with 28 stations distributed across the main rivers of the basin, as well as a quality and consistent historical data series available to the general public. However, few studies are reported with emphasis on this specific area. Overall, there are still few studies on the capacity of removal/retention of substances in natural wetland areas of tropical and subtropical regions, especially in South America [
12]. The fact is that the number of parameters and the monitoring frequency have been continuously expanded, generating water quality data from an important transition region between the Pantanal and Cerrado Biomes, as it is also located in the transition region from the Plateau to the Pantanal Plain. Therefore, analyzing the relationship between land use and water quality by jointly considering the temporal dimension with the spatial dimension is essential to understand the dynamics of water quality and land use in the region’s natural wetlands.
The objective of this study was to analyze the correlation between land use and land cover and the quality of surface waters in a natural wetland area in the MRW, in Mato Grosso do Sul, Brazil. For that it was identified which quality parameters were the most relevant (key factors) and the dynamics of substance concentrations along the main watercourses of the basin. Additionally, it aimed to analyze the influence of natural wetland areas on water quality, considering how these areas are capable of improving water quality.
4. Conclusions
The study demonstrated changes in land use and land cover patterns and water quality standards in the MRW from 2005 to 2018. It is possible to observe changes in land cover such as the reduction of native vegetation, as well as a decrease in livestock and an increase in agricultural plantations. Additionally, there was variability in the quantity of wetland areas in the MRW, characterizing multi-year cycles of drought and flood. Water quality parameters such as BOD, DO, pH, and Turbidity as well as the WQI, remained within the permissible reference limit for classification in Class 2, except for Thermotolerant Coliforms and Total Phosphorus, which showed values above the allowed range. The study also revealed a trend of increasing BOD values, indicating a higher input of organic matter into the water bodies of the basin, although in all the studied Influence Areas, the WQI was considered "GOOD" or "EXCELLENT". This indicates, overall, satisfactory levels of water quality for the basin. On the other hand, although Total Phosphorus exceeded the allowed limit, there was a gradual decreasing trend over the study period.
When analyzing correlations in the MRW, it was observed that when using a general approach for the entire basin, the indices showed low correlation. However, when correlations were analyzed by Influence Area (sub-basins), the correlation values increased, both for positive and negative correlations, indicating that policies for water quality management can be developed in terms of sub-basins, smaller planning units.
Comparing the behavior of water quality parameters concerning the wetland areas of the MRW, there is a significant improvement in WQI, and for Thermotolerant Coliforms and Turbidity. However, there is a downward trend in DO values in wetland areas and an increase in BOD. For pH, it is observed that values tend to increase, making the water more alkaline, rarely exceeding the range of 7.5. Regarding TP, TN, and TS, it was found that the values showed variability over the years, with periods of increase or decrease. Considering the climatic seasonality in the wetland areas of the MRW, it was observed that DO, Turbidity, and WQI values are higher in the driest times of the year, and COD, TN, and TP presented higher concentrations in the rainy season.
Wetland areas play important hydrological and ecological functions for flora and fauna, and especially for the water quality of the region. They are effective in the removal, storage, and cycling of nutrients, such as Nitrogen and Phosphorus, as well as improving and controlling quality parameters downstream. Therefore, they are highly relevant areas for defining strategies for watershed management and conservation, environmental protection, and the provision of ecosystem services.
Figure 1.
Miranda River Watershed with selected Water Quality Monitoring Stations (Sub-basins: 1-17) and the natural wetland area.
Figure 1.
Miranda River Watershed with selected Water Quality Monitoring Stations (Sub-basins: 1-17) and the natural wetland area.
Figure 2.
Monitoring stations in the natural wetland area of the Miranda river watershed, in the Pantanal Biome (IA: Influence Areas).
Figure 2.
Monitoring stations in the natural wetland area of the Miranda river watershed, in the Pantanal Biome (IA: Influence Areas).
Figure 3.
Land use and land cover data in the Miranda river watershed between the period of 2005 to 2018 (IA: Influence Area).
Figure 3.
Land use and land cover data in the Miranda river watershed between the period of 2005 to 2018 (IA: Influence Area).
Figure 4.
Percentage of compliance with water quality parameters according to Brazilian legislation in the Miranda Basin between the period from 2005 to 2018.
Figure 4.
Percentage of compliance with water quality parameters according to Brazilian legislation in the Miranda Basin between the period from 2005 to 2018.
Figure 5.
Variation of the Average of the Water Quality Index (WQI) in the Miranda River Watershed subdivided into 17 Influence Areas.
Figure 5.
Variation of the Average of the Water Quality Index (WQI) in the Miranda River Watershed subdivided into 17 Influence Areas.
Figure 6.
Correlogram of land use and land cover with water quality in the Miranda river watershed. Note: Class: Agri=Agriculture, Hidro=Water Body, Not_Veg=Non-Vegetated, Past=Pasture, Wet=Wetland Areas, Urb=Urban Area, Veg_Nat=Native Vegetation. Water quality parameters: DO=dissolved oxygen, BOD= biochemical oxygen demand, Therm. Coli.= Thermotolerant coliforms, TN= total nitrogen, TP=total phosphorus, TS=total solids, Tu=turbidity, WQI=water quality index.
Figure 6.
Correlogram of land use and land cover with water quality in the Miranda river watershed. Note: Class: Agri=Agriculture, Hidro=Water Body, Not_Veg=Non-Vegetated, Past=Pasture, Wet=Wetland Areas, Urb=Urban Area, Veg_Nat=Native Vegetation. Water quality parameters: DO=dissolved oxygen, BOD= biochemical oxygen demand, Therm. Coli.= Thermotolerant coliforms, TN= total nitrogen, TP=total phosphorus, TS=total solids, Tu=turbidity, WQI=water quality index.
Figure 7.
Correlogram of land with water quality by influence area in the Miranda river watershed. Note: IA: Influence Areas. Water quality parameters: DO=dissolved oxygen, BOD= biochemical oxygen demand, Therm. Coli.= Thermotolerant coliforms, TN= total nitrogen, TP=total phosphorus, TS=total solids, Tu=turbidity, WQI=water quality index.
Figure 7.
Correlogram of land with water quality by influence area in the Miranda river watershed. Note: IA: Influence Areas. Water quality parameters: DO=dissolved oxygen, BOD= biochemical oxygen demand, Therm. Coli.= Thermotolerant coliforms, TN= total nitrogen, TP=total phosphorus, TS=total solids, Tu=turbidity, WQI=water quality index.
Figure 8.
Correlogram of wetland areas with water quality by influence area in the Miranda river watershed. Note: Class: Agri=Agriculture, Hidro=Water Body, Not_Veg=Non-Vegetated, Past=Pasture, Umide=Wetland Areas, Urb=Urban Area, Veg_Nat=Native Vegetation. Water quality parameters: DO=dissolved oxygen, BOD= biochemical oxygen demand, Therm. Coli.= Thermotolerant coliforms, TN= total nitrogen, TP=total phosphorus, TS=total solids, Tu=turbidity, WQI=water quality index.
Figure 8.
Correlogram of wetland areas with water quality by influence area in the Miranda river watershed. Note: Class: Agri=Agriculture, Hidro=Water Body, Not_Veg=Non-Vegetated, Past=Pasture, Umide=Wetland Areas, Urb=Urban Area, Veg_Nat=Native Vegetation. Water quality parameters: DO=dissolved oxygen, BOD= biochemical oxygen demand, Therm. Coli.= Thermotolerant coliforms, TN= total nitrogen, TP=total phosphorus, TS=total solids, Tu=turbidity, WQI=water quality index.
Figure 9.
Spatial variation of the average Thermotolerant Coliforms (MPN/1000mL) at the monitoring stations in the Pantanal (IA07, IA14, IA15, IA16, IA17).
Figure 9.
Spatial variation of the average Thermotolerant Coliforms (MPN/1000mL) at the monitoring stations in the Pantanal (IA07, IA14, IA15, IA16, IA17).
Figure 10.
Spatial variation of the average Dissolved Oxygen (mg/L) at the monitoring stations in the Pantanal (IA07, IA14, IA15, IA16, IA17).
Figure 10.
Spatial variation of the average Dissolved Oxygen (mg/L) at the monitoring stations in the Pantanal (IA07, IA14, IA15, IA16, IA17).
Figure 11.
Spatial variation of the measurement of Biochemical Oxygen Demand BOD (mg/L) at the monitoring stations in the Pantanal (IA07, IA14, IA15, IA16, IA17).
Figure 11.
Spatial variation of the measurement of Biochemical Oxygen Demand BOD (mg/L) at the monitoring stations in the Pantanal (IA07, IA14, IA15, IA16, IA17).
Figure 12.
Spatial variation of the mean Total Nitrogen (mg/L) at the monitoring stations of the Pantanal (IA07, IA14, IA15, IA16, IA17).
Figure 12.
Spatial variation of the mean Total Nitrogen (mg/L) at the monitoring stations of the Pantanal (IA07, IA14, IA15, IA16, IA17).
Figure 13.
Spatial variation of the mean Total Phosphorus (mg/L) at the monitoring stations of the Pantanal (IA07, IA14, IA15, IA16, IA17).
Figure 13.
Spatial variation of the mean Total Phosphorus (mg/L) at the monitoring stations of the Pantanal (IA07, IA14, IA15, IA16, IA17).
Figure 14.
Spatial variation of Total Solids (mg/L) at the monitoring stations of the Pantanal (IA07, IA14, IA15, IA16, IA17).
Figure 14.
Spatial variation of Total Solids (mg/L) at the monitoring stations of the Pantanal (IA07, IA14, IA15, IA16, IA17).
Figure 15.
Spatial variation of the average Turbidity (NTU) at the monitoring stations of the Pantanal (IA07, IA14, IA15, IA16, IA17).
Figure 15.
Spatial variation of the average Turbidity (NTU) at the monitoring stations of the Pantanal (IA07, IA14, IA15, IA16, IA17).
Figure 16.
Spatial variation of the average pH at the monitoring stations of the Pantanal (IA07, IA14, IA15, IA16, IA17).
Figure 16.
Spatial variation of the average pH at the monitoring stations of the Pantanal (IA07, IA14, IA15, IA16, IA17).
Figure 17.
Spatial variation of the Water Quality Index (WQI) at the monitoring stations of the Pantanal (IA07, IA14, IA15, IA16, IA17).
Figure 17.
Spatial variation of the Water Quality Index (WQI) at the monitoring stations of the Pantanal (IA07, IA14, IA15, IA16, IA17).
Table 1.
Quality parameters and reference values for Class 2 classification Deliberation CECA No. 36/2012.
Table 1.
Quality parameters and reference values for Class 2 classification Deliberation CECA No. 36/2012.
Parameter |
Unit |
Analytical Method |
Class 2 Limit Value |
pH |
- |
4500–H [15] |
6-9 |
Dissolved Oxygen (DO) |
mg/L |
ASTM D888-12 [16] |
>5.0 |
Biochemical Oxygen Demand (5-day BOD Test) |
mg/L |
5210-D [15] |
≤5.0 |
Fecal coliforms |
MPN. 100 m/L |
9221–E [15] |
≤200 |
Total Nitrogen (TN) |
mg/L |
NBR 13796 [17] |
- |
Total Phosphorus (TP) |
mg/L |
4500 P–B, E [15] |
≤0.1 |
Total Solids (TS) |
mg/L |
2540 – C [15] |
- |
Turbidity (Tu) |
NTU |
2130-B [15] |
≤100 |
Water Quality Index (WQI) |
- |
Calculated |
- |
Table 2.
Existing water quality monitoring stations in the natural wetland, in the Pantanal Biome of the Miranda river watershed.
Table 2.
Existing water quality monitoring stations in the natural wetland, in the Pantanal Biome of the Miranda river watershed.
Position |
Station code [19] |
Name (Influence Areas) |
River |
Start |
00MS23AQ2284 |
IA07 |
Aquidauana |
Start |
00MS23MI1292 |
IA14 |
Miranda |
Intermediate |
00MS23AQ2000 |
IA15 |
Aquidauana |
Intermediate |
00MS23MI2147 |
IA16 |
Miranda |
Mouth |
00MS23MI2000 |
IA17 |
Mouth of Miranda |