NKFIH KH130337

NKFIH KH130337 Granulometric analysis of recent Saharan dust - Summaries

Summary of the research and its aims for experts
In the course of the proposed research, complex grain size and grain shape analysis of samples extracted from the annually emitted hundreds of millions of tonnes of Saharan mineral dust will be achieved. The collection of samples is based on our previous studies on possible source areas (e.g. Western Sahara; Atlas foreland: Chott Melrhir and Chott Jerid), relatively close proximity of sources (Canary Islands; Antikhytera; Crete), and remote regions (Alps; Carpathian Basin).
By using the automated static image analysis granulometry, we are able to obtain direct grain size and grain shape information on the investigated minerals completed with light transmittance and Raman-spectra (correlated with mineralogy) of the particles. The exact determination of size and shape parameters of irregularly shaped, anisotropic mineral particles is already a key issue in many scientific fields (e.g. sedimentology, geomorphology, atmospheric physics and chemistry, aerobiology), however the accurate answers to these questions are still missing. As part of our proposed project both the size and shape of particles, both their dependence on origin, transport and sedimentation will be explored. So, we could also monitor the changes of environmental conditions by analyses of these parameters.

What is the major research question?
In our research, we will determine by direct methods the particle size and particle shape, main mineralogical character and optical properties of dust material deflated from Saharan source areas. The synoptic meteorological background, transport mechanism and mode of deposition, as well as correlations of these processes with granulometric characteristics would also be defined.
We will also investigate the further possible application of the size and shape parameters gained by automated static image analysis in a deeper understanding of dust transport and deposition processes. Since both the geological and geomorphological characteristics of source areas, the mode and distance of transportation both the deposition mechanisms have an effect on the sampled mineral particles, we will apply multivariate mathematical-statistical methods to determine the weighting scores of different factors in formation of individual granulometric fingerprint. As a result of environmental changes, each of these factors may change over time, so it is more important to know the effects of all these factors.

What is the significance of the research?
Our preliminary results have shown that theoretically obtained, remote sensed and indirectly measured grain size data of Saharan are contradictory and incomparable with the direct results of our automated static image analyses. Consequently, the usage of these previous data causes problem in (1) the determination of quantity of transported dust material and (2) the mass of deposited mineral particles; (3) the deeper understanding of mechanisms of dust transportation; and (4) the determination of effects of dust on climatic, cloud physical and other environmental processes.
Precise and accurate granulometric properties serve as important input data for dust transport and deposition models to help partially eliminate the identified underestimation of deposition of numerical simulations.
Changes in atmospheric circulations through the effects of current climate change can be well observed in the Sahara dust transport processes, especially in the case of North African dust transport to Europe and, consequently, to Hungary. Arctic amplification due to increased warming in Arctic regions influences high-altitude jet-streams, which have higher north-south amplitudes and slower east-west propagation leading to more frequent development of cold drops and longer stationary meridional flow situations. As a consequence, intense and prolonged dust storm situations develop, the seasonal distribution of which also changes; it has been observed also in recent years that the number of episodes in winter and at the beginning of spring is increasing.
Similar examination of aeolian dust particles has not been done so far, and we have only partly dealt with such determination of granulometry of recent dust material. That is why we will be the first to provide detailed direct grain size and grain shape information on Saharan dust particles.

Summary and aims of the research for the public
Billions of tons of mineral dust that is emitted every year through dust storms into the atmosphere are important part of our climatic and environmental system: mineral dust affects the amount of solar radiation reaching the surface; has an effect on aerobiological processes, cloud formation, soil development, health situation. Nevertheless, the level of our scientific knowledge of dust cannot be considered satisfactory. We do not have sufficiently accurate data on (1) the amount; (2) grain size; and (3) mineral composition of Saharan dust transported towards Europe. The inadequate appearance of these three factors makes both the reconstruction and interpretation of past events and estimation of direction and potential effects of future changes difficult and inaccurate.
In the course of our research, we would like to analyse the grain size and previously not investigated grain shape of Saharan dust from source to distant settling regions. All these will be completed with mineralogical analyses, time series derived from satellite images and measurements, synoptic meteorological background and air-mass trajectories, providing fundamental new data about Saharan dust events affecting our country as well as about key properties of mineral dust material. All these data and information would lead to a deeper understanding of the potential impacts of current climate change and to more accurate knowledge if past environmental processes.

Identification of Saharan dust particles in Pleistocene dune sand- paleosol sequences of Fuerteventura (Canary Islands)

New paper: 

Varga, Gy, Roettig, C.-B. (2018). Identification of Saharan dust particles in Pleistocene dune sand-paleosol sequences of Fuerteventura (Canary Islands). Hungarian Geographical Bulletin 67 (2), pp. 121-141.

DOI: https://doi.org/10.15201/hungeobull.67.2.2

Keywords: Saharan dust, Canary Islands, grain size, grain shape, automated image analysis


