Next Article in Journal
Nitrogen Removal Characteristics and Comparison of the Microbial Community Structure in Different Anaerobic Ammonia Oxidation Reactors
Previous Article in Journal
Water Price Prediction for Increasing Market Efficiency Using Random Forest Regression: A Case Study in the Western United States
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Light Absorption Budget in a Reservoir Cascade System with Widely Differing Optical Properties

by
Nariane Bernardo
1,
Enner Alcântara
2,*,
Fernanda Watanabe
1,
Thanan Rodrigues
3,
Alisson do Carmo
1,
Ana Carolina Campos Gomes
1 and
Caroline Andrade
1
1
Department of Cartography, São Paulo State University—UNESP, Presidente Prudente, SP 19060-100, Brazil
2
Department of Environmental Engineering, São Paulo State University—UNESP, São José Dos Campos, SP 12245-000, Brazil
3
Federal Institute of Education, Science and Technology of Pará State—IFPA, Castanhal, PA 68740-970, Brazil
*
Author to whom correspondence should be addressed.
Water 2019, 11(2), 229; https://doi.org/10.3390/w11020229
Submission received: 9 October 2018 / Revised: 15 January 2019 / Accepted: 23 January 2019 / Published: 29 January 2019
(This article belongs to the Section Water Quality and Contamination)

Abstract

:
Aquatic systems are complex systems due to the environmental pressures that lead to water quality parameter changes, and consequently, variations in optically active compounds (OAC). In cascading reservoir systems, such as the Tietê Cascade Reservoir System (TCSR), which has a length of 1100 km, the horizontal gradients are expressive due to the filtration process that is caused by the sequence of dams affecting the light absorption throughout the cascade. Our new observations showed that colored dissolved organic matter (CDOM) dominate two reservoirs; non-algae particles (NAP) dominate one, and phytoplankton dominates the other. The variability of light absorption along the cascade indicates the influence of watershed dynamics in the reservoirs as much as the flow driven by previous reservoirs. Despite the effect of the variability of light absorption, light absorption by phytoplankton strongly affects the total absorption in the four reservoirs in TCSR. The results obtained in this work may enable a better understanding of how the gradient pattern changes primary production and indicates a challenge in retrieving OAC concentrations using a bio-optical model for an entire cascade composed of different optical environments.

1. Introduction

Water quality monitoring has been efficiently performed via remotely sensed data, using empirical and analytical models, and when applied to time-series images, it provides synoptic scale observations [1,2,3,4]. The success of remote sensing applications to aquatic systems depends on the understanding of underwater light interactions with the optically active compounds (OACs) [5]. These interactions are expressed by inherent optical properties (IOPs), the absorption and scattering of light, and the OACs are represented by chlorophyll-a (Chl-a), suspended particulate matter (SPM), and colored dissolved organic matter (CDOM) [5].
Much attention has been devoted to characterizing aquatic systems and determining the spatial variability of IOPs in the ocean [5,6,7,8,9], coastal environments [10,11,12], and bays [13,14]. For inland waters, efforts have been made to characterize Australian [15], African [16], and Brazilian reservoirs [3,17] to identify the optical differences between open oceans and freshwater systems. The newest studies have also been developed to classify the inland water systems [18,19,20] and overcome the challenge of bio-optical modeling in such environments [21].
Bio-optical modeling is used to create a remotely sensed-based algorithm, which can be applied to satellite images and used to monitor water quality variability over time. In inland waters with specially linked water systems, such as cascading reservoirs, a high variability in optical properties prevents the use of only one model to estimate the water quality parameters. For example, Rodrigues et al. (2017a) studied two reservoirs with widely differing optical properties and tried to estimate the SPM concentration using only one adjusted model. However, this was not possible; each reservoir needed a local model [22]. An understanding of how light is absorbed and how a connected system influences not only the absorption but also the water quality parameters is mandatory for operational monitoring from space. Until now, most studies have investigated and proposed the development of new bio-optical models; however, few studies have aimed to understand light absorption variability and how it has been conducted on cascading systems.
In contrast to open oceans, inland waters are optically complex environments, due to the different influences of OACs [18]. Moreover, inland waters are influenced by internal processes and environmental pressures, such as nutrient enrichment, wastewater discharges, and other non-natural disturbances [23,24,25]. In cascade reservoir systems, the existence of sequential dams modifies the water retention time of continuous rivers and impacts the filtration process, which changes the OACs concentrations and compositions from upstream to downstream, as predicted by the Cascade Reservoir Continuum Concept (CRCC) [26]. Relying on the OAC’s variability in TCSR, a relevant scientific question is raised: what is the behavior of the IOPS throughout cascading reservoirs, since the IOPs are related to the composition and morphology of OACs? We hypothesize that a strong relationship would exist between the concentration of OACs and the absorption properties in aquatic continuum systems, and both would be affected similarly by autochthonous and allochthonous strengths throughout the course of the Tietê River.
Considering its 1100 km of extension, the Tietê Cascade System Reservoir (TCSR) environment is a challenging and relevant region to investigate, since its water attends to fishery activities, navigation, hydroelectric productivity, and supplies industrial, agriculture, and domestic demands [26,27]. The TCSR’s waters are deeply affected by different non-natural contamination sources, such as discharges from pasture waste, wastewater from urban centers, and other agricultural activities developed in nearby areas, including sugarcane and citrus crops [4,28]. As consequence, Barra Bonita (BB) and Bariri (BAR) reservoirs are considered eutrophic environments [29], Ibitinga (IBI) is considered a meso-to-eutrophic [30,31], and Nova-Avanhandava (NAV) is considered an oligo-to-mesotrophic environment [3].
Understanding the optical features throughout long rivers is crucial for a better understanding of the riverine transfers and their potential role in global biogeochemical cycles. Therefore, the present data are suitable to contribute to the understanding of the dynamics of this type of aquatic system, aiming to gain insights into the optical information of South American aquatic continuum systems by sampling the absorption coefficients from non-algal particles (NAP), chlorophyll-a (Chl-a), and colored dissolved organic matter (CDOM) in TCSR. The specific goals of this work were to analyze the optical properties in an aquatic continuum river, identify whether there was a downstream gradient in spectral absorption, and verify the controls on the light absorption budget. Considering the optical characteristics as well known, remote sensing products can be better used for continuous monitoring and establishing effective control strategies to improve water quality, mainly for larger rivers such as the Tietê River.

2. Materials and Methods

2.1. Study Area

The Tietê River runs from west to east for 1100 km in São Paulo State, passing through the capital, and flows into the Paraná River, which is a tributary of the Plata River. There are six hydroelectric reservoirs placed along the Tietê Cascade System of Reservoirs (TCSR) from upstream to downstream. They are: BB (22.50° S; 48.51° W), BAR (22.15° S; 48.73° W), IBI (21.65° S; 48.98° W), Promissão (21.48° S; 49.83° W), NAV (22.11°S; 50.17° W), and Três Irmãos (20.68° S; 51.32° W) (Figure 1—Datum WGS 84, geographic coordinate system). BB, BAR, IBI, and NAV reservoirs are responsible for more than 90% of the hydroelectric productivity from the entire cascade, and present flooded areas of 310 km2, 63 km2, 114 km2, and 210 km2, respectively. NAV, IBI, and BAR are run-of-river reservoirs, whilst BB is a storage system.

2.2. Field Survey

Eight field campaigns were made in TCSR, with two field studies per reservoir (see details of the sampled locations in Figure 1d–g); therefore, each dataset was labeled with regard to the respective locations and time order: BB1, BB2, NAV1, NAV2, BAR1, BAR2, IBI1, and IBI2. The field studies were carried out between 5–9 May 2014 in BB1; 13–16 October 2014 in BB2; 15–18 August 2016 in BAR1; 23–24 June 2017 in BAR2; 19–23 July 2016 in IBI1; 21–22 June 2017 in IBI2; 28 April to 2 May 2014 in NAV1; and 24–26 September 2014 in NAV2. The date for each fieldwork was established with regard to the Landsat-8 overpass, to have in situ measurements that were coincident with the satellite data for future applications. Furthermore, all of the field measurements were carried out during the dried periods (Figure 1c), which did not allow us to make seasonal interpretations. Principal component analysis was performed using the standard deviation image, which was retrieved from the temporal series of the Operational Land Imager (OLI) (onboard Landsat-8). The first component of the standard deviation (SD) image was then used to determine the sampling locations within each reservoir (more details are given in [32]).
Physical parameters, such as temperature (°C), turbidity (NTU), and depth (m) were measured in triplicate in all of the sampling locations. The Secchi disk depth (ZSD, m−1) was achieved with a black and white disk with a 30-cm diameter. Furthermore, water samples were collected near the surface (depth = 0.3 m) in five-liter black plastic bottles that were rinsed three times with the water from each sampling spot. To preserve the original features of the samples, the bottles were stored in the dark at low temperatures until laboratory analysis, where absorption spectra measurements were made. The SPM and Chl-a concentrations were established following the American Public Health Association (APHA) (1998) and Golterman et al. (1978) protocols [33,34], respectively.

2.3. Absorption Coefficients

