Next Article in Journal
Towards Establishing Empirical Rainfall Thresholds for Shallow Landslides in Guangzhou, Guangdong Province, China
Next Article in Special Issue
Short-Term Effect of Eucalyptus Leachates on Green Food Webs in Headwaters
Previous Article in Journal
A Tariff Model for Reclaimed Water in Industrial Sectors: An Opportunity from the Circular Economy
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Exploring Key Determinants of the Periphytic Diatom Community in a Southern Brazilian Micro-Watershed

1
Graduate Program in Agricultural Engineering (PGEAGRI), Western Paraná State University, Cascavel 85819-110, PR, Brazil
2
Botany Department, Federal University of Paraná, Curitiba 80060-000, PR, Brazil
3
Graduate Program in Ecology and Evolutionary Biology (EBIO), University of Colorado Boulder, Boulder, CO 80309, USA
4
Graduate Program in Agricultural Energy Engineering (PPGEA), Western Paraná State University, Cascavel 85819-110, PR, Brazil
5
Graduate Program in Conservation and Natural Resources Management (PPRN), Western Paraná State University, Cascavel 85819-110, PR, Brazil
*
Author to whom correspondence should be addressed.
Water 2022, 14(23), 3913; https://doi.org/10.3390/w14233913
Submission received: 22 October 2022 / Revised: 17 November 2022 / Accepted: 28 November 2022 / Published: 1 December 2022
(This article belongs to the Special Issue Freshwater and/or Brackish Diatoms: Ecology and Bioindication)

Abstract

:
Associating anthropogenic effects with variations in biodiversity is key to understanding how anthropogenic impacts are extrapolated in public supply micro-watersheds. The structure and dynamics of metacommunities in aquatic environments depend not only on the river network itself, but on a multitude of factors. Therefore, we associate the density and species richness of diatoms, assessed in a micro-watershed, with the following driver factors: local environmental variables, spatial variables, landscape characteristics, and the historical community, comparing their possible dispersal routes. Variance partitioning was performed using partial RDA models, with prior selection of predictor variables, to estimate the relative role of each predictor in the diatom community. The small scale of the micro-watershed resulted in a small spatial gradient, reflecting in the low variation in community richness across sampling stations. However, temporal heterogeneity associated with fluctuating precipitation throughout the year may cause temporal variation in the relative abundance of species. This pattern is a result of the supply of resources that increases biodiversity over time, as it allows the coexistence of species that alternate between dominance and persistence. Thus, even on a small scale and during one year of sampling, we demonstrated that predictors of different natures act together to explain diatom communities in micro-watersheds.

1. Introduction

Water resources have been considered the most important of all natural resources [1]. Despite being essential for human survival, whether for consumption or other applications, drinking water is poorly distributed across the planet [2,3]. In addition, territories with high water availability such as Brazil are facing a generalized water crisis, with an imminent reduction in the level of water resources [4]. In the Brazilian scenario, these consequences are maximized by population expansion, affecting the generation of electricity and, above all, the potability of water, especially in rivers used for public supply [5]. By providing environmental and ecosystem services central to human society, natural resources are indispensable in maintaining the communities present in lotic environments and all biodiversity [6]. However, anthropic activities are modifying natural characteristics and threatening ecosystems associated with water bodies [7]. Some important discussions in ecological studies aim to understand which factors explain species variation in aquatic communities [8,9]. Possible dispersal routes, environmental characteristics, and landscape configuration are often considered predictors of metacommunity distribution patterns [10,11].
Although it is agreed that the aforementioned factors act at different spatial scales [12], the relative importance of these factors is likely to change [13]. On a regional scale, the landscape tends to become homogenized as a consequence of the history of this territorial cutout [14]. At finer scales, local environmental factors (abiotic variables and biological interactions) are probably more important in explaining environmental filtering over smaller territorial extents [15]. In addition, the smaller the territorial extent, the greater the cumulative anthropic pressure on aquatic ecosystems due to the significant portion of land use that mixes agricultural and urban purposes, which contribute to reduced biodiversity [16]. High population densities imply the entry of domestic, industrial, and agricultural effluents into streams, loss of riparian vegetation, conversion of forested areas to agricultural use, and increases in impermeable areas in urban settlements [17]. The sum of these factors not only simplify the landscape, but also affect aquatic communities, altering limnological characteristics and impairing water quality [18]. With this in mind, linking anthropogenic effects to biodiversity variation is critical to understanding how anthropogenic impacts are extrapolated in public supply micro-watersheds (see [19]).
Among aquatic communities, periphytic diatoms are commonly used to assess the effects of environmental and spatial predictors [20], since they are algae that easily adhere to some substrate through mucilage production, which can restrict their dispersal and maximize the influence of these factors [21]. Furthermore, these microorganisms respond efficiently to changes in water characteristics, as well as in climate, geology, and land use in micro-watersheds [22,23]. Despite this, few studies use periphytic diatoms as a model of metacommunities in lotic environments [24].
The development of computational tools and increased data availability allows for better comparisons of the relative role of different predictors in structuring metacommunities. Nevertheless, studies comparing landscape descriptors, local descriptors, and spatial configuration for diatoms are rare [25,26,27]. Spatial configuration is commonly used as an indicator of dispersal routes [28]. However, dispersal of microalgae can occur via multiple vectors: by air [29], by flow-independent water (e.g., [30]), and by flow-dependent water (e.g., [31]). Considering each of these different dispersal hypotheses may more clearly reveal the relative role of predictors of metacommunities [32].
In Brazil, research involving dispersal routes of diatoms related to environmental and spatial filtering are usually developed in reservoirs addressing phytoplankton communities. Bortolini et al. [33] demonstrated that connectivity and dispersal by water flow in a reservoirs–river–floodplain gradient maximize dispersal rates of phytoplankton communities, contributing to increased biomass and homogenization of neighboring communities. Zorzal-Almeida et al. [34] observed that environmental and spatial factors together have a strong influence on variation between planktonic and surface sediment diatom communities in tropical reservoirs, indicating mass effects as a strong structuring predictor. Oliveira et al. [35] also evaluated the relative importance of local and spatial environmental processes on planktonic diatoms in streams in central Brazil, indicating that both act together in affecting metacommunity structuring. Medeiros et al. [36], however, evaluated the response of the periphytic diatom community through environmental predictors along a linear gradient using a distance-by-watercourse matrix. There is greater interest in the study of reservoirs due to their applicability; however, it is also crucial to understand the dynamics of river metacommunities belonging to micro-watersheds. Although these ecosystems are heavily impacted by human actions, they are still essential for public supply, and biomonitoring them provides information on water quality, consequently helping to maintain social health.
Indeed, it is agreed that the structure and dynamics of metacommunities in aquatic environments depend not only on the river network itself, but on a multitude of factors [37]. Understanding how these factors reflect on diatom community structure can help to understand ecological processes in river systems, which are important for the maintenance of the entire ecosystem, but suffer severe threats [38]. Our aim was to compare possible determinants of the periphytic diatom community in a micro-watershed: (i) determinants generated to represent different dispersal routes, (ii) determinants referring to local microhabitat filtering, and (iii) determinants related to landscape-scale environmental filtering. Thus, diatom community structure was explained both by considering the relative density of species (representing dominance and rarity relationships), and by the simple presence and absence of species (representing local occurrence capacity). Finally, we also considered the community composition from a previous period to explain the diatom community of the following period, thereby considering previous colonization as one of the determinants of the metacommunity (see [39]).

2. Materials and Methods

2.1. Study Area

Cascavel city (24°57′21″ S and 53°27′19″ W) is located in the western region of the State of Paraná, a subtropical area with a mean annual temperature of 21 °C [40]. The Cascavel River (24°32′ and 25°17′ S; 53°05′ and 53°50′ W) has a drainage area of 17.50 km2 divided in similar percentages between agricultural/urban areas and reserve or forest areas [41] (Figure 1). The main riverbed is 17.85 km long, with the main headwaters located in urban areas [42].

2.2. Sampling Design

Eight sampling stations were selected in the micro-watershed of the Cascavel River (Table 1), considering the land use and occupation for agricultural, industrial, extractive activities, and urban settlements, following the head–mouth direction. Samplings were performed quarterly (one per season) during 2018, totaling four samplings for each community, with communities 1, 2, 3, and 4 being associated with March, June, August, and December, respectively.

2.3. Landscape and Climate Variables

The land use and occupation measurements were performed by manual classification in QG is software (available in open licensing), using images obtained from the Sentinel 2A satellite with a spatial resolution of 10 m. The study area was classified according to five main classes: Forest, Mining, Water, Rural, and Urban. The divisions were based on the Land Use Technical Manual [43], which provides information about the types of land use included in each class (Figure 2). Climatological variables such as temperature, air humidity, wind speed and direction, solar radiation, and rainfall were provided by the Paraná Meteorological Institute (SIMEPAR).

2.4. Local Environmental Variables

The physical and chemical variables such as temperature, electrical conductivity, dissolved oxygen, pH, and turbidity were measured in situ by a HORIBA multiparameter probe, model U-5000.
We collected water samples for physical–chemical analysis by filling 5 L polyethylene gallons, and also 100 ml of water using dark glass bottles for microbiological analysis, both cooled and kept in the dark until their use in the laboratory. The oxygen consumption due to chemical oxidation (COD) and organic matter (BOD), concentrations of total Kjeidahl nitrogen (NT), total dissolved phosphorus (TDS), orthophosphate (PO4), nitrate (N-NO3), ammonium (N-NH4) and chlorophyll a were analyzed by the Environmental and Agronomic Analyses Laboratory (ACQUASOLLUS-Campo Mourão), following Standard Methods [44].
The superficial sediment was collected following the protocol of the Environmental Company of the State of São Paulo and the National Water Agency [45]. Samples with approximately 500 g of superficial sediment (0–5 cm depth) were collected utilizing plastic bags. Three sediment samples were collected at each sampling station and mixed to constitute a composite sample. In the laboratory, the samples were dried at room temperature in a closed place, protected from sunlight, wind, and excessive light. Then, the sediments were sieved using a 63 µm mesh to perform metal analyses. The sediment processing and analyses were performed in the laboratory of Unioeste—Cascavel campus. Aluminum (Al), barium (Ba), cadmium (Cd), calcium (Ca), chromium (Cr), cobalt (Co), copper (Cu), iron (Fe), lead (Pb), magnesium (Mg), manganese (Mn), nickel (Ni), potassium (K), sodium (Na), and zinc (Zn) were read, according to U.S. Environmental Protection Agency method USEPA 3050B [46], analyzing the silt/clay fraction of the soils using atomic absorption spectrometry (AAF) or inductively coupled plasma atomic emission spectrometry (ICP-AES).

2.5. Dispersal Route Indicators

