Aeolian Dust Preserved in the Guliya Ice Cap (Northwestern Tibet): A Promising Paleo-Environmental Messenger

: Asian aeolian dust is a primary factor in Northern Hemisphere atmospheric dynamics. Predicting past and future changes in atmospheric circulation patterns relies in part on sound knowledge of Central Asian dust properties and the dust cycle. Unfortunately for that region, data are too sparse to constrain the variation in dust composition over time. Here, we evaluate the potential of a Tibetan ice core to provide a comprehensive paleo-atmospheric dust record and thereby reduce uncertainties regarding mineral aerosols’ feedback on the climate system. We present the ﬁrst datasets of the mineralogical, geochemical, and Sr-Nd isotope composition of aeolian dust preserved in pre-Holocene layers of two ice cores from the Guliya ice cap (Kunlun Mountains). The composition of samples from the Summit (GS; 6710 m a.s.l.) and Plateau (GP; 6200 m a.s.l.) cores reveals that the characteristics of the dust in the cores’ deepest ice layers are signiﬁcantly different. The deepest GS layers reveal isotopic values that correspond to aeolian particles from the Taklimakan desert, contain a mix of ﬁne and coarse grains, and include weathering-sensitive material suggestive of a dry climate at the source. The deep GP layers primarily consist of unusual nodules of well size-sorted grey clay enriched in weathering-resistant minerals and elements typically found in geothermal waters, suggesting that the dust preserved in the oldest GP layers originates from a wet and possibly anoxic source. The variability of the dust composition highlighted here attests to its relevance as a paleo-environmental messenger and warrants further exploration of the particularly heterogenous Guliya glacial dust archive.


