remote sensing Toward Sentinel-2 High Resolution Remote Sensing of Suspended Particulate Matter in Very Turbid Waters: SPOT4 (Take5) Experiment in the Loire and Gironde Estuaries

: At the end of the SPOT4 mission, a four-month experiment was conducted in 2013 to acquire high spatial (20 m) and high temporal (5 days) resolution satellite data. In addition to the SPOT4 (Take5) dataset, we used several Landsat5, 7, 8 images to document the variations in suspended particulate matter (SPM) concentration in the turbid Gironde and Loire estuaries (France). Satellite-derived SPM concentration was validated using automated in situ turbidity measurements from two monitoring networks. The combination of a multi-temporal atmospheric correction method with a near-infrared to visible reflectance band ratio made it possible to quantify SPM surface concentration in moderately to extremely turbid waters (38–4320 g·m − 3 ), at an accuracy sufficient to detect the maximum turbidity zone (MTZ) in both estuaries. Such a multi-sensor approach can be applied to high spatial resolution satellite archives and to the new ESA Sentinel-2 mission. It offers a promising framework to study the response of estuarine ecosystems to global changes at unprecedented spatio-temporal resolution.


Introduction
Due to the Water Framework Directive and Marine Strategy Framework Directive, European state members are required to monitor the environmental status of their inland and coastal waters.Water quality monitoring is all the more needed in major estuaries hosting large urban areas and important industrial activities where the transfer of suspended material from terrestrial to oceanic ecosystems is significantly affected by land-use changes [1] and anthropogenic contamination [2,3].Amongst other environmental issues, suspended particulate matter (SPM) transport and accumulation have direct and indirect effects on aquatic ecosystems and human activities.Direct effects mainly concern the modification of river morphology due to erosion and sedimentation, locally resulting in bed sediment washout, mudflat and sandbar moves, or silting-up of harbours and navigation channels [4].Indirect effects include changes in light attenuation and primary production [5], oxygen saturation [6], alteration of fish nursery function [7], transport of SPM-attached pollutants [8], modification of nutrients pathways [9] and biogeochemical cycles [10].
In estuaries, SPM concentration is characterized by a high degree of variability due to the concomitant intrusion of marine waters and freshwater discharge.This results in intermingled physical and chemical processes such as sediments deposition and resuspension during tidal cycles [11], turbulence-induced particle accumulation [12] and flocculation [13].Measurement is needed at high temporal and high spatial resolution to accurately document SPM variation.Though the deployment of automated sensors makes it possible to study SPM temporal variation over a range of scales from hours to years [14], it is restricted to the geographical location of the instrumented stations in instrumented estuaries, and most estuaries are not instrumented [6].Due to its global coverage and ability to zoom in specific areas, satellite remote sensing provides synoptic observations of water reflectance in the visible and near infrared (NIR) spectral regions, which in turn can be used to map SPM surface concentration in turbid waters [15] at worldwide scale.
Landsat and SPOT high spatial resolution data have been widely used to detect SPM concentration in estuaries and nearshore areas.Since the early observations of the turbid discharge of the Mississippi waters in the Louisiana Bight [16], a variety of Landsat and SPOT images have been analysed to map SPM surface concentration in the world's main estuaries, including the Amazon [17], Gironde [18], Mekong [19], Rhone [20], and Yangtze [21,22] estuaries.Due to their low acquisition frequency, Landsat and SPOT observations are however temporally limited to document SPM variability in highly dynamical coastal areas.With a revisit time of 1-3 days, MODIS and MERIS data can be used more systematically to quantify surface SPM dynamics in large estuaries [23][24][25], but their spatial resolution (250 m -1 km) is not sufficient to detect small scale SPM variations in narrow estuaries and nearshore areas.
Launched in June 2015 by the European Space Agency (ESA), Sentinel-2 is the first satellite mission to systematically observe the world's coastal zone at both high spatial (10-60 m) and temporal resolution (5 days).With 13 spectral bands from the visible to the short wave infrared (SWIR) spectral regions, the Multispectral Instrument (MSI) on board Sentinel-2 will provide continuity and improvements to SPOT and Landsat image archives.To prepare the analysis of Sentinel-2 data, the French Space Agency (Centre National d'Etudes Spatiales, CNES, Paris, France) implemented the SPOT4 (Take5) experiment in 2013, providing a four-month time series of images simulating the revisit frequency, spatial resolution and large field of view of Sentinel-2 data.SPOT4 (Take5) data have been used for a variety of applications, from land and habitat studies to monitoring of inland and coastal waters.We present here an analysis of the SPOT4 (Take5) images acquired over two major estuaries of the French Atlantic coast, namely the Gironde and Loire estuaries.Satellite data were analysed in conjunction with automated in situ turbidity measurements, in order to evaluate the accuracy of the SPM algorithms previously proposed for both study sites [26].The Doxaran et al. (2003) algorithm has previously been validated in the Gironde estuary at moderate resolution using MODIS data [23], and already applied to SPOT and Landsat images [18], but a thorough evaluation of the approach to high resolution satellite images has never been done, neither in the Gironde, nor in the Loire estuary.In the present study, we used the opportunity of the SPOT4 (Take5) experiment to evaluate SPOT4 and Landsat SPM data processed using the Doxaran et al. (2003) algorithms [26], and to study the spatial distribution of SPM surface concentration at high resolution in both estuaries.

