Dynamics of the Estuarine Turbidity Maximum Zone from Landsat-8 Data: The Case of the Maroni River Estuary, French Guiana

: The estuarine turbidity maximum (ETM) zone occurs in river estuaries due to the e ﬀ ects of tidal dynamics, density-driven residual circulation and deposition / erosion of ﬁne sediments. Even though tropical river estuaries contribute proportionally more to the sediment supply of coastal areas, the ETM in them has been hardly studied. In this study, surface suspended particulate matter (SPM) determined from OLI (Operational Land Imager)-Landsat 8images was used to gain a better understanding of the spatio-temporal dynamics of the ETM of the tropical Maroni estuary (located on the Guianas coast, South America). A method to estimate the remotely-sensed ETM location and its spatiotemporal evolution between 2013 and 2019 was developed. Each ETM was deﬁned from an envelope of normalized SPM values > 0.6 calculated from images of the estuary. The results show the inﬂuence of the well-marked seasonal river discharge and of tides, especially during the dry season. The ETM is located in the middle estuary during low river-ﬂow conditions, whereas it shifts towards the mouth during high river ﬂow. Neap–spring tidal cycles result in a push of the ETM closer to the mouth under spring-tide conditions or even outside the mouth during the rainy season. An increase in SPM, especially since 2017, coincident with an extension of the ETM, is shown to reﬂect the periodic inﬂuence of mud banks originating from the mouth of the Amazon and migrating along the coast towards the Orinoco (Venezuela). These results demonstrate the advantages of ocean color data in an exploratory study of the spatio-temporal dynamics of the ETM of a tropical estuary, such as that of the Maroni.


Introduction
The trapping of suspended particles in river estuaries can often lead to the formation of highly turbid concentrations, sometimes culminating in the generation of an estuarine turbidity maximum (ETM). The ETM commonly occurs near the upstream limits of the salinity intrusion and can be the result (c) View of the estuary depicting locations where river discharge (blue dot), wind data (green dot), tide levels (point: Les Hattes, red dot) and wave data (purple dot) were collected. S1 and S2 correspond to the satelliteimage areas used to estimate SPM evolution. Location of Kilometric Points (KP) 10 and 20 (yellow dots), which represent distances, in kilometers, along the estuary starting from KP 0. (d) Bathymetry of the Maroni estuary from two surveys conducted in 2016 and 2017 [36]. Bathymetry in (a) is from GEBCO bathymetry data and in (c) from NOAA Bathymetric Data Viewer.
The salinity intrusion length in the Maroni varies from the dry to the wet season, reaching ∼12 and ∼5 km upstream from the mouth (corresponding to Kilometric Point (KP) 0; Figure 1d), respectively, with strong stratification during the rainy season marked by a salty mass of water during high tides [15]. This study also showed a high seasonality in the estuary's turbidity and (c) View of the estuary depicting locations where river discharge (blue dot), wind data (green dot), tide levels (point: Les Hattes, red dot) and wave data (purple dot) were collected. S1 and S2 correspond to the satellite-image areas used to estimate SPM evolution. Location of Kilometric Points (KP) 10 and 20 (yellow dots), which represent distances, in kilometers, along the estuary starting from KP 0. (d) Bathymetry of the Maroni estuary from two surveys conducted in 2016 and 2017 [36]. Bathymetry in (a) is from GEBCO bathymetry data and in (c) from NOAA Bathymetric Data Viewer.
Although the Maroni and several other Guiana estuaries are large estuaries in their own right, the periodic influence of Amazon mud is expected to induce patterns in ETM distribution and dynamics that can be quite challenging to unravel as field approaches systematically entail considerable logistical difficulties [31]. Earth observation methods can overcome to a certain extent the problems of spatial and temporal relevance of in-situ estuarine monitoring, as they can provide a synoptic perspective. Remote sensing data from satellite imagery have, indeed, enabled advances in understanding of the dynamics of the large mud-bank and inter-bank systems [29,37,38]. Recently-launched high-resolution sensors, such as Landsat 8 (16 days of temporal resolution and 15-60 m of spatial resolution) and Sentinel 2 (5 days of temporal resolution and 10-60 m of spatial resolution), present a spatial scale that can improve significantly the study of the small-scale spatial variability of the high inshore mud-concentrations of the Guianas [29]. The estimation of suspended particulate matter (SPM) concentration from ocean Remote Sens. 2020, 12, 2173 4 of 23 color remote sensing can, therefore, potentially provide relevant information on mud dynamics in estuaries such as the Maroni river.
The objective of this work was to assess the potential of such high-resolution sensor information in contributing to the understanding of the spatio-temporal dynamics of the Maroni estuary under a phase of significant mud-bank influence that has prevailed since 2011 [36]. Due to the complex dynamics that characterize the ETM zone and the challenges of identifying SPM depth concentration gradients, and given the opportunities for the observation of the surface of the water column offered by satellite data, we defined in this study a number of ETM situations from an envelope of normalized SPM values. From this, we propose an exploratory approach that can lay the foundations for more detailed remote sensing-based studies of the ETM. In the frame of this study, SPM variability in the Maroni waters was analyzed from Landsat 8 scenes for the period between April 2013 and June 2019 and confronted with dynamic drivers deemed as important in influencing the ETM. The spatio-temporal evolution of the ETM was evaluated and described based on regional climate seasonality and changes in neap-spring tidal cycles, and the identified patterns are discussed below.

