Statistical distribution models for migratory fish in Jacuí basin , South Brazil

The aim of the present study was to identify the distribution patterns of migratory fishes in the Jacuí river basin (Rio Grande do Sul, South Brazil), proposing a statistical model of presumed distribution based on geomorphologic environmental data. Through maps of occurrence probability, we hope to contribute to decisions regarding basin environmental management. The analyzed species were: Salminus brasiliensis (dourado), Leporinus obtusidens (piava), Prochilodus lineatus (grumatã) and Pimelodus pintado (pintado). Samples were made through interviews with fishermen and local inhabitants, covering the main channel and tributaries of the rivers Jacuí, Taquari-Antas, Vacacaí, Vacacaí-Mirim, Pardo, Pardinho, Sinos, and Caí. The sampling program resulted in 204 interviews, being 187 considered as valid in 155 different sampling points. The probability of migratory fish occurrence was adjusted through the LOGIT routine of the Idrisi Andes Software: P = e(b0 + b1 . altitude + b2 . basin area) . (1 + e(b0 + b1 . altitude + b2 . basin area))-1, where P is the occurrence probability of the species (0-1) and b0, b1 and b2 are the equation parameters. Model accuracy, for estimating presence, ranged from 82% to 93%. Pimelodus pintado was cited to occur in 121 points among the 155 sampled (78.06%), Prochilodus lineatus in 72 (46.45%), L. obtusidens in 62 (40.00%) and S. brasiliensis in 58 (37.42%). Equation parameters were estimated (± standard error) as follow: S. brasiliensis: b0 = -2.8762 ± 0.2597; b1 = 1.3028 ± 0.0332; b2 = 1.1487 ± 0.0301; L. obtusidens: b0 = -0.8364 ± 0.2213; b1 = -1.5564 ± 0.0462; b2 = 0.9947 ± 0.0206; Prochilodus lineatus: b0 = 0; b1 = -1.3067 ± 0.0544; b2 = 0.8128 ± 0.0177; Pimelodus pintado: b0 = -0.9487 ± 0.3688; b1 = -0.8269 ± 0.0496; b2 = 0.9255 ± 0.0304.


Introduction
More than 15% of the Neotropical fish fauna is composed by migratory species (Carolsfeld et al., 2004), which may need to cover several kilometers to stimulate gonadal maturation (Godoy, 1987).The watershed area for maintaining most of these species could reach 40.000 km 2 (Godoy, 1987).Spawning use to happen upstream the adult feeding areas, and stream current carries the eggs and larvae to areas where they will develop.In these sites, juveniles will feed and grow until they become big enough to join the main stock (Pitcher & Hart, 1982).Malabarba (1989) presents a list of species of freshwater fish of Patos Lagoon drainage, which includes the Jacuí River.Following Agostinho & Júlio (1999) and Vazzoler et al. (1997), the migratory species for Jacuí basin are: Salminus brasiliensis respect, predicting species occurrence through a modeling approach based on geographical information system (GIS) represents a useful methodological tool (Pearce & Ferrier, 2000b).
Generalized linear models (GLM), implemented within a GIS are especially useful as it can be applied with many types of predictors (continuous, binary, qualitative, and ordinal), and the distribution of the species should comprise just presence and absence data (Syartinilia & Tsuyuki, 2008), information available through interviews.These models are then applied to extrapolate the probability of occurrence of species across the entire region of interest (Pearce & Ferrier, 2000b).
Although several ecological factors are known to influence fish distribution patterns, access to such data is not easy or even accurate at large areas.In a large regional scale, the use of simple factors derived from GIS, as altitude and basin area, seems to be enough for understanding distributional patterns along the rivers (Matthews et al., 1992;Pouilly et al., 2006).
The present study aims to identify the longitudinal historical distribution pattern of long distance migratory fish in the hydrographic basin of the Jacuí River, which already has an installed hydroelectric production of 1.428,674 MW, and a potential, still not installed one, of 271,004 MW.Apart from presenting distribution maps for long distance migratory fish, probabilistic distribution models for each species will also be proposed, using altitude and basin area as main presence predictors.

