Trophic dynamics of methylmercury and trace elements in a remote Amazonian Lake.

Information on pollutant trophodynamics can be crucial for public health, as contaminated food consumption may lead to deleterious effects. This study was performed in Puruzinho Lake, a remote body of water in the Brazilian Amazon from which a riparian human population obtains an important part of its animal protein intake. Samples from 92 individuals, comprising 13 species and four trophic guilds (iliophagous, planktivorous, omnivorous, and piscivorous fish) were analysed for the determination of trace elements (Fe, Cr, Mn, Ni, Zn, Ca, Sr, Cd, Sn, Tl and Pb) and methylmercury concentrations. Samples from the same individuals had already been analysed for stable isotope (SI) measurements (δ13C and δ15N) in a previous investigation and the SI data have been statistically treated with those generated in this study for the evaluation of trophic dynamics of contaminants. Methylmercury was the only analyte that biomagnified, presenting TMF values of 4.65 and 4.55 for total and resident ichthyofauna, respectively. Trace elements presented either trophic dilution or independence from the trophic position, constituting a behaviour that was coherent with that found in the scientific literature. The similarity between Ni behaviour through the trophic web to that of essential elements contributes to the discussion on the essentiality of this metal to fish. Considering the Non-cancer Risk Assessment, the calculated Target Hazard Quotient (THQ) values were higher than 1.0 for all analysed individuals for methylmercury, as well as for only one individual for nickel. No other analyte rendered THQ values higher than 1.0.

Bacterial Communities Along Environmental Gradients in Tropical Soda Lakes.

Soda lake environments are known to be variable and can have distinct differences according to geographical location. In this study, we investigated the effects of different environmental conditions of six adjacent soda lakes in the Pantanal biome (Mato Grosso do Sul state, Brazil) on bacterial communities and their functioning using a metagenomic approach combined with flow cytometry and chemical analyses. Ordination analysis using flow cytometry and water chemistry data from two sampling periods (wet and dry) clustered soda lakes into three different profiles: eutrophic turbid (ET), oligotrophic turbid (OT), and clear vegetated oligotrophic (CVO). Analysis of bacterial community composition and functioning corroborated this ordination; the exception was one ET lake, which was similar to one OT lake during the wet season, indicating drastic shifts between seasons. Microbial abundance and diversity increased during the dry period, along with a considerable number of limnological variables, all indicative of a strong effect of the precipitation-evaporation balance in these systems. Cyanobacteria were associated with high electric conductivity, pH, and nutrient availability, whereas Actinobacteria, Alphaproteobacteria, and Betaproteobacteria were correlated with landscape morphology variability (surface water, surface perimeter, and lake volume) and with lower salinity and pH levels. Stress response metabolism was enhanced in OT and ET lakes and underrepresented in CVO lakes. The microbiome dataset of this study can serve as a baseline for restoring impacted soda lakes. Altogether, the results of this study demonstrate the sensitivity of tropical soda lakes to climate change, as slight changes in hydrological regimes might produce drastic shifts in community diversity.

Empirical evidence for resilience of tropical forest photosynthesis in a warmer world.

Tropical forests may be vulnerable to climate change1-3 if photosynthetic carbon uptake currently operates near a high temperature limit4-6. Predicting tropical forest function requires understanding the relative contributions of two mechanisms of high-temperature photosynthetic declines: stomatal limitation (H1), an indirect response due to temperature-associated changes in atmospheric vapour pressure deficit (VPD)7, and biochemical restrictions (H2), a direct temperature response8,9. Their relative control predicts different outcomes-H1 is expected to diminish with stomatal responses to future co-occurring elevated atmospheric [CO2], whereas H2 portends declining photosynthesis with increasing temperatures. Distinguishing the two mechanisms at high temperatures is therefore critical, but difficult because VPD is highly correlated with temperature in natural settings. We used a forest mesocosm to quantify the sensitivity of tropical gross ecosystem productivity (GEP) to future temperature regimes while constraining VPD by controlling humidity. We then analytically decoupled temperature and VPD effects under current climate with flux-tower-derived GEP trends in situ from four tropical forest sites. Both approaches showed consistent, negative sensitivity of GEP to VPD but little direct response to temperature. Importantly, in the mesocosm at low VPD, GEP persisted up to 38 °C, a temperature exceeding projections for tropical forests in 2100 (ref. 10). If elevated [CO2] mitigates VPD-induced stomatal limitation through enhanced water-use efficiency as hypothesized9,11, tropical forest photosynthesis may have a margin of resilience to future warming.

