Printer Friendly

Periphytic diatom as bioindicators in urban and rural streams/Diatomaceas perifiticas como bioindicadores em corregos urbanos e rurais.


Over the last century, humans have caused drastic environmental changes through deforestation and the conversion of land for agricultural use (WINEMILLER et al., 2008). Urbanization stands out as a major modification to an ecosystem, by reducing soil permeability and causing hydrographic alterations, increasing nutrient concentrations and changing the morphology of waterways (WALSH et al., 2005).

With the increased human influence on hydrographic basins, the identification of factors that affect aquatic ecosystems and how these factors influence aquatic communities is a great challenge (PAN et al., 2004). The structure of aquatic communities is determined by processes that operate at multiple spatial scales, such as differences in habitat, current flow, nutrient availability and luminosity, which directly affect the organisms (TISON et al., 2005).

Biological communities in rivers and streams are important components in the evaluation of water quality (WHITTON; KELLY, 1995). Diatoms are important autotrophs in streams and rivers (STOERMER; SMOLL, 1999), and are sensitive to the physical and chemical variations of a body of water (WINTER; DUTHIE, 2000a and b). The composition and relative abundance of diatoms are determined by the preferences and tolerance of the species (LANGE-BERTALOT, 1979; VAN DAM et al., 1994; POTAPOVA; CHARLES, 2003). Variations in these attributes occurs in both the spatial and temporal scales and are ruled by differences in physical and chemical properties of water (STEVENSON; PAN, 1999), which has led to the growing use of diatoms in studies that monitor the water quality of rivers and streams (STOERMER; SMOLL, 1999; SOININEN et al., 2004). Variations may also be related to human activity in water basins, such as urbanization and agriculture (FORE; GRAFE, 2002).

Although several indices have been developed to evaluate the water quality using diatoms (VAN DAM et al., 1994; STEVENSON; PAN, 1999), only few studies have been conducted to assess the responses from diatoms along basins influenced by urbanization and by agricultural practices (LOBO et al., 1996; GOMEZ; LICURSI, 2001; JUTTNER et al., 2003), especially in tropical regions. In Brazil, bioindication studies of water quality using diatoms were developed from the 80's, mainly in South and Southeastern regions, focusing on lotic environments (LOBO; TORGAN, 1988; LOBO et al., 1996, 2002, 2004a, b and c; LOBO; CALLEGARO, 2000; BURLIGA et al., 2004; HERMANY et al., 2006; SALOMONI et al., 2006; DUPONT et al., 2007, among others).

The objectives of this study were: i) to evaluate differences in the structure of periphytic diatom assemblages of streams from an urban and an rural watersheds and to correlate the differences in the assemblages with physical and chemical variables; ii) to analyze the temporal variability in the periphytic diatom assemblage in these streams; and iii) to select indicator species in each environment. Therefore, the following hypotheses were tested: a) the chemical influence of water on assemblages of periphytic diatoms in an urban environment is greater than on assemblages from a rural environment; b) the spatial variation of the assemblage of periphytic diatoms (urban and rural streams) is more important than the temporal variation (month of sampling).

Material and methods

Study area

The Pirapo river basin is located in the northeast corner of Parana State, within the upper Parana river system. The region has a drainage area of approximately 5,076 [km.sup.2], situated in the Serra Geral Formation, which is composed of basic igneous rocks, such as basalt (MINEROPAR, 2006), and a predominance of red nitosols. Waters from this basin are used for supply, development of agricultural activities and ecological tourism for most of the cities in this region. One city that receives the water from this basin is Maringa that is on the water divide of the Pirapo river basin and the Ivai river basin, and the headwaters of several streams are located within the urban perimeter of Maringa.

The Nazare Stream micro-basin has a drainage area of 867.928 hectares, is located in the urban area of Maringa (Figure 1) and contains 34.8% impervious surfaces, which are mostly located in residential and industrial areas. In this microbasin, the main industrial activities are related to metallurgy, plastic and petroleum products (KUHL et al., 2010). The banks are steep throughout its length. At the mouth and at the intermediate region, riparian forest is present on the right bank. At the headwater, the riparian forest is more developed and is present on both banks.

The Remo Stream micro-basin has a drainage area of 792.325 hectares, is located in the rural area of the municipality of Maringa (Figure 1) and has 0.5% impervious surfaces. Farming practices are the main activities performed in this microbasin (crop rotation: maize, soybean and wheat). The Remo Stream has riparian vegetation on both banks along its course.

Abiotic variables

The sampling to determine the abiotic variables was performed concurrently with the sampling of biotic variables. Data of physical and chemical conditions of water, such as pH (DIGIMED DM2), electrical conductivity (DIGIMED DM3, [micro]S [cm.sup.-1]), dissolved oxygen, water temperature (YSI 55/12FT, ml [L.sup.-1] and [degrees]C, respectively) and flow (FLO-MATE 2000, Marsh McBirdey, m [s.sup.-1]) were all measured in the field with portable analytical equipment. The total nitrogen concentration, orthophosphate, biochemical oxygen demand ([BOD.sub.5], mg O2 [L.sup.-1]), chemical oxygen demand (COD, mg O2 [L.sup.-1]) and oil (mg [L.sup.-1]) in water samples were analyzed by technicians from the Sanitation and Agrochemistry Laboratories from the State University of Maringa.

Sampling and collection of diatoms

Six samples of periphytic diatoms were collected from each stream every two months from July 2007 to June 2008. Each sample consisted of three pebbles (composite sample). This substrate was selected because it was the most abundant substrate and was present in both streams. The side of the rock opposite to the direction of the water current was scraped with the aid of a brush and a blade, and the material was fixed with 4% formalin (1:1 ratio). Measurements of the shaved surface area were taken with a caliper. The material was oxidized and Hyrax was used to mount the permanent slides. Slides with control material were deposited in the Herbarium of the State University of Maringa, in Maringa, Parana State (HUEM).

Identification and counting of the diatom species was performed using a light microscope (Olympus CX31). Individuals were identified and counted until recording a minimum of 600 valves, as recommended by Kobayasi and Mayama (1982), and a 90% counting efficiency was used in accordance with Pappas and Stoermer (1996). The concentration of cells.[cm.sup.-2] was estimated by multiplying the number of valves from each taxon by a conversion factor according to the methods of Hermany et al. (2006).

