Multi-Scale Analysis of Grain Size in the Component Structures of Sediments Accumulated along the Desert-Loess Transition Zone of the Tengger Desert and Implications for Sources and Aeolian Dust Transportation

: Aeolian sediments accumulated along the desert-loess transition zone of the Tengger Desert include heterogeneous textures and complex component structures in their grain-size distributions (GSD). However, the sources of these aeolian sediments have not been resolved due to the lack of large reference GSD sample datasets from adjacent regions that contain various types of sediments; such datasets could be used for fingerprinting based on grain-size properties. This lack of knowledge hinders our understanding of the mechanism of aeolian dust releases in these regions and the effects of forcing of atmospheric circulations on the transportation and accumulation of sediments in this region. In this study, we employed a multi-scale grain-size analysis method, i.e., a combination of the single-sample unmixing (SSU) and the parametric end-member modelling (PEMM) techniques, to resolve the component structures of sediments that had accumulated along the desert-loess transition zone of the Tengger Desert. We have also analyzed the component structures of GSDs of various types of sediments, including mobile and fixed sand dunes, lake sediments, and loess sediments from surrounding regions. Our results demonstrate that the patterns observed in coarser fractions of sediments (i.e., sediments with a mode grain size of >100 µ m) from the transition zone match well with the patterns of component structures of several types of sediments from the interior of the Tengger Desert, and the patterns seen in the finer fractions (i.e., fine, medium, and coarse silts with a modal size of <63 µ m) were broadly consistent with those of loess sediments from the Qilian Mountains. The deflation/erosion of loess from the Qilian Mountains by wind was the most important mechanism underlying the production of these finer grain-size fractions. The East Asia winter monsoon (EAWM) played a key role in transportation of the aeolian dust from these source regions to the desert-loess transition zone of the desert.

Single-sample unmixing (SSU) and parametric end-member modelling (PEMM) of GSDs are two common numerical decomposition techniques used to reveal potential grain-size subpopulations/end-members, and they reveal the component structures of GSDs at different scales [16].Subpopulations defined by the SSU method are derived from individual GSDs and therefore provide process-related information on individual GSDs and the associated between-sample variability, while end-members defined by the PEMM method are derived from the average component structure shared by all GSDs and therefore represent general states of the grain-size components and the associated between-member variability [16].However, the combination of these two different methods to investigate the component structure of GSDs has seldom been applied to the grain-size analysis of aeolian sediments.
Studies on aeolian sediments along the desert-loess transition zone of the Tengger Desert have revealed great variability in grain-size characteristics [17][18][19][20][21] and significant complexity in the component structures of GSDs [22,23] due to significant spatial-temporal variations in both wind direction and wind strength across the region and to complex local topographical and hydrological processes.Guan et al. [20] analyzed 56 surface aeolian samples from the desert-loess transition zone of the Tengger Desert and suggested that the coarse-grain-size fraction (61.8-185.4µm) had a greater degree of sensitivity to the expansion/retraction of the aeolian system compared to finer fractions.They suggested that the coarse fraction originated from the interior of the Tengger Desert and that finer fractions might have more complicated sources.Peng et al. [22] analyzed a large number of GSDs of aeolian sediments along the transition zone using the SSU method and suggested that these samples consist of at least five distinct components, each representing a different potential source, including both local and distant materials from the upwind Tengger Desert, Qilian Mountains, Hexi Corridor, and Gobi desert system in the Alxa Plateau.However, due to the lack of data that could serve as a basis for fingerprinting, specifically data on the various types of sediments from adjacent regions, their grain-size properties, the conclusions of Peng et al. [22] remain speculative and therefore require validation.Liu et al. [23] analyzed the component structures of various types of surface sediments across the Tengger Desert using the PEMM method.They demonstrated that the abundance of different end-members varied between different sediment types and used end-members with coarser grain sizes to indicate aeolian activity.However, the vast majority of samples analyzed by Liu et al. [23] were collected from the interior of the desert.In conclusion, at this time, no study has been conducted to investigate the similarities/differences between the component structures of aeolian sediments within the desert and those along the transition zone and to reveal the potential linkage between them.The resultant knowledge gap hinders our understanding of the sources and transportation processes of aeolian sediments accumulated along the transition zone of the desert.
In this study, we applied the SSU and PEMM methods to GSDs of aeolian sediments collected from wide regions in both the interior and the transition zone of the Tengger Desert to reveal the potential linkage between their component structures.The degrees of variability in grain-size characteristics, structural compositions, and end-members between various types of sediments were investigated and compared in a systematic manner, and the implications for dust release from adjacent regions and for the sources and transportation of aeolian dust along the transition zone are discussed.

