Atmospheric Conditions for Uplift and Dust Transport in the Latitudinal 10 ◦ North–20 ◦ North Band in Africa

: Desert aerosols suspended in the atmosphere are a very marked fact in West Africa with estimates of 400 to 1000 million tons produced annually and concentrations exceeding 50 µ g · m 3 in Burkina. In Bamako, the daily dust concentration can go up to reach 504 µ g/m 3 . The Sahara and the Sahel are recognized as the primary desert aerosol producing regions. Source areas continue to be discovered as the desert advances. Previous studies have mainly focused on the spatial and temporal variability of aerosols. The current question is: What makes an area a source of dust emission? Our study brings together all the climatic parameters of the 10–20 band, as well as the soil types and their characteristics; it reveals 4 soils characteristic of ﬁne sandy semi-arid soils in Chad. The Ouadaï plateau in Chad was identiﬁed as a source area for dust emissions. We noted for JFM (January, February, March) that the strongest wind intensities were located mainly towards Chad for average rmaximum temperatures around 34.7 ◦ C. The statistical study reveals a correlation of 66.8% between direct and indirect links between the climatic factors of the 10–20 band and the source area. The presence of vortexes throughout the year and a vertical wind proﬁle that is among the strongest in the 10–20 band, this gradient is strongly localized in the grid “10 ◦ North, 20 ◦ North and 20 ◦ East, 30 ◦ East” next to the Kapka massif. The study shows that the AEJ (African Easterly Jet) proﬁle, which is a strong wind, associated with the harmattan circulation, allows the transport of aerosols from Ouadaï to the West African coast. In Senegal, a signiﬁcant deposition was observed.


Introduction
The presence of desert aerosols suspended in the atmosphere is a very prominent fact in West Africa. The emission processes of these aerosols result from complex interactions between wind speed and soil roughness [1]. The Sahara and the Sahel are recognized as the first regions producing desert aerosols [2][3][4][5]. Several types of instruments have been used over time to quantify and identify the source areas [6][7][8][9]. Correlations between wind, temperature and precipitation during the African monsoon season have been established [10]. However, there are other parameters that can have a significant direct or indirect influence 2 of 19 on the lifting, transport and deposition of fine particles in the Sahelian zone. In the literature, the focus has been mainly on the spatial and temporal variability of aerosols. Authors such as [11], have summarized the source areas of dust emissions. A large majority of them are located in the 10 • -20 • North African band (hereafter 10-20 band). Source zones continue to be discovered with the advance of the desert and the climatic changes which result in a warming of the planet. We have studied all the climate parameters available in the 10-20 band in relation with this issue. These parameters include the volumetric water level of the soil which plays a role on evaporation as well as the type of soil, the evapotranspiration, which is a factor associated with the transpiration of the vegetation cover, and the geopotential at 1000 hPa, which is the gravitational potential at a point of the earth. These parameters associated with the vertical velocity, the vorticity at 1000 hPa and the winds at 10 m will be the catalysts of the wind erosion of the soils at the base of the uplift of the particles. The study also includes total precipitation and temperature at 2 m, which in turn will play an essential role on the soil volumetric water level and soil type.
The objective of this work is to make a climatic and statistical study of all the parameters having a direct or indirect link with the Mechanism of Uprising, Transport, and Deposition of Aerosols (MSTDA). The relevance of this study lies in the fact that beyond establishing a source zone of dust emission, we show what the factors are at the origin of this emission beyond the simple wind erosion of the soil.

Materials and Methods
Our study area is the scene of climatic instability, caused by the superposition of four types of climate. These are the Saharan, Sahelian, Sahelo-Sudanese, Sudano-Sahelian climates, as well as forest areas. This is the band of 10 • North and 20 • North latitudes (10-20 band) bounded by the African continent. The Sahelian zone is one of the largest and most important dust emission regions [11]. Emissions that can be related to parameters are evaporation, evapotranspiration, geopotential at 1000 hpa, total precipitation, soil type, temperature at 2 m, wind at 10 m, vertical velocity, soil volumetric water, vorticity at 1000 hpa and specific humidity. Parameters are directly extracted from the ERA5 platform, which is the fifth generation of ECMWF reanalysis (European Centre for Medium-Range Weather Forecasts) for global climate and weather for the last 4-7 decades, headquarters: Reading, UK. Currently, data are available from 1979 onwards. ERA5 (European Environment Agency) Copenhague, Danemark replaces the ERA-Interim reanalysis [12,13]. The data were regridded to a regular lat-lon grid of 0.25 degrees for the reanalysis and 0.5 degrees for the uncertainty estimation. There are four main subsets: hourly and monthly products, both on pressure levels (upper air fields) and single levels (atmospheric, ocean wave and land surface quantities).

