Spatiotemporal Dynamics of Suspended Sediments in the Negro River, Amazon Basin, from In Situ and Sentinel-2 Remote Sensing Data

Monitoring suspended sediments through remote sensing data in black-water rivers is a challenge. Herein, remote sensing reflectance (Rrs) from in situ measurements and Sentinel-2 Multi-Spectral Instrument (MSI) images were used to estimate the suspended sediment concentration (SSC) in the largest black-water river of the Amazon basin. The Negro River exhibits extremely low Rrs values (<0.005 sr−1 at visible and near-infrared bands) due to the elevated absorption of coloured dissolved organic matter (aCDOM at 440 nm > 7 m−1) caused by the high amount of dissolved organic carbon (DOC > 7 mg L−1) and low SSC (<5 mg L−1). Interannual variability of Rrs is primarily controlled by the input of suspended sediments from the Branco River, which is a clear water river that governs the changes in the apparent optical properties of the Negro River, even at distances that are greater than 90 km from its mouth. Better results were obtained using the Sentinel-2 MSI Red band (Band 4 at 665 nm) in order to estimate the SSC, with an R2 value greater than 0.85 and an error less than 20% in the adjusted models. The magnitudes of water reflectance in the Sentinel-2 MSI Red band were consistent with in situ Rrs measurements, indicating the large spatial variability of the lower SSC values (0 to 15 mg L−1) in a complex anabranching reach of the Negro River. The in situ and satellite data analysed in this study indicates sedimentation processes in the lower Negro River near the Amazon River. The results suggest that the radiometric characteristics of sensors, like sentinel-2 MSI, are suitable for monitoring the suspended sediment concentration in large tropical black-water rivers.


