# Infrared optical signature reveals the source-dependency and along-transport evolution of dust mineralogy as shown by laboratory study

Claudia Di Biagio, Jean-Francois Doussin, Mathieu Cazaunau, Edouard Pangui, Juan Cuesta, Pasquale Sellitto, Milagros Ródenas, Paola Formenti

<figure><img src="image_1.png" /></figure>

## To cite this version:

Claudia Di Biagio, Jean-Francois Doussin, Mathieu Cazaunau, Edouard Pangui, Juan Cuesta, et al. Infrared optical signature reveals the source-dependency and along-transport evolution of dust mineralogy as shown by laboratory study. Scientific Reports, 2023, 13 (1), pp.13252. 10.1038/s41598-023-39336-7. [hal-04192330](https://hal.science/hal-04192330v1)

HAL Id: hal-04192330

https://hal.science/hal-04192330v1

Submitted on 31 Aug 2023

HAL is a multi-disciplinary open access archive for the deposit and dissemination of scientific research documents, whether they are published or not. The documents may come from teaching and research institutions in France or abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est destinée au dépôt et à la diffusion de documents scientifiques de niveau recherche, publiés ou non, émanant des établissements d'enseignement et de recherche français ou étrangers, des laboratoires publics ou privés.

<figure><img src="image_2.png" /></figure>
<figure><img src="image_3.png" /></figure>

Check for updates

OPEN

# Infrared optical signature reveals the source-dependency and along-transport evolution of dust mineralogy as shown by laboratory study

Claudia Di Biagio<fnref n="1" />✉, Jean-François Doussin², Mathieu Cazaunau², Edouard Pangui², Juan Cuesta², Pasquale Sellitto²,³, Milagros Ródenas⁴ & Paola Formenti<fnref n="1" />

Advancing knowledge of the mineralogical composition of dust is key for understanding and predicting its climate and environmental impacts. The variability of dust mineralogy from one source to another and its evolution during atmospheric transport is not measured at large scale. In this study we use laboratory measurements to demonstrate that the extinction signature of suspended dust aerosols in the 740–1250 cm⁻<fnref n="1" /> atmospheric window can be used to derive dust mineralogy in terms of the main infrared-active minerals, namely quartz, clays, feldspars and calcite. Various spectral signatures in dust extinction enable to distinguish between multiple global sources with changing composition, whereas modifications of the dust extinction spectra with time inform on size-dependent particles mineralogy changes during transport. The present study confirms that spectral and hyperspectral infrared remote sensing observations offer great potential for elucidating the size-segregated mineralogy of airborne dust at regional and global scales.

Mineral dust is amongst the most abundant and widespread aerosol specie on Earth<fnref n="1" />⁻³. Dust aerosol affects the climate system in multiple ways, including direct effects on atmospheric shortwave and longwave radiative budget⁴⁻⁸, indirect effect on liquid and ice cloud formation and properties⁹⁻¹², contribution to biogeochemical cycles by acting as a source of nutrients for the ecosystems¹³,¹⁴, participation to atmospheric chemical reactions¹⁵,¹⁶, and contribution to air quality degradation and impacts on human health¹⁷,¹⁸. As demonstrated by extensive scientific work in the last decades, the strength and sign of these different effects depend on the dust mineralogical composition, i.e. the abundance, relative proportions and state of mixing of the different minerals composing the aerosols, which includes mainly silicates in the form of clays (kaolinite, illite, smectite, chlorite,...), quartz, and feldspars (orthose, albite...), carbonates (calcite, dolomite), sulfates (gypsum) and iron and titanium oxides¹⁹⁻²¹. Indeed, different minerals show different features in terms of their spectral absorption and scattering properties, capacity to act as cloud condensation or ice nuclei particles, chemical reactivity, or solubility²²⁻²⁸. The knowledge of the dust mineralogical composition is fundamental for assessing its broad role in the Earth's climate system and the environment.

Observations show that the mineralogy of dust aerosols is far from being homogeneous in the atmosphere³,¹⁹,²⁰,²⁹⁻³⁴. First, due to the diverse soil mineralogy of the different source areas, the composition of dust varies with the region of emission, and this is both at the global, regional and local scales¹⁹,³⁵⁻³⁷. The dust mineralogical composition also changes with size: quartz, feldspars, and calcium- rich species are generally more abundant in the coarse mode component, while clays and iron oxides rely mostly on the fine fraction below 2 µm²⁰,³⁸. As a consequence of its size-dependence, the dust mineralogy modifies during transport because of the progressive loss of coarse particles due to gravitational settling³³,³⁹,⁴⁰. Mixing of dust plumes with different origins and history may additionally occur in the atmosphere, further affecting the mineralogy of suspended dust.

<footnote n="1">¹Université Paris Cité and Univ Paris Est Creteil, CNRS, LISA, F-75013 Paris, France. ²Univ Paris Est Creteil and Université Paris Cité, CNRS, LISA, F-94010 Créteil, France. ³Istituto Nazionale di Geofisica e Vulcanologia, Osservatorio Etneo, Catania, Italy. ⁴EUPHORE Labs., Fundación CEAM, 46980 Paterna, Valencia, Spain. ⁵email: claudia.dibiagio@lisa.ipsl.fr</footnote>
As of today, available information on the dust mineralogical composition, its size - dependency and its spatio-temporal changes in the atmosphere remains still sparse and mostly limited to the specific conditions sampled during intensive field campaigns. The capability to get regional and global mapping of airborne dust mineralogy is still missing, a crucial limitation for properly developing and validating the representation of dust in Earth System Models, as well as for constraining its regional and global climate forcing<sup>7,26,41</sup>.

Due to their global coverage, frequency of observations, and multiyear duration, hyperspectral satellite measurements in various spectral ranges are currently used for characterizing desert dust abundance from daily to multi-year scale, in terms of horizontal and vertical distributions (for example<sup>42–46</sup>). Because the different minerals composing the fine and coarse fractions of dust show different spectral absorption signatures, these satellite observations can also be potentially used to detect the presence of diverse minerals and reconstruct their relative proportions in the dust aerosols<sup>28,47,48</sup>. Therefore, they offer the potential to provide highly valuable knowledge on size-resolved dust mineralogy and its spatio-temporal pattern. Ultraviolet-visible (UV-Vis, ~300–800 nm) and Visible to Short Wavelength Infrared (VSWIR, ~350–2500 nm) regions contain the absorption signatures of iron oxides (hematite, goethite), clays and carbonates, more abundant in the dust fine-fraction, whereas the longwave infrared region (LWIR, ~2.5–25 µm) is more sensitive to the presence of coarse-sized silicates and carbonates characterized by specific lattice vibrational-rotational transitions in this spectral region, including quartz and feldspars to which VSWIR radiation is not sensitive<sup>23,28,48,49</sup>.

Up to now, efforts aimed at quantitatively exploiting satellite observations to derive dust mineralogy from their optical signature have focused on the UV-Vis and VSWIR regions of the spectrum, therefore mostly looking at dust small-sized minerals<sup>50–52</sup>. Using the Deep Space Climate Observatory (DSCOVR) Earth Polychromatic Imaging Camera (EPIC), Go et al.<sup>50</sup> developed a retrieval algorithm to determine the total iron-oxide mass fractions in dust and its speciation in hematite and goethite minerals. Sanwlani and Das<sup>52</sup> use MODIS (Moderate Resolution Imaging Spectroradiometer) spectral optical depth observations at discrete wavelengths to reconstruct the main mineralogical features of the detected dust plumes. The EMIT mission (Earth Surface Mineral Dust Source Investigation)<sup>51</sup>, which integrated the International Space Station (ISS) in June 2022, will employ VSWIR reflectance observations to retrieve new data on the mineralogy of source soils (clays, iron oxides, sulfates, carbonates) from dust-emitting regions worldwide.

Hyperspectral observations in the LWIR, as those performed by IASI (Infrared Atmospheric Sounding Interferometer)<sup>53</sup> and the IASI-NG (Next Generation)<sup>54</sup> instruments, can complement information obtained from UV-Vis and VSWIR sensors by detecting coarse-sized mineral components, including silicates, in particular quartz, clays, and feldspars, and carbonates. In the infrared atmospheric window region located at 8–13.5 µm (740–1250 cm<sup>-1</sup>, i.e. the spectral region of the infrared spectrum with relatively low absorption by atmospheric gases) the LWIR dust signature arises from the overlapping contribution of quartz, clays and feldspars above 900 cm<sup>-1</sup>, and separate bands for calcite (876 cm<sup>-1</sup>) and quartz (800 and 777 cm<sup>-1</sup>) below 900 cm<sup>-1</sup><sup>123,28</sup>. Past literature based on laboratory transmission, reflectance, and emission measurements on pellet, thin film and soil samples has demonstrated that the LWIR spectra of mixtures of minerals in this spectral region can be reconstructed as the linear combination of the spectra of its single components, an approach named linear spectral mixing (LSM)<sup>48,55–57</sup>, which is used to investigate the surface composition of the Earth and Mars<sup>58,59</sup>. The LSM analysis of LWIR spectra of natural dust mixtures can be used to retrieve quantitatively its mineralogical composition. A first successful attempt of application of an LSM-like method to atmospheric IASI thermal infrared spectra has been conducted for retrieving information on the coarse dust composition for a case study over the Gobi desert<sup>60</sup>.

In the present study we explore the potential of LWIR spectral extinction signal of airborne mineral dust to inform about its mineralogy at the global-scale level. We use data of LWIR high-resolution extinction spectra, mineralogical composition and size distribution obtained from laboratory chamber experiments on suspended natural dust aerosols from nineteen source regions worldwide with contrasting mineralogy<sup>61</sup>. The dust aerosols, generated with a system that mimics the saltation and sandblasting of soils during wind erosion, are injected in the CESAM large simulation chamber<sup>62</sup> and aged for a time between 60 and 120 min under the effect of gravitational settling. The progressive loss of coarse particles during aging time allows to reproduce the temporal changes of dust size distribution and mass load as observed in the real atmosphere from emission to medium-to long-range transport<sup>61,63,64</sup>, comparable to what is observed for African dust from source areas towards the Mediterranean and the Atlantic<sup>20,32,39,40,65–67</sup>. The LSM analysis is applied to the dust extinction spectra measured during the experiments. A library of reference extinction spectra for single minerals composing dust including silicates (phyllosilicates (mica, serpentine/kaolinite, smectite, chlorite, and talc/pyrophyllite groups) and tectosilicates (feldspar and silica groups) families), carbonates, diatomite, and sulfates, constructed from data available in the literature for suspended particles and compressed pellets is used in the analysis (see Methods).

We take advantage of this unique dataset to test the capability of the LSM approach to provide consistent information on dust aerosol mineralogical composition, its size- and source-dependency, and its modifications due to size-selective coarse depletion. The bulk mineralogy and mass concentration derived from X-ray diffraction (XRD) and X-ray fluorescence (XRF) on dust particles collected on filters over the whole duration of each experiment is compared to the retrieved mineralogical composition and reconstructed mass concentration of the suspended dust from LSM analysis. Thanks to the specificity of the spectral features in the 740 to 1250 cm<sup>-1</sup> (8–13.5 µm) range, it is possible to reproduce the high-resolution LWIR extinction spectra of natural dust mixtures starting from a limited number of single mineral references. A quantitative estimate of the relative amounts of phyllosilicates, quartz, feldspars, carbonates, as well as minor components as sulfates or diatomite, can be derived in agreement with XRD analysis and literature observations. Based on their unique extinction signature, the mineralogical composition of samples of dust from different sources can be retrieved. Modifications of the dust extinction spectra with residence time in the chamber are shown to inform on modifications in clays partitioning and variations in feldspars content with evolving size distribution. Realistic estimation of the total
dust mass concentration and its temporal evolution is provided by LSM. The analysis performed in the present study demonstrates the possibility to retrieve consistent mineralogical signatures and mineral quantification from the measurements of high-resolution extinction spectra of dust in the LWIR spectral window and to follow the evolution of mineralogy with size-selective gravitational settling as occurring during large scale dust transport. These results support new ways of exploiting existing and future satellite measurements and suggest that LWIR remote sensing spectroscopy can help in significantly advance our knowledge of dust mineralogy at a global scale by mapping diversity of origin and properties.

## Results

**LSM fitting of experimental dust spectra.** Results of LSM analysis applied to chamber experimental data are illustrated (Fig. 1) for three experiments with dust aerosols from Tunisia, Australia, and Kuwait, showing different extinction intensities and spectral features. The comparison between measured and modelled spectra, as taken at the peak of the dust injection in CESAM and after 20 and 60 min of aging is reported for the three samples (Fig. 1a-c). The spectral extinction versus LSM fitting is also illustrated for all the remaining sixteen samples in this study (Supplementary Fig. 1). As seen in Fig. 1a-c, the difference between measured and modelled spectra (residual) is low, meaning that the model provides a very good fit of the measured spectra at the different times of the experiment. The Normalized Root Mean Square Error is low, typically between 1 and 5%. The modelled spectra (which is given by the sum of minerals contribution, i.e. the portion of the extinction spectrum reconstructed as sum of single mineral extinctions, plus a baseline offset term) very well reproduce band positions and intensities. The minerals contribution (in red in Fig. 1a-c) accounts for most of the modelled signal and reproduces the different absorption bands in the spectra. The polynomial baseline (in dotted black in Fig. 1a-c), which can assume both positive and negative values, shows a smoothed spectral variation and contributes to a fraction of the modelled signal that is usually less than 20% at 1000 cm⁻¹ but increases up to 50% towards low wavenumbers.

**Minerals abundance retrieved from LSM analysis.** The LSM procedure allows to resolve both overlapping and distinct mineral absorbing features and the contribution of different species to dust mass concentration is reconstructed as a function of time. Results are illustrated for samples from Tunisia, Australia and Kuwait (Fig. 1d-f; reporting absolute mass concentrations), and summarized in Table 1 (as percent mineral contributions) for the nineteen dust samples of this study. Globally, realistic mineral abundances are simulated, and sample-to-sample differences are identified, in overall agreement with literature-based knowledge.

Phyllosilicates represent between 36 and 91% in mass for all samples at the peak of injection in CESAM, values which are consistent with field investigations of African and Asian dusts¹⁹,²⁰,³¹,⁶⁸. Illite dominates the clay speciation for North African and Chinese samples, in line with observations²⁹,³⁰,³⁶, whereas montmorillonite is dominant for Mauritania, Ethiopia, Kuwait, Patagonia, Namib-2 and Australia samples. Kaolinite contributes significantly to clay speciation for Northern African and Sahelian samples, but it is almost absent in all the other regions. Quartz is lower than 5% for most samples but reaches up to 21% content for Sahelian, Gobi, Patagonia and Australia samples at injection peak. Quartz-rich samples in Sahel and Asia are also reported in the literature¹⁹,³⁰,³¹. Feldspars are present in variable amounts ranging between 1 and 26% in all samples, with the highest values for Niger and Bodélé. The geographical variability in dust feldspars content is in agreement with soil composition showing the highest feldspars content in the Sahel⁶⁹. However, available field observations have shown a much lower contribution (<2%) of feldspars from Niger and Bodélé sources¹⁹. Calcite is found to be the highest in Tunisia, Morocco, Saudi Arabia, Atacama and Namib-1 samples (33 to 46% at the peak of injection in CESAM) and less than 13% for the other samples, with exception of Niger, Mali, Patagonia, Namib-2 and Australia for which no calcite is detected. Geographical distribution of calcite between analyzed samples well reflects the global distribution of this mineral in soil samples⁶⁹ and field aerosol observations²⁰,²⁹⁻³¹. Dolomite has a negligible contribution to reconstructed mineralogy, in line with previous mineralogical analysis in Northern Africa¹⁹,²⁰, while anhydrite is between 0 and 8% for all samples at the peak of injection. Field observations of sulfates are limited to crystalline gypsum, shown to contribute a few percent to the dust mass¹⁹,²⁰. Diatomaceous Earth has no contribution to all samples, with the exception of Bodélé (6% at peak) and a small contribution to Ethiopia (2%). The Bodélé depression is a diatomite-rich area, and aerosols from this zone contain resuspension of detrital diatoms³⁵. We are capable therefore to identify this unique feature of Bodélé features from our dataset and analysis.

**Temporal evolution of the mineralogical composition from LSM analysis.** Results reported in Fig. 1d-f and Table 1 illustrate the evolution of LSM reconstructed mineralogy with a lifetime in the chamber up to more than 2 h aging that corresponds to a gravitational processing equivalent to that potentially occurring over medium-to-long-range transport. The effective diameter ($D_{eff}$) of dust, also reported in Fig. 1d-f and Table 1, varies in the range 3.1-10.9 µm for fresh conditions (injection peak) to 1.1-4.0 µm for aged dust, values which reproduce the size-selective changes observed in the field within 1 to 5 days of transport³⁹,⁶⁵,⁶⁶,⁷⁰. Results indicate that the dust mineralogy retrieved from LSM changes from the beginning to the end of each experiment. However, no systematic tendencies are identified for all samples. The content of phyllosilicates, calcite, quartz and feldspars is alternatively increasing or decreasing for the various samples for fresh and aged conditions. A closer look at the data however evidences some consistent features. In terms of phyllosilicates speciation, it can be noticed that when kaolinite is present, it tends to increase systematically from fresh to aged conditions, with a maximum of 32% increase for Libya (that is the % contribution of kaolinite to total dust mass goes from 17 to 49% from fresh to aged conditions). As a consequence, kaolinite, which is never the dominant clay specie for fresh dust as discussed in the previous paragraph, becomes the dominant clay specie for aged Libya, Niger
<figure><img src="image_4.png" /><figcaption>**Figure 1.** Illustration of Linear Spectral Mixing (LSM) analysis and comparison to online measurements. (a, b, c) Measured and modelled extinction spectra (expressed as absorbance in log10 convention) for (a) Tunisia, (b) Australia, and (c) Kuwait samples at the peak of the dust injection in the CESAM chamber and after 20 and 60 min aging time. The modelled spectra from LSM analysis is the sum of the single minerals contributions (Sum Minerals, plotted in red) and the baseline offset term (Baseline, plotted in black dotted) (see Methods). The Residual (difference between measured and modelled signal) is also shown (plotted in grey). Examples of baseline and residuals are reported at the bottom of each plot in reference only to the spectra at 60 min. The horizontal black line represents the 0 value. (d, e, f) Temporal evolution of the total mass concentration of dust aerosols and minerals contributions obtained from LSM modelling and comparison with the total mass concentration of dust in CESAM obtained from online size distribution data. Mass concentrations obtained both considering full size distribution (size total) and size cut at 5 µm diameter, i.e. integrating only diameters below 5 µm, (size cut 5 µm) are reported. The size total data are plotted against the right y-axis while the size cut 5 µm data are plotted against the left y-axis. The effective diameter ($D_{eff}$) calculated from online size distribution data at the beginning and at the end of the experiments is also indicated. Note that the peak of dust injection is reached more rapidly for Kuwait sample (panel f) compared to Tunisia and Australia (panels d-e) which causes a sharper increase of the mass concentration curve at the initial stages of the experiment.</figcaption></figure>

and Bodélé dust, in good agreement with field observations¹⁹. This could potentially suggest that the fresh dust in CESAM is representative of the first stages of dust emissions in atmosphere, a condition which has been rarely sampled in the field²⁰,⁷¹. Montmorillonite shows on the contrary a tendency to decrease (down to 30% variation) going from fresh to aged conditions, whereas illite increases with aging for most samples and almost systematically when kaolinite is not present. While the speciation between different clay species seems to be affected by gravitational settling, the total amount of phyllosilicates is observed to either remain almost constant with time
<table><thead><tr><th rowspan="2">Sample name</th><th colspan="4">Phyllosilicates (%)</th><th rowspan="2">Quartz (%)</th><th rowspan="2">Feldspars (%)</th><th rowspan="2">Calcite (%)</th><th rowspan="2">Dolomite (%)</th><th rowspan="2">Anhydrite (%)</th><th rowspan="2">Diatomaceous Earth (%)</th><th rowspan="2">D<sub>eff</sub> (μm)</th></tr><tr><th>Total</th><th>Illite (%)</th><th>Kaolinite (%)</th><th>Montmorillonite (%)</th></tr></thead><tbody><tr><td>Tunisia - fresh</td><td>46 ± 1</td><td>34 ± 2</td><td>1 ± 1</td><td>9 ± 2</td><td>5 ± 1</td><td>10 ± 1</td><td>35 ± 1</td><td>0 ± 0</td><td>2 ± 0</td><td>0 ± 0</td><td>3.5 ± 0.9</td></tr><tr><td>Tunisia - aged (62)</td><td>41 ± 2</td><td>24 ± 3</td><td>0 ± 0</td><td>13 ± 1</td><td>5 ± 0</td><td>13 ± 1</td><td>31 ± 1</td><td>0 ± 0</td><td>3 ± 0</td><td>0 ± 0</td><td>2.5 ± 0.6</td></tr><tr><td>Morocco - fresh</td><td>48 ± 0</td><td>34 ± 0</td><td>5 ± 1</td><td>10 ± 1</td><td>5 ± 0</td><td>2 ± 1</td><td>41 ± 1</td><td>0 ± 0</td><td>4 ± 0</td><td>0 ± 0</td><td>4.5 ± 1.1</td></tr><tr><td>Morocco - aged (100)</td><td>48 ± 3</td><td>36 ± 4</td><td>11 ± 2</td><td>2 ± 1</td><td>4 ± 1</td><td>1 ± 1</td><td>39 ± 4</td><td>0 ± 0</td><td>8 ± 1</td><td>0 ± 0</td><td>2.6 ± 0.7</td></tr><tr><td>Libya - fresh</td><td>83 ± 1</td><td>43 ± 4</td><td>17 ± 4</td><td>24 ± 7</td><td>2 ± 0</td><td>6 ± 1</td><td>6 ± 1</td><td>0 ± 0</td><td>3 ± 0</td><td>0 ± 0</td><td>9.5 ± 2.4</td></tr><tr><td>Libya - aged (184)</td><td>84 ± 5</td><td>34 ± 8</td><td>49 ± 5</td><td>1 ± 1</td><td>2 ± 1</td><td>6 ± 6</td><td>4 ± 3</td><td>0 ± 0</td><td>4 ± 2</td><td>0 ± 0</td><td>1.9 ± 0.5</td></tr><tr><td>Algeria - fresh</td><td>82 ± 1</td><td>27 ± 2</td><td>5 ± 1</td><td>50 ± 3</td><td>4 ± 0</td><td>12 ± 1</td><td>1 ± 0</td><td>0 ± 0</td><td>1 ± 0</td><td>0 ± 0</td><td>8.5 ± 2.1</td></tr><tr><td>Algeria - aged (120)</td><td>82 ± 3</td><td>41 ± 8</td><td>19 ± 2</td><td>22 ± 7</td><td>6 ± 2</td><td>10 ± 3</td><td>0 ± 0</td><td>1 ± 1</td><td>3 ± 1</td><td>0 ± 0</td><td>3.0 ± 0.8</td></tr><tr><td>Mauritania - fresh</td><td>81 ± 2</td><td>23 ± 4</td><td>0 ± 0</td><td>58 ± 6</td><td>3 ± 1</td><td>5 ± 2</td><td>11 ± 0</td><td>0 ± 0</td><td>0 ± 0</td><td>0 ± 0</td><td>7.0 ± 1.8</td></tr><tr><td>Mauritania - aged (82)</td><td>88 ± 4</td><td>32 ± 5</td><td>0 ± 0</td><td>56 ± 5</td><td>8 ± 2</td><td>2 ± 3</td><td>2 ± 3</td><td>0 ± 0</td><td>0 ± 0</td><td>0 ± 0</td><td>2.8 ± 0.7</td></tr><tr><td>Niger - fresh</td><td>56 ± 2</td><td>43 ± 2</td><td>12 ± 1</td><td>1 ± 1</td><td>22 ± 1</td><td>21 ± 3</td><td>0 ± 0</td><td>0 ± 0</td><td>2 ± 0</td><td>0 ± 0</td><td>5.1 ± 1.3</td></tr><tr><td>Niger - aged (138)</td><td>67 ± 4</td><td>27 ± 3</td><td>40 ± 1</td><td>0 ± 0</td><td>22 ± 2</td><td>11 ± 6</td><td>0 ± 0</td><td>0 ± 0</td><td>1 ± 1</td><td>0 ± 0</td><td>1.1 ± 0.3</td></tr><tr><td>Mali - fresh</td><td>75 ± 3</td><td>33 ± 4</td><td>7 ± 2</td><td>35 ± 6</td><td>7 ± 2</td><td>14 ± 3</td><td>0 ± 0</td><td>0 ± 0</td><td>5 ± 1</td><td>0 ± 0</td><td>4.1 ± 1.0</td></tr><tr><td>Mali - aged (56)</td><td>80 ± 10</td><td>36 ± 9</td><td>16 ± 4</td><td>27 ± 9</td><td>11 ± 5</td><td>3 ± 6</td><td>0 ± 0</td><td>0 ± 0</td><td>7 ± 1</td><td>0 ± 0</td><td>3.3 ± 0.8</td></tr><tr><td>Bodélé - fresh</td><td>44 ± 8</td><td>31 ± 4</td><td>4 ± 3</td><td>7 ± 4</td><td>17 ± 3</td><td>26 ± 4</td><td>0 ± 0</td><td>0 ± 0</td><td>8 ± 2</td><td>6 ± 1</td><td>8.4 ± 2.1</td></tr><tr><td>Bodélé - aged (138)</td><td>57 ± 4</td><td>12 ± 6</td><td>27 ± 4</td><td>9 ± 5</td><td>19 ± 12</td><td>12 ± 3</td><td>4 ± 4</td><td>0 ± 0</td><td>0 ± 0</td><td>8 ± 5</td><td>3.3 ± 0.8</td></tr><tr><td>Ethiopia - fresh</td><td>76 ± 1</td><td>28 ± 2</td><td>0 ± 0</td><td>48 ± 2</td><td>9 ± 0</td><td>9 ± 1</td><td>4 ± 1</td><td>0 ± 0</td><td>0 ± 0</td><td>2 ± 0</td><td>8.4 ± 2.1</td></tr><tr><td>Ethiopia - aged (132)</td><td>75 ± 9</td><td>22 ± 9</td><td>0 ± 0</td><td>53 ± 6</td><td>6 ± 5</td><td>7 ± 10</td><td>9 ± 2</td><td>0 ± 0</td><td>0 ± 0</td><td>4 ± 1</td><td>3.3 ± 0.8</td></tr><tr><td>Saudi Arabia - fresh</td><td>56 ± 1</td><td>41 ± 2</td><td>3 ± 0</td><td>12 ± 3</td><td>1 ± 0</td><td>9 ± 1</td><td>33 ± 0</td><td>0 ± 0</td><td>0 ± 0</td><td>0 ± 0</td><td>8.5 ± 2.1</td></tr><tr><td>Saudi Arabia - aged (132)</td><td>44 ± 1</td><td>33 ± 3</td><td>8 ± 2</td><td>3 ± 1</td><td>0 ± 0</td><td>9 ± 1</td><td>44 ± 1</td><td>0 ± 0</td><td>0 ± 0</td><td>0 ± 0</td><td>3.3 ± 0.8</td></tr><tr><td>Kuwait - fresh</td><td>75 ± 2</td><td>34 ± 5</td><td>0 ± 0</td><td>40 ± 7</td><td>1 ± 1</td><td>8 ± 2</td><td>13 ± 1</td><td>0 ± 0</td><td>5 ± 0</td><td>0 ± 0</td><td>5.4 ± 1.4</td></tr><tr><td>Kuwait - aged (60)</td><td>70 ± 7</td><td>42 ± 9</td><td>0 ± 0</td><td>28 ± 5</td><td>3 ± 2</td><td>9 ± 6</td><td>11 ± 3</td><td>0 ± 0</td><td>6 ± 1</td><td>0 ± 0</td><td>3.1 ± 0.8</td></tr><tr><td>Gobi - fresh</td><td>54 ± 1</td><td>46 ± 2</td><td>0 ± 0</td><td>7 ± 1</td><td>15 ± 0</td><td>13 ± 2</td><td>12 ± 0</td><td>0 ± 0</td><td>6 ± 0</td><td>0 ± 0</td><td>3.1 ± 0.8</td></tr><tr><td>Gobi - aged (78)</td><td>61 ± 1</td><td>52 ± 1</td><td>0 ± 0</td><td>8 ± 1</td><td>15 ± 0</td><td>2 ± 1</td><td>14 ± 0</td><td>0 ± 0</td><td>9 ± 0</td><td>0 ± 0</td><td>1.5 ± 0.4</td></tr><tr><td>Taklimakan - fresh</td><td>65 ± 1</td><td>42 ± 1</td><td>0 ± 0</td><td>0 ± 0</td><td>4 ± 0</td><td>15 ± 1</td><td>13 ± 2</td><td>0 ± 0</td><td>3 ± 0</td><td>0 ± 0</td><td>7.3 ± 1.8</td></tr><tr><td>Taklimakan - aged (124)</td><td>72 ± 4</td><td>64 ± 5</td><td>0 ± 0</td><td>0 ± 0</td><td>3 ± 0</td><td>5 ± 4</td><td>16 ± 2</td><td>0 ± 0</td><td>4 ± 1</td><td>0 ± 0</td><td>2.8 ± 0.7</td></tr><tr><td>Arizona - fresh</td><td>72 ± 1</td><td>27 ± 8</td><td>0 ± 0</td><td>30 ± 4</td><td>4 ± 1</td><td>12 ± 0</td><td>7 ± 1</td><td>0 ± 0</td><td>4 ± 0</td><td>0 ± 0</td><td>8.5 ± 2.1</td></tr><tr><td>Arizona - aged (72)</td><td>70 ± 4</td><td>46 ± 7</td><td>0 ± 0</td><td>19 ± 4</td><td>3 ± 1</td><td>16 ± 8</td><td>7 ± 3</td><td>0 ± 0</td><td>4 ± 2</td><td>0 ± 0</td><td>3.9 ± 1.0</td></tr><tr><td>Atacama - fresh</td><td>57 ± 1</td><td>29 ± 5</td><td>0 ± 0</td><td>28 ± 6</td><td>1 ± 0</td><td>7 ± 1</td><td>34 ± 1</td><td>0 ± 0</td><td>2 ± 0</td><td>0 ± 0</td><td>6.9 ± 1.7</td></tr><tr><td>Atacama - aged (128)</td><td>59 ± 3</td><td>49 ± 4</td><td>5 ± 1</td><td>5 ± 2</td><td>2 ± 1</td><td>5 ± 2</td><td>32 ± 3</td><td>0 ± 0</td><td>3 ± 1</td><td>0 ± 0</td><td>2.7 ± 0.7</td></tr><tr><td>Patagonia - fresh</td><td>55 ± 1</td><td>16 ± 3</td><td>0 ± 0</td><td>39 ± 2</td><td>20 ± 1</td><td>17 ± 1</td><td>0 ± 0</td><td>1 ± 0</td><td>8 ± 1</td><td>0 ± 0</td><td>10.9 ± 2.7</td></tr><tr><td>Patagonia - aged (64)</td><td>59 ± 1</td><td>14 ± 5</td><td>0 ± 0</td><td>45 ± 4</td><td>23 ± 1</td><td>11 ± 2</td><td>0 ± 0</td><td>0 ± 0</td><td>7 ± 1</td><td>0 ± 0</td><td>3.6 ± 0.9</td></tr><tr><td>Namib-1 - fresh</td><td>36 ± 0</td><td>23 ± 1</td><td>0 ± 0</td><td>13 ± 1</td><td>0 ± 0</td><td>10 ± 0</td><td>46 ± 1</td><td>0 ± 0</td><td>8 ± 0</td><td>0 ± 0</td><td>9.4 ± 2.4</td></tr><tr><td>Namib-1 - aged (82)</td><td>36 ± 2</td><td>25 ± 2</td><td>0 ± 0</td><td>10 ± 1</td><td>0 ± 0</td><td>11 ± 2</td><td>43 ± 1</td><td>0 ± 0</td><td>11 ± 0</td><td>0 ± 0</td><td>4.0 ± 1.0</td></tr><tr><td>Namib-2 - fresh</td><td>91 ± 0</td><td>9 ± 2</td><td>0 ± 0</td><td>82 ± 2</td><td>0 ± 0</td><td>9 ± 0</td><td>0 ± 0</td><td>0 ± 0</td><td>0 ± 0</td><td>0 ± 0</td><td>9.4 ± 2.4</td></tr><tr><td>Namib-2 - aged (124)</td><td>84 ± 5</td><td>8 ± 5</td><td>0 ± 0</td><td>77 ± 2</td><td>0 ± 0</td><td>16 ± 5</td><td>0 ± 0</td><td>0 ± 0</td><td>0 ± 0</td><td>0 ± 0</td><td>1.9 ± 0.5</td></tr><tr><td>Australia - fresh</td><td>76 ± 1</td><td>12 ± 4</td><td>3 ± 1</td><td>62 ± 5</td><td>8 ± 1</td><td>14 ± 1</td><td>0 ± 0</td><td>0 ± 0</td><td>2 ± 0</td><td>0 ± 0</td><td>10.4 ± 2.6</td></tr><tr><td>Australia - aged (122)</td><td>76 ± 6</td><td>23 ± 5</td><td>0 ± 0</td><td>54 ± 3</td><td>6 ± 1</td><td>12 ± 8</td><td>0 ± 0</td><td>0 ± 0</td><td>6 ± 2</td><td>0 ± 0</td><td>2.8 ± 0.7</td></tr></tbody></table>
**Table 1.** Mineralogy for fresh and aged dust samples retrieved from linear spectral mixing (LSM) analysis. Mineralogy retrieved for each sample at the peak of injection in the CESAM simulation chamber (labelled as fresh) and at the end of the experiment (labelled as aged). Percent mass fraction for each mineral is derived by using as total mass concentration the one derived from LSM analysis. For aged conditions the number of minutes of aging in the chamber is also indicated in parenthesis in the first column. Fresh and aged values are the average (± standard deviation) of five consecutive spectra. The total phyllosilicates mass includes the contribution of illite, kaolinite, montmorillonite, chlorite, chamosite, serpentine and talc. The illite, kaolinite and montmorillonite contributions to phyllosilicates are reported in the table. The total feldspars includes albite, oligoclase and orthoclase. The effective diameter D<sub>eff</sub> calculated from online size distribution measurements during experiments is also reported. Its uncertainty is estimated at 25%<sup>61,63</sup>.

(8 cases), or to slightly decrease (4 cases) or increase (7 cases) showing however maximum variation within 12%. One consistent feature observed is that the increase/decrease of phyllosilicates is opposite to feldspars. On the contrary, quartz content remains mostly unchanged between fresh and aged conditions with variations lower than 5% which do not co-occur with those of phyllosilicates or other minerals. Calcite behaves similarly, with variations going up to 10% in both directions and not clearly correlated to other mineral changes.

**Comparison of mineralogy retrieved from LSM and measured by X-ray diffraction (XRD).** The comparison between LSM and XRD derived experiment-averaged mineralogy for the nineteen analyzed samples in this study (Fig. 2 and Supplementary Fig. 3) confirms that the main features of dust mineralogy are well reproduced by LSM. Phyllosilicates are in good agreement, both in terms of absolute values (in the range 40–90% content) and sample-to-sample variability when comparing the two techniques. For calcite and feldspars, the LSM also well reproduces the sample-to-sample variability as obtained in XRD, but tends to overestimate the absolute content of these minerals. Calcite from LSM analysis, in particular, shows the highest values for Northern African (Tunisia, Morocco), Saudi Arabia, Atacama, and Namib-1 samples. Quartz is identified in almost all samples, in line with XRD observations, with the exception of Namib-1 and Namib-2 samples for which no quartz is retrieved from LSM. These samples show amongst the lowest quartz content from XRD analyses as well. Contrarily to calcite and feldspars, quartz is significantly underestimated in LSM compared to XRD, with the only exception of Mauritania and Ethiopia.

**Comparison of total mass concentration retrieved from LSM and measured by size distribution and X-ray fluorescence (XRF).** The total mass of dust aerosols reconstructed from LSM (sum up of the minerals contributions; Fig. 1d-f) increases as the starting of the dust injection in CESAM, maximizes the peak of injection (reaching up to tenths of mg m⁻³) and reduces with residence time in the chamber (down to values less than 1 mg m⁻³ after 120 min aging), following the intensity changes in measured spectra. The lowest detected mass concentration for single minerals is around 5 µg m⁻³, a value which can be assumed as detection limit in the present analysis. The temporal evolution of the total mass concentration from LSM compares well with the behavior derived from size distribution measurements from particle sizers connected online to the chamber (Fig. 1d-f), suggesting that the dynamic of the dust population is properly captured by LSM inversion. Nevertheless, the absolute value is up to a factor 10 lower that the total mass from size distribution data (as observed for all samples). The size/LSM mass ratio is related to the fraction of coarse particles in the dust and decreases for reducing the coarse component with aging in the chamber. It should be noted that the mass concentration from size distribution is highly uncertain because of the low concentrations of coarse particles and size-selective particle losses along sampling tubes<sup>61</sup>. Cutting the size distribution at 5 µm in diameter (Fig. 1d-f) would significantly improve the comparison, suggesting that the infrared observations are more sensitive to particles in this size range.

In order to have a further reference for the comparison of the mass concentration, we consider the one calculated from XRF elemental composition during CESAM experiments<sup>64</sup>. This estimate, based on mass determination of the main elements composing dust<sup>72</sup>, provides a robust evaluation, in excellent agreement with gravimetric measurements<sup>19,64</sup>. The comparison between LSM and XRF mass concentrations (Fig. 3) supports a very good agreement between the two techniques. The linear correlation between the two variables is excellent, as indicated by the correlation coefficient (R<sup>2</sup>=0.90; 0.98) and the reduced chi-square (χ<sup>2</sup><sub>red</sub>=1.7; 0.3) (the two set of values are obtained by considering or not three outlier points in the dataset, see Fig. 3). The slope of 0.92 ± 0.12 (with the three outlier points included) and 0.69 ± 0.08 (outliers excluded) indicates an average 8 to 31% difference when fitting LSM and XRF data, which suggests a moderate but systematic underestimation of the mass concentration by LSM.

Discussion

In this study, we perform a deconvolution analysis of the measured infrared extinction spectra of suspended natural dust aerosols to separate the contributions from its single constituent minerals and use these contributions to reconstruct the particles’ mineralogical composition and their spatio-temporal variability. Nineteen samples with contrasting mineralogy originated from worldwide sources are considered.

The results of the LSM inversion (mineral relative abundances, sample-to-sample variability) are in good agreement with both literature-based knowledge of dust mineralogical features and concurrent XRD measurements. Nevertheless, differences are found when comparing absolute values, in particular for quartz, calcite and
<figure><img src="image_5.png" /><figcaption>**Figure 2.** Comparison of mineralogy retrieved from LSM and measured by X-Ray Diffraction (XRD). Data represent the comparison (percent content in mass) for (a) phyllosilicates, (b) calcite, (c) feldspars, and (d) quartz. The XRD data are measured on dust samples collected on filters and represent mineralogy integrated over the duration of each experiment⁶¹. Uncertainty on XRD data is obtained from XRD single mineral calibration and total mass estimation⁶¹. For comparison purposes, and since only the crystalline phase is visible to XRD, the % content of minerals from LSM is rescaled removing from the total mass the contribution of non-crystalline minerals such as anhydrite and diatomaceous earth, and also that of ammonium sulfate. The LSM data are averaged over the same interval of XRD filter sampling. The uncertainty on LSM average values is calculated as the quadratic combination of mass concentration uncertainty (±27%, obtained with the error propagation formula considering ±19% uncertainty on both single mineral and total dust mass from LSM analysis) and the standard deviation over the average period.</figcaption></figure>

feldspars between LSM and XRD. The XRD and LSM are looking at different matter properties, i.e., light-diffraction linked to crystalline structure in XRD versus chemical bonds leading to rotational and vibrational
<figure><img src="image_6.png" /><figcaption>**Figure 3.** Comparison of mass concentration retrieved from LSM and measured by X-Ray Fluorescence (XRF). The XRF data are measured on dust samples collected on filters and represent mass concentration integrated over the duration of each experiment⁶⁴. The LSM data are averaged over the same interval of XRF filter sampling, and the XRF concentrations are corrected for particle losses along sampling lines⁶¹,⁶⁴. Data points in blue evidence three outliers in the dataset, corresponding to Saudi Arabia, Namib-1 and Patagonia dust. The uncertainty on XRF mass concentration takes into account loss correction uncertainties as detailed in Supplementary Methods. The uncertainty on LSM is the quadratic combination of LSM total mass concentration uncertainty (± 19%) and the standard deviation over the average period. The results of the linear fit are reported including also statistical indicators of the goodness of fit (correlation coefficient, R², and reduced chi square, χ²<sub>red</sub>). Two linear fits are reported: Fit1 (blue) considering all data points, and Fit 2 (red) removing the three outlier points in blue from the dataset. The grey shaded area represents the range of Fit 1 uncertainty.</figcaption></figure>

transitions in LWIR, and differences should be expected, as discussed in previous literature⁴⁸,⁴⁹,⁵⁷. The tendencies in LSM-XRD differences in our study are in line with those found for sedimentary rocks investigations⁵⁷.

The total mass concentration reconstructed from LSM is consistent with concurrent estimates from validated compositional-based analyses, which support the potential to recover both single minerals and total atmospheric dust mass loadings from measured LWIR optical signatures. However, despite an excellent linear correlation between LSM and XRF data, a potential systematic underestimation is identified from the analysis.

The LSM approach, i.e. its capacity to provide accurate mineral abundances, mineral relative proportions and total mass concentration, is dependent on the realism of the single mineral spectra in terms of the shape, intensity and position of absorption bands. Possible differences in size, morphology, and chemical state between reference and real minerals can affect the analysis. To this aim, and because of the well-known artefacts of the pellet technique<sup>23,73</sup>, it is important to have spectra for suspended minerals instead of compressed thin layer samples. Despite a number of reference spectra for suspended minerals is nowadays available for the main species composing dust, such as those used in the present analysis<sup>74-77</sup>, the ensemble remains limited. Also, we notice that most of reference spectra in the literature, including those used in this study, are obtained for submicron particles, while in nature most of these minerals extend up to several tenths of microns. In this regard, our analysis shows that the currently available library of reference minerals permits to well reproduce absorption bands, however the quality of the fit is not always perfect, and a non-negligible baseline contribution is needed to properly reproduce the spectra, especially at low wavenumbers. This could suggest that some mineral bands are not perfectly represented, and the fact that this occurs at low wavenumbers, where most of the interaction should come from coarse-sized particles, suggests a possible misrepresentation of submicron reference spectra compared to natural conditions. When considering the ensemble of reference spectra in the present analysis, it emerges that the dataset that more frequently contributes to LSM reconstruction and dominates the retrieved dust mass concentration is the one from Mogili et al.<sup>77</sup>, which is also the only one for which single mineral sizes extend up to super-micron diameters. The D<sub>eff</sub> estimated for illite, kaolinite, montmorillonite, quartz and calcite in Mogili et al.<sup>77</sup> is between 1.0 and 2.9 µm, which is at the lower bound but comparable to the values reported for our dust samples (Table 1), and significantly larger than the D<sub>eff</sub> for the same minerals in additional references considered in our analysis, which instead remain below 0.6 µm<sup>74,75</sup> (Supplementary Table 3). These results support the critical impact of the size distribution in reference spectra on the quality of the retrieval and reconstructed mineralogy. In agreement with previous studies<sup>49</sup>, our work suggests therefore that one limiting point for LSM application is the representativeness of reference spectra for single minerals, a dataset that should be consolidated and extended. Revising (or building up a new) database of single mineral mass extinction spectra that would include both the fine and the coarse mineral size distribution is needed. Spectra should be acquired on suspended particles and at high-resolution (2 cm<sup>-1</sup> or better) in order to well record the characteristic features of the minerals which are key for their spectral speciation.
Despite limitations discussed in the previous paragraphs, the present analysis clearly proves that it is possible to retrieve the mineralogical composition of dust aerosols from worldwide sources and follow its changes during transport due to gravitational settling based on their LWIR extinction signature. Results from LSM deconvolution do not evidence any systematic trend in the observed dust mineralogy changes with time, with some of the samples showing modifications in clays partitioning and variations in feldspars content with time, and others not. In general terms, the fact that mineralogy does not evolve in a unique way for the different samples in our analysis suggests that the size-dependent composition of dust is possibly not the same between different sources, therefore leading to different temporal evolutions of the mineralogy due to gravitational settling. As a matter of fact, size-segregated compositional analyses of dust remains limited to a few case studies<sup>20,32</sup> and low-level knowledge exists. Analyses conducted for Northern African dust in Morocco and Cape Verde<sup>20,32</sup> show that the dust mineralogy varies with size, but the most important variations are below 0.5 and above 20 µm, while the composition is more uniform in between these limits. This kind of size-segregated mineralogy could explain why in some case we do not detect significant changes for changing size, knowing that our investigated size range is located mostly in the 0.5–20 µm limit.

The results of this study, show that the LSM analysis applied to infrared extinction spectra can be used to retrieve global-scale features of dust mineralogy, therefore pave the way to a new perspective of spectroscopic investigation of dust aerosols. In particular, the outcomes of this study are relevant for advancing the exploitation of past, current and next-generation hyperspectral satellite missions, as IASI and IASI-NG, in order to complement the retrievals in the UV or VSWIR, such as those from the new EMIT mission<sup>51</sup>. The increasing capabilities of current and future satellite instruments, in combination with new analyses approaches such as those presented in this paper, can allow to go further in the characterization of dust aerosols with remote sensing, thus providing new capabilities for quantifying their broad impacts on climate, air quality, and the environment. Combining UV, VSWIR and LWIR satellite measurements is expected to finally provide access to a comprehensive global scale size-resolved mineralogical information. Note that while the analysis provided here focuses on the mid-infrared region, the same considerations can be potentially extended to the far-infrared, and be of relevance also for the future FORUM satellite mission (Far-Infrared Outgoing Radiation Understanding and Monitoring), that will provide observations of atmospheric radiation up to 100 µm<sup>78</sup>, where spectral signatures of dust minerals can be found<sup>28,79-81</sup>. In addition, it is of relevance for ground-based infrared remote sensing observations, for example, but not limited to, those provided by NDACC (Network for the Detection of Atmospheric Composition Changes). Note that the lower limit of dust signal for which the LSM approach has been successfully applied in the present study is about 0.005-0.01 at 10 µm (as absorbance, a.u., that is 60-130 Mm<sup>-1</sup> extinction coefficient accounting for the chamber optical pathlength). For 1 km deep uniform dust layer these extinctions would correspond to an aerosol optical depth of 0.06-0.1, a value that can be considered as the minimum to obtain the required signal-to-noise ratio to retrieve mineralogical information from a measured signature.

Using remote sensing observations to map dust mineralogy can provide new knowledge on the composition of dust laden air masses as they are emitted and transported across the globe, including areas not yet studied up to now because of difficult access (e.g. remote desert and oceanic areas), but also can provide insight on new potential source areas which are expected to progressively emerge in consequence of climate change-induced desertification. Indeed, because the LSM approach is sensitive to both crystalline, non-crystalline and amorphous compounds, it can reveal of high relevance for the study of high latitude and volcanic dust aerosols, composed of a significant amorphous and non-crystalline fraction<sup>82,83</sup>.

Global-scale mineralogical data can provide fundamental new input for improving Earth System Models performances to reproduce the dust atmospheric cycle. While soil mineralogical datasets are more and more integrated into model schemes<sup>41,84</sup>, including in perspectives those from the EMIT mission<sup>51</sup>, new global observational data are fundamental for either model assimilation and validation exercises, in particular for a better evaluation of dust climatic-relevant properties. Finally, retrieving mineralogical information on a large scale may also provide further understanding of the mechanisms behind dust variability and source identification, and can be also used for tracking of nutrients inputs to oceans and the ecosystems, fundamental for biogeochemical forcing attribution<sup>13</sup>.

As a concluding remark, the analysis presented in this study is based on laboratory extinction spectra measured on a spectral range and resolution comparable to those that can be detected from satellite hyperspectral sensors. It must be mentioned, however, that, compared to laboratory measurements, aerosol atmospheric observations from satellite are affected by the signal-to-noise ratio of current-generation infrared radiation measurements that can partly hide the aerosol signature. In the case of nadir-viewing instruments, possibly inhomogeneous vertical distribution of aerosols and the limited vertical sensitivity of the observations can also be a limiting factor in isolating the aerosol spectral signature. Additionally, for all types of real satellite observations, interference of infrared-absorbing gaseous species, clouds and the surface is also to be considered, from a full-radiative-transfer perspective. In this regard, studies focusing on the application of the LSM approach to real LWIR remote sensing observations are strongly needed.

## Methods

**Experimental data of infrared extinction spectra and physiochemical properties of natural suspended dust aerosols.** Experimental dust data analysed in this work are obtained in the 4.2 m³ stainless-steel CESAM chamber (French acronym for Experimental Multiphasic Atmospheric Simulation Chamber)<sup>62</sup> and fully described in<sup>61,63,64,85</sup>. The objective of the experiments is to investigate the spectral optical properties of mineral dust aerosols and their global-scale variability, going from UV-Vis to LWIR (0.37–16 µm, 625–27,027cm<sup>-1</sup>). To this aim a set of 19 natural soil samples with contrasted mineralogy and being representatives of the diversity of the surface composition at the global scale are selected based on concurrent information
from a global soil mineralogy database⁶⁹. The identified samples cover the four continents and both hemispheres and include Northern Africa (Algeria, Libya, Mauritania, Morocco, Tunisia), Sahel (Chad, Mali, Niger), Eastern Africa and Middle East (Ethiopia, Kuwait, Saudi Arabia), Eastern Asia (China), North America (Arizona), South America (Chile, Argentina), Southern Africa (Namibia), and Australia. Their list and origin are provided in Supplementary Table 1.

Mineral dust aerosols of these samples are generated by mechanical shaking of the natural parent soils, a procedure able to mimic the sandblasting of soils due to saltation during wind erosion as occurring in arid and semi-arid regions. This technique produces dust aerosol with realistic mineralogical composition and size distribution with respect to atmospheric conditions⁶¹,⁶³,⁸⁵. As discussed in Di Biagio et al.⁶¹, since the same protocol is applied for generating dust, differences in the size distribution and mass concentration for the various samples are related to differences in soil characteristics and its attitude to generate aerosols, in particular the coarse fraction. The dust aerosol produced by the mechanical shaking is injected into the CESAM volume and left in suspension for a time varying between 60 and 120 min while measuring the evolution of its physico-chemical and spectral optical properties by a suite of instruments both measuring in situ and connected online to the chamber, as detailed in Supplementary Table 2. The only aging process occurring during the 120 min dust suspension in CESAM is gravitational settling, intended to reproduce the progressive loss of coarse particles during medium-to-long-range transport of plumes in the real atmosphere. As discussed by Di Biagio et al.⁶¹, the dust number size distribution, measured by Scanning Mobility Particle Sizer (SMPS) and Optical Particle Counters (OPCs), extends up to tenths of µm at the peak of the injection and evolves rapidly during the experiments: particularly the coarse mode above 5 µm diameter rapidly decreases, and the observed dynamic well reproduces the size-selective changes observed in the field within 1 to 5 days of transport³⁹,⁷⁰. The mass concentration of dust aerosols in CESAM (derived from number size distribution assuming a density of 2.5 g cm⁻³) reaches up to tenths of mg m⁻³ at the peak of the injection, in line with observations during dust storms close to source regions²⁰, and reduces to values less than 1 to 2 mg m⁻³ after 120 min aging in CESAM (due to the combined effect of gravitational deposition and dilution due to instrumental sampling), comparable to what measured after several days of transport from African sources to the Mediterranean and the Atlantic⁶⁵,⁶⁶. These observations indicate therefore that a 2 h aging experiment in CESAM reproduces the temporal changes of size distribution and mass load of dust aerosols as observed in the real atmosphere from emission to medium-to-long-range transport. As a proxy of the size distribution, the effective diameter (D<sub>eff</sub>) defined as D<sub>eff</sub> = (Σ<sub>i</sub> πD<sub>i</sub><sup>3</sup>dN<sub>i</sub>) / (Σ<sub>i</sub> πD<sub>i</sub><sup>2</sup>dN<sub>i</sub>) is calculated for the different experiments. All chamber experiments are conducted at ambient laboratory temperature and dry conditions (relative humidity < 2%). Two experiments per soil type are usually conducted in order to test repeatability. The chamber is evacuated and kept at a pressure of 3·10⁻⁴ hPa between each experiment, and manually cleaned between experiments with different soils to avoid carry-over contaminations. Background concentrations of aerosols in the chamber are less than 2.0 µg m⁻³, which is about a factor of 1000 below operating conditions.

Concurrently with the size distribution and mass concentration, the evolution of the LWIR extinction spectrum of suspended dust aerosols is monitored online along the full duration of each experiment by means of an in situ Fourier Transform Infrared spectrometer (FTIR) interfaced with a multi-pass cell⁶²,⁸⁷ providing a total optical path length of (182.5 ± 0.5) m within the chamber. The dust extinction spectrum is retrieved from FTIR-measured absorbance between 625 and 5000 cm⁻¹ (2–16 µm wavelength range) at 2 cm⁻¹ resolution over 2 min intervals. For each experiment, the FTIR background spectrum is acquired immediately before dust injection in CESAM. Note that the FTIR spectra are degraded at 8 cm⁻¹ resolution in the analysis to match reference spectra resolution (Supplementary Methods).

The bulk mineralogical composition of the dust aerosols in the size fraction < 10.6 µm diameter is obtained by X-ray diffraction (XRD) on samples collected on filters over the whole duration of each experiment. Analysis of XRD spectra provided the mass proportion of clays, quartz, calcite, dolomite, gypsum, and feldspars in the aerosols. Additionally, X-ray fluorescence (XRF) is performed on the filter samples to retrieve the dust elemental composition, then used to reconstruct the dust mass concentration during experiments⁶⁴,⁷².

**LSM analysis of experimental extinction spectra.** The LWIR extinction spectra measured in the CESAM chamber for natural dust aerosols mixtures are mathematically deconvoluted to represent them as the LSM of its single constituent minerals⁴⁸,⁵⁶,⁵⁷. Under the hypothesis of external mixing, it is possible to write the spectral extinction of dust as a function of time as:

$$ \text{Ext}_{\text{Mod}}(\lambda, t) = \sum_{i=0}^{n_{\text{minerals}}} C_i(t) \text{Ext}_i(\lambda) + \text{Ext}_b(\lambda, t) \quad (1) $$

$$ \text{Residual}(\lambda, t) = \text{Ext}_{\text{Meas}}(\lambda, t) - \text{Ext}_{\text{Mod}}(\lambda, t) \quad (2) $$

where C<sub>i</sub>(t) are the coefficients representing the contribution of the mineral i to extinction spectra at the time t of the experiment, Ext<sub>i</sub>(λ) is the reference extinction spectra for mineral i, Ext<sub>b</sub>(λ,t) is the offset term (polynomial baseline) of the equation, and Residual(λ,t) is the difference between measured, Ext<sub>meas</sub>(λ,t), and modelled, Ext<sub>mod</sub>(λ,t), spectra.

Using the estimated coefficients C<sub>i</sub>(t), the mass concentration (MC, µg m⁻³) contribution of each mineral to the dust mixture as a function of time can be calculated as:
$$MC_i(t) = C_i(t) \cdot MC_{i,ref} \cdot \frac{x_{i,ref}}{x_{CESAM}} \quad (3)$$

where $MC_{i,ref}$ is the mass concentration for the measured reference, and $x_{i,ref}$ and $x_{CESAM}$ are the pathlengths used for transmission spectroscopic measurements in the reference studies and in CESAM, respectively. The uncertainty on the derived $MC_i$ is estimated to be 19% (90% confidence interval) (Supplementary Methods).

To determine the coefficients $C_i(t)$ and the $Ext_b(\lambda,t)$ in Eq. (1) an automatic least-squares fitting procedure is applied. The algorithm takes as input the measured extinction spectra and reference mineral spectra and provides the $C_i(t)$ and $Ext_b(\lambda,t)$ values that allows the best overlap between experimental and modelled data. The procedure calculates the baseline locally. Only a constraint for positive values is set to the calculation of $C_i(t)$ in the retrieval, despite in some isolated cases specific minerals have to be removed from the analysis since their inclusion provides unrealistic results. Calculations are performed at the same 2 min resolution of the FTIR spectra acquisition and restricted in the 740-1475 cm⁻¹ range (6.8 to 13.5 µm), that is in the region where the most intense and spectrally varying dust signatures are found²,²⁸. Note that above 1250 cm⁻¹ (<8 µm) the spectral signature of main dust minerals weakens considerably (Supplementary Fig. 2) with the only exception of carbonates and also that in this spectral region CESAM extinction spectra show additional noise due to possible residual contaminations by water vapor lines due to trace levels of H₂O during experiments. Results are shown and discussed in the range 740 to 1250 cm⁻¹ (8-13.5 µm, corresponding to the infrared atmospheric window) in the paper.

The Root Mean Square Error (RMSE) of the residuals, calculated as $\sqrt{\frac{\sum_{\lambda} (Ext_{mod}(\lambda) - Ext_{meas}(\lambda))^2}{N}}$, with N the number of data points, is used to provide quality assessment of the LSM fitting. The Normalised Root Mean Square Error (NRMSE) is the ratio of the RMSE by the average $Ext_{mod}(\lambda)$.

As two experiments per soil type were conducted, the LSM was repeated twice for each sample. The comparison of the results in the two cases demonstrates the repeatability of the whole procedure, both in retrieved minerals percent contributions and in temporal dynamics. However, we want to stress that the LSM analysis is conducted independently at each time step without any priori information. This means that reproducibility of the results and the robustness of the spectral analysis is verified at each time step retrieval for each experiment.

**Library of reference extinction spectra.** A library of reference extinction spectra for single minerals composing dust including silicates (phyllosilicates (mica, serpentine/kaolinite, smectite, chlorite, and talc/pyrophyllite groups) and tectosilicates (feldspar and silica groups) families), carbonates, and sulfates was constructed based on data available in the literature. For library constitution, we privileged datasets corresponding to measurements on suspended aerosol particles instead of compressed pellets, which are better representative of natural airborne aerosol samples and whose spectra are not affected by well-known artefacts including shift of absorption bands' peak and modification of absorption band depth and profile⁷³. Several studies fulfilled this requirement, allowing to collect data for major clay species (illite, kaolinite, montmorillonite), quartz, carbonates (calcite, dolomite), feldspars (albite, oligoclase), sulfates (anhydrite) and diatomaceous Earth (amorphous non-crystallized silica material)⁷⁴⁻⁷⁷. For other minerals, for which literature data on suspended particles are not available, we referred to classical pellet measurements, including references for five phyllosilicates (chamosite, chlorite, montmorillonite, serpentine, talc)⁷³ and tectosilicates (orthoclase; Caltech mineral spectroscopy server). All reference mineral spectra were acquired for submicron-sized particles generally limiting to D<2 µm, with the only exception of Mogili et al.⁷⁷ that extended the range of measured diameters up to 7.5 µm. The spectral resolution was 8 cm⁻¹ for the majority of the datasets. A summary of reference spectra used and conditions for their acquisition are provided in Supplementary Table 3 and detailed in Supplementary Methods. Supplementary Fig. 2 shows the reference spectra for different minerals used in this study. Note that for expected dominant minerals in dust (illite, kaolinite, montmorillonite, quartz, and calcite) various reference spectra from different literature studies are included in the library in order to evaluate the impact of different datasets in the process. Conversely, iron oxides are discarded, because of their relatively flat spectra which we tested to have a nil contribution to the deconvolution analysis of absorption peaks while potentially biasing the estimate of the baseline.

In addition to dust-composing minerals, we also included in the library the spectra of H₂O and CO₂, to account for possible contribution to the dust measured spectra. Also, ammonium sulfate aerosol spectra, presenting absorption bands at 7.0 and 8.9 µm⁸⁸, is included as control non-dust specie to evaluate the deconvolution process. Control ammonium sulfate aerosol is obtained to have a negligible contribution (typically <2%) to the total dust reconstructed mass for all samples.

**Data availability**

The simulation chamber experiments' data (LWIR spectra from FTIR, and total mass concentrations and effective diameters from size distribution) that support the findings of this study are available through the Database of Atmospheric Simulation Chamber Studies (DASCS) of the EUROCHAMP Data Centre (https://data.eurochamp.org/data-access/chamber-experiments/) with the identifiers listed in Supplementary Table 4. Reference spectra for single minerals are published in⁷⁴⁻⁷⁷. Data from Dorschner et al.⁷³ are obtained from the Heidelberg-Jena-St.Petersburg-Database of Optical Constants (HJPDOC, https://www2.mpia-hd.mpg.de/HJPDOC/index.php; last access March 2023). Data from the Caltech Mineral Spectroscopy server are retrieved at http://minerals.gps.caltech.edu/files/Infrared_MIR/Minerals_From_JK/Index.htm; last access March 2023. The HITRAN (high-resolution transmission molecular absorption database) data are downloaded at https://hitran.org/about/ (last access March 2023). The XRD mineralogical data in Fig. 2 are from Di Biagio et al.⁶¹ (in their supporting information). The XRF mass concentration data in Fig. 3 are published in Caponi et al.⁶⁴.
## Code availability

The Linear Spectral Mixing analysis was performed using ANIR POLWIN, a freely available routine developed at CEAM Foundation and distributed by the Eurochamp datacenter: https://data.eurochamp.org/anasoft/ (last access March 2023).

Received: 29 March 2023; Accepted: 24 July 2023

Published online: 15 August 2023

## References

1. Adebiyi, A. et al. A review of coarse mineral dust in the earth system. *Aeolian Res.* **60**, 100849 (2023).

2. Goudie, A. S. & Middleton, N. J. Saharan dust storms: Nature and consequences. *Earth Sci. Rev.* **56**, 179–204 (2001).

3. Knippertz, P. & Stuut, J. B. W. Mineral dust: A key player in the earth system. in (2014).

4. Di Biagio, C., Balkanski, Y., Albani, S., Boucher, O. & Formenti, P. Direct radiative effect by mineral dust aerosols constrained by new microphysical and spectral optical data. *Geophys. Res. Lett.* **47**, e2019GL086186 (2020).

5. di Sarra, A. et al. Shortwave and longwave radiative effects of the intense Saharan dust event of 25–26 March 2010 at Lampedusa (Mediterranean Sea). *J. Geophys. Res. Atmos.* **116**, D23 (2011).

6. Kok, J. F. et al. Smaller desert dust cooling effect estimated from analysis of dust size and abundance. *Nat. Geosci.* **10**, 274–278 (2017).

7. Miller, R. L., Tegen, I. & Perlwitz, J. Surface radiative forcing by soil dust aerosols and the hydrologic cycle. *J. Geophys. Res. Atmos.* **109**, D4 (2004).

8. Slingo, A. et al. Observations of the impact of a major Saharan dust storm on the atmospheric radiation balance. *Geophys. Res. Lett.* https://doi.org/10.1029/2006GL027869 (2006).

9. Cziczo, D. J. et al. Clarifying the dominant sources and mechanisms of cirrus cloud formation. *Science* **340**, 1320–1324 (2013).

10. Froyd, K. D. et al. Dominant role of mineral dust in cirrus cloud formation revealed by global-scale measurements. *Nat. Geosci.* **15**, 177–183 (2022).

11. Hoose, C. & Möhler, O. Heterogeneous ice nucleation on atmospheric aerosols: A review of results from laboratory experiments. *Atmos. Chem. Phys.* **12**, 9817–9854 (2012).

12. Rosenfeld, D., Rudich, Y. & Lahav, R. Desert dust suppressing precipitation: A possible desertification feedback loop. *Proc. Natl. Acad. Sci.* **98**, 5975–5980 (2001).

13. Jickells, Boyd, P. & Hunter, K. A. Biogeochemical impacts of dust on the global carbon cycle. in *Mineral Dust: A Key Player in the Earth System* (eds. Knippertz, P. & Stuut, J.-B. W.) 359–384 (Springer Netherlands, 2014). https://doi.org/10.1007/978-94-017-8978-3_14.

14. Jickells, et al. Global iron connections between desert dust, ocean biogeochemistry, and climate. *Science* **308**, 67–71 (2005).

15. Gaston, C. J. Re-examining dust chemical aging and its impacts on earth's climate. *Acc. Chem. Res.* **53**, 1005–1013 (2020).

16. Usher, C. R., Michel, A. E. & Grassian, V. H. Reactions on mineral dust. *Chem. Rev.* **103**, 4883–4940 (2003).

17. Goudie, A. S. Desert dust and human health disorders. *Environ. Int.* **63**, 101–113 (2014).

18. Querol, X. et al. Monitoring the impact of desert dust outbreaks for air quality for health studies. *Environ. Int.* **130**, 104867 (2019).

19. Formenti, P. et al. Mapping the physico-chemical properties of mineral dust in western Africa: Mineralogical composition. *Atmos. Chem. Phys.* **14**, 10663–10686 (2014).

20. Kandler, K. et al. Size distribution, mass concentration, chemical and mineralogical composition and derived optical parameters of the boundary layer aerosol at Tinfou, Morocco, during SAMUM 2006. *Tellus B Chem. Phys. Meteorol.* **61**, 32–50 (2009).

21. Scheuvens, D., Schütz, L., Kandler, K., Ebert, M. & Weinbruch, S. Bulk composition of northern African dust and its source sediments — A compilation. *Earth-Sci. Rev.* **116**, 170–194 (2013).

22. Atkinson, J. D. et al. The importance of feldspar for ice nucleation by mineral dust in mixed-phase clouds. *Nature* **498**, 355–358 (2013).

23. Di Biagio, C. et al. Variability of the infrared complex refractive index of African mineral dust: Experimental estimation and implications for radiative transfer and satellite remote sensing. *Atmos. Chem. Phys.* **14**, 11093–11116 (2014).

24. Harrison, A. D. et al. Not all feldspars are equal: a survey of ice nucleating properties across the feldspar group of minerals. *Atmos. Chem. Phys.* **16**, 10927–10940 (2016).

25. Journet, E., Desboeufs, K. V., Caquineau, S. & Colin, J.-L. Mineralogy as a critical factor of dust iron solubility. *Geophys. Res. Lett.* https://doi.org/10.1029/2007GL031589 (2008).

26. Li, L. et al. Quantifying the range of the dust direct radiative effect due to source mineralogy uncertainty. *Atmos. Chem. Phys.* **21**, 3973–4005 (2021).

27. Marcotte, A. R., Anbar, A. D., Majestic, B. J. & Herckes, P. Mineral dust and iron solubility: Effects of composition, particle size, and surface area. *Atmosphere* **11**, 533 (2020).

28. Sokolik, I. N. & Toon, O. B. Incorporation of mineralogical composition into models of the radiative properties of mineral aerosol from UV to IR wavelengths. *J. Geophys. Res. Atmos.* **104**, 9423–9444 (1999).

29. Formenti, P. et al. Recent progress in understanding physical and chemical properties of African and Asian mineral dust. *Atmos. Chem. Phys.* **11**, 8231–8256 (2011).

30. Jeong, G. Y. Bulk and single-particle mineralogy of Asian dust and a comparison with its source soils. *J. Geophys. Res. Atmos.* **113**, D2 (2008).

31. Jeong, G. Y. Mineralogy and geochemistry of Asian dust: Dependence on migration path, fractionation, and reactions with polluted air. *Atmos. Chem. Phys.* **20**, 7411–7428 (2020).

32. Kandler, K. et al. Ground-based off-line aerosol measurements at Praia, Cape Verde, during the Saharan mineral dust experiment: Microphysical properties and mineralogy. *Tellus B Chem. Phys. Meteorol.* **63**, 459–474 (2011).

33. Kandler, K. et al. Differences and similarities of central Asian, African, and Arctic dust composition from a single particle perspective. *Atmosphere* **11**, 269 (2020).

34. Ryder, C. L. et al. Coarse-mode mineral dust size distributions, composition and optical properties from AER-D aircraft measurements over the tropical eastern Atlantic. *Atmos. Chem. Phys.* **18**, 17225–17257 (2018).

35. Bristow, C. S., Hudson-Edwards, K. A. & Chappell, A. Fertilizing the amazon and equatorial Atlantic with West African dust. *Geophys. Res. Lett.* https://doi.org/10.1029/2010GL043486 (2010).

36. Caquineau, S., Gaudichet, A., Gomes, L. & Legrand, M. Mineralogy of Saharan dust transported over northwestern tropical Atlantic Ocean in relation to source regions. *J. Geophys. Res. Atmos.* **107**, 4–1 (2002).

37. Formenti, P. et al. Regional variability of the composition of mineral dust from western Africa: Results from the AMMA SOP0/DABEX and DODO field campaigns. *J. Geophys. Res. Atmos.* **113**, D23 (2008).

38. Pye, K. *Aeolian dust and dust deposits* (Academic Press, 1987).

39. Maring, H., Savoie, D. L., Izaguirre, M. A., Custals, L. & Reid, J. S. Mineral dust aerosol size distribution change during atmospheric transport. *J. Geophys. Res. Atmos.* https://doi.org/10.1029/2002JD002536 (2003).
40. Reid, E. A. et al. Characterization of African dust transported to Puerto Rico by individual particle and size segregated bulk analysis. *J. Geophys. Res. Atmos.* **108**, D19 (2003).

41. Scanza, R. A. et al. Modeling dust as component minerals in the community atmosphere model: Development of framework and impact on radiative forcing. *Atmos. Chem. Phys.* **15**, 537–561 (2015).

42. Cuesta, J. et al. Three-dimensional pathways of dust over the Sahara during summer 2011 as revealed by new Infrared Atmospheric Sounding Interferometer observations. *Q. J. R. Meteorol. Soc.* **146**, 2731–2755 (2020).

43. Klüser, L., Kleiber, P., Holzer-Popp, T. & Grassian, V. H. Desert dust observation from space – application of measured mineral component infrared extinction spectra. *Atmos. Environ.* **54**, 419–427 (2012).

44. DeSouza-Machado, S. G., Strow, L. L., Hannon, S. E. & Motteler, H. E. Infrared dust spectral signatures from AIRS. *Geophys. Res. Lett.* https://doi.org/10.1029/2005GL024364 (2006).

45. Capelle, V. et al. Evaluation of IASI-derived dust aerosol characteristics over the tropical belt. *Atmos. Chem. Phys.* **14**, 9343–9362 (2014).

46. Sokolik, I. N. The spectral radiative signature of wind-blown mineral dust: Implications for remote sensing in the thermal IR region. *Geophys. Res. Lett.* **29**, 7–1 (2002).

47. Chudnovsky, A., Ben-Dor, E., Kostinski, A. B. & Koren, I. Mineral content analysis of atmospheric dust using hyperspectral information from space. *Geophys. Res. Lett.* https://doi.org/10.1029/2009GL037922 (2009).

48. Sadrian, M. R., Calvin, W. M., Engelbrecht, J. P. & Moosmüller, H. Spectral characterization of parent soils from globally important dust aerosol entrainment regions. *J. Geophys. Res. Atmos.* **128**, e2022JD037666 (2023).

49. Sadrian, M. R., Calvin, W. M. & McCormack, J. Contrasting mineral dust abundances from X-ray diffraction and reflectance spectroscopy. *Atmos. Meas. Tech.* **15**, 3053–3074 (2022).

50. Go, S. et al. Inferring iron-oxide species content in atmospheric mineral dust from DSCOVR EPIC observations. *Atmos. Chem. Phys.* **22**, 1395–1423 (2022).

51. Green, R. O. et al. The earth surface mineral dust source investigation: an earth science imaging spectroscopy mission. in: *2020 IEEE Aerospace Conference* 1–15 (2020). https://doi.org/10.1109/AERO47225.2020.9172731.

52. Sanwlani, N. & Das, R. Understanding haze: Modeling size-resolved mineral aerosol from satellite remote sensing. *Remote Sens.* **14**, 761 (2022).

53. Clerbaux, C. et al. Monitoring of atmospheric composition using the thermal infrared IASI/MetOp sounder. *Atmos. Chem. Phys.* **9**, 6041–6054 (2009).

54. Crevoisier, C. et al. Towards IASI-New Generation (IASI-NG): impact of improved spectral resolution and radiometric noise on the retrieval of thermodynamic, chemistry and climate variables. *Atmos. Meas. Tech.* **7**, 4367–4385 (2014).

55. Ramsey, M. S. & Christensen, P. R. Mineral abundance determination: Quantitative deconvolution of thermal emission spectra. *J. Geophys. Res. Solid Earth* **103**, 577–596 (1998).

56. Thomson, J. L. & Salisbury, J. W. The mid-infrared reflectance of mineral mixtures (7–14 µm). *Remote Sens. Environ.* **45**, 1–13 (1993).

57. Thorpe, M. T., Rogers, A. D., Bristow, T. F. & Pan, C. Quantitative compositional analysis of sedimentary materials using thermal emission spectroscopy: I. Application to sedimentary rocks. *J. Geophys. Res. Planets* **120**, 1956–1983 (2015).

58. Bell, J. F., Farrand, W. H., Johnson, J. R. & Morris, R. V. Low abundance materials at the mars pathfinder landing site: an investigation using spectral mixture analysis and related techniques. *Icarus* **158**, 56–71 (2002).

59. Metternicht, G. I. & Fermont, A. Estimating erosion surface features by linear mixture modeling. *Remote Sens. Environ.* **64**, 254–265 (1998).

60. Alalam, P., Deschutter, L., Al Choueiry, A., Petitprez, D. & Herbin, H. Aerosol mineralogical study using laboratory and IASI measurements: Application to East Asian deserts. *Remote Sens.* **14**, 3422 (2022).

61. Di Biagio, C. et al. Global scale variability of the mineral dust long-wave refractive index: A new dataset of in situ measurements for climate modeling and remote sensing. *Atmos. Chem. Phys.* **17**, 1901–1929 (2017).

62. Wang, J. et al. Design of a new multi-phase experimental simulation chamber for atmospheric Photosmog, aerosol and cloud chemistry research. *Atmos. Meas. Tech.* **4**, 2465–2494 (2011).

63. Di Biagio, C. et al. Complex refractive indices and single-scattering albedo of global dust aerosols in the shortwave spectrum and relationship to size and iron content. *Atmos. Chem. Phys.* **19**, 15503–15531 (2019).

64. Caponi, L. et al. Spectral- and size-resolved mass absorption efficiency of mineral dust aerosols in the shortwave spectrum: A simulation chamber study. *Atmos. Chem. Phys.* **17**, 7175–7191 (2017).

65. Denjean, C. et al. Size distribution and optical properties of mineral dust aerosols transported in the western Mediterranean. *Atmos. Chem. Phys.* **16**, 1081–1104 (2016).

66. Weinzierl, B. et al. Microphysical and optical properties of dust and tropical biomass burning aerosol layers in the Cape Verde region—an overview of the airborne in situ and lidar measurements during SAMUM-2. *Tellus B Chem. Phys. Meteorol.* **63**, 589–618 (2011).

67. Weinzierl, B. et al. The saharan aerosol long-range transport and aerosol–cloud-interaction experiment: Overview and selected highlights. *Bull. Am. Meteorol. Soc.* **98**, 1427–1451 (2017).

68. Klaver, A. et al. Physico-chemical and optical properties of Sahelian and Saharan mineral dust: in situ measurements during the GERBILS campaign. *Q. J. R. Meteorol. Soc.* **137**, 1193–1210 (2011).

69. Journet, E., Balkanski, Y. & Harrison, S. P. A new data set of soil mineralogy for dust-cycle modeling. *Atmos. Chem. Phys.* **14**, 3801–3816 (2014).

70. Ryder, C. L., Highwood, E. J., Lai, T. M., Sodemann, H. & Marsham, J. H. Impact of atmospheric transport on the evolution of microphysical and optical properties of Saharan dust. *Geophys. Res. Lett.* **40**, 2433–2438 (2013).

71. Panta, A. et al. Insights into the single particle composition, size, mixing state and aspect ratio of freshly emitted mineral dust from field measurements in the Moroccan Sahara using electron microscopy. *Atmos. Chem. Phys. Discuss.* **2022**, 1–40 (2022).

72. Lide, D. R. : C. Handbook of chemistry and physics 1991–1992. (1992).

73. Dorschner, J., Friedemann, C. & Gürtler, J. Laboratory spectra of phyllosilicates and the interstellar 10-micrometer absorption band. *Astron. Nachrichten* **299**, 269–282 (1978).

74. Hudson, P. K., Gibson, E. R., Young, M. A., Kleiber, P. D. & Grassian, V. H. Coupled infrared extinction and size distribution measurements for several clay components of mineral dust aerosol. *J. Geophys. Res. Atmos.* **113**, 789–963 (2008).

75. Hudson, P. K., Young, M. A., Kleiber, P. D. & Grassian, V. H. Coupled infrared extinction spectra and size distribution measurements for several non-clay components of mineral dust aerosol (quartz, calcite, and dolomite). *Atmos. Environ.* **42**, 5991–5999 (2008).

76. Laskina, O., Young, M. A., Kleiber, P. D. & Grassian, V. H. Infrared extinction spectra of mineral dust aerosol: Single components and complex mixtures. *J. Geophys. Res. Atmos.* https://doi.org/10.1029/2012JD017756 (2012).

77. Mogili, P. K., Yang, K. H., Young, M. A., Kleiber, P. D. & Grassian, V. H. Environmental aerosol chamber studies of extinction spectra of mineral dust aerosol components: Broadband IR-UV extinction spectra. *J. Geophys. Res. Atmos.* https://doi.org/10.1029/2007JD008890 (2007).

78. Palchetti, L. et al. FORUM: Unique far-infrared satellite observations to better understand how earth radiates energy to space. *Bull. Am. Meteorol. Soc.* **101**, E2030–E2046 (2020).
79. Glotch, T. D., Rossman, G. R. & Aharonson, O. Mid-infrared (5-100 µm) reflectance spectra and optical constants of ten phyllosilicate minerals. *Icarus* **192**, 605-622 (2007).

80. Glotch, T. D. & Rossman, G. R. Mid-infrared reflectance spectra and optical constants of six iron oxide/oxyhydroxide phases. *Icarus* **204**, 663-671 (2009).

81. Long, L. L., Querry, M. R., Bell, R. J. & Alexander, R. W. Optical properties of calcite and gypsum in crystalline and powdered form in the infrared and far-infrared. *Infrared Phys.* **34**, 191-201 (1993).

82. Baldo, C. et al. Distinct chemical and mineralogical composition of Icelandic dust compared to northern African and Asian dust. *Atmos. Chem. Phys.* **20**, 13521-13539 (2020).

83. Deguine, A. et al. Retrieval of refractive indices of ten volcanic ash samples in the infrared, visible and ultraviolet spectral region. *J. Aerosol. Sci.* **167**, 106100 (2023).

84. Menut, L. et al. Modelling the mineralogical composition and solubility of mineral dust in the Mediterranean area with CHIMERE 2017r4. *Geosci. Model. Dev.* **13**, 2051-2071 (2020).

85. Di Biagio, C., Formenti, P., Styler, S. A., Pangui, E. & Doussin, J.-F. Laboratory chamber measurements of the longwave extinction spectra and complex refractive indices of African and Asian mineral dusts. *Geophys. Res. Lett.* **41**, 6289-6297 (2014).

86. Fratini, G., Ciccioli, P., Febo, A., Forgione, A. & Valentini, R. Size-segregated fluxes of mineral dust from a desert area of northern China by eddy covariance. *Atmos. Chem. Phys.* **7**, 2839-2854 (2007).

87. White, J. U. Long optical paths of large aperture. *J. Opt. Soc. Am.* **32**, 285-288 (1942).

88. Toon, O. B., Pollack, J. B. & Khare, B. N. The optical constants of several atmospheric aerosol species: Ammonium sulfate, aluminum oxide, and sodium chloride. *J. Geophys. Res.* **1896-1977**(81), 5733-5748 (1976).

## Acknowledgements

This work has received funding from the European Union's Horizon 2020 research and innovation program through the EUROCHAMP-2020 Infrastructure Activity under grant agreement no. 730997. It was supported by the French national program LEFE/INSU (Les Enveloppes Fluides et l'Environnement/Institut National des Sciences de l'Univers), and by the OSU-EFLUVE (Observatoire des Sciences de l'Univers-Enveloppes Fluides de la Ville à l'Exobiologie) through dedicated research funding to the RED-DUST project. CNRS-INSU is gratefully acknowledged for supporting the CESAM chamber as a national facility as part of the French ACTRIS Research Infrastructure as well as the AERIS data center (www.aeris-data.fr) for distributing and curing the data produced by the CESAM chamber through the hosting of the EUROCHAMP datacenter (https://data.eurochamp.org). This work has been supported by the Centre National d'Études Spatiales (CNES) via TOSCA/FORUM grant. C. D. B. was supported by the CNES and by the CNRS via the Labex L-IPSL, which is funded by the ANR (grant no. ANR-10-LABX-0018). Fundación CEAM is partly supported by Generalitat Valenciana. The authors wish to acknowledge Prof. Paul Kleiber for help in the reference mineral database constitution and its use.

## Author contributions

J.F.D., C.D.B. and P.F. conceived the study, performed the data analysis, and discussed the results. C.D.B., M.C., and E.P. performed the chamber experiments with contributions by J.F.D. and P.F. M.R. developed the deconvolution software and validated the use of it. J.C. and P.S. contributed to the discussion of the results. C.D.B. led the writing of the manuscript. J.F.D., P.F., M.R., J.C., and P.S. contributed to editing and review of the manuscript.

## Competing interests

The authors declare no competing interests.

## Additional information

**Supplementary Information** The online version contains supplementary material available at https://doi.org/10.1038/s41598-023-39336-7.

**Correspondence** and requests for materials should be addressed to C.D.B.

**Reprints and permissions information** is available at www.nature.com/reprints.

**Publisher's note** Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

<figure><img src="image_7.png" /></figure>

**Open Access** This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.

© The Author(s) 2023