Effects of Small-Scale Gold Mining Tailings on the Underwater Light Field in the Tapaj ó s River Basin , Brazilian Amazon

Artisanal and Small-scale Gold Mining (ASGM) within the Amazon region has created several environmental impacts, such as mercury contamination and changes in water quality due to increased siltation. This paper describes the effects of water siltation on the underwater light environment of rivers under different levels of gold mining activities in the Tapajós River Basin. Furthermore, it investigates possible impacts on the phytoplankton community. Two field campaigns were conducted in the Tapajós River Basin, during high water level and during low water level seasons, to measure Inherent and Apparent Optical Properties (IOPs, AOPs), including scattering (b) and absorption (a) coefficients and biogeochemical data (sediment content, pigments, and phytoplankton quantification). The biogeochemical data was separated into five classes according to the concentration of total suspended solids (TSS) ranging from 1.8 mg·L−1 to 113.6 mg·L−1. The in-water light environment varied among those classes due to a wide range of concentrations of inorganic TSS originated from different levels of mining activities. For tributaries with low or no influence of mining tailings (TSS up to 6.8 mg·L−1), waters are relatively more absorbent with b:a ratio of 0.8 at 440 nm and b660 magnitude of 2.1 m−1. With increased TSS loadings from mining operations (TSS over 100 mg·L−1), the scattering process prevails over absorption (b:a ratio of 10.0 at 440 nm), and b660 increases to 20.8 m−1. Non-impacted tributaries presented a critical depth for phytoplankton productivity of up to 6.0 m with available light evenly distributed throughout the spectra. Whereas for greatly impacted waters, attenuation of light was faster, reducing the critical depth to about 1.7 m, with most of the available light comprising of red wavelengths. Overall, a dominance of diatoms was observed for the upstream rivers, whereas cyanobacteria prevailed in the low section of the Tapajós River. The results suggest that the spatial and temporal distribution of phytoplankton in the Tapajós River Basin is not only a function of light availability, but rather depends on the interplay of factors, including flood pulse, water velocity, nutrient availability, and seasonal variation of incoming irradiance. Ongoing research indicates that the effects of mining tailings on the aquatic environment, described here, are occurring in several rivers within the Amazon River Basin.


Introduction
The increase in the gold price over the past 10 years has stimulated small-scale gold mining in Africa, Asia, and South America, resulting in an annual production of around 400 tonnes of gold by roughly 16 million miners [1].Despite the contribution of approximately 12% of total global gold production [2], artisanal gold mining is known to have many negative environmental impacts [3], and most of these activities are illegal [4].Mercury contamination [2,5,6], geomorphological changes [7], and water siltation are among the primary direct impacts of small-scale gold mining [8].Water siltation caused by gold mining is commonly reported throughout the world (such as Indonesia [2] in Asia, Ghana [9] in Africa, and French Guyana [10] and Peru [11] in South America) due to the fact that most of the mining activities occur in rivers or at their margins [12].In the Brazilian Amazon, for example, sediments from mining tailings in streams and rivers may vary between one and two tonnes per gram of gold produced [13].One of the major impacts of water siltation is the increase in light attenuation by suspended particles, which can directly affect phytoplankton productivity by limiting the availability of Photosynthetic Active Radiation (PAR), and indirectly contribute to changes in biodiversity within rivers, including plankton [10], macroinvertebrates [14], and fish [15].Specific to the Amazon Basin, a major reduction in phytoplankton density has been documented in Batata Lake due to a gain in suspended matter caused by mining activity [16].Guenther and Bozelli [17] have attributed this reduction to a decrease in growth rates caused by high light attenuation due to increased suspended matter derived from mining activities.Moreover, Tudesque et al. [10] have reported changes in phytoplankton community with a rise in the proportion of benthic diatoms due to increase in water siltation in streams impacted by small-scale gold mining in French-Guiana.
While a few studies have shown the impacts on underwater light conditions caused by mining tailings [16,17], the quantification of the impacts on the spectral underwater light environment in Amazonian waters is still lacking.Understanding the underwater light environment, more specifically the total scalar irradiance, E o , is fundamentally important since it comprises the total spectral energy available for photosynthesis in the interval between 400 and 700 nm, which in turn allows for comprehension of the phytoplankton absorption efficiency [18].For example, an underwater light field rich in blue-green light will favor phytoplankton groups that have pigments that absorb light in the blue-green spectra [19].
This research focuses on defining the impact of gold mining tailings in the light environment of the Tapajós River Basin, and associated phytoplankton communities' diversity.Specifically, this study provides the following: (1) an evaluation of the effects of sediment loading on the underwater light environment of impacted and non-impacted rivers; and (2) evaluation of the relationship between light attenuation caused by water siltation and the critical depth for phytoplankton communities.To address these objectives, biogeochemical and in situ optical measurements, including Inherent Optical Properties (IOPs) and Apparent Optical Properties (AOPs), were collected during two field campaigns addressing high and low water conditions.Despite the strong emphasis on underwater light fields, this paper also provides a partial discussion on the distribution of phytoplankton communities and other important factors, such as hydrological regime, current velocity, and depth.

Study Area
The Tapajós River Basin includes the largest small-scale gold mining area in the world (28,186 km 2 ), which has been intensively mined, either using water-jets to remove top soil layers, or using small boats called 'balsas' that take the sediment from the bottom of the rivers using suction and gravity processes [20,21].Currently, more than 300 small-scale mines (most without mining licenses) with the participation of more than 50,000 miners produce gold within the Artisanal Gold Mining District.This district was created in 1983 to support gold miners and overlaps three main sub-basins: Novo, Crepori, and Tocantinzinho (abbreviated Tocantins) (Figure 1).The lower section of the Tapajós River Basin covers about 130,370 km 2 (Figure 1) and drains mostly lixiviated Pre-Cambrian rocks, which results in water that is transparent/greenish with low amounts of suspended solids, so-called "clearwaters" [24,25].The river basin can generally be The lower section of the Tapajós River Basin covers about 130,370 km 2 (Figure 1) and drains mostly lixiviated Pre-Cambrian rocks, which results in water that is transparent/greenish with low amounts of suspended solids, so-called "clearwaters" [24,25].The river basin can generally be separated into two geomorphological sections: the upstream riverine section (lotic system) from headwaters down to the Aveiro City region; and the downstream section (semi-lentic system-locally called Tapajós Lake), from Aveiro City to the mouth of the river where it runs into the Amazon River (Figure 1).
The downstream river section is relatively wide (10-15 km) with low water speed, creating a semi-lacustrine environment, able to support phytoplankton growth.Chlorophyll-a (chl-a) concentrations vary from 5.0 to 25.1 µg•L −1 [26,27], and cyanobacteria blooms related to urban-derived diffuse nutrient input at the margins of the Tapajós mouth (Santarém) have been reported [28].In this section of the river, total suspended sediments (TSS) are composed primarily of fine matter discharged from upstream [29][30][31] reaching concentrations up to 4.1 mg•L −1 during rising and receding water levels [32].In the Tapajós, for example Costa et al. [32] observed, on average, lower a cdom absorption (a cdom (440) = 3.2 m −1 ) than values measured in black water rivers of the Amazon (8.6 m −1 ).As a result of the low TSS and a cdom , the diffuse attenuation coefficient (K d PAR) was reported at around 2.0 m −1 in both rising and falling periods, whereas for black and white Amazonian rivers, K d increases to 3.7 and 5.1 m −1 , respectively [32].
Information about biogeochemistry [21,33] and hydrological optics of the upstream section of the river is rare.As opposed to the semi-lentic conditions of the Tapajós Lake section, the upstream section starting from Aveiro City (Figure 1) is narrower (2-4 km), with a strong advection current [30].In such systems, plankton development is likely suppressed by the low residence time and turbulence of running water [34].As a result of low photosynthetic production, the suspended sediment content is composed mostly of allochthonous matter, such as quartz, plant debris, and clay clusters, rather than autochthonous matter [33,35].