Mapping carbon and nitrogen isotopic composition of fingernails to demonstrate a rural-urban nutrition transition in the Center-West, Northeast, and Amazon regions of Brazil.

OBJECTIVE: The main objective of this study is to investigate diet patterns among rural and urban populations of the Center-West, Northeast, and Amazon regions of Brazil through the carbon and nitrogen isotopic composition of fingernails, recognizing that the extent of market integration is a key driver of food consumption. MATERIALS AND METHODS: In the Center-West, Northeast, and Amazon regions of Brazil, fingernails were sampled in clusters encompassing a major city, town, and rural village. A total of 2,133 fingernails were analyzed. Fingernails were clipped by donors using fingernail clippers. In the laboratory, samples were cleaned then weighed in small tin capsules before being isotopically analyzed for carbon and nitrogen. RESULTS: The overall mean δ13 C and δ15 N were -19.7 ± 2.8‰ and 10.6 ± 1.1‰, respectively. In the more remote villages, where access to food markets is more challenging, lower δ13 C prevails, suggesting that Brazilian staple foods (rice, beans, and farinha) still dominate. In areas with easier access to food markets, δ13 C values were higher, suggesting a change to a diet based on C4 plants, typical of a Brazilian supermarket diet. The variability among inhabitants in the same location expressed by a significant inverse correlation between δ13 C and δ15 N fingernail values suggested that "market integration" does not affect everyone equally in each community. DISCUSSION AND CONCLUSION: The nutrition transition has not yet reached some remote villages in these regions of Brazil or that the nutrition transition has not yet reached all residents of these remote villages. On the other hand, in several villages there is a considerable adherence to the supermarket diet or that some residents of these villages are already favoring processed food.

Identifying the composition of commercial Brazilian cat food with stable isotopes of carbon and nitrogen.

Considering the increasing pet owner's concern about the food their pets are consuming, in this study we investigated the origin of the main ingredients in wet and dry foods produced in Brazil using stable isotope ratios of carbon and nitrogen. We concluded that chicken and pork seem to be the dominant ingredients in most of the samples, with larger proportions in wet cat food. Even in pet foods showing 'beef' as the main ingredient on the label, we found a low proportion of bovine products in both wet and dry cat foods. Comparing the contribution of plant-derived products (C3 and C4 plants) and animal-derived products (chicken-pork, bovine and fish), approximately 21 % of cat foods had more than 30 % of ingredients with plant origin in their composition. The high amount of plant-derived products in cat foods found here raises the question whether this should be mentioned on package labels.

Carbon and Nitrogen Isotope Ratios of Food and Beverage in Brazil.

Several previous studies on targeted food items using carbon and nitrogen stable isotope ratios in Brazil have revealed that many of the items investigated are adulterated; mislabeled or even fraud. Here, we present the first Brazilian isotopic baseline assessment that can be used not only in future forensic cases involving food authenticity, but also in human forensic anthropology studies. The δ13C and δ15N were determined in 1245 food items and 374 beverages; most of them made in Brazil. The average δ13C and δ15N of C3 plants were -26.7 ± 1.5‱, and 3.9 ± 3.9‱, respectively, while the average δ13C and δ15N of C4 plants were -11.5 ± 0.8‱ and 4.6 ± 2.6‱, respectively. The δ13C and δ15N of plant-based processed foods were -21.8 ± 4.8‱ and 3.9 ± 2.7‱, respectively. The average δ13C and δ15N of meat, including beef, poultry, pork and lamb were -16.6 ± 4.7‱, and 5.2 ± 2.6‱, respectively, while the δ13C and δ15N of animal-based processed foods were -17.9 ± 3.3‱ and 3.3 ± 3.5‱, respectively. The average δ13C of beverages, including beer and wine was -22.5 ± 3.1‱. We verified that C-C4 constitutes a large proportion of fresh meat, dairy products, as well as animal and plant-based processed foods. The reasons behind this high proportion will be addressed in this study.

Evolutionary diversity is associated with wood productivity in Amazonian forests.