Possible dispersal routes were generated considering the following methodologies:
(1)
To test “overland”, we considered the Euclidean distances resulting from the Principal Coordinate of Neighbor Matrices (PCNM, [47]), an ordination matrix generated using the coordinates of the sampling stations, creating a spatial connection network with links at all evaluated sites without directionality (i.e., the sampling stations were connected by two directions, upstream–downstream and downstream–upstream).
(2)
To test asymmetric eigenvector mapping “AEM”, the latitude and longitude of the sites were used to relate the points, considering whether the sites had connectivity or not [48].
(3)
To test “watercourse”, the distance matrix was generated by drawing a dendritic system over a real river map with Quantum GIS software. Thus, we modeled a spatial network with directional links over the watercourse (considering that the sites were connected only as a function of the river water flow).

2.6. Periphytic Community

Three substrates (stones), three meters apart from each other, were collected at all sampling stations, making a composite sample in order to evaluate the periphytic diatom community. The biofilm containing the diatoms was scraped from all the rocks using a toothbrush and distilled water. The density and relative abundance of taxa were estimated based on standard scraping of the stones (area equal to 9 cm²), considering the down-stream face of each substrate. Samples were preserved in 1% Lugol’s acetic solution (Bicudo and Menezes 2017) and stored at the Herbário da Universidade Estadual do Oeste do Paraná (UNOPA). Permanent slides were prepared with oxidized subsamples according to Simonsen’s [49] technique, modified by Moreira-Filho & Valente-Moreira [50], with known volume (1.0 mL) and Naphrax® resin (refractive index = 1.74) as the mounting medium. Quantitative analyses were performed in light microscopy (LM) at 1000× magnification using an Olympus B × 34 microscope equipped with a DP 71 capture and phase contrast camera. Individuals were identified and counted according to Kobayasi & Mayama [51], (1982), ensuring a minimum efficiency of 90% [52]. Diatom valves were counted in an optical microscope (1000×), as oxidation tended to separate frustules (Moro and Bicudo, 2002), including at least 50% of the full-size broken ones. Cells in pleural view were identified when possible. We followed the recommendations of Battarbee [53] and Lobo et al. [54] for counting under light microscopy, calculating densities per cm2. Abundant (valve density exceeds the average density for the sample) and dominant species (valve density is over 50% of the total density) were determined according to Lobo and Leighton (1986).

2.7. Data Analysis

Environmental variables that were not collinear and previously standardized by standard deviation were subjected to Principal Component Analysis (PCA) (see [55]) in order to characterize the sampling sites. Spatial (between sampling sites) and temporal (between periods) variation in species density and occurrence, as well as in environmental variables, were assessed using non-parametric permutational multivariate analysis (PERMANOVA, [56]). A Bray–Curtis matrix was used for species density after Hellinger transformation, a Sorensen matrix for occurrence, and a Euclidean distance matrix with environmental data previously standardized by standard deviation (see [55]).
The diatom community matrix was related to matrices with the following factors: local environmental variables, spatial variables (connectivity by dispersal routes), landscape characteristics, and the community from the previous period. The effect of previous community colonization was considered in the explanation by using matrices summarized in a principal coordinate analysis applied to the Bray–Curtis or Sorensen matrix (abundance or occurrence, Principal co-ordinates analysis (PCoA)), [57], with the diatom community from the previous period considered as a predictor of diatoms from the following period (except for the first collection period, see also [39]).
In PCoA analysis, data are divided into components to maximize the linear correlation between data points in a dissimilarity matrix. Through a “coordinate transform”, the observed values in the sample stations are replaced with newly derived coordinates. Thus, we reduce the dimensionality of a dataset by discarding coordinates that may not capture a threshold of variance in the community data. This technique preserves the global structure of the data while projecting it in low dimension. The predictor variables, i.e., environmental and landscape variables (standardized by standard deviation), were also previously scaled using a PCoA applied to the Euclidean distance matrices. The PCoA was used in order to improve the correlations of the predictor matrices with the community matrices [55,58].
Variance partitions were then performed using partial RDA models (pRDA, [55]), with prior selection of predictor variables (i.e., environmental, landscape, dispersal hypotheses, and previous community) through a stepwise approach [48,59], seeking to estimate the relative role of each predictor in the diatom community. Models were generated separately for each period, considering relative densities or presence and absence.
All analyses were performed using the R Development Core Team language [60] and its environment for computational statistics, along with the vegan [61], adespatial [62] and spdep [63] packages.

3. Results

The characterization of the sampling stations regarding the landscape variables pointed out the differences between these environments and the dominance of urban and agricultural areas in the Cascavel River micro-watershed. Sampling stations 2, 3, 4 and 7 are predominantly urban, stations 5 and 6 are predominantly used for agriculture and livestock, and stations 1 and 8 have a larger area of forest than the others (Figure 2). The urban influence area strongly influenced the sampling stations in the headwaters of the micro-watershed, although they have natural vegetation, which was also observed for points under the rural extension influence area (Figure 2).
Throughout 2018, the Cascavel River micro-watershed suffered great interference from seasonal variation, showing high accumulated precipitation in the first sampling period when compared to the others (Figure 3).
The environmental variables’ spatial variation (F = 2.07; R² = 0.38; p = 0.001) and temporal variation (F = 3.21; R² = 0.26; p = 0.001). The Principal Component Analysis (PCA) summarized 31.91% of the total variability of the sampled data in the first two axes (Figure 4). The scores showed the separation of the environments according to the sampling period and the land use at the sampling points. We highlight the interference of the macroscale in the March sampling period (the highest accumulated precipitation among the evaluated periods—557.2 mm in 30 days). This period was associated with higher values of Al, Cu, Nitrate and Chlorophyll a. In contrast, June had a lower association with metals in the surface sediment, as well as higher COD and total phosphorus values. Additionally, the analysis grouped stations 1, 4, and 7 (all first-order streams, and heavily impacted by urbanization) with higher values of Na, Zn, Cr, and Pb.
We found 118 infrageneric taxa, distributed in 37 genera. Eunotia Ehrenberg, Pinnularia Ehrenberg (11 spp. each), Gomphonema Ehrenberg and Navicula Bory (9 spp. each) showed high representative taxa richness among all sampling stations (see the means of relative abundance for each species according to sampling stations in Appendix A). The variation in species density for the four samplings showed spatial variation (F = 1.82; R² = 0.35; p = 0.001), but not temporal variation (F = 1.17; R² = 0.11; p = 0.168). In contrast, there was no spatial variation in species occurrence (F = 1.02; R² = 0.23; p = 0.4), but occurrences varied in time (F = 1.98; R² = 0.18; p = 0.001). However, some species differentiate the sampling stations: Achnanthidium modestiforme (Lange-Bertalot) Van de Vijver (S1, S2, S4, S6, and S8, with the greatest relative abundance at S4), Cymbopleura naviculiformis (Auerswald) Krammer (S3), Eunotia yberai Frenguelli (S4), Navicula cryptocephala Kützing (S1, S2, S3, S6, and S8, with the greatest relative abundance at S2 and S8), Navicula cryptotenella Lange-Bertalot (S1, S2, S3, S5, and S7, with the greatest relative abundance at S5), Navicula salinicola Hustedt (rare at S1, common at S2), and Spicaticribra kingstonii Johansen, Kociolek & Lowe (S3 and S8). We also observed the dominance of Eunotia botulitropica Wetzel & Costa (June, S7), Eunotia veneris (Kützing) De Toni (December, S1, with relative abundance higher than 50%,), and Gomphonema parvulum (Kützing) Kützing (December, S3). Species variation between sampling stations is described in the Appendix A. The explanatory power of the selected predictor variables was relatively low, and most of the variability in diatom community structure remained unexplained in pRDA (see fraction containing residues [U] in Figure 5).
Figure 5 illustrates the importance of considering different dispersal hypotheses in explaining the structure of the periphytic diatom metacommunity in the Cascavel River. Overall, the predictors generated by AEM had the greatest ability to explain the metacommunity (Figure 5).
In general, the species density matrix was better explained than the abundance matrix. The June community was best explained by the predictors, with 18% of the variance explained by the AEM spatial filters (Figure 5). Overland explained the March community best, with 8% explaining species occurrence and 7% explaining species density. Watercourse was best able to explain the dispersal of the community in December 2018 (2%), also associated with the environmental filters and the historical community (Figure 5). Finally, the historical community becomes influential in explaining species dispersal starting in August, with the June community being influenced by the other factors, but not by the historical community (March community), which was impacted by the high accumulated precipitation (Figure 5).

4. Discussion