Automated static image analysis and newly introduced evaluation techniques were applied in this paper to identify Saharan dust material in the unique sand-paleosol sequence of Fuerteventura (Canary Islands). Measurements of ~50,000 individual mineral particles per samples provided huge amount of granulometric data on the investigated sedimentary units. In contrast to simple grain size and shape parameters of bulk samples, (1) parametric curve-fitting allowed the separation of different sedimentary populations suggesting the presence of more than one key depositional mechanisms. Additional (2) Raman-spectroscopy of manually targeted individual particles revealed a general relationship among grain size, grayscale intensity and mineralogy. This observation was used to introduce the (3) intensity based assessment technique for identification of large number of quartz particles. The (4) cluster and (5) network analyses showed that only joint analysis of size, shape and grayscale intensity properties provided suitable results, there is no specific granulometric parameter to distinguish Saharan dust due to their irregular shape characteristics. The presented methods allowed the separation of Saharan dust-related quartz grains from local sedimentary deposits, but due to the lack of robust granulometric characterization of coarsest fractions and due to the diverse geochemical properties of North African sources, exact volumetric amount of deposited dust material and sedimentation rates could not be determined from these data.


Samples were taken from 24 silty units considered as paleo-surfaces of stable geomorphic periods with reduced sand movements and relatively enhanced Saharan dust influence, additional dune sand and sand sheet samples were also investigated as references for intense sand transportation intervals. Detailed description of the units and stratigraphic analysis of selected sites can be found in the works of Faust, D. et al., (2015) and Roettig, C.-B. et al., (2017; under review). Air-dried and 2 mm sieved samples were measured by Malvern Morphologi G3-ID instrument in the Laboratory for Sediment and Soil Analysis (Geographical Institute, Research Centre for Astronomy and Earth Sciences, Hungarian Academy of Sciences).
The applied automated static image analysis technique is a new, innovative mode of grain size and shape analyses completed with chemical identity assessments of Raman spectrometry. In contrast to widely used laser diffraction measurements, image analysis provides direct observational data of particle size, and due to the automatic measurement technique large number of particles are characterized allowing us a more robust and objective granulometric description of particles compared to manual microscopic approaches.

7 mm3 of mineral particles per samples were dispersed by 4 bar compressed air onto a glass slide with 60 s settling time. The used 20× objective lens provide a 960× magnification, suitable for detailed characterization of particles in the size range between fine silt and fine sand fractions. Two-dimensional imaging was completed with the usage of additional vertical focal planes, two additional layers were applied above and two other ones below the focus, equivalent to a total of 27.5 µm.
The captured high-resolution grayscale images of ~50,000 individual mineral particles were automatically analysed by the device software to get a raw granulometric data-matrix. Each row of the table represents one sedimentary particle (with its own identity number), while the columns are various size and shape parameters, completed with light transmissivity data and Raman correlation scores.
Circle-equivalent (CE) diameter is the key size descriptor, calculated as the diameter of a circle with the same area as the projected two-dimensional image of a given mineral grain. Beside several other simple size properties (e.g. length, width, perimeter, sphere-equivalent volume), various shape parameters are derived from these sizes. Aspect ratio is the ratio of width and length, circularity describes the proportional relationship between circumference of a circle equal to the projected area of the particle and perimeter. Convexity (and solidity) parameters are measures of edge roughness by using the ratio of particle and convex hull perimeter (and area). Circularity and convexity values are also suitable to filter out stacked particles and aggregated particles, in this study particles with lower than 0.65 circularity and convexity values were excluded from further calculations.
Intensity mean and standard deviation parameters are determined from the grayscale images as a results of light transmissivity of particles. These values are dependent on mineralogy, particle thickness, chemical homogeneity and surface roughness (for detailed description of the method, see: Varga, Gy. et al., 2018). Intensity values together with chemical identity analyses of the build-in Raman spectrometer provide useful additional information for separation of granulometrically similar particles.

Identification of Saharan dust material

Based on the fact that the Saharan dust deposited at Fuerteventura is mainly (1) silt-sized and (2) contains a lot of quartz particles (regarded as exotic in the basaltic and carbonate-rich environment of the island), these two deterministic factors were evaluated separately to identify North African dust particles. Three different assessment methods were applied to determine the amount of Saharan dust material of the samples.
An indirect approach was applied to theoretically discriminate the silt-sized sedimentary subpopulations which were mathematically separated. The polymodal grain size distribution curves were partitioned into several unimodal Weibull-distributions by applying parametric curve-fitting technique (Sun, D. et al., 2002, 2004; Varga, Gy. et al., in press). According to the applied parametric curve fitting technique the polymodal particle size curves can be interpreted as sum of several, in this case three overlapping Weibull-functions which represent three sediment populations. According to published data on recent dust events from the area (Criado, C. et al., 2003; Menéndez, I. et al., 2007; von Suchodoletz, H. et al., 2009) and measurements of other far-travelled North African dust material (Varga, Gy. et al., 2016), the subpopulation with smallest particles are regarded as the product of long-ranged dust transport.
Raman-spectroscopy (at 785 nm wavelength with 3µm spot) was also applied to directly identify the quartz grains as an indicator of Saharan dust contribution. The acquired spectra of targeted particles were compared to Raman spectral reference libraries using KnowitAll® software from Bio-Rad to identify the minerals present.

