Analysis of the Transport of Aerosols over the North Tropical Atlantic Ocean Using Time Series of POLDER/PARASOL Satellite Data

: The time series of total, ﬁne and coarse POLAC/PARASOL aerosol optical depth (AOD) satellite products (2005–2013) processed by the POLAC algorithm are examined to investigate the transport of aerosols over the North Tropical Atlantic Ocean, a region that is characterized by signiﬁcant dust aerosols events. First, the comparison of satellite observations with ground-based measurements acquired by AERONET ground-based measurements shows a satisfactory consistency for both total AOD and coarse mode AOD (i.e., correlation coefﬁcients of 0.75 and bias ranging from − 0.03 to 0.03), thus conﬁrming the robustness and performance of POLAC/PARASOL data to investigate the spatio-temporal variability of the aerosols over the study area. Regarding ﬁne mode aerosol, POLAC/PARASOL data present a lower performance with correlation coefﬁcient ranging from 0.37 to 0.73. Second, the analysis of POLAC/PARASOL aerosol climatology reveals a high contribution of the coarse mode of aerosols (AOD c between 0.1 and 0.4) at long distance from the African sources, conﬁrming previous studies related to dust transport. The POLAC/PARASOL data were also compared with aerosol data obtained over the North Tropical Atlantic Ocean from MACC and MERRA-2 reanalyses. It is observed that the total AOD is underestimated in both reanalysis with a negative bias reaching − 0.2. In summary, our results thus suggest that satellite POLAC/PARASOL observations of ﬁne and coarse modes of aerosols could provide additional constraints useful to improve the quantiﬁcation of the dust direct radiative forcing on a regional scale but also the biogeochemical processes such as nutrient supply to the surface waters.


Introduction
Composed of fine solid particles or liquid droplets suspended in air, aerosols constitute an important part of the atmosphere. Through their microphysical and optical properties, they interact with solar and terrestrial radiation, alter cloud amount and radiative properties (e.g., [1][2][3][4][5]), fertilize ocean and land, and regulate carbon uptake (e.g., [6,7]) with ensuing impacts on climate (e.g., [8][9][10]). The sensitivity of radiative forcing to the variations of the aerosol optical properties defines the impact of aerosols on climate change and must be taken into account in climate scenarios (e.g., [11]). Yet, it is still difficult to obtain an accurate estimate of the radiative impact of aerosols. The scarcity of information on aerosol optical properties is one of the greatest uncertainties in the evaluation of climate forcing (e.g., [12][13][14][15]). Map of the study area (North Tropical Atlantic Ocean), which is divided into 3 sub-regions: the West zone corresponds to the North East Brazilian zone, the East zone corresponds to the Africa zone and the Centre zone is located between these two zones. The color scale indicates the average value of the aerosol optical depth τ a at 550 nm as derived from the POLAC/PARASOL satellite sensor time-series (2005)(2006)(2007)(2008)(2009)(2010)(2011)(2012)(2013). The geographical location of the AERONET ground-based stations that are used in this study is also reported in the figure.

POLDER/PARASOL
The PARASOL satellite (CNES), which was operational from 2005 to 2013 within the A-Train constellation, carried the third generation of the POLDER sensor (POLDER-3). POLDER-3 is a passive optical imaging radiometer and polarimeter, which consists of a rotating wheel carrying filters that allow measurements at 9 wavelengths (443, 490, 565, 670, 763, 765, 865, 910 and 1020 nm). Three of the channels, namely 490 nm, 670 nm, 865 nm measure the polarization of the incident light. In addition to its ability to measure the polarization state of light, the POLDER-3 sensor is characterized by its wide field of view telecentric optics ( 114 • ) which allowed observation of the entire Earth's surface every 1 or 2 days at least. Interestingly, such a wide field of view also enabled observation of a given ground target along the satellite track at various viewing zenith angles, typically up to 16 directions of observation, thus permitting the measurements of the radiation coming from a target at different scattering angles. Note that due to signs of aging (notably its star tracker), the PARASOL satellite orbit moved off the A-train's track in December 2009 and its orbit was lowered 4 km below the A-Train from 2009. Since then, PARASOL orbit was continuously lowered until the end of the PARASOL mission. Such modification of the PARASOL orbit caused the temporal window of PARASOL to be delayed about 1 h at the end of 2010. A total number of 2939 PARASOL images acquired over the study area has been exploited here.
Over the oceans, the PARASOL operational algorithm [61] assumes non-absorbing particles, spherical or non-spherical particles. The size distribution of aerosols is modelled as a combination of two lognormal distributions, one in the fine mode (sub-micron size) and one in the coarse mode (effective radius typically larger than 1 µm [45,46]). POLDER/PARASOL products have already been assessed in Bréon et al. [41]. In Bréon et al. [41]'s study, an evaluation of the accuracy of different aerosol products derived from 5 sensors (POLDER, MODIS, MERIS, SEVIRI and CALIPSO) was made via a statistical comparison with AERONET data during a period of 5 years (March 2005-July 2010) on a global scale. It included the distinction between total and fine mode AOD. Their method revealed a good performance from POLDER with a correlation with AERONET around 0.9 for total AOD and 0.77 for fine mode AOD. It also showed that POLDER and MODIS retrievals are of similar quality over the ocean, even if the polarization-based retrieval using POLDER data offers a better fine mode estimate than MODIS. Finally, Bréon et al. [41] reported that the other satellite products are of lesser quality, strengthening our choice to use POLDER measurements.

