Synergy between Low Earth Orbit (LEO)—MODIS and Geostationary Earth Orbit (GEO)—GOES Sensors for Sargassum Monitoring in the Atlantic Ocean

: Since 2011, massive stranding of the brown algae Sargassum has regularly affected the coastal waters of the West Caribbean, Brazil, and West Africa, leading to heavy environmental and socio-economic impacts. Ocean color remote sensing observations as performed by sun-synchronous satellite sensors such as MODIS (NASA), MERIS (ESA), or OLCI (ESA/Copernicus) are used to provide quantitative assessments of Sargassum coverage through the calculation of indices as the Alternative Floating Algae Index (AFAI). Sun-synchronous sensors usually provide at best one daytime observation per day of a given oceanic area. However, such a daily temporal revisit rate is not fully satisfactory to monitor the dynamics of Sargassum aggregation due to their potentially signiﬁcant drift over the course of the day as a result of oceanic currents and sea surface wind stress. In addition, the sun glint and the presence of clouds limit the use of low earth orbit observations, especially in tropical zones. The high frequency sampling provided by geostationary sensors can be a relevant alternative approach in synergy with ocean color sun-synchronous sensors to increase the temporal resolution of the observations, thus allowing efﬁcient monitoring of Sargassum dynamics. In this study, data acquired by a geostationary satellite sensor located at 36,000 km from Earth, namely GOES-16 (NASA/NOAA), which was primarily designed for meteorology applications, are analyzed to investigate the Sargassum dynamics. The results demonstrate that a GOES-16 hourly composite product is appropriate to identify Sargassum aggregations using an index commonly used for vegetation monitoring, namely NDVI (Normalized Difference Vegetation Index). It is also shown that GOES hourly observations can signiﬁcantly improve the simulated drift obtained with a transport circulation model, which uses geostrophic current, wind, and waves. This study thus highlights the signiﬁcant relevance of the effective synergy between sun-synchronous and geostationary satellite sensors for characterizing the Sargassum dynamics. Such a synergy could be summarized as follows:


Introduction
Unprecedented massive stranding of the brown macro-algae Sargassum fluitans and S. natans (hereinafter Sargassum) has been observed along the coastline of French Guyana, the Antilles Islands, and the Caribbean Sea over the last decade. These pelagic algae float in the upper surface layer in offshore waters and can cluster into large aggregations that are transported by currents over long distances across the Atlantic Ocean. Thus, satellite data are highly convenient to determine the algae spatial distribution. A spectral index called Maximum Chlorophyll Index (MCI) was first defined by Gower et al. [1] for the Medium Resolution Imaging Spectrometer (MERIS) satellite sensor (300 m resolution). Such an index relies on the spectral peak of the water leaving radiance at 709 nm that is mostly induced by the Sargassum algae optical signature. It was shown that such a peak indicates the occurrence of a high surface concentration of chlorophyll a, and allows extensive areas of pelagic vegetation to be detected [1,2]. The MCI index was previously used to study the spatial distribution of Sargassum in the Gulf of Mexico and in the western Atlantic waters [3,4]. The Sentinel-3/OLCI (Ocean and Land Color Instrument) sensor was launched in 2016 to ensure the continuity of MERIS observations. OLCI was thus designed to provide spatial resolution and spectral bands similar to MERIS [5].
Hu [6] proposed a Floating Algae Index (FAI) to analyze the data acquired by the Moderate Resolution Imaging Spectroradiometer (MODIS) with the objective of detecting and tracing blooms of Ulva prolifera macroalgae species in the Yellow Sea near Qingdao, China [7,8]. Since FAI was proposed based on the vegetation red-edge reflectance observed between 675 and 750 nm, such an index could be relevant to detect any floating vegetation including Sargassum [9]. However, both Sargassum slicks and clouds show high FAI values. Considering the difficulty of performing a robust cloud masking, Wang and Hu [10] suggested the Alternative Floating Algae Index (AFAI), which better separates Sargassum from clouds, using data measured in alternative spectral bands, namely 667 nm, 748 nm, and 869 nm.
High spatial resolution sensors which typically have a resolution of between 10 m and 60 m such as Sentinel-2/MSI (MultiSpectral Instrument) are also relevant to detect Sargassum. Ody et al. [11] recently adapted the AFAI index to MSI spectral features using the bands at 665 nm, 833 nm, and 940 nm, thus providing the so-called Modified Floating Algae Index (MFAI).
In addition to Sargassum detection, some algae indices allow the Sargassum fractional cover within the pixels to be quantified and converted into biomass per unit area. Wang and Hu [10] proposed a linear relationship between the AFAI and the fraction of pixels covered by Sargassum. To estimate the AFAI value for a pixel covered by 100% of Sargassum, they used bucket measurements and radiative transfer simulations. Wang et al. [12] extended this approach to estimate biomass using field measurements.
Ocean color sun-synchronous sensors are then able to provide indices that can be linked to quantitative parameters such as fractional coverage and biomass of Sargassum. However, these sensors provide at best daily retrievals and sometimes fewer in very cloudy conditions, as it frequently happens over most of the Tropical Atlantic Ocean. Thus, the temporal variation of Sargassum and hence their dynamics remain difficult to evaluate. Yet, Sargassum are known to drift rapidly (as much as 17 km over the course of one day) due to ocean currents and sea surface wind stress [11]. Geostationary satellite observations could be relevant to monitor Sargassum thanks to their high temporal revisit rate, typically every 15 to 30 min. The analysis of data acquired by the Korean geostationary sensor GOCI together with Lagrangian particle modeling demonstrated the feasibility of tracing the trajectories of pelagic Sargassum aggregation using this type of satellite orbit [13]. GOCI radiometric sensitivity is specifically dedicated to measuring low level signals such as those encountered in oceanic waters. Unfortunately, the GOCI orbit is not designed to observe the Atlantic Ocean. Therefore, there is currently a lack of geostationary sensors designed for ocean color applications to observe areas in the Atlantic Ocean where Sargassum strongly proliferate [13].
Remote Sens. 2021, 13, 1444 3 of 16 GOES-16 ABI is a geostationary sensor that provides data over the Atlantic Ocean. Although it is not designed for observing ocean targets, it may be possible to use it to monitor Sargassum patterns. GOES provides only 3 spectral bands in the visible to nearinfrared domain, i.e., blue, red, and Near Infra-Red (NIR) with a revisit rate of 15 min (32 images per day). GOES spatial resolutions are 1 km, 500 m, and 1 km for the blue, red, and NIR wavelengths respectively. Thus, the AFAI cannot be calculated using GOES data because the required bands are missing. However, the well-known NDVI index [14,15], which is commonly used to monitor land vegetation, can be calculated.
In this study, the objective is to demonstrate the feasibility of quantifying the temporal variation of Sargassum from high-temporally resolved data acquired by a geostationary sensor designed for weather forecast applications, namely GOES-16 ABI. GOES is used in synergy with the MODIS sun-synchronous sensor to characterize Sargassum dynamics.
The paper is organized as follows: The study area, namely the eastern part of Saint Lucia Island in the western area of the North Atlantic Ocean, is first described (Section 2). The main features of the satellite data and the proposed methodology using GOES data in synergy with MODIS observations are then outlined (Section 2). The spatial variation of MODIS AFAI Level-2 products and the time variation of NDVI product derived from GOES are presented in Section 3. Ways in which the quality of GOES raw data could be improved to enhance its potential to observe Sargassum are discussed in Section 4. Finally, the analysis of GOES data for improving the forecast of Sargassum transport and stranding is performed by comparing it with estimated drift velocity from an ocean circulation model (Section 4).