In this study, we showed that the structure of the diatom metacommunity of a micro-watershed is partially dependent on variables related to local environmental variables, landscape variables, dispersal hypothesis, and historical community. More importantly, we show that the relative importance of these factors varies temporally and depending on how the dispersal hypothesis is generated (see also [19,32,39]). The micro-watershed is located in a predominantly agricultural region of Brazil. The much of the city is structured around open, flat fields or fields with some barriers, mainly characterized by areas used for soybean and corn monocultures [64,65]. In addition, our characterization of land use and occupation of the micro-watershed area showed that the headwaters are surrounded by the urban center of the city. This combination of macro influence with direct influence at the sampling stations summarized the impacts along the micro-watershed, especially in the rainy season, when greater leaching occurs [16,66]. This may be associated with the community not having variation in richness between sampling stations, since diffuse pollution causes homogeneity in the spatial structuring of niches.
The landscape transition between urban and agricultural areas is reflected in different levels of anthropic interference on lotic systems [67]. When comparing micro-watersheds in Paraná State, Peres et al. [68] found that the Cascavel River micro-watershed has high environmental fragility, which refers to the sensitivity and resilience of an ecosystem. This high fragility occurs due to the effect of urban and agricultural effluents present in the river fluvial network in urban areas, which alter the limnological characteristics and promote environmental degradation, reducing forest areas and riparian vegetation [24,69,70]. Medeiros et al. [36], when evaluating the driving factors of the metacommunity in a linear design in the same micro-watershed, recorded the real effect of the continuous river theory, highlighting the effects of urbanization and agriculture on environmental variables, which had the greatest effect on the community.
Our results pointed to several interesting pieces of information about the spatial variability of the diatom community along the micro-watershed. Overall, no single spatial factor explained most of the variation in species density or occurrence in this model. Instead, we highlight the complexity to describe the variability of the metacommunity, as the variance partition did not show a significant percentage of association (see also [71]). The spatial autocorrelation structure without flow connection (overland) can be useful for modeling community abundance, which can actively move both upstream and downstream, overcoming dispersal not only by the lotic environment of connected rivers, but also by other means such as ectozoochory, involving adherence to fish and bird feathers, for instance [29,72].
The size of the spatial extents also interferes with the influence of local and landscape variables on community assembly, based on environmental heterogeneity (see more in [19,36,73,74,75]). The small scale of the micro-watershed probably results in a small spatial gradient that is reflected in the low variation in community richness depending on the sampling station. In fact, spatial and temporal variation in communities have not always been observed (i.e., see results for occurrence matrices), but the relative role of predictors was quite different between periods. This indicates that even with relative temporal stability in the communities, ecological relationships are poorly predicted (see in [76]). Spatial extent is proportional to the importance of niche selection, and the strength of these scale effects depends on the configuration of the environment, dispersal capacity, and niche breadth [74]. However, temporal heterogeneity associated with fluctuating precipitation throughout the year may cause temporal variation in the relative abundance of species. Passy et al. [71] associate this pattern with competition theory, which posits that temporal heterogeneity in resource supply increases biodiversity because it allows species that alternate between dominance and persistence to coexist.
Spatial variables generated using AEM resulted in better estimates, especially when spatial factors were the main drivers of metacommunities. Studies have already shown that it is necessary to consider the connectivity of river networks to model the dispersal routes of metacommunities, as they are important in explaining their effect on community formation [32,48]. However, when considering predictors of hydrologic connectivity across the watercourse, the percentages of total variation explained by spatial and shared effects were reduced, likely reflecting the small spatial gradient [77].
Studies show that environmental gradients may be poorly evident at small scales due to the mass effect [19,20]. This may reflect the high connectivity between sampling stations, which promotes facilitated community dispersal, interfering with the filtering of predictors [77,78]. The influence of metacommunity predictors is more difficult to distinguish at small scales due to the likely effects of niche versus neutral processes on community composition [79]. Possible causes of the complexity to explain these models could be stochastic colonization, ecological drift, and biological interactions or the homogenization of landscape effects itself [19,21,32,73].
Although landscape characteristics, especially the presence of riparian vegetation, are important in providing a favorable environment for community species richness, the results suggest that factors other than land cover may contribute to community restructuring. Previous studies point out that in urban streams, increased surface and subsurface connectivity to drainage networks can lead to water quality impacts, and the presence of forested areas can mitigate these stress pathways [80,81]. Landscape variables end up being suitable driver factors in small-scale models, helping to structure a multidimensional model.
The historical community as a predictor of the metacommunity became significant from the August community onwards. This may be related to the strong influence of precipitation on the March community, which may have needed a period to stabilize again. Foets et al. [82] pointed out that communities took one to two months to re-establish themselves after a significant change in the environment. However, even if the community responds to this temporal variation, in our study, the predominance of both the richness and density of Eunotia specimens is a consequence of the spatial homogeneity, which highlights environmental conditions that favor the occurrence of this genus (acidic and oligotrophic waters) [83,84]. Diatom species such as Eunotia rhomboidea Hustedt, E. botulitropica, Fragilaria pectinalis (Müller) Lyngbye, Gomphonema parvulum, Achnanthidium macrocephalum (Hustedt) Round & Bukhtiyarova, and Sellaphora nigri (De Notaris) Wetzel & Ector did not have specific environmental filtering, occurring throughout the entire micro-watershed; these are species commonly reported in polluted urban environments with high concentrations of solids and nutrients [85,86,87,88].
The deterioration of forested areas, which support the conservation of micro-watersheds within large agricultural extensions, is reflected in the increased concentration of nutrients and metals associated with fertilizers and pesticides. In addition, the sparse riparian vegetation increases siltation and turbidity in the riverbeds. These changes tend to favor medium to large species with cosmopolitan characteristics and competitive advantage, such as Gomphonema lagenula Kützing, Gomphonema exilissimum (Grunow) Lange-Bertalot & Reichardt, and Gomphonema parvulum [85,87,89,90,91,92,93]. These taxa are commonly registered as tolerant to pollution, occurring in conditions of elevated nutrient and light availability [94,95,96]. This environmental filtering supports the homogenizing effect of agriculture, which has been highlighted in other studies [93,97].
Another important factor is that the micro-watershed belongs to a region with soil classified as distropheric red latosol, predominated by Fe and Al ions, and acidic conditions [98]. These and other metals are carried to the aquatic environments, especially in the absence of riparian vegetation, situations in which silting facilitates the accumulation of sediments [99]. These conditions may be harmful to diatoms (see references in [100]), leading to changes in the community structure [101,102]. Thus, the streams studied here are characterized by an acidic pH and high conductivity [103], contributing to the occurrence of species belonging to the genus Eunotia [104].
Some of the effects of urbanization could be observed in stations with sparse and degraded riparian vegetation, such as S2, S4, and S7, reflecting higher concentrations of metals in the surface sediment (Figure 4), especially in March when high accumulated precipitation was observed. The representativeness of some species also helped differentiate these environments, such as Navicula salinicola on S2, Achnanthidium modestiforme and Eunotia yberai on S4 (also see the previously described Nupela semifasciata Amaral, T.Ludwig et Bueno species in [105]), and Eunotia botulitropica on S7 (dominant). These conditions also likely contributed significantly to higher leaching of metals into the surface sediment (Figure 4), reflecting intensified human activities [106,107]. These sampling stations were represented by species such as Navicula salinicola (S2) and Navicula cryptocephala (S2). Therefore, diatom communities are being affected directly or indirectly by a point or diffuse source of pollution. We suggest that other factors need to be considered in future studies, such as the bioavailability of these metals, as well as changes in diatom morphology. On a regional scale, spatial factors can be significant in the distribution of diatoms [108,109], considering land use throughout the micro-watershed [110]. In micro-watersheds, branching rivers contribute significantly to the dispersal dynamics of populations, especially when tributaries are not directly connected to the main river course [31,38]. In these cases, the contribution of the connection can be drastically limited in regions where native vegetation is suppressed, mainly due to anthropic actions [111]. Thus, the temporal maintenance of biodiversity is mainly due to asynchrony, that is, ideal conditions for different species, alternated by the seasonality of climatic events [71,112].
The diatom communities assessed along the micro-watershed presented a variety of growth forms, due to macroinfluence and to the evolutionary process itself, providing dominant profile characteristics of the benthic periphyton biomass [113,114]. The evolutionary aspect of growth may have influence on the comparison results of the methods evaluated, since its dispersal on a small scale (such as the studied micro-watershed), has reached an equilibrium that surpasses the dispersal only by the lotic environment of connected rivers, causing a mass effect [115]. However, it is agreed that time-replicated metacommunity studies are scarce; and that long-term standardized data are essential to clarify the determinants of ecological metacommunities with complex models [116]. Thus, it is necessary that monitoring of ecological communities, such as the one in this study, be continued. Even on a small scale and during one year of sampling, we demonstrated that predictors of different natures act together to explain diatom communities in micro-watersheds.
Due to the limitations of our study, we suggest that future works include a monitoring of the Cascavel River, increasing the number of sampling sites and including data from a larger historical series. It is also suggested to compare these results with other regional micro-watersheds to identify a pattern of anthropic behavior and its influence on lotic environments.

Author Contributions

Conceptualization, G.M. and A.A.P.; methodology, G.M., A.A.P., M.W.W.A., R.G. and M.C.P.; software, G.M. and A.A.P.; validation, S.C.S., T.A.V.L., N.C.B. and R.R.d.R.; formal analysis, G.M.; investigation, G.M., M.W.W.A., R.G. and M.C.P.; resources, S.C.S., N.C.B. and R.R.d.R.; data curation, G.M. and A.A.P.; writing—original draft preparation, G.M.; writing—review and editing, G.M., M.W.W.A., R.G. and M.C.P.; visualization, G.M.; supervision, A.A.P., S.C.S., T.A.V.L., N.C.B. and R.R.d.R.; project administration, R.R.d.R.; funding acquisition, R.R.d.R., S.C.S. and N.C.B. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by CAPES and Fundação Araucária considering both student and researcher scholarships, and other financial resources.

Acknowledgments

The authors wish to thank the Paraná Meteorological Institute (SIMEPAR), for the precipitation data. The authors acknowledge CAPES and Fundação Araucária for continuous support considering both student and researcher scholarships, and other financial resources. We are also grateful to two anonymous reviewers and the editor Milena Milosavljevic for valuable suggestions to improve the manuscript.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A