Study Area, Samples, and Grain-Size Measurements
The Tengger Desert, located in the drylands of northern China, is one of the major proximal deserts upwind of the Chinese Loess Plateau (CLP).It provides vast amounts of aeolian dust, which influence the regional environments [17,24,25].Variability in aridity and winds, together with complex local topographical and hydrological processes, has led to significant temporal and spatial heterogeneity in the accumulation of aeolian sediments around the desert [24].The desert, located in the margin of the East Asia summer monsoon (EASM), is a sedimentary basin facilitating sediment accumulation by carrying clastic sediments, clay-sized particles, and fine silts from the mountains around the desert to piedmonts and foreland basins [24].The climate of the area is dominated by the EASM in the summer and by the East Asia winter monsoon (EAWM) in the winter.In this study, a total of 486 samples of various types were collected from sites representing wide regions in and surrounding the desert (Figure 1), including (1) sediments from mobile, semi-mobile, and fixed/semi-fixed dunes; interdune depressions; nabkhas; and interdune lakes from the interior of the desert [23]; (2) loess from the northern piedmont of the eastern Qilian Mountains in the southwest of the desert; and (3) aeolian sand and sandy loess from the desert-loess transition zone at the southern margin of the desert [22,24,25].
Atmosphere 2024, 15, x FOR PEER REVIEW 3 of 15 and winds, together with complex local topographical and hydrological processes, has led to significant temporal and spatial heterogeneity in the accumulation of aeolian sediments around the desert [24].The desert, located in the margin of the East Asia summer monsoon (EASM), is a sedimentary basin facilitating sediment accumulation by carrying clastic sediments, clay-sized particles, and fine silts from the mountains around the desert to piedmonts and foreland basins [24].The climate of the area is dominated by the EASM in the summer and by the East Asia winter monsoon (EAWM) in the winter.In this study, a total of 486 samples of various types were collected from sites representing wide regions in and surrounding the desert (Figure 1), including (1) sediments from mobile, semi-mobile, and fixed/semi-fixed dunes; interdune depressions; nabkhas; and interdune lakes from the interior of the desert [23]; (2) loess from the northern piedmont of the eastern Qilian Mountains in the southwest of the desert; and (3) aeolian sand and sandy loess from the desertloess transition zone at the southern margin of the desert [22,24,25].Grain-size measurements were conducted at Hunan University of Science and Technology.Bulk samples (~0.65 g) were pretreated with 10% HCl to remove carbonate and 30% H2O2 to remove organic matter.Then, the samples were further treated with 10% Grain-size measurements were conducted at Hunan University of Science and Technology.Bulk samples (~0.65 g) were pretreated with 10% HCl to remove carbonate and 30% H 2 O 2 to remove organic matter.Then, the samples were further treated with 10% (NaPO 3 ) 6 in an ultrasonic bath and subsequently measured using a Malvern Mastersizer 2000 particle size analyzer.Repeated measurements suggest that the measured GSDs are characterized by a relative standard deviation of less than 3%.

Single-Sample Unmixing
The SSU method optimizes the following residual sum of squares (RSS) to obtain the optimal set of parameters describing a GSD: where x i is the i-th grain-size level, y i and ŷi represent respectively the i-th measured and fitted grain-size volume percentages (a GSD is defined as the variation in y i as a function of x i for a series of i indices ranging from 1 to 101); n is the total number of measured grain-size levels; θ j is the parameters of the j-th subpopulation; g(x i , θ j ) is a probability density function describing a grain-size subpopulation with a set of parameters θ j at a grain-size level x i ; and k is the total number of subpopulations used for unmixing a GSD.
In this study, we used the transformed skew normal function to describe a grain-size subpopulation as follows [16]: The parameters of a subpopulation are θ = {x m , v m , α, ω}.The transformed function contains two free parameters specifying the peak position (x m ) and magnitude (v m ) of a grain-size subpopulation and can be more easily initialized and constrained during the application of the SSU method [16].The fitting quality of the model is measured with the figure-of-merit statistic [27]: The degree of overlap between two adjacent grain-size subpopulations is calculated as follows [22]: where X 1 and X 2 are the grain-size levels corresponding to the maximum volume percentages of two adjacent subpopulations; δ 1 is the half-width on the right side of the first subpopulation; and τ 2 is the half-width on the left side of the second subpopulation.The un-normalized proportion p (which was later normalized to sum up to 100% to obtain the abundance), the mean grain size µ and associated standard deviation σ for a subpopulation are calculated as follows: )) The unmixing results of a GSD contain several subpopulations, each of which has its own parameters: p, µ, and σ.As a result, the unmixed results for a large number of GSDs will contain large sets of parameters for individual subpopulations from a large number of samples.We used the variance-included finite mixture model [28] to classify the pooled unmixed mean grain sizes (µ) by taking the standard deviation (σ) of subpopulations into account [22].The optimal number of clusters was determined using the Bayesian information criterion of Schwartz [29].The unmixing and statistical analysis were implemented using the efficient and flexible R program "ussGSD" [22] which is freely available at the GitHub repository https://github.com/pengjunUCAS/ussGSD,accessed on 1 December 2023.