Brief Overview of the POLAC Algorithm
To better exploit the multidirectional and polarization measurements made by POLDER-3, the POLAC atmospheric correction algorithm was developed (Harmel and Chami, 2011). The POLAC method uses as inputs the measurements of the multidirectional and polarized properties of the top of the atmosphere signal to retrieve atmospheric parameters, such as aerosol optical depth, fine and coarse modes of aerosol, and the water-leaving radiance over the ocean. To do so, the POLAC algorithm relies on the invariance of the POLDER/PARASOL radiance and polarization data to the concentration of in-water constituents (such as phytoplankton in open ocean waters), especially in the visible part of the spectrum such as 490 and 550 nm [50,52]. The POLAC algorithm is also built using an optimization scheme, which is based on the minimization of the Euclidian distance between the POLDER/PARASOL directional data, and on forward radiative transfer simulations in the coupled atmosphere-ocean system. Note that the POLAC method was preferred in this study to the so-called GRASP algorithm [62], which is also able to retrieve the aerosols properties from multidirectional and polarized satellite data, because of its consideration of the oceanic hydrosols optical properties during the inversion process.
For the purpose of this study, the Look Up Table (LUT), which is used to describe the aerosol models within POLAC algorithm, has been updated since the first version of POLAC algorithm published in 2011. In particular, the absorption of the aerosols (fine and coarse modes) is taken into account. The non-spherical aerosol model that has been proposed by Volten et al. [63] and has been employed by Herman et al. [61], is used here as well. Note that such a non spherical model could be representative of dust-like aerosols. Table 1 reports the aerosol models used in the paper. The POLAC algorithm does not take into account potential variations of the vertical profile of aerosols. The vertical distribution of aerosols is modeled as exponentially decreasing with respect to the height, with 2-km scale height for both fine mode and coarse mode aerosols. Gordon [64] showed that the impact of the vertical distribution of aerosols (absorbing or not) on the top of atmosphere satellite reflectance is practically null in the near infrared spectrum independently of the aerosol amount and model. Since the POLAC algorithm uses the radiance measured by PARASOL in the near infrared band to derive the aerosol optical properties [43], the impact of the assumption using a fixed value of 2 km for the aerosol height scale should thus be weak. However, the POLAC algorithm also uses information measured in a visible band, namely 490 nm, coupled with the signal measured in the near infrared band. Duforêt et al. [65] showed that the vertical distribution of strongly absorbing aerosols could induce a maximum difference in the top-of-atmosphere unpolarized radiance up to 7% at 490 nm for a turbid atmosphere (aerosol optical depth value of 0.6) when the aerosol height scale value is fixed at 2 km (as in POLAC algorithm) and when the aerosols are located at an altitude of 4-5 km (such as dust-like aerosols) (see Table 1 in Duforêt et al. [65]). Such a difference could potentially lead to errors in the retrieval of aerosol properties by the POLAC algorithm. First, it should be highlighted that the POLAC algorithm uses the polarized radiance (i.e., Stokes parameter Q and U) at 490 nm, which is much weaker in magnitude than the unpolarized top-of-atmosphere radiance. Second, the radiance measured in the near infrared band remains the most sensitive band to derive aerosol properties at first order. Third, the validation of the POLAC-derived aerosol optical depth, which has been performed using ground-based AERONET measurements located at various contrasted sites, namely Capo-Verde (West Africa coast, dust-like aerosols), Lampedusa (Mediterranean Sea, urban aerosols) and Tahiti (Pacific Ocean, maritime aerosols), showed a satisfactory agreement as it could be observed in Harmel [66] (see omtab.obs-vlfr.fr/fichiers_PDF/Harmerl_PhD_09.pdf, Figure 4.12a). As a result, the POLAC derived aerosol products are reliable, at least for conditions where aerosols are not strongly absorbing, despite the assumptions made on the vertical distribution of aerosols.
Harmel [66] showed that the intrinsic uncertainty associated with the POLAC inversion of PARASOL signal is estimated at 0.5%, based on simulations. Then, Harmel and Chami [43] especially highlighted that the benefit of using multi-directional observations of a given target, as POLAC algorithm does, consists in providing reliable estimates of aerosol type, including the Ångström exponent that is calculated between 490 nm and 865 nm, with an uncertainty lower than 4%. POLAC is not the only algorithm to retrieve the fine and coarse fractions of POLDER-3 aerosol products over ocean. The operational level 2 POLDER-3 algorithm developed by Herman et al. [61] also discriminates these two modes over ocean and is employed by in Tanré [45], Tanré et al. [46] and more recently in Formenti et al. [67].
Our study is based on the use of AOD daily products retrieved by the POLAC algorithm during the lifetime of POLDER/PARASOL sensor, namely between March 2005 and September 2013. These data are derived at 550 nm at the ICARE Data and Service Centre (http://www.icare.univ-lille1.fr). It would have been interesting to compare these POLAC/PARASOL products to the operational level 2 POLDER-3 algorithm fine and coarse mode AOD retrievals products. However, such a work is beyond the scope of the current study.
In the rest of the paper, the notation "PARASOL AOD" is used for the total aerosol optical depth (i.e., fine mode AOD + coarse mode AOD) as observed by the POLDER/PARASOL sensor and retrieved by POLAC algorithm ; the notation "PARASOL AOD f " is used to refer to the aerosol optical depth of the fine mode ; the notation "PARASOL AOD c " is used to refer to the aerosol optical depth of the coarse mode. The fact that the fine and coarse mode of aerosol could be distinguished here in PARASOL AOD is relevant for this study to identify the type of aerosols over the Atlantic Ocean and thus, to quantify the spatial distribution and transport of different types of aerosols.
Note that the 6 km resolution POLAC/PARASOL product was projected with conservative re-gridding onto a 25 km grid to correctly perform comparisons with MODIS observations and to investigate the relevance of MACC and MERRA-2 reanalysis products. Table 1. Microphysical properties of the fine and the coarse modes of the spherical aerosol models used in this study by POLAC algorithm for the processing of PARASOL time-series data. The mean radius r i , and the standard deviation σ i are given in µm. The microphysical properties of the non-spherical model are not indicated in the table since its optical properties were directly obtained from laboratory measurements by Volten et al. [63].

AERONET Ground Based Measurements
AERONET (AErosol RObotic NETwork) is an internationally federated ground-based instrument network of nearly 300 sites regularly reporting data, and over a thousand automatic sun and sky scanning Cimel (Cimel is NASA-AERONET's exclusive supplier of automatic Sun Sky Lunar photometers (CIMEL CE318-T) operating in near real time and providing aerosol optical and columnar microphysical properties.) radiometers, measuring and characterizing aerosol properties including aerosol optical depth [68]. Direct sun measurements are carried out in several spectral bands between 340 and 1640 nm to compute the total column atmospheric aerosol optical depth. Away from the Sun, measurements scan the sky and provide measurements of the radiation scattered down to the radiometer at 4 wavelengths, namely 0.44, 0.67, 0.87 and 1.02 µm, and over a wide range of scattering angles [69]. Then, sun and sky measurements are used to retrieve aerosol size parameters including the fraction of the fine and coarse mode of aerosol following the methodology of Dubovik et al. [70]. The AERONET inversion algorithm provides improved aerosol retrievals by fitting the entire measured field of radiances (i.e., the angular distribution of radiances and sun radiance) in the solar almucantar to a radiative transfer model. The inversion allows retrieval of the wavelength-dependent aerosol complex index of refraction and the bimodal particle size distribution in the size range [0.05-15 µm], which enables calculations of aerosol properties [70]. The AOD processing includes the spectral de-convolution algorithm (SDA) described in O'Neill et al. [71]. This algorithm provides fine (sub-micron size) and coarse (super-micron size) AOD at 500 nm. The algorithm fundamentally assumes that the coarse mode Angström exponent and its derivative are close to zero. Its advantage is based on the fact that it produces useful indicators of aerosol size discrimination at the frequency of extinction measurements. In this paper, the level 2.0 quality-assured AERONET SDA retrievals (Version 2) derived at 500 nm are used. In order to compare AERONET AOD products to other data sets at the same wavelength, AERONET derived AOD products at 500 nm are interpolated to determine AERONET derived AOD products at 550 nm following Angström power law: where τ 550 is the AOD or AOD f at 550 nm, τ 500 is the AOD or AOD f at 500 nm and α is the Angström exponent of the AOD or AOD f at 500 nm. The AERONET AOD products at 550 nm are referred to as "AERONET AOD" for total aerosol optical depth, "AERONET AOD f " for the fine mode aerosol optical depth and " AERONET AOD c " for the coarse mode aerosol optical depth. To ensure their quality, the AERONET products have undergone pre-field and post-field calibration as provided by AERONET quality assurance program ; they have been cloud cleared and manually inspected as well. The typical overall uncertainty in the AERONET AOD data for field instruments is about ±0.01 to ±0.02 [72,73]. To verify the accuracy of PARASOL AOD values, near-continuous and long-term data sets of AERONET AOD, AERONET AOD f and AERONET AOD c as provided by the following four AERONET stations (https://aeronet. gsfc.nasa.gov/) from 2005 to 2013 were processed: Santa Cruz Tenerife (28.47 (Figure 1) are used. Note that continental AERONET stations located over the study area, such as Dakar, Bambey or Dahkla, cannot be used since PARASOL AOD products derived by POLAC method are delivered over the ocean only. Cape Verde and Guadeloupe AERONET stations are also not used for evaluating POLAC retrievals because of the lack of available data. Therefore, the POLAC aerosol products may not be sufficiently accurate to compare with these AERONET sites, including the Dakar site which has one of the longest data records. Note also that the AERONET stations are spatially distributed so as to cover the various sub-regions of the study area.
From 2005 to 2010, AERONET AOD data were selected every day during a 3-hour time window centered on the approximate revisit time of PARASOL satellite, namely at 14:30 UTC in the East zone and at 17:30 UTC in the West zone. From 2011 to 2013, the approximate revisit time of PARASOL satellite is shifted to 3:30 pm in the East zone and to 6:30 pm in the West zone to compensate the continuing loss in orbit of PARASOL satellite as mentioned in Section 2.2. The AERONET AOD data are then daily averaged prior to comparing with daily PARASOL AOD products.

MACC Reanalysis
Since 2003, the European Center for Medium-Range Weather Forecasts (ECMWF) performed a reanalysis of global atmospheric composition, the so-called "Monitoring Atmospheric Composition and Climate (MACC) Reanalysis". Using the aerosol size as a proxy for aerosol origin, MACC reanalysis includes aerosol optical depths that are divided into three classes: anthropogenic, mineral-dust and sea-salt. The MACC reanalysis database is constructed by assimilating satellite data at several spectral bands (including 550 nm) into the ECMWF Integrated Forecast System (IFS) global model. Current MACC aerosol reanalysis are performed using an aerosol data assimilation system developed by Morcrette et al. [74] and Benedetti et al. [75] which relies on MODIS AOD measurements and a 4D-Var approach. The spatial resolution of the MACC reanalysis data set is about 80 km on 60 vertical levels from the surface up to approximately 65 km. The MACC-II VAL sub-project validated the outputs of the reanalysis but several quality issues remain such as those regarding the sea salt aerosol mixing ratio above freshwater or changes in biomass burning emissions for example (https: //software.ecmwf.int) [76,77]. In this study, MACC total aerosol optical depth monthly mean at 550 nm, hereafter referred to as MACC AOD, is compared with the corresponding PARASOL AOD measurements over the time period 2005-2013. The fact that MACC reanalysis is able to distinguish between aerosol origins is important to gain understanding on the fine mode and coarse mode aerosol transport over the NTA ocean.

MERRA-2 Reanalysis
The Modern-Era Retrospective analysis for Research and Applications, Version 2 (MERRA-2), developed by NASA's Global Modeling and Assimilation Office (GMAO), is the latest atmospheric reanalysis of the satellite era (1980 onward) produced using the Goddard Earth Observing System (GEOS-5) atmospheric model, data assimilation system version 5.12.4 [54,78,79] and the three-dimensional variational data analysis (3DVAR) Gridpoint Statistical Interpolation (GSI) meteorological analysis scheme [80,81]. GEOS-5 model spatial resolution is roughly 50 km in the latitudinal direction with 72 hybrid-eta layers from the surface to approximately 65 km. Most products are delivered on the standard 0.5 • × 0.625 • latitude by longitude grid [82]. In comparison with the original first version of MERRA, MERRA-2 includes an aerosol analysis as described in Randles et al. [82,83].To assimilate AOD at 550 nm, an analysis splitting technique is used in which a two-dimensional analysis is achieved using error covariances derived from innovation data [82]. Then, horizontal increments are projected vertically and across species via an ensemble method [54]. New AOD observations are assimilated from several sources, including: MERRA-2 reduces spurious jumps and trends associated with changes in the meteorological observing system [84]. Using recent developments at GMAO in the field of modeling and data assimilation, MERRA-2 is now able to combine assimilation systems for land, ocean, atmosphere and chemistry. As an example, MERRA-2 includes aerosol data assimilation that can interact with atmospheric radiative processes and thus provides a multidecadal reanalysis in which aerosol and meteorological observations are jointly assimilated within a global data assimilation system [54]. In order to compare with PARASOL AOD data, different types of MERRA-2 monthly mean products at 550 nm were used here: the Total Aerosol Extinction AOD (hereafter referred to as "MERRA-2 AOD"), the Dust Extinction AOD-PM 2.5 and the Sea Salt Extinction AOD-PM 2.5 (PM 2.5, i.e., the mass concentration of aerosols with diameters less than 2.5 µm). Note that even though dust and sea salt aerosols mostly contribute to the coarse mode AOD, they could be also present within the fine mode fraction. In that sense, unlike MACC reanalysis, the fine mode of AOD due to primary dust and sea-spray are available in the MERRA-2 reanalysis products. To directly compare fine mode AOD to the PARASOL AOD f product, it is necessary to sum both MERRA-2 fine modes aerosol products, namely the Dust Extinction AOT-PM 2.5 and the Sea Salt Extinction AOT-PM 2.5 products, to get the total fraction of the fine mode of aerosols. Such a sum is hereafter referred to as "MERRA-2 AOD f ". The fraction of coarse mode of MERRA-2 AOD products, hereafter referred to as "MERRA-2 AOD c ", is calculated by subtracting MERRA-2 AOD f to MERRA-2 AOD product.
As for POLAC/PARASOL product, MACC (80 km spatial resolution) and MERRA-2 (50 km spatial resolution) AOD products used in this study were projected onto the same 25 km grid in conservative projection. Such a projection aims at making consistent the analysis of MACC data set, MERRA-2 data set and PARASOL data set.

Comparison with AERONET Measurements
PARASOL AOD products (total, fine and coarse modes) (Section 2.2) are evaluated with respect to AERONET measurements (Section 2.3) from 2005 to 2013. To compare with AERONET ground-based measurements, the highest spatial resolution of PARASOL AOD products (i.e., 6 km) has been used. Figure 2 shows the monthly time series of the total AOD difference between PARASOL observations and corresponding ground-based measurements acquired by the four selected AERONET stations ( Figure 1 3) are observed. Such seasonal variability is likely due to mineral dust advection during large Saharan dust outbreaks occurring in summer since we can observe the same but less pronounced seasonal variability in AERONET AOD c (Figure 4). At Ragged Point, the most westerly station of the study area, the seasonal variability remains observable despite an attenuation illustrated by an annual amplitude around 0.1 (Figure 2g). Unlike other stations in the East zone, the almost total lack of high daily AOD values over 0.4 (not shown) leads to the observed lower AOD variability, which is probably due to dry and wet removal processes of dust aerosols during the transport across the Atlantic Ocean. The seasonal variability is not the only temporal variability experienced by the aerosol optical depth. Because of the high frequency and intensity of aerosol production in North Africa associated with high temporal variability of meteorological parameters (precipitations, wind,..), AERONET AOD shows an important daily variability as illustrated by the high AERONET standard deviation values compared to corresponding AERONET AOD (Table 2).
Second, the comparison between the two data sets presented by Figure 2 and Table 2 show that, at all stations, PARASOL AOD values are in a good agreement with the AERONET AOD temporal variability. Yet, some discrepancies are observed with most of the biases ranging from −0.2 to 0.2, which is a consistent range based on previous studies [85,86]. Outside the bias range [−0.2; 0.2], PARASOL biases are mostly positive, which may be explained by an occasional failure in the detection of clouds [87]. In Calhau, discrepancies between PARASOL and AERONET retrievals in the SON period, as reported in Table 2, are obtained from measurements only over one season (SON 2012) thus limiting the significance of the results. Figures 3 and 4 present similar comparisons as Figure 2, but for the fine and coarse modes of aerosols respectively. Note that such a study is unusual as most of studies referring to spatial and temporal variability of aerosols are focused on the evaluation of the total AOD without making any specific distinctions between the fine and coarse modes of aerosols. While the comparison of PARASOL AOD c with AERONET measurements shows a satisfactory agreement with low biases at all stations ( Figure 4, Table 2), PARASOL AOD f presents weaker values. An overestimation of AOD f is observed in PARASOL ( Figure 3, Table 2), with high biases [−0.2; 0.2] compared to AERONET AOD f value range [0; 0.2]. This is particularly noticeable at Ragged Point (Figure 3g,h) where daily mean PARASOL AOD f can be twice as high as daily mean AERONET AOD f ( Table 2). The cause of the overestimation of the AOD f observed by POLDER/PARASOL might be due to the fact that fine mode aerosols are more difficult to observe due to their size, which implies an increase of errors in satellite observations. Figure 3 and Table 2 also show that values of AERONET AOD f are much lower than total AOD values ( Figure 2): AOD f maxima are typically lower than 0.2 (Figure 3a,c,e,g). On the contrary, AERONET AOD c values are high, typically in the range between 0 and 0.4 (Figure 4), and can even reach total AOD values at some points. These observations point out the lower contribution of the fine mode and the opposite high contribution of the coarse mode to the total extinction. Such a result is consistent with the study made by Schutz [88] that revealed significant concentrations of coarse mode of dust close to source regions. As a consequence of the reduced number of the fine mode of aerosols, the seasonal and daily variabilities of AERONET AOD f is significantly reduced (Figure 3), while those of the coarse mode still exist but are less pronounced (average of maximum values around 0.3, Figure 4) than total AOD variabilities. Figure 5 shows the results from a statistical analysis that has been performed between PARASOL and AERONET for the total, fine and coarse AOD. It should be highlighted that the relatively small number of coincident AERONET and PARASOL data could be explained by the constraints used to ensure a consistent comparison. In particular, only AERONET AOD data that are measured within a 3-h time window centered on the approximate revisit of PARASOL satellite have been selected, as mentioned in Section 2.3, rather than all AERONET data available over the course of the day.  PARASOL observations of both AOD and AOD c exhibit satisfactory correlation coefficients (greater than 0.74) and low biases ranging from −0.03 to 0.03 at all AERONET stations (Figure 5a,c,d,f,g,i,j,l). Such results indicate that the location of the stations, from East to West zone, has a low impact on PARASOL scores and the discrepancies shown in Figure 5 do not appear increase with increasing AOD. Note that the higher number of AOD c observations (120) compared to AOD and AOD f observations (108) in La Laguna (Figure 5d-f) is due to the lack of Angström exponent observations for the total and the fine mode (coarse mode Angström exponent is set to 0 for the La Laguna site), while observations of AOD products exist.  In parallel, PARASOL AOD f presents the lowest performance with correlation coefficients between 0.37 and 0.73, and biases ranging from 0.00 to 0.04 (Figure 5b,e,h,k). At Tenerife, the low correlation coefficient (0.40) associated with zero bias can be explained by the compensation between both negative and positive bias which tends to get an average close to 0. Figure 5 also shows that, while PARASOL AOD f biases are positive, PARASOL AOD c are negative, a trend that has also been seen in Table 2. This may lead to a slight underestimation of the coarse mode AOD and an overestimation of the fine mode AOD from PARASOL. Different studies have already shown that in situ observations often lead to significantly larger mean particle sizes than retrievals from remote sensing instruments, thus impacting the coarse mode size distribution [89][90][91][92]. It may explain the slight underestimation of PARASOL AOD c . The cause of this discrepancy is still unknown. Kleidman et al. [69] suggested that it may come from some ambiguity in the definition of the radius-threshold for fine and coarse mode between the sunphotometer and the satellite products [69].
Despite the possible accumulation of several error sources such as the difference in the spatial resolution between ground-based and satellite data sets (6 km for the PARASOL AOD against few meters for the AERONET AOD) and such as occasional errors in ground based measurements and the PARASOL cloud detection, PARASOL shows a good performance in the observation of the total and coarse mode AOD over the NTA ocean. The satisfactory performance of PARASOL total AOD is also consistent with Harmel [66] on POLAC inversion made on PARASOL AOD (Section 2.2.2). Our analysis therefore highlights the ability of PARASOL to detect the contribution of coarse aerosols, and to a lesser extent, fine aerosols to the total AOD which are innovative contributions over such a region.  Figure 6 and Table 3 compare the AOD, AOD f and AOD c between PARASOL observations and the Level-3 MODIS Atmosphere Monthly Global Product MYD08_M3 collected from the Aqua platform (https://modis-atmos.gsfc.nasa.gov/products/monthly), averaged from 2005 to 2013. Note that the use of MISR satellite sensor products [93], which has a finer spatial resolution (4.4 km) than MODIS would have been interesting. However, the goal of this paper is not to compare all the size-resolved AOD products and the use of MODIS data was preferred here, in particular because the MACC and MERRA-2 reanalysis aerosol datasets used in this paper are based on the assimilation of MODIS data. Since PARASOL AOD products derived by POLAC method are not delivered over the land, this study only uses MODIS products over the ocean. It also should be highlighted that the Multi-Angle Implementation of AtmosphericCorrection (MAIAC), which is designed to work with MODIS data to facilitate deriving both aerosol and land surface reflectance products, is a relevant method [94] that has not been used here since our study deals with data collected over the ocean surface. In this study, PARASOL coarse mode aerosols are defined by an effective radius larger than 1 µm while MODIS coarse mode aerosols are defined by an effective radius larger than 0.98 µm (https://modis.gsfc.nasa.gov/data/atbd/atbd_mod02.pdf). Our results show that PARASOL and MODIS data agree reasonably well. For both of them, the AOD decreases by a factor of ∼1.5 (Table 3) from the coast of North Africa westward to the Caribbean Sea, presumably due to wet and dry depositions of mineral dust during transport and horizontal spreading. However, a few differences between PARASOL and MODIS do exist. First, the transport of aerosols over the NTA ocean seems to be less strong in MODIS observations with a horizontal dispersion of aerosols that is not as extended as PARASOL transport (Figure 6a,d) which is illustrated by total AOD values equal to 0.18 for MODIS against 0.23 for PARASOL in the West sub-region (Table 3). This can be explained by the lower amount of coarse mode aerosols observed by MODIS, especially in the West sub-region (Figure 6c,f, Table 3), which implies a decrease of MODIS total AOD over the study area, and more specifically over the West sub-region. Second, the highest AOD values derived from MODIS off the coast of Mauritania (Africa) are observed further southward with PARASOL observations, both for AOD and AOD f products, and to a lesser extent in AOD c . Usually, important features of AOD along the coast of the Gulf of Guinea are often attributable to biomass burning aerosols, but large amounts of dust aerosols, transported from sources further inland, could also be present. It may indicate that biomass burning aerosol is not well observed using POLDER/PARASOL sensor, as indicated by the largest bias between PARASOL AOD f and MODIS AOD f found in winter (0.03) in regard to other seasons. Perhaps, this is because of a less strict PARASOL cloud mask which leads to overestimated AOD values by PARASOL. The area off the coast of Mauritania is also subjected to a coastal upwelling where strong winds blow across the surface of the ocean thus making surface waters diverge horizontally [95]. In such area, the differences observed between POLDER/PARASOL and MODIS sensors could possibly come from discrepancies between the inversion algorithms used for each sensor, such as the wind retrieval algorithm. In particular, a different whitecap coverage (a parameter that is function of the wind speed), the stability of the lower atmosphere and changes in the sea-surface temperature could impact the observed apparent microwave temperature and short-wave albedo of the sea surface [96], thus leading to different observations of AOD between both satellite sensors. Despite the overestimation of PARASOL AOD f observed in the comparison with AERONET measurements, the fine mode of PARASOL shows a good consistency with MODIS AOD f (Figure 6b,e and Table 3). This may indicate the difficulty to derive size information from the spectral and directional reflectance [41]. Kalashnikova and Kahn [36] also point out that low optical depths (below 0.5), that are predominant in fine mode AOD, combined with assumptions made in the satellite retrieval algorithms about ocean surface optical properties (reflectance and emissivity) [97], may lead to an overestimation of the AOD in MODIS. This could be the case as well in PARASOL.
Finally, the overall satisfactory comparison between PARASOL AOD observations and MODIS dataset confirms that POLDER/PARASOL sensor is relevant to study the transport of aerosol over the NTA ocean.

Seasonal Variability and Climatology of AOD Observed by PARASOL
Aerosol production in North Africa and exportation over the NTA ocean is characterized by a pronounced seasonal evolution [98,99]. To analyze the seasonal variability, Figure 7 presents the climatology of PARASOL products (AOD, AOD f , AOD c ) at 25 km spatial resolution for each season and averaged over the whole time series 2005-2013. At seasonal scale, the production of different types of aerosols coming from different regions of Africa leads to both a spatial and temporal variability of the amount of aerosols in the atmosphere [56]. Thus, seasonal changes in spatio-temporal aerosol activity are correlated with seasonal changes in local meteorological conditions providing favorable atmospheric conditions for aerosol emission (e.g., [100]). Most of the aerosol variability comes from spatio-temporal dust production and, to a lesser extent, from biomass burning aerosols production. For the latter, even if vegetation fires could occur over the course of the entire year in Sub-Saharan Africa, fires are infrequent during the wet season. Therefore, the largest amount of biomass burning aerosols are emitted from the Gulf of Guinea in equatorial Africa during the dry season [101] and thus, part of the winter plume over the NTA ocean is due to smoke (Figure 7a). This variability is mainly driven by the oscillation of the ITCZ leading to the succession of two meteorological regimes over the Sahel during the year [102]: during summer, the Sahel is submitted to the wet monsoon flow from South-West due to the northern shift of the ITCZ. An occurence of extremely high daily concentrations of dust is induced by strong surface winds associated to Mesoscale Convective Systems that cross the Sahel. On the opposite side, in winter, the northeastern dry Harmattan winds allow a Saharan dust transport towards the Gulf of Guinea. This changeover explains the south-north migration of the main transport corridor for the aerosol outflow from Africa observed in Figure 7: from under 15 • N in winter (Figure 7a), with the plume extending to South America, to 10 • N-30 • N in summer (Figure 7g) with an export into the Caribbean and toward Florida [103][104][105]. In spring, the amount of aerosols reaches its maximum with PARASOL AOD greater than or equal to 0.6 near the coast of Africa and over 0.4 at 40 • W (Figure 7d) due to the highest values of PARASOL AOD c in this season (Figure 7f, Table 3). In summer, PARASOL AOD increases in the Centre zone and especially in the West zone, showing the highest transport of aerosol across the ocean over the year (Figure 7g, Table 3). This result from a summer transport at higher elevations (up to 7 km) within the SAL combined with higher wind speeds compared to winter transport [92,106]. In autumn, a significant decrease of PARASOL AOD appears with values lower than 0.4 over the entire domain (Figure 7j) marking the end of the summer intensive period of aerosol transport over the NTA ocean.
Consistent features are found for PARASOL AOD c with lower values (Figure 7, right column): PARASOL AOD c varies between 0.1 and 0.4 whereas PARASOL AOD can reach 0.6 according to season. High values of PARASOL AOD c that are observed over 10 • N also confirm the important part of the coarse mode in the dust aerosols coming from the SAL. Such results attest that the aerosol spatio-temporal variability is primarily under the influence of coarse dust aerosol spatio-temporal variability and show that the major contribution of the aerosol coarse mode, already observed in Section 3.1.1 at four AERONET stations, is found over the entire study area.
Regarding the fine mode of PARASOL AOD, Figure 7 (middle column) and Table 3 show values below 0.3 throughout the year and presents some differences in spatio-temporal variability with the coarse mode AOD and thus, with the total AOD. First, PARASOL AOD f maximum value is reached in winter where the aerosols originate from the equatorial Africa ( Figure 7b) ; second, PARASOL AOD f observed in summer at equatorial Africa also prevails over PARASOL AOD f that originates from the North of Africa, despite the high emission of aerosols at this period of the year (Figure 7h). Such an observation suggests that an important source of PARASOL AOD f is located in equatorial Africa. Since biomass burning aerosols are mainly found in this area, it does thus make sense to believe that biomass burning aerosols are composed of fine mode aerosols. However, it should be highlighted that the aerosols used in the POLAC algorithm are not as absorbing as biomass burning aerosols can be. As a result, the occurrence of biomass burning aerosols could alter significantly the performance of POLAC retrieval. PARASOL AOD f and PARASOL AOD c are also differentiated by their transport over the ocean. While PARASOL AOD f is fairly zero over 30 • W westward in winter (Figure 7b), PARASOL AOD c achieves to cross the Atlantic Ocean (60 • W) in spring and summer (Figure 7f,i). The much more important transport of coarse mode aerosols as observed through PARASOL AOD c product is a significant result verified in previous studies [10,30,36,40,92,107,108]. This should be correlated with the fact that a high contribution of the coarse mode relative to the total AOD was observed over the whole area, as mentioned earlier in the paper. The transport of coarse aerosols over long distances could be explained by the greater number of emissions sources of large aerosols, which statistically increase the probability of coarse mode aerosols reaching the West zone. Apart from Africa, Figure 7 presents a much smaller source of aerosols in the north of Brazil (0 • N ; 50 • W) (Figure 7, left column) which exhibits the occurrence of coarse mode aerosols only (Figure 7, right column). Because of their limited transport in this West sub-region, coarse mode aerosols emitted from Brazil have a negligible impact on high values of AOD c found in the East and Centre zones (Figure 7c,f,i,l). Aerosol coarse and fine modes vary depending on the sub-region and on the season. To study this spatio-temporal variability, Figure 8 presents the interannual variability of PARASOL AOD, AOD f and AOD c monthly means for each sub-region (East zone, Centre zone and West zone). Even though AODc has the largest impact on AOD over the whole study area, it can be equalized or surpassed by AOD f at the end of every year in the East zone (Figure 8a), due to the large amount of emitted biomass burning aerosols from the Golf of Guinea during the dry season. In the Center zone, and even more in the West zone, the influence of AOD c on AOD is emphasized by the weak transport of PARASOL AOD f , which is illustrated by the gap between AOD c and AOD f (Figure 8b,c). Thus, cases where AOD f contribution reaches AOD c values occur nine times in the East zone, two times in November 2006 and July 2011 in the Center zone ( Figure 8b) and no longer occur in the West zone (Figure 8c).

Implications for Models Simulations
The previous observations of the aerosol transport, and more specifically the aerosol coarse mode transport, across the NTA ocean is now analysed using reliable reanalysis retrievals that combine modelling and assimilation of AOD observations: MERRA-2 and MACC reanalysis. According to Bellouin et al. [109] and Randles et al. [82], the uncertainties of MACC and MERRA-2 reanalyses are close to 0.03 and 0.013 respectively for the total AOD. The uncertainties are close to 0.014 and 0.010 for MACC and MERRA-2 respectively when dealing with dust AOD f at 550 nm. Results presented below are calculated only based on coincident PARASOL/MACC/MERRA-2 observations. Figure 9 presents the spatial distributions of AOD, AOD f and AOD c derived from PARASOL (Figure 9a-c) and from MERRA-2 reanalysis (Figure 9d-f). The averaged AOD from 2005 to 2013 as derived from MACC reanalysis is also reported (Figure 9g). In parallel with Figure 9, Table 3 compiles seasonal and annual averages by sub-regions (East, Center and West) of the three AOD products (total, fine, coarse) from PARASOL, MERRA-2 and MACC.
The comparison with MERRA-2 reanalysis shows a significant difference with PARASOL observations: smaller amounts of aerosols are provided by MERRA-2 over the whole domain (Figure 9a,d, Table 3). Such a result could be due to several factors. First, MERRA-2 presents a lower average content of aerosols over the East zone (0.26 for MERRA-2 against 0.33 for PARASOL, Table 3) which may be due in part to an underestimation of winter fine mode aerosol production, i.e dust fine mode fraction and biomass burning aerosol under 10 • N (Figure 9e) (winter MERRA-2 AOD f : 0.10, winter PARASOL AOD f : 0.16, Table 3). Second, Figure 9d presents a much stronger AOD longitudinal gradient with a faster decrease in the MERRA-2 AOD value from East to West, compared to PARASOL (Figure 9a), thus highlighting a considerably more limited MERRA-2 transport. Indeed, the underestimation of MERRA-2 aerosol content in the East zone cannot entirely be the cause of the low MERRA-2 AOD values found further west since MERRA-2 AOD negative bias, in regard to PARASOL AOD, increases with the distance from african sources (−0.07 in the East zone; −0.09 beyond 30 • W). Note that when MERRA-2 dust content is higher than PARASOL content (in summer), MERRA-2 transport of dust remains smaller. A highly underestimated MERRA-2 coarse mode transport, observed in Figure 9f, could explain the origin of this limited transport of aerosol.  In parallel, it should be highlighted that since MODIS AOD observations have been assimilated in MERRA-2 reanalysis, a meridional shift of the MERRA-2 maxima of total AOD, AOD f , and to a lesser extent AOD c , is expectedly observed in comparison with PARASOL AOD products. However, it is worth noting that MERRA-2 underestimates the aerosol fine mode fraction under 10 • N (Figure 9e) whereas MODIS sensor could observe it ( Figure 6e). As for the coarse mode, the limited transport of MERRA-2 AOD c (Figure 9f) may be partly explained by the weaker AOD c transport in MODIS observations as that observed with PARASOL (Figure 6c,f). However, MERRA-2 process of data assimilation could also potentially be involved, owing to MERRA-2 AOD c bias with MODIS AOD c ( Table 3).
Compared to MERRA-2, MACC reanalysis presents a transport of aerosols across the ocean (Figure 9g) closer to the one observed in PARASOL (Figure 9a). This is illustrated by a constant bias around −0.4/−0.3 in each zone in regard to PARASOL AOD (Table 3). Such better performance occurs despite MACC lower aerosol activity in the East zone where a weaker amount of dust ( Figure 9g) is observed, compared to MERRA-2 and PARASOL (Figure 9a,d). Although MACC and MERRA-2 reanalysis differ on several aspects, such as the dust retrieval or the aerosol transport across the ocean, they converge on the estimated content of aerosol coming from equatorial Africa. However, with weaker values than those observed by PARASOL and MODIS, the estimated production may be underestimated. MACC reanalysis also provides the advantage of distinguishing between different types of aerosols (mineral dust, anthropogenic aerosol and sea salt aerosol) whose spatial distributions are shown in Figure 10. These simulations indicate that, even if sea salt AOD has generally much lower values compared to dust (ranging from 0.10 to 0.20) (Figure 10c), AOD reaches the value 0.10 in the central part of the Atlantic (Figure 10d) thus highlighting the non-negligible role of sea salt aerosols in this area in terms of contribution to the total extinction (Figure 10a). The contribution of anthropogenic aerosols, notably biomass burning aerosols, is limited to a more restricted area that is located close to the source of aerosol production (Figure 10b). The simulated limited transport of fine mode aerosols strengthens the assumption previously made by Adams et al. [110]. Our results confirm the important role of dust coarse aerosols in the field of aerosols transport study. In summary, the comparison of the aerosol transport between PARASOL and model reanalysis outputs revealed that models tend to underestimate the aerosol transport over the NTA ocean, which is probably due to an underestimation of the coarse mode transport of aerosols.

Discussion
Over the NTA ocean, the quarterly variations of a long time series of satellite observations obtained by PARASOL processed using the new POLAC atmospheric correction algorithm provide a satisfactory overview of the seasonal variability of aerosols, allowing the distinction between total, fine and coarse mode aerosols. Driven by the seasonal shift of the ITCZ and its impact on the SAL export [111], the aerosol seasonal variability is mainly induced by mineral dust advection over the ocean as shown by the large contribution of PARASOL AOD c relative to the total AOD and by their location over 10 • N. The high concentration of coarse mode aerosols close to the source of production and its decrease during transatlantic transport to the Caribbean are well observed by PARASOL. It strengthens the result obtained from previous works related to dust transport (e.g., [23][24][25]92]); some of which use other satellite sensors like MODIS and IASI [8,41,112]. Based on the performance of POLAC/PARASOL observations of AOD and AOD c assessed through a comparison with ground-based AERONET measurements, this study thereby confirms the enhanced long-range transport of coarse particles, such as mineral dust, over our study area.
Note that the coexistence of mineral dust and clouds could exist. However, since our satellite and AERONET data were analysed for cloud-free conditions only, it is not possible in this study to examine the consequences of such a coexistence on the comparison between observations and models. This is one limitation of our approach.

Comparison with Previous Experimental Campaigns and Satellite Data
Several field campaigns have been realized on dust aerosol transport across the Atlantic Ocean to the Caribbean. Under the Saharan Aerosol Long-Range Transport and Aerosol-Cloud-Interaction Experiment (SALTRACE) project, ground-based and airborne in situ remote sensing observations in Barbados, Puerto Rico and Cabo Verde were conducted from spring 2013 to summer 2014 [92].
Weinzierl et al. [92] exhibit the major contribution of dust to the AOD in the Caribbean and detect the unexpected presence of 10-30 µm particles even after 5 days of transport and more than 4000 km. Based on Denjean et al. [113] showing that modal peak diameter of the volume size distribution remained constant on both sides of the Atlantic, Weinzierl et al. [92] suggest that after 2-3 days from uplift gravitational settling is almost ineffective. These conclusions are consistent with Maring et al. [30] who, from comparisons between surface based measurements at sites in the Canary islands (July 1995) and in Puerto Rico (July 2000), conclude that some atmospheric processes partially counteract gravitational settling of large dust particles. Velasco-Merino et al. [107] also evaluated the impact of the transport over the Atlantic Ocean on Saharan dust optical and microphysical properties based on AERONET data. They show a decrease of about 30% in the coarse-mode concentration of the size distributions but no significant change in the coarse-mode effective radius during the transport. All these studies present total AOD values similar to our results with maximum peaks reaching 0.6 in summer near the west coast of Africa and around 0.3 in the Caribbean.
To characterize aerosol transport over the Atlantic Ocean, several other studies rely on satellite sensors and their inherent advantages over ground-based measurements due to their spatial and temporal coverage. Among them, Kalashnikova and Kahn [36], Yu et al. [35] and Prospero et al. [10] have described the significant dust transport across the tropical North Atlantic over the course of one a year. They all show consistent seasonal variations of trans-Atlantic dust transport but present notable differences in the magnitude of dust AOD, with PARASOL dust transport higher than MODIS and much higher than MISR [35]. Our study stands apart by distinguishing the fine mode of AOD, which may include the contribution of several types of aerosols, from the coarse mode of AOD, which is likely dominated by dust [16,17]. It should be mentioned that those previous studies also did not perform such an extensive evaluation as the one presented in this paper since their investigations did not concern such a long period of time series or such a large spatial extent.

Possible Causes of the Long-Range Transport of Coarse Mode Aerosols
One of the most important factors that determine the transport of coarse mode aerosols is the vertical distribution (e.g., [114]). During summer, coarse mode aerosol transport across the Atlantic Ocean mainly occurs in the SAL. As shown by CALIOP [34,103,110,114,115], the SAL is located above the top of the marine boundary layer and typically 3-4 km high. Thus, the transport of the coarse mode could cover great distances. This transatlantic transport ends in the Caribbean where aerosols are transferred from the SAL to the marine boundary layer, by gravitational settling of mainly coarse mode particles and through convective and turbulent downward mixing [92]. This makes it different from, for example, fine biomass burning aerosols which are less influenced by the trade wind since the latter restricts their transport [116] over the ocean. The seasonal meteorological regime may explain the long-distance transport of the coarse mode aerosol by causing the elevation of the dust layers over the West Africa at higher levels during summer than during winter [117]. Indeed, while long-range dust transport occurs in summer at high altitude within the SAL above the trade winds inversion (~1.5 km above sea level), the winter time long-range dust transport takes place in a low layer within the trade winds [118]. However, Figure 7c shows that coarse mode aerosol transport exceeds the fine mode aerosol transport even during winter. Such differences in long-distance transport must be correlated to aerosol fine mode or coarse mode attribution.

Long-Range Transport of Coarse Mode Aerosols in Reanalysis
Our study also highlights the difficulty for model reanalysis to properly estimate the aerosol optical depth over the study area, despite the use of data assimilation. Kim et al. [115] reached a similar conclusion by presenting a multimodel analysis and comparison with remote sensing data. Contrary to satellites which show general agreement, Kim et al. [115] revealed that models display considerable discrepancy with each other, with significant differences between the models and observation as well as among the models. In our study, such discrepancy is illustrated by an underestimation of the aerosol transport over the ocean, more or less important regarding the type of model reanalysis (MACC or MERRA-2). Even though MACC AOD values are slightly lower than observations, MACC reanalysis shows a satisfactory transport of aerosols over the NTA ocean (Figure 9g) while it is not the case for MERRA-2 reanalysis. It seems that the inability of MERRA-2 to reproduce a long distance transport as high as measured by POLDER/PARASOL comes from a significant underestimation of the coarse mode aerosol transport over the ocean (Figure 9f). The underestimation occurs despite the fact that assimilation of accurate MODIS observations, of both AOD and AOD c , are performed in MERRA-2 reanalysis (Figure 6d). It might suggest a weakness either in MERRA-2 process of data assimilation or in MERRA-2 wet and dry deposition parameterization. Such a deficiency in the evaluation of coarse mode particles along the trans-Atlantic transect has already been highlighted by Ansmann et al. [119] and Kok [120] in climate models. According to Ansmann et al. [119] based on SALTRACE campaign observations, an important removal of dust is observed in a regional model (SKIRON) and two global atmospheric models (MACC/CAMS, NMMB/BSC-Dust), as in Weinzierl et al. [92], where the main process of vertical exchange in the well-stratified SAL relies on particle settling through sedimentation. Gravitational settling depending on the simulated fine and coarse dust fractions [119], this result is coherent with Kok [120]. Indeed, the latter shows an underestimation of the coarse size fraction in climate models which may result from the sampling of dust specific size classes (PM10, PM2.5), as it is performed in MERRA-2. However, these possible causes cannot be the only explanations of the underestimation of the coarse mode particles along the trans-Atlantic transect in models. Mechanisms than retain these particles in the atmosphere during long-range transport must be involved [119].
The higher performance of MACC reanalysis, relatively to MERRA-2, for deriving the AOD is consistent with recent studies [82,83]. Randles et al. [83] suggested that MERRA-2 absorption is overestimated in areas far from aerosol sources (e.g the free troposphere) because of excessive amounts of black carbon aerosol and because the probability of precipitation of aerosols that reach the free atmosphere is almost null.
Yet, it is crucial to accurately assimilate the coarse mode of aerosols in models reanalysis. Aerosols with diameters around 0.05-10 µm are of great interest because of their direct interaction with solar and infrared radiation. Their overestimation or underestimation in climate models can likely be the causes of errors in the magnitude, and possibly, the sign of the modeled aerosol radiative forcing, depending on local parameters such as the surface albedo. A wrong consideration of aerosol coarse mode could result in erroneous numerical weather forecasting and regional climate predictions in areas which are affected by aerosol emission or transport. As an example, the consequences resulting from errors in the estimation of the coarse/fine mode of aerosols may lead to: • a possible error on the radiative forcing exerted by dust aerosols at the top of the atmosphere in case of erroneous retrievals of fine mode aerosols: less fine aerosols could lead to less radiative cooling; • a possible underestimation of warming at the top of the atmosphere in the longwave spectral range since more coarse dust aerosols lead to a higher warming; • a possible important change in the ability of dust (aged) aerosols to act as aerosol-cloud condensation nuclei.

Conclusions
The POLDER/PARASOL ability to measure the polarization state of light was exploited here to analyze a time series (from 2005 to 2013) of the amount of fine and coarse modes aerosols over the NTA ocean using a dedicated algorithm, namely the POLAC. In addition to the consistency between PARASOL data and MODIS observations, it was shown that POLAC/PARASOL products satisfactorily compared with AERONET ground-based measurements, with correlation coefficients over 0.75 and bias ranging from -0.03 and 0.03 for PARASOL AOD and PARASOL AOD c products respectively. Daily PARASOL aerosol products were then used to analyse the spatio-temporal variability of the total, fine mode and coarse mode aerosols. It clearly underlined a significant contribution of the dust coarse mode relative to the total AOD (from 15% to 65%) over the whole study area, which is consistent with previous studies using in situ observations from experimental campaigns or satellite data. Such significant transport of coarse dust mode aerosols observed by PARASOL and MODIS dataset is not reproduced by MERRA-2 reanalysis which showed a negative bias (around −0.06 beyond 30 • W). At this stage, it is difficult to identify the possible causes in the differences observed in the dust coarse mode between PARASOL and MERRA-2 reanalysis. However, as described recently by Weinzierl et al. [92] and Ansmann et al. [119] using SALTRACE observations, the daytime convective mixing within the SAL would allow a significant fraction of coarse mode aerosols to reach the Caribbean by balancing the sedimentation. Then, one reason could be related to a deficient reproduction of the dynamics in the SAL layer in models. Our study also confirms that satellite polarimetric and multidirectional daily observations, such as the POLDER/PARASOL sensor or the forthcoming "Multi-viewing Multi-channel Multi-polarization Imaging (3MI)" sensor (ESA/EUMETSAT) and the PACE sensor (NASA), should be helpful for providing additional constraints to evaluate the representation of fine/coarse dust aerosols in global and regional models.