Gironde Estuary
The Gironde estuary, formed by the confluence of the Garonne and Dordogne rivers (Figure 1), has the largest surface area of any estuary in Europe.Its length is about 75 km, and its width ranges from 3-11 km.The Garonne and Dordogne rivers exhibit watershed areas of 57,000 km 2 and 24,000 km 2 , respectively.The Gironde's flow rate averages 1100 m 3 •s −1 [27].The Garonne and the Dordogne respectively contribute 65% and 35% of freshwater inputs in the Gironde estuary.The freshwater discharge in the Garonne ranges from less than 100 m 3 •s −1 to more than 4000 m 3 •s −l , whereas it varies between 200 and 1500 m 3 •s −1 in the Dordogne [28].The Gironde's morphology is typical of wave and tide dominated estuaries.Due to tidal amplitude ranging from 2-5 m, the Gironde estuary is classified as macrotidal.It presents a well-developed maximum turbidity zone (MTZ).From a remote sensing point of view, the MTZ is defined here as the zone where suspended particulate matter (SPM) surface concentration exceeds 1000 g•m −3 .A MTZ has been both observed in the downstream and upstream sectors of the Gironde estuary.During summer, it can move upstream in the Garonne and reside around Bordeaux, about 30 km upstream the Dordogne and Garonne convergence [29,30].

Loire Estuary
With a length of 1012 km and a watershed area of 117,000 km 2 , the Loire is the largest river in France.The Loire estuary is about 100 km long, stretching from a freshwater portion upstream of Nantes to the outer estuary downstream of Saint-Nazaire (Figure 2).The Loire's flow rate averages 850 m 3 •s −1 .It varies from less than 300 m 3 •s −l during summer droughts to more than 4000 m 3 •s −l during winter floods.Estuarine water residence time varies from less than one day during floods to 10 days during droughts [31].The Loire estuary is macrotidal, with 4 m average tidal amplitude at Saint-Nazaire.The intertidal area covers more than 30 km 2 , most of which are mudflats colonized by benthic microalgae [32].
As in many other macrotidal estuaries, the Loire estuary is characterized by large variations in SPM concentration, which varies from about 50 g•m −3 in the moderately turbid freshwater sector to a zone of maximum turbidity where SPM surface concentration exceeds 1000 g•m −3 [9].In the Loire estuary, the MTZ is 20-50 km long, depending on river flow and tidal current [33].In the MTZ, SPM surface concentration typically reaches 2000 g•m −3 during spring tide, but decreases to 100 g•m −3 during neap tides, when suspended particles are temporarily trapped just above sediment bed, creating fluid mud layers where bottom SPM concentration can exceed 20 kg•m −3 [34].The total mass of SPM transported from the river to the estuary is about one million tons every year [2].