Parametric End-Member Analysis
The PEMM method tries to minimize the difference between the measured GSDs (X) and the modelled GSDs ( X). X and X are m × n matrices (m is the number of GSDs and n the number of grain-size levels).X can be described by the matrix product of abundances (M) and end-members (B) [30], as follows: where M is a m × k matrix (k is the number of end-members) of numbers that are nonnegative and constrained such that the sum of the entries in each row is equal to 100%.B is a k × n matrix that can be described as follows: where φ j is the parameters of the j-th subpopulation and f (x i , φ j ) is a probability density function describing a grain-size subpopulation with a set of parameters φ j at grain-size level x i .The transformed skew normal function describing an end-member is as follows [16]: The parameters of a subpopulation are φ = {x m , α, ω}.D is an intermediate quantity described by Equation (7).In this study, the PEMM method was implemented using the R function "PEMM" [16].The difference between the measured GSDs and the fitted values was measured using the mean of the angles (and the mean of squared Pearson correlation coefficients) of the row vectors of matrices X and X.

Grain-Size Characteristics
We summarized the volume percentages of clay, silt, and sand fractions for different types of sediments using a ternary plot (Figure 2).The plot shows that for samples of mobile and semi-mobile dunes, and for most samples of interdune depressions, sand content was greater than 96%, silt content was lower than 3.4%, and clay content was lower than 0.6%.For most samples of fixed/semi-fixed dunes and nabkhas, sand content was greater than 80%, silt content was lower than 15%, and clay content was lower than 5%.Samples from interdune lake surfaces had greater contents of silt and clay and lower content of sand compared to samples of mobile and semi-mobile dunes, interdune depressions, fixed/semifixed dunes and nebkhas.By contrast, most samples of loess from the Qilian Mountains had sand content less than 15%, silt content greater than 75%, and silt content less than 10%.For samples from the transition zone, the content of sand (from 3% to 99%) and silt (from less than 1% to 80%) varied widely, while the content of clay was less than 17%, suggesting that these samples demonstrated the greatest variability in the contents of clay, silt, and sand compared to other types of sediments.)) The parameters of a subpopulation are φ = {xm, α, ω}.D is an intermediate quantity described by Equation (7).In this study, the PEMM method was implemented using the R function "PEMM" [16].The difference between the measured GSDs and the fitted values was measured using the mean of the angles (and the mean of squared Pearson correlation coefficients) of the row vectors of matrices X and  .