Material and Methods
Eight field expeditions were carried out, covering the most important rivers draining to the Jacuí basin: Taquari-Antas, Vacacaí/Vacacaí-Mirim, Pardo/Pardinho, Sinos and Caí, as their main tributaries.The samplings were performed between February and July of 2007.Expeditions were carried out making use of the road network and its crossing points with the rivers, resulting in 155 different sampling points (Fig. 1).
At each sampling site, local fishermen and oldest residents were searched and interviewed.With the aid of a board containing a picture of interest species, each interviewed person was asked about their knowledge of the fishes, and then, their presence or absence in that specific river segment (present or historical).In case of more than one interview at the same sampling site, only the most frequent information or consensus was registered.The project totalized 204 interviews, 187 registered as useful in 155 different points.
Also, a survey in the main Brazilian fish collections was done by using the SIBIP/NEODAT III database (MCP, MNRJ, MZUSP, UFRGS); besides literature information (Petry & Schulz, 2006).This survey was extended through the analysis of Environmental Impact Reports (EIA-RIMA) produced for the licensing power-dams (UHE) in the Jacuí and Taquari-Antas basins (Dona Francisca; 14 de Julho, Castro Alves, Monte Claro).
Data analysis and cartographic products were performed  Loureiro, 2008 (pintado).The spotted catfish of the Patos basin was known until recently as Pimelodus maculatus Lacépede, 1803, a species with wide distribution in neotropics and maybe a species complex.After Azpelicueta, Lundberg & Loureiro (2008) the Patos spotted catfish was reviewed and described as Pimelodus pintado, with few information concerning migratory or reproductive biology but maybe with similar biological patterns as Pimelodus maculatus.
Nevertheless, including Pimelodus pintado as a migratory fish species is still controversial.Braun (2005) suggests that the species seems to reproduce in little channels, tributaries of Patos Lagoon.This information is reinforced by the occurrence of Pimelodus pintado in Barros lagoon, which receives water from only small creeks of first order, suggesting no dependence of large scale migration for reproductive purpose.However, considering the importance of this species for artisanal fishery, Pimelodus pintado will be included among the compulsory migrant species of the Jacuí River basin.
A river is a complex system, with continuous changing structure from upstream to lower areas.The species richness increases downwards (Reyes-Gavilán et al., 1996) and in higher altitudes both physical and chemical conditions are more stressful, with a lower number of species showing physiological and ecological capacities to survive in such places (Matthews & Styron-Júnior, 1981).
As a general rule, the distributional limit of a species is driven by several factors, which combine or interact to create the patterns found in nature (Hall et al., 1992).The altitude, for example, seems to be broadly related to hydrodynamics and morphology of a river.With the combined factors such as water temperature, conductivity, pH, current speed, declivity, and the presence of pools and waterfalls, altitude is capable of creating an environmental gradient that exerts influence on the species richness in each region (Pouilly et al., 2006).Araújo & Tejerina-Garro (2009), found that pH, water velocity, channel width and water temperature strongly affected fish assemblages in upper Paraná tributaries.Using geomorphologic data, Súarez & Petrere-Júnior (2007) suggested altitude as the most important structuring factor to fish community in the Iguatemi River basin.Apart from altitude, the water volume drained to a specific point seems to exert a large influence on the fish distribution patterns (Garutti, 1988).When the altitude decreases, there is an increase in the width and deep of the river, allowing the survival of bigger fish species (Vázquez et al., 1979).
Brazil has approximately 4360 dams, not taking into account homemade barriers and/or non-registered ones (Silveira & Cruz, 2005).Currently, 90% of the Brazilian electric power production is generated by hydroelectric turbines (Petrere-Júnior et al., 2002).Hydroelectric plants have been interfering in fish migration all around the world.Knowing the longitudinal distribution of migratory fishes is an important tool to support decisions related to dams planning in a regional scale.In this using Idrisi Andes software (Clarck Labs, 2006) and a digital elevation model (DEM, radar altitude in pixel of 92.6 per 92.6 m) for the Brazilian official reference system (SAD 69;Labgeo, 2006).
The hydrographic matrix was obtained through four basic steps: (1) DEM homogenization by applying the Filter Min 3x3 pixels; (2) defining pathways with monotonically decreasing attitudes through Pit Removal tool; (3) by applying the Runoff tool to estimate the upstream catchment area for each pixel, correcting pixel number for the real area (km²) by multiplying pixel values by 0.00857476; (4) converting each pixel with catchment area smaller then 10 to zero and bigger then 10 into one, and then creating a hydrographic matrix with watershed bigger then 10 km².Hydrographic matrix with fish presence, for each species, was created through Pathway tool, using recorded presence (GPS coordinates) as Target Image and the DEM as Cost Surface.Mask for data analysis was obtained by using also the Pathway tool, with upstream sampled coordinates as Target Image and the DEM as Cost Surface.
The occurrence probability for each migratory fish species throughout the hydrographic matrix was estimated using the Multinomial Logistic Regression routine (Logisticreg, adjust by maximum likelihood through the Newton-Raphson algorithm): P (X = 1) = e (b0 + b1 .altitude + b2 .basin area) .(1 + e (b0 + b1 .altitude + b2 .basin area) ) -1 , where P (X = 1) is the occurrence probability of the species and b 0 , b 1 and b 2 are the equation parameters.The measured variables, the measurement scale and data source is described in Table 1.
Data analysis was performed in two steps.First, the probability model for each species was estimated using data from 175 interviews from 150 different points, including literature, EIA-RIMA and museum data.A residual map was generated showing differences among observed presence (0/1) and predicted presence probability.Sampling points with residuals larger then 0.5 were revisited for new interviews.During this new field trip, additional 29 interviews and five new points were introduced into de original data matrix for new model processing.Adherence between informed distribution and presumed distribution were estimated as the percentage of the river segments with presence/absence informed as correct, according to the presumed distribution model (considered as presumed presence probabilities higher than 0.5 and as presumed absence probabilities lower than 0.5).
Considering that the Logisticreg function of the Idrisi Andes software does not estimate the standard errors of the adjusted parameters, 10 independent adjustments were carried out for each species through random selection of 10% of the pixels from de sample matrix.The standard errors were calculated as the standard deviation of the different estimates for each parameter, obtaining the significance thorough the Wald statistic, the ratio between the estimated parameter and its standard deviation, calculating P through the Z distribution (Tabachnick & Fidell, 1996). in 58 (37.42%).
Is shown in Fig. 2 the cumulative frequency of informed occurrence as a function of altitude and catchment area.Results were similar for all the species but Pimelodus pintado.As identified, Prochilodus lineatus, L. obtusidens and S. brasiliensis have an altitudinal threshold of about 280-300 m whereas Pimelodus pintado was recorded in altitudes exceeding 500 m (one informed occurrence at 680 m).Occurrences, concerning basin area, were recorded for all species in river segments with catchments over 20 km².
The estimated parameters for the probability model for each species are shown in Table 2. Table 3 presents the adherence between informed distribution and presumed distribution.
For each analyzed species, Figs.
Considering the effects of dams as barriers to fish migratory movements, Fig. 7 presents the cumulative distribution of large and small power dams of the Jacuí basin (n = 24, data from National Electric Power Agency -ANEL) according to an altitudinal gradient.As could be identified, around 50% of registered dams are present above the altitudinal threshold of 280-300 m. Figure 8 presents the geographical distribution of power dams (yellow spots) in the Jacuí basin according to altitudinal categories of under 300 m (red river segments) and over 300 m (green river segments).Inspecting the spatial distribution, it becomes clear that most of the constructed power dams are located above or near the altitudinal threshold of 280-300 m of upper migratory fish distribution, suggesting that most operating power dams presents little effect as barriers for upstream movements for migratory fish.