Methods
The methodology of this study consisted of two components (Figure 2): (1) acquisition of in situ optical (Table 1) and biogeochemical data at high (2011) and low (2012) water levels to define the in-water light environment; and (2) evaluation of the critical depth for phytoplankton productivity based on the available light.separated into two geomorphological sections: the upstream riverine section (lotic system) from headwaters down to the Aveiro City region; and the downstream section (semi-lentic system-locally called Tapajós Lake), from Aveiro City to the mouth of the river where it runs into the Amazon River (Figure 1).
The downstream river section is relatively wide (10-15 km) with low water speed, creating a semi-lacustrine environment, able to support phytoplankton growth.Chlorophyll-a (chl-a) concentrations vary from 5.0 to 25.1 μg•L −1 [26,27], and cyanobacteria blooms related to urban-derived diffuse nutrient input at the margins of the Tapajós mouth (Santarém) have been reported [28].In this section of the river, total suspended sediments (TSS) are composed primarily of fine matter discharged from upstream [29][30][31] reaching concentrations up to 4.1 mg•L −1 during rising and receding water levels [32].In the Tapajós, for example Costa et al. [32] observed, on average, lower acdom absorption (acdom (440) = 3.2 m −1 ) than values measured in black water rivers of the Amazon (8.6 m −1 ).As a result of the low TSS and acdom, the diffuse attenuation coefficient (Kd PAR) was reported at around 2.0 m −1 in both rising and falling periods, whereas for black and white Amazonian rivers, Kd increases to 3.7 and 5.1 m −1 , respectively [32].
Information about biogeochemistry [21,33] and hydrological optics of the upstream section of the river is rare.As opposed to the semi-lentic conditions of the Tapajós Lake section, the upstream section starting from Aveiro City (Figure 1) is narrower (2-4 km), with a strong advection current [30].In such systems, plankton development is likely suppressed by the low residence time and turbulence of running water [34].As a result of low photosynthetic production, the suspended sediment content is composed mostly of allochthonous matter, such as quartz, plant debris, and clay clusters, rather than autochthonous matter [33,35].

Methods
The methodology of this study consisted of two components (Figure 2): (1) acquisition of in situ optical (Table 1) and biogeochemical data at high (2011) and low (2012) water levels to define the inwater light environment; and (2) evaluation of the critical depth for phytoplankton productivity based on the available light.

Inherent Optical Property
Symbol Unit Formula

Apparent Optical Property
Upwelling radiance  1 for sample site locations).Sample sites were distributed along the Tapajós River and main tributaries, which correspond to known locations of mined (Crepori, Tocantins, and Novo) and non-mined tributaries (Jamanxim River).As a side note, logistically, fieldwork in this remote region of the Amazon was extremely difficult due to the size and weight of the optical sensors and other field and analytical equipment, as well as issues related to security due to local land conflicts.

Biogeochemical Data
For each sampled site (n = 40), two water samples (at least 500 mL sample volume) were taken at a depth of 0.3 m and at Z 1% (depth where downwelling irradiance reaches 1% of surface irradiance) to determine TSS concentrations according to the gravimetric method [36].The surface (0.3 m) and depths at Z 1% were chosen to represent the water column where phytoplankton productivity can occur [34].For each water sample, triplicates of pre-weighed GF/F (0.7 µm) filters were preserved in cold, dark environments until laboratory analysis was conducted.After TSS analysis, half of those filters were used for particulate organic carbon (POC) determination according to the high-temperature combustion method with TOC-V analyzer (Shimadzu Inc., Guelph, ON, Canada) [37].The other filters were used in determining the percentage of the organic content present by applying the Loss On Ignition (LOI) technique [36].Water samples (500 mL) were also filtered with GF/F (0.7 µm) for determination of chl-a and auxiliary pigments by HPLC (High-Performance Liquid Chromatography) [36].The pigments on the filters were extracted with 90% acetone, centrifuged, and analyzed in HPLC (Dionex P680) within 24 h after extraction to avoid pigment degradation [38].
Discrete surface water samples were taken from a subset (n = 15 including both campaigns) of the total sampled locations in mined and non-mined tributaries for quantitative and qualitative analyses of phytoplankton populations using the microscopy settling technique [39].The units, corresponding to cells, colonies, and filaments, were enumerated to at least 100 specimens of the most frequent species at 400× [40].Specific phytoplankton biovolume (mm 3 L −1 ) was estimated according to Hillebrand et al. [41].

Optical Data
Inherent optical properties (IOP) measurements were only acquired during the high water level period (April 2011, n = 27) due to logistical problems (previously explained).Total beam attenuation c(λ) and total absorption a(λ) coefficients were measured in situ with a WetLabs (Philomath, Oregon) ac-S instrument at 80 wavelengths from 390 to 750 nm.Total beam attenuation coefficient c(λ) is the sum of the total absorption a(λ) and total scattering coefficient b(λ).The absorption coefficient is defined as the sum of the absorption coefficients by water a w (λ), particulate material a p (λ), and coloured dissolved organic matter a cdom (λ) [42].Similarly, the total scattering is the sum of its component scattering coefficient of water b w (λ) and particulate material b p (λ), with the general assumption that scattering due to CDOM is negligible [42].
The ac-S output was calibrated using ultra-clean water from the Barnstead E-pure water purification system before and after field campaigns to minimize measurement deviation caused by sensor transportation.The output was also corrected for temperature and scattering in the absorption tube using the manufacturer protocol [43].The a cdom (λ) coefficient over PAR was acquired after filtering sampled water with a Pall 0.2 µm sterile filter.Spectral measurements of the filtered water were taken in situ with the Ocean Optics UV-4000 equipment [44].The measured absorbance was firstly corrected for blank offset measured with deionized water and secondly converted to a cdom (λ) according to Kirk [18].Attached in the same cage along with ac-S, the Environmental Characterization Optics-Backscattering (ECO-BB) by Wetlabs© measured particle backscattering (b bp (λ)) at three wavelengths (470, 532, and 660 nm).For each wavelength, b bp (λ) was determined as follows: where b bp (λ) is the particulate backscattering coefficient, β b (117 • , λ) is the volume scattering function given a wavelength (λ) for a specific angle (117 • ), and x is 1.1 as determined by Boss et al. [45].

