Spatial variability of surface suspended particulate matter (SPM) concentration in the Gironde estuary and their relationships with environmental forcing are investigated through high spatial resolution multispectral data collected from July 2013 to August 2018 by the Operational Land Imager (Landsat-8/OLI) and MultiSpectral Instrument (Sentinel-2/MSI). A principal component analysis using the T-mode orientation is applied to the discontinuous multiannual time series composed of 41 remotely sensed images. The three first principal components (PC1, PC2, and PC3) explain 65.7% of the total variance. The SPM distribution associated with PC1 and PC2 exhibits a privileged along-estuary direction of the oscillation modes, while the spatial patterns of PC3 are clearly dominated by lateral oscillations opposing channels and shoals. The main environmental factors affecting the SPM distribution are identified by the analysis of their temporal patterns. The tidal range and the daily river discharge control the spatial patterns of PC1 and PC2, while the tidal cycles and the wind speed are significantly correlated with PC3. Furthermore, the analysis in the along-estuary and lateral directions shows marked longitudinal and transverse SPM gradients and a strong control of bathymetry on the SPM spatial distribution. For the first time, we highlight that the maxima of surface residual turbidity are located on the shoals regardless of the environmental (tidal and hydrological) conditions. Compared to previous studies, usually based on single-point in-water column measurements, these results provide a novel and complementary description of the spatial variability of SPM. They are useful to validate sediment transport numerical models, but also may improve our understanding of suspended sediment dynamics in estuarine systems governed by an estuarine turbidity maximum.Read less <