Higher levels of taxonomic and evolutionary diversity are expected to maximize ecosystem function, yet their relative importance in driving variation in ecosystem function at large scales in diverse forests is unknown. Using 90 inventory plots across intact, lowland, terra firme, Amazonian forests and a new phylogeny including 526 angiosperm genera, we investigated the association between taxonomic and evolutionary metrics of diversity and two key measures of ecosystem function: aboveground wood productivity and biomass storage. While taxonomic and phylogenetic diversity were not important predictors of variation in biomass, both emerged as independent predictors of wood productivity. Amazon forests that contain greater evolutionary diversity and a higher proportion of rare species have higher productivity. While climatic and edaphic variables are together the strongest predictors of productivity, our results show that the evolutionary diversity of tree species in diverse forest stands also influences productivity. As our models accounted for wood density and tree size, they also suggest that additional, unstudied, evolutionarily correlated traits have significant effects on ecosystem function in tropical forests. Overall, our pan-Amazonian analysis shows that greater phylogenetic diversity translates into higher levels of ecosystem function: tropical forest communities with more distantly related taxa have greater wood productivity.

Evidence for Snake Venom Plasticity in a Long-Term Study with Individual Captive Bothrops atrox.

Variability in snake venom composition has been frequently reported and correlated to the adaptability of snakes to environmental conditions. Previous studies report plasticity for the venom phenotype. However, these observations are not conclusive, as the results were based on pooled venoms, which present high individual variability. Here we tested the hypothesis of plasticity by influence of confinement and single diet type in the venom composition of 13 adult specimens of Bothrops atrox snakes, maintained under captivity for more than three years. Individual variability in venom composition was observed in samples extracted just after the capture of the snakes. However, composition was conserved in venoms periodically extracted from nine specimens, which presented low variability restricted to the less abundant components. In a second group, composed of four snakes, drastic changes were observed in the venom samples extracted at different periods, mostly related to snake venom metalloproteinases (SVMPs), the core function toxins of B. atrox venom, which occurred approximately between 400 and 500 days in captivity. These data show plasticity in the venom phenotype during the lifetime of adult snakes maintained under captive conditions. Causes or functional consequences involved in the phenotype modification require further investigations.

Seasonal and drought-related changes in leaf area profiles depend on height and light environment in an Amazon forest.

Seasonal dynamics in the vertical distribution of leaf area index (LAI) may impact the seasonality of forest productivity in Amazonian forests. However, until recently, fine-scale observations critical to revealing ecological mechanisms underlying these changes have been lacking. To investigate fine-scale variation in leaf area with seasonality and drought we conducted monthly ground-based LiDAR surveys over 4 yr at an Amazon forest site. We analysed temporal changes in vertically structured LAI along axes of both canopy height and light environments. Upper canopy LAI increased during the dry season, whereas lower canopy LAI decreased. The low canopy decrease was driven by highly illuminated leaves of smaller trees in gaps. By contrast, understory LAI increased concurrently with the upper canopy. Hence, tree phenological strategies were stratified by height and light environments. Trends were amplified during a 2015-2016 severe El Niño drought. Leaf area low in the canopy exhibited behaviour consistent with water limitation. Leaf loss from short trees in high light during drought may be associated with strategies to tolerate limited access to deep soil water and stressful leaf environments. Vertically and environmentally structured phenological processes suggest a critical role of canopy structural heterogeneity in seasonal changes in Amazon ecosystem function.

Evidence for snake venom plasticity in a long-term study with individual captive Bothrops atrox

Variability in snake venom composition has been frequently reported and correlated to the adaptability of snakes to environmental conditions. Previous studies report plasticity for the venom phenotype. However, these observations are not conclusive, as the results were based on pooled venoms, which present high individual variability. Here we tested the hypothesis of plasticity by influence of confinement and single diet type in the venom composition of 13 adult specimens of Bothrops atrox snakes, maintained under captivity for more than three years. Individual variability in venom composition was observed in samples extracted just after the capture of the snakes. However, composition was conserved in venoms periodically extracted from nine specimens, which presented low variability restricted to the less abundant components. In a second group, composed of four snakes, drastic changes were observed in the venom samples extracted at different periods, mostly related to snake venom metalloproteinases (SVMPs), the core function toxins of B. atrox venom, which occurred approximately between 400 and 500 days in captivity. These data show plasticity in the venom phenotype during the lifetime of adult snakes maintained under captive conditions. Causes or functional consequences involved in the phenotype modification require further investigations.

Evidence for snake venom plasticity in a long-term study with individual captive Bothrops atrox