Critical Depth for Photosynthesis
Given that phytoplankton cells can utilize irradiance from all directions, the quantity and quality (spectral distribution) of the light available for photosynthesis in the aquatic environment, underwater total scalar irradiance, E o (0 − , λ) (λ = 400-700 nm), is defined as follows [18]: where E u (0 − , λ) and E d (0 − , λ) are upwelling and downwelling irradiance, respectively.
In order to compare the light availability amongst the five classes and minimize the effect of variations in atmospheric conditions during each measurement, scalar irradiance (E o ) was normalized by downward irradiance measured above water (E d + ) resulting in E o , as follows: Once the E o (0 − , 400-700 nm) was calculated, the critical depth for phytoplankton was estimated for the different classes of water.Critical depth, (Z c ), is the depth at which the photosynthetic production by a phytoplankton community is equal to the energy loss on processes such as respiration.In other words, given optimum conditions of nutrient concentration, Z c is the depth at which irradiance is the minimum necessary for photosynthesis.Z c is calculated as follows [18]: where E o (0 − , PAR) is the sub-surface scalar irradiance modeled by Hydrolight; K o (PAR) is the average scalar attenuation coefficient for PAR from the surface to Z 1% (similar to Equation ( 6)), and E c (PAR) is the species-specific compensation irradiance.To calculate Z c , we adopted E c (PAR) based on the freshwater phytoplankton minimum light requirement described by Deblois et al. [47].The authors reported that Chlamydomonas sp.(Chlorophyta) (identified in this study) shows a growth rate close to zero when exposed to irradiance of 14 µEm −2 s −1 , whereas other species, such as Aulacoseira granulate (diatom) and Cryptomonas obovata (Cryptophyta), present a growth rate of up to 0.2 divisions/day when exposed to the same light regime.For freshwater phytoplankton communities, a minimum E c of 14 µEm −2 s −1 was defined as a general compensation irradiance as described in Deblois et al. [47] and therefore used as input to calculate Z c in Equation ( 9).Furthermore, analyses of critical depth must take into account the depth of the mixing layer (Z m ).If the depth of the mixed layer is smaller than or equal to Z c , then there is sufficient light to drive net production and phytoplankton growth [48].On the other hand, if Z m is greater than Z c , light can become a limiting factor for net production in deeper layers.In river systems, where high current velocity and turbulence mixes the water column homogeneously, Z m usually coincides with the total water depth (Z) [49].In turbid and deep rivers such as the Amazon, phytoplankton is only exposed to light for a short period within the upper layers.This low light exposure can reduce phytoplankton productivity, as indicated by Dustan [50].For oceans and lakes, thermal stratification can occur, reducing the mixed layer to shallower depths than Z c and favoring phytoplankton growth and bloom conditions [49].

Data Analysis
The analyses of collected data were conducted in four steps: (1) Sampled points were organized into five classes of TSS concentration, corresponding to different levels of mining influence from low to high mining tailing influence; (2) TSS classes were submitted to non-parametric test to verify the significance difference between them at p < 0.10; (3) Once the TSS classes were defined and tested, a full optical and biogeochemical characterization were conducted including phytoplankton, pigments, IOPs, and AOPs; and (4) In order to verify the contribution of discharged sediment derived from mining activities to the optical properties (IOPs and AOPs), Pearson's correlations between sediment content and optical properties were performed (p < 0.05).

Biogeochemical Data
The biogeochemical data were separated into five preliminary classes according to the concentration of total suspended matter aiming to facilitate the interpretation of the results (Table 2).The definition of classes was based on prior knowledge of the level of mining impact on the sampled rivers and statistical tests for significant differences on TSS (non-parametric Mann-Whitney U Tests; p < 0.10, Appendix A).The non-parametric test confirmed the significant difference among the five classes of TSS (Table 2), which displays considerable changes in their spatial and seasonal distribution along the Tapajós River and main tributaries (Figure 3).classes of TSS (Table 2), which displays considerable changes in their spatial and seasonal distribution along the Tapajós River and main tributaries (Figure 3).Class 1 includes upstream Jamanxim (Jamanxim 1), Tapajós-Upstream, and the Tapajós Lake section of the river (Figure 3), corresponding to a total of 17 samples with TSS concentrations of up to 5.0 mg•L −1 .Although all samples belong to clear water [25,52], their hydrodynamic characteristics vary from lotic upriver to a semi-lentic system at the Tapajós Lake region.For that reason, they will be presented and discussed separately.
Sampled during high water period, upstream Jamanxim River site showed acdom (440) of 3.2 m −1 and maximum values of TSS of 3.6 mg•L −1 with 37% of organic content (POC = 17%).Tapajós-Upstream (Figure 3) sampled during low water period showed similar suspended sediment and organic content concentrations as those of upstream Jamanxim, but higher average chl-a concentration (1.1 μg•L −1 ).At the lower section of the Tapajós Lake, TSS values of 2.0 mg•L −1 were measured during the low water level period and approximately 3.0 mg•L −1 in the high water period.The organic content in the Tapajós Lake was consistently higher than that in upstream tributaries and Tapajós-Upstream during both seasons, with an average of 50% (POC ~ 13%).Differences were not observed for acdom (440).chl-a followed values observed at upstream Tapajós River, with concentrations up to 1.4 μg•L −1 .
Class 2 includes 15 sites with low impact mining activities, in which TSS concentrations vary from 5.0 to 10.0 mg•L −1 , and comprises water samples taken at the Jamanxim River, after the confluence of Novo (Jamanxim 2) and Tocantins rivers (Jamanxim 3), acquired during the high water period, and samples of the Tapajós River at Itaituba (Tapajós-Itaituba), sampled during both periods (Figure 3).TSS was relatively higher (6.3 mg•L −1 on average) than in Class 1 (3.3 mg•L −1 ), followed by decreased organic and POC content to 24% and 9%, respectively.On the other hand, acdom (440) (3.2 m −1 ) values were similar for both classes.The average chl-a concentration within Class 2, however, was higher (1.1 μg•L −1 ) at the Tapajós-Itaituba river section than at the tributaries (0.1 μg•L −1 in average).Class 1 includes upstream Jamanxim (Jamanxim 1), Tapajós-Upstream, and the Tapajós Lake section of the river (Figure 3), corresponding to a total of 17 samples with TSS concentrations of up to 5.0 mg•L −1 .Although all samples belong to clear water [25,52], their hydrodynamic characteristics vary from lotic upriver to a semi-lentic system at the Tapajós Lake region.For that reason, they will be presented and discussed separately.
Sampled during high water period, upstream Jamanxim River site showed a cdom (440) of 3.2 m −1 and maximum values of TSS of 3.6 mg•L −1 with 37% of organic content (POC = 17%).Tapajós-Upstream (Figure 3) sampled during low water period showed similar suspended sediment and organic content concentrations as those of upstream Jamanxim, but higher average chl-a concentration (1.1 µg•L −1 ).At the lower section of the Tapajós Lake, TSS values of 2.0 mg•L −1 were measured during the low water level period and approximately 3.0 mg•L −1 in the high water period.The organic content in the Tapajós Lake was consistently higher than that in upstream tributaries and Tapajós-Upstream during both seasons, with an average of 50% (POC ~13%).Differences were not observed for a cdom (440) .chl-a followed values observed at upstream Tapajós River, with concentrations up to 1.4 µg•L −1 .
Class 2 includes 15 sites with low impact mining activities, in which TSS concentrations vary from 5.0 to 10.0 mg•L −1 , and comprises water samples taken at the Jamanxim River, after the confluence of Novo (Jamanxim 2) and Tocantins rivers (Jamanxim 3), acquired during the high water period, and samples of the Tapajós River at Itaituba (Tapajós-Itaituba), sampled during both periods (Figure 3).TSS was relatively higher (6.3 mg•L −1 on average) than in Class 1 (3.3 mg•L −1 ), followed by decreased organic and POC content to 24% and 9%, respectively.On the other hand, a cdom (440) (3.2 m −1 ) values were similar for both classes.The average chl-a concentration within Class 2, however, was higher (1.1 µg•L −1 ) at the Tapajós-Itaituba river section than at the tributaries (0.1 µg•L −1 in average).Class 3 includes tributaries subjected to gold mining at medium impact level with TSS concentration between 10.0 and 20.0 mg•L −1 , and is represented by samples taken in the Tocantins River during the high water period, and in the Tapajós River after receiving the Crepori River inflow (named Tapajós-Crepori) during low water level periods.Overall, a cdom (440) were similar to those from previous classes, and organic matter decreased to 13% (POC = 4%), indicating an increased influence of inorganic sediment input from mining tailings to the TSS when compared to Classes 1 and 2.
Class 4 and Class 5 represent rivers subjected to high level of mining activity and comprise samples of the Crepori River at the high and low water periods, respectively (Figure 3).During the high water level period, TSS concentration of approximately 33.0 mg•L −1 was observed, whereas in the low water period, TSS increased to 111.0 mg•L −1 .Values for a cdom (440) , however, were similar in both periods (2.1 and 2.5 m −1 , respectively).The organic matter content for Class 4 (12%) was not different from that observed in Class 3, but it decreased to 2% for Class 5 due mainly to POC (Table 2).Regarding chl-a, Crepori River presented relatively higher chl-a (0.7 µg•L −1 ) concentration during low water level when compared to samples acquired during high water level (0.1 µg•L −1 ).

Phytoplankton and Pigments
Although samples were limited, phytoplankton microscopy analysis indicated a higher biovolume (mm 3 •L −1 ) during low water period compared to the high water period (Figure 4).Further, the upstream tributaries showed lower phytoplankton biovolume and dominance of diatoms and cryptomonads when compared to the Tapajós Lake section where most of the phytoplankton biovolume was higher (up to 6 mm 3 •L −1 ) and dominated by cyanobacteria cells (Figure 4).
During the low water level, total pigment concentration varied between 1.0 and 2.0 µg•L −1 for all samples, with no clear difference in pigment contents and concentrations between the Tapajós Lake and the main channel and tributaries (e.g., Crepori River).The total pigment concentration in the Crepori River was higher during low water level period than at high water period, similar to the biovolume data.Overall, the relative contribution of auxiliary pigments remained consistent throughout the samples, showing chl-b (14%), zeaxanthin (14%), and fucoxanthin (4%) as the main contributors, thus indicating the presence of chlorophyte, cyanobacteria, and diatoms, confirmed with microscopy analysis.

Inherent Optical Properties (IOPs)
The absorption spectra by coloured dissolved organic matter (a cdom ) has a typical power-law function, decreasing from blue towards red wavelengths (Figure 5a), and similar a cdom (440) magnitudes (~2.7 m −1 ) among the different classes of water, except for Class 4, which showed slightly lower a cdom (440) (2.0 m −1 ) (Table 2).

Inherent Optical Properties (IOPs)
The absorption spectra by coloured dissolved organic matter (acdom) has a typical power-law function, decreasing from blue towards red wavelengths (Figure 5a), and similar acdom (440) magnitudes (~2.7 m −1 ) among the different classes of water, except for Class 4, which showed slightly lower acdom (440) (2.0 m −1 ) (Table 2).Generally, ap, bp, and bbp showed a larger variability than acdom, and similar to the TSS trend (i.e., the IOPs magnitude increased from Class 1 to Class 4).As a result, TSS was significantly correlated (p < 0.05) to bp and bbp, with correlation coefficients as high as 0.99 for both (Table 3).5b), and bp (660) from 2.7 for Class 1 to 20.8 m −1 for Class 4 (Figure 5c).Further, the contribution of bp (660) to total attenuation varied from 34% for Class 1 to >80% for Classes 4 and 5 (Figure 5d).The backscattering coefficient, bbp (660), also exhibited a similar trend, increasing from 0.10 m −1 for Class 1 to 0.60 m −1 for Class 4 (Table 2).Generally, a p , b p , and b bp showed a larger variability than a cdom , and similar to the TSS trend (i.e., the IOPs magnitude increased from Class 1 to Class 4).As a result, TSS was significantly correlated (p < 0.05) to b p and b bp , with correlation coefficients as high as 0.99 for both (Table 3).5b), and b p (660) from 2.7 for Class 1 to 20.8 m −1 for Class 4 (Figure 5c).Further, the contribution of b p (660) to total attenuation varied from 34% for Class 1 to >80% for Classes 4 and 5 (Figure 5d).The backscattering coefficient, b bp (660) , also exhibited a similar trend, increasing from 0.10 m −1 for Class 1 to 0.60 m −1 for Class 4 (Table 2).

Apparent Optical Properties (AOPs) and Underwater Light Field
The spectra for diffuse attenuation of downwelling irradiance (K d ) was higher at blue wavelengths and decreased towards green and red wavelengths (Figure 6), similar to the IOPs spectra.Considering the PAR range, Class 1 presented a K d (PAR) of 1.9 m −1 , and Class 5 values were as high as 9.6 m −1 .Not surprisingly, K d (PAR) showed a strong linear correlation with TSS for all wavelengths, for example, at 660 nm (r = 0.96), 560 nm (r = 0.94), and 440 nm (r = 0.84, p < 0.05, n = 40).Again, similar correlations were not observed for chl-a and a cdom (440) data (Table 4), thus indicating the prevalent role of inorganic suspended particulates in the AOPs.

Apparent Optical Properties (AOPs) and Underwater Light Field
The spectra for diffuse attenuation of downwelling irradiance (Kd) was higher at blue wavelengths and decreased towards green and red wavelengths (Figure 6), similar to the IOPs spectra.Considering the PAR range, Class 1 presented a Kd (PAR) of 1.9 m −1 , and Class 5 values were as high as 9.6 m −1 .Not surprisingly, Kd (PAR) showed a strong linear correlation with TSS for all wavelengths, for example, at 660 nm (r = 0.96), 560 nm (r = 0.94), and 440 nm (r = 0.84, p < 0.05, n = 40).Again, similar correlations were not observed for chl-a and acdom (440) data (Table 4), thus indicating the prevalent role of inorganic suspended particulates in the AOPs.Overall, (0.3m, λ) showed an average magnitude around 1.0 with a spectral shape increasing from blue to red wavelengths (Figure 7a).Values higher than 1.0 can occur when the upwelling irradiance, Eu (0.3m, λ), contributes significantly to the total scalar irradiance [54].At sub-surface, (0.3m, λ) for Class 1 was about 0.5 in the blue region, and approximately 1.0 at green and red wavelengths.In more turbid waters (Classes 4 and 5), short wavelengths decreased quickly at sub-surface.Whereas, green and mainly red wavelengths showed (0.3m) values up to 1.7 of incoming light, thus comprising most of the (0.3m, PAR) at this depth (Table 2).At a depth of 2.0 m, light intensity for all classes is mostly concentrated within the green and red wavelengths; blue light is fully attenuated, even for Class 1 (Figure 7).At this depth, Class 5 presented only 0.1% of the incoming PAR light, increasing towards Class 4 (0.4%), Class 3 (10.1%),Class 2 (13.5%), and Class 1 (18.0%).
The depth in which (PAR) is 1% (Z1%) varied from 6.0 m for Class 1 to average of 5.4, 4.3, 3.6, 2.8, and 1.7 m for Classes 2 to 5, respectively.The wavelength in which light penetrates the most at Z1% shifted from 620 nm (Classes 1 and 2) to 660 nm in turbid waters (Classes 4 and 5) (Figure 7c).
The critical depth, Zc, for Chlamydomonas sp., used here is representative of the species in our samples, indicating that for water of Class 1, enough light was available for photosynthesis up to 6.0 m depth (Figure 7).As water turbidity increased, a decrease in Zc was observed (Table 2).Note that  Overall, E o (0.3m, λ) showed an average magnitude around 1.0 with a spectral shape increasing from blue to red wavelengths (Figure 7a).Values higher than 1.0 can occur when the upwelling irradiance, E u (0.3m, λ) , contributes significantly to the total scalar irradiance [54].At sub-surface, E o (0.3m, λ) for Class 1 was about 0.5 in the blue region, and approximately 1.0 at green and red wavelengths.In more turbid waters (Classes 4 and 5), short wavelengths decreased quickly at sub-surface.Whereas, green and mainly red wavelengths showed E o (0.3m) values up to 1.7 of incoming light, thus comprising most of the E o (0.3m, PAR) at this depth (Table 2).At a depth of 2.0 m, light intensity for all classes is mostly concentrated within the green and red wavelengths; blue light is fully attenuated, even for Class 1 (Figure 7).At this depth, Class 5 presented only 0.1% of the incoming PAR light, increasing towards Class 4 (0.4%), Class 3 (10.1%),Class 2 (13.5%), and Class 1 (18.0%).
The depth in which E o (PAR) is 1% (Z 1% ) varied from 6.0 m for Class 1 to average of 5.4, 4.3, 3.6, 2.8, and 1.7 m for Classes 2 to 5, respectively.The wavelength in which light penetrates the most at Z 1% shifted from 620 nm (Classes 1 and 2) to 660 nm in turbid waters (Classes 4 and 5) (Figure 7c).
The critical depth, Z c , for Chlamydomonas sp., used here is representative of the species in our samples, indicating that for water of Class 1, enough light was available for photosynthesis up to 6.0 m depth (Figure 7).As water turbidity increased, a decrease in Z c was observed (Table 2).Note that within the upper 0.5 m, the scalar irradiance available can be higher than the E max , which is the maximum level of irradiance that yields photosynthesis.
Remote Sens. 2017, 9, 861 14 of 22 within the upper 0.5 m, the scalar irradiance available can be higher than the Emax, which is the maximum level of irradiance that yields photosynthesis.

Discussion
Water siltation caused by ASGM in the Tapajós River Basin over the past 40 years has been demonstrated by Lobo et al. [51,55], using a historical Landsat imagery archive.From 1978 to 1993, the introduction of low-budget mechanization associated with very high gold prices resulted in the large expansion of mining area and high water siltation levels.After a general recession during the 1990s and early 2000s, in the past 12 years, ASGM has intensified encouraged by high gold prices, resulting in expansion of mining sites and increase of TSS input to the rivers.The current gold rush in this region raises the concern about the aquatic impacts caused by the discharge of mining tailings into the rivers.
To investigate this issue, the present study provides a full characterization of the optical changes due to increased TSS from mining tailings for the Tapajós River and tributaries.Among the main findings, measured optical properties indicated that the incoming light (PAR) is rapidly attenuated with increased TSS concentrations caused by ASGM activity.Non-impacted tributaries presented a critical depth for phytoplankton productivity of up to 6.0 m, while for greatly impacted waters, attenuation of light was rapid, reducing Zc to about 1.7 m, with most of the available light comprising red wavelengths.