Introduction
Airborne mineral dust (i.e., microparticles) has the largest mass emission rate of all aerosol types [1] and constitutes 70% of the global aerosol mass burden [2].Dust particles, emitted by aeolian erosion of arid and semi-arid soils, affect the physical and chemical properties of the atmosphere, thereby influencing Earth's climate system [3].Aeolian microparticles change the atmospheric radiation balance through their absorption and scattering properties, modify cloud microphysics by acting as condensation nuclei, influence marine biochemistry and oceanic uptake of atmospheric CO 2 , fertilize terrestrial ecosystems, and impact atmospheric chemistry.All of these processes are directly linked to the particles' size, morphology, geochemical and mineralogical composition, and their atmospheric transport pathways.In turn, the parameters controlling mineral dust emissions (i.e., aridity, vegetation and snow cover, soil types, wind velocity, precipitation, and topographic features) are very sensitive to climate variations [4] and involve feedbacks that are recognized as major uncertainties of climate change [5].
Characterizing mineralogical and geochemical fingerprints of aeolian dust and determining their change through time enables the extraction of paleo-environmental and paleo-synoptic information from aeolian deposits, and thereby reduces uncertainty regarding the role of dust in the climate cycle.However, despite the valuable paleo-environmental message carried by aeolian dust, paleo-dust records remain globally scarce and underexplored [6].Surprisingly, a paleo-dust record is also lacking in NW Tibet, one of the largest atmospheric dust source regions in the Northern Hemisphere.This research provides the first mineralogical and geochemical information on aeolian microparticles deposited on the Guliya ice cap (NW Tibet) during the Late Pleistocene and demonstrates the paleoenvironmental value that paleo-dust records extracted from the Guliya ice cores likely hold.
Ice cores are commonly used as terrestrial archives of past atmospheric microparticle deposition (e.g., [7]).Based on the principle that the particle size at the deposition site becomes finer with distance from the source, the source-to-sink relative distance, paleo-winds strength, and aridity are inferred from the particles' size and concentration, as measured in ice.In polar ice, for instance, the concentration of microparticles is much higher during glacial than interglacial periods [8][9][10].However, this is not always the case in tropical ice cores.For example, lower dust concentration in the Sajama (Bolivia) core for the last glacial period reflects the wetter conditions and a more extensive snow cover [11].Not only do the particles' physical properties reflect past atmospheric dynamics and hydroclimate; the mineralogical, isotopic, and elemental composition of the dust provide environmental information regarding the source of windborne material.Physical and geochemical characteristics thereby complement each other in constraining the geographical source of the dust, its distance to the deposition site and the atmospheric transport conditions (e.g., wind strength).
Trace element (TE) ratios and rare earth element (REE) abundances and distribution patterns are other source-specific parameters of mineral particles that can be conserved during long-range windborne transport.REEs are particularly well suited for investigating mineral aerosol characteristics that are controlled by provenance and pedogenetic processes [12][13][14].Because their variability results from changes in the mineral composition of aeolian dust [15], elemental ratios complement the panel of geochemical tools as diagnostic source tracers.
The combination of strontium ( 87 Sr/ 86 Sr) and neodymium ( 143 Nd/ 144 Nd) isotopic ratios is another powerful source-area fingerprint for dust as these ratios depend on the isotopic composition and age of the rock from which they are derived [16][17][18][19][20][21][22][23].The isotopic compositions of the ice core dust samples compared to those of Potential Source Areas (PSA) samples indicate variations in the mineral aerosols' geographical provenance, which may be used to resolve past climate changes [17,24,25].For example, the first investigations of the isotopic fingerprint of airborne dust from ice cores suggested that Greenland is influenced by proximal sources during warm periods and by remote Asian sources (e.g., the Gobi Desert) in colder periods [26,27], until subsequent studies excluded the Gobi Desert as a potential dust source [28].The most recent and extensive investigations of the main Greenland dust sources during glacial periods support the Taklimakan and Tengger deserts, and to a lesser extent Central Europe as the more likely contributors [29][30][31].On a seasonal scale, the dust source fluctuation between the Taklimakan desert in spring and Inner Mongolia in autumn is linked to the Arctic Oscillation phase and to the position of the Westerly jet.Various studies [22,32,33] provided new evidence of a latitudinal gradient in the Sr and Nd isotopic composition of cryoconite, and snowpack dust collected on five glaciers of the northeast Tibetan Plateau (TP), from which they infer dust atmospheric transport pathways.On an event-based time scale, aeolian dust from the Gobi Desert and the Chinese Loess Plateau (CLP) was also identified by its Sr-Nd isotopic signature in the snow of Mt.St. Elias (Yukon Territory, Canada), following a dust storm event that occurred in Asia and was transported across the North Pacific Ocean by clouds [34].
High elevation glaciers on the TP, dominating one of Earth's most arid regions, are repositories of tropospheric microparticles that have experienced mixing and gravitational sorting during transport.Even though local (i.e., coarser grain size fraction) soil dust is a non-negligible source, the particles entrapped in these high-altitude ice fields are substantially sorted by grain size with longer range transport enhancing sorting, making long range transported dust a prevalent deposit on TP glaciers [33,35].The major dust storm events, during which particles rise above the boundary layer and are then transported for thousands of kilometers, are most important to climate systems.Therefore, the high altitude Guliya ice cap (western Kunlun Mountains, 35.12 to 35.22 • N, 81.23 to 81.35 • E; maximum elevation: 6710 m a.s.l.; Figure 1) situated between the 400 to 470 hPa pressure levels in the Westerlies, represents a favorable deposition site for dust transported long distances by large storms.
Geosciences 2022, 12, x FOR PEER REVIEW 3 of 29 snow of Mt.St. Elias (Yukon Territory, Canada), following a dust storm event that occurred in Asia and was transported across the North Pacific Ocean by clouds [34].
High elevation glaciers on the TP, dominating one of Earth's most arid regions, are repositories of tropospheric microparticles that have experienced mixing and gravitational sorting during transport.Even though local (i.e., coarser grain size fraction) soil dust is a non-negligible source, the particles entrapped in these high-altitude ice fields are substantially sorted by grain size with longer range transport enhancing sorting, making long range transported dust a prevalent deposit on TP glaciers [33,35].The major dust storm events, during which particles rise above the boundary layer and are then transported for thousands of kilometers, are most important to climate systems.Therefore, the high altitude Guliya ice cap (western Kunlun Mountains, 35.12 to 35.22° N, 81.23 to 81.35° E; maximum elevation: 6710 m a.s.l.; Figure 1) situated between the 400 to 470 hPa pressure levels in the Westerlies, represents a favorable deposition site for dust transported long distances by large storms.A growing body of research has characterized the Sr-Nd isotopic signature of Chinese [26,31,[36][37][38] and Indian [39,40] deserts' aeolian material, making it possible to identify Asian sub-regions, or desert sources, of dust preserved in Tibetan glaciers.The remote continental location of Tibetan glaciers ensures that the marine influence on their dust Sr-Nd isotopic composition is negligible compared to that of polar ice sheets [22].However, riverine and lacustrine contributions to the Sr isotopic composition of glacial dust remains possible through the deposition of carbonate from paleo-lacustrine sediments onto the glaciers [41,42].These evaporitic minerals containing calcium carbonate (CaCO3) have lower 87 Sr/ 86 Sr ratios than those associated with aluminosilicates and can impede the accurate isotopic fingerprinting of the dust source.For this reason, the carbonate fraction is commonly removed from the ice core insoluble mineral dust by weak acid leaching.Despite their potential as archives of Asian dust provenance, to date, only three shallow ice cores from the Tibetan Plateau (Dunde ice cap, [35]), Tien Shan (Miaoergou glacier, [43]) A growing body of research has characterized the Sr-Nd isotopic signature of Chinese [26,31,[36][37][38] and Indian [39,40] deserts' aeolian material, making it possible to identify Asian sub-regions, or desert sources, of dust preserved in Tibetan glaciers.The remote continental location of Tibetan glaciers ensures that the marine influence on their dust Sr-Nd isotopic composition is negligible compared to that of polar ice sheets [22].However, riverine and lacustrine contributions to the Sr isotopic composition of glacial dust remains possible through the deposition of carbonate from paleo-lacustrine sediments onto the glaciers [41,42].These evaporitic minerals containing calcium carbonate (CaCO 3 ) have lower 87 Sr/ 86 Sr ratios than those associated with aluminosilicates and can impede the accurate isotopic fingerprinting of the dust source.For this reason, the carbonate fraction is commonly removed from the ice core insoluble mineral dust by weak acid leaching.Despite their potential as archives of Asian dust provenance, to date, only three shallow ice cores from the Tibetan Plateau (Dunde ice cap, [35]), Tien Shan (Miaoergou glacier, [43]) and the Himalaya (East Rongbuk glacier, [44]), referred to collectively as the Third Pole region, have been analyzed for Sr and Nd isotopes.Other Third Pole glacier dust provenance studies using TEs and REEs have been conducted on surface snow or samples from shallow ice cores, and thus provide information on modern aerosols only [33,35,[43][44][45][46][47][48][49][50].A deeper temporal perspective on aeolian dust composition and source is lacking for Tibet, especially in the northwestern regions.The first geochemical investigation of a Central Asian ice core was conducted by [45] on a shallow core from Inilcheck glacier in Tien Shan (Kyrgyzstan), which was found to accumulate loess during deposition events.Calcite and gypsum originating from various sites constituted the average background mineralogical composition of Inilcheck dust.Ref. [35] examined the isotopic signature in discrete samples from a shallow ice core from the Dunde ice cap (Qilian Mountains, Northern China) and found 87 Sr/ 86 Sr and ε Nd values within the same range as those of the Taklimakan and Qaidam sands.They conclude that, under the present climate, Dunde dust originates from the Tarim and the Qaidam Basins, where it is entrained and transported over long distances by westerly winds.Similarly, for the Miaoergou glacier (Tien Shan, north of Tarim Basin), an isotopic investigation was conducted on six discrete contemporary ice core samples (dated between 1956 to 2000 AD) and conclusions were drawn from a relatively short time interval [43].Using REE and Sr-Nd isotopes, these authors documented that since 1956 AD the Taklimakan and the Gobi deserts were the main sources of aeolian dust deposited on Miaoergou.These findings were corroborated by [33] using the Sr and Nd isotopic compositions of dust from cryoconites at the surface of Miaoergou glacier.
In this study, the mineralogical, elemental, and radiogenic isotopic signatures of the aeolian dust sampled from discrete deep sections of two ice cores (retrieved from the Summit and Plateau sites) from the Guliya ice cap are characterized to assess their potential in tracing the long-term variability of paleo-hydro-climatic conditions prevailing in NW Tibet, and to identify the lithological and geographical source area(s) for the dust.The sampling approach was guided by visual characteristics of dust in the ice (e.g., color, layer feature, orientation).We also investigated discrete dust samples in ice sections preceding and succeeding several stadial/interstadial transitions (i.e., the cold/warm transitions within the Last Glacial Period) recorded in the 1992 Guliya ice core [11], for which we expect to observe a pronounced variability in the dust parameters.The results are interpreted considering the compositions of aeolian deposits discussed in the literature.

The Guliya Ice Cap and Surrounding Deserts
The Guliya ice cap is a ~200 km 2 ice field located in the Kunlun Mountains and situated under the atmospheric dominance of the continental westerly flow which is particularly intense during winter (Figure 1).Indian Summer Monsoon incursions can occasionally and indirectly modulate the atmospheric circulation pattern over NW Tibet (Figure 1).This location makes the Guliya site particularly relevant for the study of the Westerlies' intensity, which can entrain and transport mineral aerosols from northern (e.g., Gurbantünggüt, Taklimakan) and western deserts (e.g., Tajik loess desert) to the very high-altitude surface of the glaciers in the region and to the free troposphere (>6000 m a.s.l.).Ultimately, transport can be achieved over long distances to Eastern China and the remote Pacific region [51,52].Additionally, with the influence of the Indian Summer Monsoon and a more southwestern circulation during warmer periods, the Guliya ice cap could potentially be the sink for mineral dust transported from the sandy Sahara, the Arabian Peninsula, the Thar Desert in Rajasthan [39,40], and the southern TP.
The Kunlun Mountains ring the southern edge of the Taklimakan sandy desert in the Tarim Basin and circumscribe the arid northwestern TP region.Moraines and loess deposits are present on the northern slopes of the Kunlun Mountains and represent a proximal source of mineral aerosols to the Guliya glacier.The loess deposits in the Kunlun Mountains are among the highest in the world.They are derived from the Taklimakan Desert dust entrained by surface northeasterly winds and deposited on the windward slopes of the Kunlun Mountains [52].Tibetan Plateau soils and Quaternary deposits (e.g., talus fans, alluvial fans, and paleo-lake beds) on the southern side of the Kunlun Mountains also represent potential sources of aeolian silt particles to the Guliya ice cap [53].The vast desert areas surrounding the Guliya ice cap and the very low snow accumulation on the ice cap (230 mm w.e.yr −1 at the Plateau site, [54]) favor very high dust concentrations in the Guliya ice core (1650-1992 A.D. median dust concentration = 194 × 10 3 particles mL −1 , [55]).The high dust content of the core presents a great advantage in terms of sampling resolution for geochemistry as compared to cleaner polar ice cores for which more voluminous ice samples are required to collect a sufficient mass of insoluble material for the Sr and Nd isotope measurements (typically 200 ng for Sr and 50 ng for Nd).Arid conditions at the Summit of the Guliya ice cap are further exacerbated by very low snow precipitation relative to sublimation.The average annual ice layer thickness at the Summit (125 mm ice eq.yr −1 ) is less than half of that on the Plateau (262 mm ice eq.yr −1 ), a lower site on the ice cap [54].  3Cl data from a deep ice core extracted from the Plateau site in 1992 suggest that the deepest 20 m of GP may be older than 500 kyr [11].If correct, GP would be the oldest ice core archive yet recovered from the mid-latitudes.

The Guliya
The top 110 kyr, which is the better chronologically constrained part of the GP core, is the only mid-latitude climate archive comparable to the deep Greenland and Antarctic ice cores in terms of age, dating precision, and diversity of paleo-environmental information [56,57].For instance, the δ 18 O curve of the 1992 Plateau core indicates that temperature fluctuations during the last 110 kyr (Figure 2) were consistent with glacial cycles recorded in Northern Hemisphere ice cores and North Atlantic Ocean sediments [11,58,59].The chronology of the 2015 GS core also is under development, the discussion of which is beyond the scope of this paper.However, to provide some context to the work presented here, we note that the bottom layers of GS extend into the deglaciation (~16 kyr).The paleo-environmental interpretation of the GP and GS dust geochemistry presented here does not take the sample's age into consideration.For this reason, and despite the reliable original timescale developed for the upper 200 m of the 1992 GP core [11], the results of this study are discussed in terms of sample depth rather than sample age (Figure 2).

Stratigraphic Observations
During the processing of the ice core sections in the freezer facilities of the BPCRC, a stratigraphic description of GS and GP allowed the distinction of different types of dust within both ice cores.The most illustrative example is a progressive transition with depth from the light brown silt-sized particles to the grey nodules dominating the bottom 20 m of GP core.Both types of dust are observed in ice layers around 300 m depth (Figure A1).These prominent features in the ice core stratigraphy suggest that the source area for dust, the transport pathways of mineral aerosols, the climate at the dust source, or englacial conditions, have changed markedly over the time period represented by the Guliya ice cores.This was explored further by characterizing mineralogical, elemental (TE and REE), and Sr-Nd isotopic signatures, and thus aiding the identification of the geographical source(s) of the different types of dust deposited on the Guliya ice cap.For the GS core, the deeper layers with more dust layers visible were a natural first choice for sampling and allowed a preliminary comparison of the bottom ice layers at the two sites (Summit and Plateau).

Sample Preparation and Analysis
Two types of samples are used in this study from both GP and GS: (1) dust particles collected on filters from melted ~1 m ice core sections, hereafter referred to as "dust filters" (Section 2.3.1), and (2) acid-leached ice core subsamples from ~1 m ice sections, hereafter referred to as "ice samples" (Section 2.3.2) (Figure A2).Subsamples refer to the ~2 to 10 cm long individual ice pieces cut along depth from each ~1 m long ice sections (Table 1).[11] illustrating the climatic periods recorded in the Guliya Plateau (GP) ice core (lower panel) as defined by [11].The vertical bars mark the 1 m ice core sections sampled for microparticle concentration and size, and trace element concentrations in ice corresponding to "warm" (purple/red) or "cold" (pink/blue) periods.The bars are increasingly wider with age (and depth) as 1 m of ice comprises more time toward the bottom of the glacier due to the increased thinning of successive ice layers with depth.Note that the vertical bars representing T290 and T292 (brown), and T305 and T306 (grey) are juxtaposed and therefore appear as two bars instead of four on the time axis.The crossed circles indicate the depth at which dust filtrate from 1 m ice sections were collected on filters (i.e., dust filter samples) for geochemical and isotopic analyses (see Table 1).

Stratigraphic Observations
During the processing of the ice core sections in the freezer facilities of the BPCRC, a stratigraphic description of GS and GP allowed the distinction of different types of dust within both ice cores.The most illustrative example is a progressive transition with depth from the light brown silt-sized particles to the grey nodules dominating the bottom 20 m of GP core.Both types of dust are observed in ice layers around 300 m depth (Figure A1).These prominent features in the ice core stratigraphy suggest that the source area for dust, the transport pathways of mineral aerosols, the climate at the dust source, or englacial conditions, have changed markedly over the time period represented by the Guliya ice cores.This was explored further by characterizing mineralogical, elemental (TE and REE), and Sr-Nd isotopic signatures, and thus aiding the identification of the geographical source(s) of the different types of dust deposited on the Guliya ice cap.For the GS core, the deeper layers with more dust layers visible were a natural first choice for sampling and allowed a preliminary comparison of the bottom ice layers at the two sites (Summit and Plateau).[11] illustrating the climatic periods recorded in the Guliya Plateau (GP) ice core (lower panel) as defined by [11].The vertical bars mark the 1 m ice core sections sampled for microparticle concentration and size, and trace element concentrations in ice corresponding to "warm" (purple/red) or "cold" (pink/blue) periods.The bars are increasingly wider with age (and depth) as 1 m of ice comprises more time toward the bottom of the glacier due to the increased thinning of successive ice layers with depth.Note that the vertical bars representing T290 and T292 (brown), and T305 and T306 (grey) are juxtaposed and therefore appear as two bars instead of four on the time axis.The crossed circles indicate the depth at which dust filtrate from 1 m ice sections were collected on filters (i.e., dust filter samples) for geochemical and isotopic analyses (see Table 1).
Five dust filters (4 from GP and 1 from GS), and sixteen ice samples (14 from GP and 2 from GS) were selected for elemental geochemistry analysis (Table 1, Figure 2).The GP samples range in age from the Holocene to the Eemian (and possibly older).The bulk 87 Sr/ 86 Sr and 143 Nd/ 144 Nd isotopic compositions of these five dust filters were determined to constrain the geographical source of the brown and the grey dust observed at different depths of GS and GP.Two additional dust filters (T50 from GS and T309 from GP, Table 1) were analyzed for bulk mineralogy only, to also help elucidate provenance.It is important to note that for this preliminary study, collecting a sufficient amount of insoluble material and being able to readily compare the Guliya samples to those from the Dunde ice cap (then the only other Tibetan ice core dust analyzed for Sr-Nd isotopes [35]) were two primary concerns.Consequently, dust has not been treated with any acid to remove the carbonate fraction, nor has it been size-sorted to target the clay-size fraction.Therefore, the results discussed here represent the bulk of the dust collected on filters.The entire GS and GP ice cores, including the ice samples used for this study, were previously analyzed for oxygen isotopes (δ 18 O), major ions, and insoluble microparticles (size and concentration) in the ICPRG's Class 100 Clean Room.However, the full sample resolution records of these data will not be presented in this manuscript.Measurements of particle size distribution and concentration were performed with a Multisizer IV © Coulter counter (aperture diameter: 30 µm; 14 size bins: 0.63 to 16 µm, [54]) on diluted aliquots of melted ice samples.Note that each dust filter and ice sample ID number (T-XXX) refers to the cylinder (or "Tube") in which the corresponding ~1 m ice section was stored following extraction from the glacier.Therefore, the sample ID number is similar, but not identical to, the sample's average depth (Table 1).

Dust Filters
A diluted portion of the melted ice samples was consumed by the Coulter counter analyses [54].From the remaining undiluted portion, particles were filtered onto a cellulose Millipore filter (pore diameter 0.22 µm), stored in a Petri dish, and each filter was kept in the Class 100 clean room until analysis.Each filter collected the bulk fraction (>0.22 µm) of insoluble material contained in 1 m of ice core.The different types and colors of dust observed in the layers of the ice cores could still be visually distinguished among the filters (Figure A1 and [60] there in).Although the ages of the samples are not considered in this study, it is important to note that the increased layer thinning with depth makes one meter of ice comprise a longer time period as it is sampled deeper along the core.Conversely, in the upper part of the ice cores, time is relatively more stretched, making transitions between dust types more "diluted" in time, and thus less conspicuous.In this preliminary study, the four dust filters from GP correspond to top depths (m) of 223 (T225), 232 (T233), 290 (T292) and 305 (T308), and the dust filter from the GS core corresponds to 48 m depth (T49) (Table 1, Figure 2).

Mineralogy
The two additional dust filters, T50 from GS and T309 from GP were first analyzed directly using powder X-ray diffraction (XRD) to determine qualitative bulk mineralogy of all grain sizes present in the samples.Filter mounts were prepared by adhering them with double-sided adhesive tape to glass slides and mounting the slides onto XRD sample holders.Small portions of each filter were used to prepare samples for scanning electron microscopy and energy dispersive X-ray spectroscopy (SEM-EDXS) by pressing aluminum stubs with carbon adhesive tape directly onto the filter samples.These first steps are important for identifying trace minerals in the filters that may be sensitive to water or alcohol washing, which are required for the separation of grains that include only the aeolian size fraction of material transported to the Guliya ice cap.Because the filter mounts produced a high background in XRD scans that made data interpretation difficult, a second set of mounts were prepared by isopropyl washing of dust from the filters onto zerobackground holders.The resulting XRD scans from this second sample preparation method yielded better results (see Section 3.1.3).XRD allows identification of major and minor minerals present, while SEM-EDXS yields complementary micro-analytical data to support mineral identification by XRD, to detect trace minerals, and to provide morphological information.Particular attention was focused on the clay minerals (commonly found within the fine fraction, <2 µm) which tend to adsorb trace metals (e.g., [61]).

Sr-Nd Isotope and Elemental Geochemistry Analysis by MC-ICPMS and ICPMS
Particulates were remobilized from each dust filter and then centrifuged in a conical tube with ultrapure water.After centrifugation (30 min), the clear supernatant was discarded using a syringe, and the remaining water was evaporated in a class 100 laminar flow box.Extracted particles were dissolved in an HF-HNO 3 mixture in Teflon hexvials that were heated to 100-120 • C for three days.The full acid digestion allows the determination of total elemental concentrate ons in the dust filter samples (contrary to a fraction of the total TEs determined by acid leaching of ice samples (see Section 2.3.2)).Following sample dissolution, Sr and REEs were isolated on 3 mL quartz columns using conventional ion-exchange chromatography primary columns packed with AG 50 W-X12 (200-400 mesh) resin.Strontium was extracted with 3.5 N HCl and the REE fraction was collected in 6 N HCl.Neodymium (Nd) was then separated from the other REEs on Teflon columns using 800 µL of LnSpec resin as the cation exchange medium.The Sr-Nd isotopes were analyzed on a Nu Plasma Multi-Collector Inductively Coupled Plasma Mass Spectrometer (MC-ICPMS) at the University of Florida (UF).Sr isotopes were analyzed using wet plasma and the time resolved analysis (TRA) method developed by [62].The long-term average 87 Sr/ 86 Sr value of TRA-measurements of NIST 987 (reference value: 87 Sr/ 86 Sr = 0.710248, [63]) at UF is 0.710244 (±0.000024, 2σ) and was used to normalize the measured data.Procedural blanks were less than 100 pg Sr or four orders of magnitude less than sample concentrations.Samples for Nd isotopic analyses were also analyzed using the TRA method.Samples were aspirated through a DSN-10 desolvating nebulizer system.The average of measured 143 Nd/ 144 Nd ratios for the standard JNdi-1 (reference value: 143 Nd/ 144 Nd = 0.512115, [64]) is 0.512115 (±0.000014, 2σ) and was used to normalize the 143 Nd/ 144 Nd isotopic ratios of samples.The long-term 2σ external reproducibility of JNdi-1 analyses (±0.000014) is equivalent to ±0.27 ε Nd units, where ε Nd is calculated using the equation: ε Nd (0) = (( 143 Nd/ 144 Nd) sample /( 143 Nd/ 144 Nd) CHUR − 1) × 10000 (1) and the isotopic composition of the Chondritic Uniform Reservoir (CHUR) is ( 143 Nd/ 144 Nd) CHUR = 0.512638 [65].Procedural blanks were approximately 15 pg Nd or four orders of magnitude less than sample concentrations.The determination of the geochemical composition of the dust filters also entailed the analysis of 19 TEs (Li, Ti, V, Cr, Co, Ni, Cu, Zn, Ga, Rb, Sr, Zr, Nb, Cs, Ba, Nd, Ta, Pb, U) and 16 REEs (La, Ce, Pr, Nd, Sm, Eu, Gd, Tb, Dy, Ho, Er, Tm, Yb, Lu) including Sc and Y on a Thermo Scientific Element 2 ICPMS.Samples were diluted 1:1 with 5% HNO 3 spiked with 8 parts per billion (ppb) rhenium (Re) and rhodium (Rh) and analyzed using in-house standards.Long-term measurement reproducibility was checked against the AGV1 standard.Relative standard error between repeated measurements of AGV1 was less than 8% for all elements.

Ice Samples
The 16 selected ice samples (two from GS and 14 from GP) were cut into a total of 158 subsamples (5 to 10 cm long) with a carbon steel band saw and were transferred to a Class 100 clean room where they were decontaminated by triple washing with ultra-pure water.The ice samples were melted at room temperature in pre-cleaned LDPE containers and acidified with ultra-pure HNO 3 (optima grade) to a final concentration of 2% (v/v).Following acidification, the samples were kept in the clean room for 30 days to allow the complete leaching of TEs by the acid, which releases a fraction of the total TEs in the sample [66].The subsequent analysis by ICPMS (Thermo Scientific Element 2 housed at BPCRC) coupled with a desolvation system for sample introduction (Apex Q, ESI), allowed the determination of 29 major and trace elements (Li, Na, Mg, Al, Ti, V, Cr, Mn, Fe, Co, Ni, Cu, Zn, Ga, As, Rb, Sr, Nb, Mo, Ag, Cd, Sn, Sb, Cs, Ba, Tl, Pb, Bi, and U).All concentrations in the Guliya samples were several times higher than the instrumental detection limits determined by [55,67].The measurements' accuracy and precision values were similar to those reported by [55,67]: measured concentrations of the TMRain-95 standard solution fall within the uncertainty limit of the certified concentrations and the precision is within 4-10% for all TEs except Nb and Sn (17%), and Ag (22%).

Enrichment Factor
The crustal enrichment factor (EF) is defined as the ratio of the concentration of a given trace element (TE) to that of an element (X) that is stable and abundant in rock and soil, typically Fe, Al, or Ti.This ratio (TE/X) spl is normalized to the same ratio obtained from a crustal reference (TE/X) ref : EF = (TE/X) spl /(TE/X) ref (2) The ice samples and dust filters selected for this preliminary study have been sampled from layers that are too old to have been contaminated by the deposition of anthropogenic emission products.Therefore, instead of highlighting the anthropogenic contribution to the sample's elemental budget, the EF is used here to extract subtle compositional differences among samples.The preparation (i.e., acid leaching) and analysis by ICPMS of the Guliya ice samples followed the same procedures as those employed by [55,68] to characterize the PSA soil samples collected from the vicinity of Guliya glacier.Consequently, this set of PSA soils is the preferred crustal reference used in the calculation of EF for the ice samples.However, for the digested dust filter samples, the Upper Continental Crust (UCC) composition [69] is the chosen crustal reference because UCC concentrations were determined by X-ray fluorescence spectrometry and atomic absorption spectroscopy [70,71], which provides total elemental concentrations comparable to those obtained after the full acid digestion of the dust filter samples.The dust concentration shown in Figure 3 corresponds to the median concentration of discrete samples or "subsamples" (n: number of discrete subsamples; Table 1) analyzed in each ice sample (i.e., each ~1 m ice section) of interest.Samples from GS display higher dust median concentrations than any of the GP ice samples investigated, except for T309 in the deep GP layers.This result is expected considering the lower annual snow accumulation at the Summit compared to the Plateau site.However, because the age of the GS samples is not yet determined precisely, their dust concentration cannot be readily compared to those of the GP samples.Generally, in GP, higher dust concentrations are observed in ice layers corresponding to cold periods when compared to warm periods occurring before and after them, which corroborates findings from the Guliya 1992 core [11,72].This variability observed on an orbital timescale is comparable to that on a seasonal timescale for Guliya, where the highest dust storm frequency occurs during the coldest and driest season (i.e., late winter-early spring), when the Westerlies intensify [73,74].An enhanced westerly circulation under a dry climate favors the transport of mineral particles over long distances, their size sorting while in route, and ultimately, their deposition onto high-altitude glaciers.This process likely contributed to the higher proportion of fine particles (<4 µm) observed in the GP samples for cold periods compared to those from warm periods (Figure 4a), even if both types of samples have a similar diameter size mode (8.4 µm (T225) and 8.6 µm (T233), respectively, Figure 4b).The more extensive snow cover around the Guliya ice cap during cold periods would inhibit the entrainment of dust from local sources, thereby further increasing the relative contribution of distant dust sources and finer dust.

Results and Discussion
Towards the bottom of GP, the ice is compressed, and 1 m of ice comprises more time, thereby smoothing out the changes in dust concentration between warm and cold periods (e.g., T262 to T278 during MIS5; Figure 3).For instance, this is characterized by the smaller 75th percentile of samples T273 (9 × 10 5 part.mL −1 ) and T292 (8.7 × 10 5 part.mL −1 ) compared to T308 (5.3 × 10 6 part.mL −1 ).The largest variation in particle counts in the deepest part of the core coincides with the stratigraphic transition from mostly diffuse brown dust to mostly grey dust.The particle count is high when the grey aggregates (or nodules) occur in a sample taken from an ice section otherwise described as "clean" relative to the dustier upper layers of the core (Figure A1).The distribution of grey nodules is heterogeneous within the "clean" bottom ice layers.Therefore, the bottom GP ice samples can either contain grey nodules with numerous microparticles or are "clean ice" without nodules and yielding a low dust concentration.Towards the deepest part of GS, the volume size distribution of particles measured in T49 is bimodal with a primary mode in the fine grain silt-size at 10.6 µm (11.4% of dust volume), and a secondary mode in the clay-size range at 1.8 µm, which represents only 7.22 % of the dust volume in that sample (Figure 4b).The difference in particle size distribution between the bottom of GS (bimodal distribution) and the bottom of GP (predominantly clay-size fraction) is one of the most striking observations of spatial variability in the Guliya dust record.A wide range in grain sizes is observed on the T50 dust filter (Figure 5).Comparing the deepest layers of the two Guliya core sites, Figure 5 shows that the dust deposited on T50 (~48 m depth on GS) contains a higher proportion of silt and fine sand particles than on the plateau of the ice cap (i.e., T309 of GP).The subangular grain shape observed with SEM (Figure 5a), and the reduced size sorting suggest a more proximal source of aeolian sediments, or a mix of particles from proximal and remote sources.The angularity of the larger grains also supports a glacial origin.The occurrence in T50 of Ca-pyribole (pyroxene or amphibole), clinochlore (Mg-rich chlorite), subhedral calcite (Figures 5 and 6) and halite suggests relatively weak chemical weathering in the dust source region.This

Mineralogy of the Deep GS Dust Filter (T50): Dry Conditions and Proximal Influence
A wide range in grain sizes is observed on the T50 dust filter (Figure 5).Comparing the deepest layers of the two Guliya core sites, Figure 5 shows that the dust deposited on T50 (~48 m depth on GS) contains a higher proportion of silt and fine sand particles than on the plateau of the ice cap (i.e., T309 of GP).The subangular grain shape observed with SEM (Figure 5a), and the reduced size sorting suggest a more proximal source of aeolian sediments, or a mix of particles from proximal and remote sources.The angularity of the larger grains also supports a glacial origin.The occurrence in T50 of Ca-pyribole (pyroxene or amphibole), clinochlore (Mg-rich chlorite), subhedral calcite (Figures 5 and 6) and halite suggests relatively weak chemical weathering in the dust source region.This implies that the climate was particularly dry or that a stronger insolation would have favored the sublimation of the ice surface and crystallization of salts.Illitic clay dominates the clay-sized fraction of T50, whereas quartz, K-feldspar and biotite compose the bulk of the silt fraction, which is typical of a granitic mineralogy.The size distribution and its mineral composition make the dust in the deep layers of the GS core resemble that of the Taklimakan desert sands, which among central Asian deserts, contain the highest amount of calcite [36,75,76].
Geosciences 2022, 12, x FOR PEER REVIEW 12 of 29 implies that the climate was particularly dry or that a stronger insolation would have favored the sublimation of the ice surface and crystallization of salts.Illitic clay dominates the clay-sized fraction of T50, whereas quartz, K-feldspar and biotite compose the bulk of the silt fraction, which is typical of a granitic mineralogy.The size distribution and its mineral composition make the dust in the deep layers of the GS core resemble that of the Taklimakan desert sands, which among central Asian deserts, contain the highest amount of calcite [36,75,76].Geosciences 2022, 12, x FOR PEER REVIEW 12 of 29 implies that the climate was particularly dry or that a stronger insolation would have favored the sublimation of the ice surface and crystallization of salts.Illitic clay dominates the clay-sized fraction of T50, whereas quartz, K-feldspar and biotite compose the bulk of the silt fraction, which is typical of a granitic mineralogy.The size distribution and its mineral composition make the dust in the deep layers of the GS core resemble that of the Taklimakan desert sands, which among central Asian deserts, contain the highest amount of calcite [36,75,76].

Mineralogy of the Deep GP Dust (T309): Grey Clay from a Weathered Source
The GP ice section T309 corresponds to a depth where small (1-5 mm) to large (<1.5 cm) scattered nodules of grey dust are the only type of dust visible to the unaided eye (Figure A1).Upon melting the ice samples, these nodules disintegrated into very fine particles, yielding the highest microparticle concentration (median concentration = 1.6 × 10 6 ± 2.10 7 part.mL −1 , Figure 3) and the finest particle size distribution (diameter median = 5.9 µm; diameter mode = 3.8 µm) among ice samples analyzed by Coulter counter (Figure 4).
The SEM analysis performed on the T309 dust filter (Figure 5b) confirms that the grey dust is mostly composed of clay-sized grains (<2 µm), with fewer silt-sized particles (<50 µm) of K-feldspar and bio-pyribole (biotite, amphibole, or pyroxene).The small size fraction could indicate a more efficient mineral sorting during long-range atmospheric transport of the microparticles to the ice cap, suggesting that the grey clay originates from a distant location.However, within the ice, the grey clay is observed as nodules.In a relatively wet environment, clay aggregates could have been entrained and transported from proximal clay deposits surrounding the glacier.Alternatively, nodules could be post depositional features.In some sections, dust nodules are concomitant to diffuse horizontal dust layers, which could argue that the grey clay was deposited already in the nodular form, or that the diffuse layers result from the smearing of a nodule along a shear plane.The aggregation of clay particles with time and depth [60] is another possible process that could explain the presence of the clay nodules (Figure A1).
The X-ray scan of T309 (Figure 6) indicates that the grey clay is composed of illite, mica (biotite/muscovite), chlorite and vermiculite.SEM imagery (Figure 5b) shows that biotite is occasionally physically weathered (rounded folia).Compared to other filter samples, the grey and green color of the dust layers visible in the deep sections of the GP core could reflect the occurrence of more ferrous (Fe 2+ )-bearing minerals such as Fe-Ti oxides (Figure 5b) [77], or a larger amount of chlorite or other Fe-bearing phyllosilicates.
Heavy REE-bearing phosphates, such as xenotime and monazite, are also found in the silt and clay-sized fraction (1-2 µm) of the bulk T309 sample, identified by EDX microanalysis (Figure A3).Potentially REE-bearing coarse fluorapatite minerals (20 µm) found in T309 are a common accessory mineral in various rocks, nevertheless they could also possibly indicate the hydrothermal nature of the dust source.Clay to fine silt-sized needles of iron and titanium oxides (ilmenite, spinel, rutile or anatase), which are resistant to weathering, are detected in higher proportions than in other samples (Figures 5b and 6).Titanium oxides are derived from the weathering of igneous and metamorphic rocks and sediments.They can also form in soils from the weathering of less-resistant minerals (e.g., biotite, amphibole) [78].The higher occurrence of small titanium oxide needles or equant-shaped grains in the deepest GP ice layers provides additional evidence that the source of aeolian particles to NW Tibet was subject to more substantial weathering prior to 110 kyr BP.

Trace Element Geochemistry
Dust filters.Elemental EFs for dust filter samples are shown in Figure 7. EFs do not covariate with the concentration of the crustal element of reference Ti in the filter samples (Table A1) and thus the influence of Ti on EF can be ruled out.The Sr EF in GS T49 (EF UCC Sr = 0.74) is notably higher than those of the four GP filters analyzed (EF UCC Sr = 0.22 (T225); 0.35 (T233); 0.27 (T292) and, 0.22 (T308), Table A1).This is consistent with the previously described mineralogical composition of the adjacent T50 dust filter, which includes calcite and halite.As strontium abundance in Chinese loess is controlled by carbonate [79], the higher Sr enrichment in T49 is likely associated with the larger amount of calcite in the dust, indicative of mixing with calcite-rich Taklimakan desert dust.The higher Sr EF in dust from "warm" intervals (T233) compared to that of dust filters from "cold" GP sections (T225) could also reflect a larger input of calcite-rich Taklimakan desert dust during warm periods.
pegmatites and released in geothermal waters and sediments ( [87] and reference therein).For example, [87], and [88] studied thermal springs and salt lakes in Tibet and found that many geothermal fields produce B, Li, and Cs-rich waters.Northern China is one of the regions of the world with the highest fluoride concentrations in its groundwater [85] (Figure A4).Thus, it is possible that the Cs, F, and Li-rich clay found in the deep part of the GP core originated from clay sediments of salt lakes or rivers fed by geothermal waters [89].For the two types of dust in GP, the elemental analysis of the brown (T292) and the grey (T308) dust filters shows that the former is enriched in Sr; while the latter is enriched in Li, Rb, Cs and Pb, similar to Southern Tibetan loess [80,81] and to Tibetan geyserite deposits [82].The high abundance of Pb in the grey dust could suggest an affiliation of Pb to phosphate minerals, zircons or to K-feldspar, which is commonly found in NW China granites for instance [83].Particle size measurements previously showed that the grey dust microparticles are much smaller (average particle diameter <2 µm) than those of the brown dust (Figures 4 and 5).The T308 enrichment in Cs (EF UCC Cs = 6) relative to the other filters is consistent with the presence of grey clay at the bottom of the core because clay minerals, such as vermiculite, tend to adsorb and sequester Cs + and Rb + ions [80,84].This adsorption mechanism by clay minerals and the release of ions in the melted ice sample, could also partially explain the higher F − concentrations measured in Guliya ice samples (Figure A4 and [85] there in) with high content of grey clay particles compared to those with brown dust [86].Studies of trace elements in granites show that Cs, Li, and F (among others) can be expelled at the end of crystallization of granitic magma to form pegmatites and released in geothermal waters and sediments ( [87] and reference therein).For example, refs.[87,88] studied thermal springs and salt lakes in Tibet and found that many geothermal fields produce B, Li, and Cs-rich waters.Northern China is one of the regions of the world with the highest fluoride concentrations in its groundwater [85] (Figure A4).Thus, it is possible that the Cs, F, and Li-rich clay found in the deep part of the GP core originated from clay sediments of salt lakes or rivers fed by geothermal waters [89].
Ice samples.Elemental analyses performed on GS ice samples (Figure 8a) and the GP "warm" (Figure 8b) and "brown dust" (Figure 8c) ice samples confirm the relative enrichment in both Sr and soluble elements such as Na, Nb, Mo, Sn and U (Table A2).Considering the high solubility of Mo in oxic systems and under alkaline conditions in calcareous soils [90], its enrichment in "warm" (i.e., 18 O-enriched) and "brown dust" ice samples (EF PSA Mo = 1.87 and 1.35, respectively) would indicate that Mo is desorbed from alkaline dust contained in the ice.The Mo desorption from the dust and its dissolution in solution may have occurred during the melting phase of the sample preparation (see Section 2.3.2.).Conversely, the noticeably lower Mo enrichment observed in "cold" (i.e., 18 Odepleted) (EF PSA Mo = 1) and "grey dust" ice samples (EF PSA Mo = 1.12) (Figure 8b) suggests that these two types of samples contain acidic dust in which Mo is adsorbed more efficiently on mineral surfaces (e.g., Fe-oxides) [91], and thus is less easily released into the liquid ice sample.While Mo enrichment could reflect alkaline conditions in the ice, confirmation of Mo enrichment as a reliable proxy for dry climate conditions (with weaker chemical weathering) requires the analysis of sufficient additional Guliya ice core samples to achieve a statistically significant result.
(EFPSA Mo = 1.87 and 1.35, respectively) would indicate that Mo is desorbed from alkaline dust contained in the ice.The Mo desorption from the dust and its dissolution in solution may have occurred during the melting phase of the sample preparation (see Section 2.3.2.).Conversely, the noticeably lower Mo enrichment observed in "cold" (i.e., 18 O-depleted) (EFPSA Mo = 1) and "grey dust" ice samples (EFPSA Mo = 1.12) (Figure 8b) suggests that these two types of samples contain acidic dust in which Mo is adsorbed more efficiently on mineral surfaces (e.g., Fe-oxides) [91], and thus is less easily released into the liquid ice sample.While Mo enrichment could reflect alkaline conditions in the ice, confirmation of Mo enrichment as a reliable proxy for dry climate conditions (with weaker chemical weathering) requires the analysis of sufficient additional Guliya ice core samples to achieve a statistically significant result.
Enrichments in Ag and to a lesser extent Tl in the ice containing grey clay (Figure 8c) seem to be characteristic of Chinese aeolian soils, as they were previously observed by [68] and [55] in samples collected, respectively, across the TP and in the vicinity of the Guliya ice cap.The relatively higher enrichment in Cs of the ice sample with grey dust compared to the ice with brown dust confirms the results obtained for the filter analysis: the GP grey clay is composed of Cs-hosting minerals such as vermiculite (Figure 6) [84].  1O-depleted (cold) and 18 O-enriched (warm) GP sections, (c) GP ice sections containing "brown" or "grey" dust.  1O-depleted (cold) and 18 O-enriched (warm) GP sections, (c) GP ice sections containing "brown" or "grey" dust.EF is calculated using Fe and the local PSAs [55] as the crustal references.Here, Fe is the preferred crustal element of reference to illustrate the Ti enrichment variability among the various types of ice samples.
Enrichments in Ag and to a lesser extent Tl in the ice containing grey clay (Figure 8c) seem to be characteristic of Chinese aeolian soils, as they were previously observed by [55,68] in samples collected, respectively, across the TP and in the vicinity of the Guliya ice cap.The relatively higher enrichment in Cs of the ice sample with grey dust compared to the ice with brown dust confirms the results obtained for the filter analysis: the GP grey clay is composed of Cs-hosting minerals such as vermiculite (Figure 6) [84].

REE and Sr and Nd Isotopic Composition of Dust Filter Samples: Provenance Keys
The five dust filters were analyzed for REE and Sr-Nd isotopes (Table 1).Their signature for their chondrite normalized REE patterns and for two size-independent provenance tracers (i.e., δEu N and (La/Yb) N , [76]) are given, respectively, in Figures 9 and 10.Chondritic normalization has been chosen over shale (PAAS) normalization to create more spread-out patterns and more readable figures.The chondrite normalized REE distribution patterns are generally homogeneous (Figure 9 and [92,93] there in) with a negative Eu anomaly (δEu N , Table A3), a steep light REE enrichment (LREE includes La to Eu), and a relatively flat heavy REE profile (HREE includes Gd to Lu) characteristic of Chinese loess [36,76,94,95].The ratio of La N to Yb N is a measure of the overall enrichment in LREE (Figure 10).EF is calculated using Fe and the local PSAs [55] as the crustal references.Here, Fe is the preferred crustal element of reference to illustrate the Ti enrichment variability among the various types of ice samples.

REE and Sr and Nd Isotopic Composition of Dust Filter Samples: Provenance Keys
The five dust filters were analyzed for REE and Sr-Nd isotopes (Table 1).Their signature for their chondrite normalized REE patterns and for two size-independent provenance tracers (i.e., δEuN and (La/Yb)N, [76]) are given, respectively, in Figure 9 and Figure 10.Chondritic normalization has been chosen over shale (PAAS) normalization to create more spread-out patterns and more readable figures.The chondrite normalized REE distribution patterns are generally homogeneous (Figure 9 and [92,93] there in) with a negative Eu anomaly (δEuN, Table A3), a steep light REE enrichment (LREE includes La to Eu), and a relatively flat heavy REE profile (HREE includes Gd to Lu) characteristic of Chinese loess [36,76,94,95].The ratio of LaN to YbN is a measure of the overall enrichment in LREE (Figure 10).[36] and to (c) moraines of west Kunlun Mountains [92] represented by the shaded area.REE abundances in the C1 chondrite [93] are used for normalization.[76]; West Kunlun moraines [92]; and Taklimakan sand and West Taklimakan loess [36].The N subscript in δEuN and (La/Yb)N indicates that values are normalized to the C1 chondrite [94].
To track the geographical source of strontium and neodymium, 87 Sr/ 86 Sr and ƐNd results for the bulk Guliya dust are matched to values reported in the literature for silt-sized surface dust samplesthe last glacial period as reported by [37].However, until more evidence is gathered, the CLP situated due east of the Kunlun Mountains is not considered a potential source of aeolian dust to Guliya ice cap because wind patterns are assumed to be more westerly such that part of the CLP material originates from the NW Chinese deserts and the Kunlun moraines [36,92].Furthermore, drawing conclusions regarding the differential influence on the Summit versus the Plateau sites would require sampling the dust material in ice layers of the same age in both cores, which was not feasible for this study as the GS time scale was under development.When compared with data from other Chinese glaciers (Figure A5) and glaciers around the world (Figure A6, and [100] there in), the Sr-Nd isotopic composition of the T49 GS dust filter resembles that of modern dust collected from Tien Shan glaciers, which was found to originate from the northern boundary of the Taklimakan desert and, to a lesser extent, from the Gurbantünggüt desert [48].The Sr-Nd isotopic values for the GP samples (except for T308) are homogeneous and indicate that either (1) the mineral dust deposited on the Guliya Plateau was well-mixed in the atmosphere, or (2) the time encompassed by 1 m of ice in the deep layers of the Guliya glacier is so compressed that dust from multiple sources are contained in a single sample (filter).The potential multiple sources of GP dust could then be similar to those of the dust collected from the Miaoergou shallow core (Eastern Tien Shan), with isotope values overlapping those of Guliya (Figures A5 and A6) and that were interpreted as a mixture of Taklimakan and, Gobi deserts material [43].It is then possible that the GP dust also originates from these two sources even though it does not plot along mixing lines (not shown) between the fields defined by Taklimakan sand and loess and the Gobi Desert.
The GP samples differ from the Dunde ice cap's modern dust samples [35].ƐNd for T225 (−9.55 ± 0.23) is the closest value to the Dunde signature, suggesting that, under cold climate conditions, some northern desert material from the Tarim and/or Qaidam Basins could also be entrained on the surface of Guliya.The 87 Sr/ 86 Sr of the grey dust (T308) is  [94] and the Central Asian deserts: Chinese Loess Plateau (CLP) loess; Taklimakan fine sand; Northeast deserts (Mu Us, Hobq, and Badain Jaran deserts); Tibetan soils, and Thar desert [76]; West Kunlun moraines [92]; and Taklimakan sand and West Taklimakan loess [36].The N subscript in δEu N and (La/Yb) N indicates that values are normalized to the C1 chondrite [93].
Among the five dust filter samples, T233 ("warm" GP) exhibits the highest (La/Yb) N ratio (10.87) and the lowest δEu N (0.60), which implies that fewer HREE-carrier minerals (e.g., clays) and relatively more plagioclases compose the mineralogy of T233.δEu N for T233 (Figure 10) falls within range of Eu anomalies measured in CLP [76] and West Taklimakan loess [36], while its REE pattern (Figure 9) matches that of West Taklimakan loess and West Kunlun Moraine [92].This suggests that the T233 dust deposited during MIS 5a on the glacier plateau could be a mix of loess and moraine material originating west of the ice cap.However, caution is required when comparing the geochemistry of late Pleistocene Guliya dust samples with that of modern moraine sediments and surface loess deposits, as these are possibly much younger than the deeper layers of the Guliya ice cap.For instance, [80] reports the lack of loess deposits older than 13 kyr in southern Tibet due to their erosion and recycling during glacial-interglacial cycles of the Quaternary.Thick (up to 700 m) and old (~800 kyr) loess deposits border the northern flanks of the western Kunlun Mountains [96], but to our knowledge, no Sr-Nd isotopic data have been reported for the Late Pleistocene layers of these sequences.This supports the strong relevance of the oldest Guliya glacial dust as an invaluable aeolian archive in the region.The Summit dust (T49) displays the second highest enrichment in LREE with a (La/Yb) N ratio of 8.83 compared to 8.54 for T292, 8 for T225 and 6.93 for T308.Because LREEs are less stable during weathering [14,97], their slightly higher enrichment in T49 is consistent with the presence of soluble and weathering-sensitive minerals (e.g., calcite; see Section 3.1.2)found in T50.The limited chemical weathering that affected T49 dust indicates dry climatic conditions prevailing at the source and in the Guliya region at the time of its deposition onto the ice cap surface.Figure 9 provides additional information on the possible dust source.The HREE part of the T49 REE profile falls within range of that of the Taklimakan desert sand and the West Kunlun moraines (Figure 9a), while the LREE part and its moderate Eu anomaly (δEu N = 0.64) are more typical of the West Taklimakan loess (Figure 9b).Dust filter T308 ("grey dust" from GP) exhibits the least negative Eu anomaly (δEu N = 0.71), which could result from the combined influence of vermiculite possessing a positive Eu anomaly [14], and of the smaller amount of plagioclase detected in the grey clay (Figure 6).Eu-hosting plagioclase minerals are characterized by a negative Eu anomaly as they release Eu 2+ during weathering [98].Therefore, less plagioclase could weaken the negative Eu anomaly of the grey clay.The reduced L/HREE fractionation in T308 is also the smallest reported here (L/HREE = 2.61 and (La/Yb) N = 6.93,Table A3) and is linked to the relatively high abundance of weathering resistant phases such as Ti-oxides, zircons or phosphates in which HREE are preferentially present (Figure A3) [14,99].The REE abundances in T292 ("brown dust" GP) and in T225 ("cold" GP) are lower than those of T49, T233 and T308 (Figure 9).Their REE patterns are quasi-identical and their Eu anomalies (δEu N = 0.67 and 0.65, respectively) are typical of the average upper continental crust [94] and of fairly well mixed sediments.
To track the geographical source of strontium and neodymium, 87 Sr/ 86 Sr and ε Nd results for the bulk Guliya dust are matched to values reported in the literature for siltsized surface dust samplesthe last glacial period as reported by [37].However, until more evidence is gathered, the CLP situated due east of the Kunlun Mountains is not considered a potential source of aeolian dust to Guliya ice cap because wind patterns are assumed to be more westerly such that part of the CLP material originates from the NW Chinese deserts and the Kunlun moraines [36,92].Furthermore, drawing conclusions regarding the differential influence on the Summit versus the Plateau sites would require sampling the dust material in ice layers of the same age in both cores, which was not feasible for this study as the GS time scale was under development.When compared with data from other Chinese glaciers (Figure A5) and glaciers around the world (Figure A6, and [100] there in), the Sr-Nd isotopic composition of the T49 GS dust filter resembles that of modern dust collected from Tien Shan glaciers, which was found to originate from the northern boundary of the Taklimakan desert and, to a lesser extent, from the Gurbantünggüt desert [48].The Sr-Nd isotopic values for the GP samples (except for T308) are homogeneous and indicate that either (1) the mineral dust deposited on the Guliya Plateau was well-mixed in the atmosphere, or (2) the time encompassed by 1 m of ice in the deep layers of the Guliya glacier is so compressed that dust from multiple sources are contained in a single sample (filter).The potential multiple sources of GP dust could then be similar to those of the dust collected from the Miaoergou shallow core (Eastern Tien Shan), with isotope values overlapping those of Guliya (Figures A5 and A6) and that were interpreted as a mixture of Taklimakan and, Gobi deserts material [43].It is then possible that the GP dust also originates from these two sources even though it does not plot along mixing lines (not shown) between the fields defined by Taklimakan sand and loess and the Gobi Desert.
The GP samples differ from the Dunde ice cap's modern dust samples [35].ε Nd for T225 (−9.55 ± 0.23) is the closest value to the Dunde signature, suggesting that, under cold climate conditions, some northern desert material from the Tarim and/or Qaidam Basins could also be entrained on the surface of Guliya.The 87 Sr/ 86 Sr of the grey dust (T308) is unusually high ( 87 Sr/ 86 Sr = 0.72425, Table 2) and unmatched among Chinese glacier dust signatures reported to date.This reveals the uniqueness of the material preserved in the bottom of the Guliya Plateau core which could reflect its old age and the unknown environmental setting at that time in northwestern Tibet.As the Sr isotope ratio is easily altered by grain size sorting [101,102], it is also possible that the higher 87 Sr/ 86 Sr value partly reflects the finer grain size fraction composing the material in T308.The origin of this grey clay could be determined by matching its Sr-Nd signature to that of the finest fraction (< 5 µm) of potential source area samples (Figure 11) as T308 dust is naturally size-sorted to the clay fraction (Figure 4).Its signature falls into the same range as those of the clay-sized CLP loess deposited between 20 and 50 kyr BP [30] and of a surface sample from the Qaidam desert [37], which likely indicate a common origin.However, other sources for T308 cannot be excluded until the analysis of the Guliya dust is rigorously constrained to its silicate fraction and its clay-size fraction.We intend to apply those constraints to future isotopic analyzes of Guliya dust samples.clay could be determined by matching its Sr-Nd signature to that of the finest fraction (< 5 μm) of potential source area samples (Figure 11) as T308 dust is naturally size-sorted to the clay fraction (Figure 4).Its signature falls into the same range as those of the clay-sized CLP loess deposited between 20 and 50 kyr BP [30] and of a surface sample from the Qaidam desert [37], which likely indicate a common origin.However, other sources for T308 cannot be excluded until the analysis of the Guliya dust is rigorously constrained to its silicate fraction and its clay-size fraction.We intend to apply those constraints to future isotopic analyzes of Guliya dust samples.  8Sr/ 86 Sr of the five Guliya dust filters (squares) compared to Sr and Nd isotopic arrays of (<75 μm) dust sampled in Chinese deserts: Gurbantunggut, Qaidam [37], Gobi [26], Taklimakan sand and Tarim loess [36,37], Kunlun moraines [92], CLP loess [37], Northeast deserts including Mu Us, Hobq and Badain Jaran deserts [36]; and Indian deserts Thar [39] and Ghaggar river sediments [103].The inset identifies the dust filter samples (squares) shown in the main figure and individual data points for the Chinese deserts' arrays within the rectangle.Data points for the clay-size fraction of surface Qaidam sand (<5 μm, [37]) and last glacial CLP loess (<10 μm, [30]) are also shown (crossed circles) as they can be compared to the naturally size sorted T308 sample (fine grey clay).

GP Grey Clays: A Possible Biological Influence
The grey color of the clay trapped within the deepest layers of the GP ice core (T309) is associated with the presence of phyllosilicates (chlorite, biotite, vermiculite) and oxides containing reduced iron (Fe 2+ ), while the reddish-brown colors of the shallower diffuse dust layers are characteristic of sediments containing ferric (oxidized) iron (Fe 3+ ) [77].In pedology, a grey color is generally used as an indicator of anaerobic (oxygen-free) conditions, which occurs when soils are water-saturated (e.g., wetland) and anaerobic bacteria are present [104].During respiration, these microorganisms utilize as nutrients the oxidized forms of elements such as N, Fe, Mn, and S, which become chemically reduced.Therefore, nitrate (NO 3 − ), manganese (Mn 4+ ), iron (Fe 3+ ) and sulfate (SO 4 2− ) are depleted by anaerobic respiration and, respectively, converted to NH 4 + , Mn 2+ , Fe 2+ , and HS − .These additional indications that the grey dust originates from water-saturated soils or wetland soils could support our hypothesis of a wetter or a warmer climate at the time the deepest Guliya ice layers were deposited.Nonetheless, the post-depositional microbial activity could also account for the reduction of Fe in clay minerals within the ice.The fact that the grey clay nodules and the diffuse layers of brown dust are occasionally found within the same ice horizon (Figure A1), suggests that anaerobic bacteria are present on the grey clay particles only.A thorough microbiological analysis, like that performed by [105], combined with further mineral analyses of the oldest layers of the GP core would be instrumental in determining whether iron in the nodular grey clay has been reduced in situ by anaerobic microorganisms.No melt features (e.g., clear, non-bubbly ice) were observed around the clay nodules, which precludes the possible cryoconite nature of the nodules.

Conclusions
This preliminary physical, mineralogical, and geochemical characterization of dust particles from deep layers of the Guliya ice cap provides a rare glimpse of the aeolian dust composition that predates the Last Glacial Maximum.The comparison between the deep layers of the Summit and the Plateau sites and between stadial and interstadial transition layers of the Plateau core reveals the high variability of the Guliya dust record and demonstrates its value as an indicator of paleo-environmental conditions prevailing in NW Tibet since the Late Pleistocene.
The reduced size sorting and the angular shape of coarse grains on dust filters from the deep layers of the GS core indicate a mix of proximally and distally sourced aeolian material.The silt fraction is dominant, and its granite-like mineralogy also includes Ca-pyribole (amphibole), chlorite (clinochlore), and100 calcite and halite.This composition is consistent with the enrichment in soluble Sr and Na of the ice samples containing these dust layers and is indicative of weak chemical weathering and dry climatological conditions at the time the dust was entrained, transported, and deposited onto the ice cap.The higher (La/Yb) N also reflects the moderate weathering experienced by the deeper (older) Summit dust.The Summit dust REE pattern and ε Nd vs. 87 Sr/ 86 Sr signature indicate that its source is likely in the Taklimakan desert, even though only a size-resolved isotopic analysis could allow the differentiation of distal and proximal provenances.
The investigation of samples from ice sections preceding and succeeding stadial/interstadial (cold/warm) transitions shows that the dust from "warm" GP ice shares some of the characteristics of the GS dust.In particular, its enrichment in Sr and in other soluble TEs (i.e., Na, Nb, Mo, Sn, and U) and its more fractionated L/HREE are indications that weathering-sensitive minerals have been sourced and transported to the glacier surface under a drier climate.However, contrary to the GS dust, the REE pattern and the ε Nd vs. 87 Sr/ 86 Sr signature of the "warm" GP dust do not point to a distinct geographical source but rather to a mix of West Taklimakan loess and West Kunlun moraine material.This preliminary study of the dust from warm/cold transitions also reflects the larger contribution of clay-size particles to the total dust volume during cold periods.This, and the upper continental crust-like REE pattern of the dust from a "cold" GP section, suggests that the ice cap then receives a bigger proportion of well-mixed, long-range-transported dust.This mix could possibly originate from the north and the Tarim and Qaidam Basins as suggested by its isotopic signature, similar to that of the modern dust from Dunde ice cap.
This work also describes the dust from the GP ice core layers older than 110 kyr, which is present in various unusual stratigraphical forms: diffuse horizontal brown or greenish layers, discrete scattered particles, and grey nodules.The geographical source of the grey dust nodules cannot yet be confidently ascertained as its more radiogenic 87 Sr/ 86 Sr value could only be approximately matched to the fine CLP loess deposit and fine Qaidam surface sands.However, its preliminary geochemical characterization already provides important paleo-environmental information and keys for further interpretation of the GP ice core archives.The grey dust is well size-sorted and composed mainly of phyllosilicates (illitic clay, biotite, chlorite, vermiculite), Ti-oxides (anatase, rutile, or brookite), K-feldspar, and trace xenotime and monazite, which are resistant to weathering.The abundance of Fe/Ti-oxides and the inherent enrichment in HREE are additional evidence that the bottom mineral dust experienced intense chemical weathering, suggesting wet conditions at the source.At this depth in the core, the ice enrichment in Rb, as well as in Cs, Li, and F may be explained by the presence of vermiculite and its high adsorption capacity of these elements in their ionic form.It is therefore suggested that geothermal waters, highly concentrated in Cs, F, and Li, may have contributed to the leaching of the grey dust before its deposition on the Guliya Plateau more than 110 kyr ago.
These initial results are encouraging and demonstrate the importance of characterizing additional samples from these two cores to develop higher resolution Guliya dust records.Investigating the temporal variations of isotopic and elemental geotracers in the Guliya ice will afford the opportunity to gain new insight into regional paleo-synoptic conditions, including their variability and associated environmental changes.Future research on compositional changes in the dust will also provide a fundamental reference for biological investigation of microbial populations within the ice and their evolution over different time periods.the sketch mark the location of diffuse brown or greenish dust layers.Dots of various sizes symbolize dust particles and nodules.The grey dust nodules are larger within the deepest 5 m of the core.Note that the size of particles and nodules is not at scale in the sketch.The arrows point to the top depth of the 6 deepest 1 m ice sections discussed in this paper.The color difference is still noticeable among dust filters, as shown here in small photos of T292 and T308 filters.Photos of deep ice core sections show the grey nodules (T311, left) and horizontal diffuse layers (T312, right).Bubbles are present throughout the core and extend all the way to the bedrock, confirming the ice's meteoritic nature.The nodules' increasing size with depth could reflect growth by aggregation with time [60].Smaller and much less conspicuous grey nodules have been observed in shallower parts of the GP core, which are not discussed in this paper.   Sr/ 86 Sr of the five Guliya dust filters (brown squares) compared to Sr and isotopic arrays of dust sampled in surface snow and shallow ice cores from Tien Shan and Pa [43,48], Tanggula, East Kunlun, Qilian [48], Dunde [35], and Everest (Himalaya) [44].   Sr/ 86 Sr of the five Guliya dust filters (brown squares) compared to Sr and Nd isotopic arrays of dust sampled in surface snow and shallow ice cores from Tien Shan and Pamirs [43,48], Tanggula, East Kunlun, Qilian [48], Dunde [35], and Everest (Himalaya) [44].[43,48], Tanggula, East Kunlun, Qilian [48], Dunde [35], and Everest (Himalaya) [44].), Miaoergou (China [43]), Everest (Nepal [44]), Mt.St. Elias (Canada [34]), GISP2 (Greenland [26]), Dye 3 (Greenland [22]) and Talos Dome (Antarctica [100]).These ice cores cover different time scales.
Table A1.TE enrichment factor of dust filters relative to the UCC [69].Ti is the crustal element of reference and its concentration (ppm) in each sample is displayed in italics.