Study Area
The study area is located in the Caribbean Sea in the North Atlantic Ocean. The frequent influx of Sargassum that reached the French Antilles in the past, especially between 2015 and 2018, represents a major threat to the ecology, health, and economy of these islands [16]. These islands can each receive an average value of 150,000 m 3 of seaweed during major stranding episodes [17]. All the bays located on the east coasts are affected by the Sargassum stranding. The west coasts can also be affected by Sargassum presence due to ocean currents. Our specific study area is located west of Martinique and Saint Lucia Islands (13-17 • N and 60-63 • W, white box in Figure 1). The selected temporal period for the data analysis is the month of August 2018, which corresponds to a time where Sargassum was abundant as revealed by MODIS data [11,18]. GOES-16 ABI is a geostationary sensor that provides data over the Atlantic Ocean. Although it is not designed for observing ocean targets, it may be possible to use it to monitor Sargassum patterns. GOES provides only 3 spectral bands in the visible to nearinfrared domain, i.e., blue, red, and Near Infra-Red (NIR) with a revisit rate of 15 min (32 images per day). GOES spatial resolutions are 1 km, 500 m, and 1 km for the blue, red, and NIR wavelengths respectively. Thus, the AFAI cannot be calculated using GOES data because the required bands are missing. However, the well-known NDVI index [14,15], which is commonly used to monitor land vegetation, can be calculated.
In this study, the objective is to demonstrate the feasibility of quantifying the temporal variation of Sargassum from high-temporally resolved data acquired by a geostationary sensor designed for weather forecast applications, namely GOES-16 ABI. GOES is used in synergy with the MODIS sun-synchronous sensor to characterize Sargassum dynamics.
The paper is organized as follows: The study area, namely the eastern part of Saint Lucia Island in the western area of the North Atlantic Ocean, is first described (Section 2). The main features of the satellite data and the proposed methodology using GOES data in synergy with MODIS observations are then outlined (Section 2). The spatial variation of MODIS AFAI Level-2 products and the time variation of NDVI product derived from GOES are presented in Section 3. Ways in which the quality of GOES raw data could be improved to enhance its potential to observe Sargassum are discussed in Section 4. Finally, the analysis of GOES data for improving the forecast of Sargassum transport and stranding is performed by comparing it with estimated drift velocity from an ocean circulation model (Section 4).

Study Area
The study area is located in the Caribbean Sea in the North Atlantic Ocean. The frequent influx of Sargassum that reached the French Antilles in the past, especially between 2015 and 2018, represents a major threat to the ecology, health, and economy of these islands [16]. These islands can each receive an average value of 150,000 m 3 of seaweed during major stranding episodes [17]. All the bays located on the east coasts are affected by the Sargassum stranding. The west coasts can also be affected by Sargassum presence due to ocean currents. Our specific study area is located west of Martinique and Saint Lucia Islands (13-17° N and 60-63° W, white box in Figure 1). The selected temporal period for the data analysis is the month of August 2018, which corresponds to a time where Sargassum was abundant as revealed by MODIS data [11,18].