Results
In order to understand the MTSDA in the Sahelian 10-20 band, the parameters directly and indirectly related to AUVI (ultraviolet index aerosol) and dust deposition have been studied in this section. The study of these parameters was conducted over a period from 1978 to 2008. Figure 1 shows the spatial distribution of wind intensity and direction. For JFM (January, February, March), the strongest wind intensities are located mainly towards Sudan, Ethiopia and Chad. The FIT (intertropical front) [14] is very easily recognizable in this figure. It is represented by the light band. This is consistent with the work of [15][16][17]. The equations explaining the atmospheric dynamics with the geostrophic wind allow a quantitative approach. In the case where the Rossby number (which is a dimensionless number used in fluid dynamics. It represents the relationship between the forces of inertia and the forces due to rotation which characterize the movement of a fluid in a rotating frame) is very small, much less than 1, we are close to a geostrophic equilibrium between the Coriolis and pressure forces (Equation (1)): where → v g is the geostrophic wind, ρ is the air density, f is the Coriolis parameter, p is the pressure at given altitude z,k is the orientation of z taken positive upward.
where ⃗ is the geostrophic wind, is the air density, is the Coriolis parameter, is the pressure at given altitude z, is the orientation of z taken positive upward.
The Coriolis force is orthogonal to the velocity vector, and opposite to the pressure force, the geostrophic wind is itself orthogonal to the pressure gradient and therefore parallel to the isobars. In the northern hemisphere, the low pressures are to the left of the wind (Buys Ballot's Law), to the right in the southern hemisphere. The geostrophic wind speed is proportional to the pressure (or z) gradient. When the Rossby number is large (i.e., small scale), the wind is accelerated from high to low pressure. The speed will be limited at the surface by friction.  The Coriolis force is orthogonal to the velocity vector, and opposite to the pressure force, the geostrophic wind is itself orthogonal to the pressure gradient and therefore parallel to the isobars. In the northern hemisphere, the low pressures are to the left of the wind (Buys Ballot's Law), to the right in the southern hemisphere. The geostrophic wind speed is proportional to the pressure (or z) gradient. When the Rossby number is large (i.e., small scale), the wind is accelerated from high to low pressure. The speed will be limited at the surface by friction. Figure 2 shows the temperature distribution at two meters from the ground in the 10-20 band. It can be seen that the temperatures range from 17 • C to 34.7 • C. From January to March, the average monthly temperature from 1979 to 2008 is around 17 • C. From April onwards, we note average maximum temperatures around 34.7 • C. These observations are close to those of the studies of [18,19]. Figure 2 shows the temperature distribution at two meters from the ground in the 10-20 band. It can be seen that the temperatures range from 17 °C to 34.7 °C. From January to March, the average monthly temperature from 1979 to 2008 is around 17 °C. From April onwards, we note average maximum temperatures around 34.7 °C. These observations are close to those of the studies of [18,19].
Air temperature is an input to most surface temperature models. Often air temperature measurements are available. Nevertheless, to make predictions, it is useful, even necessary, to have simple and reliable models, as these make it possible to limit the input data for the models [20][21][22].  Air temperature is an input to most surface temperature models. Often air temperature measurements are available. Nevertheless, to make predictions, it is useful, even necessary, to have simple and reliable models, as these make it possible to limit the input data for the models [20][21][22]. Figure 3 shows that there is almost no rainfall from January to April in the Sahelian zone. A timid start to the rains can be observed from March onwards [23]. The region gradually shifts from the continental trade wind regime to the monsoon regime, thus characterizing the real start of the rains over the Sahel Onset, as also shown in the work of [24][25][26].
Atmosphere 2022, 13, x FOR PEER REVIEW 5 of 20 Figure 3 shows that there is almost no rainfall from January to April in the Sahelian zone. A timid start to the rains can be observed from March onwards [23]. The region gradually shifts from the continental trade wind regime to the monsoon regime, thus characterizing the real start of the rains over the Sahel Onset, as also shown in the work of [24][25][26].   Figure 4 shows the vertical wind speed, also called vertical velocity. From January to March, there is a strong fluctuation of vertical wind gradient. This gradient is strongly localized in the 10-20 • North; 20-30 • East grid, which in our case covers part of the Bodélé and Banizoumbou depression. Towards the month of June, we note a strong attenuation of this gradient with a tendency of the unidirectional downward profile [27,28]. and Banizoumbou depression. Towards the month of June, we note a strong attenuation of this gradient with a tendency of the unidirectional downward profile [27,28].
The difference between the geostrophic thermal wind at pressure p1 and p2 can be written as: Here, pressure p2 is conventionally located at a higher altitude p1. In the Northern Hemisphere, the thermal wind is directed parallel to the isotherms, with the high temperatures to the right.  The difference between the geostrophic thermal wind at pressure p 1 and p 2 can be written as: Here, pressure p 2 is conventionally located at a higher altitude p 1 . In the Northern Hemisphere, the thermal wind is directed parallel to the isotherms, with the high temperatures to the right. Figure 5 shows the vortex distribution at 1000 hPa. In the 10-20 band, there is a significant presence of vortices, a clear decrease from March onwards. From January to March, a very visible pattern formed by a succession of vortices in the Sahelian band is observed. Figure 5 shows the vortex distribution at 1000 hPa. In the 10-20 band, there is a significant presence of vortices, a clear decrease from March onwards. From January to March, a very visible pattern formed by a succession of vortices in the Sahelian band is observed.
The vorticity equation is given by the following equation:  The vorticity equation is given by the following equation: Figures 6 and 7 show evaporation and evapotranspiration, respectively. According to these two figures, the Sahelian zone is subject to high evaporation from January to June, and at the same time, to very low potential evapotranspiration. Similar observations have been made in [29,30].
Atmosphere 2022, 13, x FOR PEER REVIEW 8 of 20 Figures 6 and 7 show evaporation and evapotranspiration, respectively. According to these two figures, the Sahelian zone is subject to high evaporation from January to June, and at the same time, to very low potential evapotranspiration. Similar observations have been made in [29,30].    Mathematical approach to geopotential: Hydrostatic equilibrium causes the pressure to always decrease with altitude. A locally high pressure must therefore correspond to a high altitude of the isobaric surfaces. The geopotential is a concept mainly used in meteorology. It represents a potential energy of gravity.
The pressure force acting on a surface S is normal to that surface and its norm is pS. For a parcel of air of volume δV = δxδyδz (Figure 8).
mosphere 2022, 13, x FOR PEER REVIEW 10 Mathematical approach to geopotential: Hydrostatic equilibrium causes the pressure to always decrease with altitude. A cally high pressure must therefore correspond to a high altitude of the isobaric surfa The geopotential is a concept mainly used in meteorology. It represents a potential en of gravity.
The pressure force acting on a surface is normal to that surface and its norm is For a parcel of air of volume = (Figure 8). The total pressure force in given by Equation (3): In Figure 9, the horizontal pressure force is given by Equation (4): . Equivalence between pressure and altitude differences: points A and B are at the s altitude, A and C at the same pressure. The pressure at B is therefore higher than at C. Figure 10 shows that in January, the geopotential at 1000 hPa remains relatively with values around 11,173.7 m 2 s −2 compared to the North African part, beyond 20° No From April to June, the lowest geopotential values of the continent are found in the 1 band, similar observations were made by [31]. The total pressure force in given by Equation (4): In Figure 9, the horizontal pressure force is given by Equation (5): Mathematical approach to geopotential: Hydrostatic equilibrium causes the pressure to always decrease with altitude. A cally high pressure must therefore correspond to a high altitude of the isobaric surfac The geopotential is a concept mainly used in meteorology. It represents a potential ener of gravity.
The pressure force acting on a surface is normal to that surface and its norm is For a parcel of air of volume = (Figure 8). The total pressure force in given by Equation (3): In Figure 9, the horizontal pressure force is given by Equation (4): Figure 9. Equivalence between pressure and altitude differences: points A and B are at the sa altitude, A and C at the same pressure. The pressure at B is therefore higher than at C. Figure 10 shows that in January, the geopotential at 1000 hPa remains relatively l with values around 11,173.7 m 2 s −2 compared to the North African part, beyond 20° Nor From April to June, the lowest geopotential values of the continent are found in the 10band, similar observations were made by [31]. . Equivalence between pressure and altitude differences: points A and B are at the same altitude, A and C at the same pressure. The pressure at B is therefore higher than at C. Figure 10 shows that in January, the geopotential at 1000 hPa remains relatively low with values around 11,173.7 m 2 s −2 compared to the North African part, beyond 20 • North. From April to June, the lowest geopotential values of the continent are found in the 10-20 band, similar observations were made by [31]. Figure 11 shows the amount of water contained in the soil and Table 1 gives the types of soils associated with each volumetric properties. From latitude 10 degrees North from January to June, it can be seen that the volumetric amount of water in the soil rotates around values below 0.  Figure 11 shows the amount of water contained in the soil and Table 1 gives the types of soils associated with each volumetric properties. From latitude 10 degrees North from  Table 2 summarizes soil types, Figures 12 and 13 show the spatial distribution of low vegetation cover and soil type, respectively. It is noted that the 10-20 band contains the highest low vegetation cover and a sandy, semi-arid to arid soil type by location. These observations were also made by [33,34]. In Figure 13, this parameter is the soil texture (or classification) used by the ECMWF Integrated Forecasting System (IFS) land surface scheme to predict soil water holding capacity in soil moisture and runoff calculations. It is derived from the root zone data (30-   Mathematical approach to soil volumetric water level: The volumetric soil water level is associated with the soil water holding. In the case of our study, we consider level 1. Equation (6) gives the mass water content W: M w : water mass in the soil M S : dry soil mass Volumetric water content θ of the soil is given by Equation (7): V w : volume of soil water V t : total soil volume Table 2 summarizes soil types, Figures 12 and 13 show the spatial distribution of low vegetation cover and soil type, respectively. It is noted that the 10-20 band contains the highest low vegetation cover and a sandy, semi-arid to arid soil type by location. These observations were also made by [33,34]. In Figure 13, this parameter is the soil texture (or classification) used by the ECMWF Integrated Forecasting System (IFS) land surface scheme to predict soil water holding capacity in soil moisture and runoff calculations. It is derived from the root zone data (30-100 cm below the surface) of the FAO/UNESCO Digital Soil Map of the World, DSMW (FAO, 2003). The seven soil types are shown in Table 2.  Tropical organic 0 Indicates a non-terrestrial point. This parameter does not vary over time.        In Figure 15, 66.8% of the information is held by the first three dimensions of our principal component analysis. Knowing that eigenvalues can make it possible to determine the number of principal axes, the difference between the other remaining eigenvalues being relatively small, the first two axes explain the relationships between the varia- In Figure 15, 66.8% of the information is held by the first three dimensions of our principal component analysis. Knowing that eigenvalues can make it possible to determine the number of principal axes, the difference between the other remaining eigenvalues being relatively small, the first two axes explain the relationships between the variables well. So with Figure 15, we have a confirmation of the good representativeness of our variables. In Figure 15, 66.8% of the information is held by the first three dimensions of o principal component analysis. Knowing that eigenvalues can make it possible to det mine the number of principal axes, the difference between the other remaining eigenv ues being relatively small, the first two axes explain the relationships between the var bles well. So with Figure 15, we have a confirmation of the good representativeness of o variables.   Figure 16 shows the contribution of relative vorticity, evapotranspiration, vertical velocity, specific humidity and dust deposition to the first dimension. This figure shows the average expected contribution. If the contribution of the variables is uniform, the expected threshold value would be 10%. For a given component, a variable with a contribution greater than this threshold is considered to have a significant contribution to the factorial axis. It can be seen that the deposit on each month, with the exception of December, contributes well to the formation of the analysis in principal component. This confirms the significance of the results given by the correlation circle.
Atmosphere 2022, 13, x FOR PEER REVIEW 16 of 20 Figure 16 shows the contribution of relative vorticity, evapotranspiration, vertical velocity, specific humidity and dust deposition to the first dimension. This figure shows the average expected contribution. If the contribution of the variables is uniform, the expected threshold value would be 10%. For a given component, a variable with a contribution greater than this threshold is considered to have a significant contribution to the factorial axis. It can be seen that the deposit on each month, with the exception of December, contributes well to the formation of the analysis in principal component. This confirms the significance of the results given by the correlation circle. With regard to the contribution of individuals, Figure 17 shows that grids 5, 2, 4 and 7 contribute to dimensions 1, 2 and 3, while grids 1, 3 and 6 contribute less significantly in these three dimensions.
As for Figure 17 showing the variables that contribute to the second dimension and are not involved in the first dimension, there is the wind at 10 meters and the volumetric soil water level. With regard to the contribution of individuals, Figure 17 shows that grids 5, 2, 4 and 7 contribute to dimensions 1, 2 and 3, while grids 1, 3 and 6 contribute less significantly in these three dimensions.
With regard to the contribution of individuals, Figure 17 shows that grids 5, 2, 4 and 7 contribute to dimensions 1, 2 and 3, while grids 1, 3 and 6 contribute less significantly in these three dimensions.
As for Figure 17 showing the variables that contribute to the second dimension and are not involved in the first dimension, there is the wind at 10 meters and the volumetric soil water level.  As for Figure 17 showing the variables that contribute to the second dimension and are not involved in the first dimension, there is the wind at 10 meters and the volumetric soil water level.
The climatic observation shows that the uplift took place mainly on and around the 20 • East-15 • North coordinate point in Chad, more precisely on the Ouaddaï Plateau, which is surrounded by Bodélé, the Djourab Erg, the Mourdi Depression, the Kapka Massif, the Guéra Massif and the Chad Basin, as shown in Figure 18. the principal component analysis on the different parameters reveals why the uplift takes place in this zone through the correlation established between the different climatic parameters. The climatic observation shows that the uplift took place mainly on and around the 20° East-15° North coordinate point in Chad, more precisely on the Ouaddaï Plateau, which is surrounded by Bodélé, the Djourab Erg, the Mourdi Depression, the Kapka Massif, the Guéra Massif and the Chad Basin, as shown in Figure 18. the principal component analysis on the different parameters reveals why the uplift takes place in this zone through the correlation established between the different climatic parameters.

Discussion
According to Figure 13 showing soil types in the 10-20 band, specifically in the Ouaddaï plateau, soil type indices range from 4 to 4.2. This characterizes the fine, semi-arid to arid soil types, which are caused by low rainfall; rainfall is only observed in July, August and September, and does not exceed 5692 millimeters. This is associated with low evapo-

Discussion
According to Figure 13 showing soil types in the 10-20 band, specifically in the Ouaddaï plateau, soil type indices range from 4 to 4.2. This characterizes the fine, semiarid to arid soil types, which are caused by low rainfall; rainfall is only observed in July, August and September, and does not exceed 5692 millimeters. This is associated with low evapotranspiration and a high index of low vegetation cover as shown in Figures 7 and 12.
A semi-arid to arid soil is also caused by a high evaporation of the volumetric water level of the soil which is observed in Figures 6 and 11. Similar results are found in the work of [35]. High evaporation values are associated with a rather high temperature index at 2 m throughout the year as can be seen in Figure 2. This remark is very close to those of [36] especially with the strong evaporation in the Sahel.
Ground-level winds at the Ouaddaï plateau are strongest in the 10-20 band from January to April and from September to December. This favors the saltation and sandblasting process. Since the geopotential is the gravitational potential of a locality at a given altitude, according to Figure 10, the Ouaddaï Plateau is subject to a relatively low geopotential throughout most of the year. This is associated with the presence of vortexes throughout the year and a vertical wind profile that is among the strongest in the 10-20 band just off the Kapka massif as shown in Figures 4 and 5. Thus, it can be said that particle uplift in this area is a very favorable phenomenon.
When uplift occurs on the Ouaddaï plateau, the heaviest particles with a weight of more than 20 µm fall back to the area. However, particles with a diameter of less than 20 µm manage to reach 850 hPa, at which point they enter the Harmatan circulation, which is a low-level north-northeast wind. Particles with a diameter less than or equal to 10 µm (PM10) reach the lower part of the Atmospheric Mid-Layer Jet AEJ (African Easterly Jet), which is a very strong wind, and the PM10 that reaches it is transported to the coastal country Senegal. The results obtained on the role of the wind in the uplift of particles concurred with those of [37].
The AEJ mid-layer winds weaken as they approach the coast. Aerosols are transported in the AEJ. As they approach the coast, they become heavier combined with the specific moisture brought by the westerly winds from the Atlantic Ocean and the winds from the Gulf of Guinea which are also loaded with moisture and reinforced by the Saint Helena High. As these aerosols become heavier, they lose altitude and fall back down to around 900 hPa. This is the atmospheric pressure level at which winds in Senegal become weak according to studies by [38]. Consequently, aerosols found there tend to settle. Figure 19 summarizes all the mechanisms studied.
Atmosphere 2022, 13, x FOR PEER REVIEW 18 of 20 band just off the Kapka massif as shown in Figures 4 and 5. Thus, it can be said that particle uplift in this area is a very favorable phenomenon. When uplift occurs on the Ouaddaï plateau, the heaviest particles with a weight of more than 20 µm fall back to the area. However, particles with a diameter of less than 20 µm manage to reach 850 hPa, at which point they enter the Harmatan circulation, which is a low-level north-northeast wind. Particles with a diameter less than or equal to 10 µm (PM10) reach the lower part of the Atmospheric Mid-Layer Jet AEJ (African Easterly Jet), which is a very strong wind, and the PM10 that reaches it is transported to the coastal country Senegal. The results obtained on the role of the wind in the uplift of particles concurred with those of [37].
The AEJ mid-layer winds weaken as they approach the coast. Aerosols are transported in the AEJ. As they approach the coast, they become heavier combined with the specific moisture brought by the westerly winds from the Atlantic Ocean and the winds from the Gulf of Guinea which are also loaded with moisture and reinforced by the Saint Helena High. As these aerosols become heavier, they lose altitude and fall back down to around 900 hPa. This is the atmospheric pressure level at which winds in Senegal become weak according to studies by [38]. Consequently, aerosols found there tend to settle. Figure 19 summarizes all the mechanisms studied. Figure 19. Schematic map summarizing the uplift process from the Ouaddaï plateau in Chad to Senegal.

Conclusions
The results obtained allowed the definition of the aerosol cycle from the Ouaddaï plateau to Senegal. The study showed that ground-level winds on the Ouaddaï Plateau favor the saltation and sandblasting process. The Ouaddaï Plateau is subject to a relatively low geopotential throughout most of the year. There is a presence of vortexes throughout Figure 19. Schematic map summarizing the uplift process from the Ouaddaï plateau in Chad to Senegal.

Conclusions
The results obtained allowed the definition of the aerosol cycle from the Ouaddaï plateau to Senegal. The study showed that ground-level winds on the Ouaddaï Plateau favor the saltation and sandblasting process. The Ouaddaï Plateau is subject to a relatively low geopotential throughout most of the year. There is a presence of vortexes throughout the year and a vertical wind profile that is among the strongest in the 10-20 band just off the Kapka massif. Particle uplift in this area is a phenomenon that is very conducive. Particles with a diameter of 10 micrometers or less (PM10) reach the lower part of the Atmospheric Mid-Layer Jet (AEJ), which is a very strong wind, and the PM10 that reaches it is transported to the coastal country Senegal. Institutional Review Board Statement: Ethical review and approval was waived for this study because it did not involve humans and animals.