Maroni River and Estuary
The Maroni River is located in French Guiana, South America (Figure 1), and is subject to a humid tropical-equatorial climate. The region is influenced by seasonal variations caused by the north-south movement of the Intertropical Convergence Zone (ITZC) with a well-marked rainy season from the end of December to July and a dry season from August to December. The rainy season is usually interrupted by a one-month period with dry conditions around March.
The Maroni River is 612 km long and lies astride French Guiana and Suriname (Figure 1a,b), and its mouth lies approximately 700 km northwest of the mouths of the Amazon. Its catchment covers an area of 66,184 km 2 , cut in crystalline rocks of the Guiana Shield, and it has a mean water discharge of ∼774 m 3 /s during the dry season and 2400 m 3 /s during the rainy season. Rousseau et al. [34] showed that the Maroni's seasonality index (SI), which is the ratio between the highest monthly discharge divided by the lowest one, is much higher than that of the Amazon, thus highlighting the high seasonality of the freshwater input. This is probably related to the lower dampening effect on water discharge resulting from the much smaller catchment of the Maroni, compared to the Amazon (catchment size of 3.6 million km 2 and mean discharge of 200,000 m 3 ). Discharge, however, is not erratic in this equatorial setting and increases gradually as the rainy season sets in (see Section 2.3). The Maroni presents a large funnel-shaped shallow estuary mouth rich in sand deposits, indicating the importance of fluvial bedload supply in the system [36].
Tides in the Maroni estuary are semidiurnal and low mesotidal with a range from up to 2.5 m during spring tides [22] to up to 1.1 m during neap tides. The estuary shows a converging plan shape, with tidal range amplification at about 8 km from the mouth [15]. Ross et al. [16] identified an ebb channel at the west bank of the estuary (Suriname) and a flood channel at the east bank (French Guiana) (Figure 1c). A mud bank started impinging on the mouth of the Maroni in 2011 (Figure 1d) [36], and it is still migrating across the estuary.
Onshore-directed northeast trade winds characterize the rainy season, whereas offshore southwest trade winds are dominant when the ITZC is at its northern position during the dry season [27]. The strongest winds in the year occur from January to April [39]. The interaction of the North Brazil Current (NBC), the seasonally variable trade winds and surface gravity waves result in the advection of Amazon mud in a low-salinity plume along the coast to the northwest [40,41]. This migration affects the coastal waters off French Guiana, Suriname, Guyana and Venezuela making this the longest muddy coast in the world.
The salinity intrusion length in the Maroni varies from the dry to the wet season, reaching ∼12 and ∼5 km upstream from the mouth (corresponding to Kilometric Point (KP) 0; Figure 1d), respectively, with strong stratification during the rainy season marked by a salty mass of water during high tides [15]. This study also showed a high seasonality in the estuary's turbidity and identified an ETM during the rainy season at the river mouth [15]. Two large tidal creek networks debouch on the river's east bank near the mouth, and a series of shoals and vegetated islands have accumulated between KP 15 and KP 26.5 (Figure 1c). The total suspended solid concentration in French Guiana rivers is low [42], and the Maroni has been identified as among the world's rivers with the lowest concentrations [43], no doubt a reflection of the highly crystalline basement rocks and still largely forested nature of the catchments. However, Gallay et al. [44] identified a notable 239% increase in the Maroni River's SPM after 2009 due to an upsurge in anthropogenic activities, notably gold mining and deforestation.
The Maroni River banks host the twin towns of St. Laurent du Maroni (French Guiana) and Albina (Suriname) about 20 km upstream from the mouth. These are among the fastest growing towns in the Guianas, notably St. Laurent du Maroni, where port facilities capable of handling large ocean-going vessels are being planned, providing a further justification for a better understanding of the mud dynamics at the mouth of the river.