Satellite Data
MODIS is a satellite sensor launched by NASA [19] in 1999, on board the Terra satellite and in 2002, on board the Aqua satellite. MODIS allows the entire Earth to be observed every 1 to 2 days. MODIS acquires data in 36 spectral bands ranging from 0.4 µm to 14.4 µm.
Observations are performed at various spatial resolutions (2 bands at 250 m, 5 bands at 500 m, and 29 bands at 1 km). In this study, the bands at 667 nm, 748 nm, and 869 nm at a 1-km resolution of the Terra satellite are exploited to calculate AFAI. A composite RGB image acquired on 8 August 2018 14:35 UTC is shown in Figure 2a. That image covers the Islands from the North to the South: Montserrat, Guadeloupe, Dominica, Martinique, and Saint Lucia.
RGB image acquired on 8 August 2018 14:35 UTC is shown in Figure 2a. That image covers the Islands from the North to the South: Montserrat, Guadeloupe, Dominica, Martinique, and Saint Lucia.
The Advanced Baseline Imager (ABI) is the primary instrument on the GOES-R Series for imaging Earth's weather, oceans, and environment. GOES-ABI is a multi-channel passive imaging radiometer designed to observe the western hemisphere and provides relevant information for the analysis of land surfaces, atmosphere, and cloud cover. GOES-ABI observes the Earth in 16 large spectral bands ranging from 0.47 µ m to 13.3 µ m, including two channels in the visible spectrum, four channels in the near-infrared range, and 10 channels in the infrared range. GOES-ABI improves by a factor of 3, 4, and 5 the spectral, spatial, and temporal resolutions provided by the previous generation of GOES sensors, respectively. ABI concurrently produces a full disk every 15 min. In the visible domain, ABI's spatial resolution is 1 km for the blue band (470 nm) and the NIR (860 nm) band and 500 m for the red band (640 nm). The Signal-to-Noise Ratio (SNR) is around 1200, 500, and 800 for the blue, red, and NIR bands, respectively, which is quite high compared to other meteorological sensors [20,21]. A false-color NIR-G-B composite acquired on 8 August 2018 at 14:30 UTC is shown in Figure 2b. That image covers the same area as MODIS at almost the same time. Clouds are located at the same place in both images.  The Advanced Baseline Imager (ABI) is the primary instrument on the GOES-R Series for imaging Earth's weather, oceans, and environment. GOES-ABI is a multi-channel passive imaging radiometer designed to observe the western hemisphere and provides relevant information for the analysis of land surfaces, atmosphere, and cloud cover. GOES-ABI observes the Earth in 16 large spectral bands ranging from 0.47 µm to 13.3 µm, including two channels in the visible spectrum, four channels in the near-infrared range, and 10 channels in the infrared range. GOES-ABI improves by a factor of 3, 4, and 5 the spectral, spatial, and temporal resolutions provided by the previous generation of GOES sensors, respectively. ABI concurrently produces a full disk every 15 min. In the visible domain, ABI's spatial resolution is 1 km for the blue band (470 nm) and the NIR (860 nm) band and 500 m for the red band (640 nm). The Signal-to-Noise Ratio (SNR) is around 1200, 500, and 800 for the blue, red, and NIR bands, respectively, which is quite high compared to other meteorological sensors [20,21]. A false-color NIR-G-B composite acquired on 8 August 2018 at 14:30 UTC is shown in Figure 2b. That image covers the same area as MODIS at almost the same time. Clouds are located at the same place in both images.

Environmental Data
The Sargassum aggregation dynamic is linked to environmental data such as ocean current and temperature. The ocean current velocity and direction are obtained from OSCAR provided by NASA. The "Third Degree Sea Surface Velocity" product is a five-day temporal composite. Such a product is a direct computation of global surface currents using satellite sea surface height, 10 m wind, and surface temperature [22]. Currents are calculated using a quasi-steady geostrophic model together with an eddy viscosity-based wind-driven geostrophic component and a thermal wind adjustment.

Methodology
MODIS Level-1B data (i.e., Top-of-the-Atmosphere (TOA) radiance) were provided by NASA (ladsweb.modaps.eosdis.nasa.gov, accessed on 5 April 2021). They were used to compute the AFAI following Wang and Hu's approach [10] (Equation (1)): where λ 1 = 667 nm, λ 2 = 748 nm, and λ 3 = 869 nm, which correspond to MODIS bands 13, 15, and 16. R(λ) is the corrected reflectance derived from MODIS observations. The TOA radiance was corrected for molecular (Rayleigh) scattering and gaseous absorption using NASA's Ocean Color Science Software (OCSSW-SeaDAS) package. The correction for the aerosol effects was not performed in this study to make sure that THE Sargassum optical signature WAs not removed from the data analysis. This is because the water reflectance is not necessarily close to zero in the NIR band when Sargassum are present, thus preventing the application of the black pixel atmospheric correction method commonly used to derive aerosol properties over open ocean waters [25]. Clouds were flagged using the OCSSW mask based on a threshold technique using a value of 0.030 at 2130 nm. Additional masking for clouds, sun glint, and cloud shadows was applied following Wang and Hu [10]. While the AFAI is very sensitive to Sargassum optical signature when using MODIS data, it is also sensitive to the water spectral reflectance, to potential residual sun glint contamination, and to the presence of aerosols. Thus, the presence of Sargassum cannot just be determined using a constant threshold value of AFAI. Therefore, the deviation of AFAI with respect to the residual AFAI has also been taken into account to improve detection. The residual AFAI is estimated based on surrounding Sargassum-free water. The AFAI background (AFAI bg ) can be represented as a sum of three components (Equation (2)): where AFAI 0 retains the large-scale variation of AFAI (including potential residual atmospheric effects), δAFAI w is the deviation due to local characteristics of Sargassum-free ocean reflectance, and δAFAI geom is the deviation induced by the local variation of viewing geometry.
Similarly, to Wang and Hu [10], the background AFAI that corresponds to Sargassumfree water was estimated here. A median filter in a 51 × 51 window was applied to the AFAI to estimate the AFAI background in the vicinity of each Sargassum pixel. The positive deviation of AFAI with respect to the background (δAFAI) was then computed to determine the presence of Sargassum.
The Level-1 TOA reflectance data measured by GOES-ABI were provided by the AERIS/ICARE Data and Services Center (www.icare.univ-lille.fr, Villeneuve d'Ascq Cedex, France, accessed on 7 April 2021). These data were used to compute NDVI because the GOES bands do not allow AFAI to be calculated. The extraction of the region of interest was operated for the study area shown in Figure 1, namely 13-17 • N and 60-63 • W. The NIR data were over-sampled by a factor of 2 to reach a 500-m resolution using a bilinear method [26]. Then, a mask on the cloud and on the land was generated using a threshold value in the NIR band (NIR reflectance >0.13). Four successive acquisitions were averaged each hour to derive an hourly product. Such a method not only improves the SNR and the data quality but also decreases the amount of cloudy pixels. Indeed, because of the speed at which clouds travel, they change position within an hour. As the cloud masking procedure is applied before averaging, cloud-free observations are favored in the hourly product.
NDVI (Equation (3)) is computed for all the pixels of the hourly composite image. Such an index detects the increase of the vegetation reflectance located between 640 nm (λ 1 ) and 860 nm (λ 2 ), including the red-edge region of the spectrum located between 700 and 750 nm: where f is the Sargassum fractional coverage, R s the Sargassum reflectance, and R w the Sargassum-free water reflectance. Similarly, to the AFAI calculation, the background NDVI (NDVI bg ) was also obtained by applying a median filter in a 51 × 51 window to the NDVI. The positive deviation of NDVI, noted δNDVI, is calculated by subtracting the background NDVI (NDVI bg ) from NDVI (Equation (5)): Such a subtraction corrects the NDVI for potential residual sun glint contamination and atmospheric effects, assuming the atmosphere is homogenous over the considered scene.
Unlike AFAI and because of its definition in Equation (3), NDVI is non-linearly linked to Sargassum sub-pixel fractional coverage. NDVI is rapidly saturated by the presence of Sargassum, therefore, NDVI cannot be used as quantitative information to derive Sargassum fractional coverage. An interesting synergistic strategy could be to primarily detect and quantify the abundance of Sargassum using AFAI calculated from daily MODIS data, and secondarily use the NDVI derived from GOES data to monitor the Sargassum dynamics.
The flowchart of the overall methodology used to combine MODIS and GOES data is shown in Figure 3. Remote Sens. 2021, 13, x FOR PEER REVIEW 7 of 17 Figure 3. Flowchart of the overall methodology used to exploit the synergy between MODIS and GOES data for the characterization of Sargassum occurrence and dynamics.

AFAI Deviation from MODIS
The δAFAI product derived from the MODIS data on 8 August 2018 at 14:35 UTC allows a large area of Sargassum in the western part of Saint Lucia Island to be detected ( Figure 4). Note that the grey areas correspond to clouds, glint, or land. The relationship proposed by [10] allows δAFAI to be linked with the fraction of coverage "f" of Sargassum (Equation (6)). δAFAI(f) = K * f (6) where K = 4.41 × 10 −2 . A variation of δAFAI from 0 to 10 −3 means that the fraction of coverage varies between 0 and 2.3%. Such low fractions of coverage can be estimated due to high MODIS radiometric sensitivity.

AFAI Deviation from MODIS
The δAFAI product derived from the MODIS data on 8 August 2018 at 14:35 UTC allows a large area of Sargassum in the western part of Saint Lucia Island to be detected ( Figure 4). Note that the grey areas correspond to clouds, glint, or land. The relationship proposed by [10] allows δAFAI to be linked with the fraction of coverage "f" of Sargassum (Equation (6)).
where K = 4.41 × 10 −2 . A variation of δAFAI from 0 to 10 −3 means that the fraction of coverage varies between 0 and 2.3%. Such low fractions of coverage can be estimated due to high MODIS radiometric sensitivity.

NDVI Deviation from GOES Hourly Product
δNDVI calculated from GOES hourly product acquired on 8 August 2018 shows that the Sargassum aggregations detected with MODIS can be easily distinguished over the entire study area between 13:30 and 18:30 UTC, which corresponds to ca. 10:30-15:30 local solar time ( Figure 5). Before 13:30 and after 18:30 UTC, the Sargassum aggregation cannot be detected because NDVI is significantly altered by the sensor noise, as a result of the weak amplitude of the upward oceanic radiation at high sun zenith angles. A movie illustrating the time variation of NDVI over the course of the day is provided as a Supplementary Materials of this paper. Figures 4 and 5 show the significant potential benefit obtained from the synergy between sun-synchronous (MODIS) and geostationary sensors (GOES); MODIS is highly appropriate for the identification and the quantification of the Sargassum coverage while GOES is relevant to monitor the temporal dynamics of Sargassum.

NDVI Deviation from GOES Hourly Product
δNDVI calculated from GOES hourly product acquired on 8 August 2018 shows that the Sargassum aggregations detected with MODIS can be easily distinguished over the entire study area between 13:30 and 18:30 UTC, which corresponds to ca. 10:30-15:30 local solar time ( Figure 5). Before 13:30 and after 18:30 UTC, the Sargassum aggregation cannot be detected because NDVI is significantly altered by the sensor noise, as a result of the weak amplitude of the upward oceanic radiation at high sun zenith angles. A movie illustrating the time variation of NDVI over the course of the day is provided as a Supplementary Materials of this paper.

NDVI Deviation from GOES Hourly Product
δNDVI calculated from GOES hourly product acquired on 8 August 2018 shows that the Sargassum aggregations detected with MODIS can be easily distinguished over the entire study area between 13:30 and 18:30 UTC, which corresponds to ca. 10:30-15:30 local solar time ( Figure 5). Before 13:30 and after 18:30 UTC, the Sargassum aggregation cannot be detected because NDVI is significantly altered by the sensor noise, as a result of the weak amplitude of the upward oceanic radiation at high sun zenith angles. A movie illustrating the time variation of NDVI over the course of the day is provided as a Supplementary Materials of this paper.

GOES Noise Filtering
Although NDVI applied to GOES data is a relevant index to monitor aggregations, its radiometric sensitivity is not necessarily well adapted for ocean targets, and will induce some noise in the data as observed in Figure 5 (North-South and West-East striping). The first approach that was adopted here to reduce the noise in the data was to use the high revisit rate of GOES data to calculate an hourly product instead of using the 15-min product.
In addition to hourly averaging, another approach to decrease (or remove) the noise is to use a more sophisticated filtering technique. Many filtering techniques have been developed to overcome noise effects. One of the most efficient techniques is the "non local means" filter method (hereinafter noted as NLM) developed by Buades et al. [27]. Local mean filtering, by nature, would have affected the very fine structure we would like to preserve. Unlike "local mean" filters, which take the mean value of a group of pixels surrounding a target pixel to smooth the data, the NLM filtering technique considers a mean of all pixels over the full area, which are weighted by the similarity of these pixels in terms of neighborhood average with the pixel of interest. Typically, a high similarity leads to a high weight value and vice-versa. Such a filtering method leads to a much greater post-filtering clarity, while maintaining a suitable level of details in the image as compared to "local mean" algorithms [28]. The application of the NLM method shows that the noise is significantly reduced over the whole scene ( Figure 6). The NLM filtering allows the threshold value to be decreased to improve Sargassum aggregation detection.
its radiometric sensitivity is not necessarily well adapted for ocean targets, and will induce some noise in the data as observed in Figure 5 (North-South and West-East striping). The first approach that was adopted here to reduce the noise in the data was to use the high revisit rate of GOES data to calculate an hourly product instead of using the 15-min product.
In addition to hourly averaging, another approach to decrease (or remove) the noise is to use a more sophisticated filtering technique. Many filtering techniques have been developed to overcome noise effects. One of the most efficient techniques is the "non local means" filter method (hereinafter noted as NLM) developed by Buades et al. [27]. Local mean filtering, by nature, would have affected the very fine structure we would like to preserve. Unlike "local mean" filters, which take the mean value of a group of pixels surrounding a target pixel to smooth the data, the NLM filtering technique considers a mean of all pixels over the full area, which are weighted by the similarity of these pixels in terms of neighborhood average with the pixel of interest. Typically, a high similarity leads to a high weight value and vice-versa. Such a filtering method leads to a much greater post-filtering clarity, while maintaining a suitable level of details in the image as compared to "local mean" algorithms [28]. The application of the NLM method shows that the noise is significantly reduced over the whole scene ( Figure 6). The NLM filtering allows the threshold value to be decreased to improve Sargassum aggregation detection.

δNDVI Threshold Detection
The application of both the hourly averaging and NLM filtering makes the water and the Sargassum modes more readily distinguishable. A threshold value can be used for δNDVI to separate Sargassum from background. Sargassum aggregations are then identified and tracked from one hourly image to the next. The threshold value of 0.003 δNDVI was chosen by visualizing the histogram of the δNDVI. For the case where Sargassum are not present, the histogram of δNDVI consists only of one pronounced mode. In the case where Sargassum are present, the δNDVI distribution shows a second mode

δNDVI Threshold Detection
The application of both the hourly averaging and NLM filtering makes the water and the Sargassum modes more readily distinguishable. A threshold value can be used for δNDVI to separate Sargassum from background. Sargassum aggregations are then identified and tracked from one hourly image to the next. The threshold value of 0.003 δNDVI was chosen by visualizing the histogram of the δNDVI. For the case where Sargassum are not present, the histogram of δNDVI consists only of one pronounced mode. In the case where Sargassum are present, the δNDVI distribution shows a second mode overlapping the upper tail of the Sargassum-free ocean mode. The 0.003 threshold is optimally chosen to retain the non-ambiguous Sargassum optical signature (above 0.003) while rejecting most (99%) of the Sargassum-free ocean (below 0.003).

Performance of GOES Algorithm for Sargassum Identification
The comparison of MODIS δAFAI and GOES δNDVI over the same area (white box in Figure 6b) and at the same time (8 August 2018, 16:30 UTC) shows that GOES can detect Sargassum corresponding to maximum δAFAI above 8.10 −4 (Figure 7a) i.e., a fractional coverage ca. 2% (Equation (6)). A so-called confusion matrix was computed to compare the Sargassum detection using GOES (i.e., δNDVI > 0.003) with the Sargassum detection performed using MODIS (δAFAI > 1.79 × 10 −4 ) (Figure 7). Taking MODIS as the truth reference, the True Negative and True Positive rates of GOES detections are equal to 99% and 21% respectively and the False Positive and False Negative rates of GOES detections are equal to 1% and 79%, respectively. These results show that the NDVI is reliable in detecting the highest fractional coverage of Sargassum. The 2D scatterplot (Figure 7c) shows a quasi-linear relationship between δNDVI and δAFAI, over the δAFAI > 1.79 × 10 −4 range. This detection technique still satisfactorily retrieves the main patterns of Sargassum aggregations over the entire scene. This confirms the usefulness of GOES data to detect and track the largest Sargassum aggregations and study their dynamics. The δNDVI detection is only validated by comparison with the δAFAI product. In the current study, in situ measurements were not available, however in this area, the presence of Sargassum was confirmed by numerous beaching reports (http://sargassummonitoring.com/, accessed on 5 April 2021).
the True Negative and True Positive rates of GOES detections are equal to 99% and 21% respectively and the False Positive and False Negative rates of GOES detections are equal to 1% and 79%, respectively. These results show that the NDVI is reliable in detecting the highest fractional coverage of Sargassum. The 2D scatterplot (Figure 7c) shows a quasilinear relationship between δNDVI and δAFAI, over the δAFAI> 1.79·10 −4 range. This detection technique still satisfactorily retrieves the main patterns of Sargassum aggregations over the entire scene. This confirms the usefulness of GOES data to detect and track the largest Sargassum aggregations and study their dynamics. The δNDVI detection is only validated by comparison with the δAFAI product. In the current study, in situ measurements were not available, however in this area, the presence of Sargassum was confirmed by numerous beaching reports (http://sargassummonitoring.com/, accessed on 5 April 2021).

Impact of the Hourly GOES Product on the Noise Filtering and on the NDVI Value
The first approach adopted here was to compute an hourly product instead of using the 15-min product. The higher the sun zenith angle, the weaker the SNR. Such a technique decreases the level of noise in the data. When a given number of n measurement values are averaged, the noise decreases by √ and then the SNR is multiplied by √ . Here, because four images are averaged, the SNR is then multiplied by 2 (provided that the given pixel is cloud free for the four images). It also decreases the overall cloud cover within the area. However, this noise filtering approach has an impact on the Sargassum aggregation detection area. Since the drift velocity was estimated at 1 m.s -1 from GOES data, the drift distance over the course of one hour is about 3600 m, which corresponds to 7.2 pixels at a 500-m resolution. Thus, some pixels containing Sargassum and located on the aggregation edges may be averaged with Sargassum-free pixels. However, as GOES is used to estimate the velocity and direction of Sargassum aggregations, not to quantify their abundance, there is no downside to use such hourly averaged data.

Impact of the Hourly GOES Product on the Noise Filtering and on the δNDVI Value
The first approach adopted here was to compute an hourly product instead of using the 15-min product. The higher the sun zenith angle, the weaker the SNR. Such a technique decreases the level of noise in the data. When a given number of n measurement values are averaged, the noise decreases by √ n and then the SNR is multiplied by √ n. Here, because four images are averaged, the SNR is then multiplied by 2 (provided that the given pixel is cloud free for the four images). It also decreases the overall cloud cover within the area. However, this noise filtering approach has an impact on the Sargassum aggregation detection area. Since the drift velocity was estimated at 1 m·s −1 from GOES data, the drift distance over the course of one hour is about 3600 m, which corresponds to 7.2 pixels at a 500-m resolution. Thus, some pixels containing Sargassum and located on the aggregation edges may be averaged with Sargassum-free pixels. However, as GOES is used to estimate the velocity and direction of Sargassum aggregations, not to quantify their abundance, there is no downside to use such hourly averaged data.

Sargassum Aggregation Tracking
The high temporal resolution of GOES sensor is well adapted to track individual Sargassum aggregations, and then to investigate their drift speed and direction.
Hourly GOES Sargassum detections of 8 August 2018 were superimposed on the same plot in Figure 8. Each detection is represented with several dots scattered across the extent of the corresponding pixel. The number of dots is proportional to the δNDVI value (0 for δNDVI ≤ 0.003 to 15 for δNDVI ≥ 0.01), i.e., the dot density in a pixel corresponds to the Sargassum abundance. The color of the dots corresponds to the observation time, such that detections retrieved at different observation times in the same pixel can be distinguished. The aggregations retrieved at a given acquisition time are revealed by dot patterns of the same color and tracked in the figure.
plot in Figure 8. Each detection is represented with several dots scattered across the extent of the corresponding pixel. The number of dots is proportional to the NDVI value (0 for NDVI ≤0.003 to 15 for NDVI ≥0.01), i.e., the dot density in a pixel corresponds to the Sargassum abundance. The color of the dots corresponds to the observation time, such that detections retrieved at different observation times in the same pixel can be distinguished. The aggregations retrieved at a given acquisition time are revealed by dot patterns of the same color and tracked in the figure. Persistent and rather compact aggregations were first tracked. Their start position (13:30, blue) and end position (18:30, yellow) was determined visually for each selected aggregation. Then, the direction and speed value of about 20 aggregations were estimated with an accuracy that depends on the size of each aggregation. Accuracy was lower than or equal to 2 km in 5 h, which corresponds to about 0.1 m·s −1 . Although the method is manual, the main circulation pattern retrieved with that method, as illustrated by Figure  9. Persistent and rather compact aggregations were first tracked. Their start position (13:30, blue) and end position (18:30, yellow) was determined visually for each selected aggregation. Then, the direction and speed value of about 20 aggregations were estimated with an accuracy that depends on the size of each aggregation. Accuracy was lower than or equal to 2 km in 5 h, which corresponds to about 0.1 m·s −1 . Although the method is manual, the main circulation pattern retrieved with that method, as illustrated by Figure 9. The purple and red arrows show the direction and speed of the Sargassum aggregations as estimated from the GOES time series over the course of 5 h. The retrieved motion shows two main features. First, aggregations are subjected to one cyclonic eddy centered at 62.2° N-3.6° W (Figure 9a). Second, other aggregations move northwestwards ( Figure  9a). The day after (9 August 2018, Figure 9b), aggregations remain organized along an eddy on the western part, thus confirming the real occurrence of the eddy observed the day before. This aggregation appears to shift about 8.5 km to the southwest during one day, corresponding to a velocity of 0.1 m·s -1 . The eddy seems to strengthen with the diameter increasing by 75% to approximately 120 km. In the meantime, in the eastern part, the The purple and red arrows show the direction and speed of the Sargassum aggregations as estimated from the GOES time series over the course of 5 h. The retrieved motion shows two main features. First, aggregations are subjected to one cyclonic eddy centered at 62.2 • N-3.6 • W (Figure 9a). Second, other aggregations move northwestwards (Figure 9a). The day after (9 August 2018, Figure 9b), aggregations remain organized along an eddy on the western part, thus confirming the real occurrence of the eddy observed the day before. This aggregation appears to shift about 8.5 km to the southwest during one day, corresponding to a velocity of 0.1 m·s −1 . The eddy seems to strengthen with the diameter increasing by 75% to approximately 120 km. In the meantime, in the eastern part, the velocity of northwestward current also increased in magnitude (arrows).
Our analysis also points out the feasibility of tracking aggregations from one day to another by estimating the Sargassum drift distance. Table 1 provides the drift distance and velocity for four different Sargassum aggregations that were identified on 8 and 9 August 2018 at 16:30 UTC (Figure 10). Although a significant drift (white arrows) is observed, it remains possible to locate the aggregation. Table 1. Estimation of drift for the four Sargassum aggregations indicated in Figure 10.  Monitoring aggregations from one day to the next is possible by visual analysis. Such an analysis can help to improve forecasting of the short-term drift, hence the stranding location and time.

Consistency of Observed Drift with Satellite Derived Current and Temperature Data
Despite its coarse temporal (5 days) and spatial (about 40 km) resolutions, the OS-CAR product reveals the occurrence of an eddy (Figure 11) near the eddy that was identified using GOES observations (Figure 9a). A northwestward jet also appears in the western area (Figure 11a). The oceanic circulation observed in Figure 9 is then in agreement with the OSCAR-current speed data ( Figure 11). Finally, it is interesting to note that a temperature anomaly can be observed by GHRSST data (Figure 11b,c) in the specific area (3.5° N, 62.4° W) on 8 and 9 August 2018, which is consistent with the occurrence of an eddy there. Monitoring aggregations from one day to the next is possible by visual analysis. Such an analysis can help to improve forecasting of the short-term drift, hence the stranding location and time.

Consistency of Observed Drift with Satellite Derived Current and Temperature Data
Despite its coarse temporal (5 days) and spatial (about 40 km) resolutions, the OSCAR product reveals the occurrence of an eddy (Figure 11) near the eddy that was identified using GOES observations (Figure 9a). A northwestward jet also appears in the western area (Figure 11a). The oceanic circulation observed in Figure 9 is then in agreement with the OSCAR-current speed data ( Figure 11). Finally, it is interesting to note that a temperature anomaly can be observed by GHRSST data (Figure 11b

Pertinence of GOES Sargassum Product for Transport Models
Systematic comparison of observed Sargassum drift with the components acting on the drift (current, wind, and waves) will help to improve drift modeling and the subsequent drift forecast for a short-term period, and thus the stranding in location and time. There are several unknowns in the drift parameters, such as wind effect and Stokes contribution. Comparison of observed drift with the different components can help to better constrain these parameters.
Here, the comparison of the observed drift with modeled aggregations drift ( Figure  12) reveals some differences. The eddy observed with GOES is absent in the modeled aggregations drift velocity. The modeled advection jet towards the northwest direction is consistent with the Sargassum aggregations motion observation, but with a southwestward shift of hydrodynamics structures. The inflexion of the current direction is encountered further west in the model velocity drift. In addition, the modeled drift velocity is underestimated in comparison to GOES observed velocity; typically, the velocity ranges from 0.3 to 0.4 m·s −1 (61.4° N-4° W) for the model vs. 0.6 to 0.85 m·s −1 for GOES observations. Finally, even if these differences are small, their impact will grow

Pertinence of GOES Sargassum Product for Transport Models
Systematic comparison of observed Sargassum drift with the components acting on the drift (current, wind, and waves) will help to improve drift modeling and the subsequent drift forecast for a short-term period, and thus the stranding in location and time. There are several unknowns in the drift parameters, such as wind effect and Stokes contribution. Comparison of observed drift with the different components can help to better constrain these parameters.
Here, the comparison of the observed drift with modeled aggregations drift ( Figure 12) reveals some differences. The eddy observed with GOES is absent in the modeled aggregations drift velocity. The modeled advection jet towards the northwest direction is consistent with the Sargassum aggregations motion observation, but with a southwestward shift of hydrodynamics structures. The inflexion of the current direction is encountered further west in the model velocity drift. In addition, the modeled drift velocity is underestimated in comparison to GOES observed velocity; typically, the velocity ranges from 0.3 to 0.4 m·s −1 (61.4 • N-4 • W) for the model vs. 0.6 to 0.85 m·s −1 for GOES observations. Finally, even if these differences are small, their impact will grow when simulating the Sargassum transport. Consequently, these differences have to be reduced as much as possible, or at least, estimated.

Conclusions
In this paper, the synergy between a sun-synchronous (MODIS) and a geostationary (GOES) sensor for studying the spatio-temporal variation of Sargassum in the North Atlantic was investigated. Well-defined radiometric indices based on the reflectance at specific wavelengths were used for each sensor to derive Sargassum occurrence and abundance. While MODIS AFAI could be linked to the Sargassum fractional coverage and to the biomass, it suffered from heavy cloud masking and low revisit frequency. Current geostationary satellite such as GOES-16, with an hourly revisit, are well suited to identify and monitor Sargassum aggregations but with a weaker sensitivity. A "non local mean" filtering method was satisfactorily used to further reduce noise level. NDVI and AFAI were compared for the same acquisition time and showed a satisfactory agreement. Geostationary data allow individual Sargassum aggregations to be tracked over the course of several hours per day. The data were shown to be highly relevant for deriving their transport speed and direction. It was shown that Sargassum could drift between 40 and 60 km over the course of 19 h. The comparison of the remotely sensed drift with the drift velocity modeled from a global ocean circulation model showed a significant discrepancy. A systematic comparison of observed and simulated drift will help to better constrain drift components. This study demonstrated that the use of geostationary data for Sargassum detection and tracking is a promising approach to improve Sargassum drift modeling and forecast of Sargassum transport and stranding. The Sargassum detections from GOES and MODIS for the same acquisition time were compared. While GOES only retrieves 21% of the MODIS detections, the main Sargassum aggregations with high δAFAI values were retrieved with a satisfactory agreement. The forecast of Sargassum transport is of primary importance to support the need of stakeholders and end-users that make decisions in link with socio-economic interests for the Antilles Islands.  One way of improving transport model is to evaluate the wind effect in the modeling drift. The effect of wind was determined from large scale analysis [18,23] and non-local analysis. Moreover, here the wind is estimated from a large-scale model [18,23], which may not be accurate in the region. Finally, Stokes velocity has to be taken into account when modeling Sargassum drift. However, due to their density, Sargassum aggregations may dampen high frequency waves. Large-scale models of Stokes drift do not take into account this effect and thus the Stokes velocity contribution may be reduced by Sargassum. Therefore, wind and the Stokes effect on Sargassum drift are poorly known and these images may help to better investigate these effects.

Conclusions
In this paper, the synergy between a sun-synchronous (MODIS) and a geostationary (GOES) sensor for studying the spatio-temporal variation of Sargassum in the North Atlantic was investigated. Well-defined radiometric indices based on the reflectance at specific wavelengths were used for each sensor to derive Sargassum occurrence and abundance. While MODIS AFAI could be linked to the Sargassum fractional coverage and to the biomass, it suffered from heavy cloud masking and low revisit frequency. Current geostationary satellite such as GOES-16, with an hourly revisit, are well suited to identify and monitor Sargassum aggregations but with a weaker sensitivity. A "non local mean" filtering method was satisfactorily used to further reduce noise level. NDVI and AFAI were compared for the same acquisition time and showed a satisfactory agreement. Geostationary data allow individual Sargassum aggregations to be tracked over the course of several hours per day. The data were shown to be highly relevant for deriving their transport speed and direction. It was shown that Sargassum could drift between 40 and 60 km over the course of 19 h. The comparison of the remotely sensed drift with the drift velocity modeled from a global ocean circulation model showed a significant discrepancy. A systematic comparison of observed and simulated drift will help to better constrain drift components. This study demonstrated that the use of geostationary data for Sargassum detection and tracking is a promising approach to improve Sargassum drift modeling and forecast of Sargassum transport and stranding. The Sargassum detections from GOES and MODIS for the same acquisition time were compared. While GOES only retrieves 21% of the MODIS detections, the main Sargassum aggregations with high δAFAI values were retrieved with a satisfactory agreement. The forecast of Sargassum transport is of primary importance to support the need of stakeholders and end-users that make decisions in link with socioeconomic interests for the Antilles Islands.