To identify the periphytic diatom taxa, the following taxonomic texts were used: Krammer and Lange-Bertalot (1991a and b), Rumrich et al. (2000), Lange-Bertalot (2001), Krammer (2002), Metzeltin and Lange-Bertalot (1998, 2007), Metzeltin et al. (2005) and Tremarin et al. (2009, 2010).

Data analysis

To identify possible significant spatial and temporal differences between the means of the attributes (richness, evenness, Shannon-Wiener diversity index and density) from the periphytic diatom assemblages, analysis of variance (ANOVA) was used, and the temporal variability was controlled through blocks (ANOVA; period factor: block; Local factor: rural and urban). Assumptions of normality and homoscedasticity were checked using Shapiro-Wilk and Levene tests, respectively. When the assumptions were not met, data were log transformed. When the ANOVA was significant, a Tukey's test was employed to determine the level that has varied.

In order to summarize the assemblage structure of periphytic algae, a non-metric multidimensional scaling (NMS) was used (KRUSKAL, 1964). The Sorensen distance was calculated, and the general procedure for NMS was followed according to previously published methods (MCCUNE; GRACE, 2002). One-hundred permutations were performed, and the standard deviation was used as the stability criterion ([less than or equal to] 0.005, stress above 100 interactions). This analysis was run with the matrix of abundance data (log transformed to remove the effect of high values) in different sampling locations and periods.

For testing significant differences between each location and period (summarized by NMS), this study used a multi-response permutation procedure (MRPP), which is a non-parametric method used to test for multivariate differences between predefined groups (ZIMMERMAN et al., 1985). The significance of the null hypothesis that the locations and periods were not different was tested by Monte Carlo randomization (based on 10,000 permutations).

To determine indicator species (IndVal), the methods presented by Dufrene and Legendre (1997) were followed, by which the abundance and frequency of occurrence of the species in each group was used as the input data and indicator values calculated for each species (MCCUNE; GRACE, 2002). The species with p < 0.05 (based on 10,000 permutations) from the Monte Carlo test were considered to be an indicator species.

Environmental variables were summarized by a principal component analysis (PCA). To determine the components that should be retained for interpretation, the Broken Stick criterion was employed. According to this model, only axes with an eigenvalue greater than the eigenvalues generated by randomization should be interpreted (MCCUNE; GRACE, 2002). The abiotic data, except pH, were log transformed for PCA.

Relationships between the multivariate analyses (environmental variables and structure of periphytic diatom assemblage) were examined by a Procrustes analysis (PERES-NETO; JACKSON, 2001). In this analysis, the two matrices are compared using a logarithm that minimizes the residual sum of squares between two matrices (ROHLF; SLICE, 1990). The resulting value of [m.sup.-2] is the best fit, which assumes that it describes the degree of association between the matrices.

NMS, MRPP, IndVal and PCA analyses were performed using the program PC-Ord[R] 4.0 (MCCUNE; MEFFORD, 1999). Procrustes statistics were calculated with the PROTEST[R] program (JACKSON, 1995). The ANOVA was computed using the software Statistica[TM] 7.0. The level of significance was p < 0.05.


Table 1 shows the variation of physical and chemical variables of the studied streams. Water temperature and pH had no expressive differences between the urban and rural streams. Conductivity and total nitrogen presented higher values in the urban stream. Considering the orthophosphate, higher concentrations were observed in the urban stream, except in periods 1 and 2. In periods 3 and 6, values of this variable were similar between streams. Concentrations of dissolved organic carbon were greater in the rural stream, except in periods 4 and 5. Values of biochemical oxygen demand were also higher in the rural stream, with exception of periods 4 and 6. For oxygen and flow, values were slightly higher in the rural stream.

PCA has determined that 53.3% total variability of the abiotic data was in the first two axes (Table 2). In PCA 1 (30.63%), there was a clear separation between the streams. There were positive correlations with dissolved oxygen, chemical oxygen demand, biochemical oxygen demand and flow. A negative correlation was detected with conductivity and total nitrogen. In PCA 2 (22.67%), there were fluctuations in abiotic factors along the sampled periods, especially at the urban stream (Table 2 and Figure 2). The positively correlated variables were pH and dissolved oxygen, while water temperature, orthophosphate, chemical oxygen demand and biochemical oxygen demand showed a negative correlation.

In total, 135 species were identified, and 69 taxa were common to both streams. One hundred and twenty four were detected in the rural stream, being 55 exclusive in this stream. Eighty taxa were observed in the urban stream and 11 species were exclusively found in this stream. Significant differences were detected for species richness between the streams (ANOVA, F = 28.92; p = 0.000009; Figure 3a), however, there was no significant difference in the time scale for this attribute (F = 0.4523; p = 0.80; Figure 3b).

The higher mean evenness was recorded in the rural stream (ANOVA, F = 5.274; p = 0.02; Figure 4a). However, no significant difference was found in evenness on the time scale (F = 0.519; p = 0.75; Figure 4b).

The Shannon-Wiener diversity index had higher values in the rural stream (ANOVA, F = 12.38; p = 0.001; Figure 5a); however, there was no significant difference in diversity between the sampled periods (F = 0.38; p = 0.85; Figure 5b).

The highest mean values of density were observed in the urban stream (ANOVA, F = 13.002; p = 0.001; Figure 6a); but no significant difference in density was registered on the time scale (F = 1.28; p = 0.29; Figure 6b).

The structure of periphytic diatom assemblages summarized by a NMS showed a difference between locations and periods (Figure 7). After 31 interactions, the stability criterion was achieved with a final stress of 13.73 (Monte Carlo test, p = 0.009), and three axes were retained for interpretation. The proportion of variance represented by each axis, which was based on the distance between r2 within the ordination space and distances in the original space, was 0.512 for axis 1, 0.158 for axis 2 and 0.218 for axis 3, all totaling 0.888.

The spatial scale was identified as the main pattern of the structure of the periphytic diatom assemblages, which was determined by plotting axes 1 and 2 (Figure 7a), axes 1 and 3 (Figure 7b) and axes 2 and 3 (Figure 7c), corroborating the result of the MRPP (p = 0.00000). The formation of a group for the rural stream and another group for the urban stream was related to the increased conductivity and total nitrogen values recorded for the urban stream, and to higher dissolved oxygen values, COD, [BOD.sub.5] and flow in the rural stream. These results are confirmed by the matrix correlation test (Procrustes) between the first three axes of PCA ordination, and the first three axes of NMS ordination. The adjusted value for the distribution of periphytic diatoms was [m.sup.-2] = 0.7607 and p = 0.0001, which supported statistically the influence of abiotic variables on the spatial and temporal distribution of the assemblages in the urban and rural streams.

A clear temporal variability was not observed for attributes analyzed separately (Figures 3b, 4b, 5b and 6b). But when considered the community structure (richness and density simultaneously) a significant difference was detected (Figures 7a, b and c--MRPP, p = 0.0006).

The analysis of indicator species (IndVal) has identified diatom species strongly associated with urban and rural streams (IndVal, p [less than or equal to] 0.05; Table 3). Thirty-two taxa were indicators in the rural stream, and 17 taxa were indicators in the urban stream.


The assemblages of the periphytic diatoms studied were different between the urban stream and the rural stream. Analyses performed using attributes from this assemblage, such as richness, evenness, diversity and total density, indicated that the spatial variation (urban and rural streams) was more significant than temporal variation (month of sampling).

Studies on streams that drain urban areas have shown a decline in the richness of diatom species, and the reduction in species richness is associated with organic pollution (SONNEMAN et al., 2001). Common effects from pollution are the reduction of species diversity, increased density and an increase of tolerant species (JUTTNER et al., 2003; NDIRITU et al., 2006). In addition, under conditions of intermediate nutrient concentrations, the periphytic diatoms may have increased diversity (JUTTNER et al., 2003). In this study, higher concentrations of nutrients in the urban stream have affected periphytic diatoms. This stream presented a high density, but attributes such as richness, evenness and diversity were reduced, owing nutrient enrichment and human disturbance. In contrast, richness had higher values in the rural stream when compared to those recorded for the urban stream.

Studies that involve periphytic diatom assemblages have reported the influence of urban and rural environments on the structure and distribution of these assemblages. Winter and Duthie (1998) observed that the differentiation between streams basin were related to changes in temperature, [BOD.sub.5], total phosphorus and suspended solids. Winter and Duthie (2000a) showed that the differentiation of periphytic diatom assemblages along an urban-rural gradient was related to differences in the concentrations of total phosphorus and nitrogen. Ndiritu et al. (2006) found differences in the assemblages of periphytic diatoms in areas with smaller towns and subsistence-farming activities when compared to assemblages found in large urban centers and areas with intensive agriculture, which was associated with changes in the composition of species and was related to the intensity of pollution. In this study, the difference in periphytic diatom assemblages between urban and rural streams were related to the increased conductivity and total nitrogen values recorded for the urban stream, the higher dissolved oxygen values, COD, [BOD.sub.5] and flow in the rural stream (Procrustes--[m.sup.-2] = 0.7607 and p = 0.0001). In urban landscapes, impervious surface areas intensify the surface runoff, which carries nutrients (WALSH et al., 2005; PAUL; MEYER, 2008) and contaminants to the streams. Deficiencies in waste treatment systems, in addition to illicit discharges of effluents into urban streams, cause an increase in the concentration of nutrients, especially nitrogen (PAUL; MEYER, 2008). High nitrogen concentrations contribute to the excessive development of algae (PORTER et al., 2008), which may lead to a reduction in the concentration of dissolved oxygen (KANNEL et al., 2007). These conditions were observed in the urban stream, which presented severe conditions for the maintenance of species sensitive to increased concentrations of nitrogen and conductivity.

Increased conductivity values were also considered to be one of the main effects of urbanization to streams (SONNEMAN et al., 2001; WALKER; PAN, 2006), and this result can explain great part of the variability between periphytic diatom assemblages (POTAPOVA; CHARLES, 2003). The influence of this variable on periphytic diatoms has been previously reported in several studies (SONNEMAN et al., 2001; SOININEN et al., 2004; WALKER; PAN, 2006; WINEMILLER et al., 2008).

Factors that are frequently indicated as important for the distribution and abundance of periphytic algae, such as phosphorous (WINTER; DUTHIE, 2000a), pH (PAN et al., 1996) and temperature (DENICOLA, 1996), had no influence on the separation of periphytic diatom assemblages between the streams studied; however, they did contribute to the temporal distribution. As observed in PCA, the temporal variation of these abiotic factors did not show a similar pattern between the urban and rural streams because they reflected the activities developed in their micro-basins, and these micro-basins do not show differences in soil or climate. Changes in the structure of periphytic diatom assemblages were expected in response to seasonal changes in abiotic variables (STEVENSON; PAN, 1999). Nevertheless, changes in physical and chemical variables caused by human activities in the micro-basins of the studied streams were more important for the structuring of the assemblages of the periphytic diatoms than the seasonal variation of abiotic factors. Although both agriculture and urbanization affect rivers and streams, pollution from urban areas is more intense than the pollution detected in rural areas (KANNEL et al., 2007).

Species that had a greater indicator potential for the urban stream have already been identified by other authors as tolerant to nutrient enrichment, such as Fallacia monoculata (POTAPOVA; CHARLES, 2007) and Mayamaea atomus var. permitis (FORE; GRAFE, 2002). Species that are frequently associated with conditions of organic pollution were also indicators for the urban stream, such as Achnanthes lanceolata (VAN DAM et al., 1994), Amphora montana (LOBO et al., 2002; VAN DAM et al., 1994), Cyclotella meneghiniana (LOBO et al., 2002; POTAPOVA; CHARLES, 2007), Eolimna minima (RIMET, 2009), Navicula tenelloides (VAN DAM et al., 1994; SOININEN et al., 2004), Nitzschia cf. inconspicua (VAN DAM et al., 1994), Pinnularia gibba (LOBO et al., 2002) and Sellaphora seminulum (LOBO et al., 2002; SALOMONI et al., 2006; RIMET, 2009). Caloneis bacillum was reported by Van Dam et al. (1994) as being typical of mesotrophic environments, and this species was an indicator for the urban environment in the streams studied. Gomphonema lagenula was an indicator for the urban stream in this study and was an indicator of nutritionally poor environments in a study by Chessman et al. (2006), and our results suggest a greater range of tolerance to nutrient concentrations for this species. Achnanthidium minutissimum tolerates conditions that range from oligotrophic to eutrophic (VAN DAM et al., 1994; POTAPOVA; CHARLES, 2007). Here, this species was recorded as an indicator of the urban stream, which was more affected when compared to the rural stream.