The third applied technique was based on the grayscale intensity mean values of particles, the relatively high values were used as a proxy for quartz grains as it was found that there is a strong correlation between light transmissivity and chemical identity (especially in this special case of an environment characterized with the overwhelming majority of carbonate and quartz particles).

Cluster and networks analysis techniques were also applied to differentiate various mineral particle populations based on their general normalized shape (aspect ratio, circularity, convexity, solidity) and grayscale intensity (mean, standard deviation) values. Hierarchical cluster trees were created by using the Euclidean distance pairs of the selected parameters of separated quartz and carbonate size fractions (fine, medium, coarse silt and sand).
For network analysis 192×192 [(24 samples × 2 minerals × 4 size fractions) × (24×2×4)] matrix was compiled, where coefficient of determination was calculated for each pair of records based on the normalized shape and grayscale intensity parameters. This matrix was transformed into an adjacency matrix with values of 0 if r2<0.99 and 1 if r20.99, in this way all of the similar mineral grains were coupled and the whole database can be handled as a network or a finite graph, where the similar records (nodes) are connected (edges) to each other. The Gephi network visualization software was used to analyse the compiled network by applying the ForceAtlas2 continuous graph layout algorithm (Jacomy, M. et al. 2014).

Granulometric characterization of paleosols in loess series by automated static image analysis

An automated image analysis method is proposed here to study the size and shape of siliciclastic sedimentary particles of paleosols of Central European loess sequences. 

New paper: Varga et al. (2018) Sedimentary Geology 370, pp. 1-14. DOI: j.sedgeo.2018.04.001

Several direct and indirect measurement techniques are available for grain size measurements of sedimentary mineral particles. Indirect techniques involve the use of some kind of physical laws, however, all requirements for calculations are in many cases not known. Even so, the direct manual microscopic observation and measurement of large, representative number of grains is time-consuming and sometimes rather subjective. Therefore, automated image analyses techniques provide a new and perspective way to analyse grain size and shape sedimentary particles. 
Here we test these indirect and direct techniques and provide new granulometric data of paleosols. Our results demonstrate that grain size data of the mineral dust samples are strongly dependent on shape parameters of particles, and shape heterogeneity was different of the different size classes. Due to the irregular grain shape parameters, uncertainties have arisen also for the sizes. 
In this paper we present a possible correction procedure to reduce the differences among the results of the laser diffraction and image analysis methods. By applying new correction factors, results of the two approaches could be get closer but the most definite factor, the unknown thickness of particles remained a problem to solve. The other presented method to assess the uncertain 3rd dimension of particles by their intensity-size relationships makes us able to reduce further the deviations of the two sizing methods.

Granulometry - 101

Interpretation of sedimentary (sub)populations extracted from grain size distributions of Central European loess-paleosol series

Varga, Gy., Újvári, G., Kovács, J. (in press). Interpretation of sedimentary (sub)populations extracted from grain size distributions of Central European loess-paleosol series. Quaternary International

Grain size proxies of aeolian dust deposits have widely been applied in environmental and sedimentary studies. However, large body of research papers are not taking into consideration that a complex grain size distribution curve cannot be an indicator of a single one environmental factor (e.g. wind speed/strength, transportation distance, aridity).
The aim of the present paper is to discuss the main differences of frequently used statistical methods and to provide possible interpretations of the results by applying these various approaches on the high-resolution loess-paleosol profile of Dunaszekcső, South Hungary (Central Europe). Beside single statistical descriptors (mean, median, mode) of grain size and simple indices of size-fraction ratios (U-ratio, Grain Size Index), some more complex algorithms were also used in our paper. The applied parametric curve-fitting, end-member modelling and hierarchical cluster analysis techniques are using the whole spectrum of the measured grain size distributions and provide a more reliable and more representative results even in case of small scale variations.
According to our findings, approaches which provide direct linkage among simple statistical descriptors and single atmospheric or other environmental elements are rather oversimplified as properties aeolian dust deposits are influenced by the integrated effects of several concurrent processes. Differences of more complex decomposition methods arise from the different approach and scope. End-members are determined from the unmixing based on the covariance structure of the whole grain size data-series of the section, while the parametric curve-fitting is based on the one-by-one deconvolution of the grain size distribution curves. End-members of loess-paleosol samples are regarded as representation of the average dust grain size distribution of various temporal sediment clusters of seasonal or other short-term intervals, while (sub)populations by parametric curve-fitting are proposed to illustrate process-related elements of background and dust storm depositional components for each sample. Results of cluster analysis represent similar grouping conditions as end-member modelling with a reduced sedimentary and genetically meaning.

How to fit Weibull functions?

You can separated the different sedimentary sub-populations by using a simple MS Excel software with SOLVER. The difference between the measured and modelled values (sum of two/three Weibull distributions) has to be minimized.