Discussion
In the last twenty years, the use of multivariate statistic in the modeling of species distribution has increased through the application of a large variety of techniques.Particularly, the regression models have been broadly used to predict distribution, abundance and habitat preference of species (Brito et al., 1999).However, when this information is joined with Geographic Information Systems (GIS), it is possible to incorporate ecological factors directly obtained from remote images to prediction models.Also, indirect gradients constituted by variables that do not show direct physiological action, as altitude, declivity and topographic position, are easier to be obtained than environmental parameters of direct gradient, such temperature or pH, which often introduce spatial uncertainties due to the lack of data and interpolation bias (Guisan et al., 1999).
Although in the present work we have opted to model distribution based just in two parameters of indirect gradient, altitude and catchment area, the predicted distribution, when confronted with the informed distribution, always produced an adherence ratio higher than 80%.The combined effect of both variables suggest that altitude is the main factor limiting migratory fish distribution for areas comprising the Serra Geral (altitudinal areas) whereas catchment area limits distribution in the plains of the Central Depression of Rio Grande do Sul.This result suggests that even a relatively small set of geomorphologic variables are able to show an important predictive power, maybe due to colinearity effects with associated direct parameter, such as water temperature, river dimension (deep, width), current speed and turbidity.
Nevertheless, even considering adherence over than 80% as satisfactory, the resulted error should not be neglected.Among the factors contributing to the observed residues are: (1) information mismatch; (2) deficiency in determining the exact point of species disappearance into a river system; (3) influence of other environmental factors, such as barrier downstream; and (4) the intrinsic variability related to any natural phenomena.
Although we have tried to minimize the probability of mistaken information, through a resample process of discrepant interviews (validation), this process was limited to river segments in which residuals (both positive and negative) were higher than 0.5 for at least two species, due to the large geographical area under analysis.Also, new interviews resulted mostly into confirmation of previous data and no significant change was obtained.
On the other hand, the methodological approach used in this work, statistic operation of images instead of direct numeric data, presents advantages and disadvantages that need closer examination.The main advantage consists in increasing the available data set.Pearce & Ferrier (2000b) identified a sharp increase in predictive performance of GLM as sample size is increased from 50 to 250 sites.By applying conventional statistic software, each point is an independent sample, with measured attributes.In the present work, we have 155 different points spread in the Jacuí basin.However, by using Geographic Information System (GIS) and statistic operation of images, additional information could be incorporated.When the presence of a species is registered in adjacent points of a river, it is reasonable to suppose that the species do also occurs along the whole segment between these two points.By using image operation, this information can be retrieved, increasing the sampling effort.In the other hand, when a species is not registered in a sampling point, is reasonable to consider that it will be absent from that point to upstream.
With the aid of this procedure, a matrix of around 40 thousand information points (pixels) was created, with presence/absence data for each species, as well as the information of altitude and catchment area.This protocol allowed us to perform independent adjustments for each species, using in each only 10% of the sampled points (approximately four thousand points), in a bootstrap process to estimate the confidence intervals for the adjusted coefficients.
The process, however, could incorporate bias at an inverse proportion to the density of the sample matrix.When a species disappears between two sampling points, this analysis protocol considers it as absent from the last point with informed presence.If the points are distant apart, a considerable segment of probable presence will be computed as absence.The bigger the density of sampled points in the distribution limits, the smaller the errors from the analysis process.This kind of possible bias was the reason driving our priority for sampling in the upstream segments of main rivers instead of the lower areas.
The error margin could even have suffered the influence of other factors not included in the model, as the presence of downstream barriers.Moreover, the incorporation of natural barriers into the present model was presently unviable.Although a lot of natural barriers have been identified in loco, the DEM precision by using pixel of 92 per 92 meters and altitude information coming from a radar sensor, appeared to be inadequate for the identification of potentially significant small accidents, so that we opted to not include this variable in the proposed model.Also, to include information related to dam presence brings other problems as dam age and the presence and effectiveness of transposition facilities.So, the final distributional models of the present work intend to identify the historical distribution of migratory fish along the Jacuí basin, do not including the subsequent interference of constructed dams or barriers.This historical distribution, however, could be biased (underestimated) if interviewed residents were younger than operating power dams.Large power-dams in the basin are operating since 1952 in the Caí River (UHE Bugres, 544 m altitude), 1962 in the upper Jacuí (UHE Jacuí, 270 m altitude) and 2004 in the Taquari-Antas (UHE Monte Claro, 150 m altitude) allowing to identify undisturbed fish presence from old residents for most of the basin.Although a lot of small units are much older and widespread, they are located mainly in small tributaries and upper altitudes, with small probability of fish occurrence.By interpreting each species distribution pattern, it is interesting to identify that Pimelodus pintado presents a relative low accuracy in relation to estimated absence (48%).This low adherence can be related to two factors: (1) possible intentional introduction of the species in upstream areas, or (2) the natural transposition of barriers through dam flooding.The introduction of the Pimelodus pintado was informed by residents in the surroundings of the Salto System.According to them, individuals of Pimelodus pintado were introduced in the flooded area of the Divisa Dam (Caí River).Also residents of the upper Jacuí reported that the species started to be captured only after the construction of Maia Filho and Passo Real Dams, suggesting possible transposition of natural barriers due to lake flooding.
By observing the species presumed distribution maps, it is possible to verify that a considerable superposition has occurred.The altitude coefficients obtained for the four species showed a negative influence of this parameter over the analyzed species distribution.With higher altitudes the smaller habitat heterogeneity and the dynamic instability determine conditions in which few species are able to tolerate (Bistoni & Hued, 2002).In rivers of central Argentina, distribution of the species according to an altitudinal gradient was also observed: Pimelodus albicans reaches altitudes between 600 and 800 m; Prochilodus platensis (= P. lineatus), between 400 and 600 m; and Salminus maxillosus (= S. brasiliensis), between 200 and 400 m (Bistoni & Hued, 2002).In the present study, we observed that Pimelodus pintado (680 m) and Prochilodus lineatus (297 m) occurred in lower altitudes than their co-generic species from Argentina, although S. brasiliensis (292 m) presented the same altitudinal pattern of S. maxillosus (= S. brasiliensis).Even so, in both cases the genus maintained the same altitudinal sequence, with Pimelodus reaching upstream areas, suggesting the existence of a phylogenetic component in the habitat distribution pattern.
The upstream catchment area appears to be extremely significant to explain the distribution pattern of the four species, with positive coefficients, showing the preference of the species for downstream segments.When reported for small upstream segments, migratory fish presence was frequently associated with unusual floods and not as a regular occurrence.However, direct comparisons among the occurrence probabilities of the species should not be done taking into account each coefficient independently.The nature of the multiple regression models, due to the variables multicolinearity, may result in difficulties to correctly identify the individual effects of each coefficient (Hise et al., 1983).However, S. brasiliensis and Prochilodus lineatus, which presented almost the same altitude coefficients, can be directly compared in relation to the catchment area coefficient.Despite suffering similar influence from the altitude over their distribution, the S. brasiliensis seems to be more affected by water volume, requiring larger river dimension to fix its occurrence.
Available information suggest that S. brasiliensis and Prochilodus lineatus need at least 100 km of unobstructed river segments to complete gonad development and the spawning process (Agostinho et al., 2004).There is no information concerning Pimelodus pintado, although Pimelodus maculatus needs shorter segments to spawn (Godoy, 1967), apparently being a species which is less affected by artificial barriers.
Concerning the usefulness of distribution predictability, Generalized Linear Models (GLM) could be interpreted by two distinct ways (Pearce & Ferrier, 2000a).Reliable predictions may be used as if each predicted probability can be an accurate estimate of the likelihood of detecting the species at a given area.For example, an estimated 0.8 probability means that the species has a real likelihood of presence of 80%.On the other hand, a model with good discrimination ability is the one that can only discriminate between occupied (p > 0.5) and unoccupied sites (p < 0.5), but estimated probabilities may be not accurate (Pearce & Ferrier, 2000a).
The proposed models (standard GLM), imply that the effect of each environmental variable is linear (multiplied by a constant), that is, an increase of the value of the variable from, for example, 100 to 200 has the same effect than an increase from 500 to 600, despite the sigmoid pattern of the probability distribution.Nevertheless, the effect of an environmental variable may follow a more complex pattern, emulating a linear effect only on a specific range of values.Some species, for example, may occur in the middle segment of large rivers, and may be rare upstream or downstream.For this kind of distribution neither altitude nor basin area would have linear effects from zero to high values; they may present an almost linear effect inside a more restricted range.This implies that if the behavior of any predictor in not clearly known, the likelihood of presence predicted by the GLM should be used only for discriminatory purpose: the species is believed to be present when likelihood of presence is above 0.5 and absent with lower values.
Finally, the Geographic Information System (GIS) is a tool with increasing importance for environmental management.In this respect, all products of this work have the objective to contribute for an integrated environmental planning of the Jacuí River Hydrographic Basin.Moreover, the presumed distribution models can be useful to subsidize decisions related to the construction of new power dams.In this respect our data suggests that power dam planning, as potential barriers for fish migration, should consider the altitudinal threshold of 280-300 m as upstream limits of migratory fish in the Jacuí basin.Also, the proposed distribution models do not intend to substitute in loco sampling for verification of species occurrence.We intended, more than anything, to provide to the licensing agencies an independent tool to evaluate the environmental field data related to licensing of new dam.
Fig. 2. Inverse cumulative frequency of informed occurrence of migratory fish of Jacuí basin according to an upstream gradient of altitude (m) and basin catchment (km²).