Introduction
As important components of continental landforms, large rivers play a significant role in continental erosion, and they are the main agents for transferring erosion products to the oceans [1,2]. Thus, monitoring sediment transport in large rivers contributes to the understanding of the biogeochemical cycles in order to quantify soil erosion and the formation of landscapes [3][4][5]. However, systematic and regular measurements of sediments by gauging stations in large rivers are not available and they require indirect estimates [1], especially in the tropical region. Therefore, Earth observation satellites have helped to obtain information on suspended sediment in different aquatic systems over the The total suspended sediment discharge values range from 1 to 7 Mt per year, primarily from the Guyana Shields and savannah areas, with an erosion rate of less than 10 ton. km 2 year −1 [49,53]. The Branco River, which is the main tributary of the basin in terms of discharge and suspended sediment, has an average annual flow of 3000 m 3 s −1 (in Caracaraí station) and it contributes about 50% of the suspended flow in the lower Negro River [53,54]. The mean annual suspended sediment concentration of the Negro and the Branco rivers, based on in-situ data, is ~5 mg L −1 and ~30 mg L −1 , respectively [49,[55][56][57]. The soils in the basin are covered by dense equatorial forest, poor in nutrients, and formed by the weathering of crystalline rocks [58].
The Rio Negro Basin is located in the geomorphological domains of the Guiana Shields, which consist of highly eroded Precambrian lands with little influence on the The total suspended sediment discharge values range from 1 to 7 Mt per year, primarily from the Guyana Shields and savannah areas, with an erosion rate of less than 10 ton. km 2 year −1 [49,53]. The Branco River, which is the main tributary of the basin in terms of discharge and suspended sediment, has an average annual flow of 3000 m 3 s −1 (in Caracaraí station) and it contributes about 50% of the suspended flow in the lower Negro River [53,54]. The mean annual suspended sediment concentration of the Negro and the Branco rivers, based on in-situ data, is~5 mg L −1 and~30 mg L −1 , respectively [49,[55][56][57]. The soils in the basin are covered by dense equatorial forest, poor in nutrients, and formed by the weathering of crystalline rocks [58].
The Rio Negro Basin is located in the geomorphological domains of the Guiana Shields, which consist of highly eroded Precambrian lands with little influence on the dissolved inorganic carbon flows that are carried by the rivers [23]. In this basin, old rocks predominate (granites), and there is not a significant presence of carbonates [59].
The climate in this basin is tropical without a dry season (Af), according to the Köppen classification. The average annual rainfall is 2000 mm, with values exceeding 3500 mm in the Northwest portion of the basin. The annual variability of the water level of the lower Negro River reaches 10 m, with a flood period exceeding seven months. In the study area, the high water and low water periods occur, respectively, between June and July and October and November.
Sioli [60] classified the primary types of water in the Amazonian rivers into three types: rivers with white, black, or clear water. Black-water rivers are typical in the Negro River basin, with the presence of acidic water (pH close to 4), low electrical conductivity (<20 µS), dissolved organic carbon concentration (DOC) over 8 mg L −1 , and suspended sediment concentration of less than 10 mg L −1 [25]. With an average chlorophyll-a concentration around 1 µg L −1 , the high DOC concentration increases light attenuation in the water column, limiting primary phytoplankton production in the Negro River [61,62]. Figure 2 presents a flowchart of the methodology that was used in this study. Data acquisition occurred between 2016 and 2019 at eight sample stations that were visited in different water level periods ( Figure 1b). Note that the following nomenclature was used to differentiate rivers and lakes: BCO, RN, and LAP stand for Branco River, Negro River, and Lake Apacú, respectively. Three stations are located near the confluence region of the Negro (stations RN1 and RN2) and the Branco Rivers (station BCO). The stations RN1 and RN2 are located upstream and downstream the mouth of the Branco River, respectively. Five stations are situated upstream (RN3), inside (RN4, RN5, and LAP), and downstream (RN6) of the Anavilhanas Archipelago. Inside the Anavilhanas Archipelago, the station LAP is located on Lake Apacú, and stations RN4 and RN5 are situated on channels at the right and left banks of Negro River, respectively. Each sample station corresponds to a fixed coordinate. dissolved inorganic carbon flows that are carried by the rivers [23]. In this basin, old rocks predominate (granites), and there is not a significant presence of carbonates [59]. The climate in this basin is tropical without a dry season (Af), according to the Köppen classification. The average annual rainfall is 2000 mm, with values exceeding 3500 mm in the Northwest portion of the basin. The annual variability of the water level of the lower Negro River reaches 10 metres, with a flood period exceeding seven months. In the study area, the high water and low water periods occur, respectively, between June and July and October and November.

In situ Data Collection
Sioli [60] classified the primary types of water in the Amazonian rivers into three types: rivers with white, black, or clear water. Black-water rivers are typical in the Negro River basin, with the presence of acidic water (pH close to 4), low electrical conductivity (< 20 µS), dissolved organic carbon concentration (DOC) over 8 mg L −1 , and suspended sediment concentration of less than 10 mg L −1 [25]. With an average chlorophyll-a concentration around 1 µg L −1 , the high DOC concentration increases light attenuation in the water column, limiting primary phytoplankton production in the Negro River [61,62]. Figure 2 presents a flowchart of the methodology that was used in this study. Data acquisition occurred between 2016 and 2019 at eight sample stations that were visited in different water level periods (Figure 1(b)). Note that the following nomenclature was used to differentiate rivers and lakes: BCO, RN, and LAP stand for Branco River, Negro River, and Lake Apacú, respectively. Three stations are located near the confluence region of the Negro (stations RN1 and RN2) and the Branco Rivers (station BCO). The stations RN1 and RN2 are located upstream and downstream the mouth of the Branco River, respectively. Five stations are situated upstream (RN3), inside (RN4, RN5, and LAP), and downstream (RN6) of the Anavilhanas Archipelago. Inside the Anavilhanas Archipelago, the station LAP is located on Lake Apacú, and stations RN4 and RN5 are situated on channels at the right and left banks of Negro River, respectively. Each sample station corresponds to a fixed coordinate.

Radiometry
Measurements of the remote sensing reflectance, R rs in sr −1 , were performed based on above-water radiometry acquisitions using three TriOS RAMSES spectroradiometers operating in the spectral range 350-950 nm, which is widely used in large Amazonian rivers [18,38,44]. Sensor configuration and data collection geometry were set to follow the protocol that was suggested by Mobley [63], i.e., using a viewing angle of 40 • and a relative azimuth from the Sun of 135 • . The main principle of above-water radiometry is to correct for the contribution of the surface reflection (mostly skylight) to the measured signal to properly obtain the water-leaving radiance [64]. This step is especially critical for highly absorbing water, where exiting radiances remain low and potentially yield large uncertainty on the retrieved R rs .
Hundreds of R rs were calculated with consideration of the impacts of polarisation of light [65] and the spectral dependence of the skylight refection on the water surface [66] in order to reduce potential errors and uncertainties. Hyperspectral R rs values were retrieved through equation (1), in which geometrical dependencies are ignored for the sake of brevity: where L t is the total radiance (water + surface) sensor pointing downward, L s is the sky radiance sensor pointing upward, and E s is the total above-water downward spectral irradiance. The so-called surface reflectance factor, ρ (rho), converts the skylight measurement into the contribution of the surface to the total radiance. Here, ρ was taken from radiative transfer computations through the software OSOAA [67] when considering the impacts of the aerosol load in the atmosphere and the water roughness modulated by the wind speed following the procedure that was described in Gilerson et al., [66]. Wind speed, ws, was measured during data acquisition and the aerosol optical thickness, aot, was set from the Copernicus Atmosphere Monitoring Service (CAMS) dataset [68]. Uncertainties that were attached to the retrieved R rs values were taken as the standard deviation calculated over the hundreds of individual acquisitions gathered for each sample station. Note that each acquisition sequence lasted more than 10 minutes, enabling statistically robust retrievals of R rs with lower variance [69]. This time of measurement produced a dataset that was suitable for examining the stability and variability of R rs [64].

Water Quality
A total of nine measurements at the water surface were collected in order to determine the suspended sediment concentration (SSC) in each sample station and in different hydrological periods. The SSC was determined following the protocol of HYBAM observation service [13], where the water sample is filtered with a 0.45 µm porosity filter (cellulose acetate membranes) and then weighed on an analytical balance. The uncertainties attached to the SSC measurements were taken as the standard deviation calculated over the nine measurements at the surface collected in each sample station.
Water samples (triplicates in each station) were filtered with microfiber glass membranes (47 mm diameter and 0.7 µm porosity), decontaminated for organic compounds (calcination at 450 • C for 4 h), stored in dark vials, refrigerated, and processed on the Shimadzu TOC-VCPH total carbon analyzer to determine the concentration of the dissolved organic carbon (DOC).
Absorption estimates of the coloured dissolved organic material (aCDOM) were performed using the TriOS VIPER hyperspectral spectrophotometer, which operates with a 10 cm cuvette in the range between 350 and 750 nm with a resolution of 2 nm. The water sample was filtered in a microfiber glass membrane (47 mm diameter and 0.7 µm porosity), stored in dark vials, and then refrigerated. Milli-Q water was used as a reference.

Models for Estimating SSC
A wide variety of methods are currently available for relating radiometric data with SSC, such as the use of single band [34,36,40], band ratio [18,44,70], semi-analytical methods [33,71], machine learning [72], and artificial neural networks [73]. Previous studies indicate that empirical models that estimate SSC as a function of R rs in the visible and near-infrared (VNIR) bands perform well in single-band adjusted linear regressions [18,36,42,74] and with the NIR and Red bands ratio [18,38,43].
In this study, we developed empirical models that were based on regression analysis between the SSC and the R rs . The principal characteristic of inorganic suspended sediment is its large impact on the scattering coefficient and, in turn, the increase of the water-leaving radiance due to their greater refractive index [9,31]. Thus, the regression approach was chosen because the considered SSC values range from low to moderate. For such a range, R rs increases nearly linearly before reaching asymptotic regime (or "saturation" reflectance) for greater values of SSC, typically >70 mg L −1 [75]. The Orthogonal Distance Regression (ODR) approach [76] was used in addition to the commonly used Ordinary Least Square (OLS) method to better account for the inherent uncertainties of each dataset (R rs , SSC). The primary advantage of ODR over OLS is its ability to fully take the uncertainty attached to both R rs and SSC into account [77,78].
The in situ R rs were analysed as a function of spectral shape and magnitude. Pearson correlation coefficients (r) were calculated to describe the relationship between SSC, DOC, and R rs , at different wavelengths and as a function of the hydrological period. The spatiotemporal variability of the SSC was analysed by the R rs , together with the DOC and the water level of the Negro River. This optically active constituent was not considered in this study due to the low concentration of chlorophyll-a from the Negro River (≈1 µg L −1 ).
The estimates of those uncertainties were assigned to the standard deviation calculated from the several R rs acquisitions and SSC samples of each station, respectively, as highlighted in the previous section. For each empirical model using the OLS method, the Root Mean Square Error (RMSE, Equation (2)) and Mean Absolute Percentage Error (MAPE, Equation (3)) were used to assess the performance in conjunction with the coefficient of determination (R 2 ).
where x obs and x i are the in situ and estimated SSC values, respectively; n is the number of samples.

Satellite Data and Processing
A partially cloud-free dataset of Sentinel-2A and Sentinel-2B satellites images acquired during 2017 was used to map the spatial and temporal dynamics of the SSC in the central Anavilhanas Archipelago area (stations RN4, RN5 and LAP in Figure 1). Sentinel-2 MSI data were analysed to quantify performances to remotely monitor anabranching systems that are defined by numerous narrow sections and small ponds that are only accessible through high-resolution satellite data. Note that Sentinel-2 is also well suited for such analyses because of the large number of spectral bands. Available acquisition dates of satellite Sentinel-2 images used in this study can be seen in the Supplementary material Table S1. It was not possible to assess the relationship between water quality data and satellite Sentinel-2 data due to large time lag between satellite and in situ measurements (±10 days).
In situ hyperspectral R rs was used to simulate the multispectral characteristics of the Sentinel-2 MSI sensor, while taking their specific band widths into account. The reflectance of the Sentinel-2 MSI bands was simulated based on their relative response functions. The simulation of one spectral band consists of applying the relative response functions of the multispectral sensors on the in situ hyperspectral data [69]. Table 1 presents the spectral configuration of the evaluated bands and the spectral response function of MSI sensor is available at [79]. The Sentinel-2 MSI images with Level-1C (L1C) were accessed from the Copernicus Open Access Hub (https://scihub.copernicus.eu) (accessed on 24 January 2021) that was provided by European Space Agency (ESA). Sentinel-2 satellite images were frequently contaminated by high sunglint signals (undesired reflection of sunlight on the water surface) due to the low latitude location of the Negro River. Consequently, corrections for the sunglint in addition to the atmosphere itself were mandatory in maximising the number of exploitable images and pixels within those images. For such a purpose, the algorithm GRS, standing for Glint Removal of Sentinel-2-like images [81] was used to fully process the L1C images and provide the spectral values of R rs .
The GRS algorithm was specifically developed to handle and correct for the direct sunlight reflected by the water surface and potentially reaching the sensor (i.e., sunglint signal) of Sentinel-2 missions, which is, a nadir or near-nadir viewing sensor with Shortwave Infrared (SWIR) bands. The GRS processor consists of three modules to correct for (i) gaseous absorption, (ii) diffuse light from sky and its reflection by the air-water interface, and (iii) the sunglint signal to retrieve the water-leaving signal at the water surface level.
First, the gaseous absorption (mainly CO 2 , H 2 O and O 3 ) correction is performed with the SMAC software [82]. For each pixel, the diffuse radiance component is reconstructed for the given viewing geometry (i.e., sensor and Sun viewing angles and relative azimuth) and prescribed aerosol parameter values (i.e., aerosol optical thickness and its spectral dependency) from pre-computed look-up tables (LUT). Atmospheric pressure, gas concentrations, and aerosol parameters are retrieved from bilinear interpolation within the grid of the CAMS dataset [68,83].
The sunglint correction is performed based on an estimation of the bidirectional reflectance distribution function (BRDF) of the rough water surface from the SWIR bands (i.e.,~1610 and~2200 nm). The sunglint signal that is obtained in the SWIR is then extrapolated toward the NIR and visible bands. The estimation of the sunglint radiance is based on the fact that the water body is virtually totally absorbing; the water absorption coefficient in the SWIR is several orders of magnitude greater than that in the NIR. Once corrected for the atmosphere diffuse radiance, the remaining radiance in the SWIR is interpreted as the pure surface component of the signal and then translated into BRDF. This BRDF in the SWIR is extrapolated to the other bands when considering the spectral variation of the refractive index of water and its important consequences onto the spectral variation of the sunglint signal; see Harmel et al. [81] for details. The water-leaving component at the water surface level is eventually obtained after division by the total transmittance (i.e., diffuse + total transmittances) that is calculated for the bimodal aerosol model from the LUT. Table 2  water period). After confluence with the Branco River, the Negro River at station RN2 presented an increase of SSC (average 6.18 mg L −1 ) with higher values during the highwater period. In the Anavilhanas Archipelago (stations RN4, RN5 and LAP), the mean SSC was 3.41 mg L −1 , with the lowest concentrations being observed in the high-water period. In the station RN1, the average DOC concentration of Negro River was 10.64 mg L −1 , ranging from 9.54 mg L −1 in low water period to 11.26 mg L −1 in high water period. The lowest concentration of DOC was observed in the Branco River during the low water period (4.20 mg L −1 ), with a mean value of 6.98 mg L −1 . Close to the Anavilhanas Archipelago, the mean value of DOC concentration was 9.50 mg L −1 and it ranged from 6.47 mg L −1 (station LAP) to 12.74 mg L −1 (station RN3). The dissolved organic carbon of the Negro River presented small spatiotemporal variability, with a coefficient of variation (CV) of 19%, lower than the suspended sediment (CV = 73%). Figure 3 shows the variation of light absorption coefficients (aCDOM) for eight samples stations at the Negro and Branco Rivers and a floodplain lake at the Anavilhanas Archipelago. The aCDOM at 440 nm was higher (> 7 m −1 ) in the stations that are dominated by black water rivers (rich in dissolved organic matter) and lower (< 2 m −1 ) at the station BCO (a clear water river).      Figure 4 presents the in situ R rs data of the Negro and Branco Rivers during the high and low water period. The shape and magnitude of R rs have clear differences between both rivers (Figure 4a,b). Despite the low R rs in the Negro River, spatial and seasonal differences could be observed upstream (Figure 4c Figure 4 presents the in situ Rrs data of the Negro and Branco Rivers during the high and low water period. The shape and magnitude of Rrs have clear differences between both rivers (Figure 4a Figure 4b). Despite the low Rrs in the Negro River, spatial and seasonal differences could be observed upstream (Figure 4c,d), inside (Figure 4e-g), and downstream (Figure 4h) of the Anavilhanas Archipelago.  Typical reflectance of the Negro River can be seen in Figure 4a, with R rs < 0.002 sr −1 across the in situ spectrum of station RN1 with no large variation between the high and low water periods. During the high water period, the in situ R rs of the Branco River (Figure 4b) presented higher values in the Red band (≈0.020 sr −1 between 640 and 690 nm), while, in the low water period, the R rs peak was in the green band (≈0.017 sr −1 , between 560 and 590 nm). This pattern clearly represents the seasonal variation of the SSC of this river ( Table 2). The in situ R rs of the Branco River showed absorption features at 610 and 675 nm, with a minor peak of R rs at 810 nm. The Negro River at station RN2 (Figure 4c) presented R rs with a high intensity in the Red band (between 630 and 700 nm), as a result of the contribution of suspended sediment from the Branco River. Therefore, with higher SSC, the reflectance tends to increase in the direction of longer wavelengths, especially in the Red and NIR bands, thereby modifying the spectral shape of R rs .

Relationships between Rrs and Water Quality
Upstream, inside, and downstream of the Anavilhanas Archipelago (Figure 4d-h), the R rs show a rapid increase between 500 and 650 nm, followed by a peak at 690 nm, reduction between 700 and 730 nm, and a second peak in 800 nm. These peaks in the Red band have different intensities, with higher values of R rs being observed upstream of the Anavilhanas (station RN3), in the left bank channel (station RN5), and in Lake Apacú (station LAP), sample stations with higher SSC than observed in station RN4 ( Table 2).
The in situ R rs observed in the Branco River represents an environment with moderate SSC, while the Negro River R rs represent an environment that is dominated by light absorption due to dissolved organic matter, with reduced reflectance in the visible (between 400 and 760 nm) and a slight increase in NIR (between 760 and 900 nm). Therefore, between the stations BCO and RN6, a reduction of up to 90% in R rs at 665 nm was observed, which is related to the reduction of SSC in the reach between the mouth of the Rio Branco and upstream from the mouth of the Rio Negro.
We observed that SSC was significantly correlated with DOC (r = −0.68, p < 0.05) ( Figure 5). The magnitude of the correlation coefficient was high during the high water level period (r = −0.83, p < 0.05). The SSC was not significantly correlated with DOC (r = −0.52, p > 0.05) in the low water period.
Anavilhanas (station RN3), in the left bank channel (station RN5), and in Lake Apacú (station LAP), sample stations with higher SSC than observed in station RN4 ( Table 2).
The in situ Rrs observed in the Branco River represents an environment with moderate SSC, while the Negro River Rrs represent an environment that is dominated by light absorption due to dissolved organic matter, with reduced reflectance in the visible (between 400 and 760 nm) and a slight increase in NIR (between 760 and 900 nm). Therefore, between the stations BCO and RN6, a reduction of up to 90% in Rrs at 665 nm was observed, which is related to the reduction of SSC in the reach between the mouth of the Rio Branco and upstream from the mouth of the Rio Negro.
We observed that SSC was significantly correlated with DOC (r = -0.68, p < 0.05) ( Figure 5). The magnitude of the correlation coefficient was high during the high water level period (r = −0.83, p < 0.05). The SSC was not significantly correlated with DOC (r = −0.52, p > 0.05) in the low water period.

Multispectral MSI simulated bands correlation with SSC and DOC
The Pearson coefficients correlation (r) between SSC and Sentinel-2 MSI simulated bands indicated a positive relationship at different hydrological periods (Figure 6b,c). During the high water period, the correlation was higher, ranging from 0.77 to 0.98 (p < Figure 5. Relationship between DOC and SSC including low water period (red) and high water (blue) period. The shaded areas indicate the uncertainty.

Multispectral MSI Simulated Bands Correlation with SSC and DOC
The Pearson coefficients correlation (r) between SSC and Sentinel-2 MSI simulated bands indicated a positive relationship at different hydrological periods (Figure 6b,c). During the high water period, the correlation was higher, ranging from 0.77 to 0.98 (p < 0.05, Figure 6b). For all of the datasets (low water and high water periods), the r ranged between 0.67 (at 865 nm) and 0.92 (at 665 nm) (Figure 6a).
Conversely, the correlations between DOC and Sentinel-2 MSI simulated bands indicated a negative relationship. During the high water period, the correlation coefficients ranged from −0.75 to −0.87 (p < 0.05, Figure 6b) and, in the low water period, the correlation ranged from −0.60 to −0.64 (p < 0.05, Figure 6c). For all the DOC datasets, the correlation ranged from −0.62 (at 560 nm) to −0.66 (at 705 nm) ( Figure 6).
In general, the visible and red-edge bands showed high correlations with SSC while the NIR (Band 8a) and the NIR/Red ratio (Band 8a/Band 4) showed lower correlations. This evaluation indicates that the visible bands have better performance for retrieval SSC of Sentinel-2 MSI data on the Negro River.

SSC Estimations Based in Sentinel-2 MSI Simulated Bands
Several VNIR Sentinel-2 MSI simulated bands were tested based on the observed correlations ( Figure 6), and Figure 7 shows the results of the SSC as a function of R rs in the lower Negro River basin. Using the OLS method, the models with the best fit were obtained while using simulated band 3 (R rs at 560 nm) and band 4 (R rs at 665 nm) of Sentinel-2 MSI (R 2 > 0.80, RMSE < 1.50 mg L −1 and percentage error < 17%) ( Table 3). Good results were also obtained with the band 5 (R 2 > 0.80, RMSE ≈ 1.50 mg L −1 and percentage error < 17%). The lowest performance was obtained using the ratio of NIR and Red bands (R 2 < 0.15, RMSE > 5 mg L −1 and percentage error > 60%). coefficients ranged from −0.75 to −0.87 (p < 0.05, Figure 6b) and, in the low water period, the correlation ranged from −0.60 to −0.64 (p < 0.05, Figure 6c). For all the DOC datasets, the correlation ranged from −0.62 (at 560 nm) to −0.66 (at 705 nm) ( Figure 6).
In general, the visible and red-edge bands showed high correlations with SSC while the NIR (Band 8a) and the NIR/Red ratio (Band 8a/Band 4) showed lower correlations. This evaluation indicates that the visible bands have better performance for retrieval SSC of Sentinel-2 MSI data on the Negro River.

SSC estimations based in Sentinel-2 MSI simulated bands
Several VNIR Sentinel-2 MSI simulated bands were tested based on the observed correlations ( Figure 6), and Figure 7 shows the results of the SSC as a function of Rrs in the lower Negro River basin. Using the OLS method, the models with the best fit were obtained while using simulated band 3 (Rrs at 560 nm) and band 4 (Rrs at 665 nm) of Sentinel-2 MSI (R 2 > 0.80, RMSE < 1.50 mg L −1 and percentage error < 17%) ( Table 3). Good   Figure 8 presents the SSC observed versus the SSC derived from Sentinel-2 MSI simulated band 4. It was observed that the SSC derived from simulated band 4 indicated a goodness of fit (R 2 = 0.86 for OLS and ODR methods) and, for SSC values lower than 6 mg L −1 , the OLS method overestimated the observed values. Figure 8 showed that the main advantage of the ODR method was presenting a lower tendency (less-biased slope values) as compared with those that were obtained using the OLS method (decrease 36%, from 1.41 to 0.90). Therefore, it is emphasized that the ODR method (which considers the associated uncertainties) is more suitable than the OLS method for estimating SSC. results were also obtained with the band 5 (R 2 > 0.80, RMSE ≈ 1.50 mg L −1 and percentage error < 17%). The lowest performance was obtained using the ratio of NIR and Red bands (R 2 < 0.15, RMSE > 5 mg L −1 and percentage error > 60%).   Figure 8 presents the SSC observed versus the SSC derived from Sentinel-2 MSI simulated band 4. It was observed that the SSC derived from simulated band 4 indicated a goodness of fit (R 2 = 0.86 for OLS and ODR methods) and, for SSC values lower than 6 mg L −1 , the OLS method overestimated the observed values. Figure 8 showed that the main advantage of the ODR method was presenting a lower tendency (less-biased slope values) as compared with those that were obtained using the OLS method (decrease 36%, from 1.41 to 0.90). Therefore, it is emphasized that the ODR method (which considers the associated uncertainties) is more suitable than the OLS method for estimating SSC.

SSC Spatiotemporal Variability from Satellite Imagery
Inside the Anavilhanas, it was possible to detect a relationship between the Rrs at the Red band (Sentinel-2 MSI band 4) and the water level of the Negro River (Figure 9(a)). The Rrs at the Red band of station RN4 were progressively lower than that observed in station RN5. This difference may be related to the higher SSC at the station RN5, a region that is influenced by the input of suspended sediments from the Branco River that slightly modifies the Negro River AOPs during the high water period.
At the station RN5, the maximum SSC occurs near the flood peak of the Negro River, while, at station RN4, it was not possible to detect a clear behaviour in the SSC and water level (Figure 9b). Except in the data that were collected during November 2017 (see Figure 9b), the SSC value was greater in the channel on the left bank (station RN5) than the channel on the right bank (station RN4). Conversely, the DOC was always higher in station RN4 as compared to that in station RN5.

SSC Spatiotemporal Variability from Satellite Imagery
Inside the Anavilhanas, it was possible to detect a relationship between the R rs at the Red band (Sentinel-2 MSI band 4) and the water level of the Negro River (Figure 9a). The R rs at the Red band of station RN4 were progressively lower than that observed in station RN5. This difference may be related to the higher SSC at the station RN5, a region that is influenced by the input of suspended sediments from the Branco River that slightly modifies the Negro River AOPs during the high water period. The relationship obtained with ODR method using band 4 in Section 3.2.2 was applied to Sentinel-2 MSI images over year 2017. This satellite dataset is used to compare different water levels during the annual cycle. Figure 10 presents the estimated SSC distributions based on Sentinel-2 MSI band 4 images during the rising, high, falling, and low water periods of 2017. In this reach of the Anavilhanas Archipelago, the Negro River has two primary channels and several lakes with different shapes and sizes. The Sentinel-2 MSI band 4 images captured the SSC spatiotemporal variability well. The high detail, the predominance of suspended sediment in the floodplain lakes and channels close to the left bank of the Anavilhanas Archipelago are very noticeable due to the influence of the suspended sediment from Branco River. At the station RN5, the maximum SSC occurs near the flood peak of the Negro River, while, at station RN4, it was not possible to detect a clear behaviour in the SSC and water level (Figure 9b). Except in the data that were collected during November 2017 (see Figure 9b), the SSC value was greater in the channel on the left bank (station RN5) than the channel on the right bank (station RN4). Conversely, the DOC was always higher in station RN4 as compared to that in station RN5.
The relationship obtained with ODR method using band 4 in Section 3.2.2 was applied to Sentinel-2 MSI images over year 2017. This satellite dataset is used to compare different water levels during the annual cycle. Figure 10 presents the estimated SSC distributions based on Sentinel-2 MSI band 4 images during the rising, high, falling, and low water periods of 2017. In this reach of the Anavilhanas Archipelago, the Negro River has two primary channels and several lakes with different shapes and sizes. The Sentinel-2 MSI band 4 images captured the SSC spatiotemporal variability well. The high detail, the predominance of suspended sediment in the floodplain lakes and channels close to the left bank of the Anavilhanas Archipelago are very noticeable due to the influence of the suspended sediment from Branco River. Figure 11a shows an overview of the SSC over the Anavilhanas Archipelago in the high and low water periods of the Rio Negro in 2017, highlighting the predominance of higher values close to the left bank of the archipelago. Figure 11b presents the cross profile of SSC between the left and right banks of the Negro River. The difference between the satellite estimates of SSC in the left and right banks were higher than 50% in the high and low water periods.   Figure 11b presents the cross profile of SSC between the left and right banks of the Negro River. The difference between the satellite estimates of SSC in the left and right banks were higher than 50% in the high and low water periods.

Suspended Sediment Transport in an Anabranching River
The lower Negro River is an old and complex multichannel system due to the presence of the Anavilhanas Archipelago [84], with a very low suspended sediment concentration and under the strong influence of backwater effects of the Amazon River [85]. This hydro-sedimentary behaviour allows for the storage of water and sediments in this archipelago, which are primarily sourced from the Branco River [49]. In the results presented (Table 2), the Negro River downstream from the mouth of the Branco River, (station RN2) in comparison to the upstream sample station (RN1), showed an increase in the SSC of 39%. On the other hand, the mean SSC difference between the sample stations RN2 and RN6 (250 km distance) shows a reduction of 75%, which indicated suspended sediment retention in the Anavilhanas Archipelago. Despite the high volume of water,

Suspended Sediment Transport in an Anabranching River
The lower Negro River is an old and complex multichannel system due to the presence of the Anavilhanas Archipelago [84], with a very low suspended sediment concentration and under the strong influence of backwater effects of the Amazon River [85]. This hydrosedimentary behaviour allows for the storage of water and sediments in this archipelago, which are primarily sourced from the Branco River [49]. In the results presented (Table 2), the Negro River downstream from the mouth of the Branco River, (station RN2) in comparison to the upstream sample station (RN1), showed an increase in the SSC of 39%. On the other hand, the mean SSC difference between the sample stations RN2 and RN6 (250 km distance) shows a reduction of 75%, which indicated suspended sediment retention in the Anavilhanas Archipelago. Despite the high volume of water, this sedimentation process is related to the very low velocity of the Negro River (≈0.50 m s −1 ), the backwater effect from the Amazon River, and the low erosion process in the Anavilhanas Archipelago [49,52].
In the channels of the Negro River without influence of suspended sediment from the Branco River (sample stations RN1, RN4 and RN6 in Figure 10), the low SSC values registered during the high water period could be related to the dilution process that was caused by the larger volume of water available. The hydro-sedimentary dynamics of the Negro River are similar to the process observed at the confluence between the Negro and Amazon rivers [38] and in the Amazon floodplain [86,87].
In the Amazon River and its main tributaries, the concentration of dissolved organic carbon decreases as the concentration of suspended sediment increases, as highlighted by Moreira-Turcq et al. [25]. In the lower Negro River, we observed that SSC was found to correlate inversely with DOC, primarily in the high water period ( Figure 5). The water of the Negro River is dominated by dissolved organic matter [23,25] that causes intense light absorption due to coloured dissolved organic matter (CDOM). The Negro River exhibits the higher values of aCDOM when compared to those of Amazon Rivers and its floodplains lakes [46,88] or the reservoirs in the Brazilian semiarid region [89], while demonstrating the same magnitude as those of the Estonian lakes dominated by CDOM [90], as observed in Figure 3.
Even in environments with strong absorption due to CDOM, the results that are presented here (Figure 7) indicate that models with a single spectral band are best suited for water with low suspended sediment concentrations. Previous studies also reported good results for SSC estimation of inland waters using single bands [33,34,36,41,42]. As expected, a higher R rs in the Negro and Branco Rivers is related to an increase of light scattering by mineral particles [9,18].
In the central area of Anavilhanas (see stations RN4 and RN5 in Figure 1), the water discharge flows through two large channels near to the right and left banks of the archipelago. However, more SSC in the channel on the left bank of the Negro River was observed from in situ (sample station RN5 in Table 2) and Sentinel-2 data (Figure 11a). For example, in the high water period, the difference in suspended sediment discharge between the stations RN4 and RN5 can be over ten times [49]. Inside Anavilhanas Archipelago, the highest amplitude of SSC was observed in the Lake Apacú (station LAP in Table 2), with the minimum and maximum values being obtained during the high and low water periods, respectively. During the low water period, part of the water that was stored in the Anavilhanas lakes throughout the high water period flows to the Negro River channels, and this dynamic may influence the process that increases the concentration of the SSC in some channels and lakes of Anavilhanas Archipelago, as observed in Figures 10 and 11.
Despite the high aCDOM values, the Sentinel-2 MSI data could confirm the spatial heterogeneity of the SSC estimated with in situ R rs data. During the high and low water periods, the preferred route of suspended sediment inside Anavilhanas can be observed in great detail (Figure 11a). The Sentinel-2 MSI images that were obtained in 2017 allowed for us to see that, even 90 km from its mouth, the suspended sediment from the Branco River influences the R rs in the Negro River at the Anavilhanas Archipelago, a dynamic that produces differences in the composition of the local fauna and flora [56,91,92]. The SSC cross profile that is retrieved by Sentinel-2 data (see Figure 11b) details this biogeochemical gradient, showing a tendency of reduced suspended sediment concentration between the left bank (station RN5), lakes (station LAP), and right bank (station RN4) of the Anavilhanas Archipelago.

Use of Passive Satellite Remote Sensing for SSC Monitoring in the Negro River
Multispectral sensors, like Sentinel-2 MSI, were not developed for remote sensing of water; nevertheless, the results here indicate that these spectral configurations are suitable for monitoring SSC in black-water rivers with complex channel pattern, like the Negro River. In this study, the largest percentage error was obtained with the model that uses the NIR and NIR/Red bands ratio, which may be a result of the low NIR response over the Negro River (Table 3). Despite the good performance of the green band (at 560 nm), this band onboard in orbital sensors can be strongly impacted by atmospheric constituents in shorter wavelengths and by sunglint [48,93].
We showed that full consideration of uncertainties attached to the in situ dataset changes the parameterisation of the satellite bands in terms of suspended sediment concentration. A comparison between the OLS (no uncertainty consideration) and ODR methods demonstrated that an estimation of SSC with Sentinel-2 MSI simulated band 4 is less underestimated with the latter (see Figure 8). It is recommended that efforts to fully characterise bias and uncertainties attached to the in situ data and, in particular, for the radiometry measurements are pursued [64][65][66]. The good adjustment that is obtained between SSC and R rs in the Red band (at 665 nm) is influenced by the presence of suspended sediment from the Branco River.
The lower Negro River is a complex anabranching system with channel width varying from 230 m to 2800 m. Despite this environmental configuration, the Sentinel-2 MSI images thoroughly captured the temporal and spatial variability of the SSC in the Anavilhanas channels and lakes, highlighting, in detail, the predominance of suspended sediments in the lakes and channels near the left bank channel of Anavilhanas.
Although the combination of orbital sensors enables higher temporal resolution [17], the probability of cloud-free imaging in the study area is low throughout a hydrological year [94]. However, multispectral sensors onboard unmanned aerial vehicles can partially overcome these limitations, since they operate at low altitude, and with greater flexibility in revisit time and image acquisition geometry. Conversely, the use of commonly used watercolour sensors with higher temporal resolutions, such MODIS, OLCI, and NPP/VIIRS, allows for us to retrieve a time series of the SSC in the Negro River reaches with channel width >2000 m by using a similar approach that is applied for other large rivers of the Amazon basin [43][44][45][46]. However, the satellite sensors with higher spatial resolutions generally have a greater capacity to detect higher SSC values [95].
The Negro River basin has the lowest coverage in the sediment measurement network that is managed by the Brazilian National Water Agency. Without sedimentometric stations around the Anavilhanas Archipelago, it is difficult to assess temporal variations of the SSC and the relationships between extreme climatic events (El Niño/La Niña) or the temporal variations that will occur due to the Bem Querer dam [96] in the Branco River and subsequently reduce the amount of sediment that arrives in the Anavilhanas Archipelago. The suspended sediment concentration dataset collected in this study before and after they arrive at the archipelago indicate that the Anavilhanas Archipelago retains suspended sediments and acts as sink of fine particles in the lower Negro River (Table 2). Thus, the use of passive satellite data could be an important tool for monitoring suspended sediment transport in areas of great ecological and socioeconomic relevance, such as the Anavilhanas Archipelago.

Conclusions
This is the first study that has utilised Sentinel-2 MSI calibrated data for assessing the suspended sediment concentration (SSC) in a large black-water Amazonian river. The waters in the Negro River basin exhibit an absorption budget that is dominated by coloured dissolved organic material. We demonstrated that the so-called black-water rivers, such as the Negro River, exhibit very low suspended sediment concentrations (≈5 mg L −1 ) and a high coloured dissolved organic material absorption (aCDOM > 7 m −1 at 440 nm). We show that remote sensing reflectance at 665 nm (Red band) can be used in order to monitor suspended sediment concentration in the Negro River with high accuracy (R 2 > 0.86 and error less than 30% in adjusted models). Thus, passive sensors that are used in satellites operating in the 650-690 nm range can be used to remotely monitor suspended sediment concentrations in this large river. The orthogonal distance regression method demonstrated that estimation of SSC using simulated bands of Sentinel-2 MSI results in less-biased slope values when compared to those obtained using ordinary least square regression.
Our study also highlighted that the suspended sediment from the Branco River causes changes in the apparent optical properties on the low course of Negro River. Indeed, the remote sensing reflectance pattern that was detected with the Sentinel-2 MSI images in the channels and lakes of the Anavilhanas Archipelago was not only similar to the in situ data, but it also showed higher suspended sediment concentration in channels and lakes that were influenced by the Branco River, which is the primary clear water river in the Negro River basin. The suspended sediment transport characterised in this study indicates sedimentation processes in an anabranching reach of the Negro river upstream of the Amazon River.
Therefore, the in situ reflectance dataset from this study contributed to the knowledge of the apparent optical properties in Amazonian rivers. However, further research must be carried out to better assess the inherent optical properties that are based on multispectral and hyperspectral approaches to exploit the full information content of the measured spectra. This will be particularly important for the monitoring of the Negro River basin, including impacts of extreme climate events and potential changes in sediment transport due to deforestation, mining, or the construction of a hydropower dam on this basin. Furthermore, the authors thank the anonymous reviewers who provided valuable suggestions to improve the manuscript.