Epilithic diatom species identified from Cascavel River microbasin.
Table A1. Analysis of variance (Anova) to compare the means of relative abundance of each diatom species according to the sampling stations. Species that had significantly different means across sampling stations, p < 0.05, have p-value highlighted in bold.
Table A1. Analysis of variance (Anova) to compare the means of relative abundance of each diatom species according to the sampling stations. Species that had significantly different means across sampling stations, p < 0.05, have p-value highlighted in bold.
Sampling Station S1S2S3S4S5S6S7S8Fp-Value
Achnanthidium catenatum (J. Bily & Marvan) Lange-Bertalotmean---7.5%0.9%---0.2200.9660.478
standard deviation---14.9%1.8%---
Achnanthidium exiguum (Grunow) D.B.Czarneckimean-------1.1%0,.261.0000.455
standard deviation-------2.2%
Achnanthidium macrocephalum (Hustedt) Round & Bukhtiyarovamean1.1%-3.2%2.7%---2.8%0.2591.1980.341
standard deviation2.2%-3.2%5.3%---3.3%
Achnanthidium minutissimum (Kützing) Czarneckimean5.8%3.6%3.0%5.5%-7.8%5.0%9.1%0.1720.7140.661
standard deviation8.7%4.4%3.7%7.1%-9.1%6.0%9.0%
Achnanthidium modestiforme (Lange-Bertalot) Van de Vijvermean1.9%1.4%-11.1%-1.4%-2.2%0.4773.1300.017
standard deviation2.6%2.9%-9.9%-2.8%-4.3%
Achnanthidium tropicocatenatum G.C.Marquardt, C.E.Wetzel & Ectormean0.4%0.7%3.6%-----0.2741.2910.297
standard deviation0.7%1.5%6.0%-----
Actinella hermes-moreirae Ruwer, Ludwig & Rodriguesmean-0.2%------0.2261.0000.455
standard deviation-0.5%------
Aulacoseira ambigua (Grunow) Simonsenmean-0.5%4.4%-0.3%--2.8%0.3792.0890.085
standard deviation-1.0%3.2%-0.6%--5.6%
Aulacoseira granulata (Ehrenberg) Simonsenmean-2.0%------0.2261.0000.455
standard deviation-3.9%------
Aulacoseira pusilla (Meister) Tuji & Houkimean2.1%-1.3%-----0.2301.0260.439
standard deviation4.3%-1.6%-----
Aulacoseira tenella (Nygaard) Simonsenmean--1.4%-0.8%--0.8%0.1790.7470.635
standard deviation--2.8%-1.6%--1.5%
Brachysira brebissonii Rossmean-0.2%------0.2261.0000.455
standard deviation-0.5%------
Brachysira microcephala (Grunow) Compèremean--1.4%-----0.2261.0000.455
standard deviation--2.9%-----
Caloneis hyalina Hustedtmean---1.0%0.9%---0.2000.8580.552
standard deviation---2.1%1.8%---
Caloneis westii (Smith) Hendeymean------0.7%-0.2261.0000.455
standard deviation------1.4%-
Craticula riparia (Hustedt) Lange-Bertalotmean-0.2%------0.2261.0000.455
standard deviation-0.5%------
Craticula submolesta (Hustedt) Lange-Bertalotmean---2.6%----0.2261.0000.455
standard deviation---5.2%----
Cymbopleura naviculiformis (Auerswald) Krammermean--0.7%-----0.4652.9790.021
standard deviation--0.9%-----
Discostella stelligera (Cleve & Grunow) Houk & Kleemean3.5%0.2%13.7%1.9%-1.0%1.6%7.2%0.3561.8970.115
standard deviation7.0%0.5%17.1%3.7%-1.9%3.1%0.8%
Encyonema minutiforme Krammermean-----2.1%--0.2261.0000.455
standard deviation-----4.2%--
Encyonema neomesianum Krammermean--0.5%--0.7%--0.2020.8670.546
standard deviation--1.0%--1.4%--
Encyonema perpusillum (Cleve-Euler) Mannmean0.4%-------0.2261.0000.455
standard deviation0.7%-------
Encyonema silesiacum (Bleisch) Mannmean-0.2%-0.5%-0.6%-2.0%0.2591.1990.341
standard deviation-0.5%-1.0%-1.1%-3.2%
Eunotia bilunaris (Ehrenberg) Schaarschmidtmean3.4%-7.0%-0.3%5.8%4.2%0.3%0.3301.6860.160
standard deviation6.8%-4.0%-0.6%7.2%6.4%0.5%
Eunotia botulitropica Wetzel & Costamean11.3%2.0%4.9%0.5%11.9%6.3%16.8%0.5%0.2311.0290.437
standard deviation16.6%2.9%7.3%1.1%9.6%11.4%23.6%1,1%
Eunotia georgii Metzeltin & Lange-Bertalotmean----1.4%---0.2261.0000.455
standard deviation----2.7%---
Eunotia intricans Lange-Bertalot & Metzeltinmean-----0.4%--0.2261.0000.455
standard deviation-----0.9%--
Eunotia kruegeri Lange-Bertalotmean0.7%0.8%-1.6%2,4%---0.1680.6930.677
standard deviation1.4%1.5%-3.2%4.8%---
Eunotia meridiana Metzeltin & Lange-Bertalotmean3.5%5.6%6.8%2.0%6.1%3.5%0.8%2.2%0.1350.5350.799
standard deviation4.2%4.9%9.2%2.3%10.4%4.2%1.6%4.3%
Eunotia rhomboidea Hustedtmean7.0%5.6%2.8%1.9%15.8%3.1%5.8%4.8%0.3331.7160.153
standard deviation5.2%4.4%4.1%3.7%14.7%3.6%3.9%6.0%
Eunotia subarcuatoides Alles, Nörpel & Lange-Bertalotmean-----0.7%--0.2261.0000.455
standard deviation-----1.4%--
Eunotia tropico-arcus Metzeltin & Lange-Bertalotmean-----0.3%0.2%-0.2010.8600.551
standard deviation-----0.6%0.5%-
Eunotia veneris (Kützing) De Tonimean19.6%-1.6%-1.4%-0.2%0.7%0.2391.0750.409
standard deviation36.7%-3.1%-2.7%-0.5%1.4%
Eunotia yberai Frenguellimean---2.0%----0.4652.9830.021
standard deviation---2.3%----
Fragilaria gracilis Østrupmean--2.4%-0.3%0.9%--0.3621.9430.107
standard deviation--2.9%-0.6%1.8%--
Fragilaria pectinalis (Müller) Lyngbyemean-1.2%1.3%--2.7%0.8%1.0%0.1520.6140.739
standard deviation-2.4%1.5%--5.4%1.6%2.0%
Fragilaria tenera var. tenera (Smith) Lange-Bertalotmean-----2.0%-1.4%0.2020.8660.547
standard deviation-----3.9%-2.8%
Fragilaria sp.mean--0.5%--3.1%--0.2180.9570.483
standard deviation--1.0%--6.3%--
Fragilariforma javanica (Hustedt) Wetzel, Morales & Ectormean----1.0%---0.2261.0000.455
standard deviation----2.0%---
Frustulia acidophilissima Wydrzycka & Lange-Bertalotmean1.4%--0.5%1.0%1.5%1.7%0.5%0.1070.4100.887
standard deviation2.7%--1.1%2.0%2.9%3.4%1.1%
Frustulia crassinervia (Brébisson ex W.Smith) Lange-Bertalot & Krammermean1.2%--1.5%5.4%--0.7%0.3311.6980.157
standard deviation0.9%--3.1%7.2%--1.4%
Frustulia guayanensis Metzeltin & Lange Bertalotmean-0.8%0.4%-0.7%1.7%--0.1710.7060.667
standard deviation-1.5%0.8%-1.4%3.4%--
Frustulia neomundana Lange-Bertalot & Rumrichmean------0.8%-0.2261.0000.455
standard deviation------1.6%-
Frustulia saxonica Rabenhorstmean-----0.5%-0.5%0.2000.8580.552
standard deviation-----1.0%-1.1%
Frustulia undosa Metzeltin & Lange-Bertalotmean0.2%-------0.2261.0000.455
standard deviation0.4%-------
Frustulia vulgaris (Thwaites) De Tonimean----0.5%1.0%--0.2060.8900.530
standard deviation----1.0%2.1%--
Frustulia weinholdii Hustedtmean0.4%------0.9%0.3681.9930.098
standard deviation0.7%------1.1%
Geissleria punctifera (Hustedt) Metzeltin, Lange-Bertalot e García-Rodríguezmean-0.2%--0.5%--1.1%0.1950.8300.573
standard deviation-0.5%--1.0%--2.3%
Gogorevia exilis(Kütz.) Kulikovskiy and Kociolekmean-------1.4%0.2261.0000.455
standard deviation-------2.8%
Gomphonema angustatum (Kützing) Rabenhorstmean-------1.6%0.2261.0000.455
standard deviation-------3.3%
Gomphonema exilissimum (Grunow) Lange-Bertalot & Reichardtmean--0.6%--6.3%-1.5%0.2090.9040.520
standard deviation--1.1%--12.5%-3.1%
Gomphonema guaraniarum Metzeltin & Lange-Bertalotmean-----1.1%-0.8%0.2010.8640.548
standard deviation-----2.1%-1.5%
Gomphonema graciledictum Reichardtmean--0.8%----0.3%0.2070.8940.527
standard deviation--1.6%----0.7%
Gomphonema lagenula Kützingmean3.3%7.0%3.8%1.0%-2.7%0.8%4.5%0.2020.8690.544
standard deviation4.1%7.4%6.7%2.1%-4.2%1.6%7.3%
Gomphonema naviculoides Smithmean0.4%-0.5%--1.3%--0.1950.8290.574
standard deviation0.7%-1.0%--2.7%--
Gomphonema obtusatum (Kützing) Grunowmean--0.4%1.1%----0.2090.9080.517
standard deviation--0.8%2.1%----
Gomphonema parvulum (Kützing) Kützingmean5.9%13.0%19.7%2.8%9.6%22.7%3.1%10.3%0.2361.0600.418
standard deviation5.8%9.4%33.2%3.2%8.0%11.7%6.3%11.4%
Gomphonema pumilum (Grunow) Reichardt & Lange-Bertalotmean-2.2%0.7%--0.9%2.7%-0.2080.9000.522
standard deviation-2.9%0.8%--1.8%5.4%-
Halamphora montana (Krasske) Levkovmean-0.7%-1.0%-0.3%3.9%0.7%0.1850.7810.610
standard deviation-1.5%-2.1%-0.6%7.8%1.4%
Hantzschia amphioxys (Ehrenberg) Grunowmean-0.2%---0.3%--0.2000.8590.552
standard deviation-0.5%---0.6%--
Humidophila contenta (Grunow) Lowe et al.mean0.4%4.8%-1.5%--7.4%-0.3091.5330.204
standard deviation0.9%8.9%-3.1%--8.5%-
Humidophila subtropica (Metzeltin, Lange-Bertalot & García-Rodríguez) Lowe et al.mean-1.5%------0.2261.0000.455
standard deviation-2.9%------
Iconella tenuissima (Hustedt) D.Kapustin & Kulikovskiymean---1.6%-0.5%-0.3%0.3421.7790.138
standard deviation---2.0%-1.0%-0.7%
Luticola acidoclinata Lange-Bertalotmean-0.2%--2.8%1.9%-0.9%0.2621.2140.333
standard deviation-0.5%--3.8%3.9%-1.1%
Luticola goeppertiana (Bleisch) Mannmean--0.4%-1.5%--2.5%0.1940.8260.576
standard deviation--0.8%-2.9%--5.0%
Luticola hustedtii Levkov, Metzeltin & Pavlovmean0.2%-------0.2261.0000.455
standard deviation0.4%-------
Luticola permuticoides Metzeltin e Lange-Bertalotmean-0.7%---1.0%--0.2020.8660.546
standard deviation-1.4%---1.9%--
Navicula angusta Grunowmean-0.2%--1.1%---0.2150.9380.496
standard deviation-0.5%--2.1%---
Navicula cryptocephala Kützingmean1.5%7.5%2.3%--0.7%-4.3%0.4753.1020.018
standard deviation1.5%6.8%2.8%--1.4%-3.7%
Navicula cryptotenella Lange-Bertalotmean1.3%0.2%1.5%-10.0%-1.2%-0.4723.0630.019
standard deviation2.0%0.5%2.2%-10.4%-1.5%-
Navicula eichhorniaephila Manguin ex Kociolek & Reviersmean---0.5%----0.2261.0000.455
standard deviation---1.1%----
Navicula leptostriata E.G.Jørgensenmean--0.4%-----0.2261.0000.455
standard deviation--0.8%-----
Navicula radiosa Kützingmean-2.0%0.3%-0.5%0.6%--0.1860.7830.608
standard deviation-4.1%0.7%-1.0%1.3%--
Navicula rostellata Kützingmean-----1.1%-0.3%0.2150.9390.496
standard deviation-----2.3%-0.5%
Navicula salinicola Hustedtmean0.2%3.6%------0.4312.5960.038
standard deviation0.4%4.4%------
Navicula symmetrica Patrickmean-1.7%-----1.4%0.2010.8600.551
standard deviation-3.4%-----2.8%
Neidium essequiboanum Metzeltin & Krammermean-----0.5%-0.8%0.2030.8760.540
standard deviation-----1.0%-1.7%
Neidium gracile f. aequale Hustedtmean0.4%-------0.2261.0000.455
standard deviation0.7%-------
Neidium iridis (Ehrenberg) Clevemean----0.5%---0.2261.0000.455
standard deviation----1.0%---
Neidium tenuissimum Hustedtmean0.7%-----1.4%-0.2050.8860.533
standard deviation1.4%-----2.7%-
Nitzschia amphibia Grunowmean0.7%-------0.2261.0000.455
standard deviation1.4%-------
Nitzschia brevissima Grunowmean---0.5%----0.2261.0000.455
standard deviation---1.0%----
Nitzschia palea (Kützing) Smithmean1.6%4.5%0.9%2.1%1.0%0.6%-1.3%0.2681.2550.313
standard deviation2.6%3.7%1.0%4.1%2.0%1.1%-2.1%
Nitzschia perminuta Grunow in Van Heurckmean0.4%-------0.2261.0000.455
standard deviation0.7%-------
Nupela pardinhoensis Bes, Torgan & Ector in Bes et al.mean2.4%2.6%-3.9%2.4%2.1%8.3%2.1%0.2230.9820.467
standard deviation2.8%3.6%-3.5%4.8%4.2%10.1%3.1%
Nupela praecipuoides Tremarin & Ludwigmean1.1%--2.0%4.3%1.9%5.4%7.3%0.2841.3630.266
standard deviation2.1%--2.9%6.5%3.9%4.5%8.7%
Nupela semifasciata Amaral, T.Ludwig et Buenomean2.7%2.7%-2.4%1.5%2.6%8.3%2.4%0.1410.5640.777
standard deviation5.4%5.4%-3.5%2.9%3.9%14.4%4.9%
Nupela wellneri (Lange-Bertalot) Lange-Bertalotmean------1.1%-0.2261.0000.455
standard deviation------2.3%-
Orthoseira roeseana (Rabenhorst) Pfitzermean----2.0%---0.2261.0000.455
standard deviation----4.1%---
Pinnularia brauniana (Grunow) Studnickamean--0.4%-----0.2261.0000.455
standard deviation--0.8%-----
Pinnularia divergens W.Smithmean0.4%-1.4%0.5%----0.1950.8320.572
standard deviation0.7%-2.9%1.1%----
Pinnularia graciloides var. latecapitata Metzeltin & Krammermean-----0.6%--0.2261.0000.455
standard deviation-----1.1%--
Pinnularia laucensis Lange-Bertalotmean0.4%---0.7%1.5%-0.5%0.1720.7120.662
standard deviation0.7%---1.4%3.0%-1.1%
Pinnularia microstauron var. rostrata Krammermean0.4%-------0.2261.0000.455
standard deviation0.7%-------
Pinnularia obscura Krasskemean----2.0%---0.2261.0000.455
standard deviation----4.1%---
Pinnularia subanglica Krammermean-0.8%------0.2261.0000.455
standard deviation-1.5%------
Pinnularia subgibba Krammermean---0.5%----0.2261.0000.455
standard deviation---1.0%----
Pinnularia subinterrupta Krammer & Schroetermean----0.3%---0.2261.0000.455
standard deviation----0.6%---
Pinnularia sp.mean1.4%0.2%------0.2170.9510.488
standard deviation2.7%0.5%------
Pinnularia sp.2mean0.2%-------0.2261.0000.455
standard deviation0.4%-------
Pinnularia tabellaria Ehrenbergmean1.1%-------0.2261.0000.455
standard deviation2.1%-------
Placogeia kriegeri (K.Krasske) Bukhtiyarovamean---0.5%-0.3%--0.2040.8800.537
standard deviation---1.0%-0.6%--
Placoneis elginensis (W.Gregory) E.J.Coxmean----0.7%---0.2261.0000.455
standard deviation----1.4%---
Placoneis hambergii (Hustedt) K.Brudermean-------1.1%0.2261.0000.455
standard deviation-------2.2%
Planothidium frequentissimum (Lange-Bertalot) Lange-Bertalotmean----0.7%---0.2261.0000.455
standard deviation----1.4%---
Psammothidium hustedtii (Krasske) Mayamamean------0.2%-0.2261.0000.455
standard deviation------0.5%-
Rhopalodia gibberula var. vanheurckii O. Müllermean--0.4%-----0.2261.0000.455
standard deviation--0.8%-----
Sellaphora laevissima (Kützing) D.G.Mannmean----0.5%---0.2261.0000.455
standard deviation----1.0%---
Sellaphora nigri (De Notaris) Wetzel & Ectormean9.1%9.3%-18.8%6.0%1.0%14.5%4.4%0.3241.6410.172
standard deviation8.0%10.7%-20.9%5.6%1.9%11.5%5.4%
Sellaphora sassiana (Metzeltin & Lange-Bertalot)mean--2.4%1.5%----0.2020.8700.544
standard deviation--4.8%3.1%----
Sellaphora saugerresii (Desmazières) Wetzel & Mannmean-7.9%1.2%10.1%0.3%-3.0%3.7%0.1590.6490.712
standard deviation-15.9%1.8%20.2%0.6%-6.0%5.1%
Sellaphora sp.mean---1.8%----0.2261.0000.455
standard deviation---3.6%----
Sellaphora sp.2mean-0.1%------0.2261.0000.455
standard deviation-0.2%------
Spicaticribra kingstonii Johansen, Kociolek & Lowemean--1.0%----0.3%0.6165.5090.001
standard deviation--0.7%----0.5%
Stauroneis gracilis Ehrenbergmean-----0.6%--0.2261.0000.455
standard deviation-----1.1%--
Surirella angusta Kützingmean0.2%--1.5%----0.2190.9610.481
standard deviation0.4%--3.1%----
Surirella roba Leclercqmean-0.4%--0.3%---0.2000.8590.551
standard deviation-0.7%--0.6%---
Ulnaria ulna (Nitzsch) Compèremean0.2%0.2%0.1%1.5%---2.1%0.2471.1260.380
standard deviation0.4%0.5%0.3%3.1%---3.1%