Fig. 3 .
Fig. 3. Salminus brasiliensis.(a) Hydrographic map of presumed occurrence (blue) through interviews (yellow spots), literature and museums (purple spots) over sampled hydrography (red).(b) Probability map of the presumed occurrence through the adjusted logistic model to the whole hydrographic basin with a catchment area higher than 10 km 2 .(c) Distribution map of residuals applied to sampled hydrographic basin.Red segments represent negative residuals less then -0.34, indicating predicted presence but with no registered occurrence.Yellow segments represent positive residuals greater then 0.34, indicating registered occurrence do not predicted by the model.Hydrographic maps generated through a digital elevation model (DEM, radar altitude in pixel of 92.6 per 92.6 m) for the Brazilian official reference system (SAD 69; Labgeo, 2006).

Fig. 4 .
Fig. 4. Prochilodus lineatus.(a) Hydrographic map of presumed occurrence (blue) through interviews (yellow spots), literature and museums (purple spots) over sampled hydrography (red).(b) Probability map of the presumed occurrence through the adjusted logistic model to the whole hydrographic basin with a catchment area higher than 10 km 2 .(c) Distribution map of residuals applied to sampled hydrographic basin.Red segments represent negative residuals less then -0.34, indicating predicted presence but with no registered occurrence.Yellow segments represent positive residuals greater then 0.34, indicating registered occurrence do not predicted by the model.Hydrographic maps generated through a digital elevation model (DEM, radar altitude in pixel of 92.6 per 92.6 m) for the Brazilian official reference system (SAD 69; Labgeo, 2006).