Variability in snake venom composition has been frequently reported and correlated to the adaptability of snakes to environmental conditions. Previous studies report plasticity for the venom phenotype. However, these observations are not conclusive, as the results were based on pooled venoms, which present high individual variability. Here we tested the hypothesis of plasticity by influence of confinement and single diet type in the venom composition of 13 adult specimens of Bothrops atrox snakes, maintained under captivity for more than three years. Individual variability in venom composition was observed in samples extracted just after the capture of the snakes. However, composition was conserved in venoms periodically extracted from nine specimens, which presented low variability restricted to the less abundant components. In a second group, composed of four snakes, drastic changes were observed in the venom samples extracted at different periods, mostly related to snake venom metalloproteinases (SVMPs), the core function toxins of B. atrox venom, which occurred approximately between 400 and 500 days in captivity. These data show plasticity in the venom phenotype during the lifetime of adult snakes maintained under captive conditions. Causes or functional consequences involved in the phenotype modification require further investigations.

Isoprene emission structures tropical tree biogeography and community assembly responses to climate.

The prediction of vegetation responses to climate requires a knowledge of how climate-sensitive plant traits mediate not only the responses of individual plants, but also shifts in the species and functional compositions of whole communities. The emission of isoprene gas - a trait shared by one-third of tree species - is known to protect leaf biochemistry under climatic stress. Here, we test the hypothesis that isoprene emission shapes tree species compositions in tropical forests by enhancing the tolerance of emitting trees to heat and drought. Using forest inventory data, we estimated the proportional abundance of isoprene-emitting trees (pIE) at 103 lowland tropical sites. We also quantified the temporal composition shifts in three tropical forests - two natural and one artificial - subjected to either anomalous warming or drought. Across the landscape, pIE increased with site mean annual temperature, but decreased with dry season length. Through time, pIE strongly increased under high temperatures, and moderately increased following drought. Our analysis shows that isoprene emission is a key plant trait determining species responses to climate. For species adapted to seasonal dry periods, isoprene emission may tradeoff with alternative strategies, such as leaf deciduousness. Community selection for isoprene-emitting species is a potential mechanism for enhanced forest resilience to climatic change.

Ecosystem heterogeneity and diversity mitigate Amazon forest resilience to frequent extreme droughts.

The impact of increases in drought frequency on the Amazon forest's composition, structure and functioning remain uncertain. We used a process- and individual-based ecosystem model (ED2) to quantify the forest's vulnerability to increased drought recurrence. We generated meteorologically realistic, drier-than-observed rainfall scenarios for two Amazon forest sites, Paracou (wetter) and Tapajós (drier), to evaluate the impacts of more frequent droughts on forest biomass, structure and composition. The wet site was insensitive to the tested scenarios, whereas at the dry site biomass declined when average rainfall reduction exceeded 15%, due to high mortality of large-sized evergreen trees. Biomass losses persisted when year-long drought recurrence was shorter than 2-7 yr, depending upon soil texture and leaf phenology. From the site-level scenario results, we developed regionally applicable metrics to quantify the Amazon forest's climatological proximity to rainfall regimes likely to cause biomass loss > 20% in 50 yr according to ED2 predictions. Nearly 25% (1.8 million km2 ) of the Amazon forests could experience frequent droughts and biomass loss if mean annual rainfall or interannual variability changed by 2σ. At least 10% of the high-emission climate projections (CMIP5/RCP8.5 models) predict critically dry regimes over 25% of the Amazon forest area by 2100.

Partitioning controls on Amazon forest photosynthesis between environmental and biotic factors at hourly to interannual timescales.

Gross ecosystem productivity (GEP) in tropical forests varies both with the environment and with biotic changes in photosynthetic infrastructure, but our understanding of the relative effects of these factors across timescales is limited. Here, we used a statistical model to partition the variability of seven years of eddy covariance-derived GEP in a central Amazon evergreen forest into two main causes: variation in environmental drivers (solar radiation, diffuse light fraction, and vapor pressure deficit) that interact with model parameters that govern photosynthesis and biotic variation in canopy photosynthetic light-use efficiency associated with changes in the parameters themselves. Our fitted model was able to explain most of the variability in GEP at hourly (R2  = 0.77) to interannual (R2  = 0.80) timescales. At hourly timescales, we found that 75% of observed GEP variability could be attributed to environmental variability. When aggregating GEP to the longer timescales (daily, monthly, and yearly), however, environmental variation explained progressively less GEP variability: At monthly timescales, it explained only 3%, much less than biotic variation in canopy photosynthetic light-use efficiency, which accounted for 63%. These results challenge modeling approaches that assume GEP is primarily controlled by the environment at both short and long timescales. Our approach distinguishing biotic from environmental variability can help to resolve debates about environmental limitations to tropical forest photosynthesis. For example, we found that biotically regulated canopy photosynthetic light-use efficiency (associated with leaf phenology) increased with sunlight during dry seasons (consistent with light but not water limitation of canopy development) but that realized GEP was nonetheless lower relative to its potential efficiency during dry than wet seasons (consistent with water limitation of photosynthesis in given assemblages of leaves). This work highlights the importance of accounting for differential regulation of GEP at different timescales and of identifying the underlying feedbacks and adaptive mechanisms.

Evolutionary heritage influences Amazon tree ecology.

Lineages tend to retain ecological characteristics of their ancestors through time. However, for some traits, selection during evolutionary history may have also played a role in determining trait values. To address the relative importance of these processes requires large-scale quantification of traits and evolutionary relationships among species. The Amazonian tree flora comprises a high diversity of angiosperm lineages and species with widely differing life-history characteristics, providing an excellent system to investigate the combined influences of evolutionary heritage and selection in determining trait variation. We used trait data related to the major axes of life-history variation among tropical trees (e.g. growth and mortality rates) from 577 inventory plots in closed-canopy forest, mapped onto a phylogenetic hypothesis spanning more than 300 genera including all major angiosperm clades to test for evolutionary constraints on traits. We found significant phylogenetic signal (PS) for all traits, consistent with evolutionarily related genera having more similar characteristics than expected by chance. Although there is also evidence for repeated evolution of pioneer and shade tolerant life-history strategies within independent lineages, the existence of significant PS allows clearer predictions of the links between evolutionary diversity, ecosystem function and the response of tropical forests to global change.

Mercury biomagnification and the trophic structure of the ichthyofauna from a remote lake in the Brazilian Amazon.

The present study assesses mercury biomagnification and the trophic structure of the ichthyofauna from the Puruzinho Lake, Brazilian Amazon. In addition to mercury determination, the investigation comprised the calculation of Trophic Magnification Factor (TMF) and Trophic Magnification Slope (TMS), through the measurements of stable isotopes of carbon (δ13C) and nitrogen (δ15N) in fish samples. These assessments were executed in two different scenarios, i.e., considering (1) all fish species or (2) only the resident fish (excluding the migratory species). Bottom litter, superficial sediment and seston were the sources used for generating the trophic position (TP) data used in the calculation of the TMF. Samples from 84 fish were analysed, comprising 13 species, which were categorized into four trophic guilds: iliophagous, planktivorous, omnivorous and piscivorous fish. The δ13C values pointed to the separation of the ichthyofauna into two groups. One group comprised iliophagous and planktivorous species, which are linked to the food chains of phytoplankton and detritus. The other group was composed by omnivorous and piscivorous fish, which are associated to the trophic webs of phytoplankton, bottom litter, detritus, periphyton, as well as to food chains of igapó (blackwater-flooded Amazonian forests). The TP values suggest that the ichthyofauna from the Puruzinho Lake is part of a short food web, with three well-characterized trophic levels. Mercury concentrations and δ13C values point to multiple sources for Hg input and transfer. The similarity in Hg levels and TP values between piscivorous and planktivorous fish suggests a comparable efficiency for the transfer of this metal through pelagic and littoral food chains. Regarding the two abovementioned scenarios, i.e., considering (1) the entire ichthyofauna and (2) only the resident species, the TMF values were 5.25 and 4.49, as well as the TMS values were 0.21 and 0.19, respectively. These findings confirm that Hg biomagnifies through the food web of Puruzinho Lake ichthyofauna. The migratory species did not significantly change mercury biomagnification rate in Puruzinho Lake; however, they may play a relevant role in Hg transport. The biomagnification rate (TMS value) in Puruzinho Lake was higher than the average values for its latitude, being comparable to TMS values of temperate and polar systems (marine and freshwater environments).

Developing Cost-Effective Field Assessments of Carbon Stocks in Human-Modified Tropical Forests.

Across the tropics, there is a growing financial investment in activities that aim to reduce emissions from deforestation and forest degradation, such as REDD+. However, most tropical countries lack on-the-ground capacity to conduct reliable and replicable assessments of forest carbon stocks, undermining their ability to secure long-term carbon finance for forest conservation programs. Clear guidance on how to reduce the monetary and time costs of field assessments of forest carbon can help tropical countries to overcome this capacity gap. Here we provide such guidance for cost-effective one-off field assessments of forest carbon stocks. We sampled a total of eight components from four different carbon pools (i.e. aboveground, dead wood, litter and soil) in 224 study plots distributed across two regions of eastern Amazon. For each component we estimated survey costs, contribution to total forest carbon stocks and sensitivity to disturbance. Sampling costs varied thirty-one-fold between the most expensive component, soil, and the least, leaf litter. Large live stems (≥10 cm DBH), which represented only 15% of the overall sampling costs, was by far the most important component to be assessed, as it stores the largest amount of carbon and is highly sensitive to disturbance. If large stems are not taxonomically identified, costs can be reduced by a further 51%, while incurring an error in aboveground carbon estimates of only 5% in primary forests, but 31% in secondary forests. For rapid assessments, necessary to help prioritize locations for carbon- conservation activities, sampling of stems ≥20cm DBH without taxonomic identification can predict with confidence (R2 = 0.85) whether an area is relatively carbon-rich or carbon-poor-an approach that is 74% cheaper than sampling and identifying all the stems ≥10cm DBH. We use these results to evaluate the reliability of forest carbon stock estimates provided by the IPCC and FAO when applied to human-modified forests, and to highlight areas where cost savings in carbon stock assessments could be most easily made.

Linking canopy leaf area and light environments with tree size distributions to explain Amazon forest demography.

Forest biophysical structure - the arrangement and frequency of leaves and stems - emerges from growth, mortality and space filling dynamics, and may also influence those dynamics by structuring light environments. To investigate this interaction, we developed models that could use LiDAR remote sensing to link leaf area profiles with tree size distributions, comparing models which did not (metabolic scaling theory) and did allow light to influence this link. We found that a light environment-to-structure link was necessary to accurately simulate tree size distributions and canopy structure in two contrasting Amazon forests. Partitioning leaf area profiles into size-class components, we found that demographic rates were related to variation in light absorption, with mortality increasing relative to growth in higher light, consistent with a light environment feedback to size distributions. Combining LiDAR with models linking forest structure and demography offers a high-throughput approach to advance theory and investigate climate-relevant tropical forest change.

Intraspecific isotopic niche variation in broad-snouted caiman (Caiman latirostris).

Strategies to minimise intraspecific competition are common in wild animals. For example, individuals may exploit food resources at different levels of the food chain. Analyses of stable isotopes are particularly useful for confirming variations in an intraspecific niche because the chemical composition of animals tends to reflect both the food consumed and the habitats occupied by the species. However, studies using this methodology to investigate neotropical crocodilians are scarce. This study aimed to verify the existence of ontogenetic and sexual niche variation in broad-snouted caiman in a silvicultural landscape in Brazil through the use of carbon and nitrogen stable isotopes. The isotopic ratios of carbon and nitrogen were determined in claw samples collected from 24 juveniles, 8 adults, and 16 hatchlings of C. latirostris. We identified a discrete ontogenetic variation in the isotopic niche and sexual difference only for juveniles. These results may indicate differences in the exploitation resources and a consequent reduction in competition between age classes.

Influence of land use changes on water chemistry in streams in the State of São Paulo, southeast Brazil.

Streamwater is affected by several processes in the watershed including anthropogenic activities that result in changes in water quality as well as in the functioning of these stream ecosystems. Therefore, this work aims to evaluate the concentration of major ions (Ca(2+), Mg(2+), Na(+), K(+), NH4(+), NO3(-), NO2(-), Cl(-), SO4(2-), PO4(3-), HCO3(-)) in streams in the state of São Paulo (southeast Brazil). The sampling sites are located at undisturbed (ombrophilus dense forest, semideciduous forest and savanna - cerrado) and disturbed areas (pasture, urbanization and sugar cane crops). Streamwater chemistry varied according to land use change and, in general, was higher in disturbed sites. Streams located in undisturbed sites at Ribeira de Iguape/Alto Paranapanema watershed (streams 1, 2 and 3) seem to be regulated by soil characteristics, as the disturbed streams located at the same watershed covered by pasture (stream 7) showed high concentration for the most of the variables. Exception to streams located at Pontal do Paranapanema watershed where both disturbed (stream 8) and undisturbed streams (stream 4 and 5) presented similar patterns for almost all variables measured.