Particle Size Analysis of African Dust Haze over the Last 20 Years: A Focus on the Extreme Event of June 2020

: Over the last decades, the impact of mineral dust from African deserts on human health and climate has been of great interest to the scientiﬁc community. In this paper, the climatological analysis of dusty events of the past 20 years in the Caribbean area has been performed using a particulate approach. The focus is made on June 2020 extreme event dubbed “Godzilla”. To carry out this study, different types of data were used (ground-based, satellites, model, and soundings) on several sites in the Caribbean islands. First, the magnitude of June 2020 event was clearly highlighted using satellite imagery. During the peak of this event, the value of particulate matter with an aerodynamic diameter of less than 10 µ m ( PM 10) reached a value 9 times greater than the threshold recommended by the World Health Organization in one day. Thereafter, the PM 10, the aerosol optical depth, and the volume particle size distribution analyses exhibited their maximum values for June 2020. We also highlighted the exceptional characteristics of the Saharan air layer in terms of thickness and wind speed for this period. Finally, our results showed that the more the proportion of particulate matter with an aerodynamic diameter of less than 2.5 µ m ( PM 2.5) in PM 10 increases, the more the inﬂuence of sea salt aerosols is signiﬁcant.


Introduction
In the literature, it is well known that Africa is the world's largest source of mineral dust [1,2]. Every year, millions of tons of mineral dust are transported from the African continent towards the Atlantic ocean in the Saharan Air Layer (SAL), i.e., dry and stable air masses confined between two inversion layers at an altitude of 1 to 5 km [3,4].
Dust particles play a significant role in the air quality and climate system [5,6]. Indeed, many epidemiological studies have already demonstrated their harmful effects on health. Respiratory and cardiovascular diseases are the most cited [5,[7][8][9][10]. Recently, Domínguez-Rodríguez et al. [11] showed that an inhalation of desert dust may result in a 2% increase in cardiovascular mortality risk. Regarding the climate, dust particles influence Earth's radiation budget directly by scattering, absorbing, and re-emitting radiation, and indirectly by impacting the atmospheric dynamics [12][13][14][15]. Recent research highlighted that African dust influences air temperature behavior in the atmospheric boundary layer [16,17]. Furthermore, dust particles may also inhibit the development of tropical cyclones [18,19].
In boreal summer, the Caribbean basin is frequently impacted by African dust outbreaks [20][21][22]. During this period, there is a continuous alternation between African Easterly Waves (AEWs) and dust plumes [23,24]. Recently, in June 2020, an extremely intense dust plume dubbed "Godzilla", which reached the Atlantic ocean [11,25,26], aroused the attention of the scientific community. This event is the biggest dust storm in half a century [26]. The dust plume was so large that it covered the Caribbean Sea and darkened the skies in several states of the United States [11]. Thus, it is crucial to study the characteristics of this event in order to better understand its impact on the Caribbean basin.
In this work, the granulometric distribution of the dusty events that have occurred over the last 20 years in the Caribbean islands of Puerto Rico (PR), Guadeloupe (GPE), and Barbados (BAR) (see Figure 1) was studied and compared to the extreme event of June 2020. To carry out this study, different types of data were used (ground-based, satellites, model, and soundings) on several sites in order to describe the dust events in terms of seasonality, intensity, frequency, and particle size behavior. All these measures are complementary to each other and essential to fully analyze the magnitude of an event. The remainder of this paper is organized as follows. Section 2 describes the databases used and the applied methods. Section 3 presents the achieved results and discuss them. Finally, a conclusion and an outlook for future studies are given in Section 4.

Satellites Images
Our analysis begins with a seasonal satellite overview of Aerosols Optical Depth (AOD) over the tropical North Atlantic. For this purpose, there is a wealth of product data derived from MODIS observations used to describe the characteristics of land, oceans, and atmosphere for local and global scales. In order to observe the seasonal behavior of mineral aerosols transported from African coast to the Caribbean basin, the AOD "Dark Target algorithm" satellite image processing is selected to improve AOD detection from MODIS on earth [27]: (i) spectral band at 550 nm, (ii) satellite mode Terra, and (iii) spatial resolution at 1 • . Thus, the AOD monthly mean is computed from the daily averaged data. The AOD images have been generated for the monthly periods from 2000 to 2020 (data available on the GIOVANNI (NASA) website: https://giovanni.gsfc.nasa.gov/giovanni/ (accessed on 20 January 2021).

Optical Data
The particle size distribution is studied using optical data reflecting the state of the atmospheric column. Caribbean (PR, GPE, and BAR) and African (Capo Verde and Dakar) sites are equipped with CIMEL sun photometer model CE-318-1 and supplied AOD measurements (stored as 15 min averages) in wavelengths of 1020-870-675-440-936-500-340-380 nm [28]. Sun photometer radiance measurements can be inverted to determine aerosol optical properties. This study is based on AERONET data at wavelength 440 nm: • AOD: measure of the extinction of the solar beam by particles in the atmosphere. • Ångström Exponent (AE) (440-870 nm): basically, referring to the particle size: the smaller the particle, the larger the value. • The Volume Particle Size Distribution (VPSD), corresponding to the daily concentration [dV(r)/dln(r) (µm 3 µm −2 )] retrieved for 22 logarithmically equidistant radius (r) in the range of sizes 0.05 ≤ r ≤ 15 µm. VPSD is an inversion product derivative of direct and diffuse radiation measurements [29,30].
To analyze the most data possible, only measurements at 1.5 quality level (cloud free) are selected. Table 1 presents the number of daily data used for this study. The uncertainty on data provided by AERONET is ± 0.01 for optical thickness measurements [31]. The surface concentration of particulate matter with an aerodynamic diameter of less than 10 and 2.5 µm, typically termed PM10 and PM2.5, respectively, is measured by the Thermo-scientific Tapered Element Oscillating Microbalance (TEOM) particle monitors 1400ab and 1400-FDMS (Filter Dynamics Measurement System). French territories (GPE) and the United States (PR) use the same instrument. Continuous measurement of PM10 and PM2.5 is reported to 15 min averages, the latter being used to compute hourly and then daily data. The measurement accuracy is ±1.5 µg/m 3 for 1 h and ±0.5 µg/m 3 for 24 h. The reliability of the measurements carried out in GPE is controlled by the Central Air Quality Monitoring Laboratory (LCSQA), and The Environmental Protection Agency (EPA) is in charge of controlling data from PR. Data from PM10 and PM2.5 are used to quantify air quality level, while characterizing the spatial and temporal dimension of natural or anthropogenic pollution events. Regarding BAR, surface dust concentration measurements are collected on a ground base. These data, not used in our work, are comparable to PM10 [24]. Table 2 presents the PM database considered for this study. A recent study demonstrated the consistency of data from several sites in the Caribbean area [21]. Indeed, dust events deteriorate air quality of the islands in the same way. Here, data from PR have been selected due to the regularity and completeness of PM10 and PM2.5 database between 2000 and 2020.

Visibility Data
Classically, visibility is defined as the greatest horizontal distance at which a black object of suitable dimensions and located near the ground can be seen and recognized when observed against a background scattering of hydrometeors (rain, snow, fog, and mist) or lithometeors (dust processes) [32]. In order to quantify the impact of the extreme event of June 2020 on visibility, data from the Meteorological Office of GPE (Météo France) are used. This measurement is made on the international airport of Pôle Caraïbes at Abymes (16.263 • N, −61.515 • E). Measured at hourly basis, these data are converted into daily data.

Back Trajectories
The Hybrid Single Particle Lagrangian Integrated Trajectory (HYSPLIT) version 4 (PC Windows-based HYSPLIT, Unregistered Version Download) [33,34] has been used to determine the origin of dusty air masses, analogously to studies carried out in Europe and the Caribbean [35][36][37][38]. Kalnay  This model was used to identify the origin of the air mass responsible for the June 2020 event.

Soundings
The thermodynamic structure of air is represented by the type of Skew-T diagram showing the description of humidity as a function of temperature, wind speed, and direction. Thus, the vertical profiles of temperature and humidity (mixing ratio) are shown at the following standard pressures: 1000, 925, 850, 700, 500, 400, and 300 hPa using atmospheric radiosonde measurements. The atmospheric soundings are measured daily at 00 UTC (8 p.m., local time) and 12 UTC (8 a.m., local time) at the meteorological stations of PR (code: 78526 TJSJ, San Juan) and BAR (code: 78954 TBPB, Grantley Adams). These soundings are available on Wyoming website (http://weather.uwyo.edu/upperair/sounding.html) (accessed on 15 December 2020). They are used to determine the characteristics of the SAL (thickness and wind) during the extreme event of June 2020.

Atmospheric Context of Caribbean Zone
The seasonal transport of mineral aerosol from African desert sources over the North Atlantic Ocean significantly impacts the atmospheric background of the Caribbean area [40]. Indeed, this zone is mainly characterized by the dust season and the non-dust season, highlighting two distinct thermodynamic and granulometric profiles [41]. In this work, the general atmospheric state of the Lesser Antilles was examined in comparison to dust episodes that occurred in the last twenty years, notably the extreme one in June 2020 [25]. Figure 2 presents an overview of the seasonal redistribution of mineral aerosols between the African coasts and Central America. The monthly behavior of mineral aerosols transport to the Caribbean islands over the last two decades shows that (i) the low dust season is from September to May (Figure 2a-c); nevertheless, this period also coincides to the high dust season of the Northern part of South America (Amazonia) [20,21], and (ii) the high dust season corresponds to June-July-August (Figure 2d).
At this stage, note that June is particularly marked by more frequent and more intense dusty events [16,42]. Figure 3 illustrates the magnitude of the phenomenon of June 2020 on the North Atlantic Ocean scale. Compared to Figure 2d, which gives the trend of the high dust season over 2 decades, we observe that the maximum AOD intensity corresponds to the area extending from Capo Verde to the middle of the Atlantic Ocean in Figure 3. It can be seen that the Caribbean Basin was strongly impacted by this exceptional event. According to Francis et al. [25], the largest AOD was recorded during this event at Capo Verde and over large swaths of the Atlantic Ocean due to the high content of dust suspended in the atmosphere.
Subsequently, the classical dust episodes compared to June 2020 were studied using different types of measurements located at several sites in the Caribbean.

AOD Versus AE
In this section, our analysis deals with optical data (AOD versus AE distribution) to determine the nature of the particles in the Caribbean zone. In the insular context, sea salt particles compose the main atmospheric component [43,44]. The latter is characterized by varied sized particles where the coarse mode dominates without significant impact on AOD associated with visibility [21]. However, during dust events, the AOD values increase while the AE decreases, highlighting the coarse mode of the particles [45,46]. On the other hand, anthropogenic pollution generates fine particles [47]. Previous studies established thresholds related to AOD versus AE values reflecting marine, dust, continental, and biomasses burning (BB) particles [21,29,30,48]: • AOD ≥ 0.2 combined with AE ≤ 0.6 related to dust particles. • AOD < 0.2 combined with AE < 2 related to marine particles (background atmosphere). • AOD ≥ 0.15 combined with AE > 0.6 related to continental/BB particles (anthropogenic pollution). Figure 4 shows the classification of particle types according to these known thresholds for the Caribbean zone. According to the results presented in Table 3, we notice that from 1996 to 2021, among the 9044 daily data collected, 64.89% are associated with the background atmosphere, 22.11% cases are related to dust episodes, and 3.47% to infrequent events of industrial pollution (continental/BB).  These results are consistent with previous studies [21]. This implies that the Caribbean area constitutes a privileged field of study, free from the impact of anthropogenic pollution due to its low industrial activity [49,50]. The main contributor of mineral pollution remains the seasonal episode of desert dust [20,21,51]. For June 2020, 68.42% of daily cases are associated to dust episodes. Later in this study, we will focus on the extreme event of June 2020.

Monthly Volume Particles Size Distribution (Vpsd) Analysis from 1996 to 2021
The size distribution of the particles leads to the characterization of their composition and properties. The predominance of the fine/coarse mode results in a concentration peak for a specific radius which testifies to their nature. As a reminder, the radius 2.24 µm is deeply associated with desert dust aerosols [21,46,52,53], while the radius 3.86 µm corresponds to sea salt aerosols [21,52].
For African sites (see Figure 5a), all months show the same behavior with a peak concentration at radius 2.24 µm, i.e., the desert dust signature. This is not surprising, as African dust sources (Mauritania, Sahara, and Bodélé) are active all year round, but their intensities fluctuate according to the seasons. Indeed, the maximum activity of the Sahara/Mauritania area is in winter, while the maximum activity of Bodélé is in summer. The cycle of activation related to dust sources has been highlighted in many studies [54][55][56][57][58][59][60][61], to mention a few. Over the 26 studied years, the maximum peak concentration corresponds to June and July with a concentration of 0.33 and 0.25 µm 3 µm −2 , respectively. The other months have an average concentration ranging from 0.13 to 0.20 µm 3 µm −2 . December has the minimum concentration with 0.1 µm 3 µm −2 . June's concentrations in 2020 stand out. Indeed, this latter is clearly above the seasonal value with 0.83 µm 3 µm −2 , i.e., the highest concentration. June 2020 strongly influences the overall average concentration for the months of June. For Caribbean sites (see Figure 5b), June and July show the highest peak of concentration at radius 2.24 µm as expected. May and August show intermediate concentration rates but are still dominated by desert dust particles. The other months (September to April) are mainly marked by marine particles (background atmosphere) with a peak of concentration at radius 3.86 µm. June and July correspond to the maximum peak with a concentration of 0.14 and 0.11 µm 3 µm −2 , respectively. The concentration in May and August is approximately 0.08 µm 3 µm −2 , while the other months have an average concentration ranging from 0.02 to 0.05 µm 3 µm −2 . In the Caribbean basin, only June 2020 stands out. Indeed, the concentration of June 2020 is approximately 1.5 times higher than the seasonal norm with 0.19 µm 3 µm −2 . It is important to emphasize that concentrations of the Caribbean are 4 times lower than concentrations observed on the African coasts. The decrease in volume particle size distribution is due to the dust particles deposition by dry or wet process during transport over the Atlantic. Overall, these results are in agreement with Velasco et al. [46] study. This analysis allowed to quantify the intensity of June 2020 event between Africa (source) and the Caribbean area (receptor).

Surface Dust Analysis
In addition to the optical data describing the column aerosols measurement, ground measurements were also analyzed. Thus, we examine the PM10 concentrations measured by the air quality networks of Caribbean islands reflecting the pollution phenomenon impact on a local scale. The Pearson correlation coefficient between AOD and PM10 in the West Indies is R = 0.8 [21,46], i.e., better than those found in Europe [62], the United States [63], and China [64,65]. The PM10 thresholds inform the population of the air quality level and warn them of critical thresholds for human health. PM2.5, being more harmful, is closely monitored. However, mineral dust aerosols are mainly associated to coarse particles, i.e., PM10. In the Caribbean Basin, desert dust events are detected from PM10 concentration superior or equal to 35 µg.m −3 , i.e., dusty cases [21]. The classically threshold used in literature to detect dust episodes and alert the population is 50 µg.m −3 . The latter corresponds to the European legislation (Directive 2008/50/EC) [66], which has proved to be unsuitable for the the Caribbean area [21].  From one year to another, one can notice a strong heterogeneity in the statistical parameters values, i.e., median, mean, standard deviation, skewness, and kurtosis. One of the reasons that may explain this variability is the activation of dust sources in Africa [67]. Indeed, many sources of dust in Africa become more active in summer [68]. These activations are related to the development and movement of AEWs in collaboration with extra-tropical disturbances. In addition, the variations of soil properties in these source areas (rainfall, vegetation cover, or land use) coupled with the climate processes that affect them will act to modulate transport to receptor sites in the Western Atlantic [69].  14 9 In Figure 6, the presence of outliers is not systematic in the box-plot for each year. June 2020 stands out with a value at 453 µg.m −3 the 23rd related to PM10 peak concentration of the so-called "Godzilla" event. This value is 9 times higher than the 50 µg.m −3 recommended by the World Health Organization over a 24 h period [70].

PM10 Statistical Analysis
In Table 4, June 2020 exhibits the highest values for mean, standard deviation, skewness, and kurtosis. This clearly shows the exceptional criterion of this month. During this month, PM10 daily concentrations are not high, except from the 20th to the 29th when [PM10] ≥ 50 µg.m −3 for 9 consecutive days. Apart from the intensity of a dust event, it is important to emphasize its duration. Taking the example of June 2012, which has the second highest average, we see that [PM10] ≥ 35 µg.m −3 for 21 days out of 30. However, among the 14 days when [PM10] ≥ 50 µg.m −3 , the longest event is 5 consecutive days between the 1st and the 5th. For large-scale phenomena such as African dust, it seems that the more intense the event is, the longer its duration will be.
In order to quantify the impact of June 2020 event on a meteorological parameter in the Caribbean, the visibility has been used. The Pearson correlation coefficient computed between PM10 and visibility is R = −0.8, i.e., strong anti-correlation. Figure 7 illustrates the behavior of PM10 concentrations and visibility data in GPE for June 2020.
During the first 18 days of the month there is no dusty event. Thus, PM10 average values are almost constant with small standard deviations highlighting the value of the background atmosphere. By contrast, due to the lack of turbidity, visibility average values and their standard deviations are high. From the 19th to the 24th, one can clearly see that the passage of the event has a significant impact on visibility. Due to the persistence and intensity of the event, the standard deviations values are, respectively, minimum for the visibility and maximum for PM10. On the 22nd, the higher value of PM10 (317 µg.m −3 ) corresponds to the weaker value of visibility (9 km). During the strong dusty event from 19th to 23th June 1994 in GPE, the visibility also decreased to 9 km [71]. Note that the June 2020 event occurred during the same period as the June 1994 event.
Between GPE and PR, there is a one day lag for the peak of the event at PR, however the intensity is greater at PR. For the lag, this is due to the fact that PR is further from the African coast than GPE. Regarding the intensity, the mechanisms explaining this difference will be analyzed in a future study.

Relationship between PM2.5/PM10 Ratio and VPSD
Before analyzing the relationship between the PM ratio and the VPSD, it is essential to first investigate PM2.5 and PM10 behavior. From 2000 to 2020, the Pearson correlation coefficient between PM2.5 and PM10 is R = 0.73. Figure 8 illustrates the monthly average concentrations for PM2.5 and PM10 highlighting a representative year over two decades. Overall, one can observe that PM2.5 and PM10 concentrations seem to follow the same temporal pattern. The highest PM2.5 value corresponds to the highest PM10 value (or vice versa). The impact of the dusty episodes in May and June clearly stands out. For May, PM2.5 and PM10 concentrations are, respectively, 8.8 and 33.7 µg.m −3 ; in June, they are, respectively, 10.3 and 48.2 µg.m −3 . It is relevant to emphasize that the standard deviations are weaker for the low dust season. Therefore, it can be reasonably considered that these two quantities are influenced by common factors. Nevertheless, we underline that from August to November, the PM2.5 − PM10 curves seem to diverge slightly. The origin of this phenomenon will be analyzed in a future study.
Analysis of PM2.5 proportion in PM10 (PM2.5/PM10 ratio) leads to improve the characterization of the dust events by giving it a health significance. Some studies carried out in the Mediterranean [72,73] and in Asia [74] highlighted fine particle emissions in highly industrialized areas. For the Caribbean zone, we deepened our analysis on PM2.5/PM10 by refining the following ratios: • A1: ratio < 0.2; • A2: 0.2 ≤ ratio < 0.4; • B: 0.4 ≤ ratio < 0.6; • C: 0.6 ≤ ratio < 0.8; • D: ratio ≥ 0.8. Figures 9-11 show three kinds of particle analysis: (i) PM2.5 versus PM10 sorted by previously defined ratio groups, (ii) the corresponding VPSD, and (iii) AE versus AOD for these ratios. This multiple approach allows us to determine the nature of the particles based on the known classification "VPSD" and "AE versus AOD", while evaluating the potential health risk linked to the proportion of PM2.5 in PM10. Indeed, the higher the PM2.5, the more harmful the health effects, especially for people with respiratory diseases and asthma [75,76]. Our multi-analysis shows two specific situations from 2000 to 2020: all cases (Figures 9a-11a) and dusty cases ([PM10] ≥ 35 µg.m −3 ; Figures 9b-11b).  For C and D (ratio ≥ 0.6), PM10 values are low contrary to PM2.5. However, due to the few cases retrieved over the study period, the achieved results are not significant.   The VPSD curves in Figure 10a depict the main dominance of desert dust, i.e., most of peaks of concentration are at radius 2.24 µm. Only the concentration levels stand out from each other. A1 and A2 show the highest values with, respectively, 0.080 and 0.060 µm 3 µm −2 . B and D correspond to approximately 0.045 µm 3 µm −2 . Then, C displays a more extended distribution with a low peak at 0.020 µm 3 µm −2 . Ultimately, B-C-D ratios proportionally have the largest part of PM2.5 in PM10 and also exhibit the lowest VPSD concentrations. This shows that the increase in VPSD concentration is more influenced by PM10 than by PM2.5. Compared to A1-A2 groups, the B-C-D group represents few cases. Note that a number of cases associated with the D ratio occur mainly outside the high dust season: September-November and February-May. Cases associated with the B ratio occur throughout the year at low frequencies. From our point of view, we think that B-C cases could be linked to the same type of phenomenon, i.e., a weak dust episode coupled with an influence of anthropogenic origin. Finally, most cases of the C ratio are intermittent and occurred during the same period: March-April-May 2017. This reflects a particular and non-seasonal event. According to the literature, this behavior could be explained by the smoke plumes which come up from Central Africa and the Amazon/Orinoco, mainly during the February to April season, just before the ITCZ rains move north [77].   [21]. In this analysis, only A1, A2, and B are concerned with pure dust episodes. The coefficient of determination between PM2.5 and PM10 remains high with R 2 = 0.7 and 0.8, respectively, for A1 and A2. Nevertheless, B displays R 2 = 0.3 with very few cases. As regards the VPSD graph (see Figure 10b), the peaks at 2.24 µm are more significant: 0.125, 0.100, and 0.060 µm 3 µm −2 , respectively, for A1, A2, and B. Figures 4 and 11a show that the case rates associated with dust, marine, and continental/BB aerosols are equivalent to, respectively, 22%, 65%, and 3%. The results from Figure 11b are presented in Table 5. Pure daily dust episodes account for more than 51.75% of cases, and dusty events with a marine dominance represent 46.36%. On the other hand, continental/BB cases are less significant with approximately 2%. As the atmospheric background is mainly dominated by marine aerosols, the results seem consistent. One can assume that the more the part of PM2.5 increases, the more important the marine influence is. Thus, A2 is more marked by the influence of marine aerosols than A1. This reinforces the postulate that dusty episodes are more significantly related to PM10 than to PM2.5.

Thermodynamic Analysis of The Atmosphere
Back trajectories from the sites located in San Juan (PR) and Grantley Adams (BAR) provide information on the origin of the dust-laden air masses that generated the extreme natural pollution episode of 23 June 2020. These two stations are representative of the general thermodynamic structure of the Caribbean area. From the HYSPLIT model (see Figure 12a), one can notice that the back trajectories come from Africa, and they indicate significant air uplift occurred 9 days before (the 14th at ≈6000 m of altitude) in the source area border between Mali and Algeria, i.e., the main source referenced in the bibliography [43,71,[78][79][80]. According to Cavalieri et al. [81] and Adams et al. [82], 6000 m is the maximum altitude corresponding to the mechanism for lifting particles. This result is consistent with Francis et al. [25] study using satellite imagery. In this study, the processes behind the lifting and transport of this intense dust storm have been identified. Due to an anomalous atmospheric circulation pattern in the mid-latitudes associated with a circum-global wavetrain, a subtropical high-pressure system to the west of the Saharan Heat Low developed [25]. This generated a continuous dust emissions from the Sahara the 14 to 19 June 2020.
The dust plumes are subsequently driven westward in high-altitude wind shears, thus crossing the Atlantic Ocean to the Caribbean Basin via the SAL. The latter presents thermodynamic characteristics associated with dry and hot air [3,41]. The soundings in Figure 12b,c illustrate the temperature and the humidity profile of a dusty atmosphere. On the right curve, it can be seen the sudden temperature inversion around 1600 m is associated with a strong drying out of the layer at the same altitude on the left curve. It is the base of the SAL. The second temperature inversion marks the top of the SAL, i.e., ≈5000 m. In other words, the thickness of the SAL is approximately 3400 m above the Caribbean area during this event. One can notice that the top of the SAL (≈5000 m) is much higher than what has already been reported in literature, i.e., ≈3400-3600 m [36,77,82]. In addition, the soundings also display the intensity and direction of the wind related to the altitudes of dust transport. The average value of wind speed measured in this layer is 15 m.s −1 with a northeast average direction. In Africa, Francis et al. [25] found a wind speed of 20 m.s −1 in the SAL. The difference between both values can be attributed to the distance between the sites. Nevertheless, due to the magnitude of this event, these values are greater than the 10 m.s −1 usually found in the SAL [36,71,77].

Conclusions
To conclude, this climatological analysis of the dusty events in the Caribbean area over the last 20 years highlighted the extreme intensity of the episode in June 2020 dubbed "Godzilla". In addition to the recent meteorological study made by Francis et al. [25], we propose a particulate approach of this event through data measured from the surface (PM10, PM2.5) and the atmospheric column (AOD, AE, and VPSD). The peculiarity of the Caribbean Basin is the fact that both approaches correlate, which makes it a privileged field of study. Thus, the thresholds reported in our recent work [21] have been confirmed and different types of particles are characterized: sea salt aerosols (background atmosphere), mineral aerosols (dust events), and anthropogenic particles (rare case). The determination of the thresholds relating to the size of the particles provides information on their nature and their properties.
At synoptic scale, the magnitude of June 2020 event was clearly observed using satellite imagery. With 453 and 317 µg.m −3 , respectively, measured at PR and GPE for PM10 concentrations, this event reached daily values never recorded in the last 20 years, even for volcanic eruption (150 µg.m −3 ). Due to this extreme event, the PM10 and the atmospheric column (AOD and VPSD) analysis exhibited their maximum values for June 2020. From a thermodynamic point of view, we highlighted the exceptional characteristics of the SAL in terms of thickness and wind speed. This episode had a strong impact on the environment, as is evident from the visibility factor reported during the event. Therefore, it was crucial to characterize this phenomenon in terms of the size, concentration, and nature of the particles.
To the best of the authors' knowledge, the influence of marine aerosols in dusty events has never been assessed before in the Caribbean Basin. In this study, this topic was analyzed in order to better understand their possible interaction with PM2.5 which are more harmful than PM10. Our results showed that the more the proportion of PM2.5 in PM10 increases, the more the influence of sea salt aerosols is significant. In order to more precisely quantify the proportion of PM2.5 in PM10, a transfer entropy analysis should be performed in a future study. Funding: The authors declare that they have not received any fund for the present paper. The paper is the sole work of the authors and is not a part/product of any project.

Institutional Review Board Statement: Not applicable.
Informed Consent Statement: Not applicable.