Grain-Size Characteristics
We summarized the volume percentages of clay, silt, and sand fractions for different types of sediments using a ternary plot (Figure 2).The plot shows that for samples of mobile and semi-mobile dunes, and for most samples of interdune depressions, sand content was greater than 96%, silt content was lower than 3.4%, and clay content was lower than 0.6%.For most samples of fixed/semi-fixed dunes and nabkhas, sand content was greater than 80%, silt content was lower than 15%, and clay content was lower than 5%.Samples from interdune lake surfaces had greater contents of silt and clay and lower content of sand compared to samples of mobile and semi-mobile dunes, interdune depressions, fixed/semi-fixed dunes and nebkhas.By contrast, most samples of loess from the Qilian Mountains had sand content less than 15%, silt content greater than 75%, and silt content less than 10%.For samples from the transition zone, the content of sand (from 3% to 99%) and silt (from less than 1% to 80%) varied widely, while the content of clay was less than 17%, suggesting that these samples demonstrated the greatest variability in the contents of clay, silt, and sand compared to other types of sediments.The relationships between standard deviation (in unit Φ), skewness, and mean grain size (in unit Φ) among various types of sediments differed (Figure 3).For samples from mobile and semi-mobile dunes and interdune depressions, the variation in sorting (standard deviation) and skewness as a function of the mean size was not significant.The mean sizes were smaller than 3 Φ.The standard deviation and skewness of most samples were centered at 0.5 Φ and 0, respectively, suggesting good sorting and symmetry.For samples The relationships between standard deviation (in unit Φ), skewness, and mean grain size (in unit Φ) among various types of sediments differed (Figure 3).For samples from mobile and semi-mobile dunes and interdune depressions, the variation in sorting (standard deviation) and skewness as a function of the mean size was not significant.The mean sizes were smaller than 3 Φ.The standard deviation and skewness of most samples were centered at 0.5 Φ and 0, respectively, suggesting good sorting and symmetry.For samples from fixed/semi-fixed dunes and nebkhas, the sorting and skewness were positively correlated with the grain size; the sorting increased and the skewness decreased as the mean size became coarser.For samples from the interdune lake surfaces, the sorting was positively correlated with the mean size.However, the variation in skewness as a function of mean size did not follow a uniform pattern.For almost all samples of loess from the Qilian Mountains, the sorting was positively correlated with the mean size, but the skewness was negatively corrected with the mean size.Samples from the transition zone demonstrated the greatest variability in the patterns of correlation between the mean size and the other two statistics.
from fixed/semi-fixed dunes and nebkhas, the sorting and skewness were positiv related with the grain size; the sorting increased and the skewness decreased as th size became coarser.For samples from the interdune lake surfaces, the sorting w tively correlated with the mean size.However, the variation in skewness as a fun mean size did not follow a uniform pattern.For almost all samples of loess from th Mountains, the sorting was positively correlated with the mean size, but the skewn negatively corrected with the mean size.Samples from the transition zone demon the greatest variability in the patterns of correlation between the mean size and th two statistics.

Structural Compositions
The SSU results for representative GSDs of various types of sediments are pr in Figure 4.The FOM values were smaller than 5%.For GSDs unmixed with mo one subpopulation, the minimum degree of overlapping between adjacent subp tions was greater than 0.7, indicating the unmixing results were acceptable.The d tions of the unmixed grain-size means of various types of sediments are shown as plified) radial plot [31] in Figure 5.For samples from mobile dunes, the pooled un results could be clustered only as one group (Figure 5a).Samples from semi-mobil and interdune depressions could be best described by two clusters (Figure 5b,c).S of fixed/semi-fixed dunes, nebkhas, and loess from the Qilian Mountains were de by four clusters (Figure 5d,e,g).Samples from the interdune lake surfaces and the tion zone formed five distinct clusters (Figure 5f,h).

Structural Compositions
The SSU results for representative GSDs of various types of sediments are presented in Figure 4.The FOM values were smaller than 5%.For GSDs unmixed with more than one subpopulation, the minimum degree of overlapping between adjacent subpopulations was greater than 0.7, indicating the unmixing results were acceptable.The distributions of the unmixed grain-size means of various types of sediments are shown as a (simplified) radial plot [31] in Figure 5.For samples from mobile dunes, the pooled unmixing results could be clustered only as one group (Figure 5a).Samples from semi-mobile dunes and interdune depressions could be best described by two clusters (Figure 5b,c).Samples of fixed/semi-fixed dunes, nebkhas, and loess from the Qilian Mountains were described by four clusters (Figure 5d,e,g).Samples from the interdune lake surfaces and the transition zone formed five distinct clusters (Figure 5f,h). .SSU results of grain-size distribution for various types of sediments, using the skew nor mal distribution."SN0" is the transformed skew normal distribution; "mrsl" is the minimum degree of overlapping between adjacent subpopulations; "RSS" and "FOM" are the residual sum of square and the figure-of-merit value, respectively.Note that the x-axis is plotted on a log scale for bette visualization.SSU results of grain-size distribution for various types of sediments, using the skew normal distribution."SN0" is the transformed skew normal distribution; "mrsl" is the minimum degree of overlapping between adjacent subpopulations; "RSS" and "FOM" are the residual sum of squares and the figure-of-merit value, respectively.Note that the x-axis is plotted on a log scale for better visualization.
We compared the clustering results, including the mean grain size and the normalized abundance of individual clusters, between various types of sediments; the results are presented in Figure 6.Almost all types of sediments (except loess from the Qilian Mountains) predominantly contained a coarse subpopulation in the size range 100-250 µm.The abundance of this coarse fraction decreased gradually as the sediment types changed from mobile and semi-mobile dunes to interdune depressions, fixed/semi-fixed dunes, and nebkhas.In addition, samples from fixed/semi-fixed dunes, nebkhas, interdune lake surfaces, and desert-loess transition zone contained an extremely coarse fraction in the interval of 500-1000 µm.Samples of fixed/semi-fixed dunes and nebkhas demonstrated very similar structural compositions.The structural compositions of samples from interdune lake surfaces and transition zone were also relatively consistent.Samples from fixed/semi-fixed dunes, nebkhas, interdune lake surfaces, loess from the Qilian Mountains, and the transition zone contained two subpopulations in the size range 2-32 µm.In addition, samples from interdune lake surfaces, loess from the Qilian Mountains, and transition zone contained a clay fraction.

End-Member Variations
The distributions of end-members of various types of sediments are presented in Figure 7.The optimized mean angles between the measured and modelled GSDs for various types of sediments were smaller than 5 • , indicating good fit quality.Sediments of mobile and semi-mobile dunes were modelled using four closely overlapping endmembers, suggesting that the PEMM algorithm had an intrinsic tendency to generate extra end-members only within grain-size intervals associated with larger volume percentages as the number of end-members increased, as pointed out by Peng et al. [16].Loess from the Qilian Mountains also yielded four end-members, but they were significantly much finer.By contrast, the remaining types of sediments yielded at least five end-members.
Figure 4. SSU results of grain-size distribution for various types of sediments, using the skew nor mal distribution."SN0" is the transformed skew normal distribution; "mrsl" is the minimum degree of overlapping between adjacent subpopulations; "RSS" and "FOM" are the residual sum of square and the figure-of-merit value, respectively.Note that the x-axis is plotted on a log scale for bette visualization.Comparisons of the modal size and the corresponding abundance of end-members of various types of sediments are presented in Figure 8. Sediments from mobile and semimobile dunes and interdune depressions contained large fractions of fine-to-coarse sand and those from fixed/semi-fixed dunes were dominated by a fine-sand fraction.Sediments from nebkhas, interdune lake surfaces, and the desert-loess transition zone contained non-negligible fractions of silt and various sands.The modal numbers of end-members decreased gradually between samples from mobile dune and sediments from interdune lake surfaces.The modal sizes of all four end-members for loess from the Qilian Mountains represented silt fractions, suggesting that their sources were intrinsically different from those of sediments from the interior of the desert.In addition, a comparison of the SSU and PEMM results (Figures 6 and 8) suggested that the clusters of subpopulations derived by the SSU method are distributed more evenly across the grain-size scales, while the PEMM end-members fall within much narrower grain-size intervals, suggesting that SSU has a greater capacity to recognize independent grain-size components as compared to PEMM [22].
faces and transition zone were also relatively consistent.Samples from fixed/semi-fixed dunes, nebkhas, interdune lake surfaces, loess from the Qilian Mountains, and the transition zone contained two subpopulations in the size range 2-32 µm.In addition, samples from interdune lake surfaces, loess from the Qilian Mountains, and transition zone contained a clay fraction.[22] for various types of sediments.Note that the y-axis is plotted on a log scale for better visualization.

End-Member Variations
The distributions of end-members of various types of sediments are presented in Figure 7.The optimized mean angles between the measured and modelled GSDs for various types of sediments were smaller than 5°, indicating good fit quality.Sediments of mobile and semi-mobile dunes were modelled using four closely overlapping end-members, suggesting that the PEMM algorithm had an intrinsic tendency to generate extra end-members only within grain-size intervals associated with larger volume percentages as the number of end-members increased, as pointed out by Peng et al. [16].Loess from the Qilian Mountains also yielded four end-members, but they were significantly much finer.By contrast, the remaining types of sediments yielded at least five end-members.Comparisons of the modal size and the corresponding abundance of end-members of various types of sediments are presented in Figure 8. Sediments from mobile and semimobile dunes and interdune depressions contained large fractions of fine-to-coarse sand and those from fixed/semi-fixed dunes were dominated by a fine-sand fraction.Sediments from nebkhas, interdune lake surfaces, and the desert-loess transition zone contained non-