References

  1. Tang, W.; Pei, Y.; Zheng, H.; Zhao, Y.; Shu, L.; Zhang, H. Twenty years of China’s water pollution control: Experiences and challenges. Chemosphere 2022, 295, 133875. [Google Scholar] [CrossRef]
  2. Elsheikh, A.H.; Panchal, H.N.; Sengottain, S.; Alsaleh, N.A.; Ahmadein, M. Applications of Heat Exchanger in Solar Desalination: Current Issues and Future Challenges. Water 2022, 14, 852. [Google Scholar] [CrossRef]
  3. Elsheikh, A.H.; Sharshir, S.W.; Mostafa, M.E.; Essa, F.A.; Ahmed Ali, M.K. Applications of nanofluids in solar energy: A review of recent advances. Renew. Sustain. Energy Rev. 2018, 82, 3483–3502. [Google Scholar] [CrossRef]
  4. Chiquito Gesualdo, G.; Sone, J.S.; de Galvão, C.O.; Martins, E.S.; Montenegro, S.M.G.L.; Tomasella, J.; Mendiondo, E.M. Unveiling water security in Brazil: Current challenges and future perspectives. Hydrol. Sci. J. 2021, 66, 759–768. [Google Scholar] [CrossRef]
  5. Getirana, A.; Libonati, R.; Cataldi, M. Brazil is in water crisis—It needs a drought plan. Nature 2021, 600, 218–220. [Google Scholar] [CrossRef]
  6. Kakore, B.G.; Mamun, M.; Lee, S.-J.; An, K.-G. Land-Use Pattern as a Key Factor Determining the Water Quality, Fish Guilds, and Ecological Health in Lotic Ecosystems of the Asian Monsoon Region. Water 2022, 14, 2765. [Google Scholar] [CrossRef]
  7. Albert, J.S.; Destouni, G.; Duke-Sylvester, S.M.; Magurran, A.E.; Oberdorff, T.; Reis, R.E.; Winemiller, K.O.; Ripple, W.J. Scientists’ warning to humanity on the freshwater biodiversity crisis. Ambio 2021, 50, 85–94. [Google Scholar] [CrossRef]
  8. Zhang, Y.; Peng, C.; Huang, S.; Wang, J.; Xiong, X.; Li, D. The relative role of spatial and environmental processes on seasonal variations of phytoplankton beta diversity along different anthropogenic disturbances of subtropical rivers in China. Environ. Sci. Pollut. Res. 2019, 26, 1422–1434. [Google Scholar] [CrossRef]
  9. Pineda, A.; Bortolini, J.C.; Rodrigues, L.C. Effects of space and environment on phytoplankton distribution in subtropical reservoirs depend on functional features of the species. Aquat. Sci. 2022, 84, 5. [Google Scholar] [CrossRef]
  10. Zhou, S.; Wu, N.; Zhang, M.; Peng, W.; He, F.; Guo, K.; Yan, S.; Zou, Y.; Qu, X. Local environmental, geo-climatic and spatial factors interact to drive community distributions and diversity patterns of stream benthic algae, macroinvertebrates and fishes in a large basin, Northeast China. Ecol. Indic. 2020, 117, 106673. [Google Scholar] [CrossRef]
  11. Huszar, V.L.M.; Nabout, J.C.; Appel, M.O.; Santos, J.B.O.; Abe, D.S.; Silva, L.H.S. Environmental and not spatial processes (directional and non-directional) shape the phytoplankton composition and functional groups in a large subtropical river basin. J. Plankton Res. 2015, 37, 1190–1200. [Google Scholar] [CrossRef] [Green Version]
  12. Leibold, M.A.; Holyoak, M.; Mouquet, N.; Amarasekare, P.; Chase, J.M.; Hoopes, M.F.; Holt, R.D.; Shurin, J.B.; Law, R.; Tilman, D.; et al. The metacommunity concept: A framework for multi-scale community ecology. Ecol. Lett. 2004, 7, 601–613. [Google Scholar] [CrossRef]
  13. Lansac-Tôha, F.M.; Bini, L.M.; Heino, J.; Meira, B.R.; Segovia, B.T.; Pavanelli, C.S.; Bonecker, C.C.; Deus, C.P.; Benedito, E.; Alves, G.M.; et al. Scale-dependent patterns of metacommunity structuring in aquatic organisms across floodplain systems. J. Biogeogr. 2021, 48, 872–885. [Google Scholar] [CrossRef]
  14. Rodriguez, J.M.M. Geoecologia das Paisagens: Uma visão Geossistêmica da Análise Ambiental, 6th ed.; Rodriguez, J.M.M., da Silva, E.V., de Paula Brito Cavalcanti, A., Eds.; Imprensa Universitária: Fortaleza, Brazil, 2021; Volume 41, ISBN 8572821481. [Google Scholar]
  15. Stenger-Kovács, C.; Lengyel, E.; Sebestyén, V.; Szabó, B. Effects of land use on streams: Traditional and functional analyses of benthic diatoms. Hydrobiologia 2020, 847, 2933–2946. [Google Scholar] [CrossRef]
  16. Li, C.; Wang, J.; Jiang, Y.; Bai, Y.; Cheng, L. The evolution of different dissolved organic matter components and release characteristics of heavy metals in leaching process from sewage sludge under simulated rain. Environ. Sci. Pollut. Res. 2022, 29, 86651–86664. [Google Scholar] [CrossRef]
  17. Docile, T.; Rosa, D.C.O.; Figueiró, R.; Nessimian, J. Urbanisation alters the flow of energy through stream food webs. Insect Conserv. Divers. 2016, 9, 416–426. [Google Scholar] [CrossRef]
  18. Barnum, T.R.; Weller, D.E.; Williams, M. Urbanization reduces and homogenizes trait diversity in stream macroinvertebrate communities. Ecol. Appl. 2017, 27, 2428–2442. [Google Scholar] [CrossRef]
  19. Heino, J.; Melo, A.S.; Siqueira, T.; Soininen, J.; Valanko, S.; Bini, L.M. Metacommunity organisation, spatial extent and dispersal in aquatic systems: Patterns, processes and prospects. Freshw. Biol. 2015, 60, 845–869. [Google Scholar] [CrossRef]
  20. Soininen, J.; Teittinen, A. Fifteen important questions in the spatial ecology of diatoms. Freshw. Biol. 2019, 64, 2071–2083. [Google Scholar] [CrossRef] [Green Version]
  21. Soininen, J. Environmental and spatial control of freshwater diatoms—A review. Diatom Res. 2007, 22, 473–490. [Google Scholar] [CrossRef]
  22. Costa, A.P.T.; Schneck, F. Diatoms as indicators in running waters: Trends of studies on biological assessment and monitoring. Environ. Monit. Assess. 2022, 194, 695. [Google Scholar] [CrossRef] [PubMed]
  23. Zorzal-Almeida, S.; Lehmkuhl, E.A.; Marquardt, G.C.; Morais, K.S.; Trancoso, M.S.; Silva-Lehmkuhl, A.M. 177 years of diatom studies in Brazil: Knowledge, gaps, and perspectives. An. Acad. Bras. Cienc. 2022, 94, e20210959. [Google Scholar] [CrossRef] [PubMed]
  24. Hu, J.; Song, Z.; Zhou, J.; Soininen, J.; Tan, L.; Cai, Q.; Tang, T. Differences in diversity and community assembly processes between planktonic and benthic diatoms in the upper reach of the Jinsha River, China. Hydrobiologia 2022, 849, 1577–1591. [Google Scholar] [CrossRef]
  25. Maurya, P.K.; Ali, S.A.; Alharbi, R.S.; Yadav, K.K.; Alfaisal, F.M.; Ahmad, A.; Ditthakit, P.; Prasad, S.; Jung, Y.-K.; Jeon, B.-H. Impacts of Land Use Change on Water Quality Index in the Upper Ganges River near Haridwar, Uttarakhand: A GIS-Based Analysis. Water 2021, 13, 3572. [Google Scholar] [CrossRef]
  26. Khouni, I.; Louhichi, G.; Ghrabi, A. Use of GIS based Inverse Distance Weighted interpolation to assess surface water quality: Case of Wadi El Bey, Tunisia. Environ. Technol. Innov. 2021, 24, 101892. [Google Scholar] [CrossRef]
  27. Setia, R.; Lamba, S.; Chander, S.; Kumar, V.; Singh, R.; Litoria, P.K.; Singh, R.P.; Pateriya, B. Spatio-temporal variations in water quality, hydrochemistry and its controlling factors in a perennial river in India. Appl. Water Sci. 2021, 11, 169. [Google Scholar] [CrossRef]
  28. Shurin, J.B.; Cottenie, K.; Hillebrand, H. Spatial autocorrelation and dispersal limitation in freshwater organisms. Oecologia 2009, 159, 151–159. [Google Scholar] [CrossRef]
  29. Manning, F.S.; Curtis, P.J.; Walker, I.R.; Pither, J. Potential long-distance dispersal of freshwater diatoms adhering to waterfowl plumage. Freshw. Biol. 2021, 66, 1136–1148. [Google Scholar] [CrossRef]
  30. Jamoneau, A.; Passy, S.I.; Soininen, J.; Leboucher, T.; Tison-Rosebery, J. Beta diversity of diatom species and ecological guilds: Response to environmental and spatial mechanisms along the stream watercourse. Freshw. Biol. 2018, 63, 62–73. [Google Scholar] [CrossRef] [Green Version]
  31. Liu, J.; Soininen, J.; Han, B.-P.; Declerck, S.A.J. Effects of connectivity, dispersal directionality and functional traits on the metacommunity structure of river benthic diatoms. J. Biogeogr. 2013, 40, 2238–2248. [Google Scholar] [CrossRef]
  32. Padial, A.A.; Ceschin, F.; Declerck, S.A.J.; De Meester, L.; Bonecker, C.C.; Lansac-Tôha, F.A.; Rodrigues, L.; Rodrigues, L.C.; Train, S.; Velho, L.F.M.; et al. Dispersal Ability Determines the Role of Environmental, Spatial and Temporal Drivers of Metacommunity Structure. PLoS ONE 2014, 9, e111227. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  33. Bortolini, J.C.; Pineda, A.; Rodrigues, L.C.; Jati, S.; Velho, L.F.M. Environmental and spatial processes influencing phytoplankton biomass along a reservoirs-river-floodplain lakes gradient: A metacommunity approach. Freshw. Biol. 2017, 62, 1756–1767. [Google Scholar] [CrossRef]
  34. Zorzal-Almeida, S.; Soininen, J.; Bini, L.M.; Bicudo, D.C. Local environment and connectivity are the main drivers of diatom species composition and trait variation in a set of tropical reservoirs. Freshw. Biol. 2017, 62, 1551–1563. [Google Scholar] [CrossRef]
  35. De Oliveira, P.H.F.; Machado, K.B.; Teresa, F.B.; Heino, J.; Nabout, J.C. Spatial processes determine planktonic diatom metacommunity structure of headwater streams. Limnologica 2020, 84, 125813. [Google Scholar] [CrossRef]
  36. Medeiros, G.; Padial, A.A.; Wedig Amaral, M.W.; Ludwig, T.A.V.; Bueno, N.C. Environmental variables likely influence the periphytic diatom community in a subtropical lotic environment. Limnologica 2020, 80, 125718. [Google Scholar] [CrossRef]
  37. Tonkin, J.D.; Heino, J.; Altermatt, F. Metacommunities in river networks: The importance of network structure and connectivity on patterns and processes. Freshw. Biol. 2018, 63, 1–5. [Google Scholar] [CrossRef] [Green Version]
  38. Tonkin, J.D.; Altermatt, F.; Finn, D.S.; Heino, J.; Olden, J.D.; Pauls, S.U.; Lytle, D.A. The role of dispersal in river network metacommunities: Patterns, processes, and pathways. Freshw. Biol. 2018, 63, 141–163. [Google Scholar] [CrossRef] [Green Version]
  39. Wojciechowski, J.; Heino, J.; Bini, L.M.; Padial, A.A. The strength of species sorting of phytoplankton communities is temporally variable in subtropical reservoirs. Hydrobiologia 2017, 800, 31–43. [Google Scholar] [CrossRef]
  40. Alvares, C.A.; Stape, J.L.; Sentelhas, P.C.; Gonc, L.D.M.; Sparovek, G. Köppen’s climate classification map for Brazil Clayton. Meteorol. Z. 2013, 22, 711–728. [Google Scholar] [CrossRef]
  41. Cembranel, A.S.; Sampaio, S.C.; Remor, M.B.; Gotardo, J.T.; Rosa, P.M.D. Geochemical background in an oxisol. Eng. Agrícola 2017, 37, 565–573. [Google Scholar] [CrossRef]
  42. de Cascavel, P.M. Recuperação Ambiental da Bacia Hidrográfica do Rio Cascavel; FUNDETEC: Sincelejo, Colombia, 1995. [Google Scholar]
  43. Instituto Brasileiro de Geografia e Estatística. Manual Técnico de Uso da Terra; Instituto Brasileiro de Geografia e Estatística: Rio de Janeiro, Brazil, 2013; ISBN 9788524043079.
  44. APHA. Standard Methods for the Examination of Water and Wastewater, 22nd ed.; Franson, M.A.H., Ed.; Copyright: Washington, DC, USA, 2012; ISBN 978-0875530130. [Google Scholar]
  45. CETESB. ANA Guia Nacional de Coleta e Preservação de Amostras, 2nd ed.; Brandão, C.J., Botelho, M.J.C., Sato, M.I.Z., Lamparelli, M.C., Eds.; CETESB: São Paulo, Brazil, 2011.
  46. USEPA 3050B. Acid Digestion of Sediments, Sludges, and Soils; USEPA: Washington, DC, USA, 1996; pp. 1–12.
  47. Borcard, D.; Legendre, P. All-scale spatial analysis of ecological data by means of principal coordinates of neighbour matrices. Ecol. Modell. 2002, 153, 51–68. [Google Scholar] [CrossRef]
  48. Blanchet, F.G.; Legendre, P.; Borcard, D. Modelling directional spatial processes in ecological data. Ecol. Modell. 2008, 215, 325–336. [Google Scholar] [CrossRef]
  49. Simonsen, R. The Diatom Plankton of the Indian Ocean Expedition of RV “Meteor”, 1964–1965; Borntraeger: Berlin, Germany, 1974; ISBN 3443200192 9783443200190. [Google Scholar]
  50. Moreira-Filho, H.; Valente Mreira, I. Avaliacao taxonomica e ecologica das Diatomaceas (Bacillariophyceae) epifitas em algas pluricelulares obtidas nos litorais dos Estados do Parana, Santa Catarina e Sao Paulo. Bol. Mus. Bot. Munic. 1981, 1, 17. [Google Scholar]
  51. Kobayasi, H.; Mayama, S. Most pollution-tolerant diatoms of severely polluted rivers in the vicinity of Tokyo. Jpn. J. Phycol. 1982, 30, 188–196. [Google Scholar]
  52. Pappas, J.L.; Stoermer, E.F. Quantitative method for determining a representative algal sample count. J. Phycol. 1996, 32, 693–696. [Google Scholar] [CrossRef] [Green Version]
  53. Battarbee, R.W. Diatom Analysis. In Handbook of Holocene Palaeoecology and Palaeohydrology; Berglund, B.E., Ralska-Jasiewiczowa, M., Eds.; International Geological Correlation Programme; Wiley: New York, NY, USA, 1986; p. 869. ISBN 9780471906919. [Google Scholar]
  54. Lobo, E.A.; Katoh, K.; Aruga, Y. Response of epilithic diatom assemblages to water pollution in rivers located in the Tokyo Metropolitan area, Japan. Freshw. Biol. 1995, 34, 191–204. [Google Scholar] [CrossRef]
  55. Legendre, P.; Legendre, L. Numerical Ecology, Volume 24. Dev. Environ. Model. 1988, 24, 870. [Google Scholar] [CrossRef]
  56. Anderson, M.J. A new method for non-parametric multivariate analysis of variance. Austral Ecol. 2001, 26, 32–46. [Google Scholar] [CrossRef]
  57. Gower, J.C. Some distance properties of latent root and vector methods used in multivariate analysis. Biometrika 1966, 53, 325–338. [Google Scholar] [CrossRef]
  58. Tuomisto, H.; Ruokolainen, K. Analyzing or explaining beta diversity? Understanding the targets of different methods of analysis. Ecology 2006, 87, 2697–2708. [Google Scholar] [CrossRef]
  59. Radbruch, L.; De Lima, L.; Knaul, F.; Wenk, R.; Ali, Z.; Bhatnaghar, S.; Blanchard, C.; Bruera, E.; Buitrago, R.; Burla, C.; et al. Redefining Palliative Care—A New Consensus-Based Definition. J. Pain Symptom Manag. 2020, 60, 754–764. [Google Scholar] [CrossRef]
  60. Ntongani, W.A.; Andrew, S.M.; R Development Core Team. R: A Language and Environment for Statistical Computing; R Development Core Team: Vienna, Austria, 2011; Volume 1, ISBN 3900051070. [Google Scholar]
  61. Oksanen, J.; Blanchet, F.G.; Friendly, M.; Kindt, R.; Legendre, P.; Mcglinn, D.; Minchin, P.R.; O’hara, R.B.; Simpson, G.L.; Solymos, P.; et al. Package “Vegan” Title Community Ecology Package Version 2.5-7. 2020. [Google Scholar]
  62. Dray, S.; Pélissier, R.; Couteron, P.; Fortin, M.-J.; Legendre, P.; Peres-Neto, P.R.; Bellier, E.; Bivand, R.; Blanchet, F.G.; De Cáceres, M.; et al. Community ecology in the age of multivariate multiscale spatial analysis. Ecol. Monogr. 2012, 82, 257–275. [Google Scholar] [CrossRef]
  63. Bivand, R. R Packages for Analyzing Spatial Data: A Comparative Case Study with Areal Data. Geogr. Anal. 2022, 54, 488–518. [Google Scholar] [CrossRef]
  64. Cima, E.G.; Uribe-Opazo, M.A.; Da Rocha Junior, W.F.; De Souza Fragoso, R.M. A spatial analysis of western paranÁ: Scenarios for regional development. Rev. Bras. Gest. E Desenvolv. Reg. 2021, 17, 151–164. [Google Scholar]
  65. Guicho, R.; Ferreira, J.H.D.; Medeiros, G.; Siqueira, J.A.C.; de Souza, S.N.M.; Prior, M. Method for estimating the wind power micro and minigeneration applied to a city with a subtropical climate in South America. Res. Soc. Dev. 2021, 10, e105101220009. [Google Scholar] [CrossRef]
  66. da Silva, J.N.; da Cunha, L.S.; de Oliveira, N.E.C. O Uso De Parâmetros Físico-Químicos Na Delimitação De Contaminação Por Lixiviado Em Áreas Degradadas Por Resíduos Sólidos Urbanos. Uma Revisão Sistemática De Literatura Com Ênfase Em Metanálise. Rev. Bras. Geogr. Física 2022, 15, 1587. [Google Scholar] [CrossRef]
  67. De Almeida, T.P.; Macena, D.Â.; Simões, J.S.T.; Mareco, E.A.; Calciolari Rossi, R.; Favareto, A.P.A. Análise de parâmetros de qualidade da água e teste de genotoxicidade em peixes da bacia hidrográfica do rio Pirapozinho—SP, Brasil. Res. Soc. Dev. 2022, 11, e46711319309. [Google Scholar] [CrossRef]
  68. Peres, K.K.; Guicho, R.; Medeiros, G.; Amaral, M.W.W.; da Silva, T.T.; Pilatti, M.C.; Prior, M.; Bueno, N.C. Environmental fragility as an indicator of the risk of contamination by human action in watersheds used for public supply in western Paraná, Brazil. Environ. Earth Sci. 2022, 81, 486. [Google Scholar] [CrossRef]
  69. De Carvalho, L.H.W.; Laudares, S.; Libório, M.P.; Ekel, M.P.; de Oliveira Marques, R. Urban Settlement Spatial Analysis in Permanent Preservation Area of M’boicy Watershed River, Foz do Iguassu City in Brazil. Int. J. Geosci. 2016, 7, 1222–1231. [Google Scholar] [CrossRef] [Green Version]
  70. Odeh, T.; Mohammad, A.H.; Pradhanang, S.M.; Ismail, M.; Rödiger, T. GIS-based analytical modeling on evaluating impacts of urbanization in Amman water resources, Jordan. Environ. Earth Sci. 2022, 81, 160. [Google Scholar] [CrossRef]
  71. Passy, S.I.; Larson, C.A.; Jamoneau, A.; Budnick, W.; Heino, J.; Leboucher, T.; Tison-Rosebery, J.; Soininen, J. Biogeographical patterns of species richness and abundance distribution in stream diatoms are driven by climate and water chemistry. Am. Nat. 2018, 192, 605–617. [Google Scholar] [CrossRef] [Green Version]
  72. Erős, T.; Lowe, W.H. The Landscape Ecology of Rivers: From Patch-Based to Spatial Network Analyses. Curr. Landsc. Ecol. Rep. 2019, 4, 103–112. [Google Scholar] [CrossRef] [Green Version]
  73. Vilmi, A.; Karjalainen, S.M.; Hellsten, S.; Heino, J. Bioassessment in a metacommunity context: Are diatom communities structured solely by species sorting? Ecol. Indic. 2016, 62, 86–94. [Google Scholar] [CrossRef]
  74. Viana, D.S.; Chase, J.M. Spatial scale modulates the inference of metacommunity assembly processes. Ecology 2019, 100, e02576. [Google Scholar] [CrossRef] [PubMed]
  75. He, S.; Soininen, J.; Chen, K.; Wang, B. Environmental Factors Override Dispersal-Related Factors in Shaping Diatom and Macroinvertebrate Communities Within Stream Networks in China. Front. Ecol. Evol. 2020, 8, 141. [Google Scholar] [CrossRef]
  76. Ferraz Sampaio, M.; Ferreira, P.C.; Veiga Ludwig, T.A.; Padial, A.A. Temporal stability in beta diversity does not guarantee surrogacy or compositional stability in a micro-phytoplankton metacommunity. Limnetica 2023, 42, 1. [Google Scholar] [CrossRef]
  77. Marquardt, G.C.; Padial, A.A.; Bicudo, C.E.M. Variance partitioning of deconstructed tropical diatom communities in reservoirs cascade. Aquat. Sci. 2018, 80, 17. [Google Scholar] [CrossRef]
  78. Vilmi, A.; Tolonen, K.T.; Karjalainen, S.M.; Heino, J. Metacommunity structuring in a highly-connected aquatic system: Effects of dispersal, abiotic environment and grazing pressure on microalgal guilds. Hydrobiologia 2017, 790, 125–140. [Google Scholar] [CrossRef]
  79. Gilbert, B.; Lechowicz, M.J. Neutrality, niches, and dispersal in a temperate forest understory. Proc. Natl. Acad. Sci. USA 2004, 101, 7651–7656. [Google Scholar] [CrossRef] [Green Version]
  80. Kaushal, S.S.; Belt, K.T. The urban watershed continuum: Evolving spatial and temporal dimensions. Urban Ecosyst. 2012, 15, 409–435. [Google Scholar] [CrossRef]
  81. Mutinova, P.T.; Kahlert, M.; Kupilas, B.; McKie, B.G.; Friberg, N.; Burdon, F.J. Benthic Diatom Communities in Urban Streams and the Role of Riparian Buffers. Water 2020, 12, 2799. [Google Scholar] [CrossRef]
  82. Foets, J.; Wetzel, C.E.; Teuling, A.J.; Pfister, L. Temporal and spatial variability of terrestrial diatoms at the catchment scale: Controls on communities. PeerJ 2020, 8, e8296. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  83. Şahin, B.; Barınova, S. Role of Altitude in Formation of Diatom Diversity of High Mountain Protected Glacier Lakes in the Kaçkar Mountains National Park, Rize, Turkey. Environments 2022, 9, 127. [Google Scholar] [CrossRef]
  84. Poradowska, A. Diatoms (Bacillariophyta) from the Genus Eunotia and Pinnularia developing on Soils in the Open Landscape of the Low Beskids. J. Ecol. Eng. 2020, 21, 257–270. [Google Scholar] [CrossRef]
  85. Soininen, J.; Paavola, R.; Muotka, T. Benthic diatom communities in boreal streams: Community structure in relation to environmental and spatial gradients. Ecography 2004, 27, 330–342. [Google Scholar] [CrossRef]
  86. Tuji, A.; Williams, D.M. Type Examination of Japanese Diatoms Described by Friedrich Meister (1913) from Lake Suwa. Natl. Sci. Mus. 2007, 33, 69–79. [Google Scholar]
  87. Lobo, E.A.; Schuch, M.; Heinrich, C.G.; Ben, A.; Düpont, A.; Wetzel, C.E.; Ector, L. Development of the Trophic Water Quality Index (TWQI) for subtropical temperate Brazilian lotic systems. Environ. Monit Assess 2015, 187, 354. [Google Scholar] [CrossRef]
  88. Levkov, Z.; Mitic-Kopanja, D.; Reichardt, E. The Diatom Genus Gomphonema from the Republic og Macedonia; Lange-Bertalot, H., Ed.; Diatoms of Europe; Koeltz Botanical Books: Oberreifenberg, Germany, 2016; ISBN 9783946583004/3946583008. [Google Scholar]
  89. Domingues, R.B.; Guerra, C.C.; Barbosa, A.B.; Galvão, H.M. Are nutrients and light limiting summer phytoplankton in a temperate coastal lagoon? Aquat. Ecol. 2015, 49, 127–146. [Google Scholar] [CrossRef]
  90. Soininen, J.; Jamoneau, A.; Rosebery, J.; Passy, S.I. Global patterns and drivers of species and trait composition in diatoms. Glob. Ecol. Biogeogr. 2016, 25, 940–950. [Google Scholar] [CrossRef]
  91. Lange, K.; Liess, A.; Piggott, J.J.; Townsend, C.R.; Matthaei, C.D. Light, nutrients and grazing interact to determine stream diatom community composition and functional group structure. Freshw. Biol. 2011, 56, 264–278. [Google Scholar] [CrossRef]
  92. Tapolczai, K.; Bouchez, A.; Stenger-Kovács, C.; Padisák, J.; Rimet, F. Taxonomy- or trait-based ecological assessment for tropical rivers? Case study on benthic diatoms in Mayotte island (France, Indian Ocean). Sci. Total Environ. 2017, 607–608, 1293–1303. [Google Scholar] [CrossRef]
  93. Stenger-Kovács, C.; Lengyel, E.; Buczkó, K.; Padisák, J.; Korponai, J. Trait-based diatom functional diversity as an appropriate tool for understanding the effects of environmental changes in soda pans. Ecol. Evol. 2020, 10, 320–335. [Google Scholar] [CrossRef] [Green Version]
  94. Jüttner, I.; Ector, L.; Reichardt, E.; Van De Vijver, B.; Jarlman, A.; Krokowski, J.; Cox, E.J. Gomphonema varioreduncum sp. Nov., a new species from northern and western Europe and a re-examination of Gomphonema exilissimum. Diatom Res. 2013, 28, 303–316. [Google Scholar] [CrossRef]
  95. Faustino, S.B.; Fontana, L.; Cristina, E.; Bartozek, R.; Eduardo, C.; Bicudo, D.M.; Bicudo, D.D.C. Composition and distribution of diatom assemblages from core and surface sediments of a water supply reservoir in Southeastern Brazil. Biota Neotrop. 2016, 16, 1–23. [Google Scholar] [CrossRef]
  96. Marra, R.C.; Tremarin, P.I.; Algarte, V.M.; Ludwig, T.A.V. Epiphytic diatoms (Diatomeae) from Piraquara II urban reservoir, Paraná state. Biota Neotrop. 2016, 16, 1–19. [Google Scholar] [CrossRef] [Green Version]
  97. Pillsbury, R.; Stevenson, R.J.; Munn, M.D.; Waite, I. Relationships between diatom metrics based on species nutrient traits and agricultural land use. Environ. Monit. Assess. 2019, 191, 228. [Google Scholar] [CrossRef]
  98. Nystrand, M.I.; Österholm, P.; Nyberg, M.E.; Gustafsson, J.P. Metal speciation in rivers affected by enhanced soil erosion and acidity. Appl. Geochem. 2012, 27, 906–916. [Google Scholar] [CrossRef]
  99. Martins, W.A.; Martins, L.L.; de Maria, I.C.; de Moraes, J.F.L.; Pedro Júnior, M.J. Reduction of sediment yield by riparian vegetation recovery at distinct levels of soil erosion in a tropical watershed. Ciência Agrotecnologia 2021, 45, e028220. [Google Scholar] [CrossRef]
  100. Falasco, E.; Ector, L.; Wetzel, C.E.; Badino, G.; Bona, F. Looking back, looking forward: A review of the new literature on diatom teratological forms (2010–2020). Hydrobiologia 2021, 848, 1675–1753. [Google Scholar] [CrossRef]
  101. Morin, S.; Cordonier, A.; Lavoie, I.; Arini, A.; Blanco, S.; Duong, T.T.; Torne’s, E.; Bonet, B.; Corcoll, N.; Faggiano, L.; et al. Consistency in Diatom Response. In Emerging and Priority Pollutants in Rivers; Lavoie, I., Fortin, C., Arini, A., Blanco, S., Becares, E., Duong, T.T., Tornés, E., Torne’s, E., Sabater, S., Bonet, B., et al., Eds.; Springer: Berlin/Heidelberg, Germany, 2012; pp. 118–146. ISBN 9783642257223. [Google Scholar]
  102. Pandey, L.K.; Lavoie, I.; Morin, S.; Park, J.; Lyu, J.; Choi, S.; Lee, H.; Han, T. River water quality assessment based on a multi-descriptor approach including chemistry, diatom assemblage structure, and non-taxonomical diatom metrics. Ecol. Indic. 2018, 84, 140–151. [Google Scholar] [CrossRef] [Green Version]
  103. Dos Santos, H.G.; Bhering, S.B.; Bognola, I.A.; Curcio, G.R.; de Carvalho Junior, W.; Chagas, C.D.S.; Manzatto, C.V.; Aglio, M.L.D.; Silva, J.D.S.; Chaffin, C.E. Mapa de Solos do Estado do Paraná. Embrapa Solos 2007, 96, 73. [Google Scholar]
  104. Costa, L.F.; Wetzel, C.E.; Lange-Bertalot, H.; Ector, L.; Bicudo, D.C. Taxonomy and Ecology of Eunotia Species (Bacillariophyta) in Southeastern Brazilian Reservoirs; Cramer: Stuttgart, Germany, 2017; ISBN 9783443570552. [Google Scholar]
  105. Amaral, M.W.W.; Medeiros, G.; Daufenbach, V.; Suszek Gonçalves, M.; Alvim Veiga Ludwig, T.; Catarina Bueno, N. Nupela semifasciata (Bacillariophyceae), a new species from subtropical lotic environments in Western Paraná State, Brazil. Fottea 2021, 21, 152–163. [Google Scholar] [CrossRef]
  106. Ortiz-Colón, A.I.; Piñero-Santiago, L.E.; Rivera, N.M.; Sosa, M.A. HHS Public Access. J. Env. Anal Toxicol 2016, 6, 1000353. [Google Scholar] [CrossRef]
  107. Tao, Y.; Yuan, Z.; Xiaona, H.; Wei, M. Ecotoxicology and Environmental Safety Distribution and bioaccumulation of heavy metals in aquatic organisms of different trophic levels and potential health risk assessment from. Ecotoxicol. Environ. Saf. 2012, 81, 55–64. [Google Scholar] [CrossRef]
  108. Virtanen, L.; Soininen, J. The roles of environment and space in shaping stream diatom communities. Eur. J. Phycol. 2012, 47, 160–168. [Google Scholar] [CrossRef]
  109. Pajunen, V.; Luoto, M.; Soininen, J. Climate is an important driver for stream diatom distributions. Glob. Ecol. Biogeogr. 2016, 25, 198–206. [Google Scholar] [CrossRef]
  110. Zorzal-Almeida, S.; Salim, A.; Andrade, M.R.M.; de Nascimento, M.N.; Bini, L.M.; Bicudo, D.C. Effects of land use and spatial processes in water and surface sediment of tropical reservoirs at local and regional scales. Sci. Total Environ. 2018, 644, 237–246. [Google Scholar] [CrossRef] [PubMed]
  111. Gauthier, M.; Launay, B.; Le Goff, G.; Pella, H.; Douady, C.J.; Datry, T. Fragmentation promotes the role of dispersal in determining 10 intermittent headwater stream metacommunities. Freshw. Biol. 2020, 65, 2169–2185. [Google Scholar] [CrossRef]
  112. Cook, J.R. Decompressive Procedures. Vet. Clin. N. Am. Small Anim. Pract. 1992, 22, 917–921. [Google Scholar] [CrossRef] [PubMed]
  113. Vilmi, A.; Karjalainen, S.M.; Heino, J. Ecological uniqueness of stream and lake diatom communities shows different macroecological patterns. Divers. Distrib. 2017, 23, 1042–1053. [Google Scholar] [CrossRef] [Green Version]
  114. Wetzel, C.E.; de Bicudo, D.C.; Ector, L.; Lobo, E.A.; Soininen, J.; Landeiro, V.L.; Bini, L.M. Distance Decay of Similarity in Neotropical Diatom Communities. PLoS ONE 2012, 7, e45071. [Google Scholar] [CrossRef] [PubMed]
  115. Bried, J.T.; Vilmi, A. Improved detection of mass effect species assembly for applied metacommunity thinking. J. Appl. Ecol. 2022, 59, 921–926. [Google Scholar] [CrossRef]
  116. Record, S.; Voelker, N.M.; Zarnetske, P.L.; Wisnoski, N.I.; Tonkin, J.D.; Swan, C.; Marazzi, L.; Lany, N.; Lamy, T.; Compagnoni, A.; et al. Novel Insights to Be Gained From Applying Metacommunity Theory to Long-Term, Spatially Replicated Biodiversity Data. Front. Ecol. Evol. 2021, 8, 612794. [Google Scholar] [CrossRef]
Figure 1. Location of sampling stations (S) in Cascavel River micro-watershed, Cascavel City—PR, Brazil.
Figure 1. Location of sampling stations (S) in Cascavel River micro-watershed, Cascavel City—PR, Brazil.
Water 14 03913 g001
Figure 2. Area of influence containing detailed land use, and sampling stations following the upstream-downstream flow, with a diameter of 1 km.
Figure 2. Area of influence containing detailed land use, and sampling stations following the upstream-downstream flow, with a diameter of 1 km.
Water 14 03913 g002
Figure 3. Daily precipitation (mm) in Cascavel city—PR throughout the year 2018. Data provided by SIMEPAR.
Figure 3. Daily precipitation (mm) in Cascavel city—PR throughout the year 2018. Data provided by SIMEPAR.
Water 14 03913 g003
Figure 4. Principal Component Analysis (PCA) for the data evaluated at the sampling stations (S) of the Cascavel River, Cascavel city—PR, Brazil. Limnological data (Ec: electrical conductivity, DO: dissolved oxygen, Turb: turbidity, Fl: flow, COD: chemical oxidation demand, BOD: biochemical oxygen demand, NO3: nitrate, TP: total dissolved phosphorus, CLa: chlorophyll a, TS: total solids, CT: total coliforms), metals measured in the surface sediment (Cu: copper; Al: aluminum; Ni: nickel; Zn: zinc; Na: sodium; Mn: manganese; Pb: lead; Cr: chrome).
Figure 4. Principal Component Analysis (PCA) for the data evaluated at the sampling stations (S) of the Cascavel River, Cascavel city—PR, Brazil. Limnological data (Ec: electrical conductivity, DO: dissolved oxygen, Turb: turbidity, Fl: flow, COD: chemical oxidation demand, BOD: biochemical oxygen demand, NO3: nitrate, TP: total dissolved phosphorus, CLa: chlorophyll a, TS: total solids, CT: total coliforms), metals measured in the surface sediment (Cu: copper; Al: aluminum; Ni: nickel; Zn: zinc; Na: sodium; Mn: manganese; Pb: lead; Cr: chrome).
Water 14 03913 g004
Figure 5. Results of the Partial Redundancy Analysis (pRDA). Relative contributions (% explanation) of environmental (E), spatial (S), historical community (HC) and landscape (L) that explain the variation in relative abundance (A) and Occurrence (O) of periphytic diatom metacommunities. We used possible dispersal routes “overland”, “AEM”, and “Watercourse” as spatial predictors (see methods). U = unexplained component. Values < 0.5% were not shown. Explanation percentages significantly different from a null model (p < 0.05) are represented by “*”.
Figure 5. Results of the Partial Redundancy Analysis (pRDA). Relative contributions (% explanation) of environmental (E), spatial (S), historical community (HC) and landscape (L) that explain the variation in relative abundance (A) and Occurrence (O) of periphytic diatom metacommunities. We used possible dispersal routes “overland”, “AEM”, and “Watercourse” as spatial predictors (see methods). U = unexplained component. Values < 0.5% were not shown. Explanation percentages significantly different from a null model (p < 0.05) are represented by “*”.
Water 14 03913 g005
Table 1. Location with geographic coordinates of the collection points along the Cascavel—PR river, Brazil.
Table 1. Location with geographic coordinates of the collection points along the Cascavel—PR river, Brazil.
Sampling SitesGeographic Coordinates
LatitudeLongitude
S124°57′32.93″ S53°26′33.32″ W
S224°58′35.81″ S53°26′6.85″ W
S325°0′13.36″ S53°26′19.12″ W
S424°58′17.29″ S53°24′24.75″ W
S524°59′5.99″ S53°25′23.83″ W
S624°59′18.90″ S53°25′32.20″ W
S724°59′48.82″ S53°26′42.08″ W
S825°0′13.36″ S53°26′19.12″ W
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Medeiros, G.; Padial, A.A.; Amaral, M.W.W.; Guicho, R.; Pilatti, M.C.; Sampaio, S.C.; Ludwig, T.A.V.; Bueno, N.C.; dos Reis, R.R. Exploring Key Determinants of the Periphytic Diatom Community in a Southern Brazilian Micro-Watershed. Water 2022, 14, 3913. https://doi.org/10.3390/w14233913

AMA Style

Medeiros G, Padial AA, Amaral MWW, Guicho R, Pilatti MC, Sampaio SC, Ludwig TAV, Bueno NC, dos Reis RR. Exploring Key Determinants of the Periphytic Diatom Community in a Southern Brazilian Micro-Watershed. Water. 2022; 14(23):3913. https://doi.org/10.3390/w14233913

Chicago/Turabian Style

Medeiros, Gabriela, André Andrian Padial, Mailor Wellinton Wedig Amaral, Ricardo Guicho, Maria Clara Pilatti, Silvio Cesar Sampaio, Thelma Alvim Veiga Ludwig, Norma Catarina Bueno, and Ralpho Rinaldo dos Reis. 2022. "Exploring Key Determinants of the Periphytic Diatom Community in a Southern Brazilian Micro-Watershed" Water 14, no. 23: 3913. https://doi.org/10.3390/w14233913

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