The IOPs are represented by a(λ), scattering coefficients (b(λ)), and the volume scattering function (VSF). To analyze the optical composition of the cascade system in terms of IOPs, only the absorption spectra measured in situ were assessed in this work. The total absorption spectra at(λ) that resulted from the contribution of each OAC spectrum (Equation (1)) can be mathematically described as follows:
a t ( λ ) = a c d o m ( λ ) +   a ϕ ( λ ) + a n a p ( λ ) + a w ( λ ) ,
where at(λ) is the total absorption spectrum, acdom(λ) is the absorption spectrum of CDOM (Equation (2)), aϕ(λ) is the absorption spectrum of phytoplankton (Equation (3)), anap(λ) is the absorption spectrum of non-algal particles (NAP, Equation (4)), and aw(λ) is the absorption spectrum of pure water from Smith and Baker et al. (1981) [35].
All of the absorption coefficients were derived using absorbance readings (ABS) on a SHIMADZU UV-2600 UV–VIS spectrophotometer (SHIMADZU, Japan) with a one-nm spectral resolution, ranging from 280 to 800 nm. The acdom(λ) was established after water samples were filtered through a nylon membrane with a 47-mm diameter and 0.22-μm pore size. The filtrate was placed in quartz cuvettes, with an optical path of 10 cm (L = 0.1 m). Then, the spectral ABS readings were taken to derive acdom(λ) after filtration and acclimation to room temperature (Equation (2)):
a c d o m ( λ ) = 2.303 × A B S ( λ ) L ,
where L is the optical path length of the cuvette, and ABS is the spectral absorbance reading within 280 to 800 nm of the interval. Baseline data were obtained by using Milli-Q water placed in a 10-cm quartz cuvette. The conversion factor (2.303) between the base 10 logarithm and natural logarithm was used, and was applied to the baseline correction by subtracting the residual absorbance, i.e., the average value between 700–750 nm. The acdom spectra were established by the methodology proposed by Bricaud et al. (1981) [6], and the fit of the exponential function was also computed to achieve the Scdom results [9].
To compute anap and aϕ, water samples were filtered through fiberglass filters (Whatman GF/F with a 0.7 pore size and 47-mm diameter), placed inside the integrated sphere module using a blank filter as reference, and then wetted with drops of Milli-Q water during the laboratory analysis [16]. The transmittance and reflectances were read to compute the particulate absorption (ap), following the transmittance-reflectance method [36,37].
The bleaching of the organic fraction of SPM was performed using hypochlorite solution (NaClO) at 10%. The bleached filters were washed with Milli-Q water to prevent the eventual contamination of the remaining NaClO, and new transmittances and reflectances readings were made to compute anap. The aϕ was achieved by subtracting anap from ap (Equations (3) and (4)):
a ϕ ( λ ) = a p ( λ )   a n a p ( λ )
a n a p ( λ ) = a n a p ( λ 0 ) × e S n a p ( λ λ 0 )
where Snap is the slope of the exponential decreasing function considering a spectral range (the present study uses 400 to 700 nm), which mathematically describes the relative steepness of the curve [38], and λ 0 is the reference wavelength, which is usually adopted near 443 nm for acdom and anap. To avoid systematic errors due to different methodologies, all of the absorption values were established following specific protocols for each OAC.

3. Results

3.1. Water Quality Scenery

The descriptive statistics of water quality parameters measured in situ are summarized in Table 1. In BB (both campaigns), particulate inorganic matter (PIM) presented the highest variation, with CV = 78.8% and 37.3%, respectively. BAR1 showed the highest variation in Chl-a (CV = 80.5%), while BAR2 presented the highest variation for PIM (CV = 42.9%). In IBI (both campaigns), the larger variability resulted from Chl-a (CV = 86.0% and 67.8%), and for NAV (both campaigns), the larger variability was obtained in PIM (CV = 76.8% and 65.8%, respectively). The lowest variability resulted from temperature for all of the field studies, with CV comprising between 1.90–15.4%.
The TCSR showed pH values ranging between 5.50 (IBI2) and 10.1 (BB2), and the temperature ranged from 21.1 °C to 39.4 °C. The highest SPM and Chl-a concentrations varied between three orders of magnitude in TCSR, whereas PIM varied between two orders of magnitude. BB2 reached the maximum concentrations (SPM = 44.0 mg·L−1 and Chl-a = 797.8 mg·m−3), whilst NAV1 reached the lowest SPM concentration (SPM = 0.1 mg·L−1). The lowest Chl-a concentration was found in IBI1 (Chl-a = 1.4 mg·m−3). The full TCSR dataset demonstrated high correlations between Chl-a and SPM concentrations, with r = 0.95 (p-value < 0.05). It is important to note that the Chl-a and particulate organic matter (POM) fractions showed r = 0.94 (p-value < 0.05, plot is not depicted here), whilst the SPM and POM fractions showed r = 0.99 (p-value < 0.05). The SPM and PIM concentrations also showed significant statistical relations (r = 0.67, p-value < 0.05).
The PIM and POM concentrations showed an organic-to-inorganic standard from upstream to downstream reservoirs. Although the average proportions are similar between all of the reservoirs (Figure S1), we observed expressive POM concentrations in BB and BAR compared to PIM concentrations, as well as the expressive PIM concentrations in NAV (Table 2) compared to POM concentrations. Details about the inorganic-to-organic relation along the cascade reservoirs are presented in the Supplementary Materials (Figure S1).
Inverse relationships were observed between water transparency (represented by ZSD) against Chl-a, as well as the ZSD against PIM, ZSD against SPM, and ZSD against POM. Through the cascade, the ZSD was slightly more influenced by SPM, with r = 0.70 (p-value < 0.05) than Chl-a, PIM, and POM (r = 0.66, 0.66, and 0.69, respectively). The minimum value of ZSD was obtained in BB2 (= 0.37 m), and the maximum was found in NAV1 (= 4.80 m). Overall, the downstream reservoirs (NAV and IBI) presented higher ZSD, and the upstream reservoirs (BB and BAR) showed lower ZSD values.
The boxplots (Figure 2) illustrated the water quality parameter gradients throughout the cascade system. High SPM concentrations were accomplished by the high Chl-a concentrations in upstream reservoirs (BB and BAR reservoirs in Figure 2a,b). The turbidity (Figure 2c) decays from upstream to downstream (from BB to NAV), following the SPM and Chl-a concentration behaviors. The ZSD is lower in reservoirs with high turbidity, and it is higher in reservoirs with low turbidity, as expected (Figure 2d). Although the reservoirs demonstrated the interconnection process, the temperature and pH did not present any detectable pattern (boxplots are not shown).

3.2. Absorption Spectra

3.2.1. Absorption by CDOM

The CDOM absorption spectral curves for the TCSR dataset (acdom(λ)), with 175 sampling spots, are depicted in Figure 3a, which followed the exponential fit, as expected.
Regarding the TCSR values, the acdom(443) ranged from 0.21 m−1 (NAV1) to 3.17 m−1 (BAR2), as shown in Table 2. The assessment of the coefficient of variation (CV) depicted BAR2 as having the highest acdom(443) variability (CV = 33.2%), and NAV1 as having the lowest variability (CV = 8.5%). To quantify the distinct decay between the fieldwork, the slopes of CDOM (Scdom) were calculated by fitting an exponential function and non-linear least squares method. Scdom comprised between 0.004 (in BB2) and 0.020 nm−1 (in NAV2) (Table 3). Although the Scdom range for BB1, NAV1, and NAV2 are similar (0.016–0.018; 0.014–0.017; and 0.016–0.020 nm−1), a paired t-test showed that the Scdom values are different at 5% of significance. The same conclusion was obtained from comparisons between BB2 versus BB1, NAV1, NAV2, and between BAR1 versus IBI1.
Regarding the frequency for all of the Scdom values, the highest frequency was found in the 0.010–0.014 nm−1 range (Figure S2a). Moreover, investigations about relationships between Scdom and acdom(443) demonstrated statistically significant relationships in BB (both field studies) and IBI2 (Figure S2b–d), with r values of −0.64, −0.87, and −0.48, respectively. Other fieldwork showed weaker relationships, such as IBI1 (r = 0.45) and NAV1 (r = −0.37). BAR (both field studies) and NAV2 did not indicate a statistical relationship between acdom(443) and Scdom.
The relations between acdom(443) versus Chl-a concentrations, and acdom(443) versus SPM concentrations, were also investigated. Considering the entire TCSR dataset, Chl-a and acdom(443) presented coefficients of correlation where −r = 0.15 (p-value < 0.05). However, taking into account each dataset separately, the relationships changes, with BB1, BAR2, IBI2, and NAV2 (Figure S3) showing significant relationships at 5% significance (r = 0.77, 0.66, 0.51, 0.61, respectively). Other field studies did not show relationships with statistical significance (BB2, BAR1, IBI1, and NAV1).
Considering SPM and acdom(443) relationships, the entire TCSR dataset showed r = 0.20 (p-value < 0.05). This scenery changes when field campaigns were individually considered. BB1, IBI1, and NAV1 showed r = 0.54, 0.49, and −0.06 (p-value < 0.05), respectively (graphs are not shown here). The other fields (BB2, BAR1, BAR2, IBI2, and NAV2) did not show significant relationships between SPM and acdom(443).

3.2.2. Absorption by Particulate Matter