Implication for Sources of Aeolian Sediments from the Desert-Loess Transition Zone
The proportions of clay, silt, and sand (Figure 2), the mean grain size, the sorting, and the skewness of the GSDs (Figure 3) are different between various types of sediments.Data points representing six types of sediments from the interior of the desert and loess from the northern piedmont of the eastern Qilian Mountains cover only certain parts of the ternary plot in Figure 2 and the scatterplot in Figure 3.By contrast, data points representing aeolian sediments from the desert-loess transition zone cover almost the whole of both plots.Thus, it can be seen that aeolian sediments from the transition zone have the greatest grain-size variability and that they inherit the grain-size characteristics of all other types of sediments.These lines of evidence imply that their GSDs formed as an admixture of various sources of sediments, a finding in accordance with conclusions drawn from previous studies [22,23].
The SSU and PEMM methods that reveal grain-size component structures on different scales were used to analyze various types of sediments (see Figures 6 and 8).For aeolian sediments from the transition zone, the first three clusters of subpopulations defined by the SSU method (with mean sizes of 0.6 µm, 3.5 µm, and 14.2 µm) are closely related to the first two clusters of samples from fixed/semi-fixed dunes and nebkhas and the first three clusters of samples from interdune lake surfaces and loess from the Qilian Mountains (Figure 6).In addition, for aeolian sediments from the transition zone, the first two end-members defined by the PEMM method (with modal sizes of 18.4 µm and 50.1 µm) can be broadly associated with certain end-members of samples from interdune lake surfaces within the Tengger Desert and loess from the Qilian Mountains (Figure 8).These results suggest that the sources of the finer fractions of aeolian sediments in the transition  The proportions of clay, silt, and sand (Figure 2), the mean grain size, the sorting, and the skewness of the GSDs (Figure 3) are different between various types of sediments.Data points representing six types of sediments from the interior of the desert and loess from the northern piedmont of the eastern Qilian Mountains cover only certain parts of the ternary plot in Figure 2 and the scatterplot in Figure 3.By contrast, data points representing aeolian sediments from the desert-loess transition zone cover almost the whole of both plots.Thus, it can be seen that aeolian sediments from the transition zone have the greatest grain-size variability and that they inherit the grain-size characteristics of all other types of sediments.These lines of evidence imply that their GSDs formed as an admixture of various sources of sediments, a finding in accordance with conclusions drawn from previous studies [22,23].
The SSU and PEMM methods that reveal grain-size component structures on different scales were used to analyze various types of sediments (see Figures 6 and 8).For aeolian sediments from the transition zone, the first three clusters of subpopulations defined by the SSU method (with mean sizes of 0.6 µm, 3.5 µm, and 14.2 µm) are closely related to the first two clusters of samples from fixed/semi-fixed dunes and nebkhas and the first three clusters of samples from interdune lake surfaces and loess from the Qilian Mountains (Figure 6).In addition, for aeolian sediments from the transition zone, the first two endmembers defined by the PEMM method (with modal sizes of 18.4 µm and 50.1 µm) can be broadly associated with certain end-members of samples from interdune lake surfaces within the Tengger Desert and loess from the Qilian Mountains (Figure 8).These results suggest that the sources of the finer fractions of aeolian sediments in the transition zone may be the Qilian Mountains and the Tengger Desert.Such finer grain-size fractions have been reported from loess sediments from various regions [7,[32][33][34][35].
There is debate on the production of these finer grain-size fractions (fine, medium, and coarse silts).Crouvi et al. [36] proposed aeolian abrasion as an important mechanism for the production of coarse silt quartz grains.Amit et al. [37] revealed a genetic association between the coarse silt fraction of the loess in CLP and the primary production of coarse silt in the proximal Mu Us, Tengger, and Badain Jaran deserts and suggested that coarse silts in the loess of CLP were produced by abrasion of quartz sand within these deserts.On the other hand, the empirical findings of Swet et al. [38] suggested that less dust was generated by sand abrasion than by the removal of clay minerals coated on sand-grain surfaces and that non-sandy soils were the major source of dust.Guan et al. [20] found that the Tengger Desert is not the main source area for the silt grain component of size 18.2-61.8µm identified in sediments from the transition zone.Liu et al. [23] identified a PEMM end-member characterized by a modal size of 28.3 µm in their GSDs, which were collected from the interior of the Tengger Desert, and suggested that the component originated from distant source areas rather than from local regions.We demonstrated that the finer grain-size fractions in aeolian sediments from the transition zone could be most easily associated with loess in the Qilian Mountains, which contained a large amount of silt (see Figures 6 and 8).Loess sediments were widely distributed along the Qilian Mountains.Their accumulation started as early as 0.83-0.85Ma ago, and they reached a maximum thickness of >200 m [39,40].Many studies [41][42][43] have suggested that sediments along the Qilian Mountains might have been a major source of silts deposited on the western part of the CLP.Considering that the transition zone is located between the desert and the CLP, the Qilian Mountains should also be an important source of the finer grains accumulated there.Accordingly, we suggest that the deflation/erosion of loess from the Qilian Mountains by wind was the most important mechanism underlying the production of these finer grain-size fractions (see Section 4.2).Interdune lake surfaces within the Tengger Desert were the second potential source for the finer grain-size fractions, although they contain less silt compared to the loess in the Qilian Mountains (Figures 6 and 8).
For all types of sediments collected from the interior of the desert, clusters of SSU subpopulations with coarser mean sizes are generally consistent with those of sediments from the transition zone, which have a mean size of 128.8 µm (fine sand) (Figure 6).For sediments from the transition zone, the fifth subpopulation cluster, with a mean size of 636 µm (coarse sand), matches the last subpopulation clusters for fixed/semi-fixed dunes, nebkhas, and sediments from interdune lake surfaces.In addition, PEMM end-members with coarser mean sizes are found in all types of sediments collected from the interior of the desert, but they are absent from loess from the Qilian Mountain (Figure 8).Grainsize components with such coarser mean sizes have been reported in many dune fields in northern China [3,4,9] and they can be transported only over short distances.These points suggest that the only source of the coarser grain-size fractions accumulated in the transition zone is the interior of the Tengger Desert, which is adjacent and located in the upwind direction.A proximal supply of coarse-grained dust upwind has been reported not only from sand/dune fields in northern China, such as the Mu Us Desert [44], Qinghai Lake basin [4], and the Kumtagh Desert [9], but also from aeolian sediments worldwide, for example, in aeolian/loess sediments from the Mediterranean [45], Egypt [46], and Israel [36].

Implications for Aeloian Dust Transportation
The first two SSU clusters, with central means of 0.6 µm and 3.5 µm (not identified by the PEMM method), have been widely reported in aeolian sediments from northern China [19,32].The enrichment of clay and very fine silt can result from pedogenesis and weathering of finer grains, such as primary loess [47].According to Pye [48], aeolian grains with a diameter <20 µm can be carried by the wind up to thousand kilometers.Their presence can be related to long-term inputs originating from both local and distant fluvial, lacustrine, and pedogenic components [24] and carried by local wind regimes and the upper-level westerlies; they can also be the result of the aggregation/adhering of fine particles to larger ones [49].
The EAWM is the main transporting force for loess sediments in the CLP [44].The desert-loess transition zone of the Tengger Desert is upwind of the CLP, and it is also controlled by the East Asian monsoon system [17,24].The SSU cluster with a mean size of 14.2 µm and the PEMM cluster with modal sizes of 18.4 µm and 50.1 µm are associated with typical loess sediments from northern China [32,34,35].Aeolian grains with diameters ranging between 20 µm and 70 µm are suspended above the surface in the short-term [48].Loess-like sediments consisting mainly of silt can be transported from the Qilian Mountains to the transition zone through the EAWM.During the winter and spring, the dominant nearsurface wind direction along the desert margin is from the northwest [19,50], which causes frequent dust storms as a result of the dry conditions that are typical during these months.Sedimentologic evidence of strengthened EAWM and enhanced release of aeolian dust has been identified in several sites around the south margin of the Tengger Desert [17,23,25].
The SSU cluster with a mean size of 128.8 µm and the PEMM clusters with modal sizes of 99.1 µm, 165.3 µm, and 334.7 µm can be transported only by extremely strong near-surface winds in saltation or creep mode.For example, Guan et al. [51] suggested that a grain-size fraction between 209 µm and 550 µm identified from the Shagou loess section in the southwest margin of the Tengger Desert is an indication that severe sandstorms occurred in the desert during the last interglacial period.It has been widely recognized that strong near-surface winds of this kind can be triggered by strong EAWM events [51,52].Lastly, considering that fluvial strata are widely distributed along the southern margin of the desert, the SSU cluster with a mean size of 636 µm may originate from the reworking of near-source intermittent fluvial sediments and/or transportation of very coarse sand from the desert by extremely strong sandstorms [24].

Conclusions
A combination of the SSU and PEMM methods has been used to resolve the component structures of sediments within and surrounding the desert-loess transition zone of the Tengger Desert.The SSU method has greater capacity to identify independent grain-size components compared to the PEMM method.Aeolian sediments accumulated along the transition zone are formed as an admixture of various types of sediments.The coarser and finer grain-size fractions resemble (inherit) characteristics and component structures of sand sediments from the interior of the Tengger Desert and loess sediments from the Qilian Mountains, respectively.The coarser grain-size fractions have proximal sources and are transported downwind from the interior of the desert through strong EAWM events and extremely strong sandstorms, while the finer grain-size fractions are mainly transported from the Qilian Mountains through continuous EAWM during the winter and spring.By contrast, significant contribution of silts from the interior of the Tengger Desert is demonstrated only in sediments from interdune lake surfaces.

Figure 1 .
Figure 1.(a) grain-size sampling sites around the Tengger Desert.Different types of sediments are shown in different colors.(b-j) representative photos of sampling sites for various types of sediments, including (b) mobile dunes, (c) semi-mobile dunes, (d) interdune depressions, (e) fixed/semifixed dunes, (f) nebkhas, (g) interdune lake surfaces, (h) loess from the Qilian Mountains, (i) and (j) sediments from the desert-loess transition zone.

Figure 1 .
Figure 1.(a) grain-size sampling sites around the Tengger Desert.Different types of sediments are shown in different colors.(b-j) representative photos of sampling sites for various types of sediments, including (b) mobile dunes, (c) semi-mobile dunes, (d) interdune depressions, (e) fixed/semi-fixed dunes, (f) nebkhas, (g) interdune lake surfaces, (h) loess from the Qilian Mountains, (i,j) sediments from the desert-loess transition zone.

Figure 2 .
Figure 2. Ternary plot showing abundances of clay, silt, and sand for various types of sediments.

Figure 2 .
Figure 2. Ternary plot showing abundances of clay, silt, and sand for various types of sediments.

Figure 3 .
Figure 3. Variations in (a) sorting and (b) skewness as a function of mean grain size for vario of sediments.

Figure 3 .
Figure 3. Variations in (a) sorting and (b) skewness as a function of mean grain size for various types of sediments.

Figure 4
Figure 4. SSU results of grain-size distribution for various types of sediments, using the skew nor mal distribution."SN0" is the transformed skew normal distribution; "mrsl" is the minimum degree of overlapping between adjacent subpopulations; "RSS" and "FOM" are the residual sum of square and the figure-of-merit value, respectively.Note that the x-axis is plotted on a log scale for bette visualization.

Figure 5 .
Figure 5. Simplified radial plots showing the classification of the means of grain-size component

Figure 4 .
Figure 4. SSU results of grain-size distribution for various types of sediments, using the skew normal distribution."SN0" is the transformed skew normal distribution; "mrsl" is the minimum degree of overlapping between adjacent subpopulations; "RSS" and "FOM" are the residual sum of squares and the figure-of-merit value, respectively.Note that the x-axis is plotted on a log scale for better visualization.

Figure 5 .
Figure 5. Simplified radial plots showing the classification of the means of grain-size component extracted from the pooled unmixing results using the finite mixture model-based Bayesian

Figure 5 .
Figure 5. Simplified radial plots showing the classification of the means of grain-size components extracted from the pooled unmixing results using the finite mixture model-based Bayesian clustering algorithm [22] for various types of sediments, including (a) mobile dunes, (b) semi-mobile dunes, (c) interdune depressions, (d) fixed/semi-fixed dunes, (e) nebkhas, (f) interdune lake surfaces, (g) loess from the Qilian Mountains, (h) sediments from the desert-loess transition zone.The gray region indicates the two-sigma range of a cluster.

Figure 6 .
Figure 6.Comparison of the means and abundances for clusters of grain-size components classified using the finite mixture model-based Bayesian clustering algorithm[22] for various types of sediments.Note that the y-axis is plotted on a log scale for better visualization.

Figure 6 .
Figure 6.Comparison of the means and abundances for clusters of grain-size components classified using the finite mixture model-based Bayesian clustering algorithm[22] for various types of sediments.Note that the y-axis is plotted on a log scale for better visualization.

Figure 7 .
Figure 7. Variations of probability density distributions of end-members for various types of sediments.Note that the x-axis is plotted on a log scale for better visualization.

Figure 7 .
Figure 7. Variations of probability density distributions of end-members for various types of sediments.Note that the x-axis is plotted on a log scale for better visualization.

Figure 8 .
Figure 8.Comparison of the modal sizes and abundances of end-members for various types of sediments.

Figure 8 .
Figure 8.Comparison of the modal sizes and abundances of end-members for various types of sediments.

1 .
Implication for Sources of Aeolian Sediments from the Desert-Loess Transition Zone