Among the indicator species for the rural stream, there are taxa that have already been observed to be tolerant to conditions that vary from oligotrophic to eutrophic environments, such as Achnanthes exigua, Hantzschia amphioxys, Navicula cryptotenella (VAN DAM et al., 1994), Ulnaria ulna (VAN DAM et al., 1994; POTAPOVA; CHARLES, 2007) and Achnanthes rupestoides, which was an indicator of oligotrophic environments in the study of Van Dam et al. (1994) and of meso-eutrophic environments in Hermany et al. (2006). Amphipleura lindheimeri was recorded by Lobo et al. (2002) as an indicator of heavily polluted environments, and Lobo et al. (2004b) showed that this species has an intermediate tolerance to pollution. Sellaphora pupula preferred less eutrophic environments in the study of Hermany et al. (2006) and polluted to very polluted rivers in Lobo et al. (2002) and in Salomoni et al. (2006). There are other taxa mentioned in the literature that are characteristic of eutrophic environments, for example, Amphora copulata, Tryblionella levidensis (VAN DAM et al., 1994) and Frustulia crassinervea (SALOMONI et al., 2006), which were indicators of the rural stream in this study. Frustulia vulgaris was determined to be sensitive to high concentrations of nutrients (KELLY et al., 1998), corroborating the results for the two streams studied herein. Meanwhile, species that showed preference for less severe conditions of this type of pollution were Fallacia ecuadoriana, Fallacia insociabilis, Geissleria neosubtropica, Luticola dapalis, Navicula lohmannii, Pinnularia sp., Placoneis constans var. symmetrica, Placoneis disparilis, Placoneis hambergii, Nupela praecipua, Sellaphora sp. 1, Sellaphora sp. 2, Stauroneis cf. kriegeri and Stenopterobia schweickerdtii. No ecological information was found with regard to these species, which suggests that they may be indicators of less affected environments.


There was a clear effect of land-use on abiotic variables of the streams, which was in turn reflected by diatoms. The spatial variation of these assemblages was greater than the temporal variation, and the human influence was more pronounced in urban stream, which showed a reduced richness and increased density due to high values of total nitrogen and conductivity. These conditions may have contributed to the highest number of tolerant species in the urban environment, suggesting that the land-use had a greater influence on the assemblage structure of periphytic diatoms than the seasonal variation of abiotic factors.

Doi: 10.4025/actascibiolsci.v36i1.18175


To the CNPq for the financial support to the project "Identification of possible bioindicators in urban aquatic ecosystems: the response of the groups of organisms to stress gradients" and to CAPES for scholarship. To PEA/UEM and Nupelia/UEM for logistic support. To Priscila I. Tremarin and Thelma A.V. Ludwig for the help in identifying diatoms. To Eder A. Gubiani for assistance with statistical analyses.


BURLIGA, A. L. M.; SCHWARZBOLD, A.; LOBO, E. A.; PILLAR, V. D. Functional types in epilithon algae communities of the Maquine River, Rio Grande do Sul, Brazil. Acta Limnologica Brasiliensia, v. 16, n. 4, p. 369-380, 2004.

CHESSMAN, B. C.; FRYIRS, K. A.; BRIERLEY, G. J. Linking geomorphic character, behaviour and condition to fluvial biodiversity: implications for river management. Aquatic Conservation: Marine and Freshwater Ecosystems, v. 16, n. 3, p. 267-288, 2006.

DENICOLA, D. M. Periphyton responses to temperature at different ecological levels. In: STEVENSON, R. J.; BOTHWELL, M. L.; LOWE, R. L. (Ed.). Algal ecology: freshwater bentic ecosystems. San Diego: Academic Press, 1996. p. 149-181.

DUFRENE, M.; LEGENDRE, P. Species assemblages and indicator species: the need for a flexible symmetrical approach. Ecological Monographs, v. 67, n. 3, p. 345-366, 1997.

DUPONT, A.; LOBO, E. A.; COSTA, A. B.; SCHUCH, M. Avaliacao da qualidade da agua do Arroio do Couto, Santa Cruz do Sul, RS, Brasil. Caderno de Pesquisa. Serie Biologia, v. 19, n. 1, p. 20-31, 2007.

FORE, L. S.; GRAFE, C. Using diatoms to assess the biological condition of large rivers in Idaho (U.S.A.). Freshwater Biology, v. 47, n. 10, p. 2015-2037, 2002.

GOMEZ, N.; LICURSI, M. The Pampean Diatom Index (IDP) for assessment of rivers and streams in Argentina. Aquatic Ecology, v. 35, n. 2, p. 173-181, 2001.

HERMANY, G.; SCHWARZBOLD, A.; LOBO, E. A.; OLIVEIRA, M. A. Ecology of the epilithic diatom community in a low-order stream system of the Guaiba hydrographical region: subsidies to the environmental monitoring of southern Brazilian aquatic systems. Acta Limnologica Brasiliensia, v. 18, n. 1, p. 9-27, 2006.

JACKSON, D. A. Protest: a procrustean randomization test of community environment concordance. Ecoscience, v. 2, n. 3, p. 297-303, 1995.

JUTTNER, I.; SHARMA, S.; DAHAL, B. M.; ORMEROD, S. J.; CHIMONIDES, P. J.; COX, E. J. Diatoms as indicators of stream quality in the Kathmandu Valley and Middle Hills of Nepal and India. Freshwater Biology, v. 48, n. 11, p. 2065-2084, 2003.

KANNEL, P. R.; LEE, S.; KANEL, S. R.; LEE, Y. S. Spatial-temporal variation and comparative assessment of water qualities of urban river system: a case study of river Bagmati (Nepal). Environmental Monitoring and Assessment, v. 129, n. 1-3, p. 433-459, 2007.

KELLY, M. G.; CAZAUBON, A.; CORING, E.; DELL'UOMO, A.; ECTOR, L.; GOLDSMITH, B.; GUASH, H.; HURLIMANN, J.; JARLMAN, A.; KAWECKA, B.; KWANDRANS, J.; LAUGASTE, R.; LINDSTROM, E. A.; LEITAO, M.; MARVAN, P.; PADISAK, J.; PIPP, E.; PRYGIEL, J.; ROTT, E.; SABATER, S.; VAN DAM, H.; VIZINET, J. Recommendations for the routine sampling of diatoms for water quality assessments in Europe. Journal of Applied Phycology, v. 10, n. 2, p. 215-224, 1998.

KOBAYASI, H.; MAYAMA, S. Most pollution tolerant diatoms of severely polluted rivers in the vicinity of Tokyo. Japanese Journal of Phycology, v. 30, n. 3, p. 188-196, 1982.

KRAMMER, K.; LANGE-BERTALOT, H. Bacillariophyceae. Achnanthaceae. Kritische Erganzungen zu Navicula (Linolatae) und Gomphonema. SuBwasserflora von Mitteleuropa, v. 2, n. 4, p. 1-437, 1991a.

KRAMMER, K.; LANGE-BERTALOT, H. Bacillariophyceae: Centrales, Fragilariaceae, Eunotiaceae. SuBwasserflora von Mitteleuropa, v. 2, n. 3, p. 1-576, 1991b.

KRAMMER, K. The genus Cymbella. In: LANGEBERTALOT, H. (Ed.). Diatoms of Europe: diatoms of the european inland waters and comparable habitats. Ruggell: A. R. G. Gantner Verlag K. G., 2002. v. 3, p. 1-584.

KRUSKAL, J. B. Nonmetric multidimensional scaling: a numerical method. Psychometrika, v. 29, n. 2, p. 115-129, 1964.

KUHL, A. M.; ROCHA, C. L. M. S. C.; ESPINDOLA, E. L. G.; LANSAC-TOHA, F. A. Rural and urban streams: anthropogenic influences and impacts on water and sediment quality. International Review of Hydrobiology, v. 95, n. 3, p. 260-272, 2010.

LANGE-BERTALOT, H. Navicula sensu stricto--10 genera separated from Navicula sensu lato Frustulia. In:

LANGE-BERTALOT, H. (Ed.). Diatoms of Europe: diatoms of the european inland waters and comparable habitats. Ruggell: A. R. G. Gantner Verlag K. G., 2001. v. 2, p. 1-526.

LANGE-BERTALOT, H. Pollution tolerance of diatoms as a criterion for water quality estimation. Nova Hedwigia, v. 64, n. 1, p. 285-304, 1979.

LOBO, E. A.; TORGAN, L. C. Analise da estrutura da comunidade de diatomaceas (Bacillariophyceae) em duas estacoes do sistema Guaiba, Rio Grande do Sul, Brasil. Acta Botanica Brasilica, v. 1, n. 2, p. 103-109, 1988.

LOBO, E. A.; CALLEGARO, V. L. M. Avaliacao da qualidade de aguas doces continentais com base em algas diatomaceas epiliticas: Enfoque metodologico. In:

TUCCI, C. E. M.; MOTTA-MARQUES, D. M. (Org.). Avaliacao e controle da drenagem urbana. Porto Alegre: Ed. Universidade/UFRGS, 2000. p. 277-300.

LOBO, E. A.; CALLEGARO, V. L. M.; BENDER, E. P. Utilizacao de algas diatomaceas epiliticas como indicadoras da qualidade da agua em rios e arroios da regiao hidrografica do Guaiba, RS, Brasil. Santa Cruz do Sul: Edunisc, 2002.

LOBO, E. A.; BES, D.; TUSDEQUE, L.; ECTOR, L. Water quality assessment of the Pardinho river, RS, Brazil, using epilithic diatom assemblages and faecal coliforms as biological indicators. Vie Milieu, v. 54, n. 2-3, p. 115-125, 2004a.

LOBO, E. A.; CALLEGARO, V. L. M.; HERMANY, G.; BES, D.; WETZEL, C. E.; OLIVEIRA, M. A. Use of epilithic diatoms as bioindicators from lotic systems in southern Brazil, with special emphasis on eutrophication. Acta Limnologica Brasiliensia, v. 16, n. 1, p. 25-40, 2004b.

LOBO, E. A.; CALLEGARO, V. L. M.; WETZEL, C. E.; HERMANY, G.; BES, D. Water quality evaluation of Condor and Capivara Streams, Municipal District of Porto Alegre, RS, Brazil, using epilithic diatom communities as bioindicators. Oceanological and Hydrobiological Studies, v. 33, n. 1, p. 77-93, 2004c.

LOBO, E. A.; CALEGARO, V. L. M.; OLIVEIRA, M. A.; SALOMONI, S. F.; SCHULER, N. S.; ASAI, K. Pollution tolerant diatoms from lotic systems in the Jacui Basin, Rio Grande do Sul, Brazil. Iheringia. Serie Botanica, v. 47, n. 1, p. 45-72, 1996.

MCCUNE, B.; GRACE, J. B. Analysis of ecological communities. Oregon: Gleneden Beach, 2002.

MCCUNE, B.; MEFFORD, M. J. Pc-ord: multivariate analysis of ecological data (Version 4). Oregon: MJM Software Design, 1999.

METZELTIN, D.; LANGE-BERTALOT, H. Tropical diatoms of South America I. About 700 predominantly rarely know or new taxa representative of the neotropical flora. In: LANGE-BERTALOT, H. (Ed.). Iconographia diatomologica: annotated diatom micrographs. Konigstein: Koeltz Scientific Books, 1998. v. 5, p. 1-695.

METZELTIN, D.; LANGE-BERTALOT, H. Tropical diatoms of South America II, Special remarks on biogeographic disjunction. In: LANGE-BERTALOT, H. (Ed.). Iconographia diatomologica: annotated diatom micrographs. Konigstein: A.R.G. Gantner Verlag K.G., 2007. v. 18, p. 1-875.

METZELTIN, D.; LANGE-BERTALOT, H.; GARCIA-RODRIGUEZ, F. Diatoms of Uruguay compared with other taxa from South America and elsewhere. In: LANGEBERTALOT, H. (Ed.). Iconographia diatomologica: Annotated Diatom Micrographs. Konigstein: A. R. G. Gantner Verlag K. G., 2005. v. 15, p. 1-736.

MINEROPAR. Mapa Geologico do Estado do Parana, 1:650.000. Curitiba: Mineropar, 2006.

NDIRITU, G. G.; GICHUKI, N. N.; TRIEST, N. Distribution of epilithic diatoms in response to environmental conditions in an urban tropical stream, Central Kenya. Biodiversity and Conservation, v. 15, n. 10, p. 3267-3293, 2006.

PAN, Y.; HERLIHY, A.; KAUFMANN, P.; WIGINGTON, J.; VAN SICKLE, J.; MOSER, T. Linkages among land-use, water quality, physical habitat conditions and lotic diatom assemblages: A multi-spatial scale assessment. Hydrobiologia, v. 515, n. 1-3, p. 59-73, 2004.