All of the sampled particulate absorption spectra, ap(λ) are plotted in Figure 3b, and demonstrated time-fold variation, with the highest magnitudes in BB2, BAR1, and BB1. The lowest magnitudes were found in BAR2 and IBI2. Considering inorganic and organic particle composition, it is worth highlighting that the phytoplankton absorption contributions to the ap(λ) curves were mainly in BB (both field studies), due to the absorption features at 443 nm and 675 nm. In contrast, the inorganic absorption contribution can be highlighted in NAV (both field studies), due to the exponential shape of the absorption curves.
The values of ap(443) for the TCSR dataset ranged from 0.10 m−1 to 11.2 m−1 (see Table 2 for details). The highest variability was found in BAR1 (CV = 65%), and the lowest variability was encountered in NAV2 (CV = 26%). In addition to the ap curves, inorganic and organic particulate contributions were also evaluated by plotting the NAP (anap(λ)) and phytoplankton (aϕ(λ)) absorption spectra (Figure 3c,d, respectively).
The absorption of blue light (aϕ(443)) decreased markedly from the upstream site BB2 (10.9 m−1) to the downstream site NAV1 (0.02 m−1) (Table 2). The anap(λ) curves did not present the same standard from upstream to downstream, with anap(443) ranging between 0.03 m−1 (NAV1) and 1.7 m−1 (BB2) (Table 2).
Two highlighted absorption features are observed in Figure 3d: one at 443 nm and another at 675 nm, which are due to Chl-a absorption. The highest magnitude of phytoplankton absorption was found in BB2, with aϕ(443) = 10.9 m−1, and aϕ(675) = 5.44 m−1 (Table 2). The lowest magnitudes were found in NAV1 (aϕ(443) = 0.02 m−1) and BAR2 (aϕ(675) = 0.004 m−1). IBI1 depicted the highest variability among all of the field sites, with CV443nm = 90% and CV675nm = 95%.
Some fieldwork (IBI2, NAV1, and NAV2) did not show statistically significant relationships between Chl-a concentration and aϕ(443) or aϕ(675). BB1, BAR1, and IBI1 presented high relationships, and did not vary the relation with wavelength (r = 0.90, r = 0.99, and r = 0.89, respectively at 443 nm and 675 nm). In BB2 and BAR2 fieldwork, relationships showed slight differences at each wavelength: at 443 nm, BB2 presented r = 0.64 and BAR2 presented r = 0.43; whereas at 675 nm, BB2 presented r = 0.54 and BAR2 presented r = 0.65 (BB2 and BAR2 plots can be found in Supplementary Information, Figure S4).
Regarding anap(λ) curves (Figure 3c), the exponential shape is a result of the isolation of residual pigment absorptions [9]. This means that there were neither phytoplankton nor accessory pigments features in the anap spectra. Among all of the field campaigns, only BAR1 and BAR2 depicted a remarkable change in magnitude for the entire spectrum dataset, demonstrating a considerable division between the two datasets (two to three times the magnitude of difference). IBI1 and IBI2 presented a slight overlap of some curves (the lowest curves for IBI1 overlap the highest curves for IBI2). Overall, BB1 and BB2, and NAV1 and NAV2 presented high similarity among the curves, except for some samples that presented the highest magnitudes.
In terms of magnitude, BB2 showed the highest values for the entire spectra, and NAV1 showed the lowest magnitude (Figure 3c). Among all of the field sites, the largest variability of anap(443) was encountered in BB2, with CV = 77%, and the smallest variability was in BAR1, with CV = 23% (Table 2).
The entire cascade presented Snap values ranging between 0.003–0.023 m−1. Statistical analysis for Snap (Table 2) demonstrated that IBI2 presented the highest variability (0.005 nm−1 < Snap < 0.022 nm−1, and CV = 39%), and the lowest variability was found in BAR (both fieldworks with CV = 4%). The frequency for the entire dataset of Snap values was analyzed and showed a normal distribution, with the highest frequency in the 0.008–0.010 m−1 range. The relationships between anap(443) and PIM concentrations were evaluated, and a statistical relationship was found only in BAR1, with r = 0.50 (p-value < 0.05).

3.3. Relative Contribution of OAC’s Absorptions in TCSR

The at curves represent the sum of all of the individual OAC’s absorption spectra with the pure water absorption curve (aw (λ), Smith and Baker, 1981), as described Equation (1). To identify each OAC absorption contribution to the total absorption without water (at-w), the average curves are plotted in Figure 4.
Regarding at-w, the BB (both field studies) dataset is markedly influenced by aϕ. BAR1 and BAR2 are distinctly different: BAR1 showed absorption features near 443 nm, 620 nm, and 675 nm, resulting from aϕ contributions, while BAR2 presented a continuous exponential decrease along the spectra (Figure 4c,d). IBI1 presented smoothed features near 443 nm and 620 nm, which were derived from aϕ contributions (Figure 4e), whilst IBI2 presented an exponential decrease (Figure 4f). NAV1 and NAV2 (Figure 4g,h) depicted an exponential decrease with anap and acdom contributions up to 670 nm, where aϕ became expressive. NAV1 showed the anap predominance, and NAV2 presented a main contribution from acdom, mainly up to 450 nm.
To compare the spectral behavior throughout the cascade, Figure 4 was reorganized in comparative graphs, as shown in Figure 5. The magnitude of aϕ(λ) (Figure 5a) coincides with the order of reservoirs, since the magnitudes followed in the sequence of BB2, BAR1, BB1, IBI1, NAV2, and the remainder fields. The anap(λ) presented varied magnitudes (Figure 5b), with BAR1 higher than BB2, BB1, IBI1, NAV1, NAV2, BAR2, and IBI2. The acdom(λ) also did not show the upstream to downstream decay, with magnitudes varying in descending order from BAR1, BAR2, IBI1, IBI2, BB2, BB1, NAV2, and NAV1 (Figure 5c).

3.4. Absorption Budget

The ternary diagram represented the OAC’s absorption contributions to the absorption itself (at-w(λ)) without water. Then, the diagrams were plotted in Figure 6 for the specific wavelengths of the Operational Land Imager sensor onboard ongoing Landsat 8 (OLI/Landsat-8): 443 nm, 482 nm, 560 nm, and 655 nm. The OLI sensor was used because its images are suitable for inland water applications due to the spatial resolution and revisit time [21].
The at-w(443) (Figure 6a) showed different OAC contributions in each reservoir. Overall, in BAR (both field studies) and IBI (both field studies), the major contributor was CDOM (averages of 53.3%, 82.74%, 64.9%, and 87.30% for BAR1, BAR2, IBI1, and IBI2). In BB (both field studies), we identified a marked contribution from phytoplankton absorption (BB1: 20.7–67.8%; BB2: 48.1–87.9%). NAV1 presented a major influence of NAP (43.2%), whereas NAV2 presented high contributions from CDOM and NAP (39.4% and 31.7%, respectively). Considering the region of Photosynthetically Active Radiation (PAR) from 400-nm to 700-nm range and including the contribution of water to the total absorption, we maintain our findings for almost all of the field works, excepting NAV2, where water presented a dominance of 41% against 30% from NAP (Table 3).
The at-w(482) (Figure 6b), showed a high contribution of aϕ to BB (BB1: 40.07% and BB2: 73.06% of averages). The acdom contribution was predominant in BAR1, BAR2, IBI1, and IBI2 (averages of 54.6%, 82.5%, 63.5%, and 86.2%, respectively). NAP was dominant in NAV1 (30.4–58.0%), whereas NAV2 presented a slightly higher aϕ contribution (37.8%) compared to anap (31.7%).
For the at-w(560) (Figure 6c), BB1 increased the dominance of the aϕ (42% on average) and anap (32% on average) contributions. BB2 presented higher aϕ contributions (31.5–95.7%), and some samples showed a slight increase in NAP contributions. BAR1 presented a high CDOM contribution (30.1–77.6%) with two samples that showed remarkable aϕ(λ) contributions (they also presented the highest Chl-a concentrations in the BAR dataset). BAR2 maintained a high contribution from CDOM (47.0–94.1%). IBI1 showed a strong influence by CDOM (41.4–90.5%), as well as IBI2 (74.3–91.8%); BAR and IBI datasets were clustered toward the CDOM apex (Figure 6c). NAV1 has anap-dominant (contributions ranging between 28–70%) with one CDOM-dominated sample (acdom = 51.6%). Overall, NAV2 conserved the high influence of anap (28.6–70.2%), with one CDOM-dominated sample (acdom = 58%).
The at-w(655) (Figure 6d) depicted a high influence by aϕ in BB1 (52.8–90.1%) and BB2 (73.2–96.1%). BAR1 was phytoplankton-dominant (30.0–95.2%), with three samples that presented meaningful acdom contributions (acdom(655) > 50%), whilst BAR2 showed the predominance of CDOM (52.8–91.3%). IBI1 and IBI2 were CDOM-dominant (with averages of 48.5% and 71.4%, respectively), while there was one phytoplankton-dominated sample in IBI1 (aϕ(655)= 58.7%), and another sample was displaced to the phytoplankton apex in IBI2 (aϕ(655) = 82.6%). NAV1 presented a homogeneous distribution among all of the OACs, with a remarkable anap contribution (30.5% on average), and a slight highlight of the aϕ contribution in some samples, which reached 60%. NAV2 was mostly dominated by NAP (34.0–79.4%).

3.5. Light Absorption Variability from Upstream to Downstream

To better represent the longitudinal variability of light absorption from upstream to downstream, the average values of aϕ, anap, acdom, and at-w at 443 nm are displayed in boxplot graphs (Figure 7). We considered 443 nm as a reference wavelength, due to its relation to dissolved organic matter concentration [39] and the marked Chl-a pigment absorption feature [4,40]; furthermore, it is the central wavelength of the first OLI spectral band.
The light absorption decreases from upstream to downstream, following a similar pattern that was observed for OACs (see Figure 6). However, the absorption is very sensitive to the concentration of OACs at the time of the water sampling. There is an increase in the light absorption with higher concentrations and a decrease with lower concentrations. The acdom(443) data showed that the BAR reservoir has more CDOM in the water than the BB reservoir, although both reservoirs had similar Chl-a concentrations. The anap(443) showed the smallest longitudinal gradient of light absorption coefficients along the cascade. The at-w(443) decreased from upstream to downstream, and showed that phytoplankton was the main contributor to at-w(443).

4. Discussion