Landsat 8 OLI Data Set
Landsat 8 images were used in this study. The Operational Land Imager (OLI) sensor has eight bands at 30 m spatial resolution and a panchromatic one at 15 m spatial resolution. The sensors present four visible bands (1-4), a near infrared band (NIR) and two bands (6 and 7) in the shortwave infrared spectral region (SWIR).
OLI level 1 data (L1T) corresponding to the period between April 2013 and June 2019 were downloaded from the Earth Explorer database (https://earthexplorer.usgs.gov/). In total, 66 Landsat 8-OLI images (Table 1) were used after images highly affected by clouds were discarded. Quasi-permanent glint contamination is present in the study area as a consequence of the characteristics of the OLI sensor, such as the solar viewing path geometry, and the low latitude of the Maroni River and estuary (5 • -6 • N). This results in the loss of contaminated pixels over the river and estuary. To recover these contaminated pixels, the images were corrected for this effect using the method proposed by Abascal-Zorrilla et al. [45].The method uses SWIR information to perform an automatic correction of this effect based on the assumption that the water-leaving radiance is negligible on this band. For application of a blind correction, several masks have to be defined, such as the land mask, the clear waters mask and the sun glint mask. Once the regions of interest are defined, the SWIR-based empirical approach was deployed following a two-step procedure: the definition of the spectral slope parameter, which allows relating the sun glint observed in the SWIR signal to the one of the NIR and visible bands to be corrected, and the estimation of the offset value corresponding to the background SWIR 2 of the areas not affected by sun glint.
ACOLITE software (acolite_win 20170718.0), developed at the Royal Belgian Institute of Natural Science (RBINS), was used to process L1T sun glint-corrected images. No external measurements are needed to perform the atmospheric correction using this software. The OLI dataset was processed with the ACOLITE SWIR processing using a fixed epsilon computed on the data prior to sun glint correction, as explained by Abascal-Zorrilla et al. [45]. The ACOLITE atmospheric correction chain used in the framework of this study was validated for French Guiana [45] and R rs data were obtained for the study area.

Post-Treatment of R rs Data
Clouds and Cloud-Shadow Detection Cloud-shadow detection is one of the limitations of high spatial-resolution sensors. In this study, water pixels (considering the detection of cloud shadows and thin clouds) of the ACOLITE sun glint-corrected images were identified using an algorithm of water-pixel extraction (WiPE) [46]. This algorithm was developed for Landsat 8-OLI and Sentinel 2-MSI and consists of a two-step procedure that considers the statistical analysis of the Rayleigh-corrected top of the atmosphere (TOA) reflectance (ρ rc (λ)). The first step uses (ρ rc (λ)) and the spectral shape analysis of each object is performed allowing water pixels to be identified and separated from cloud, vegetation, barren land and construction pixels. The Hue-Saturation-Value space is estimated during the second step, which allows a better detection of cloud shadow over waters and considering the whole image. The two-step algorithm allows the identification of thin clouds and showed a better performance over turbid waters than other algorithms [46].

Estimation of SPM
The estimation of SPM was conducted using a semi-analytical algorithm developed by Han et al. [47]. This algorithm was chosen because it covers four orders of magnitude, from clear to very turbid waters and because it was developed for several sensors, including Landsat 8, and for several water masses, including French Guiana.
SPM estimation was defined using two distinct algorithms related to turbidity levels of the water masses. A smooth procedure was applied to avoid any regularity in the turbidity levels between medium and high turbidity situations. The reader interested in the algorithm is referred to the work of Han et al. [47].

ETM Detection
Based on the different SPM attributes of the images and given that the ETM is affected by different dynamic drivers, its clear identification at the water column surface becomes, in some cases, a difficult task. To overcome this hurdle, a simple methodology was applied to the SPM values (in mg/L, estimated using the algorithm evoked in Section 2.2.3 to identify the presence of the ETM and its dynamics. 3-Definition of the ETM zone: The limits of the ETM were defined from an envelope of SPM-normalized values higher than 0.6. This > 0.6 threshold value was empirically defined following extensive work using remote sensing data on the muddy coast of the Guianas [29], and it was averaged for the different years and seasonal conditions. 4-Length and core of the ETM zone: The length was estimated to gain insight on ETM spatio-temporal variations. We use the terminology of Uncles et al. [48] to identify a "tail" and a "nose" (which represent, respectively, the downstream and upstream separation of the ETM from areas of lower turbidity). The core of the ETM was estimated from the barycenter using the geometric properties tool in QGIS software.

Data Processing
To identify the relationship between SPM values, on which ETM identification is based, and the influence of different dynamic drivers, a mean SPM value was estimated using the selected scenes without applying the previous normalization step of the methodology (see Sections 3 and 4).
Image classification as a function of the tide (neap, spring or mean) was performed as follows based on the tidal range: 1-Estimate the peak values of the tide series. 2-Estimate the mean and standard deviation values (peak std) of the identified peaks.

Water Levels, River Discharge, Wind and Wave Data
Maroni water discharge data are only available as daily means and were obtained from the database of the service Fleuves, Littoral, Aménagement et Gestion (FLAG) of the DEAL (Direction de l'environnement, de l'aménagement et du logement) (http://hydro.eaufrance.fr/). Data were collected at Apatou (~80 km upstream from KP 0; Figure 1b) between January 2013 and June 2019.
Water levels were retrieved from the tidal gauge record of Les Hattes (5 • 45 00"N, 53 • 57 00"W Figure 1b). Due to the short length of the record, which does not cover the entire study period, water level data were also retrieved from the tidal gauge of Île Royale (5 • 17 06" N, 52 • 34 48" W Figure 1a). The data are managed by the Service Hydrographique et Océanographique de la Marine (SHOM) and Direction de la Mer-Guyane (http://data.shom.fr/). Both series were compared to extract the lag difference between Île Royale and Les Hattes, in order to use the Île Royale data with the lag time corrected.
Wind velocity data (averaged every 10 min) were provided by MeteoFrance at Mana airfield (5 • 39 48" N, 53 • 45 42" W; Figure 1b) for the period 2015-2019. As our study period starts in April 2013 and the station data start in 2015, wind velocity data were also retrieved from the WAVEWATCH III (WWIII) model [49][50][51] which is a third-generation wave model developed by NOAA's National Center for Environmental Prediction (NCEP). The data concern the u and v wind components (scaled to a reference height of 10 m). WWIII monthly averaged data ( Figure 2c) turned out to be accurate when compared with in-situ data (R 2 = 0.9) ( primary wave directions (degrees true, defined as 0° coming from north and 90° coming from east). These variables were monthly-averaged for the study area ( Figure 2e). Waves in the Maroni estuary area arrive from an east to northeast direction (Figure 2d,f) in response to the predominant trade winds. High waves prevail between October and May, and smaller waves from June to September (Figure 2e). Significant wave heights have tended to exhibit a decreasing trend since 2015, with maximum values dropping from 2.3 to 1.6 m (Figure 2e).

Results and Discussion
In this section, an analysis and discussion of the evolution of the spatial patterns of the ETM, identified using the SPM-envelope methodology defined in Section 2.2.4, is performed as a function of seasonal and interannual conditions. Besides the influence of the tide and the neap-spring cycle in both the dry and the rainy seasons, the way these two seasons affect the ETM is also discussed. Other points of discussion are the influence of river discharge on the ETM and of the presence of the mud Wave data were also retrieved from WWIII for an offshore point where water depth is about 25 m (Figure 1b, wave data), well beyond the zone of coastal mud-bank influence on waves. The model provided significant wind and swell wave heights (m), the primary mean wave periods (s) and the primary wave directions (degrees true, defined as 0 • coming from north and 90 • coming from east). These variables were monthly-averaged for the study area (

Results and Discussion
In this section, an analysis and discussion of the evolution of the spatial patterns of the ETM, identified using the SPM-envelope methodology defined in Section 2.2.4, is performed as a function of seasonal and interannual conditions. Besides the influence of the tide and the neap-spring cycle in both the dry and the rainy seasons, the way these two seasons affect the ETM is also discussed. Other points of discussion are the influence of river discharge on the ETM and of the presence of the mud bank in the study area since 2011 [36]. Due to the complexity of the dynamics in the study area, the differences in SPM values for the different images and the way they influence detection and extension of the ETM in situations of different dynamic forcing are also important themes that merit discussion in this section. Finally, the question of the frequency of the images in capturing ETM dynamics is also assessed.

Seasonal and Interannual Spatial Patterns
The seasonal and inter-annual variations of the turbidity of the Maroni estuary are analyzed in terms of spatial differences. Normalized SPM values were estimated as explained in Section 2.2.4 and averaged to calculate the location and extension of the ETM under different conditions (dry and rainy season) (Tables 2 and 3) and years (between 2013 and 2019) ( Figure 3). As stated above, the outlines of the ETM have been defined using normalized values above 0.6 for the averaged seasonal and inter-annual situations. Figure 3 depicts a synthesis of the averaged annual and seasonal situations. SPM values are affected by the combined influence of tidal currents, river discharge, winds and waves, such that instant SPM spatial patterns are quite heterogeneous.   Table 1. The pink dot represents the core of the ETM.
For the years 2019 and 2013, another area of maximum turbidity, quite distinct from the classical ETM, is observed upstream (Figure 3a-Upstream turbidity). The presence of an upstream turbidity area in the Maroni River could be related to fluvial sediment influx as it occurs in the rainy season when runoff values are high [44]. It should be noted, however, that the presence of this river turbidity zone for 2019 goes with the availability of only rainy season scenes for this year (Table 1). Therefore, the situation in 2019 is representative only of that year's rainy season. Following the spatial patterns found for the rainy and the dry season (Table 3), we could expect that the ETM during the dry season of 2019 will undergo a landward "tail" and "nose" shift. Conversely to that situation, the available scenes for 2013 correspond to both seasons (Table 1), with a better representation of the dry season for that year. Despite that, a river turbidity zone is still observed upstream (Figure 3a, 2013: Upstream turbidity). This same situation is, as already mentioned, also observed in 2014, when the ETM was present all along the west bank of the river (Figure 3a, 2014), between 7.4 ± 6.1 and 31.5 ± 14.9 km upstream of KP 0 (Table 3). SPM values inside the ETM during 2014 are lower than during 2013 (not shown) and more heterogeneous, and this could be related to the lower 2014 river discharge values compared to those of 2013 ( Figure 2b). Even though the dry and the rainy season are better represented by the years 2013 and 2019, respectively, these spatial patterns, depicted in Figure 3, allow us to evaluate the evolution of the location of the ETM over a period of six years, which corresponds to about 20% of the time necessary for a mud bank to go past a given location [53]. In this case, the years 2013-2019 are within the range of the first years of encroachment of the mud bank.  Table 1. The pink dot represents the core of the ETM.
The results for each of the years are depicted in Figure 3 and Table 2. The area of the ETM (Table 2) shows an increase in 2018 and 2019, with values of 69.4 ± 19.9 and 76.6 ± 28.7 km 2 , respectively. The ETM "tail" was found near the coast (KP 0) in the last three years of the study (2017-2019) (Figure 3a,b and Table 3). The shift of the averaged ETM "tail" towards the mouth of the estuary for those years coincides with the increase in ETM extension for the years 2018 and 2019. The ETM "nose" (Figure 3a and Table 3) also experienced a migration of 9.6 km after 2014. However, none of the forcing agents ( The direct exposure of the west bank to the prevailing waves from east to northeast could explain significant mud resuspension and alongshore mobilization, and hence the absence of an ETM. In 2014, the ETM showed a length of 40.9 km (from the "tail" to the "nose"), the longest in the course of the study period. This length could be related to the low discharge values in 2014. Low discharge is known to induce landward shifts of the ETM [52].
For the years 2019 and 2013, another area of maximum turbidity, quite distinct from the classical ETM, is observed upstream (Figure 3a-Upstream turbidity). The presence of an upstream turbidity area in the Maroni River could be related to fluvial sediment influx as it occurs in the rainy season when runoff values are high [44]. It should be noted, however, that the presence of this river turbidity zone for 2019 goes with the availability of only rainy season scenes for this year (Table 1). Therefore, the situation in 2019 is representative only of that year's rainy season. Following the spatial patterns found for the rainy and the dry season (Table 3), we could expect that the ETM during the dry season of 2019 will undergo a landward "tail" and "nose" shift. Conversely to that situation, the available scenes for 2013 correspond to both seasons (Table 1), with a better representation of the dry season for that year. Despite that, a river turbidity zone is still observed upstream (Figure 3a, 2013: Upstream turbidity). This same situation is, as already mentioned, also observed in 2014, when the ETM was present all along the west bank of the river (Figure 3a, 2014), between 7.4 ± 6.1 and 31.5 ± 14.9 km upstream of KP 0 ( Table 3). SPM values inside the ETM during 2014 are lower than during 2013 (not shown) and more heterogeneous, and this could be related to the lower 2014 river discharge values compared to those of 2013 (Figure 2b). Even though the dry and the rainy season are better represented by the years 2013 and 2019, respectively, these spatial patterns, depicted in Figure 3, allow us to evaluate the evolution of the location of the ETM over a period of six years, which corresponds to about 20% of the time necessary for a mud bank to go past a given location [53]. In this case, the years 2013-2019 are within the range of the first years of encroachment of the mud bank.
During both the rainy and the dry seasons, the averaged extension and location of the ETM show similar values (Tables 2 and 3). However, although the extension and ETM "tail" and "nose" are similar, the ETM core migrated landward from the rainy to the dry season (with the exception of the year 2013) ( Table 3). Sottolichio et al. [15] showed that the salt intrusion during the rainy season could go up to~5 km upstream from the mouth of the Maroni in the rainy season and~12 km during the dry season, these values being similar to those identified by Jouanneau and Pujos [33] (~5 and~11 km, respectively) over three decades ago. Even though these are rough estimates, the variation of the salt intrusion between both seasons could explain the landward migration of the core of the ETM.

Tidal Influence
Averaged SPM-normalized values were estimated to evaluate the tidal influence on the ETM's spatial location. This influence is analyzed in Figure 4 for both seasons.
During neap tides, the ETM is located between KP 6.2 and KP 27.2 (Figure 4a,b), and increasing SPM-normalized concentrations from~0.62 ± 0.31 to 0.93 ± 0.76 are observed from the upstream part of the ETM to its downstream part. The core location of the ETM during the dry season is 14.8 ± 2.6 km compared to the mean location of 14.1 ± 6.4 km in the rainy season, thus highlighting a more stable ETM core position during the dry season (Figure 4a,b). High SPM-normalized concentrations are observed all along the ETM during the rainy season, with average values of~0.76 ± 0.39. A river turbidity maximum is observed during this season at the west bank of the river 29.5 km from KP 0 (Figure 4b). High SPM concentrations are found in this river turbidity maximum, with averaged SPM-normalized values of 0.87 ± 0.62 (Figure 4b).
During spring tides, the ETM is located between KP 0 and KP 25.3 (Figure 4c,d), and it is longer than during neap conditions (Figure 4a,b). The dry season is characterized by an ETM with a core position of 10.2 ± 5.8 km, whereas the rainy season shows a position of 9.9 ± 6.9 km. SPM-normalized concentrations exhibit a longitudinal variation, with an increase from the west in the downstream area (~0.76 ± 0.49) (between KP0 and KP10) to the east in the area surrounding the downstream islands (~0.81 ± 0.52) (islands between the KP 1 0 and KP20) (Figure 4c). High SPM-normalized concentrations are found in the coastal area of the estuary and outside the ETM, with mean values of~0.43 ± 0.26 (Figure 4c). The ETM during the rainy season has a length of 22.3 km, and up to 3.1 km of the ETM length is located in the coastal area of the estuary (after KP 0) (Figure 4d). Averaged SPM-normalized concentrations are uniform at both banks of the estuary, with a mean value of 0.72 ± 0.21, whereas averaged SPM at the coastal area of the estuary and seaward of the ETM (after KP 0) shows a value of 0.47 ± 0.27 (Figure 4d).
Tidal currents seem to play a major role in the spatial variability of the ETM during both seasons. However, river discharge seems to modify the ETM "nose" and "tail" positions, which, for both neap and spring tides, are downstream for high river discharge conditions. These results are consistent with those obtained by Sottolichio et al. [15], who showed that SPM concentrations were influenced more by tidal currents during the dry season, whereas they were influenced more by river discharge during the rainy season.
Remote Sens. 2020, 12, x 12 of 23 Tidal currents seem to play a major role in the spatial variability of the ETM during both seasons. However, river discharge seems to modify the ETM "nose" and "tail" positions, which, for both neap and spring tides, are downstream for high river discharge conditions. These results are consistent with those obtained by Sottolichio et al. [15], who showed that SPM concentrations were influenced more by tidal currents during the dry season, whereas they were influenced more by river discharge during the rainy season. During neap tides, SPM values inside the ETM zone are higher for the rainy season, 170.9 mg/L, than for the dry season, 113.76 mg/L. This difference could be explained by the inter-annual variation of river discharge, as during the dry season the decrease in discharge would allow the development of a well stratified water column that can impede resuspension [54].

River Discharge
As the mixing zone normally occurs near the limit of the salt intrusion [55,56] and the salt intrusion depends on the magnitude of river discharge and tidal range [57], the influence of river discharge on the ETM's spatial location is presented in Figure 5. During neap tides, SPM values inside the ETM zone are higher for the rainy season, 170.9 mg/L, than for the dry season, 113.76 mg/L. This difference could be explained by the inter-annual variation of river discharge, as during the dry season the decrease in discharge would allow the development of a well stratified water column that can impede resuspension [54].

River Discharge
As the mixing zone normally occurs near the limit of the salt intrusion [55,56] and the salt intrusion depends on the magnitude of river discharge and tidal range [57], the influence of river discharge on the ETM's spatial location is presented in Figure 5. manifested by low values of averaged normalized SPM where the bathymetry is deep [36], whereas values of ~0.52 ± 0.22 are present at the shallower east and west sides of the channel.
During high flow conditions (exceeding 2000 m 3 /s, as in the case of the study of Berthois and Hoorelbeck [58]), the ETM is located between the mouth of the estuary (KP 0) and the downstream zone of the islands (KP 15.4) (Figure 5b). Averaged SPM values exhibit a heterogeneous distribution (Figure 5b) as maximum values are observed at the east bank of the river (between KP 7.5 and KP 10), over an area with shallow bathymetry (Figure 1c). The rest of the ETM exhibits an almost constant SPM-averaged value (~0.67 ± 0.35). During high river discharge conditions, both the coastal and the upstream areas of the estuary (outside the ETM) exhibit high normalized SPM, with averages of 0.36 ± 0.31 and 0.55 ± 0.27, respectively. In the coastal area, high SPM values are observed in the area with the deepest bathymetry as reported by Jolivet et al. [36]. The annual variation of river discharge, which reached a minimum of 146 m 3 /s and a maximum of 6670 m 3 /s for the period of study (Figure 2b), is consistent with the seasonal variation in the location of the ETM. Under high river discharge, relatively high SPM values are observed at the upstream part of the estuary, especially at the west bank (Figure 5b). Even though the suspended sediment supply of French Guiana rivers [42], and the Maroni River in particular, is low, the sharp increase in anthropogenic activities over the last decade, notably gold mining in this river's catchment, especially along the border with Suriname [44], could explain the relatively high values along the west bank of the river during these high river flow conditions.

SPM and Variability of ETM Extension Due to External Conditions
SPM values inside the estuary are the result of the combination of different sources of dynamic forcing; nevertheless, the previous results correspond to a synthesis for different dynamic conditions.   (Figure 5a), where the bathymetry is relatively shallow (Figure 1c). Averaged SPM-normalized values inside the ETM are very homogeneous, with a mean of~0.92 ± 0.47, which corresponds to a mean value of 146.4 mg/L (estimated as the mean value of the pixels inside the ETM zone using the corresponding non-normalized SPM maps, see Section 2.2.5). In the coastal area, the influence of the channel is detected, manifested by low values of averaged normalized SPM where the bathymetry is deep [36], whereas values of~0.52 ± 0.22 are present at the shallower east and west sides of the channel.
During high flow conditions (exceeding 2000 m 3 /s, as in the case of the study of Berthois and Hoorelbeck [58]), the ETM is located between the mouth of the estuary (KP 0) and the downstream zone of the islands (KP 15.4) (Figure 5b). Averaged SPM values exhibit a heterogeneous distribution (Figure 5b) as maximum values are observed at the east bank of the river (between KP 7.5 and KP 10), over an area with shallow bathymetry (Figure 1c). The rest of the ETM exhibits an almost constant SPM-averaged value (~0.67 ± 0.35). During high river discharge conditions, both the coastal and the upstream areas of the estuary (outside the ETM) exhibit high normalized SPM, with averages of 0.36 ± 0.31 and 0.55 ± 0.27, respectively. In the coastal area, high SPM values are observed in the area with the deepest bathymetry as reported by Jolivet et al. [36].
The annual variation of river discharge, which reached a minimum of 146 m 3 /s and a maximum of 6670 m 3 /s for the period of study (Figure 2b), is consistent with the seasonal variation in the location of the ETM. Under high river discharge, relatively high SPM values are observed at the upstream part of the estuary, especially at the west bank (Figure 5b). Even though the suspended sediment supply of French Guiana rivers [42], and the Maroni River in particular, is low, the sharp increase in anthropogenic activities over the last decade, notably gold mining in this river's catchment, especially along the border with Suriname [44], could explain the relatively high values along the west bank of the river during these high river flow conditions.

SPM and Variability of ETM Extension Due to External Conditions
SPM values inside the estuary are the result of the combination of different sources of dynamic forcing; nevertheless, the previous results correspond to a synthesis for different dynamic conditions. To better analyze these differences and the influence of these dynamic conditions, SPM values under different tidal, river discharge, wind and wave conditions are presented in Figure 6. Contrasting dynamics are observed for neap tides during ebb and flood (Figure 6b,c,f). During ebb tides, the ETM is well developed (Figure 6b,f) and is located in the area between the downstream islands (KP 20and KP 0). High values of SPM are observed at the west bank of the estuary, which has been identified as the ebb channel [16], especially for 2018, when the length of the ETM was considerably longer than that of 2014 (Figure 6b,f). This is consistent with the results of Jolivet et al. [36]: the dominant direction of transport of sediments from this river is westwards along the Suriname coast. During flood, a clear ETM at the surface of the water column is absent in the estuary (Figure 6c). This would suggest a stratification of the water column for this scene of the dry season. As a consequence of this stratification, high resuspension of SPM is impeded [54] and low SPM values are observed at the surface of the water column ( Figure 6c). Therefore, flood tides are characterized by low SPM values in comparison to ebb tides, as observed by Asp et al. [59] in a river on the coast of Para State in North Brazil, also under the influence of the Amazon. No clear influence of wind velocity values is detected when the different scenes are compared ( Figure 6).
During spring tides, SPM-averaged values inside the ETM are higher than during neap tides (Figures 4 and 6) for both seasons. During the dry season, these values reach 160.6 mg/L, whereas they attain 195.6 mg/L during the rainy season. Likewise, at neap tides, ebb tides present a clearly delineated ETM located at the mouth of the estuary (Figure 6e). Even though the river discharge is 3.7 times lower in Figure 6a (Figure 6a). This could be related to the river discharge, as even though the daily value for that date is relatively low (545 m 3 /s), an increase in discharge is observed in the days prior to 24 December 2014 (Figure 2a). The increase in river flow could have enhanced suspended sediment remobilization after a period of low river discharge, thus leading to the high SPM along the river and estuary when compared to the other ebb scene (Figure 6e). SPM at the coastal area in front of the estuary is also characterized by high values (averaged SPM of 96.9 mg/L) when compared to neap tide values (averaged SPM of 44.8 mg/L) ( Figure 6). During these tidal conditions, high resuspension characterizes the sand-filled estuary mouth with its various shoals and banks (Figure 6a,e). The area surrounding the vegetated islands ( Figure 1c) exhibits relatively high flood-tide SPM values (~110 mg/L) and the flood channel is clearly identified at the mouth of the estuary (Figure 6d). In the coastal area of the estuary, high resuspension prevails, especially at the shallow west side of the flood channel [36].
The impact of the river on the mud bank is also illustrated in Figure 6. River discharge impacts the subtidal resuspension of muddy deposits at the east bank in the coastal area of the estuary (Figure 6b,c,f). This confirms the influence of large rivers such as the Maroni on mud banks, as the river discharge, even in low-flow conditions (Figure 6c,d,f), seems to generate a hydraulic-groyne effect that slows down mud-bank migration through mud fluidization [31], inducing muddy accretion on the updrift sector of the estuary mouth [36]. However, such impedance of mud-bank migration across the Maroni estuary mouth would not only impact migration rates, but could also lead to some mud being displaced up the estuary [14].  SPM dynamics in the estuary have been observed ( Figure 6) to be largely variable, with an ETM that is not always well developed (Figure 6a,c). As shown in Tables 2 and 3, the ETM extension and the locations of its "tail", "nose" and "core" exhibit relatively high standard deviations relative to their averaged values for some of the years and seasons. To enhance our analysis of the variation of the ETM extension in response to the different dynamic forcings, we confronted these averaged values to values of river water discharge, tide, wave height and wind corresponding to each image.

Medium-Term (2013-2019) Variation: The Influence of Periodic Amazon-Derived Mud Banks
Although an overall increase in SPM values in the Maroni River has been observed since 2009 [44], the SPM increase in the upstream area seems too small to justify the increase in ETM extension between 2013 and 2019 (Table 2). Moreover, river discharge seasonality alone cannot explain this extension, as mean SPM values for 2016, for example, were ~30 mg/L for all upstream Maroni stations [44]. Therefore, the sediment source is likely to be related to the mud bank present in the study area since 2011, reflecting a situation similar to that reported by Orseau et al. [14] for the smaller Mahury estuary in eastern French Guiana. This influence is observed in the increasingly more seaward position of the ETM "tail" (Figure 3b and Table 3). This seaward position could be related to the mudbank migration process. The back (trailing edge) of the bank is composed of semi-consolidated mud that is eroded and fluidized by waves, tidal currents and wind-induced geostrophic currents [60,61]. Tidal current speeds at the east bank of the Maroni estuary mouth exceed 0.75 m/s during ebb and 1 m/s during flood [36], which correspond to velocities that can initiate erosion [62,63], in addition to wave attack [27]. The erosion of the back of the bank could provide the extra amount of mud for injection into the Maroni's ETM, as also suggested by Jolivet et al. [36].
In support of the fore-going observations on the location and extension of the ETM (Tables 2 and  3), SPM values have increased downstream since 2013 (Figure 8). At station S1, located at the mouth  figure), neither of which seems to be effective in determining ETM extensions. However, given the moderate relationship between river discharge and ETM, it is clear that this dynamic agent is not the sole driver of ETM extension. For those images corresponding to a discharge of < 500 m 3 /s, the ETM extension ranges from 32.2 to 91.7 km 2 (Figure 7a). The same is observed when images with a discharge > 2000 m 3 /s are considered; the range variability of ETM extension goes from 40.3 to 63.7 km 2 (Figure 7a). In the case of the tide, a clearer tendency emerges (R = −0.31 and p = 0.015), especially when the tidal range is < 2 m (Figure 7b). As the previous section shows that both tide and flow rate highly modulate the variations and location of the ETM, its extension is compared to the flow rate for those images presenting flood and ebb conditions (Figure 7c,d). ETM extension for both tidal situations varies with the flow rate, increasing when this value does (Figure 7c,d). However, this tendency is not that clear when discharge values are higher than 2500 m 3 /s for both tidal conditions, and lower than 500 m 3 /s, especially for flood conditions (Figure 7c,d).
In the case of the images associated with ebb conditions and high river discharge, the almost constant value of~50 km 2 for the extension of the ETM could be related to both high discharge and ebb tide currents pushing the ETM seaward, inducing a limitation of its development to the mouth of the estuary. For flood conditions and low discharge, it seems that images capturing the early hours of the ebb tide depict a larger extension that those closer to low tide (figure not shown). However, it is observed that the dynamics affecting the ETM extension are really complex, as no clear relationship among them comes out, even when they are combined. This could be related to the stratification of the water column and how it affects what can be detected from the water surface. Although an overall increase in SPM values in the Maroni River has been observed since 2009 [44], the SPM increase in the upstream area seems too small to justify the increase in ETM extension between 2013 and 2019 (Table 2). Moreover, river discharge seasonality alone cannot explain this extension, as mean SPM values for 2016, for example, were~30 mg/L for all upstream Maroni stations [44]. Therefore, the sediment source is likely to be related to the mud bank present in the study area since 2011, reflecting a situation similar to that reported by Orseau et al. [14] for the smaller Mahury estuary in eastern French Guiana. This influence is observed in the increasingly more seaward position of the ETM "tail" (Figure 3b and Table 3). This seaward position could be related to the mud-bank migration process. The back (trailing edge) of the bank is composed of semi-consolidated mud that is eroded and fluidized by waves, tidal currents and wind-induced geostrophic currents [60,61]. Tidal current speeds at the east bank of the Maroni estuary mouth exceed 0.75 m/s during ebb and 1 m/s during flood [36], which correspond to velocities that can initiate erosion [62,63], in addition to wave attack [27]. The erosion of the back of the bank could provide the extra amount of mud for injection into the Maroni's ETM, as also suggested by Jolivet et al. [36].
In support of the fore-going observations on the location and extension of the ETM (Tables 2  and 3 [29], the subtidal part of the mud bank (which forms the bulk of a mud bank [14,38]) has been crossing the estuary mouth since 2017. Different features, such as the increase in SPM, extension of the ETM, and the clear seaward evolution of the ETM "tail" since 2017 (Table 3), could then be explained by the migration of the mud bank across the estuary mouth. As the evolution and spatial patterns of the ETM are mainly modified by river discharge and by tide-induced transport, these would be the principal drivers of mud intrusion into the estuary, as also suggested by Orseau et al. [14]. Landward mud ingression has also been identified as the main source of sediments in other estuaries under the influence of the Amazon River plume [59], confirming the influence of muddy deposits from the Amazon in the sediment dynamics and potential long-term infill of estuaries in northern South America. It may be expected that, during inter-bank phases separating mud-bank phases, these Amazon-influenced estuaries would have a functional regime similar to that of other tropical estuaries elsewhere, with minimal influence of outside mud-sourcing.
ingression has also been identified as the main source of sediments in other estuaries under the influence of the Amazon River plume [59], confirming the influence of muddy deposits from the Amazon in the sediment dynamics and potential long-term infill of estuaries in northern South America. It may be expected that, during inter-bank phases separating mud-bank phases, these Amazon-influenced estuaries would have a functional regime similar to that of other tropical estuaries elsewhere, with minimal influence of outside mud-sourcing.

Capacity of Satellite Images in Contributing to Analysis of Estuarine Dynamics
The spatial resolution of data from the new generation of satellite images for ocean color, such as Landsat 8 (30-15 m; used in this study), has enabled us to attempt to characterize the spatial evolution of the ETM, using for the first time six-year time series data in a tropical estuary. Although the methods used in the present study have enabled an exploratory characterization of the seasonality of the ETM and the influence of tidal currents and of an external source of mud, namely, an Amazonderived mud bank, the temporal resolution of Landsat 8 satellite images (16 days) is a limiting factor in gaining insight into small temporal-scale processes. The number of available images (Table 1) is restricted due to the high cloud cover in this equatorial area. This limited temporal coverage could affect our comprehension of the seasonal differences in the study area, as only a small number of rainy season images can be exploited (Table 1), a period when high river discharge considerably impacts ETM dynamics. The inter-annual variation of the ETM is, thus, highly hinged on the available scenes, the dry season being much better represented than the rainy season (Table 1). This, therefore, constitutes a limitation in the analysis. For 2019, only rainy season scenes were selected. An upstream river turbidity maximum was detected that year. This maximum appears also in 2013 and 2014, indicating that river influence becomes appreciable only under certain dynamic conditions ( Figure  3a and Figure 7).

Capacity of Satellite Images in Contributing to Analysis of Estuarine Dynamics
The spatial resolution of data from the new generation of satellite images for ocean color, such as Landsat 8 (30-15 m; used in this study), has enabled us to attempt to characterize the spatial evolution of the ETM, using for the first time six-year time series data in a tropical estuary. Although the methods used in the present study have enabled an exploratory characterization of the seasonality of the ETM and the influence of tidal currents and of an external source of mud, namely, an Amazon-derived mud bank, the temporal resolution of Landsat 8 satellite images (16 days) is a limiting factor in gaining insight into small temporal-scale processes. The number of available images (Table 1) is restricted due to the high cloud cover in this equatorial area. This limited temporal coverage could affect our comprehension of the seasonal differences in the study area, as only a small number of rainy season images can be exploited (Table 1), a period when high river discharge considerably impacts ETM dynamics. The inter-annual variation of the ETM is, thus, highly hinged on the available scenes, the dry season being much better represented than the rainy season (Table 1). This, therefore, constitutes a limitation in the analysis. For 2019, only rainy season scenes were selected. An upstream river turbidity maximum was detected that year. This maximum appears also in 2013 and 2014, indicating that river influence becomes appreciable only under certain dynamic conditions (Figures 3a and 7).
The recently-launched Sentinel 2 missions, with a spatial resolution of 10 m and a time resolution of five days, will allow for greater insight into the ETM, as the combination of data from this satellite and from Landsat 8 will generate a much higher temporal resolution. This higher resolution should allow for a better characterization of occasional (or semi-permanent) features such as the river turbidity maximum observed upstream in the Maroni. It should also enable us to identify the influence of different drivers of mud-bank migration and sediment resuspension, such as oceanic currents and trade-wind waves [27,28,64], as well as better characterization of resuspension processes at the sand-filled estuary mouth of the Maroni and at the back of the mud bank. These processes cannot presently be explained using the Landsat 8 series only. One current limitation of reliance on satellite data is that only surface SPM patterns are apprehended. It would, thus, also be necessary to compare these surface results with in-situ values of SPM in the water column. This comparison should enable a better explanation of the dynamics of the ETM, comparing them also to the tidal excursion between neap and spring tides and the influence of the limit of salt intrusion related to the stratification of the water column in tropical estuaries under the influence of the Amazon. It would also be necessary to compare these results with those derived from hydro-sedimentary models that provide a better description of mixing processes for different dynamic forcing conditions.

Conclusions
Understanding the spatial and temporal variability of the ETM is of primary importance in tropical estuaries where studies focusing on the sediment dynamics are rare, whereas development needs, including estuarine navigation, port development and monitoring of river-to-coast pollution, are rising with population increase. The dynamics of the Maroni estuary are representative of that of tropical-equatorial estuaries subject to strongly seasonal discharge, but this study also illustrated the idiosyncratic character of Guiana Shield estuaries, generally poor in intrinsic mud supply from their crystalline catchments, but periodically subject to extraneous mud under the influence of the Amazon River and its massive muddy discharge into the coastal ocean of northern South America. Many of these Guiana estuaries are also undergoing changes in sedimentation patterns caused by human activities and deforestation upstream, while port facilities are being actively developed. Our study illustrated the potential of time series obtained using high-resolution sensors, such as Landsat 8-OLI, for characterizing the ETM of a tropical estuary and the periodic influence of Amazon mud banks on this feature, notwithstanding the complexity of ETM detection using surface SPM values without a consideration of depth concentration gradients. Satellite data are, indeed, proving to be an effective tool for monitoring coastal areas with very turbid waters, such as the Guianas coast [27,29,38,65,66] where field monitoring is a challenging task.
SPM values were obtained using a general algorithm that allows the estimation of four ranges of SPM magnitude. The estimates were normalized and image pixels that presented a value >0.6 were retained to yield an envelope of values that was used to identify the ETM zone. This methodology allowed for the definition of the ETM under different dynamic conditions and the characterization of tropical seasonality, inter-annual variations and the influence of tides and river discharge on the ETM. This study evidenced an ETM variation influenced by the well-marked seasonal river discharge and the neap-spring tidal cycle, especially during the dry season, as observed in Sottolichio et al. [15]. As a result, the ETM is located in the estuary mouth during high river discharge conditions and it shifts inwards towards the middle estuary during low flow conditions. Neap-spring tidal range variations result in a push of the ETM closer to the mouth of the estuary under spring-tide conditions or even outside the mouth during spring tides in the rainy season.
The influence of the migration of a mud bank present at the mouth of the Maroni since 2011 was also observed. The ETM is pushed towards the coastal area of the estuary as the mud bank migrates across its mouth. An increase in SPM values is also observed within the ETM zone, as well as an increase in ETM extension. This study confirms observations made in a previous study in the area that highlighted the role of alongshore migrating Amazon-derived mud banks as purveyors of sediments that intrude into the estuaries of the Guianas [14]. Thus, an important source of mud for the Maroni estuary, a river with a catchment considered as a rather poor purveyor of mud, is likely to be offshore and related to the periodic passage of mud banks "en route" from the Amazon to the Orinoco. To determine to which extent Amazon mud scales relative to the intrinsic mud derived from the Maroni catchment in the constitution of the estuary's ETM would require a set of observations and analyses conducted during an inter-bank period, when coastal waters under the influence of the Amazon are characterized by relatively low turbidity in the vicinity of the Maroni.
The results obtained provide new perspectives for the analysis of tropical estuaries through satellite data and estimates derived thereof, such as SPM, and for the identification of the spatial and temporal evolution of the ETM in relation to tidal, seasonal and inter-annual regimes. Even though Landsat-8 presents a low revisit frequency, the methodology used and the results obtained are very promising for studying coastal and estuarine areas. The combination of data from this satellite with data from the Sentinel-2A and -2B satellites recently launched by the European Space Agency will make for considerable enhancement of knowledge on the inter-annual and seasonal variability of ETM zones. The application of this analysis to Sentinel 2 data will enable better characterization of the ETM during the rainy season and provide significant advantages in the observation of short timescale processes. Moreover, augmentation of the quantity of data would enable a determination Remote Sens. 2020, 12, 2173 20 of 23 of the influence of wind-induced waves and currents on SPM values and ETM evolution, the main drivers of coastal suspended sediment dynamics on the Guianas coast [27,64]. Better knowledge of the short-term processes could benefit from integration with in-situ measurements and hydro-sedimentary models, providing a more thorough description of mixing processes for different dynamic forcing conditions, and a valuable aid to estuarine management scenarios under increasing human pressures and climate change.