Discussion
Water siltation caused by ASGM in the Tapajós River Basin over the past 40 years has been demonstrated by Lobo et al. [51,55], using a historical Landsat imagery archive.From 1978 to 1993, the introduction of low-budget mechanization associated with very high gold prices resulted in the large expansion of mining area and high water siltation levels.After a general recession during the 1990s and early 2000s, in the past 12 years, ASGM has intensified encouraged by high gold prices, resulting in expansion of mining sites and increase of TSS input to the rivers.The current gold rush in this region raises the concern about the aquatic impacts caused by the discharge of mining tailings into the rivers.
To investigate this issue, the present study provides a full characterization of the optical changes due to increased TSS from mining tailings for the Tapajós River and tributaries.Among the main findings, measured optical properties indicated that the incoming light (PAR) is rapidly attenuated with increased TSS concentrations caused by ASGM activity.Non-impacted tributaries presented a critical depth for phytoplankton productivity of up to 6.0 m, while for greatly impacted waters, attenuation of light was rapid, reducing Z c to about 1.7 m, with most of the available light comprising red wavelengths.
4.1.Mining-Derived TSS as the Main Factor Changing the Water Optical Properties and Light Field Amongst the sampled rivers, the upstream section of the Tapajós River is naturally classified as "clearwater" [21,24,25].This class presented relatively low TSS (~5.0 mg•L −1 ), low dissolved organic matter (a cdom < 2.5 m −1 ), and low chl-a (<1.0 µg•L −1 thus resulting in a relatively deep euphotic zone (Z 1% ~6.0 m).In these waters, the suspended sediment has a considerable amount of organic matter (~30% of TSS), composed mostly of allochthonous plant debris [33].The characteristics and concentration of the suspended sediments change abruptly as the Tapajós River receives clay-rich tributaries, such as the heavily-mined Crepori River (TSS ~111.3 mg•L −1 , particulate organic matter <3%, and euphotic depth ~2.0 m; defined as Class 5).The sediment plume from the Crepori River only fully mixes with the Tapajós River waters at about 200 km downstream after passing through rapids [21,51].After these rapids, as the water velocity decreases, the fine suspended solids sink, and concentrations decrease to values similar to those of the upstream Tapajós River (see Telmer et al. [21]).Similarly to the Tapajós River, TSS at the Jamanxim River (Class 1 and Class 2) increases as it receives a sediment-rich discharge from the Novo and the Tocantins sub-basins (Class 3) subject to mining operations (Figure 3).
Seasonally, from low to high water level period, the increase of TSS is, in part, a result of the sediment input carried during rain events as the water rises [24,56].Gold mining activities are also temporally dynamic given that they intensify during the dry season [57] when access to mining pits are facilitated.This intensification combined with the lower river levels during that season, increases sediment concentrations to above 100 mg•L −1 at Crepori River (Figure 3), the most ASGM-impacted basin.Recent results of an ongoing research [58] assessing the proportion sediment production from diffused soil erosion to that of ASGM showed that, approximately, 10% of sediment concentration measured at Crepori River (above 100 mg•L −1 ) could be attributed to diffuse soil erosion throughout the year.Therefore, the remaining 90% of the TSS concentration observed can be attributed to mining tailings.
In the Tapajós region, the mining tailings flushed into the rivers are composed of an agglomeration of fine clay particles [21,59], mostly kaolinite [23].Specifically, the fine inorganic clay particles are more effective at scattering and backscattering light due to their higher refractive index compared to large inorganic particulates (silt, medium sand), organic particulates such as flocs (agglomeration of particles), and phytoplankton [60][61][62].These mine-derived clay particles play a strong role in the magnitude and spectral dependency of IOPs and AOPs compared with waters under less or no influence of mined rivers (see correlations on Tables 3 and 4).
The absorption curves for particulate matter showed typical power-law function [42]; the typical chl-a absorption peak at 675 nm, however, was not clearly observed, indicating that suspended matter, mostly of inorganic nature, and dissolved organic matter are the main optical components absorbing light in the non-impacted (Class 1) and impacted tributaries (Classes 3-5), respectively.
Measured b 660 magnitudes increased one order of magnitude from Class 1 (2.1 m −1 ) to Class 4 (20.8 m −1 ), and are comparable with other scattering coefficients reported for turbid inland waters.For example, Sun et al. [63] reported b 660 acquired in Taihu Lake, a turbid lake in China, up to 22.0 m −1 with increasing inorganic concentration (around 33.0 mg•L −1 , which corresponds to Class 4 in our study, for which b 660 is on average 20.8 m −1 ).Similarly to scattering, b b660 showed high values (0.65 m −1 for Class 3, for example), which are comparable to the average backscattering (b b700 ) of 0.80 m −1 reported for waters where TSS reached 18.9 mg•L −1 in the Curuai Lake, adjacent to the Tapajós 'Lake' [64].
Both scattering and backscattering coefficients are primarily controlled by mineral concentration (TSS) as reported by several authors [45,63,[65][66][67][68]. Increasing TSS concentration yielded positive and significant (p <0.05) correlation with b p (660) , b bp (660) , and a p (440) .Correlations of optical properties with chl-a and a cdom (440) were not as significant (see correlations on Tables 3 and 4), what is explained by the inorganic nature of the soils subject to mining in the Tapajos watershed [7,23].
The IOP of waters with high concentrations of fine sediment results in changes in the K d and consequently in the underwater light field.This is represented by the defined five water classes representing increasing TSS (Figure 7).Measured K d (PAR) values varied from 2.0 m −1 for Class 1 to 9.6 m −1 for Class 5 (Table 2), which is similar to results reported for other mine-impacted areas in New Zealand, where K d (PAR) varied on average from 1.6 to 8.4 m −1 as averaged TSS increased from 1.2 to >100.0 mg•L −1 [69].Other studies, specifically in the mine tailings of the Batata Lake, located in the Trombetas River watershed, showed that re-suspended inorganic particles derived from mining increased K d (PAR) by 35% when compared to non-impacted sections of the lake [16,17].Spectrally, the attenuation coefficient showed a similar trend; K d (440) , K d (560) , K d (660) , varying from 2.0, 1.5, 1.2, and 1.0 (Class 1) to 13.1, 8.0, and 7.0 m −1 (Class 5), respectively.These K d (440) magnitudes are comparable to values reported by Costa et al. [32] for different waters in the Amazon Basin (K d (440) up to 3.0 m −1 in clear and 15.1 m −1 in turbid waters).The authors also reported K d (660) values lower than those of the blue spectra, but still higher in white waters (9.5 m −1 , corresponding to Class 5) in comparison to clear waters (up to 1.7 m −1 , corresponding to Class 1).
The distinct spectral attenuation of the incoming light is a function of the variation in the scattering and absorption processes within the PAR range (400-700 nm).For tributaries with low or no influence of mining tailings, waters are relatively more absorbent with b:a ratio of 0.8 at the blue range, and increasing towards the red range (b:a of 20.0 at 660nm).The relative low scattering properties of these waters result in lower spectral K d and, consequently, a more homogeneous spectral light field with a significant amount of light at the green spectra (Figure 7a).On the other hand, with an increase in TSS loadings from mining operations, the scattering coefficient increases.In this case, the scattering process prevails over absorption coefficient, specifically at the green and red wavelengths, thus explaining the observed higher underwater light availability at the green and, mainly, at the red wavelengths (Figure 7).This is particularly the case when analyzing E o since both the upwelling and downwelling light fields are considered in the calculations.In turbid waters where the b:a in the blue spectra is approximately 30.0, such as the case for Class 4, E o (0.3m) values became 50% higher than E d (0 + ) (Figure 7).Overall, considering the scalar irradiance (E o ) for assessment of light availability is particularly important for highly turbid waters where increased E u contributes significantly to scalar irradiance (E o ) available for suspended phytoplankton.

Underwater Light Field and Phytoplankton
As an initial attempt to evaluate the consequences of reduced light to the phytoplankton biovolume, Z c was calculated for the studied waters considering the light requirements, E c (PAR) , for Chlamydomonas sp.(Chlorophyta) [47].However, Z c alone does not explain PAR availability for phytoplankton growth.A second parameter is needed for this analysis, which is the depth of the mixing layer (Z m ).In a river system, Z m depends on water mixing originated from bottom interaction (friction) and water flow [70,71].For the studied rivers, a deep mixing layer was considered due to the reported high turbulence and water flow.On average, water velocity varies from 1.09 ms −1 during the high water level period to 0.55 ms −1 during the low water level period for several studied tributaries, such as the Crepori, Jamanxim, and Tocantins rivers [72].For reference, water velocities over 0.5 ms −1 are considered the threshold velocity for phytoplankton growth as indicated by Salmaso et al. [73].Therefore, in both seasons, those rivers are likely to have their phytoplankton production limited by turbulent flux.
Given the high turbulence, this study assumes that the mixing layer, Z m , is equivalent to the total water depth (Z) to calculate Z m /Z c .For non-impacted tributaries, such as the upstream Tapajós and the Jamanxim rivers (Class 1), the Z m /Z c is lower than 1.3, indicating that there is enough light available in approximately 75% of the water column.On the other hand, impacted rivers, such as the Crepori (Class 4) and the Tocantins (Class 3), showed Z m /Z c of ~2.0 during the high water level period, which suggests a reduction of light availability to half of the water column.Analogous restriction of light availability has been reported in the Amazon River and it is associated with the low chl-a concentration given that nutrients are not limited in these waters [24,32,50,56].
Although the phytoplankton data available for this study is limited, the results show low values of phytoplankton biovolume and chl-a (Figure 4) for both impacted (Crepori, Tocantinzinho, and Novo rivers) and non-impacted rivers (Tapajos upstream, Jamanxin upstream rivers), except the Tapajós lake region.For the impacted and non-impacted rivers, our findings suggest that the increase in PAR attenuation due to increase in water siltation does not necessarily cause a direct reduction in light availability for phytoplankton growth.Here, we are suggesting that a combination of light availability and river turbulence, as described by Z m /Z c , control phytoplankton growth, assuming these rivers are not nutrient limited.Interestingly, seasonal factors may play a more important role in biovolume and pigment concentration than the degree of mining.For instance, phytoplankton biovolume data and chl-a for both groups of rivers (impacted and non-impacted) are higher during the low water level period, which is possibly related to other factors such as seasonal irradiance levels, nutrient availability, and water flow velocity [34].In the Amazon region, the seasonal irradiance regime shows relatively lower values during the high water period (rainy season) (~950 mol•m −2 •mo −1 ) then during the low water period (~1300 mol•m −2 •mo −1 ) [74].The higher irradiance during the low water period further explains the observed higher phytoplankton production in some of the Amazonian rivers [75].On the other hand, our data showed that during the rainy season, phytoplankton biovolume and chl-a decreased, most likely because of monthly PAR irradiance reduction [74,75].In fact, according to INMET (inmet.gov.br)data, the number of sunny hours for the period of our data collection varied from 231 during the low water level period to 126 sunny days during high water level; that is, 1.83 times higher during the low water period.

Conclusions
The ASGM expansion over the past 40 years in the study area [51,55] has directly impacted the investigated rivers with high sediment content and mercury contamination [21].To complement these previous studies, the present article investigated the effects of water siltation on the underwater light field of various tributaries within the Tapajós River Basin.The effects of TSS derived from mining activities on both the inherent and apparent optical properties were quantified.The main conclusions are as follows.
The inorganic nature of mine tailings is the main factor affecting the underwater scalar irradiance in the Tapajós River Basin.The TSS concentration varies seasonally during the year in a synergism between water level and mining activities: during low water level periods, mining activities intensify and, associated with low water volume, TSS rapidly increases, which in turn changes the optical characteristics of the water.For tributaries with low or no influence of mine tailings, waters are relatively more absorbent.On the other hand, with TSS loadings from mining operations, the scattering process prevails over the absorption coefficient at the green and red wavelengths, thus explaining why we observed higher underwater light availability at the red wavelengths in highly turbid waters.Moreover, ASGM not only affects the sediment load but also its natural seasonality since in non-impacted basins the maximum TSS is caused by diffuse erosion during the rainy season.This change in load and seasonality might affect, in the long run, biota composition of a previous clear water environment to a distinct light availability regime as the river becomes subjected to mining operations.
Although a strong light reduction due to water siltation was observed in impacted tributaries, the effects on the phytoplankton community was not clearly observed, which can be attributed to a combination of (i) minimal number of samples for proper comparison between impacted and non-impacted tributaries, and (ii) general low phytoplankton productivity in all upstream tributaries due to a combination of high current velocity (flushing) in rivers and light conditions.Therefore, due to the lack of comparable (impacted vs. non-impacted) data on phytoplankton species and concentration, the question of whether the reduction of light negatively affects phytoplankton growth requires further investigation.
In previous publications [51,55], we provided historical maps with TSS retrieved from satellite images indicating that the optical description provided in this article is directly related to remote sensing signal, and can be extended to several rivers that are subject to similar mining impacts.In fact, using remote sensing techniques, ongoing research indicates that other clear and black water rivers such as Xingu River and Amanã River, respectively, are also presenting effects of mine tailings on the aquatic environment.As such, information derived from this research and other ongoing initiatives will support control actions by IBAMA (Brazilian Environmental Protection Agency) and ICMBio (Brazilian Agency of Protected Areas) towards illegal mining activities along streams and within protected areas in the Amazon.Further, the hydrological approach by Abe et al. [58] will aid in understanding the sources of suspended sediment by quantifying the natural runoff sediment production and comparing it to the sediment concentration observed in the present article.

Figure 1 .
Figure 1.Study area.(a) Overview of the Brazilian Amazon.(b) Tapajós River Basin in the Brazilian Amazon showing the main tributaries, sample sites (see Section 4.1), deforestation [22] (light grey), and mines [23].(c) Water flow (Q) and water speed (v) are also shown for the Tapajós River at the Itaituba City region during high and low water level.

Figure 1 .
Figure 1.Study area.(a) Overview of the Brazilian Amazon.(b) Tapajós River Basin in the Brazilian Amazon showing the main tributaries, sample sites (see Section 4.1), deforestation [22] (light grey), and mines [23].(c) Water flow (Q) and water speed (v) are also shown for the Tapajós River at the Itaituba City region during high and low water level.

Figure 2 .
Figure 2. The methodology comprises (a) field campaigns for acquisition of optical properties and biogeochemical measurements to quantify underwater light field changes from non-impacted to impacted rivers; (b) assessment of light availability for phytoplankton, including critical depth analyses.

Figure 2 .
Figure 2. The methodology comprises (a) field campaigns for acquisition of optical properties and biogeochemical measurements to quantify underwater light field changes from non-impacted to impacted rivers; (b) assessment of light availability for phytoplankton, including critical depth analyses.

Figure 3 .
Figure 3. Schematic representation of TSS concentrations along the Tapajós River and its tributaries for high (April 2011) and low (September 2012) water levels.The river and its tributaries were classified according to TSS (mg•L −1 ) concentrations.TSS for Tocantins, Novo, and Jamanxim rivers during the low water season were retrieved from Landsat surface reflectance (red band) [51].The level of mining impact is an arbitrary classification considering the intensity level of mining and mining area distribution.Not to scale.

Figure 3 .
Figure 3. Schematic representation of TSS concentrations along the Tapajós River and its tributaries for high (April 2011) and low (September 2012) water levels.The river and its tributaries were classified according to TSS (mg•L −1 ) concentrations.TSS for Tocantins, Novo, and Jamanxim rivers during the low water season were retrieved from Landsat surface reflectance (red band) [51].The level of mining impact is an arbitrary classification considering the intensity level of mining and mining area distribution.Not to scale.

Figure 4 .
Figure 4. (a) Spatial distribution of phytoplankton groups (mm 3 •L −1 ) along the Tapajós River for low and high water level periods.(b) Spatial distribution of pigments concentration (μg•L −1 ) for the same seasons.The correspondent TSS class is indicated in parentheses for each sample point, except for point stations with no data available, which is indicated by asterisks.

Figure 4 .
Figure 4. (a) Spatial distribution of phytoplankton groups (mm 3 •L −1 ) along the Tapajós River for low and high water level periods.(b) Spatial distribution of pigments concentration (µg•L −1 ) for the same seasons.The correspondent TSS class is indicated in parentheses for each sample point, except for point stations with no data available, which is indicated by asterisks.

Figure 5 .
Figure 5. Spectral distribution of in situ IOPs for different water classes: absorption by particles (a) and CDOM (b); and particulate scattering (c) and backscattering (d).Gray scale curves represent the five classes of water as explained in the section above.Note that Class 5 is not shown in (b-d) due to lack of data.

Figure 5 .
Figure 5. Spectral distribution of in situ IOPs for different water classes: absorption by particles (a) and CDOM (b); and particulate scattering (c) and backscattering (d).Gray scale curves represent the five classes of water as explained in the section above.Note that Class 5 is not shown in (b-d) due to lack of data.

Figure 6 .
Figure 6.Diffuse attenuation coefficient Kd (λ), for both field campaigns for classes under different mining impacts from low impact (Class 1) to very high (Class 5).

Figure 6 .
Figure 6.Diffuse attenuation coefficient K d (λ), for both field campaigns for classes under different mining impacts from low impact (Class 1) to very high (Class 5).

Figure 7 .
Figure 7. Normalized scalar irradiance at (a) 0.3 m and (b) 2.0 m for all samples grouped by TSS.(c) Spectral profile of Z1% averaged for each class, and (d) Eo(PAR) availability from surface to bottom with depth.Compensation irradiance Ec(PAR) for Chlamydomonas sp. is indicated (thick vertical black line).The correspondent critical depth, Zc(PAR), for each class can be drawn from the intersection of Eo(PAR) with the Ec(PAR) line.Emax, is the maximum level of irradiance that yields photosynthesis.Above this point, the incident light combined with upwelling light can be harmful for phytoplankton cells.The optimum irradiance level, Ek, is between Ec (critical) and Emax.

7 .
Normalized scalar irradiance at (a) 0.3 m and (b) 2.0 m for all samples grouped by TSS.(c) Spectral profile of Z 1% averaged for each class, and (d) E o (PAR) availability from surface to bottom with depth.Compensation irradiance E c (PAR) for Chlamydomonas sp. is indicated (thick vertical black line).The correspondent critical depth, Z c (PAR), for each class can be drawn from the intersection of E o (PAR) with the E c (PAR) line.E max , is the maximum level of irradiance that yields photosynthesis.Above this point, the incident light combined with upwelling light can be harmful for phytoplankton cells.The optimum irradiance level, E k , is between E c (critical) and E max .

Table 1 .
Compilation of all optical parameters used in this paper, including respective symbol, unit, and calculation formula.
April 2011 during high water level (27 sample points) and September 2012 during low water level periods (13 sample points, no IOPs data due to loss of ac-S optical equipment in the Tapajós River at the beginning of the field work) (see Figure

Table 2 .
Values of suspended solids, chl-a, and optical properties of the water classes represented in Figure3.* IOPs for low water samples were estimated based on linear regressions established for the measurements taken during the high water period.& Parameters related to photic zone for phytoplankton productivity.Critical depth (Z c ) and depth where scalar irradiance reaches 1% of the incoming light (Z 1% ) are approximately the same values.

Table 3 .
Correlation coefficients between IOPs and suspended solids, chl-a and AOPs.Data are from high water only (n = 27).* indicates significant correlation at p < 0.05.

Table 3 .
Correlation coefficients between IOPs and suspended solids, chl-a and AOPs.Data are from high water only (n = 27).* indicates significant correlation at p < 0.05.

Table 4 .
Correlation coefficients between suspended sediment, chl-a, and acdom, with measured and modeled AOPs.Data are from high and low water (n = 40).* indicates significant correlation at p < 0.05.

Table 4 .
Correlation coefficients between suspended sediment, chl-a, and a cdom , with measured and modeled AOPs.Data are from high and low water (n = 40).* indicates significant correlation at p < 0.05.

Table A1 .
Mann-Whitney U Test for TSS concentration Marked tests are significant at p < 0.10.