Evaluating optical properties throughout inland water systems that showed extreme upstream gradients is not an easy task [41]. There many variables that affect the hydrodynamic of hydroelectrical reservoirs, such as the hydraulic residence time of water (and fluctuation of water level), land use and cover use of drainage basin, wind, rainfall rates, runoff, and internal properties, including temperature, incident solar radiation into the water surface, and nutrients, for instance [26,42]. It is already known that the absorption properties are proportional to the amount and diversity of optically active constituents within the water [5,43]; however, the relations along cascading reservoirs with substantially different characteristics were not fully investigated in the literature.
All of the fieldwork was carried out in dry seasons (except IBI1) and at similar wind speeds. According to Tundisi and Matsumura-Tundisi (1990) [44], water management leads to relatively high residence times during the dry seasons, and dry seasons normally present higher concentrations of Chl-a and suspended matter, due to high nutrient loads and light availability [39]. The rainfall rates in BB did not vary too much when compared to the other fieldwork (Figure S5a) and presented high levels of Chl-a. The main cause can be attributed to the incident solar radiation onto the surface (Figure S5b) in BB2, which is not as predominant in IBI2 and BAR2 (for field studies that were carried out in June 2017).
The water residence times in each reservoir were: 30–180 days in BB; seven to 24 days in BAR; 25–41 days in IBI; and 46 days (on average) in NAV. Similar Chl-a concentrations for BAR1/IBI1 and BAR2/IBI2 are related to the low residence time similar length (BAR is almost 70 km long, and IBI is 48 km) and comparably meanders along their course. Therefore, this suggests that the same environmental effects affect both reservoirs, since they show similar physical characteristics. They are geographically near (they are in sequence within the cascade), and they are located in drainage basins with similar developed activities [45]. The higher Chl-a concentrations in the BAR and IBI reservoirs are a consequence of local discharges from respective tributaries. It is also reasonable to suggest that the high retention time and the proximity to urban areas, which are responsible for the high discharges of nutrients and optimal conditions for phytoplankton growth, are the main causes of the high Chl-a concentrations found in BB [4]. The low nutrient concentrations in NAV are responsible for the lack of phytoplankton development as much as bottom vegetation, due to river bottom nutrients, and the deeper euphotic zone [22,27].
The most variable qualitative parameters were PIM (in BB1, BB2, BAR2, NAV1, and NAV2 field studies) and Chl-a concentration (in BAR1 and IBI1 field studies), with outliers of Chl-a concentrations in BAR1 and IBI1 (Figure 2b). Outliers in BAR1 were the samples that were located in front of the Ribeirão Grande River, which is the tributary that receives wastewater discharges from Bauru city [31], which did not have domestic sewage treatment at the time of the fieldwork (Figure 1f). The outliers in IBI1 are located in the convergence zone between the Jacaré-Guaçu River and the Claro River (Figure 1g), which are the main tributaries of the IBI reservoir. These results suggest the existence of regional nutrient inputs from tributaries, which modify the regional characteristics within the BAR and IBI reservoirs. Temperature and pH did not present gradients along the TCSR (Table 1); however, higher pH values were accomplished by higher Chl-a concentrations, which can be caused by intense phytoplankton activity, which removes CO2 from aquatic systems, and changes the pH, to values normally higher than 7.0.
Regarding the water quality parameters in TCSR, the Chl-a concentrations are highly related to POM (r = 0.94), and POM is highly related to SPM (r = 0.99). The inorganic fraction (PIM) demonstrated a weak relationship with SPM (r = 0.67); such a decrease can be the consequence of the NAV reservoir, which is an inorganic-dominated environment [3]. The variations of organic-to-inorganic particulate matter from upstream to downstream reservoirs are depicted in Figure S1a, and demonstrated changes in the composition of the water quality parameters. In this respect, it is also possible to note a qualitative change with high concentrations of SPM and Chl-a in the BB and BAR reservoirs, and low concentrations in NAV and IBI (Figure S1a,b).
Based on the organic to inorganic changes throughout the cascade, the results suggested higher aϕ(λ) magnitudes in the upstream reservoirs, and lower aϕ(λ) magnitudes in the downstream reservoirs. Such a trend was observed when the absorption spectrum curves were evaluated (Figure 5a). As opposed to the phytoplankton and NAP absorption curves, the acdom spectra were investigated due to the CDOM absorption effects on the available light for primary production [5].
The sources of colored DOM can be either autochthonous resulting from phytoplankton degradation that increases CDOM concentrations, or allochthonous resulting from carbon matter carried by the surrounding terrestrial basin drainage [46,47,48]. Considering the amount of energy sampled during the field studies (Figure S5b), we can note a trend in the relationships between high incident solar radiation and low levels of CDOM (represented by acdom(443 nm)) in NAV and BB (both field studies). The inverse is also validated: low levels of incident solar radiation and high levels of acdom(443 nm) were observed in BAR (both field studies) and IBI2. This trend is the result of the photochemical degradation of carbon organic matter due to light exposure, resulting in CDOM decreases [41,49]. The exception, IBI1, presented high levels of CDOM accompanied by high levels of incident solar radiation, which indicates that the main source of carbon organic matter was high rainfall in July 2016 (Figure 1c). These high levels of CDOM could be the main factor that replaces photodegraded CDOM with terrestrial organic matter from the drainage basin and tributaries [47].
The acdom(443) from the entire cascade ranged from 0.21 m−1 to 3.17 m−1 (Table 2), which is similar to other published studies for aquatic systems [12,16,31,50,51,52]. The highest acdom magnitudes were found in BAR1 and IBI1 (sampled curve in Figure 5c). The highest acdom curve in IBI1 was obtained at the same sampling location that presented the highest Chl-a concentration in the IBI1 dataset. Additionally, this sample is located near the Environmental Protection Area of Ibitinga, which is an area with high forest coverage and bog lake formations, resulting from wetlands of tributaries (Jacaré-Guaçu and Jacaré-Pepira Rivers, Figure 1g). More detailed information can be found in Figure S6.
The exponential shape of acdom described in terms of Scdom showed a narrow range between 0.006 nm−1 < Scdom < 0.020 nm−1 in the TCSR. The Scdom showed an inverse proportional relation to its molecular weight, represented by acdom(443) (Kirk, 1994). The low values of Scdom implied a high molecular weight of CDOM, whilst high values of Scdom implied a low molecular weight [38,53]. Such an inverse relationship between Scdom and acdom(443) can be related to the CDOM photochemical degradations from phytoplankton [54,55].
With the aim of investigating the sources of CDOM, we plotted acdom(443) against Chl-a concentrations (Figure S3). The plots can indicate whether CDOM was derived from autochthonous material (phytoplankton photobleaching) or from allochthonous material (riparian vegetation loads and runoff discharges) [56]. The results showed that at least one field campaign from each reservoir presented a significant statistical relationship between Chl-a and acdom(443) (BB1, BAR2, IBI2, and NAV2). Nevertheless, the relationships changed between the field studies, which suggested that the CDOM origin varied. Moreover, the BB2 and BAR1 fieldworks did not show statistical relationships, even with high Chl-a concentrations, demonstrating that CDOM production is not dependent on the trophic status encountered within the aquatic system; BB and BAR are considered eutrophic environments.
The Chl-a concentrations, accessory pigments, and the different physiology states of phytoplankton assemblages are the main factors that modify the intensity and spectral position of features in aϕ(λ) [9,15,57,58]. Due to the existence of absorption peaks and their spectral positions, it is possible to infer the types of phytoplankton and pigments encountered in a sample [5,58].
The aϕ(λ) peaks near 620 nm are related to phycocyanin, which is a pigment that is commonly present in cyanobacteria [59,60]. BB and BAR exhibited such features in almost all of the aϕ(λ) curves. Previous studies have related the dominance of cyanobacteria to phytoplankton assemblages in such reservoirs [26,27]. At approximately 443 nm, the accessory pigments (Chl-b,c) and Chl-a are the major contributors to aϕ(λ), whereas peaks at 675 nm result from Chl-a plus phaeophytin pigments [57]. The features near 443 nm and 675 nm are remarkably observed in the aϕ(λ) curves obtained in BB (both field studies), BAR1, and IBI1; these features are less remarkable in NAV (both field studies), and disappear in BAR2 and IBI2 (Figure 4). A shoulder near 490 nm, which is encountered in some spectral curves in BB (both field studies), BAR1, and IBI1 indicates the presence of carotenoids in phytoplankton pigments [61]. The reduction in Chl-a concentrations directly affects the peaks near 440 nm and 675 nm in the aϕ curves (Figure 4), which are mainly in BAR1 and compared to those in BAR2 and IBI2. The total absorption is also influenced by such decay, and demonstrates the relationship between phytoplankton production and the spectral absorption of available energy. The feature near 440 nm represented the pigment accessories of phytoplankton (alpha and beta carotene), while the 675-nm feature represented Chl-a, b, or c pigments. The variability of features observed in samples from different reservoirs indicate the species richness of the phytoplankton community along the cascade. Due to the different absorption peaks in aϕ(λ) from upstream to downstream caused by the high variability in Chl-a concentrations and package effects [62], it is not possible to categorize the reservoirs according to the dominance of certain types of phytoplankton assemblages per reservoir. Conversely, it is possible to identify a longitudinal gradient of aϕ(λ) magnitudes in TCSR (Figure 5a), as observed for Chl-a concentrations (Figure 2).
The anap(λ) curves showed the iron oxides features between 450–550 nm [63] in BAR, IBI, and NAV fieldwork (Figure 4), which is consistent with the sandy soil formation in the drainage basin in the BAR and IBI reservoirs [64,65], and demonstrated the influence of sand extraction activities observed in the NAV reservoir during the field studies. The highest and lowest anap(λ) average curves (Figure 5b) are related to the average PIM concentrations (Table 1), as observed in BAR1 (2.30 mg·L−1) and IBI2 (0.40 mg·L−1), respectively. The anap(λ) magnitudes depend on the concentration and type of particles [66,67]. Despite the significant differences in Snap in the fieldwork, these values are too small to hypothesize a water classification as a function of Snap, such as the classification made in African reservoirs [16].
The evaluation of the absorption contribution of each OAC to at-w made by a ternary diagram allowed for an understanding of the optical features in each reservoir. The ternary diagram is considered to be a tool for water classification, due to sampling distributions along the three axes (Case 1 or Case 2-type [68]). It is not our intention to classify the waters of the TCSR, but rather to properly understand what happens in the absorption coefficients along the cascade. Overall, the results indicated that the BAR and IBI reservoirs are CDOM-dominated environments, except for samples with high Chl-a concentrations (higher than 700 mg·m−3). The CDOM-dominance was conserved through the spectra (at 443 nm, 482 nm, 560 nm, and 655 nm) for the BAR and IBI fieldwork. BB (both field studies) is a phytoplankton-dominated environment, whilst NAV (both field studies) was identified as an NAP-dominated environment. Based on the average amount of contributions from each field, the reservoirs showed specific OAC dominances that were not wavelength dependent.
Considering 443 nm as a reference wavelength for plotting the absorption coefficient variation along the cascade (Figure 7), we observed a slight decrease in absorption from upstream to downstream, mainly in at-w(443). Furthermore, phytoplankton absorption acted as the major contributor of at-w(443), and showed the same decreasing trend encountered for OAC concentrations (Figure 2). The increasing dominance of CDOM in the BAR and IBI (Figure 6c) showed a relevant enrichment of CDOM within these reservoirs. Considering that CDOM is the colored fraction of dissolved organic matter (DOM), and DOM is an important fraction of dissolved organic carbon (DOC) in inland waters, CDOM can be used as a useful indicator of DOM and DOC [69,70]. Therefore, the CDOM dominance in BAR and IBI indicated high DOM and DOC concentrations, caused by internal production, in BAR2 and IBI2 (Figure S3), and terrestrial sources as in IBI1, due to the runoff discharges driven by watershed contributions.