Figure 1 .
Figure 1.Third Pole region with the location of the Guliya ice cap (yellow star) in the Kunlun Mountains (yellow dashed line), the main deserts and loess sub-regions.The arrows depict the trajectories of the main seasonal air masses that influence NW Tibet climate.(Source background image: ESRI, world imagery.).

Figure 1 .
Figure 1.Third Pole region with the location of the Guliya ice cap (yellow star) in the Kunlun Mountains (yellow dashed line), the main deserts and loess sub-regions.The arrows depict the trajectories of the main seasonal air masses that influence NW Tibet climate.(Source background image: ESRI, world imagery).

Figure 2 .
Figure 2. 5 m-averaged δ 18 O record vs. depth of the Guliya 2015 Summit core (GS) (upper panel); δ 18 O time series of the 1992 Plateau core[11] illustrating the climatic periods recorded in the Guliya Plateau (GP) ice core (lower panel) as defined by[11].The vertical bars mark the 1 m ice core sections sampled for microparticle concentration and size, and trace element concentrations in ice corresponding to "warm" (purple/red) or "cold" (pink/blue) periods.The bars are increasingly wider with age (and depth) as 1 m of ice comprises more time toward the bottom of the glacier due to the increased thinning of successive ice layers with depth.Note that the vertical bars representing T290 and T292 (brown), and T305 and T306 (grey) are juxtaposed and therefore appear as two bars instead of four on the time axis.The crossed circles indicate the depth at which dust filtrate from 1 m ice sections were collected on filters (i.e., dust filter samples) for geochemical and isotopic analyses (see Table1).