PAN, Y.; STEVENSON, R. J.; HILL, B. H.; HERLIHY, A. T.; COLLINS, G. B. Using diatoms as indicators of ecological conditions in lotic systems: a regional assessment. Journal of the North American Benthological Society, v. 15, n. 4, p. 481-495, 1996.

PAPPAS, J. L.; STOERMER, E. F. Quantitative method for determining a representative algal sample count. Journal of Phycology, v. 32, n. 4, p. 693-696, 1996.

PAUL, M. J.; MEYER, J. L. Streams in the urban landscape. In: MARZLUFF, J. M.; SHULENBERGER, E.; ENDLICHER, W.; ALBERTI, M.; BRADLEY, G.; RYAN, C.; SIMON, U.; ZUMBRUNNEN, C. (Ed.). Urban ecology: an international perspective on the interaction between humans and nature. New York: Springer, 2008. p. 207-231.

PERES-NETO, P. R.; JACKSON, D. A. How well do multivariate data sets match? The advantages of a Procrustean superimposition approach over the Mantel test. Oecologia, v. 129, n. 2, p. 169-178, 2001.

PORTER, S. D.; MUELLER, D. K.; SPAHR, N. E.; MUNN, M. D.; DUBROVSKY, N. M. Efficacy of algal metrics for assessing nutrient and organic enrichment in flowing waters. Freshwater Biology, v. 53, n. 5, p. 1036-1054, 2008.

POTAPOVA, M.; CHARLES, D. F. Diatom metrics for monitoring eutrophication in rivers of the United States. Ecological Indicators, v. 7, n. 1, p. 48-70, 2007.

POTAPOVA, M.; CHARLES, D. F. Distribution of benthic diatoms in U.S. rivers in relation to conductivity and ionic composition. Freshwater Biology, v. 48, n. 8, p. 1311-1328, 2003.

RIMET, F. Benthic diatom assemblages and their correspondence with ecoregional classifications: case study of rivers in north-eastern France. Hydrobiologia, v. 636, n. 1, p. 137-151, 2009.

ROHLF, F. J.; SLICE, D. Extensions of the Procrustes method for the optimal superimposition of landmarks. Systematic Biology, v. 39, n. 1, p. 40-59, 1990.

RUMRICH, U.; LANGE-BERTALOT, H.; RUMRICH, M. Diatoms of the Andes from Venezuela to Patagonia/Tierra del Fuego and two additional contributions. In: LANGE-BERTALOT, H. (Ed.). Iconographia diatomologica: annotated diatom micrographs. Ruggell: A. R. G. Gantner Verlag K. G., 2000. v. 9, p. 1-673.

SALOMONI, S. E.; ROCHA, O.; CALLEGARO, V. L.; LOBO, E. A. Epilithic diatoms as indicators of water quality in Gravatai river, Rio Grande do Sul, Brasil. Hydrobiologia, v. 559, n. 1, p. 233-246, 2006.

SOININEN, J.; PAAVOLA, R.; MUOTKA, T. Benthic diatom communities in boreal streams: community structure in relation to environmental and spatial gradients. Ecography, v. 27, n. 3, p. 330-342, 2004.

SONNEMAN, J. A.; WALSH, C. J.; BREEN, P. F.; SHARPE, A. K. Effects of urbanization on streams of the Melbourne region, Victoria, Australia. II. Benthic diatom communities. Freshwater Biology, v. 46, n. 4, p. 553-565, 2001.

STEVENSON, R. J.; PAN, Y. Assessing environmental conditions in rivers and streams with diatoms. In:

STOERMER, E. F.; SMOLL, J. P. (Ed.). The diatoms: applications for the environmental and earth sciences. Cambridge: Cambridge University Press, 1999. p. 11-40.

STOERMER, E. F.; SMOLL, J. P. Applications and uses of diatoms: prologue. In: STOERMER, E. F.; SMOLL, J. P. (Ed.). The diatoms: applications for the environmental and earth sciences. Cambridge: Cambridge University Press, 1999. p. 3-8.

TISON, J.; PARK, Y. S.; COSTE, M.; WASSON, J. G.; ECTOR, L.; RIMET, F.; DELMAS, F. Typology of diatom communities and the influence of hydroecoregions: A study on the French hydrosystem scale. Water Research, v. 39, n. 14, p. 3177-3188, 2005.

TREMARIN, P. I.; LUDWIG, T. A. V.; BERTOLLI, L. M.; FARIA, D. M.; COSTIN, J. C. Gomphonema Ehrenberg e Gomphosphenia Lange-Bertalot (Bacillariophyceae) do rio Mauricio, Parana, Brazil. Biota Neotropica, v. 9, n. 3, p. 1-20, 2009.

TREMARIN, P. I.; MOREIRA-FILHO, H.; LUDWIG, T. A. V. Pinnulariaceae (Bacillariophyceae) do rio Guaraguacu, bacia hidrografica litoranea paranaense, Brasil. Acta Botanica Brasilica, v. 24, n. 2, p. 330-348, 2010.

VAN DAM, H.; MERTENS, A.; SKINDELAM, J. A coded checklist and ecological indicator values of freshwaters diatoms from Netherlands. Netherland Journal of Aquatic Ecology, v. 28, n. 1 , p. 117-133, 1994.

WALKER, C. E.; PAN, Y. Using diatom assemblages to assess urban stream conditions. Hydrobiologia, v. 561, n. 1, p. 179-189, 2006.

WALSH, C. J.; ROY, A. H.; FEMINELLA, J. W.; COTTINGHAM, P. D.; GROFFMAN, P. M.; MORGAN, R. P. The urban stream syndrome: current knowledge and search for a cure. Journal of the North American Benthological Society, v. 24, n. 3, p. 706-723, 2005.

WHITTON, B. A.; KELLY, M. G. Use of algae and other plants for monitoring rivers. Australian Journal of Ecology, v. 20, n. 1, p. 45-56, 1995.

WINEMILLER, K. O.; AGOSTINHO, A. A.; CARAMASCHI, E. P. Fish ecology in tropical stream. In: DUDGEON. D. (Org.). Tropical stream ecology. London: Academic Press, 2008. p. 107-140.

WINTER, J. G.; DUTHIE, H. C. Effects of urbanization on water quality, periphyton and invertebrate communities in a southern Ontario stream. Canadian Water Resources Journal, v. 23, n. 3, p. 245-257, 1998.

WINTER, J. G.; DUTHIE, H. C. Epilithic diatoms as indicators of stream total N and total P concentration. Journal of the North American Benthological Society, v. 19, n. 1, p. 32-49, 2000a.

WINTER, J. G.; DUTHIE, H. C. Stream epilithic, epipelic and epiphytic diatoms: habitat fidelity and use in biomonitoring. Aquatic Ecology, v. 34, n. 4, p. 345-353, 2000b.

ZIMMERMAN, G. M.; GOETZ, H.; MIELKE, P. W. Use of an improved statistical method for group comparisons to study effects of prairie fire. Ecology, v. 66, n. 2, p. 606-611, 1985.

Received on August 7, 2012.

Accepted on March 20, 2013.

License information: This is an open-access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Carina Moresco * and Liliana Rodrigues

Programa de Pos-graduacao em Ecologia de Ambientes Aquaticos Continentais, Nucleo de Pesquisas em Limnologia, Ictiologia e Aquicultura, Universidade Estadual de Maringa, Av. Colombo, 5790, 87020-230, Maringa, Parana, Brazil. * Author for correspondence.


Table 1. Mean and standard deviation of abiotic variables
from the urban (U) and rural (R) streams. Sampling period
1 (July/07), 2 (September/07), 3 (December/07), 4 (February/08),
5 (April/08) and 6 (May-June/08). T = water temperature
([degrees]C), pH = potential of hydrogen, Cond. = conductivity
([micro]S [cm.sup.-1]), DO = dissolved oxygen (mL [L.sup.-1]),
P[O.sub.4] = orthophosphate (mg [L.sup.-1]), TN = total nitrogen
(mg [L.sup.-1]), COD = chemical oxygen demand (mg [L.sup.-1]),
BD[O.sub.5] = biochemical oxygen demand (mg [L.sup.-1])
and flow (m [s.sup.-1]).

                T                         pH

R1     16.2 ([+ or -] 0.8)        6.5 ([+ or -] 0.1)
R2     20.4 ([+ or -] 0.3)        7.4 ([+ or -] 0.2)
R3     21.9 ([+ or -] 0.7)        7.1 ([+ or -] 0.1)
R4     21.5 ([+ or -] 0.4)        7.0 ([+ or -] 0.1)
R5     20.9 ([+ or -] 0.5)        7.0 ([+ or -] 0.2)
R6     17.6 ([+ or -] 0.6)        6.8 ([+ or -] 0.1)
U1     15.3 ([+ or -] 0.2)        7.5 ([+ or -] 0.1)
U2     18.3 ([+ or -] 0.2)        7.4 ([+ or -] 0.1)
U3     22.4 ([+ or -] 0.1)        7.1 ([+ or -] 0.1)
U4     21.6 ([+ or -] 0.2)        7.0 ([+ or -] 0.1)
U5     20.8 ([+ or -] 0.2)        7.0 ([+ or -] 0.1)
U6     18.4 ([+ or -] 0.3)        6.8 ([+ or -] 0.1)

              Cond.                       DO

R1     105.6 ([+ or -] 8.7)     123.3 ([+ or -] 101.1)
R2     116.5 ([+ or -] 6.6)      8.2 ([+ or -] 0.44)
R3     112.7 ([+ or -] 4.7)      7.8 ([+ or -] 0.29)
R4     119.2 ([+ or -] 5.7)      8.4 ([+ or -] 0.35)
R5      125 ([+ or -] 5.2)       8.7 ([+ or -] 0.25)
R6      119 ([+ or -] 5.6)        8.1 ([+ or -] 0.5)
U1    264.6 ([+ or -] 46.5)      8.5 ([+ or -] 0.47)
U2    294.3 ([+ or -] 57.4)      8.1 ([+ or -] 0.62)
U3    282.6 ([+ or -] 55.4)      7.6 ([+ or -] 1.14)
U4    310.3 ([+ or -] 70.5)      7.9 ([+ or -] 0.29)
U5    297.6 ([+ or -] 59.7)      6.9 ([+ or -] 0.64)
U6    282.6 ([+ or -] 58.2)      6.9 ([+ or -] 0.42)

            P[O.sub.4]                    TN

R1    123.3 ([+ or -] 101.1)   7466.6 ([+ or -] 2138.5)
R2     86.6 ([+ or -] 45.1)     733.3 ([+ or -] 404.1)
R3     93.3 ([+ or -] 23.1)    1633.3 ([+ or -] 230.9)
R4       50 ([+ or -] 10)      1566.6 ([+ or -] 208.1)
R5     62.6 ([+ or -] 7.6)       1700 ([+ or -] 200)
R6      80 ([+ or -] 26.4)       1800 ([+ or -] 200)
U1    43.33 ([+ or -] 20.8)    7466.6 ([+ or -] 808.2)
U2      30 ([+ or -] 17.3)      5600 ([+ or -] 3803.9)
U3     96.6 ([+ or -] 40.4)    6933.3 ([+ or -] 2936.5)
U4    116.6 ([+ or -] 73.7)    7866.6 ([+ or -] 3499.0)
U5     351 ([+ or -] 513.7)     8100 ([+ or -] 3143.2)
U6     86.6 ([+ or -] 30.5)    8433.3 ([+ or -] 3197.3)

               COD                   BO[D.sub.5]

R1     10.9 ([+ or -] 14.0)      1.34 ([+ or -] 1.23)
R2     5.46 ([+ or -] 4.68)      1.83 ([+ or -] 1.02)
R3     3.65 ([+ or -] 2.51)      1.17 ([+ or -] 1.04)
R4     5.96 ([+ or -] 2.50)      1.93 ([+ or -] 0.66)
R5     3.1 ([+ or -] 0.55)       1.53 ([+ or -] 0.05)
R6     3.4 ([+ or -] 1.45)       0.93 ([+ or -] 0.05)
U1     0.48 ([+ or -] 0.25)      0.17 ([+ or -] 0.04)
U2     1.36 ([+ or -] 0.61)      0.7 ([+ or -] 0.43)
U3     1.08 ([+ or -] 0.62)      0.09 ([+ or -] 0.03)
U4     7.3 ([+ or -] 6.11)       2.53 ([+ or -] 2.57)
U5     3.6 ([+ or -] 0.55)        1.1 ([+ or -] 0.1)
U6     2.83 ([+ or -] 1.40)      1.5 ([+ or -] 0.45)


R1     0.22 ([+ or -] 0.03)
R2     0.18 ([+ or -] 0.05)
R3     0.21 ([+ or -] 0.05)
R4     0.24 ([+ or -] 0.07)
R5     0.19 ([+ or -] 0.04)
R6     0.20 ([+ or -] 0.04)
U1     0.19 ([+ or -] 0.04)
U2     0.12 ([+ or -] 0.04)
U3     0.13 ([+ or -] 0.03)
U4     0.16 ([+ or -] 0.04)
U5     0.17 ([+ or -] 0.04)
U6     0.15 ([+ or -] 0.06)

Table 2. Results from the principal component analysis (PCA)
using the matrix of abiotic variables sampled from July/2007
to May-June/2008 in urban and rural streams.

                                 Axis 1    Axis 2    Axis 3

Eigenvalues                       3.063     2.267     1.429
Broken-stick                      2.929     1.928     1.391
% of variance                     30.63    22.673    12.907
T ([degrees]C)                   0.0057    -0.3161   -0.6486
pH                               -0.0886   0.3503    -0.3459
Cond. ([micro]S-1 [cm.sup.-1])   -0.4954   -0.0973   0.1192
DO (mg [L.sup.-1])               0.3535    0.3823    0.0332
P[O.sub.4] (mg [L.sup.-1])       -0.0546   0.3823    0.0332
TN (mg [L.sup.-1])               -0.4041   -0.1062   0.5012
COD (mg [L.sup.-1])              0.3511    -0.4581   0.1154
BO[D.sub.5] (mg [L.sup.-1])      0.3184    -0.4308   0.0501
Oil (mg [L.sup.-1])              -0.2744   -0.1622   -0.376
Flow (m [s.sup.-1])              0.3944    0.0936    0.1717

Table 3. Indicator Species Analysis (IndVal) with indicator
values (IV). Only statistical significant results
(p < 0.05 for the Monte Carlo test) are shown.

                                           IV         Monte Carlo

                                      Rural   Urban        p

Achnanthes exigua Grunow               70      13       0.0066
Achnanthes rupestoides Hohn            78      15       0.0008
Amphipleura lindheimeri Grunow         89       0       0.0001
Amphora copulata (Kutzing)             89       0       0.0001
Schoeman and Archibald
Dipbneis subovalis Cleve               81       0       0.0001
Encyonema mesianum (Cholnoky) Mann     65       0       0.0002
Fallacia ecuadoriana                   44       2       0.0359
Lange-Bertalot and Rumrich
Fallacia insociabilis                  67       0       0.0001
(Krasske) Mann
Fragilaria rumpens (Kutzing)           39       0       0.0069
Frustulia crassinervia (Brebisson)     28       0       0.0473
Lange-Bertalot and Krammer
Frustulia vulgaris (Thwaites)          46       2       0.0094
De Toni
Geissleria neosubtropica               56       0       0.0002
Metzeltin, Lange-Bertalot and
Gomphonema brasiliense Grunow          78       0       0.0001
Gyrosigma acuminatum (Kutzing)         59       0       0.0008
Gyrosigma scalproides (Rabenhorst)     42       0       0.0078
Hantzschia amphioxys (Ehrenberg)       28       0       0.0486
Hippodonta capitata (Ehrenberg)        28       0       0.0462
Lange-Bertalot, Metelzin &
Luticola dapalis (Frenguelli) Mann     48       0       0.0034
Navicula cryptotenella                 67      24       0.0185
Navicula lohmannii                     89       4       0.0001
Lange-Bertalot & Rumrich
Nupela praecipua (Reichardt)           83       8       0.0002
Pinnularia sp.                         28       0       0.0442
Placoneis constans var. symmetrica     56       0       0.0006
(Hustedt) Kobayasi
Placoneis disparilis (Hustedt)         49       0       0.0018
Metelzin & Krammer
Placoneis hambergii (Hustedt)          49       9       0.0297
Sellaphora pupula (Kutzing)            79       5       0.0002
SeUaphora sp. 1                        61       0       0.0001
Sellaphora sp. 2                       76       0       0.0001
Stauroneis cf. kriegeri Patrick        64       1       0.0003

Stenopterobia schweickerdtii           56       0       0.0004
(Cholnoky) Brassac,
Ludwig & Torgan
Tryblionella levidensis Smith          79       1       0.0001
Ulnaria ulna (Nitzsch) Compere         56      11       0.0245
Sellaphora seminulum (Grunow) Mann      8      92       0.0001
Achnanthes lanceolata                   9      91       0.0001
(Brebisson ex Kutzing) Grunow
Achnanthidium minutissimum              2      91       0.0003
(Kutzing) Czarnecki
Craticula sp.                           0      89       0.0001
Nupela sp.                             13      87       0.0018
Fallacia monoculata (Hustedt) Mann      0      83       0.0001
Gomphonema lagenula Kutzing            19      81       0.0415
Amphora montana Krasske                 1      77       0.0063
Eolimna minima (Grunow)                27      73       0.0375
Navicula tenelloides Hustedt            6      69       0.0153
Navicula sp. 5                          1      69       0.0008
Caloneis bacillum (Grunow) Cleve        0      66       0.0007
Pinnularia gibba (Ehrenberg)            1      53       0.0292
Mayamaea atomus (Kutzing)               0      50       0.0028
Lange-Bertalot var. permitis
Cyclotella meneghiniana Kutzing         0      39       0.0129
Nitzschia cf. inconspicua Grunow        0      28       0.0454
Sellaphora sp.                          0      28       0.0444
COPYRIGHT 2014 Universidade Estadual de Maringa
No portion of this article can be reproduced without the express written permission from the copyright holder.
Copyright 2014 Gale, Cengage Learning. All rights reserved.

Article Details
Printer friendly Cite/link Email Feedback
Author:Moresco, Carina; Rodrigues, Liliana
Publication:Acta Scientiarum. Biological Sciences (UEM)
Date:Jan 1, 2014
Previous Article:Establishment and content of sugars and phenols in Physalis callus obtained from different explants and concentrations of bap and naa/Estabelecimento...
Next Article:Diversity and fish distribution at Rodrigo de Freitas Lagoon, Rio de Janeiro State, Brazil, using GIS/Diversidade e distribuicao de peixes na Lagoa...

Terms of use | Privacy policy | Copyright © 2019 Farlex, Inc. | Feedback | For webmasters