5. Conclusions

Understanding and interpreting bio-optical properties, such as IOPs, further the understanding of the underwater light field, and are essential to bio-optical modeling developments and accurate water quality monitoring applications. The present results allowed an analysis of the dynamics of a cascading reservoir system, showing the continuous longitudinal gradient of optical absorptions and OACs in TCSR. Variations of the magnitude of OACs affect the light absorption coefficients, mainly aϕ and at-w. From upstream to downstream, the TCSR presented a slight decrease for the NAP and CDOM absorption coefficients, and they showed more susceptibility to allochthonous effects.
acdom was the highly variable parameter within the cascade, mainly in the BAR and IBI reservoirs. The acdom spectra did not show a decaying trend along TCSR; instead, it presented low magnitudes in BB and NAV, and major magnitudes in IBI and BAR. The different magnitudes of the acdom curves among the samples collected within the same reservoir indicate additional sources of CDOM throughout the aquatic systems; in other words, the tributaries and activities in the nearby areas affected the optical properties within reservoirs as much as the water flow from previous reservoirs. NAP absorption presented a less remarkable decaying trend from upstream to downstream, with high variations in the BAR and IBI fieldwork that can result from the contribution of runoff sediments.
The longitudinal gradient of aϕ is proportional to the Chl-a concentration throughout the cascade, and directly influences at-w, despite the individual optical dominance for each reservoir. BB is a phytoplankton-dominated environment; BAR and IBI are CDOM-dominant environments, and NAV is an NAP-dominant environment. Some exceptions to the dominance of OAC were found in the specific sampling locations that showed distinct Chl-a, SPM, or acdom values compared to the entire dataset. In these cases, the highest absorption contribution within the reservoir changed accordingly to the most expressive OAC.
The insights from this work demonstrated a longitudinal gradient of absorption coefficients that were directly affected by OAC variations throughout the TCSR. Moreover, CDOM and NAP are essential to watershed contributions as much as the internal characteristics (such as residence time, incident solar radiation, and so forth). Furthermore, the land use and land cover in watersheds provide nutrient enrichments and pollution discharges within the reservoirs, and are reflected by OACs in specific locations, which affect the absorption coefficients. These local variations highlight the needs of local management to identify the main sources of higher OAC concentrations. The absorption variability found in TCSR confirms the challenges in retrieving OAC concentrations using a unique bio-optical model for the entire cascade. In contrast, our results highlighted the presence of a total absorption gradient and the strong influence of phytoplankton absorption on the total absorption along the aquatic continuum. Further analysis in a longer time series might provide temporal absorption coefficient datasets in TCSR and support strategies for bio-optical modeling in different optical conditions.

Supplementary Materials

The following are available online at https://www.mdpi.com/2073-4441/11/2/229/s1.

Author Contributions

Conceptualization, E.A. and N.B.; methodology, E.A., N.B., F.W., T.R.; software, A.d.C.; validation, N.B., A.C.C.G. and C.A.; formal analysis, N.B., A.C.C.G. and T.R.; investigation, N.B.; resources, E.A.; data curation, N.B., F.W., T.R., A.d.C., A.C.C.G. and C.A.; writing—N.B.; writing—review and editing, E.A., F.W., A.d.C. and C.A.; visualization, N.B.; supervision, funding acquisition and project administration, E.A.

Funding

This research was funded by São Paulo Research Foundation - FAPESP, grant numbers 2012/19821-1 and 2015/21586-9; and National Counsel of Technological and Scientific Development—CNPq (Project numbers: 400881/2013-6 and 472131/2012-5).

Acknowledgments