Figure 2 .
Figure 2. 5 m-averaged δ 18 O record vs. depth of the Guliya 2015 Summit core (GS) (upper panel); δ 18 O time series of the 1992 Plateau core[11] illustrating the climatic periods recorded in the Guliya Plateau (GP) ice core (lower panel) as defined by[11].The vertical bars mark the 1 m ice core sections sampled for microparticle concentration and size, and trace element concentrations in ice corresponding to "warm" (purple/red) or "cold" (pink/blue) periods.The bars are increasingly wider with age (and depth) as 1 m of ice comprises more time toward the bottom of the glacier due to the increased thinning of successive ice layers with depth.Note that the vertical bars representing T290 and T292 (brown), and T305 and T306 (grey) are juxtaposed and therefore appear as two bars instead of four on the time axis.The crossed circles indicate the depth at which dust filtrate from 1 m ice sections were collected on filters (i.e., dust filter samples) for geochemical and isotopic analyses (see Table1).
3.1.Spatial and Temporal Variability of Dust Concentration, Size, and Mineralogy 3.1.1.Dust Concentration and Particle Size Distribution in Ice Samples

Figure 3 .
Figure 3. Dust concentration (number of particles/mL) measured by Coulter counter in 1 m ice sections from the GS and GP cores.For each section, the box represents concentrations falling within the 25-75 percentile range.The horizontal bar in the box is the median dust concentration, and the full data range is represented by the whiskers.The dust concentrations of 18 O-enriched or "warm period" ice samples are shown in purple (GS) and red (GP).The dust concentrations of 18 O-depleted or "cold period" ice samples are shown in pink (GS) and blue (GP).Dust concentration data for the deepest GP ice samples are shown in brown and grey, referring to the dust's dominant color, independent of the sample's mean δ 18 O value.The circles mark in which 1 m section a dust filter was sampled and analyzed for elemental and isotopic geochemistry (crossed circles) or mineralogy (open circles).