Fig. 5 .
Fig. 5. Leporinus obtusidens.(a) Hydrographic map of presumed occurrence (blue) through interviews (yellow spots), literature and museums (purple spots) over sampled hydrography (red).(b) Probability map of the presumed occurrence through the adjusted logistic model to the whole hydrographic basin with a catchment area higher than 10 km 2 .(c) Distribution map of residuals applied to sampled hydrographic basin.Red segments represent negative residuals less then -0.34, indicating predicted presence but with no registered occurrence.Yellow segments represent positive residuals greater then 0.34, indicating registered occurrence do not predicted by the model.Hydrographic maps generated through a digital elevation model (DEM, radar altitude in pixel of 92.6 per 92.6 m) for the Brazilian official reference system (SAD 69; Labgeo, 2006).

Fig. 6 .
Fig. 6.Pimelodus pintado.(a) Hydrographic map of presumed occurrence (blue) through interviews (yellow spots), literature and museums (purple spots) over sampled hydrography (red).(b) Probability map of the presumed occurrence through the adjusted logistic model to the whole hydrographic basin with a catchment area higher than 10 km 2 .(c) Distribution map of residuals applied to sampled hydrographic basin.Red segments represent negative residuals less then -0.34, indicating predicted presence but with no registered occurrence.Yellow segments represent positive residuals greater then 0.34, indicating registered occurrence do not predicted by the model.Hydrographic maps generated through a digital elevation model (DEM, radar altitude in pixel of 92.6 per 92.6 m) for the Brazilian official reference system (SAD 69; Labgeo, 2006).

Fig. 8 .
Fig. 8. Hydrographic map of the Jacuí basin with dam distribution (yellow spots).Red river segments are below the altitudinal threshold of 300 m (green above this limit).Hydrographic maps generated through a digital elevation model (DEM, radar altitude in pixel of 92.6 per 92.6 m) for the Brazilian official reference system (SAD 69; Labgeo, 2006).

Table 1 .
Measured variables, measurement scale and data source for probability distribution models adjustment for migratory fishes at Jacuí basin (Brazil).

Table 2 .
Average, Standard Deviation, Wald and significance of the coefficients (Intercept, b 0 ; altitude, b 1 ; and basin area, b 2 ) estimated from ten independent adjustments of the logistic regression performed using the software Idrisi Andes to each of the four migratory fish species in Jacuí basin (Brazil).S.D. = Standard Deviation.

Table 3 .
Mean accuracy percentage (adherence) of the presumed occurrence model for each one of the four migratory fish species in Jacuí basin (Brazil).