The authors thank to FAPESP Projects aforementioned. The authors also thank the Coordination for the Improvement of Higher Education Personnel (CAPES) for scholarship, and Professor Edivaldo D. Velini and staff from FCA/UNESP for allowing the use of their laboratory facilities.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Lee, Z.; Carder, K.L.; Arnone, R.A. Deriving inherent optical properties from water color: A multiband quasi-analytical algorithm for optically deep waters. Appl. Opt. 2002, 41, 5755–5772. [Google Scholar] [CrossRef] [PubMed]
  2. Odermatt, D.; Gitelson, A.; Brando, V.E.; Schaepman, M. Review of constituent retrieval in optically deep and complex waters from satellite imagery. Remote Sens. Environ. 2012, 118, 116–126. [Google Scholar] [CrossRef] [Green Version]
  3. Rodrigues, T.; Alcântara, E.; Watanabe, F.; Bernardo, N.; Rotta, L.; Imai, N. Spatial and temporal variations of the inherent optical properties in a tropical cascading reservoir system. Model. Earth Syst. Environ. 2016, 2, 86. [Google Scholar] [CrossRef]
  4. Watanabe, F.; Mishra, D.R.; Astuti, I.; Rodrigues, T.; Alcântara, E.; Imai, N.; Barbosa, C. Parametrization and calibration of a quasi-analytical algorithm for tropical eutrophic waters. ISPRS J. Photogramm. Remote Sens. 2016, 12, 28–47. [Google Scholar] [CrossRef]
  5. Kirk, J. Light and Photosysnthesis in Aquactic Ecosystems, 2nd ed.; Cambridge University Press: Cambridge, UK, 1994. [Google Scholar]
  6. Bricaud, A.; Morel, A.; Louis, P. Absorption by dissolved organic matter of the sea (yellow substance) in the UV and visible domains. Limnol. Oceanogr. 1981, 26, 43–53. [Google Scholar] [CrossRef]
  7. Bricaud, A.; Babin, M.; Morel, A.; Claustre, H. Variability in the chlorophyll-specific absorption coefficients of natural phytoplankton: Analysis and parameterization. J. Geophys. Res. 1995, 100, 13321–13332. [Google Scholar] [CrossRef]
  8. Bricaud, A.; Morel, A.; Babin, M.; Allali, K.; Claustre, H. Variations of light absorption by suspended particles with chlorophyll a concentration in oceanic (case 1) waters: Analysis and implications for bio-optical models. J. Geophys. Res. 1998, 103, 31033–31044. [Google Scholar] [CrossRef]
  9. Babin, M.; Stramski, D.; Ferrari, G.M.; Claustre, H.; Bricaud, A.; Obolensky, G.; Hoepffner, N. Variations in the light absorption coefficients of phytoplankton, nonalgal particles, and dissolved organic matter in coastal waters around Europe. J. Geophys. Res. 2003, 108. [Google Scholar] [CrossRef] [Green Version]
  10. Gokul, E.A.; Palanisamy, S.; Sundarabalan, B.; Sahay, A.; Chauhan, P. Modelling the inherent optical properties and estimating the constituents’ concentrations in turbid and eutrophic waters. Cont. Shelf Res. 2014, 84, 120–138. [Google Scholar] [CrossRef]
  11. Shaju, S.S.; Minu, P.; Srokanth, A.S.; Ashraf, P.M.; Vijayan, A.K.; Meenakumari, B. Decomposition study of in vivo phytoplankton absorption spectra to identify the pigments and phytoplankton group in complex case 2 waters of coastal Arabian Sea. Oceanol. Hydrobiol. Stud. 2017, 44, 282–293. [Google Scholar] [CrossRef]
  12. Vishnu, P.S.; Shaju, S.S.; Tiwari, S.P.; Menon, N.; Nashad, M.; Joseph, C.A.; Raman, M.; Hatha, M.; Prabhakaran, M.P.; Mohandas, A. Seasonal variability in bio-optical properties along the coastal waters of Cochin. Int. J. Appl. Earth Obs. Geoinf. 2018, 66, 184–195. [Google Scholar] [CrossRef]
  13. Das, S.; Hazra, S.; Giri, S.; Das, I.; Chanda, A.; Akhand, A.; Maity, S. Light absorption characteristics of chromophoric dissolved organic matter (CDOM) in the coastal waters of northern Bay of Bengal during winter season. Indian J. Geo Mar. Sci. 2017, 46, 884–992. [Google Scholar]
  14. Ferreira, A.; Ciotti, A.M.; Giannini, M.F.C. Variability in the light absorption coefficients of phytoplankton, non-algae particles, and colored dissolved organic matter in a subtropical bay (Brazil). Estuar. Coast. Shelf Sci. 2014, 139, 127–136. [Google Scholar] [CrossRef]
  15. Campbell, G.; Phinn, S.R.; Daniel, P. The specific inherent optical properties of three subtropical and tropical water reservoirs in Queensland, Australia. Hydrobiologia 2011, 658, 233–252. [Google Scholar] [CrossRef]
  16. Matthews, M.W.; Bernard, S. Characterizing the absorption properties for remote sensing of three small optically-diverse South African reservoirs. Remote Sens. 2013, 5, 4370–4404. [Google Scholar] [CrossRef]
  17. Alcântara, E.; Watanabe, F.; Rodrigues, T.; Bernardo, N.; Rotta, L.; Carmo, A.; Curtarelli, M.; Imai, N. Field measurements of the backscattering coefficient in a cascading reservoir system: First results from Nova Avanhandava and Barra Bonita Reservoirs (São Paulo, Brazil). Remote Sens. Lett. 2016, 7, 417–426. [Google Scholar] [CrossRef]
  18. Moore, T.; Dowell, M.D.; Bradt, S.; Verdu, A.R. An optical water type framework for selecting and blending retrievals from bio-optical algorithms in lake and coastal waters. Remote Sens. Environ. 2014, 143, 97–111. [Google Scholar] [CrossRef]
  19. Shen, Q.; Li, J.; Zhang, F.; Li, J.; Li, W.; Zhang, B. Classification of several optically complex waters in China using in situ remote sensing reflectance. Remote Sens. 2015, 7, 429–440. [Google Scholar] [CrossRef]
  20. Spyrakos, E.; O’Donnell, R.; Hunter, P.D.; Miller, C.; Scott, M.; Simis, S.G.; Neil, C.; Barbosa, C.C.; Binding, C.E.; Bradt, S.; et al. Optical types of inland and coastal waters. Limnol. Oceanogr. 2018, 63, 846–870. [Google Scholar] [CrossRef]
  21. Pahlevan, N.; Lee, Z.; Wei, J.; Schaaf, C.B.; Schott, J.R.; Berk, A. On-orbit radiometric characterizations of OLI (Landast-8) for applications in aquatic remote sensing. Remote Sens. Environ. 2014, 154, 272–284. [Google Scholar] [CrossRef]
  22. Rodrigues, T.; Alcântara, E.; Watanabe, F.; Imai, N. Retrieval of Secchi Disk Depth from reservoir using a semi-analytical scheme. Remote Sens. Environ. 2017, 198, 213–228. [Google Scholar] [CrossRef]
  23. Muow, C.B.; Greb, S.; Aurin, D.; DiGiacomo, P.M.; Lee, Z.; Twardowski, M.T.; Binding, C.; Hu, C.; Ma, R.; Moore, T.; et al. Aquatic color radiometry remote sensing of coastal and inland waters: Challenges and recommendations for future satellite missions. Remote Sens. Environ. 2013, 160, 15–30. [Google Scholar] [CrossRef]
  24. Palmer, S.C.J.; Kutser, T.; Hunter, P.D. Remote sensing of inland waters: Challenges, progress and future directions. Remote Sens. Environ. 2015, 157, 1–7. [Google Scholar] [CrossRef]
  25. Shi, L.; Mao, Z.; Wu, J.; Liu, M.; Zhang, Y.; Wang, Z. Variations in spectral absorption properties of phytoplankton, non-algae particles and Chromophoric Dissolved Organic Matter in Lake Qiandaohu. Water 2018, 9, 352. [Google Scholar] [CrossRef]
  26. Barbosa, F.A.R.; Padisák, J.; Espíndola, E.L.G.; Borics, G.; Rocha, O. The cascading reservoirs continuum concept (CRCC) and its application to the river Tietê-basin, São Paulo State, Brazil. In Theoretical Reservoir Ecology and its Application; Tundisi, J.G., Straskraba, M., Eds.; International Institute of Ecology, Academy of Sciences and Backhuys Publishers: São Carlos, Brazil, 1999; pp. 425–437. [Google Scholar]
  27. Smith, W.S.; Espíndola, E.L.G.; Rocha, O. Environmental gradients in reservoirs of the medium and low Tietê River: Limnological differences trough the habitat sequence. Acta Limnol. Bras. 2014, 16, 73–88. [Google Scholar] [CrossRef]
  28. Frascareli, D.; Cardoso-Silva, S.; Mizael, J.D.; Rosa, A.H.; Pompeo, M.L.M.; Doval, J.C.L.; Moschini, C.V. Spatial distribution, bioavailability, and toxicity of metals in surface sediments of tropical reservoirs, Brazil. Environ. Monit. Assess 2018, 190, 199. [Google Scholar] [CrossRef] [PubMed]
  29. Dellamano-oliveira, M.J.; Vieira, A.A.H.; Rocha, O.; Colombo, V.; Sant’Anna, C.L. Phytoplankton taxonomic composition and temporal changes in a tropical reservoir. Fund Appl. Limnol. 2008, 171, 27–38. [Google Scholar] [CrossRef]
  30. Londe, L.R.; Novo, E.M.L.M.; Barbosa, C.; Araujo, C.A.S. Water residence time affecting phytoplankton blooms: Study case in Ibitinga Reservoir (São Paulo, Brazil) using Landsat/TM images. Braz. J. Biol. 2016, 76, 664–672. [Google Scholar] [CrossRef] [PubMed]
  31. Cairo, C.T.; Barbosa, C.C.F.; Novo, E.M.L.M.; Calijuri, M.C. Spatial and seasonal variation in diffuse attenuation coefficients of downward irradiance at Ibitinga Reservoir, São Paulo, Brazil. Hydrobiologia 2017, 784, 265–282. [Google Scholar] [CrossRef]
  32. Rodrigues, T.W.P.; Guimarães, U.S.; Rotta, L.H.; Watanabe, F.S.; Alcântara, E.; Imai, N.N. Sampling design in reservoirs based on Landsat-8/OLI images: A case study in Nova Avanhandava reservoir (São Paulo State, Brazil). Bol. Ciênc. Geod. 2016, 22, 304–323. [Google Scholar] [CrossRef]
  33. APHA (American Public Health Association); AWWA (American Water Works Association); WEF (Water Environmental Federation). Standard Methods for the Examination of Water and Wastewater, 20th ed.; APHA; AWWA; WEF: Washington, DC, USA, 1998. [Google Scholar]
  34. Golterman, H.L. Developments in Water Science 2. Physiological Limnology: An Approach to the Physiology of Lake Ecosystems; Elsevier: Amsterdam, The Netherlands, 1975. [Google Scholar]
  35. Smith, R.C.; Baker, K.S. Optical properties of the clearest natural waters (200–800 nm). Appl. Opt. 1981, 20, 177–184. [Google Scholar] [CrossRef] [PubMed]
  36. Tassan, S.; Ferrari, G.M. An alternative approach to absorption measurement of aquatic particles retained on filters. Limnol. Oceanogr. 1995, 40, 1358–1368. [Google Scholar] [CrossRef]
  37. Tassan, S.; Ferrari, G.M. Measurement of light absorption by aquatic particles retained on filters: Determination of the optical path length amplification by the ‘transmittance-reflectance’ method. J. Plankton Res. 1998, 20, 1699–1709. [Google Scholar] [CrossRef]
  38. Twardowski, M.S.; Boss, E.; Sullivan, J.M.; Donaghay, P.L. Modeling the spectral shape of absorption by chromophoric dissolved organic matter. Mar. Chem. 2004, 89, 69–88. [Google Scholar] [CrossRef]
  39. Alcântara, E.; Bernardo, N.; Watanabe, F.; Rodrigues, T.; Rotta, L.; Carmo, A.; Shimabukuro, M.; Gonçalves, S.; Imai, N. Estimating the CDOM absorption coefficient in tropical inland waters using OLI/Landsat-8 images. Remote Sens. Lett. 2016, 7, 661–670. [Google Scholar] [CrossRef]
  40. Cetinic, I. Phytoplankton: Optical Constituents of the Ocean. Available online: http://www.oceanopticsbook.info/view/optical_constituents_of_the_ocean/_phytoplankton#searchResult1 (accessed on 5 September 2018).
  41. Massicote, P.; Asmala, E.; Stedmon, C.; Markager, S. Global distribution of dissolved organic matter along the aquatic continuum: Across rivers, lakes and oceans. Sci. Total Environ. 2017, 609, 180–191. [Google Scholar] [CrossRef] [PubMed]
  42. Straškraba, M. Retention time as a key variable of reservoir limnology. In Theoretical Reservoir Ecology and Its Applications; Tundisi, J.G., Straškraba, M., Eds.; International Institute of Ecology, Brazilian Academy of Sciences and Backhuys Publishers: São Carlos, Brazil, 1999; pp. 385–410. [Google Scholar]
  43. Bukata, R.P.; Jerome, J.H.; Kondratyev, A.S.; Pozdnyakov, D.V. Optical Properties and Remote Sensing of Inland and Coastal Waters; CRC Press: Boca Raton, FL, USA, 1995. [Google Scholar]
  44. Tundisi, J.G.; Matsumura-Tundisi, T. Limnology and eutrophication of Barra Bonita reservoir, S. Paulo State, Southern Brazil. Arch. Hydrobiol. Beih. 1990, 33, 661–676. [Google Scholar]
  45. Periotto, N.A.; Tundisi, J.G. A characterization of ecosystems services, drives and values of two watersheds in São Paulo State, Brazil. Braz. J. Biol. 2018, 78, 397–407. [Google Scholar] [CrossRef] [PubMed]
  46. Massicotte, P.; Frenette, J.-J. Spatial connectivity in a large river system: Resolving the sources and fate of dissolved organic matter. Ecol. Appl. 2011, 21, 2600–2617. [Google Scholar] [CrossRef]
  47. Lambert, T.; Teodoru, C.R.; Nyoni, F.C.; Bouillon, S.; Darchambeau, F.; Massicotte, P.; Borges, A.V. Along-stream transport and transformation of dissolved organic matter in a large tropical river. Biogeosciences 2016, 13, 2727–2741. [Google Scholar] [CrossRef] [Green Version]
  48. Martins, S.; Chokmani, K.; Alcântara, E.; Ogashawara, I.; El-Alem, A. Mapping the coloured dissolved organic matter absorption coefficient in a eutrophic reservoir using remotely sensed images. Inland Waters 2018, 8, 488–504. [Google Scholar] [CrossRef]
  49. Vähätalo, A.V.; Wetzel, R.G. Photochemical and microbial decomposition of chromophoric dissolved organic matter during long (months–years) exposures. Mar. Chem. 2004, 89, 313–326. [Google Scholar] [CrossRef]
  50. Zhang, Y.; Zhang, B.; Wang, X.; Li, J.; Feng, S.; Zhao, Q.; Liu, M.; Qin, B. A study of absorption characteristics of chromophoric dissolved organic matter and particles in Lake Taihu, China. Hydrobiologia 2007, 592, 105–120. [Google Scholar] [CrossRef]
  51. Zhang, M.; Tang, J.; Dong, Q.; Song, Q.; Ding, J. Retrieval of total suspended matter concentration in the Yellow and East China Seas from MODIS imagery. Remote Sens. Envron. 2010, 114, 392–403. [Google Scholar] [CrossRef]
  52. Riddick, C.A.L.; Hunter, P.D.; Tyler, A.N.; Vicente, V.M.; Horváth, H.; Kovács, A.W.; Vörös, L.; Preston, T.; Présing, M. Spatial Variability of absorption coefficient over a biogeochemical gradient in a large and optically complex shallow lake. J. Geophys. Res. 2015, 120, 7040–7066. [Google Scholar] [CrossRef]
  53. Carder, K.L.; Steward, R.G.; Harvey, G.R.; Ortner, P.B. Marine humic and fulvic acids: Their effects on remote sensing of ocean chlorophyll. Limnol. Oceanogr. 1989, 34, 68–81. [Google Scholar] [CrossRef] [Green Version]
  54. Helms, J.R.; Stubbins, A.; Ritchie, J.D.; Minor, E.C.; Kieber, D.J.; Mopper, K. Absorption spectral slopes and slope ratios as indicators of molecular weight, source, and photobleaching of chromophoric dissolved organic matter. Limnol. Oceanogr 2008, 53, 955–969. [Google Scholar] [CrossRef] [Green Version]
  55. Matsuoka, A.; Hooker, S.B.; Bricaud, A.; Babin, M. Estimating absorption coefficients of colored dissolved organic matter (CDOM) using a semi-analytical algorithm for southern Beaufort Sea Waters: Application to deriving concentrations of dissolved organic carbon from space. Biogeosciences 2013, 10, 917–927. [Google Scholar] [CrossRef]
  56. Tilstone, G.H.; Lotliker, A.A.; Miller, P.I.; Ashraf, P.M.; Kumar, T.S.; Suresh, T.; Ragavan, B.R.; Menon, H.B. Assessment of MODIS-Aqua chlorophyll-a algorithms in coastal and shelf waters of the eastern Arabian Sea. Cont. Shelf Res. 2013, 65, 14–26. [Google Scholar] [CrossRef]
  57. Sathyendranath, S.; Lazzara, L.; Prieur, L. Variations in the spectral values of specific absorption of phytoplankton. Limnol. Oceanogr. 1987, 32, 403–415. [Google Scholar] [CrossRef] [Green Version]
  58. Mobley, C.D. Light and Water: Radiative Transfer in Natural Water; Academic Press: San Diego, CA, USA, 1994. [Google Scholar]
  59. Richardson, L.L. Remote sensing of algal bloom dynamics: A new research fuses remote sensing of aquatic ecosystems with algal accessory pigment analysis. Bioscience 1996, 46, 492–501. [Google Scholar] [CrossRef]
  60. Mishra, S.; Mishra, D.R.; Lee, Z. Bio-optical inversion in highly turbid and cyanobacteria-dominated waters. IEEE Trans. Geosci. Remote 2014, 52, 375–388. [Google Scholar] [CrossRef]
  61. Perkins, M.G.; Effler, S.W.; Strait, C.M. Phytoplankton absorption and the chlorophyll a-specific absorption coefficient in dynamic Onondaga Lake. Inland Waters 2014, 4, 133–146. [Google Scholar] [CrossRef]
  62. Watanabe, F.S.Y.; Alcântara, E.H.; Rodrigues, T.W.P.; Bernardo, N.M.R.; Rotta, L.H.S.; Imai, N.N. Phytoplankton community dynamic detection from the chlorophyll-specific absorption coefficient in productive inland waters. Acta Limnol. Bras. 2017, 29, s2179–s2975. [Google Scholar] [CrossRef]
  63. Estapa, M.L.; Boss, E.; Mayer, L.M.; Roesler, C.S. Role of iron and organic carbon in mass-specific light absorption by particulate matter from Louisiana coastal waters. Limnol. Oceanogr 2012, 57, 97–112. [Google Scholar] [CrossRef]
  64. Souza, A.D.G.; Tundisi, J.G. Hidrogeochemical comparative study of the Jaú and Jacaré-Guaçu River watersheds, São Paulo, Brazil. (English). Rev. Bras. Biol. 2000, 60, 563–570. [Google Scholar] [CrossRef]
  65. Tundisi, J.D.; Matsumura-Tundisi, T.; Pareschi, D.C.; Luzia, A.P.; von Haeling, P.H.; Frollini, E. The Tietê/Jacaré watershed: A case study in research and management. Estud. Av. 2008, 22, 159–172. [Google Scholar] [CrossRef]
  66. Woźniak, S.B.; Meler, J.; Lednicka, B.; Zdun, A.; Stoń-Egiert, J. Inherent optical properties of suspended particulate matter in the southern Baltic Sea. Oceanologia 2011, 53, 691–729. [Google Scholar] [CrossRef]
  67. Meler, J.; Ostrowska, M.; Ston-Egiert, J.; Zablocka, M. Seasonal and spatial variability of light absorption by suspended particles in the southern Baltic: A mathematical description. J. Mar. Syst. 2017, 170, 68–87. [Google Scholar] [CrossRef]
  68. Prieur, L.; Sathyendranath, S. An optical classification of coastal and oceanic waters based on specific spectral absorption curves of phytoplankton pigments, dissolved organic matter, and other particulate materials. Limnol. Oceanogr. 1981, 26, 671–689. [Google Scholar] [CrossRef]
  69. Cole, J.J.; Prairie, Y.T.; Caraco, N.F.; McDowell, W.H.; Tranvik, L.J.; Striegl, R.G.; Duarte, C.M.; Kortelainen, P.; Downing, J.A.; Middelburg, J.J.; Melack, J. Plumbing the Global Carbon Cycle: Integrating inland waters into the terrestrial carbon budget. Ecosystems 2007, 10, 175–185. [Google Scholar] [CrossRef]
  70. Zhu, W.; Yu, Q.; Tian, Y.Q.; Becker, B.L.; Zheng, T.; Carrick, H.J. An assessment of remote sensing algorithms for colored dissolved organic matter in complex freshwater environments. Remote Sens. Environ. 2014, 140, 776–778. [Google Scholar] [CrossRef]