Figure 3 .
Figure 3. Dust concentration (number of particles/mL) measured by Coulter counter in 1 m ice sections from the GS and GP cores.For each section, the box represents concentrations falling within the 25-75 percentile range.The horizontal bar in the box is the median dust concentration, and the full data range is represented by the whiskers.The dust concentrations of 18 O-enriched or "warm period" ice samples are shown in purple (GS) and red (GP).The dust concentrations of 18 O-depleted or "cold period" ice samples are shown in pink (GS) and blue (GP).Dust concentration data for the deepest GP ice samples are shown in brown and grey, referring to the dust's dominant color, independent of the sample's mean δ 18 O value.The circles mark in which 1 m section a dust filter was sampled and analyzed for elemental and isotopic geochemistry (crossed circles) or mineralogy (open circles).

Figure 4 .
Figure 4. (a) number and (b) relative volumetric size distribution of microparticles measured by Coulter counter (14 bins) in 1 m ice sections prior to filtering.3.1.2.Mineralogy of the Deep GS Dust Filter (T50): Dry Conditions and Proximal Influence

Figure 4 .
Figure 4. (a) number and (b) relative volumetric size distribution of microparticles measured by Coulter counter (14 bins) in 1 m ice sections prior to filtering.

Figure 6 .
Figure6.X-ray diffraction (XRD) patterns for the dust sampled from the deepest layers of the GS (T50) and the GP (T309) cores discussed in this study.The T50 diffractogram is shifted up to allow the identification of key features and to facilitate the visual comparison with T309.The orange and green shadings highlight minerals found in higher abundance, respectively, in T50 and in T309 (conventional mineral abbreviations from USGS.gov).

Figure 6 .
Figure6.X-ray diffraction (XRD) patterns for the dust sampled from the deepest layers of the GS (T50) and the GP (T309) cores discussed in this study.The T50 diffractogram is shifted up to allow the identification of key features and to facilitate the visual comparison with T309.The orange and green shadings highlight minerals found in higher abundance, respectively, in T50 and in T309 (conventional mineral abbreviations from USGS.gov).

Figure 6 .
Figure 6.X-ray diffraction (XRD) patterns for the dust sampled from the deepest layers of the GS (T50) and the GP (T309) cores discussed in this study.The T50 diffractogram is shifted up to allow the identification of key features and to facilitate the visual comparison with T309.The orange and green shadings highlight minerals found in higher abundance, respectively, in T50 and in T309 (conventional mineral abbreviations from USGS.gov).

Figure 7 .
Figure 7. Enrichment factor (EF) of 18 TEs analyzed in GS and GP dust filters.EF is calculated using Ti as the TE of reference (Fe was not analyzed in the filters) and the Upper Continental Crust (UCC) [69] as the crustal reference.

Figure 7 .
Figure 7. Enrichment factor (EF) of 18 TEs analyzed in GS and GP dust filters.EF is calculated using Ti as the TE of reference (Fe was not analyzed in the filters) and the Upper Continental Crust (UCC) [69] as the crustal reference.

Figure 8 .
Figure 8. TE average (± st dev) enrichment factor (EF) in ice samples from (a) GS, (b)18 O-depleted (cold) and18 O-enriched (warm) GP sections, (c) GP ice sections containing "brown" or "grey" dust.EF is calculated using Fe and the local PSAs[55] as the crustal references.Here, Fe is the preferred crustal element of reference to illustrate the Ti enrichment variability among the various types of ice samples.

Figure 9 .
Figure 9. Chondrite-normalized REE patterns of dust filters from GS and GP ice cores compared to Taklimakan fine desert sand, West Taklimakan loess [36] and to moraines of west Kunlun

Figure 9 .
Figure 9. Chondrite-normalized REE patterns of dust filters from GS and GP ice cores compared to (a) Taklimakan fine desert sand, (b) West Taklimakan loess[36] and to (c) moraines of west Kunlun Mountains[92] represented by the shaded area.REE abundances in the C1 chondrite[93] are used for normalization.

Figure A1 .
Figure A1.Simplified visual stratigraphy of the deepest 45 m of the GP ice core.Horizontal bars in the sketch mark the location of diffuse brown or greenish dust layers.Dots of various sizes symbolize dust particles and nodules.The grey dust nodules are larger within the deepest 5 m of the core.Note that the size of particles and nodules is not at scale in the sketch.The arrows point to the top depth of the 6 deepest 1 m ice sections discussed in this paper.The color difference is still noticeable among dust filters, as shown here in small photos of T292 and T308 filters.Photos of deep ice core sections show the grey nodules (T311, left) and horizontal diffuse layers (T312, right).Bubbles are present throughout the core and extend all the way to the bedrock, confirming the ice's meteoritic nature.The nodules' increasing size with depth could reflect growth by aggregation with time [60].Smaller and much less conspicuous grey nodules have been observed in shallower parts of the GP core, which are not discussed in this paper.

Figure A2 .
Figure A2.Sampling scheme of the ice and filter samples from a 1 m long ice core section of GP or GS.An ice sample is the entire 1 m section comprising subsamples which have been individually analyzed for trace elements (TE), major ions, dust size and concentration and oxygen isotopes.These parameters are not discussed at the subsample resolution in this manuscript.A filter sample corresponds to the insoluble particles (ø > 0.22 μm) collected by filtering all melted subsamples together on one filter.

Figure A2 . 29 Figure A3 .
Figure A2.Sampling scheme of the ice and filter samples from a 1 m long ice core section of GP or GS.An ice sample is the entire 1 m section comprising subsamples which have been individually analyzed for trace elements (TE), major ions, dust size and concentration and oxygen isotopes.These parameters are not discussed at the subsample resolution in this manuscript.A filter sample corresponds to the insoluble particles (ø > 0.22 µm) collected by filtering all melted subsamples together on one filter.Geosciences 2022, 12, x FOR PEER REVIEW 23 of 29

Figure A3 .
Figure A3.Spectral plot and Backscattered Electron (BSE) image of a HREE-bearing xenotime identified on dust filter T309 by EDXS trace mineralogy analysis.

Figure A3 .
Figure A3.Spectral plot and Backscattered Electron (BSE) image of a HREE-bearing xenotime i tified on dust filter T309 by EDXS trace mineralogy analysis.

Figure A4 .
Figure A4.Probability of occurrence of excessive fluoride concentrations in Asian groundwater panel) [85] and median fluoride concentration (ppb) measured in 1 m ice samples from the G ice cores discussed in this paper (right panel).

Figure A4 .
Figure A4.Probability of occurrence of excessive fluoride concentrations in Asian groundwater (left panel) [85] and median fluoride concentration (ppb) measured in 1 m ice samples from the Guliya ice cores discussed in this paper (right panel).

Figure A3 .
Figure A3.Spectral plot and Backscattered Electron (BSE) image of a HREE-bearing xenotime identified on dust filter T309 by EDXS trace mineralogy analysis.

Figure A4 .
Figure A4.Probability of occurrence of excessive fluoride concentrations in Asian groundwater (left panel) [85] and median fluoride concentration (ppb) measured in 1 m ice samples from the Guliya ice cores discussed in this paper (right panel).

Table 1 .
Sampling and analytical information.

Table 2 .
Sr-Nd isotopic composition of dust filter samples.

Table 2 .
Sr-Nd isotopic composition of dust filter samples.

Table A1 .
[69]nrichment factor of dust filters relative to the UCC[69].Ti is the crustal element of reference and its concentration (ppm) in each sample is displayed in italics.