SPOT4 (Take5) Experiment in the Gironde and Loire Estuaries
At the end of SPOT4's mission, before de-orbiting the satellite in June 2013, the Centre d'Etudes Spatiales de la BIOsphere (CESBIO) proposed a four-month experiment to increase the revisit time of SPOT4 to 5 days, thus potentially serving as a simulator for the Sentinel-2 mission (Table 1).The experiment, which was named SPOT4 (Take5) in reference to Paul Desmond and Dave Brubeck's famous jazz piece (Take Five), was successfully implemented by the CNES.The worldwide remote sensing community was invited to join the experiment and to propose observation sites for which high temporal resolution was relevant to environmental monitoring or methodology development purposes.We proposed to study two sites centered on the Gironde and Loire estuaries.For the Gironde study, the SPOT4 (Take5) spatial coverage was restricted to the southern part of the estuary, focusing on the Garonne river (Figure 1), which is a narrow (400 m wide in Bordeaux) meandering and macrotidal river.The SPOT4 (Take5) image extent does not cover the whole estuary, and is susceptible to match the records of only one automated in situ turbidity station, as described in Section 2.3.A dataset of Landsat images covering the whole Gironde estuary was therefore used to complement SPOT4 (Take5) data.Landsat images were acquired from January 2011 to December 2013.For the Loire estuary, the SPOT4 (Take5) footprint covers the whole estuary, and a total of four automated in situ turbidity stations were available for match-up purpose (Figure 2).SPOT4 (Take5) and Landsat data were processed and distributed by the Theia Land Data Centre, a French national inter-agency organization designed to foster the use of images issued from the space observation of land surface and nearshore coastal areas.Top-of-atmosphere radiance data were corrected from aerosol effects using the multi-sensor atmospheric correction and cloud screening (MACCS) method [35][36][37].Level 2A data include ortho-rectified surface reflectance, along with a mask for clouds and their shadows, and a correction of adjacency effects.In situ regional algorithms developed by Doxaran et al. (2003 [26]) to retrieve SPM surface concentration from reflectance measurements in turbid waters were then applied to SPOT4 data in the Gironde (Equation ( 1)) and Loire (Equation ( 2)) estuary, and to Landsat data in the Gironde estuary (Equation ( 3)): where Rrs is the remote sensing surface reflectance, XS1 and XS3 are SPOT4 spectral bands in the green (0.5-0.59 nm), and near infrared (XS3: 0.79-0.89nm), and LG and LNIR are Landsat spectral bands corresponding to 0.53-0.59nm, and 0.85-0.88nm, respectively.Doxaran et al. (2003 [26]) demonstrated that relations based on reflectance ratios between near infrared and visible spectral bands were robust to variations in sediment type and to changes in illumination conditions.The robustness of the approach was previously evaluated in the Gironde estuary using MODIS Aqua data.Based on 75 match-ups, the relative uncertainty was found to be around 20% over a range of SPM concentration from 10-2250 g•m −3 [23].Though the algorithm has been already applied to SPOT and Landsat data in the Gironde estuary [18], a thorough validation of its performance with high resolution satellite data is still lacking in the Gironde and Loire estuaries.In the present study, a dataset of automated in situ measurements was used to assess the accuracy of satellite-derived SPM concentration in both estuaries.

In Situ Turbidity Measurements in the Gironde Estuary
The Marel Gironde estuary network (MAGEST), operating since 2005, is composed of four MAREL stations (Mesures Automatisées en Réseau pour l'Environnement et le Littoral) [38] located in the central and fluvial estuary.Due to the spatial coverage of the SPOT4 (Take5) experiment in the Gironde estuary, only one MAGEST station was used in the present study.It is located in the fluvial section of the Gironde to the north of the city of Bordeaux (Figure 1).The water was sampled at less than 1 m below surface, and pumped to the sensors.Turbidity measurements, given in Nephelometric Turbidity Unit (NTU), were automatically performed using commercial sensors (Turbimax CUS31, Endress & Hauser, Reinach, Switzerland) covering the range 0-9999 NTU range.Data were qualitycontrolled and validated following a protocol adapted to the estuarine environment [6], including a regular calibration of the turbidity sensors.Turbidity measurements were used to estimate SPM concentration using a regional empirical relationships derived from laboratory calibration.The most recent turbidity vs. SPM concentration relationship [23], which was obtained for a range of SPM concentration up to 2000 g•m −3 , was used for this study (Equation ( 4)): A recent analysis of the 2005-2013 turbidity time series in the fluvial section showed seasonal variations from <100 NTU during droughts to >9000 NTU during floods [39].These changes can be related to MTZ migration from the Gironde estuary to the Garonne River [40].In addition to the impact of freshwater discharge on the MTZ location, surface turbidity varies during the tidal cycle: first, turbidity shows typical patterns of erosion-deposition cycle observed in macrotidal estuaries [41][42][43], and second, two distinct SPM peaks are observed during the semi-diurnal tidal cycle at mid-flood and mid-ebb due to resuspension by tidal currents [6].

In Situ Turbidity Measurements in the Loire Estuary
Along the Loire estuary, continuous turbidity measurements have been performed since 2007 in the frame of the SYVEL (Système de veille dans l'estuaire de la Loire) monitoring network operated by the GIPLE (Groupement d'Intérêt Public Loire Estuaire, Nantes, France) using MAREL stations (Figure 2).Sensors are housed inside an instrumented chamber fixed on a pier.As in the Gironde estuary, the Loire water was sampled at less than 1 m below surface, pumped to the sensors, and turbidity measurement was automatically performed every 30 min using a Turbimax CUS31 sensor (Endress & Hauser, Reinach, Switzerland) covering the range 0-9999 NTU range.Turbidity data were calibrated in SPM concentration using a regional relationship [44].

Evaluation of SPM Algorithm Performance
High-frequency in situ SPM time series from the MAGEST and SYVEL were used for match-ups with satellite-derived SPM concentration in the Gironde and Loire estuaries.For this exercise, in situ data were interpolated to the exact time of image acquisition in order to prevent from uncertainties associated with small scale temporal variation.For both study sites, match-ups were performed using four pixels, including the pixel nearest to the station's location and a longitudinal line of three pixels located 30 m away from the riverbank, in order to limit potential adjacency effects.The accuracy of satellite-derived concentration was evaluated using standard metrics, including the relative root mean square error (RRMSE) and the mean normalized bias (MNB) [45].All statistics were computed using R [46].

Satellite Data Summary
During the SPOT4 (Take5) experiment, satellite acquisitions were successfully performed every 5 days from February to early June 2013 in the Gironde and Loire estuaries.Due to poor weather conditions however, only a limited number of (partially) cloud free data was available.A total of nine SPOT4 images were selected within the four months of the Take5 experiment (Tables 2 and 3).
Selected images were acquired over a wide range of hydrological conditions.River flow varied from 581 to 1650 m 3 •s −1 in the Gironde estuary, and from 1400-3200 m 3 •s −1 in the Loire estuary.A total of 11 Landsat images acquired in 2011 and 2013 were selected and added to the Gironde satellite dataset.The addition of Landsat imagery in 2011 made it possible to extend the lower limit of hydrological conditions experienced during satellite overpasses, with river flow smaller than 100 m 3 •s −1 on several occasions.In the Gironde estuary, the addition of SPOT4 (Take5) and Landsat8 data in 2013 made it possible to study SPM concentration changes on a monthly basis from February to October 2013.For the selected satellite imagery, average tidal range respectively varied between 3.30 and 5.15 m in the Gironde estuary, and between 2.24 and 5.55 m in the Loire estuary.In the Loire estuary, SPOT4 images were acquired during both neap and spring tides, whereas Gironde's SPOT4 and Landsat cloud free images were mostly acquired during neap tides, thus probably leading to an underestimated picture of SPM concentration in the fluvial section of the estuary.The high river flows observed during the SPOT4 (Take5) experiment further contributed to flush SPM seaward and decrease SPM surface concentration in the Garonne section of the Gironde estuary.

Range of SPM Concentration during SPOT4 (Take5) Experiment
The Gironde and Loire automated turbidity sensors recorded a very wide range of SPM surface concentration during the SPOT4 (Take5) experiment, from 8.95 to 5440 g•m −3 (Figure 3).At Bordeaux, minimum SPM concentration was observed in the Gironde estuary on 19 April 2013.At this station, SPM concentrations did not exceed 859.3 g•m −3 (Figure 3a).The range of concentration recorded that year was relatively low in comparison with interannual averages [6].This is probably caused by the downstream shift of the MTZ resulting from fresh water discharges associated with important precipitations that occurred at the beginning of the year.Winter and spring 2013 were indeed very rainy, and it is likely that important land runoff and river flow contribute to an overall decrease of SPM concentration in the Garonne section of the Gironde estuary.Moreover, the location of the Bordeaux station did not allow recording the upstream intrusion of turbid waters during spring tides.  2 and 3. SPM surface concentration was approximately five times higher in the Loire estuary than in the Garonne estuary.Automated SPM time series at Donges consistently displayed a high degree of variation (Figure 3b).SPM concentration exceeded 5000 g•m −3 several times in winter 2013.Maximum SPM concentration (5.4 kg•m −3 ) was recorded on 12 March 2013 during spring tide (tidal range >6 m).The influence of tidal dynamics on turbidity time-series was clearly visible.The lunar neap-spring cycle significantly imprinted SPM variations with a regular fluctuation of about 15 days, and maximum turbid events always occurred during spring tides.At shorter time-scales, the quarter-diurnal cycle lead to dramatic changes in SPM concentration in terms of amplitude (>5000 g•m −3 ) and duration (<3 h), resulting from water height changes as well as from the alternation of sedimentation and resuspension processes caused by tidal dynamics.

SPM Match-Ups in the Gironde and Loire Estuaries
A total number of 30 data points was available for SPM match-up purpose (Figure 4).In the Gironde estuary, the Bordeaux turbidity station provided 15 match-ups.In the Loire estuary, though only five images were selected, a total of 15 pairs of concomitant in situ and satellite measurements were extracted at four in situ turbidity stations.The match-up dataset covered a very wide range of SPM concentration, from 37.9-4320 g•m −3 .The RRMSE and MNB were 220.0% and 89.3%, respectively.There was a significant linear relationship between in situ and satellite-derived SPM concentration (r-squared regression coefficient is 0.75, p-value < 10 −9 ), with a slope of 0.66.The accuracy of satellite SPM concentration retrieval was better in the 100-1000 g m −3 interval, with a RRMSE and a MNB of respectively 49.5% and 16.3%, and a slope very close to 1 (Figure 4, solid line).

SPM Spatial Distribution in the Gironde and Loire Estuary
A focus on a few SPOT4 and Landsat8 images spanning different hydrological and tidal regimes allowed us to better appraises the performance of the SPM algorithm for both satellite sensors.Examples of satellite-derived SPM surface concentration maps are displayed in Figures 5 and 6 for the Gironde and Loire estuaries, respectively.

SPM Concentration Maps in the Gironde Estuary
Around Bordeaux, in the Garonne section of the Gironde estuary, SPOT4 (Take5) data generally displayed low SPM concentrations in February 2013 (Figure 5a).Due to strong precipitation in winter and spring, enhanced land runoff increased water height and vertically diluted SPM concentrations in the estuary's fluvial section.The SPOT4 (Take5) time series of SPM concentration did not exhibit strong spatio-temporal variations, and turbidity levels were quite similar from February to May 2013.In June, a patch of turbid water was observed to the north of Bordeaux (Figure 5b), with SPM surface concentration 5 to 30 times higher than in February.This increase in SPM concentration can probably be explained by bed sediment resuspension resulting from spring tide currents and increase in river flow.
The spatial frame of the SPOT4 (Take5) images in the Gironde estuary did not allow us to observe SPM spatial distribution outside the Garonne section.The addition of Landsat data made it possible to observe the rest of the Gironde estuary, from the Garonne and Dordogne confluence to the Gironde's mouth (Figure 5c,d).In the July 2013 example, the Landsat8 image confirmed the confinement of the relatively clear waters within the upstream sections of the estuary, whereas the most turbid waters stretched from the Garonne and Dordogne confluence down to the sea (Figure 5c).In October 2013, the situation was reversed: maximum SPM concentration was found in the fluvial section, whilst surface waters were generally less turbid downstream Pauillac (Figure 5d).Longitudinal SPM transects displayed reversed gradients on 10 July and 30 October (Figure 7).In addition to the overall longitudinal changes, SPM surface concentration showed variations associated with varying sedimentary patterns.For example, in the central section of the estuary, the July 2013 example exhibited low-high turbidity strips parallel to the orientation of surface currents, whereas these patterns were far less pronounced in October 2013.In July, the strips of high turbidity were presumably linked to increased resuspension over shallow mud flats.Tidal conditions were different on 10 July (mid-ebb during spring tide) and 30 October (low tide during neap tide), thus leading to distinct tidal currents' strength and associated resuspension processes [6].The July image may also illustrate the impact of spring tides on mud bank erosion.More data are necessary to better appreciate the respective role of tidal processes on sediment dynamics at the scale of the estuary.

SPM Concentration Maps in the Loire Estuary
Two SPOT4 (Take5) images and corresponding longitudinal transects are shown to illustrate the diversity of SPM spatial distribution in the Loire estuary (Figures 6 and 8).The first example corresponds to an image acquired on 21 February 2013, three hours after low tide during neap tide (tidal range was 2.24 m, Table 3).SPM surface concentration longitudinally decreased from Nantes to the ocean off Saint-Nazaire (Figure 6a).SPM concentration did not exceed 130 g•m −3 , even in the most turbid sectors of the estuary around 1°50′W (Figure 8).SPM concentration remained above 100 g•m −3 over a 20 km section upstream of 2°05′W.In this section, SPM was homogeneously distributed between the river banks.In the outer estuary, SPM concentration was generally lower than 35 g•m −3 , but varied importantly from the northern to the southern side of the river, with maximum values up to 250 g•m −3 near shore at about 47°12′N, 2°11′W.The turbid plume spreading out the Loire was confined to the south-eastern section of the estuary.Off the estuary, the marine waters were moderately turbid, probably due to wave-induced particle resuspension.Wave height was higher than 1 m, and wind speed averaged 37 km•h −1 , with gusts at 59 km•h −1 .
The second example corresponds to an image acquired one month later on 28 March 2013, 44 min before low tide during a spring tide (tidal range is 5.55 m, Table 3).Wave height was smaller than 0.5 m, and wind speed did not exceed 30 km•h −1 .The distribution of SPM concentration was dramatically different than on 21 February 2013.In contrast with the previous situation, the most turbid waters were located downstream to Nantes (Figure 6b).Relatively low turbid waters, with SPM concentration as low as 5 g•m −3 , were observed in the fluvial section.Within a few km from Nantes, SPM concentration exhibited very large longitudinal changes, from less than 10 g•m −3 to more than 1 kg•m −3 (Figure 8).The most noticeable feature consisted in the rapid formation of a maximum turbidity zone between 2°10′ and 2°0′W.The MTZ occupied a length of 3.5 km.Downstream of Saint-Nazaire, a turbid plume spreading out of the Loire was located in the center of the estuary.Within this turbid plume, across-side variations in SPM concentration were very important at small spatial scale.They were mostly associated with bathymetry changes, and the two navigational channels were visible along the two sides of the outer Loire estuary between 2°10′ and 2°15′W (Figure 9).A patch of extremely turbid waters (SPM concentration higher than 2000 g•m −3 ) was observed at 2°50′W, probably due to the concentration of SPM in very shallow waters and to enhanced sediments erosion and resuspension by nearby mudflats.

Accuracy of Satellite-derived SPM Concentration Maps
It is generally acknowledged that the accuracy of satellite match-ups is affected not only by uncertainties in bio-optical algorithms, but also by errors in field measurements, spatial and temporal differences between sampling and satellite overflight, and atmospheric correction uncertainties.

In Situ Measurements Uncertainties
Variability in particle size, composition and inherent optical properties affects the turbidity vs. SPM concentration relationship [47].Regional calibration of turbidity measurements into SPM concentration are therefore required [48].In the Gironde estuary, the relationship between SPM concentration and turbidity was linear from 0-2000 g•m −3 [23].A polynomial fit was used to calibrate turbidity sensors in the Loire estuary.Beside the calibration of turbidity data, another source of error in match-up exercises is generally associated with temporal differences between water sampling and satellite overflight.Due to the high-frequency of in situ turbidity measurements acquired automatically in the frame of the MAGEST and SYVEL networks, it was possible in the present study to interpolate SPM time-series to the exact time of satellite images acquisition, thus preventing temporal biases.In highly dynamic estuarine and nearshore waters, SPM concentration can change dramatically within less than one hour due to tidal processes, and the use of high-frequency data is strongly recommended for satellite match-ups [49].

Spatial Heterogeneity and Match-Up Pixel Selection
An in situ experiment was carried out in the Gironde estuary at the Pauillac turbidity station to assess the impact of small-scale spatial heterogeneity on satellite turbidity retrieval over pixels adjacent to the riverbank.Water sampling was simultaneous performed just along the Pauillac embankment, and 50 m away from it during two consecutive days in order to encompass ebb and flood tidal conditions.The difference between near-and offshore SPM concentrations varied with tidal currents.A significant increase in SPM concentration was observed shoreward, with nearshore values up to two times higher than in the center of the estuary.This preliminary result highlights the importance of small-scale SPM spatial variation in nearshore waters.Small-scale spatial variability prevents the use of pixels coarser than 50 m for nearshore estuarine studies.It has also to be taken into account in match-ups exercises.
As the SPM match-up was performed using in situ data acquired from turbidity stations sampling water in the immediate vicinity of the riverbank (Figures 1 and 2), a preliminary sensitivity analysis was performed in the selection of SPOT4 match-up pixels.Several configurations in the number and location of match-up pixels were tested.The best configuration was obtained when four pixels were used, including the pixel nearest to the sampling point and three additional pixels located along a line parallel to the flow, at about 30 m from the riverbank.In the Loire estuary, due to small-scale spatial heterogeneity, the match-up accuracy decreased when the number selected pixels was higher than 5, and when the pixels selected were too far away (i.e., >40 m) from the station location.Pixels too close to the shore were removed, in order to take into account the possible contamination of marine pixels by land adjacency effect [50].In the present study, the size of SPOT4 (20 m) and Landsat (30 m) pixels was small enough to select the best configuration of match-up pixels, and it is expected that the use of higher resolution satellite data such as Sentinel-2 (10-60 m), SPOT5 (10 m), SPOT6 (6 m) or Pléiades (2 m) might improve pixel selection and match-up accuracy.Further studies are needed to implement a standard match-up protocol for the validation of remote-sensing data using land-based in situ turbidity station with nearshore water sampling.The Doxaran et al. (2003) algorithm is based on empirical relationships between SPM concentration and a NIR/VIS reflectance ratio established using a large dataset of in situ measurements [26].Two distinct relationships were developed for the retrieval of SPM concentration from SPOT data in the Gironde and Loire estuaries, and a third relationship was also proposed to derive SPM concentration from c (Figures 7 and 8 in [26]), and when used together, no bias was observed in the multi-sensor and multi-site match-up performed in the present study (Figure 4).Though the match-up was satisfactory over a very wide range of SPM concentration (38-4320 g•m −3 ), its performance was better in the 100-1000 g•m −3 range.
The Doxaran et al. (2003) relationships displayed some scatter at the lower range (see Figure 9 in [26]), and the use of an exponential fit could amplify the uncertainties.For turbid waters in the 10-500 g•m −3 range, it has been recently demonstrated that a linear relationship between SPM concentration and a NIR/VIS reflectance band ratio was sufficient to accurately determine SPM concentration in the Loire estuary and regional waters [51].In extremely turbid waters (SPM>1000 g•m −3 ), it is generally recognized that NIR/VIS Rrs ratio may saturate [26], and more data are still needed to develop robust SPM algorithms for highly turbid waters [48].Beside the SPM algorithm itself, other factors such as small-scale spatial heterogeneity (see Section 4.1.2) and atmospheric correction uncertainties (see below) also contribute to degrade match-up performance in turbid waters.

Atmospheric Correction
The reflectance at the top of the atmosphere (TOA) is the sum of the atmospheric reflectance and of the surface marine reflectance transmitted by the atmosphere.Over clear oceanic waters, the contribution of the atmosphere is computed in the NIR where the marine reflectance is assumed to be zero because of the strong pure water absorption and of the negligible contribution by SPM [52].In turbid waters, the NIR dark pixel assumption is not valid anymore [53].The contribution of SPM to the marine reflectance in the NIR, and for highly turbid waters in the SWIR between 1000 and 1250 nm [51,54], has to be taken into account to separate the atmospheric effects from the TOA reflectance.Beyond 1600 nm, it can be assumed that the SWIR marine reflectance is negligible, and a SWIR-based atmospheric correction has been recently implemented for Landsat8-OLI data over turbid waters [55].As this latter method has not been implemented for SPOT4 data, another approach was used in the present study for SPOT4 atmospheric correction.
In the present study, the MACCS atmospheric correction was used for cloud detection and atmospheric correction of satellite data [35].This method has been developed for the correction of multi-temporal images acquired using a constant viewing angle.Over land, it assumes that the surface reflectance in the blue does not change quickly with time.Changes in TOA reflectance will be detected and identified as clouds or aerosol optical thickness variations, making it possible to characterize the atmospheric reflectance, and to remove it from the TOA signal to determine surface reflectance [35].A notable advantage of the MACCS method is its ability to be applied to different sensors, as long as a sufficient number of multi-temporal data is available, and that the viewing angle remains constant for each sensor [36,37].These two conditions are met for SPOT4 (Take5) and Landsat8 imagery, and the MACCS atmospheric correction was applied to all images selected in Tables 2 and 3. A thorough evaluation of the MACCS atmospheric correction for the Gironde and Loire waters was out of scope of the present study, but as the SPM match-up did not show any bias between sites and sensors (Figure 4), it seemed that MACCS can be seamlessly applied to SPOT4, Landsat5, 7 and 8 data to study SPM concentration in turbid waters (SPM concentration > 100 g•m -3 ).

Influence of Environment on SPM Dynamics
We presented here a few examples of satellite-derived SPM distribution, and the number of SPM maps was not sufficient for a detailed analysis of SPM dynamics in response to environmental changes.The interpretation of SPM spatial distribution was mainly based on analyzing the time of satellite acquisition in the context of the tidal cycle, which is generally recognized as one of the main factors affecting SPM dynamics [6,29].River topography (i.e., mud banks location), morphology and bathymetry can also affect SPM spatial distribution [39].The high spatial resolution of SPOT4 made it possible to show the influence of river topography on SPM small-scale distribution.The tidal river is also very sensitive to river discharge variations, and it has been recently shown that the MTZ was expulsed from the Bordeaux sector when the river flow exceeded 350 m 3 •s −1 [39], which was generally consistent with the situations described here in the Gironde estuary.In the Loire estuary, winter and spring 2013 were exceptionally rainy, and the high rain falls affected the river discharge.From February to May 2013, the daily river flow was systematically higher than the decadal average (850 m 3 •s −1 ).The mean SPOT4 (Take5) river flow was 1711 m 3 •s −1 , thus exceeding the interannual average by a factor of 2. The series of SPM data presented here is however too short to analyze the effect of above-average hydrological conditions on SPM dynamics and MTZ characteristics.

Conclusions and Future Directions
In this paper, we validated the use of SPOT4 (Take5) and Landsat data to detect and quantify SPM surface concentration in the Gironde and Loire estuaries.We demonstrated that the combination of the MACCS atmospheric correction [35] with a simple NIR/VIS reflectance band ratio algorithm [26] satisfactorily quantified SPM concentration in moderately to extremely turbid waters (38-4320 g•m −3 ) at an accuracy sufficient to detect the maximum turbidity zone in two optically distinct estuaries.We presented a multi-sensor approach that can be applied to SPOT and Landsat archives, as well as to other high spatial resolution Earth Observation sensors such as the Sentinel-2 mission by the European Space Agency.
Using moderate resolution remote sensing, it has been previously demonstrated that decadal satellite time series can be used to study the turbidity response of large estuaries to global changes [22,23,25,56].
Based on a few high spatial resolution SPOT and Landsat images, we provided here a preliminary description of the fine-scale spatial distribution of SPM concentration in the Gironde and Loire estuaries during contrasted hydrological conditions.At the time we wrote the paper, the SPOT5 (Take5) experiment (April-August 2015) was currently under way, and the Sentinel-2 satellite has been successfully launched on 23 June 2015.It is expected that in a near future, the high spatial resolution (10 m) and high frequency (5 days) of the data provided by these two missions will make it possible to more systematically document SPM concentration changes in nearshore and estuarine areas, and to more accurately study the response of coastal ecosystems to environmental changes.

Figure 1 .
Figure 1.Map of the Gironde estuary showing SPOT4 (Take5) spatial coverage and the location of the MAGEST automated turbidity station used for match-up.

Figure 2 .
Figure 2. Map of the Loire estuary showing SPOT4 (Take5) spatial coverage and the location of the automated turbidity stations of the SYVEL network.

Figure 3 .
Figure 3. (a) Suspended particulate matter (SPM) concentration measured using automated turbidity station in the Gironde estuary at Bordeaux; (b) SPM concentration measured using automated turbidity station in the Loire estuary at Donges.Vertical arrows show the dates of partially cloud free SPOT4 (Take5) images listed in Tables2 and 3.

Figure 4 .
Figure 4. Suspended particulate matter (SPM) concentration match-up between in situ measurements and satellite data, respectively obtained using SPOT4 (black circles) and Landsat (black triangles) in the Gironde estuary, and SPOT4 (white circles) in the Loire estuary.Vertical lines show the standard deviation of the four pixels selected for match-up.The dashed line shows the linear fit from 37.9-4320 g•m −3 , whereas the solid line shows the linear fit between 100 and 1000 g•m −3 .

Figure 6 .
Figure 6.(a) SPOT4 suspended particulate matter (SPM) concentration in the Loire estuary during neap tide on 21 February 2013; (b) SPOT4 SPM concentration during neap tide on 28 March 2013.The crosses show the location of transect pixels used in Figure 8.

Table 2 .
List of selected satellite data for the Gironde estuary.Water height and tidal range provided by the SHOM at Bordeaux.River flow data measured at Tonneins, 100 km upstream of Bordeaux.

Table 3 .
List of selected SPOT4 (Take5) data for the Loire estuary.Water height and tidal range provided by the SHOM for the main harbour in the Loire estuary (Saint-Nazaire).River flow data measured at Montjean-sur-Loire, 60 km upstream of Nantes.