Figure 1. Brazil with São Paulo State location (a); Map of five reservoirs in TCSR (note that Três Irmãos is located after the Nova Avanhandava reservoir) and respective land cover classes of the Tietê Basin (water body, forest, shrubland, bare soil, wetlands, and urban areas by Coordenadoria de Planejamento Ambiental de São Paulo, CPLA, 2010) (b). Monthly average rainfall data from NASA’s (National Aeronautics and Space Administration) Giovanni database for the period of 2014–2017 (c); and the sampling locations detailed for each reservoir: NAV1 (n = 20) and NAV2 (n = 20) (d); BB1 (n = 20) and BB2 (n = 30) (e); IBI1 (n = 10) and IBI2 (n = 20) (f); and BAR1 (n = 30) and BAR2 (n = 18) (g). All maps are in WGS84. NAV: Nova-Avanhandava, BB: Barra Bonita, IBI: Ibitinga, and BAR: Bariri.
Figure 1. Brazil with São Paulo State location (a); Map of five reservoirs in TCSR (note that Três Irmãos is located after the Nova Avanhandava reservoir) and respective land cover classes of the Tietê Basin (water body, forest, shrubland, bare soil, wetlands, and urban areas by Coordenadoria de Planejamento Ambiental de São Paulo, CPLA, 2010) (b). Monthly average rainfall data from NASA’s (National Aeronautics and Space Administration) Giovanni database for the period of 2014–2017 (c); and the sampling locations detailed for each reservoir: NAV1 (n = 20) and NAV2 (n = 20) (d); BB1 (n = 20) and BB2 (n = 30) (e); IBI1 (n = 10) and IBI2 (n = 20) (f); and BAR1 (n = 30) and BAR2 (n = 18) (g). All maps are in WGS84. NAV: Nova-Avanhandava, BB: Barra Bonita, IBI: Ibitinga, and BAR: Bariri.
Water 11 00229 g001
Figure 2. Boxplot of water quality parameters obtained in each field campaigns: (a) SPM concentration (mg·L−1); (b) Chl-a concentrations (mg·m−3); (c) Turbidity (NTU); (d) Secchi disk depth (ZSD, meters).
Figure 2. Boxplot of water quality parameters obtained in each field campaigns: (a) SPM concentration (mg·L−1); (b) Chl-a concentrations (mg·m−3); (c) Turbidity (NTU); (d) Secchi disk depth (ZSD, meters).
Water 11 00229 g002
Figure 3. In situ absorption spectra from all of the sampling locations for (a) acdom; (b) ap; (c) anap; and (d) aϕ.
Figure 3. In situ absorption spectra from all of the sampling locations for (a) acdom; (b) ap; (c) anap; and (d) aϕ.
Water 11 00229 g003
Figure 4. Average absorption coefficients by CDOM (acdom, high dashed line), phytoplankton (aϕ, dotted line), non-algae particles (anap, low dashed line), and total non-water absorption (at-w, continuous line) for the (a) BB1; (b) BB2; (c) BAR1; (d) BAR2; (e) IBI1; (f) IBI2; (g) NAV1; and (h) NAV2. Note the differences in the Y-axis between the surveys.
Figure 4. Average absorption coefficients by CDOM (acdom, high dashed line), phytoplankton (aϕ, dotted line), non-algae particles (anap, low dashed line), and total non-water absorption (at-w, continuous line) for the (a) BB1; (b) BB2; (c) BAR1; (d) BAR2; (e) IBI1; (f) IBI2; (g) NAV1; and (h) NAV2. Note the differences in the Y-axis between the surveys.
Water 11 00229 g004
Figure 5. Average absorption spectra of (a) aϕ; (b) anap; (c) acdom from all of the field surveys and respective optically active compounds (OACs).
Figure 5. Average absorption spectra of (a) aϕ; (b) anap; (c) acdom from all of the field surveys and respective optically active compounds (OACs).
Water 11 00229 g005
Figure 6. Ternary diagrams showing the contribution (%) of absorption coefficients to at-w at (a) 443 nm; (b) 482 nm; (c) 560 nm; and (d) 655 nm from all of the field surveys (BB1, BB2, BAR1, BAR2, IBI1, IBI2, NAV1, and NAV2).
Figure 6. Ternary diagrams showing the contribution (%) of absorption coefficients to at-w at (a) 443 nm; (b) 482 nm; (c) 560 nm; and (d) 655 nm from all of the field surveys (BB1, BB2, BAR1, BAR2, IBI1, IBI2, NAV1, and NAV2).
Water 11 00229 g006
Figure 7. Boxplot of absorption coefficients at 443 nm obtained in all of the field campaigns: (a) aϕ(443)(m−1); (b) anap(443)(m−1); (c) acdom(443)(m−1); and (d) at-w(443)(m−1).
Figure 7. Boxplot of absorption coefficients at 443 nm obtained in all of the field campaigns: (a) aϕ(443)(m−1); (b) anap(443)(m−1); (c) acdom(443)(m−1); and (d) at-w(443)(m−1).
Water 11 00229 g007
Table 1. Descriptive statistics of water quality parameters in BB, BAR, IBI, and NAV from two field campaigns. Here Min, Max, SD, and CV represent the minimum, maximum, standard deviation, and coefficient of variation (in %), respectively. Zsd is the Secchi Disk Depth (m), SPM is suspended particle matter (mg·L−1), PIM is particulate inorganic matter (mg·L−1), POM is particulate organic matter (mg·L−1), Chl-a is Chlorophyll-a concentration (mg·m−3), and Temp is temperature (°C).
Table 1. Descriptive statistics of water quality parameters in BB, BAR, IBI, and NAV from two field campaigns. Here Min, Max, SD, and CV represent the minimum, maximum, standard deviation, and coefficient of variation (in %), respectively. Zsd is the Secchi Disk Depth (m), SPM is suspended particle matter (mg·L−1), PIM is particulate inorganic matter (mg·L−1), POM is particulate organic matter (mg·L−1), Chl-a is Chlorophyll-a concentration (mg·m−3), and Temp is temperature (°C).
Zsd (m)Turbidity (NTU)pHChl-a (mg·m³)SPM (mg·L−1)PIM (mg·L−1)POM (mg·L−1)Temp. (°C)
Barra Bonita 1 (n = 20)—BB1
Min0.801.667.1817.753.600.202.8024.5
Max2.3012.509.25279.8616.304.4014.726.9
Mean1.495.178.36120.447.211.106.1025.6
CV28.947.08.3258.445.2178.852.02.8
Barra Bonita 2 (n = 20)—BB2
Min0.3711.607.12263.2010.800.6010.224.5
Max0.7833.2010.1797.8044.003.8030.432.1
Mean0.5718.649.28428.7221.972.6018.228.1
CV17.1828.269.4436.0332.0537.3026.27.8
Bariri (n = 30)—BAR1
Min0.507.806.1025.673.600.901.421.1
Max1.6080.909.90709.8940.334.0036.339.4
Mean1.1616.607.90119.768.282.305.924.3
CV20.0345.8210.580.5254.7621.475.115.4
Bariri (n = 18)—BAR2
Min1.603.506.833.800.200.200.4022.0
Max3.208.807.2819.02.601.301.6023.9
Mean2.205.76.978.001.600.601.1022.8
CV10.921.91.9040.927.9042.428.81.90
Ibitinga (n = 30)—IBI1
Min1.602.825.501.371.000.300.5021.2
Max3.208.877.00119.048.102.606.0030.1
Mean2.234.296.1021.752.610.801.8023.7
CV10.9117.906.8085.9739.2035.349.69.50
Ibitinga (n = 16)—IBI2
Min1.901.856.502.500.200.200.3021.40
Max3.803.606.9613.72.201.001.9024.00
Mean2.902.476.786.641.060.400.9322.82
CV19.521.11.7767.253.561.849.83.80
Nova Avanhandava 1 (n = 20)—NAV1
Min2.291.018.502.460.100.100.2025.1
Max4.802.478.9012.562.602.200.9026.3
Mean3.151.668.606.211.010.700.9326.0
CV19.9525.401.3940.061.776.740.81.12
Nova Avanhandava 2 (n = 20)—NAV2
Min2.451.017.604.510.500.140.3023.8
Max4.652.568.3020.52.802.001.1025.6
Mean3.411.738.109.011.000.500.5024.6
CV14.118.982.2034.937.665.826.51.90
Entire Dataset (n = 174 *)
Min0.371.015.501.370.100.080.2021.1
Max4.8080.910.1797.8044.004.4043.0039.4
Mean1.1513.318.36197.5011.341.949.4025.8
CV88.767.013.980.069.753.181.010.8
* Some samples and its respective water quality parameters were not evaluated.
Table 2. Descriptive statistics of ap(443) in m−1, aϕ(443), anap(443), acdom(443), non-algae particles (NAP) Slope (Snap), and colored dissolved organic matter (CDOM) Slope (Scdom). Where Max, Min, SD, and CV are the maximum, minimum, standard deviation, and coefficient of variation (in %). The exponential slope was obtained from CDOM within the 400–700 nm of range.
Table 2. Descriptive statistics of ap(443) in m−1, aϕ(443), anap(443), acdom(443), non-algae particles (NAP) Slope (Snap), and colored dissolved organic matter (CDOM) Slope (Scdom). Where Max, Min, SD, and CV are the maximum, minimum, standard deviation, and coefficient of variation (in %). The exponential slope was obtained from CDOM within the 400–700 nm of range.
Field Sites
BB1BB2BAR1BAR2IBI1IBI2NAV1NAV2
ap(443)Min0.6903.120.700.160.230.120.100.24
Max2.9411.29.980.423.940.340.541.20
Mean1.675.151.990.300.800.210.220.52
CV38.9027.1064.7025.9054.1038.4060.4025.60
aϕ(443)Min0.292.770.340.040.060.050.020.06
Max2.6210.99.190.203.410.140.180.44
Mean1.214.671.410.120.420.090.060.25
CV52.3031.9064.7028.2089.5033.3072.8031.80
anap(443)Min0.320.230.340.090.150.060.030.14
Max0.801.700.840.260.630.200.380.78
Mean0.470.490.580.180.390.120.120.27
CV25.3076.5022.7025.9027.4048.3071.2034.00
acdom(443)Min0.450.771.120.320.720.770.260.24
Max0.971.352.463.172.232.290.540.47
Mean0.831.051.711.841.291.350.260.32
CV13.216.620.233.117.338.18.512.6
ScdomMin0.0160.0040.0100.0090.0060.0110.0140.016
Max0.0180.0160.0150.0190.0130.0150.0170.020
Mean0.0170.0120.0120.0130.0090.0130.0150.018
CV3.219.610.210.923.88.84.55.7
SnapMin0.0080.0060.0100.0090.0080.0050.0080.003
Max0.0110.0090.0120.0110.0140.0220.0110.007
Mean0.0090.0080.0110.0100.0100.0130.0090.005
CV7.909.403.904.1011.6038.606.7016.70
Table 3. Mean contribution of acdom, aϕ, anap, and aw to at within the 400–700 nm of range (PAR region).
Table 3. Mean contribution of acdom, aϕ, anap, and aw to at within the 400–700 nm of range (PAR region).
Field Sites
BB1BB2BAR1BAR2IBI1IBI2NAV1NAV2
acdom26.1%19.0%40.5%72.6%49.0%72.1%22.8%18.6%
aϕ43.0%68.9%37.4%5.2%21.3%4.7%16.7%10.9%
anap18.3%7.4%14.1%8.1%16.8%6.3%31.2%29.8%
aw12.6%4.7%8.0%14.1%12.9%16.9%29.3%40.7%

Share and Cite

MDPI and ACS Style

Bernardo, N.; Alcântara, E.; Watanabe, F.; Rodrigues, T.; Carmo, A.d.; Gomes, A.C.C.; Andrade, C. Light Absorption Budget in a Reservoir Cascade System with Widely Differing Optical Properties. Water 2019, 11, 229. https://doi.org/10.3390/w11020229

AMA Style

Bernardo N, Alcântara E, Watanabe F, Rodrigues T, Carmo Ad, Gomes ACC, Andrade C. Light Absorption Budget in a Reservoir Cascade System with Widely Differing Optical Properties. Water. 2019; 11(2):229. https://doi.org/10.3390/w11020229

Chicago/Turabian Style

Bernardo, Nariane, Enner Alcântara, Fernanda Watanabe, Thanan Rodrigues, Alisson do Carmo, Ana Carolina Campos Gomes, and Caroline Andrade. 2019. "Light Absorption Budget in a Reservoir Cascade System with Widely Differing Optical Properties" Water 11, no. 2: 229. https://doi.org/10.3390/w11020229

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop