# Greenland Ice Core Record of Last Glacial Dust Sources and Atmospheric Circulation

G. Újvári, U. Klötzli, T. Stevens, A. Svensson, P. Ludwig, T. Vennemann, S. Gier, M. Horschinegg, L. Palcsu, D. Hippler, et al.

► To cite this version:

G. Újvári, U. Klötzli, T. Stevens, A. Svensson, P. Ludwig, et al., Greenland Ice Core Record of Last Glacial Dust Sources and Atmospheric Circulation. Journal of Geophysical Research: Atmospheres, 2022, 127, 10.1029/2022JD036597. insu-03839369

HAL Id: insu-03839369

https://insu.hal.science/insu-03839369

Submitted on 4 Nov 2022

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

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

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

Distributed under a Creative Commons Attribution 4.0 International License
# JGR Atmospheres

## RESEARCH ARTICLE

10.1029/2022JD036597

**Key Points:**

* Hf-Nd isotope data demonstrate that Central Greenland last glacial dust is not sourced from mixed SE Asian and circum-Pacific volcanic dust

* Most likely sources are the Taklimakan in SE Asia, European and/or North African glacial dust deposits based on Hf-Sr-Nd and D/H isotopes

* Model simulations show that sourcing of North Greenland Ice Core Project ice core dust from Europe is plausible, the addition of North African dust cannot be excluded

**Supporting Information:**

Supporting Information may be found in the online version of this article.

**Correspondence to:**

G. Újvári,
ujvari.gabor@csfk.org

**Citation:**

Újvári, G., Klötzli, U., Stevens, T., Svensson, A., Ludwig, P., Vennemann, T., et al. (2022). Greenland ice core record of last glacial dust sources and atmospheric circulation. *Journal of Geophysical Research: Atmospheres*, 127, e2022JD036597. https://doi.org/10.1029/2022JD036597

Received 3 FEB 2022
Accepted 8 JUL 2022

**Author Contributions:**

**Conceptualization:** G. Újvári

**Data curation:** G. Újvári

**Formal analysis:** G. Újvári, T. Stevens, A. Svensson, P. Ludwig, T. Vennemann, S. Gier, M. Horschinegg, D. Hippler, J. Kovács, C. Di Biagio

**Funding acquisition:** G. Újvári, U. Klötzli

**Investigation:** G. Újvári, U. Klötzli, T. Stevens, D. Hippler

**Methodology:** G. Újvári, U. Klötzli, A. Svensson, P. Ludwig, T. Vennemann,

## Greenland Ice Core Record of Last Glacial Dust Sources and Atmospheric Circulation

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

G. Újvári¹,²,³ D, U. Klötzli³, T. Stevens⁴ D, A. Svensson⁵ D, P. Ludwig⁶ D, T. Vennemann⁷ D, S. Gier⁸ D, M. Horschinegg³, L. Palcsu⁹, D. Hippler¹⁰ D, J. Kovács¹¹,¹² D, C. Di Biagio¹³, and P. Formenti¹³ D

¹Centre for Astronomy and Earth Sciences, Institute for Geological and Geochemical Research, Eötvös Loránd Research Network, Budapest, Hungary, ²CSFK, MTA Centre of Excellence, Budapest, Hungary, ³Department of Lithospheric Research, University of Vienna, Vienna, Austria, ⁴Department of Earth Sciences, Uppsala University, Uppsala, Sweden, ⁵Physics of Ice, Climate and Earth, Niels Bohr Institute, University of Copenhagen, Copenhagen, Denmark, ⁶Institute for Meteorology and Climate Research, Karlsruhe Institute of Technology, Karlsruhe, Germany, ⁷Institute of Earth Surface Dynamics, University of Lausanne, Lausanne, Switzerland, ⁸Department of Geology, University of Vienna, Vienna, Austria, ⁹Isotope Climatology and Environmental Research Centre, Institute for Nuclear Research, Debrecen, Hungary, ¹⁰Institute of Applied Geosciences, Graz University of Technology, Graz, Austria, ¹¹Environmental Analytical and Geoanalytical Research Group, Szentágothai Research Centre, University of Pécs, Pécs, Hungary, ¹²Institute of Geography and Earth Sciences, University of Pécs, Pécs, Hungary, ¹³Université de Paris Cité and University Paris Est Creteil, CNRS, LISA, Paris, France

**Abstract** Abrupt and large-scale climate changes have occurred repeatedly and within decades during the last glaciation. These events, where dramatic warming occurs over decades, are well represented in both Greenland ice core mineral dust and temperature records, suggesting a causal link. However, the feedbacks between atmospheric dust and climate change during these Dansgaard-Oeschger events are poorly known and the processes driving changes in atmospheric dust emission and transport remain elusive. Constraining dust provenance is key to resolving these gaps. Here, we present a multi-technique analysis of Greenland dust provenance using novel and established, source diagnostic isotopic tracers as well as results from a regional climate model including dust cycle simulations. We show that the existing dominant model for the provenance of Greenland dust as sourced from combined East Asian dust and Pacific volcanics is not supported. Rather, our clay mineralogical and Hf-Sr-Nd and D/H isotopic analyses from last glacial Greenland dust and an extensive range of Northern Hemisphere potential dust sources reveal three most likely scenarios (in order of probability): direct dust sourcing from the Taklimakan Desert in western China, direct sourcing from European glacial sources, or a mix of dust originating from Europe and North Africa. Furthermore, our regional climate modeling demonstrates the plausibility of European or mixed European/North African sources for the first time. We suggest that the origin of dust to Greenland is potentially more complex than previously recognized, demonstrating more uncertainty in our understanding dust climate feedbacks during abrupt events than previously understood.

**Plain Language Summary** Abrupt climate change represents an existential threat to civilization. However, the feedbacks that modulate these abrupt changes are poorly understood, undermining our ability to predict future events. Last glacial Greenland ice core records show abrupt climate events coupled to changes in abundance of atmospheric mineral dust, but how dust impacts these events is unclear as the processes involved in dust emission remain elusive. Here we apply multiple novel tracers of dust provenance as well as regional dust cycle modeling to address this uncertainty. We show that the dominant model of mixed East Asian and Pacific volcanic sources to Greenland dust is not supported. Instead, multiple other source scenarios are plausible, demonstrating far more uncertainty in dust climate feedbacks than previously understood.

## 1. Introduction

Abrupt climate changes have occurred repeatedly and within decades in the recent geological past (Alley et al., 2003), including the last glaciation. These changes were centered on the North Atlantic, where vigorous wind systems encountered the southernmost extension of sea ice and oceanic currents with connections to the deep ocean via convection (Li & Born, 2019; Lynch-Stieglitz et al., 2007). There is growing concern that recent CO₂ forcing and warming will lead to a global cascade of rapid regime shifts in the near future (Lenton

© 2022. The Authors.

This is an open access article under the terms of the Creative Commons Attribution License, which permits use, distribution and reproduction in any medium, provided the original work is properly cited.
S. Gier, M. Horschinegg, L. Palcsu, J.

Kovács, C. Di Biagio, P. Formenti

**Project Administration:** G. Újvári, U.
Klötzli

**Resources:** G. Újvári

**Software:** G. Újvári, P. Ludwig

**Supervision:** U. Klötzli, A. Svensson, L.

Palcsu, D. Hippler, P. Formenti

**Validation:** G. Újvári, U. Klötzli, A.
Svensson, P. Ludwig, T. Vennemann,

S. Gier, M. Horschinegg, L. Palcsu, D.

Hippler, C. Di Biagio, P. Formenti

**Visualization:** G. Újvári, P. Ludwig

**Writing – original draft:** G. Újvári, T.
Stevens

et al., 2019), resulting in abrupt climate changes, which are difficult to model and predict, but would be a global threat to civilization. The reason for this lack of ability to predict tipping points is due to the limited understanding of feedbacks in the climate system during abrupt events. In particular, the role of atmospheric dust, which has both direct and indirect influence on regional and global climate through diverse physical and biogeochemical processes and feedbacks (Mahowald, 2011; Sokolik & Toon, 1996; Tegen et al., 1996), is still poorly known.

Abrupt changes in climate of the Northern Hemisphere (NH) over the last glaciation are well recorded in Greenland ice. High-resolution ice core oxygen isotope records reveal that the cold glacial climate was interrupted by numerous rapid shifts to warmer interstadial conditions, called Dansgaard-Oeschger (D-O) events, over the last glaciation (Bond et al., 1993; Dansgaard et al., 1993). Air temperatures rose by 5–16°C within decades during D-O events in Greenland (Huber et al., 2006; Kindler et al., 2014), followed by a much slower drop in temperature back to stadial conditions. This D-O type climate variability is coupled with dust concentration and particle size variations, with rapid decreases in Ca²⁺ ions (a dust concentration proxy) at the onset of interstadials and slower increases following their end (Fuhrer et al., 1999). Dust particle concentrations in the Greenland Ice Core Project (GRIP) and North Greenland Ice Core Project (NGRIP, Figure 1) ice cores were found to be higher by a factor of 100 during the Last Glacial Maximum (LGM: 26–19 ka, Clark et al., 2009) than the Preboreal warm period and typically a factor of 8 higher at the sharp transitions of D-O events, with modes of particle volume size distributions being systematically coarser over cold periods (Ruth et al., 2003; Steffensen, 1997). These abrupt variations in dust characteristics were variously accounted for by changes in transit and atmospheric residence time, transport distance, extent and intensity of the polar atmospheric cell, dust source strength, aridity and storminess (Fischer et al., 2007; Hansson, 1994; Mayewski et al., 1994; Ruth et al., 2003; Schüpbach et al., 2018; Steffensen, 1997), nonetheless the precise causal mechanisms of these changes remain elusive. Comparison of δ¹⁸Oice and dust (mass, Ca²⁺) records of the NGRIP ice core demonstrates that the average lag of dust versus δ¹⁸Oice is 1 ± 8 years for all D-O inceptions (Ruth et al., 2007). This synchrony of dust concentration and δ¹⁸Oice changes has recently been confirmed by higher resolution analyses of NGRIP and NEEM ice core data (Capron et al., 2021; Erhardt et al., 2019), indicating that the climate of continental dust sources and Greenland must have been coupled. However, the mechanisms of this coupled response over wide areas of the NH remain poorly constrained. Resolving the source(s) of Greenland ice core dust is a key step in resolving the uncertainty over this coupling, as it allows the environmental controls on dust emission to be constrained, as well as provides insight into major dust transport pathways and atmospheric circulation patterns during abrupt climate events of the Last Glacial Period. This facilitates a better understanding of feedbacks between dust and climate at decadal to millennial timescales and the role of different atmospheric patterns in dispersing dust over the NH.

Early studies of the clay mineralogy and Sr-Nd-Pb isotopic compositions of dust aerosols in the Greenland Ice Sheet Project 2 (GISP2) and GRIP ice cores suggested derivation from the SE Asian deserts (e.g., Gobi, Taklimakan) and Chinese Loess Plateau (CLP) over the LGM, with additions of circum-Pacific volcanic material in proportions of 10%-25% to account for the Sr isotope compositions of ice core dust (Biscaye et al., 1997; Svensson et al., 2000). Later work partly supported this SE Asian origin (Burton et al., 2007; Han et al., 2018; Meyer et al., 2017), and also for preindustrial aerosols in the Dye-3 ice core (Lupker et al., 2010) or modern dust collected in snow pits at the NGRIP site (Bory, Biscaye, & Grousset, 2003). However, despite the wide acceptance of this model, the sources of dust have not as yet been unequivocally constrained and other potential dust sources to Greenland, including Africa (Han et al., 2018; Meyer et al., 2017) and central Europe (Újvári et al., 2015), and admixtures of dust from sources other than SE Asia (Lupker et al., 2010; Meyer et al., 2017), were also proposed based on expanded datasets and currently are equally supported by available data. Studies of preindustrial/Holocene aerosols recovered from ice cores at lower elevation sites, located closer to the edge of the ice cap (Hans Tausen, Renland, RECAP) demonstrate dust derivation from proximal sources (Bory, Biscaye, Piotrowski, & Steffensen, 2003), including central East Greenland (Simonsen et al., 2019).

Model simulations of the LGM dust cycle are likewise diverse in simulated potential dust sources to central Greenland. SE Asia is one of the major LGM dust emission centers and aerosol sources for central Greenland ice in numerous models (De Angelis et al., 1997; Lunt & Valdes, 2001; Mahowald et al., 1999, 2006, 2011; Reader et al., 1999; Werner et al., 2002), while glaciogenic/non-glaciogenic sources in the continental USA, Alaska and Siberia are identified as additional potential contributors (Lambert et al., 2015; Mahowald et al., 2006). An active dust cycle in Europe can also be identified in regional simulations (Ludwig et al., 2021; Schaffernicht et al., 2020). Over the LGM, the Sahara Desert was the most active and largest dust source (Albani et al., 2018;
<figure><img src="image_3.png" /></figure>

**Figure 1.** Locations of the North Greenland Ice Core Project (NGRIP) and other central Greenland ice cores and potential source area samples across the Northern Hemisphere. Sites displayed on map are as follows (codes also listed in Dataset S1 and S2). Greenland: NGRIP – North Greenland Ice Core Project ice core, GRIP - Greenland Ice Core Project ice core, GISP2 – Greenland Ice Sheet Project 2 ice core, NEEM - North Greenland Eemian Ice Drilling. North America: SiC – Silver Creek (Yukon, Canada), Obe – Obert (Nebraska, USA), Pra – Prairie Lake (Nebraska, USA), Wan – Wauneta (Nebraska, USA). European sites are abbreviated, where necessary (with original codes given in parenthesis, and in Dataset S2): B – Basaharc (Bh, Hungary), Bk – Bialy Kosciol (Poland), C – Crvenka (Crv, Serbia), D – Dunaszekcső (Dsz, Hungary), KW – Krems-Wachtberg (KW-GH9, Austria), LMJ – La Motte (Jersey, Channel Islands), M – Mende (Me, Hungary), Nus – Nussloch (Germany), P – Paks (Pa, Hungary), PbE – Pegwell Bay (UK), RG – Raigorod (Lower Volga, Russia), T – Titel core (Tic-1A, Serbia), Url – Urluia (Romania), Z – Zmajevac (Zm, Croatia). North Africa: Mao – Maouna (Tunisia), MoLC – off-road site (Erg Chebbi, Morocco), S107 – off-road site (Ahaggar Mts., Algeria). Central Asia: Cha – Chasmanigar (Tajikistan), Rem – Remizovka (Kazakhstan). Northeast Asia: KHG – Khangar volcano (Kamchatka, Russia), KO – Kurile lake (Kamchatka, Russia), Sib – Tumara valley (Siberia, Russia). Southeast Asia: BD – road site (Badain Jaran desert, China) Bei – Beiguoyuan (Chinese Loess Plateau (CLP), China), Fan – Fanshan (around Beijing, China), Jb54 – dune site (Junggar basin, China), Lin – Lingtai (CLP, China), Luo – Luochuan (CLP, China), OD – road site (Otin Dag desert, China), TD/TD2 – road sites (Tengger desert, China), Td06/Td25/Td34 – overbank, stream (fan) and loess sites (Taklimakan desert, China). Regional color coding corresponds with those in Figures 2–4.

Takemura et al., 2009; Werner et al., 2002) and the shift of the Hadley cell circulation outer boundary toward the equator may have enhanced desert source areas by factors of 2–3 during the last glacial, allowing meridional dust transport by mid-latitude storm systems to polar regions (Chylek et al., 2001).

As such, despite the prevalence of the SE Asian desert/loess model for dust sources to Greenland, available dust provenance data and computer simulations provide a range of plausible scenarios. Such ambiguities mean we have limited knowledge of wider glacial dust loading and present a major impediment to understanding what specific areas were emitting dust during abrupt climate change events, the cause of these emissions, and the consequences for global climate. To gain deeper insights into this inherently complex problem, we present new Sr–Nd–Hf and D/H isotope compositional data from early LGM (25.64–24.85 ka) dust samples of the NGRIP ice core and new clay mineralogy and isotopic data from the fine (<5 and <2 µm) fractions of a large range of potential source area (PSA) samples collected from major last glacial dust hotspots over the NH (Figure 1), combined with new
regional model simulations of LGM dust transport over Europe. Since Sr isotope compositions strongly suffer from grain size and pre-treatment effects and Nd isotopes are not source diagnostic for Greenland ice core dust (Újvári et al., 2015), two additional, novel dust tracers are applied here; Hf and H isotopes. To our knowledge, hafnium isotope ratios have only once been utilized to trace sources of mineral dust from preindustrial segments in the Dye-3 ice core (Lupker et al., 2010), potentially due to analytical difficulties, yet are likely highly source diagnostic (Újvári et al., 2015, 2018). A recently established methodology developed for the combined Hf-Sr-Nd isotope measurements of low mass dust samples (2–10 mg; Újvári et al., 2021) has now made it possible to analyze the Hf isotope compositions of a relatively high number of PSA samples, together with ice core dust. In addition, the D/H isotope composition of clay structural water (δDcsw, i.e., that of hydroxyl groups of clays; Savin & Epstein, 1970) of size separated dust samples, has been analyzed for the first time in ice core dust provenance research. In general, the D/H isotope composition of naturally occurring clay minerals (kaolinite, montmorillonite, etc.) are used to reconstruct conditions under which the minerals are formed, and thus to better understand paleoclimate and pedogenetic processes (Savin & Hsieh, 1998; Sheppard & Gilg, 1996). Here, δDcsw, which is a combined signal of a mixture of clay minerals present in the dust separates, is simply used as a source fingerprint and not for making inferences about the paleoclimate. This, we will show, is feasible given that the PSAs form distinct groups in terms of δDcsw, allowing further constraints on the origin of ice core dust. Finally, we present regional model simulations results of the LGM dust cycle for Europe (including the northernmost part of Africa) with the WRF-Chem model to investigate if dust emitted from these areas may have contributed to dust deposition on central parts of the Greenland Ice Sheet, as this potential pathway has not been simulated up until today.

## 2. Materials and Methods

### 2.1. Studied Samples and Pretreatments

Four ice core samples were selected for this study from the NGRIP ice core, Greenland (Figure 1), corresponding to one of the periods with the highest dust accumulation of the LGM. When selecting the samples, we focused on the GS-3 stadial period (Figure S1 in Supporting Information S1), which is characterized by the highest dust concentrations (4–8 mg dust/kg ice; Ruth et al., 2003) and did not investigate interstadials with much lower concentrations (0.3–0.7 mg dust/kg ice) as at least 2–3 mg dust was required for the combined Sr–Nd–Hf and H isotope analyses. The sample bags 3,327, 3,326, 3,323 and 3,306–3,307 span depths of 1829.85–1818.30 m and cover the age intervals of 25,640–25604, 25,604–25567, 25,491–25453 and 24,924–24859 years GICC05 b2k (Text S1, Figure S1 in Supporting Information S1). Impurities were obtained by melting about 0.6 kg ice each (except for the combined bags 3,306–3,307: about 1 kg) at room temperature, after removing the outer ca. 3–5 mm of ice from all sides with pre-cleaned PFA chisel. Melted ice samples were collected in acid-cleaned Savillex PFA beakers and evaporated to incipient dryness at 50°C in clean laboratory environment in Vienna. Subsequently, the impurities were treated with weak (0.5 mol/L) acetic acid for 1 hr to remove carbonates and sea salt aerosols of marine origin, following Svensson et al. (2000). Both the leachates and remaining aluminosilicate fractions were analyzed for isotopic compositions.

PSA samples include modern soil, desert/dune sand, overbank, stream/fan and lake deposits and late Quaternary loess deposits, listed in Datasets S1 and S2 (also in Újvári, 2021). Loess sediments were collected from luminescence- or radiocarbon-dated sequences and most of them represent LGM or last glacial (L1) dust of the respective regions. Two pumice fall deposits of Holocene age were also investigated from Kamchatka, one of them (KHG) identified in Greenland ice cores. Long-range atmospheric transport fractionates dust and only the finest fractions (typically <10 or <5 µm) can reach Greenland from distal continental sources, as evidenced by the particle size distribution of the NGRIP ice core (Ruth et al., 2003), which is dominated by <10 µm particles (Figure S2b in Supporting Information S1). To minimize grain size and mineralogy effects on the isotopic compositions and ensure that the same particle size fractions of PSA and ice core samples are compared, the PSA samples were size separated using 5 µm hydrophobic Mitex filter and/or wet sedimentation (Stokes law) (Text S2, Figure S2 in Supporting Information S1). The <5 and <2 µm size separates were subsequently leached in 0.5 mol/L acetic acid for 1 hr and washed in 30% H₂O₂ to remove secondary carbonates and organic matter. The residues (aluminosilicate fractions) were washed with Milli-Q water and dried (50°C) for subsequent XRD or isotopic analyses.
## 2.2. Powder X-Ray Diffraction

For clay mineral analyses samples were disaggregated with diluted H₂O₂ to remove organic matter and the slurry was subsequently prepared for XRD analyses as described in detail in Moore and Reynolds (1997). The prepared clay samples were analyzed with a Panalytical PW 3040/60X'Pert PRO diffractometer (CuKα radiation, 40 kV, 40 mA, step size 0.0167, 5 s per step). The X-ray patterns were interpreted following Moore and Reynolds (1997). Smectite is identified by a broad peak at about 14 Å with Mg-saturation, which shifts to about 12 Å with K-saturation and collapses to 10 Å after heating to 550°C. Saturation of the K-saturated sample with ethylene glycol expands smectite again to 17 Å, and the Mg-saturated sample with glycerol expands to 18 Å. Chlorite was identified by peaks at 14, 7, 4.7 and 3.53 Å, which retain their position during treatments. Illite peaks at 10, 5 and 3.3 Å likewise keep their positions. Kaolinite peaks at 7 and 3.57 Å no longer appear after heating to 550°C. Mg-saturated vermiculite has a strong peak at 14 Å, and with K saturation it shifts to 10 Å. Ethylene glycol saturation of the K-saturated sample or glycerol-saturation of the Mg-saturated sample do not lead to changes in positions. The clay mineral compositions of samples were quantified after Biscaye (1965), as also described in Svensson et al. (2000). The peak areas of clay minerals in the Mg- and glycerol-saturated X-ray patterns were determined using the Panalytical X'Pert Highscore plus software.

## 2.3. Hafnium, Neodymium and Strontium Isotope Analyses

All chemical separations in this study were performed in PicoTrace class 100 clean rooms at the Department of Lithospheric Research, University of Vienna, Austria. Details on chemicals, lab wares and columns used in this study, as well as the ion-exchange chemistry of the elemental separations and purifications can be found in Újvári et al. (2021). Pretreated NGRIP dust samples and size separates of PSA samples were digested in a convection oven at 230°C using ammonium bifluoride (ABF, NH₄HF₂), and subsequently by concentrated HNO₃ and HCl (Újvári et al., 2021). Total chemistry blanks were determined to be 0.5 ng for Sr, 0.2 ng for Nd and <25 pg for Hf, all negligible in the present application. Performance of the ammonium-bifluoride method together with the column chemistry setup was thoroughly tested and demonstrated using five USGS geological reference materials, including AGV-2, BCR-2, GSP-2, RGM-2 and STM-2 (Újvári et al., 2021).

Mass spectrometry of Sr and Nd was performed at the Department of Lithospheric Research, University of Vienna, Austria using a Thermo-Finnigan Triton TI multi-collector TIMS instrument in static mode. Pure element fractions were analyzed using a Re double filament assembly. A mean ⁸⁷Sr/⁸⁶Sr ratio of 0.710259 ± 0.000002 (2SE, *n* = 20) was measured for NBS987 (ref. value: ⁸⁷Sr/⁸⁶Sr = 0.710248; Faure, 2001) and a mean ¹⁴³Nd/¹⁴⁴Nd ratio of 0.511843 ± 0.000002 (2SE, *n* = 31) for the La Jolla (ref. value: ¹⁴³Nd/¹⁴⁴Nd = 0.511858; Lugmair & Carlson, 1978) isotope standards during the analysis periods. Mass fractionation during measurement was corrected for ⁸⁶Sr/⁸⁸Sr = 0.1194, and ¹⁴⁶Nd/¹⁴⁴Nd = 0.721903, respectively. The measured ⁸⁷Sr/⁸⁶Sr and ¹⁴³Nd/¹⁴⁴Nd isotope ratios of NGRIP dust and PSA samples were finally normalized to the reference isotope ratios of 0.710248 and 0.511858 of the NBS987 and La Jolla isotope standards. Uncertainties of isotopic ratios represent 2 standard errors of the mean (2SE).

Hafnium isotope analyses were done on a Thermo Neptune Plus MC-ICP-MS equipped with an Aridus three desolvation nebulizer (flow rates of 100 µl/min) in the laboratory of the Institute of Nuclear Research, Debrecen, Hungary. A Faraday gain calibration was made prior to instrument tuning. Masses of ¹⁷²Yb, ¹⁷⁴Hf and ¹⁷⁵Lu were measured using 10¹³ Ω, while others (¹⁷⁶Hf, ¹⁷⁷Hf, ¹⁷⁸Hf, ¹⁷⁹Hf, and ¹⁸⁰Hf) using 10¹¹ Ω amplifiers. Gain calibration protocol details and Faraday cup configurations are given in Újvári et al. (2021). An exponential mass bias correction for Hf (βHf) is applied using ¹⁷⁹Hf/¹⁷⁷Hf = 0.7325 (Blichert-Toft et al., 1997; Chu et al., 2002; Griffin et al., 2006; Patchett, 1983). ¹⁷⁶Lu and ¹⁷⁶Yb interference corrections on ¹⁷⁶Hf were made using recommended values of ¹⁷⁶Lu/¹⁷⁵Lu = 0.026549 and ¹⁷⁶Yb/¹⁷²Yb = 0.58862 (Chu et al., 2002), with mass bias factors βYb = βLu = βHf. Interference corrected Hf isotope ratios and reported errors represent the mean and 2 standard error of the mean (2SE) values. Repeat analyses of the JMC-475 Hf isotope standard (10 ppb solutions) yielded mean ¹⁷⁶Hf/¹⁷⁷Hf isotopic ratios of 0.282139 ± 0.000035 (2SD, *n* = 29), 0.282141 ± 0.000019 (2SD, *n* = 4) and 0.282152 ± 0.000011 (2SD, *n* = 2) for the three analytical sessions in 2020 versus the accepted value of 0.282163 ± 0.000009 (Blichert-Toft et al., 1997). The measured ¹⁷⁶Hf/¹⁷⁷Hf isotopic ratios were normalized to the accepted ¹⁷⁶Hf/¹⁷⁷Hf isotopic ratio of JMC-475.
Nd and Hf isotopic ratios are also reported as $\epsilonNd(0) = ((^{143}Nd/^{144}Nd_{sample})^{143}Nd/^{144}Nd_{CHUR}) - 1) \times 10^4$ and $\epsilonHf(0) = ((^{176}Hf/^{177}Hf_{sample})^{176}Hf/^{177}Hf_{CHUR}) - 1) \times 10^4$ in this study using the present-day chondritic uniform reservoir (CHUR) values of 0.512630 ± 0.000011 and 0.282785 ± 0.000011 (Bouvier et al., 2008).

Uncertainties of $\epsilonNd(0)$ and $\epsilonHf(0)$ values were propagated as $\sqrt{((\partial\epsilon/\partial x)\sigma_x)^2 + ((\partial\epsilon/\partial y)\sigma_y)^2} = \sqrt{((1/y)10000\sigma_x)^2 + ((-x/y^2)10000\sigma_y)^2}$, where $\epsilon = \epsilonNd(0)$ or $\epsilonHf(0)$, $x = ^{143}Nd/^{144}Nd_{sample}$ or $^{176}Hf/^{177}Hf_{sample}$, $y = ^{143}Nd/^{144}Nd_{CHUR}$ or $^{176}Hf/^{177}Hf_{CHUR}$, $\sigma_x$ and $\sigma_y$ are uncertainties of x and y (Újvári et al., 2018). Further information on reproducibility of Sr–Nd–Hf isotope measurements can be found in Text S3 and Figures S3-S4 of Supporting Information S1.

## 2.4. Hydrogen Isotope Analysis of Clay Structural Water

Measurements of the water content and D/H ($^{2}$H/$^{1}$H) isotope compositions of bulk NGRIP dust, and fine fractions (<5 and <2 µm) of PSA samples were made using high-temperature (1,450°C) reduction methods with He-carrier gas and TC-EA attached to a gas chromatographic column and Thermo-Finnigan Delta Plus XL mass spectrometer (Bauer & Vennemann, 2014; Sharp et al., 2001). Sample mass ranged from 0.3 to 2.3 mg. To remove adsorbed surface and interlayer water, fine-grained clay separates were dehydrated in a vacuum line at 250°C for 3 hr and subsequently isolated in small break seal tubes up until analysis with a zero-blank sample tray on the TC-EA, allowing for the isotopic analysis of structurally bound water in clays only (Bauer & Vennemann, 2014). Testing with clay standards demonstrated that the TC/EA method provides accurate D/H isotopic compositions and has comparable precision and reproducibility to the Zn-based method (Bauer & Vennemann, 2014). The results are given in the standard $\delta$-notation ($\delta$D), expressed relative to V-SMOW in permil ($\permil$). The precision of the in-house kaolinite and G1 biotite standards for hydrogen isotope analyses is better than $\pm$2 $\permil$. All values were normalized using values of -125 $\permil$ for the K-17 kaolinite standard and -65$\permil$ for NBS-30. Additional information on D/H isotope measurements can be found in Text S4, Figures S5 and S6 of Supporting Information S1.

## 2.5. Dust Emission and Transport Simulations

Simulations of the LGM dust cycle over Europe at high spatial (50 km horizontal grid spacing) and temporal resolution (6 hourly output interval) were performed using the Weather and Research Forecast model with chemistry (WRF-Chem; version 4.1.2, Skamarock et al., 2019). In WRF-Chem, the GOCART (Chin et al., 2000; Ginoux et al., 2001) aerosol physics have been implemented by including algorithms for dust emissions, transport and deposition. WRF-Chem has been used in dust-only mode. Size resolved dust emissions (five dust bins: 0–2, 2–3.6, 3.6–6, 6–12 and 12–20 µm) were simulated based on the University of Cologne dust emission scheme (Shao, 2004). The dust flux into the atmosphere is triggered by saltation bombardment and aggregate disintegration, depending on the soil texture. Dust emission is then proportional to the friction velocity and the simulations are run in dust-only mode, implying that only emission, transport, and deposition is considered, without any radiative and/or cloud feedbacks. Additionally, we used the polar WRF extension (version 4.1.1, Hines et al., 2015) to account for an optimized surface energy balance and heat transfer for the Noah land surface model over sea ice and permanent ice surfaces and to update the sea ice concentration during the simulation. The potential dust source areas were determined following the approach by Ginoux et al. (2001). Other physical parametrizations applied in this study are summarized in Table S1 of Supporting Information S1 (Text S5 in Supporting Information S1). The WRF-Chem model has been successfully used in this configuration to simulate different aspects of the regional LGM climate and dust cycle over Europe (Ludwig et al., 2017, 2021; Schaffernicht et al., 2020).

The WRF model is forced by initial and 6-hourly lateral boundary conditions, including time-varying sea-ice and SST, from global LGM simulation (1.875° × 1.875° horizontal grid spacing in longitude/latitude) by the Max-Planck-Institute of Meteorology earth system model (MPI-ESM-P, Jungclaus et al., 2012), performed during the third phase of the Paleomodel Intercomparison Project (PMIP3, Braconnot et al., 2012). An overview of the capabilities of the MPI to simulate the general LGM climate conditions and atmospheric circulation in comparison with other global climate models is given in for example, Ludwig et al. (2016) and Kageyama et al. (2021). The MPI-ESM simulations are very close to the mean of the PMIP3 simulations in terms of temperature and precipitation (Figures 1e and 7e in Kageyama et al., 2021), so the MPI-ESM-P cannot be considered as an outlier either at the high or low end of the model spread. In comparison with proxy data, the MPI-ESM-P model is in
good agreement with the proxy reconstructions for temperature and outperforms some other PMIP3 models, which show larger deviations from proxy data (Figures 10b, 10d and 11 in Kageyama et al., 2021).

The dust in the WRF simulations is only sourced within the regional domain since the MPI-ESM-P does not provide any dust data in its standard PMIP3 setup and thus no tracers can be passed at the lateral boundaries. This makes the results much easier to interpret, as the dust that arrives at Greenland can only be attributed to European/ North African sources.

Several modifications have been introduced to the WRF model to account for the LGM boundary conditions. These included the implementation of an orbital routine accounting for changes in orbital parameters (Prömmel et al., 2013) and modifications of the vegetation cover and land use (CLIMAP Project Members, 1984; for implementation see Ludwig et al., 2017). Further corrections were related to ice sheets, land-sea-mask and trace gas concentration changes based on PMIP3 guidelines for the 21k experiment for consistency with the MPI-ESM-P forcing data (https://pmip3.lsce.ipsl.fr/) and the alpine ice shield (Seguinot et al., 2018). An overview of specific model settings is given in Table S2 of Supporting Information S1. The added value of regional paleoclimate model simulations from the WRF model in comparison with the outcome of the global MPI-ESM-P model have been described for example, in Ludwig et al. (2017, 2019). Furthermore, Pinto and Ludwig (2020) analyzed the statistics and characteristics of cyclones based on MPI-ESM-P driven WRF simulations. The WRF LGM dust simulation for this study comprises 32 years, where the first 2 years have been excluded from the analysis to account for a sufficient model spin-up time. The horizontal grid spacing is 50 km over the model domain, which included the northern Atlantic Ocean and Europe, with the British and Fennoscandinavian ice sheets (Figure 5). The vertical grid consists of 35 levels up to 50 hPa (~20 km height) to cover the whole troposphere. The vertical levels are irregularly distributed with height, with the highest vertical resolution in the boundary layer. The time step is set to 180 s for the simulation.

The 50 most intense dust deposition events (TOP50) over Greenland were selected based on average daily deposition rates over central Greenland (35–45 W; 70–80 N) to analyze the annual cycle of deposition. The TOP20 dust events were additionally analyzed using the LAGRANTO (Sprenger & Wernli, 2015) software package adapted for WRF model output to determine the source areas, transport pathways and transit times from the source areas toward central Greenland. For each of the TOP20 events, a set of 150 backward trajectories, starting at six equidistant levels (200m) up to 1,200 m above ground over central Greenland, were calculated and traced back to the dust source regions. The final set of trajectories per event include only those trajectories passing over a region where actually dust emission occurred at that time. This allows for an estimation of the major transport pathways and transit times of dust particles from Europe toward Greenland over the LGM.

## 2.6. Criteria of Literature Data Selection

New clay mineralogy and Sr–Nd–Hf isotope measurements of PSA samples and NGRIP dust were compared to existing data available in the literature. However, only a small subset of clay mineralogy and isotope data could be used for that purpose. With regard to XRD data, only those measurements were considered that were performed on the clay (<2 µm) fraction of sediments and where compositional data include all clay mineral groups (smectite, illite, chlorite) and kaolinite coming from (semi)quantified spectra. All literature and newly obtained XRD data are normalized to 100. Most of these clay mineral compositions were determined using the method of Biscaye (1965). In some cases, recent sediments are also included in the comparisons (see Dataset S1). Regarding the Sr–Nd–Hf isotope datasets, measurements undertaken on the aluminosilicate fractions of fine dust (<10, <5 or <2 µm) were exclusively considered (see Dataset S2). Isotopic data obtained on other size fractions or acid untreated samples were disregarded. Some of the isotope data are from modern sediments.

# 3. Results

## 3.1. Clay Mineralogy

Clay minerals from the illite group are the most abundant (x̄ ± 1SD = 54 ± 6%) phases in LGM dust samples from the GISP2 and GRIP ice core, followed by chlorite (x̄ ± 1SD = 27 ± 3%) and kaolinite (x̄ ± 1SD = 17 ± 5%), and with minor amounts of smectite (x̄ ± 1SD = 2 ± 2%) (Biscaye et al., 1997; Svensson et al., 2000; Dataset S1). Based on XRD data published in the literature and obtained in this study, PSA samples from Alaska, Siberia,
<figure><img src="image_4.png" /><figcaption>**Figure 2.** Ternary clay mineralogy diagrams of Greenland Ice Sheet Project 2 (GISP2)/Greenland Ice Core Project (GRIP) ice core dust and potential source areas in the (a) illite-smectite-kaolinite, (b) illite-kaolinite-chlorite and (c) chlorite/kaolinite-kaolinite/illite-smectite/kaolinite space. Note that the XRD data are normalized to 100 percent in each ternary plot. Dots without rims and fields in the background indicate published literature data (<2 µm fractions), while dots with rims denote newly acquired XRD data obtained in this study. Geographic regions represented by only one datum (e.g., NE Asia) is not defined by fields, but a rimless dot. Error bars on panels a) and b) represent a general 10 wt% uncertainty of XRD determination.</figcaption></figure>

Central and SE Asia have similar clay mineral compositions (Figures 2a and 2b; Dataset S1). Other PSAs show more heterogeneous clay mineralogy with samples resembling those of central Greenland ice core dust (some glacial continental North American samples and modern fluvial sediments from Western/Central/East Central Europe and North Africa) and others with smectite-richer compositions (West Africa, most of North Africa, East Central/Eastern Europe and continental US loess; Figure 2a). A crucial feature is that the GISP2 and GRIP last glacial dust samples have relatively limited variations in terms of clay mineral ratios and plot in the right corner of the chl/kao-kao/ill-sme/kao ternary diagram (Figure 2c), resulting from a smectite-poor composition. Few PSA samples have similar clay mineral ratios, but include Central Asia loess (only those from the literature), SE Asia loess and desert dust (region B), some Alaska and continental US samples and a small number of modern fluvial sediment samples from Western/East Central Europe. A multivariate statistical distance (Text S6 in Supporting Information S1), the so-called squared Mahalanobis distance (Md²) values confirm this picture (Figure S8a in Supporting Information S1) and these sources, most compatible with those of ice core dust, are characterized by Md² < 20–30.

## 3.2. Sr, Nd, Hf, and H Isotope Compositions

The aluminosilicate fractions of LGM dust samples from the NGRIP ice core have narrow ranges in Sr–Nd isotope compositions of 0.719557 – 0.720359 and –10.74 to –10.00 εNd(0). The ⁸⁷Sr/⁸⁶Sr isotope ratios are slightly more radiogenic than those measured for LGM dust from the GISP2 and GRIP ice cores. By contrast, the εNd(0) values are marginally less radiogenic, but in general the respective values overlap (Figure 3a; Dataset S2). One dust sample of the GISP2 ice core (G2) has outlying Sr–Nd isotope compositions, very likely caused by volcanic contributions. The ⁸⁷Sr/⁸⁶Sr isotope ratios of the leachates of NGRIP dust are lower (0.711551–0.712138) than those of the aluminosilicate fractions, reflecting the less radiogenic Sr isotope composition of sea salt and carbonates removed by acetic acid leaching. The opposite holds true for Nd, for which the removed fraction is relatively more radiogenic (–8.11 to –7.46) than the aluminosilicate one (Dataset S2). Since the isotope composition of the sea salt/carbonate fraction is meaningless for source discrimination of the insoluble mineral dust fraction, these are not considered further in this study. Numerous PSA samples fall in the Nd isotope compositional range of central Greenland ice core aerosols (εNd(0): about –11 to –9), while a few of them overlap or are close in Sr isotope ratios (Figure 3a). The most compatible PSA samples have Md² values below 8 (Figure S8b in Supporting Information S1) and include samples from SE Asia desert region B and Chinese loess, Northern Mongolia Plateau, Central Asian, continental North American (Nebraska), North African and Eastern/East Central European dust/loess. Some PSAs are much more (Central Asia, SE Asia desert region A, Yukon) or less
<figure><img src="image_5.png" /><figcaption>**Figure 3.** Sr-Nd (a) and Nd-Hf (c) isotope compositions of central Greenland last glacial dust and potential source area samples with mixture models in the Sr-Nd (b) and Nd-Hf (d) isotope space. Fields and rimless dots indicate literature data (list of data sources in Dataset S2), while dots with rims denote new isotopic data obtained in this study. End-members displayed on panels (b and d) are loess from Luochuan (Luo, CLP, China) and pumice fall deposit of the Khangar volcano (KHG, Kamchatka, Russia). Mixing lines are calculated with the equations given in Faure and Mensing (2005) and using elemental concentrations of Luo and KHG determined by ICP-MS and shown in Dataset S2. Measured Sr-Nd-Hf isotope compositions of artificial mixtures of KHG and Luo (proportions of 10%–90% and 30%–70%) are displayed in panel (b and d).</figcaption></figure>

(West Africa, Alaska, SE Asia desert region C) radiogenic in Nd, others are less (Egypt in North Africa, many SE Asian regions and parts of West Africa) or more radiogenic in Sr isotopic compositions (CLP, Western and Central Europe, parts of North Africa). Volcanic material from Kamchatka (samples KHG and KG) gave unra-diogenic ⁸⁷Sr/⁸⁶Sr ratios (about 0.7033) and very radiogenic Nd isotope compositions (εNd(0): +8–9). Artificial
<figure><img src="image_6.png" /><figcaption>**Figure 4.** Box/scatter-(a) and bi-plots (b–c) of hydrogen isotope compositions of the clay structural water ($\delta D_{csw}$) of North Greenland Ice Core Project (NGRIP) last glacial dust and potential source area (PSA) samples as a function of Nd-Hf isotope compositions. The NGRIP mean $\delta D_{csw}$ value is -71‰ with a $\pm 9\%$ band derived from the repeatability of 2–5 µm separates of PSA samples (see Text S4 and Figure S6 in Supporting Information S1). The $\delta D_{csw}$ of both 2 and 5 µm separates are included in boxes of panel (a) and can be found in Dataset S2. In panel (b–c) all hydrogen isotope data are displayed as measured on both 2 and 5 µm separates, where corresponding Nd and Hf isotope ratios were available from the 5 µm fractions (Dataset S2).</figcaption></figure>

10:90 mixture of two end-members (KHG tephra and Luo loess, China) yields Sr–Nd isotope ratios close to ice core dust, while the mixture of 30:70 gave compositions similar to sample G2 of the GISP2 ice core (Figure 3b; Dataset S2).

Hafnium isotope composition of the aluminosilicate fractions of LGM dust samples of the NGRIP ice core range from $\epsilon$Hf(0) of -7.06 to -4.67 (Dataset S2), overlapping with values of two dust samples (132A, B) of the Dye-3 ice core preindustrial segment (Lupker et al., 2010). These Hf isotope compositions are in agreement with those obtained from some East Central/Eastern European, North African, NE Asian and continental US samples, also indicated by Md² values <10 (Figure S8c in Supporting Information S1). Most of the SE Asian dust samples, including desert region B, Northern Mongolia Plateau and CLP/Beijing loess, are more radiogenic in Hf, with $\epsilon$Hf(0) values varying between -3 and +3 (Figure 3c). The few exceptions include aeolian deposits from Taklimakan and Tengger deserts with $\epsilon$Hf(0) values between -5.17 and -3.71 (samples TD, TD2 and Td25, Dataset S2), which show values at the higher end of ice core dust values. Central Asia loess deposits are clearly more radiogenic in Hf ($\epsilon$Hf(0): -3.03 to -1.99) than NGRIP ice core dust, likewise Yukon loess with $\epsilon$Hf(0) values of +3.82 to +3.96, reflecting volcanic dust additions to these sediments. The KHG and KG pumice fall deposits yield, as expected, highly radiogenic Hf isotope compositions ($\epsilon$Hf(0): +13.97 to +17.03) and the 10:90 and 30:70 mixtures of Chinese loess (Luo) and the KHG sample gave $\epsilon$Hf(0) values of -1.22 and 1.00 (Dataset S2) and both plot far away from NGRIP dust in the Nd–Hf isotope space (Figure 3d).

Hydrogen isotope measurements of clay structural water of two LGM dust samples of the NGRIP ice core yield $\delta D_{csw}$ values of -72.2 and -70.4 ‰ (Dataset S2). Most PSA samples have different, usually more negative $\delta D_{csw}$ values with a range from about -100 to -80 ‰, including Central Asia, SE Asian loess (CLP, around Beijing), Eastern Europe and continental North America (Nebraska loess). The most negative $\delta D_{csw}$ values (-116 to -101) are characteristic for dust samples originating from the northernmost sources (Yukon loess and NE Siberia loess), while the less negative values (-67 to -62 ‰) from North Africa and SE Asian desert region B (Figure 4a). This geographical pattern indicates a latitudinal dependence of $\delta D_{csw}$ in the studied samples and is reminiscent of the modern atmospheric water vapor and rainfall depleted in deuterium toward the poles (Craig, 1961; Dansgaard, 1964). None of the PSA samples match perfectly the NGRIP dust $\delta D_{csw}$ values, but those being close in composition include European loess, SE Asian desert sediments (region B, Taklimakan), CLP loess and North African dust. A small subset of these samples (e.g., Zm from East Central Europe or Td25 from the Taklimakan desert) plot near the NGRIP last glacial aerosols in $\epsilon$Nd(0) versus $\delta D_{csw}$ and $\epsilon$Hf(0) versus $\delta D_{csw}$ spaces (Figures 4b and 4c), and have Md² values of about 16–17 in Nd–Hf-$\delta D_{csw}$ (Dataset S2, Figure S8d in Supporting Information S1).
<figure><img src="image_7.png" /></figure>

**Figure 5.** (a) Total simulated dust deposition (g m⁻²) over Greenland accumulated over 30 model years. (b) Trajectories (black lines), total dust load (shading, mg m⁻²) and geopotential height (gpdm) at 700 hPa during an individual dust event in July. (c) Density of emission source points found from backward trajectories based on TOP20 dust deposition events over Greenland when they reach an emission area for the first time (trajectory points per 1°), based on 1,023 trajectory points in total. (d) Trajectory density (percentage of trajectory points per 2° radius) for all TOP20 dust deposition events, based on 24,404 trajectory points in total. Emission areas hatched. Numbering in panels (b–d): 1. Basaharc (Bh, Hungary), 2. Crvenka (Crv, Serbia), 3. Dunaszekcső (Dsz, Hungary), 4. Mende (Me, Hungary), 5. Paks (Pa, Hungary), 6. Titel core site (Tic, Serbia), 7. Urluia (Url, Romania), 8. Zmajevac (Zm, Croatia), 9. La Motte (LMJ, Green Island, Jersey, Channel Islands), 10. Krems-Wachtberg (KW, Austria), 11. Nussloch (Nus, Germany), 12. Bialy Kosciol (Bk, Poland), 13. Pegwell Bay site (PbE, UK). Red line depicts Last Glacial Maximum land mask based on Paleomodel Intercomparison Project3 guidelines for the 21k experiment.

### 3.3. Dust Transport Simulations

The 30 year simulations of the LGM dust cycle performed in this study with WRF-Chem demonstrates that mineral dust emitted from European sources could reach central Greenland during the LGM. Most dust events were found to be triggered by easterly winds, induced by an anticyclonic circulation around a semi-permanent high-pressure system over the Fennoscandinavian ice sheet (Ludwig et al., 2016). Dust particles were then transported toward the west and became involved in the cyclonic circulation of low-pressure systems over the North Atlantic, which finally directed dust toward Greenland, where deposition took place (Figures 5a–5c). Another dust transport route to the east of the Fennoscandinavian high-pressure center is also identified in the model simulations (Figures 5d and 6 and Figure S7 in Supporting Information S1), albeit with subordinate occurrence. Interestingly, trajectories reached back to North Africa for one of the 20 most intense (TOP20) dust events (Figure 6). Major source regions of the TOP20 dust events include the Po Plain and the NW Adriatic Sea region, the Carpathian Basin, the northern part of the East European Plain, and intense dust emission hotspots from the North European Plain to the North Sea continental shelf (Figures 5c and 6). The TOP20/TOP50 dust events
<figure><img src="image_8.png" /><figcaption>**Figure 6.** Backward trajectories for the TOP20 deposition events and dust emission areas. Event number and month indicated at the bottom left for each case. Shaded: Mean dust emission [µg m⁻² s⁻¹] between the first and last time a trajectory passed over a region with dust emission for each event.</figcaption></figure>

have two seasonal maxima during winter/early spring and late summer/early autumn (Table S3 in Supporting Information S1). Mean transit time calculated from model outputs is 4.68 days for the TOP20 dust events, with the shortest and longest transit times of 2.0 and 7.75 days (Table S4 in Supporting Information S1). The total simulated dust deposition over central Greenland for the 30 model years is ~0.53 g m⁻² (Figure 5a), translating to a dust deposition flux of 17.7 mg m⁻² yr⁻¹.
## 4. Discussion

Source discrimination of aerosols emitted from potentially multiple continental areas, entrained at various heights in the atmosphere and subsequently deposited on polar ice sheets in the past is inherently complex. The reliability of provenance interpretations largely depends on the dust source tracers applied, the general understanding of the major controls on these indicators, and the analytical uncertainties and reproducibility of tracer data. We first discuss these issues prior to considering the interpretation of these data in terms of dust source.

### 4.1. Uncertainties Related to Dust Tracers

In this study, one of these dust fingerprints used in source identification is the clay mineralogy of PSAs and ice core dust. The distribution of clay minerals on the continents is controlled by the local lithology, drainage characteristics, weathering and climate, and it is latitude dependent (Biscaye, 1965; Biscaye et al., 1997; Griffin et al., 1968). Therefore, each dust source has a characteristic spectrum of clay minerals and the abundances and abundance ratios may be source diagnostic. It is of note that the XRD intensity-based calculation of the composition of a given sample in weight percent of mineral phases is semi-quantitative, with results being affected by ±10% uncertainties for major and ±20% for minor constituents (<20 wt.% amount) and with detection limits of about 1 wt.%–5 wt.% (Moore & Reynolds, 1997).

Other traditional dust provenance indicators include the 87Sr/86Sr and 143Nd/144Nd isotope ratios of the carbonate-free, aluminosilicate fractions of dust (Grousset & Biscaye, 2005). The Sr isotope composition of aeolian dust is controlled by carbonates, feldspars and various clay minerals (Brass, 1975; Capo et al., 1998) and there is a strong coupling between Sr isotopic ratios and grain size in sediments (Chen et al., 2007; Dasch, 1969; Feng et al., 2009; Meyer et al., 2011). In ice core dust source discrimination using Sr isotopes, calcite and gypsum aerosols being entirely of marine origin are removed by acid treatment (Biscaye et al., 1997; Svensson et al., 2000), as these phases are capable of masking the Sr isotope composition of the aluminosilicate (or lithic) fraction that is of interest here (Delmonte et al., 2004). Furthermore, the <10 or <5 µm particle size fractions representative of the grain size of ice core dust are compared (Figure S2 in Supporting Information S1; Biscaye et al., 1997; Bory, Biscaye, Piotrowski, & Steffensen, 2003; Svensson et al., 2000), thereby minimizing grain size effects. However, the use of different acids (hydrochloric or acetic acid) in various concentrations and different reaction times may strongly influence the measured Sr isotope ratio due to the partial (calcite) or complete (calcite/dolomite) removal of carbonates and even a fraction of clay minerals (Asahara et al., 1999; Yokoo et al., 2004). Re-measurements of new separations of loess samples previously prepared and measured for Sr isotope ratios by Újvári et al. (2015) demonstrate that subtle deviations in grain size (Figures S2a in Supporting Information S1) may result in significant differences in 87Sr/86Sr isotope ratios (0.0032–0.0059; Figure S3a in Supporting Information S1). This implies that extreme care must be exercised when the Sr isotope compositions published in different studies are to be compared and it is important to strictly follow the same methodology (acid treatment, grain size fraction) in all these studies to reach meaningful interpretations using Sr isotopes. By contrast, 143Nd/144Nd isotope compositions, controlled by allanite, apatite, garnet, monazite, sphene, zircon and also plagioclase, amphibole, biotite and clay minerals (Aubert et al., 2001; Viers & Wasserburg, 2004), are much less or not at all affected by these issues (0.000010–0.000026 or 0.2–0.5 εNd units, Figure S3b in Supporting Information S1). These 0.2–0.5 εNd unit variations are comparable to the analytical reproducibility (about 0.3 εNd units) of Nd isotope measurements as determined by replicate analyses of three aeolian sediment samples (Figure S4 in Supporting Information S1).

Hafnium isotope ratios have only been used once, to our knowledge, in ice core dust source discrimination on the Dye-3 ice core (Lupker et al., 2010), partly because of analytical issues, low amount of ice core dust available and the scarcity of 176Hf/177Hf data from the fine fractions of PSA samples. Bulk samples of wind-blown sediments are dominated by zircons having extremely high Hf contents (Hoskin & Schaltegger, 2003) and unradiogenic 176Hf/177Hf ratios (Újvári & Klötzli, 2015), while the clay fractions contain little zircon (Újvári et al., 2018). This is further evidenced in Figure S16 of Supporting Information S1 (Text S8 in Supporting Information S1), where zircon depletion toward the finer grain size fractions of Quaternary aeolian loess, paleosol and red clay samples is clearly visible in the Zr, Hf concentration and Nd–Hf isotope data sets. The more radiogenic Hf isotopic compositions of the clay fraction (Zhao et al., 2015) are controlled by clay minerals that incorporate radiogenic Hf released from higher Lu/Hf phases (apatite, garnet, xenotime and sphene) during incongruent silicate weathering (Bayon et al., 2009; van de Flierdt et al., 2007). In this study, the ammonium-bifluoride method
was used, which allowed for a fast and complete dissolution of a large number of 5–10 mg samples for Hf–Sr–Nd isotope analyses (Újvári et al., 2021), including NGRIP ice core dust. Reproducibility of Hf isotope analyses was found to be between 0.3 and 0.8 εHf units, translating to differences of about 0.000009–0.000023 in ¹⁷⁶Hf/¹⁷⁷Hf isotopic ratios based on replicate measurements of three aeolian dust samples (Figure S4 in Supporting Information S1). These are usually within or slightly higher than analytical errors of individual measurements.

The δDcsw values of the <5 and < 2 µm fractions of ice core dust and PSA samples are thought to be controlled by the D/H isotope compositions of various hydrous silicates in the clay-mineral size fraction of these samples. The D/H isotope composition of hydrous minerals depends on the isotopic composition of ambient water and the temperature at the time of mineral formation, as well as the mineral-specific isotopic fractionation factor at the temperature of formation (under isotopic equilibrium), and finally whether the mineral has retained its original isotopic composition (Savin & Hsieh, 1998). In the case of clay minerals formed by surficial weathering, the ambient water would be meteoric water and temperatures would approximate Earth's surface temperatures. The dust samples studied are thought to be primarily made up of detrital clay minerals, largely by clays formed in-situ by chemical weathering that may have a range of formation ages and origin. This suite of clay minerals with diverse D/H isotope compositions defines δDcsw values characteristic for the different PSAs. Indeed, there is a spatial pattern in δDcsw of PSA samples, with the most negative values for dust samples from higher latitudes (Figures 4a, Dataset S2), and some PSAs have a set of δDcsw values sufficiently well-defined to be diagnostic of their source. For the source discrimination to be meaningful, the reproducibility of D/H isotope measurements must be considered. This was found to be about 1–4 ‰ within the <5 and < 2 µm fractions and mostly range between 3 and 9 ‰ when these two fractions are compared (Figure S6 in Supporting Information S1). Thus, sources with δDcsw values in the range of –71.3 ± 9 ‰ (mean ± reproducibility of NGRIP dust, Dataset S2) must be considered as potential contributors to last glacial dust of central Greenland.

## 4.2. Last Glacial Greenland Dust Sources

### 4.2.1. Low(er) Probability Single Source Regions and Source Combinations

Based on the suite of dust tracers applied here and considerations outlined above, some PSAs have a low probability of being the source of last glacial dust to central Greenland. While Alaskan loess samples are among the most compatible with central Greenland last glacial dust in terms of clay mineralogy, direct sourcing (i.e., without mixing) of dust from Alaska has been excluded by Biscaye et al. (1997) based on Sr–Nd isotopic compositions. The Sr–Nd isotope ratios here indeed support the exclusion of the Alaskan source (Dataset S2), although more Alaskan sediments would need to be analyzed to further confirm this. The Yukon loess, as a single source, is likewise incompatible with central Greenland dust, as it has much less radiogenic Sr and more radiogenic Nd isotope compositions (Zdanowicz et al., 2006), positive Hf isotope ratios and very negative δDcsw values compared to Greenland dust (Dataset S2, Figures 3a and 3c and 4). The Sr–Nd–Hf isotopic compositions of Nebraska loess fine fractions are similar to central Greenland last glacial dust (Újvári et al., 2015; Dataset S2), but these sediments are usually extremely smectite-rich and poor in chlorite and kaolinite (Újvári et al., 2015, Dataset S1). Smectite aggregation and resulting fractionation during atmospheric transport (Scheuvens et al., 2013; Singer et al., 2004) would not lead to the almost total loss of such high amounts of smectites (about 70–75 wt%) and chlorite and kaolinite enrichment during atmospheric transport can likewise be excluded, unless mixed with another chlorite and kaolinite rich dust source. Although direct derivation of last glacial ice core dust from this PSA is thus unlikely, sourcing of ice core dust from the US drylands (e.g., Nevada; Aarons et al., 2017) based on Nd–Hf isotopic compositions cannot entirely be ruled out, but this requires a more systematic study of the fine, lithic fractions of last glacial dust deposits of that region.

Although direct dust transport to Greenland from single North American sources is unlikely based on the data here, the possibility of a contribution as part of a mixture with dust from other sources must also be considered, especially because the atmosphere over Greenland was shown to be well-mixed during colder periods of the last glacial cycle (Mayewski et al., 1994). Indeed, given the previously modeled SE Asian dust transport pathway to Greenland (Andersen et al., 1998; Mahowald et al., 2011), North American glacially sourced aerosols could potentially have mixed with SE Asian dust en route to Greenland. However, given the very negative δDcsw values of these sources (SE Asia: mostly <–75 and North America: <–90 ‰) and the NGRIP dust average (–71.3 ‰), this possibility is quite unlikely. In terms of δDcsw values, a low latitude source with higher than ~–70.0 ‰ such as North Africa would be needed to reproduce NGRIP ice core dust compositions. However, North
African transport routes were northward during the glacial, as shown in our modeling results detailed later (see Figure 6), meaning that mixing of SE Asian-North American and North African dust in a single pathway is implausible. Despite this, our Monte Carlo (MC) simulations reveal that the Sr-Nd-Hf isotopic compositions of central Greenland last glacial dust is relatively well explained by mixing North African dust with Nebraska loess (Figure S9a and S9b in Supporting Information S1), in contrast to the Yukon loess, although the $\delta D_{csw}$ values are more difficult to reproduce (Figure S9c and S9d in Supporting Information S1). However, it is important to emphasize that all four isotopic compositions must be considered together in assessing plausible mixing ratios for dust admixtures. Indeed, mixing ratios that simultaneously reproduce all the isotopic compositions (Sr-Nd-Hf and $\delta D_{csw}$) of ice core dust well are very few, with 10 and 2 such combinations found out of 1,000 modeled cases for the Nebraska and Yukon loess, respectively (Table S7 in Supporting Information S1). Of those few mixing ratios, none yield modeled clay mineral compositions compatible with last glacial ice core dust; all show much more smectite, less illite and kaolinite and much less chlorite compared to ice core dust (Figure S14 in Supporting Information S1). For these reasons, and based on the currently available data, we believe that the likelihood of North American sources contributing significantly to the last glacial ice core dust deposited in central Greenland, even when mixed with other sources, is low.

Central Asian dust sources, represented by loess from Tajikistan, Kazakhstan and western China (Ili Basin), are also considered unlikely as direct single contributors of aerosols to central Greenland over the LGM due to their more radiogenic Nd-Hf isotope compositions and low $\delta D_{csw}$ values (Figures 3a and 3c and 4, and Figure S8d in Supporting Information S1). The NE Asian (Siberia) region is underrepresented by only one loess sample, which has clay mineralogy and Nd-Hf isotope compositions (Figures 2 and 3 and Figure S8a, S8b) resembling central Greenland ice core dust. However, its more radiogenic $^{87}$Sr/$^{86}$Sr isotope ratio and very negative $\delta D_{csw}$ values (down to -115 ‰; Figure 4) suggest that NGRIP last glacial dust was unlikely to be directly derived from this Siberian source.

These indications of North American/Siberian sources being at least partly incompatible with central Greenland last glacial aerosol compositions imply that despite being important in simulations, certain modeled LGM dust sources (Lambert et al., 2015; Mahowald et al., 2006; Werner et al., 2002) do not seem to contribute dust to central Greenland. A possible explanation for the incompatibility with North American sources, despite the large interior continental planes and upwind location to Greenland, could be that the North Atlantic jet stream and storm tracks were more zonally oriented and less variable at the LGM compared to today (Löfverström, 2020; Löfverström et al., 2016).

### 4.2.2. High(er) Probability Single Source Regions and Source Combinations

Derivation of last glacial ice core dust from North Africa was excluded by early studies (Biscaye et al., 1997; Svensson et al., 2000), but later considered plausible (Han et al., 2018; Meyer et al., 2017). In general, western/ central Saharan dust samples have abundant kaolinite and smectite (Dataset S1, Figure 2; Scheuvens et al., 2013), while some recent fluvial sediment samples around the High Atlas in Morocco have smectite-poor compositions similar to LGM ice core dust. There is a wide variety in isotopic compositions of North African/Saharan dust samples with $^{87}$Sr/$^{86}$Sr of 0.730 to 0.708, $\epsilon$Nd(0) of -18.5 to -4 and $\epsilon$Hf(0) values of -13.7 to +3.5, controlled by lithologies of the West African craton (Mali) to young volcanic rocks in Egypt (Abouchami et al., 2013; Grousset & Biscaye, 2005; Zhao et al., 2018). Our isotopic data, together with those presented by Zhao et al. (2018), do not exclude the possibility of sourcing of, or at least an admixture to, ice core dust from the northernmost part of Africa (Morocco/Tunisia) (e.g., sample Mao, Datasets S1-S2). This is especially intriguing in the light of our model simulation results, which include some backward trajectories over Europe reaching back to North Africa (Tunisia, Figure 6, Figure S7 in Supporting Information S1).

Out of the analyzed PSAs, some specific regions of SE Asia and Europe are the most plausible single sources of last glacial dust of central Greenland ice cores. Clay mineralogical compositions of most SE Asian loess/desert sediment samples match relatively well those of last glacial dust from the GISP2/GRIP ice cores (Figure 2, Table S5 in Supporting Information S1), usually with higher illite content of up to 77% (Taklimakan, Tengger desert) and less abundant chlorite (Dataset S1). However, few SE Asian sources are compatible with last glacial ice core dust in terms of isotope compositions. The SE Asian desert region A (northern deserts, including the Gobi) and the Northern Mongolia Plateau samples are usually less radiogenic in Sr isotope and more radiogenic in Nd isotope composition, and loess deposits from the CLP and around Beijing have higher $^{87}$Sr/$^{86}$Sr (Figure 3a). This latter discrepancy was explained by volcanic dust additions to SE Asian aerosols and mixing in 10/90%
proportions (Biscaye et al., 1997), which is tenable based on the Sr–Nd isotope data (Figure 3b). However, some dust samples from the Taklimakan and Tengger desert (desert region B) have overlapping Sr–Nd isotope compositions with GISP2/GRIP (partly by NGRIP) dust, and thus, volcanic dust admixture is not required to explain ice core dust $^{87}$Sr/$^{86}$Sr and $^{143}$Nd/$^{144}$Nd ratios.

Hafnium isotope compositions put further constraints on plausible sources, as most PSAs of SE Asia have much more radiogenic $^{176}$Hf/$^{177}$Hf isotope ratios ($\varepsilon$Hf(0): –3 to +4) than NGRIP ice core dust ($\varepsilon$Hf(0): –7.06 to –4.67). Few exceptions include Tengger (TD/TD2) and Taklimakan (Td25) desert samples with $\varepsilon$Hf(0) values of –5.17 to –3.71, being at the higher end of NGRIP ice core dust compositions (Figure 3c). Importantly, the Hf isotope ratios also demonstrate that the previously proposed mixing model of circum-Pacific volcanics and SE Asian dust is untenable. This confirms what has already been proposed by Svensson et al. (2000) that an almost constant (10%) contribution of volcanic material, inferred from the restricted range of Sr isotope ratios of central Greenland dust, seems unlikely considering the stochastic nature of volcanic eruptions. Measurements of Kamchatka volcanics (KHG and KO tephra) and artificial mixtures of the KHG tephra, which has been identified in Holocene ice of the NGRIP core (Cook et al., 2018), and loess from the Luochuan loess-paleosol sequence (Luo) indicate that an admixture of volcanic material to SE Asian dust would further increase the already radiogenic $^{176}$Hf/$^{177}$Hf isotope ratios of SE Asian dust, shifting the composition of the mixtures in the opposite direction compared to those of the NGRIP last glacial dust (Figure 3d). Thus, this mixing model can safely be excluded as an explanation for ice core dust composition.

By contrast, direct sourcing of last glacial dust of central Greenland from the Taklimakan/Tengger deserts is supported by the D/H isotope compositions. The $\delta D_{esw}$ values of NGRIP dust are different from most SE Asian dust samples (Figure 4), but within reproducibility overlap with one sample from the Taklimakan (Td25) and another from the CLP (Dataset S2). Additional support for SE Asian deserts being the major contributors of dust to the Greenland Ice Sheet over the LGM comes from the fact that they are upwind to Greenland and their high altitude, facilitating uplift of dust into the long-range transporting, westerly wind regime (Svensson et al., 2000). Further evidence is the origin of volcanic material in Greenland ice. Tephra analyses of Greenland ice cores indicate that beyond the clear dominance of relatively proximal Jan Mayen and Icelandic sources, volcanic ash from the Pacific Arc was repeatedly entrained to and deposited in Greenland, whereas no ash particles from central Europe have been identified in Greenland ice cores despite that region being very active in volcanism during the last glacial period (Abbott & Davies, 2012; Bourne et al., 2016).

An alternative hypothesis for the origin of last glacial dust in Greenland is the derivation from European sources (Újvári et al., 2015), which has recently been dismissed by Han et al. (2018) based on Sr and Pb isotopic considerations. Unfortunately, the rejection of this hypothesis was based on bulk sediment Sr–Pb isotopic data from Europe (Han et al., 2018) disregarding the fact that both Sr (Chen et al., 2007; Dasch, 1969; Feng et al., 2009; Meyer et al., 2011) and Pb isotope compositions (Feng et al., 2010) are grain size-dependent. The large variability in Sr–Pb isotope compositions between the bulk and clay fractions (e.g., $\Delta^{206}$Pb/$^{207}$Pb and $\Delta^{208}$Pb/$^{207}$Pb: 0.01–0.04; Figure S17 in Supporting Information S1) reaches and may even exceed differences in isotopic compositions between ice core dust and European loess, thereby introducing a significant uncertainty and preventing meaningful source identification. As such, the data presented in Han et al. (2018) cannot be used in testing European dust sources to Greenland. Our new data further suggest that the hypothesis of European dust transport to Greenland cannot be ignored. The clay mineralogy of most glacial loess samples from Europe are relatively smectite-rich (Figure 2). This is especially true for East Central/Eastern European loess dominated by local sources of the Carpathians/Carpathian Basin (Fu et al., 2021), also seen in modern river sediments of the Danube catchment (Martinez-Lamas et al., 2020). However, tributaries of the Danube River draining the Alps show smectite-poor compositions close to those of ice core dust (Figure 2, Figure S8a in Supporting Information S1) and some loess deposits, for instance along or near the Drava River, also have lower smectite contents of about 20%–25% (e.g., sample Zm, Dataset S1). Considering the uncertainties of semi-quantitative XRD data and that smectite fractionation during transport cannot be ruled out (see discussion above), derivation from East Central Europe is still a possibility. This is even more so because Taklimakan dust has about 20% higher illite content compared to last glacial dust from central Greenland ice cores, so the direct sourcing of ice core dust from the Taklimakan would likewise require illite loss during atmospheric transport. Another option, assuming that the particulate matter retains its original mineralogical composition during atmospheric transport, is that the dust mixed with other sources to produce the final composition of the dust in the ice cores. This situation clearly demonstrates the
ambiguities related to the use of mineral tracers alone, which in many cases is misleading (Fu et al., 2021) and does not allow for robust source discrimination. The multi-tracer approach of this study reveals this complexity and highlights the importance of using isotope tracers as well. In Sr–Nd isotopic space, some East Central and Eastern European samples plot close to central Greenland ice core dust (Figure 3a, S8b in Supporting Information S1), while the majority of samples are more radiogenic in $^{87}$Sr/$^{86}$Sr. Part of the East Central/Eastern European loess samples overlap with the NGRIP dust samples in Hf isotopic compositions, while those from Western (ice marginal) Europe are usually less radiogenic compared to ice core dust (Figure 3c, Table S5 in Supporting Information S1). In terms of D/H isotope compositions, glacial dust from Western and East Central Europe overlap the range measured for NGRIP last glacial dust (Figure 4), and some European glacial loess samples (especially sample Zm, Datasets S1–S2) are very similar to NGRIP last glacial aerosols in Sr–Nd–Hf and D/H isotope compositions (Figure S8c and S8d in Supporting Information S1).

Beyond this tracer evidence for potential European dust contributions, recent studies demonstrate higher dust mass accumulation rates in Europe than in China over the LGM (Rousseau et al., 2021) and that the variations in glacial dust deposition on centennial–millennial timescales in Europe and Greenland were synchronous within uncertainty on radiocarbon chronologies (Moine et al., 2017; Újvári et al., 2017). More southerly positions of the upper-tropospheric, eddy-driven Atlantic jet during cold, dusty stadials (Löfverström et al., 2016; Újvári et al., 2017) in combination with anticyclonic circulation over the Fennoscandian Ice Sheet induced by a semi-permanent high-pressure system (Ludwig et al., 2016), would have promoted dust transport from European sources to Greenland. At the same time, dust flux estimates based on compiled optically stimulated luminescence ages from multiple Chinese and Central Asian loess sites suggest high dust accumulation during the late LGM (23–19 ka) and much lower for the early LGM (26.5–23 ka) in these areas (Cheng et al., 2021; Kang et al., 2015), in contrast to what has been measured in Greenland ice cores (Rasmussen et al., 2014). However, a reduction or absence of dust accumulation in an area does not necessarily imply lower dust production and emission by the same area. Indeed, previous research argues that the CLP is a highly dynamic environment which may have resulted in substantial internal aeolian recycling and erosion of pre-deposited material (Kapp et al., 2015; Licht et al., 2016; Stevens et al., 2018).

Results from our high spatial resolution regional model simulation further reduce the ambiguities about the possibility of dust transport from Europe to Greenland over the LGM, and clearly demonstrate that aerosols emitted from several European glacial dust hotspots could reach the Greenland Ice Sheet in all seasons, albeit with varying seasonal occurrence (Table S3 in Supporting Information S1, Figure 6), resulting in a dust deposition flux of 17.7 mg m⁻² yr⁻¹ in central Greenland. This is, however, much lower than the GS-2.1 stadial dust flux for the 21–22 ka period (~102–203 mg m⁻² yr⁻¹, e.g., Serno et al., 2015), for which the model was run, and the GS-3 stadial peak dust deposition at 25.8 ka (498 mg m⁻² yr⁻¹), while it is comparable with the GI-3 and 2.2 interstadial values of ~20 mg m⁻² yr⁻¹. The lower simulated dust flux values may be due to the fact that the strength of European dust sources and dust emission values in the model have not been calibrated with actual dust accumulation values, for example, using observed loess sediment accumulation rates. However, this does not affect the main dust transport pathways or their characteristics such as directions and seasonal distributions. The modeled transit times of 2.0 and 7.75 days (mean of 4.68 days) of the TOP20 events are obviously much shorter than those considered in earlier studies for the LGM transport of calcium (Ca²⁺) from East Asian sources (Schüpbach et al., 2018), and would have important implications for interpreting the variability of the Ca²⁺ record of ice cores in terms of significance of drivers either related to sources (strength, storminess, aridity) or atmospheric transport (residence time, wet/dry deposition, transport distance) (Werner et al., 2002). Modeled trajectories crossing Europe and reaching back to North Africa imply that mixing of dust from these sources may have been possible, further complicating the already complex problem of ice core dust provenance Indeed, our MC simulations with two component mixtures of North African and European dust demonstrate that both the isotopic compositions (Figure S10, Table S8 in Supporting Information S1) and clay mineralogy (Figure S15 in Supporting Information S1) of last glacial dust from central Greenland ice cores can be reproduced. Further simulations indicate that the isotopic compositions of last glacial ice core dust can be explained relatively well by the mixing of dust from SE Asian desert regions and the CLP (Figure S11 in Supporting Information S1) or even three distinct PSAs of the NH, including the Taklimakan, European, and North African sources (Figures S12–S13 in Supporting Information S1). In the present study, the PSAs considered most likely by previous geochemical and modeling studies have been investigated, but we cannot exclude that others may have existed. Also, the possibility that aerosols from more than two or three distinct sources may have been mixed also cannot be excluded, but modeling of such
mixtures has not been undertaken in the present study as this would require additional dust tracer datasets from even more PSAs of the NH.

## 5. Conclusions and Outlook

We present here the most in-depth and multi-proxy analysis of last glacial dust sources to Greenland to date, supported by a regional model simulation of the dust cycle, including the computation of dust trajectories. Our datasets reveal that previous claims of well-constrained, specific Greenland dust sources are not supported. Rather, the data reveal three plausible, competing hypotheses for the origin of LGM dust recovered from central Greenland ice cores. Provided that ice core dust is from single sources, still the most likely scenario is that these aerosols are directly derived from the Taklimakan and/or Tengger desert, while other SE Asian sources are excluded. However, based on the compositional and modeling data presented, sourcing of ice core dust from Europe is also plausible, and furthermore, an admixture with North African sources cannot be ruled out. To gain more insight into Greenland ice core dust provenance, further, wide-scale investigation of many more PSA samples would be required, together with high spatial resolution dust cycle simulations for the entire NH over the LGM. These models should preferably be calibrated with observed dust accumulation rates, for example, using absolutely dated loess deposits on the continents (Kohfeld & Harrison, 2001; Mahowald et al., 2006; Stevens et al., 2016; Újvári et al., 2010, 2017), to accurately simulate terrestrial dust source strengths and deposition flux in central Greenland. In addition, further isotope geochemical studies of other sections of the central Greenland ice cores with lower dust accumulation would be worthwhile to better understand dust-climate feedbacks. In part, such studies have already been carried out (Han et al., 2018; Svensson et al., 2000) and found low variability in dust Sr-Nd-Pb isotopic compositions over time, suggesting that if Greenland had had multiple dust sources, their relative contributions have been fairly well preserved during the abrupt climate transitions of the last glaciation. However, additional information from as many isotopic dust tracers as possible is still required from specific periods of the last glaciation. As such, targeting intervals of high NH summer insolation, when a low Sahara-Sahel signal is expected ("Green Sahara" phases/humid North Africa; e.g., Tjallingii et al., 2008) and targeting lows in NH summer insolation when the opposite would be expected, could be a way of better quantifying the potential contributions from north African sources and deconvolving the dust provenance signal in the central Greenland ice core records.

## Data Availability Statement

The data that supports the findings of this study are published in the main text and Supporting Information S1 of this article and are openly available in the PANGAEA Data Repository (Újvári, 2021).

## References

Aarons, S. M., Blakowski, M. A., Aciego, S. M., Stevenson, E. I., Sims, K. W. W., Scott, S. R., & Aarons, C. (2017). Geochemical characterization of critical dust source regions in the American West. *Geochimica et Cosmochimica Acta*, **215**, 141–161. https://doi.org/10.1016/j.gca.2017.07.024

Abbott, P. M., & Davies, S. M. (2012). Volcanism and the Greenland ice-cores: The tephra record. *Earth-Science Reviews*, **115**(3), 173–191. https://doi.org/10.1016/j.earscirev.2012.09.001

Abouchami, W., Näthe, K., Kumar, A., Galer, S. J. G., Jochum, K. P., Williams, E., et al. (2013). Geochemical and isotopic characterization of the Bodélé Depression dust source and implications for transatlantic dust transport to the Amazon Basin. *Earth and Planetary Science Letters*, **380**, 112–123. https://doi.org/10.1016/j.epsl.2013.08.028

Albani, S., Balkanski, Y., Mahowald, N., Winkler, G., Maggi, V., & Delmonte, B. (2018). Aerosol-climate interactions during the Last Glacial Maximum. *Current Climate Change Reports*, **4**(2), 99–114. https://doi.org/10.1007/s40641-018-0100-7

Alley, R. B., Marotzke, J., Nordhaus, W. D., Overpeck, J. T., Peteet, D. M., Pielke, R. A., et al. (2003). Abrupt climate change. *Science*, **299**(5615), 2005–2010. https://doi.org/10.1126/science.1081056

Andersen, K. K., Armengaud, A., & Genthon, C. (1998). Atmospheric dust under glacial and interglacial conditions. *Geophysical Research Letters*, **25**(13), 2281–2284. https://doi.org/10.1029/98GL51811

Asahara, Y., Tanaka, T., Kamioka, H., Nishimura, A., & Yamazaki, T. (1999). Provenance of the north Pacific sediments and process of source material transport as derived from Rb-Sr isotopic systematic. *Chemical Geology*, **158**(3–4), 271–291. https://doi.org/10.1016/S0009-2541(99)00056-X

Aubert, D., Stille, P., & Probst, A. (2001). REE fractionation during granite weathering and removal by waters and suspended loads: Sr and Nd isotopic evidence. *Geochimica et Cosmochimica Acta*, **65**(3), 387–406. https://doi.org/10.1016/S0016-7037(00)00546-9

Bauer, K., & Vennemann, T. W. (2014). Analytical methods for the measurement of hydrogen isotope composition and water content in clay minerals by TC/EA. *Chemical Geology*, **363**, 229–240. https://doi.org/10.1016/j.chemgeo.2013.10.039

Bayon, G., Burton, K. W., Soulet, G., Vigier, N., Dennielou, B., Etobleau, J., et al. (2009). Hf and Nd isotopes in marine sediments: Constraints on global silicate weathering. *Earth and Planetary Science Letters*, **277**(3–4), 318–326. https://doi.org/10.1016/j.epsl.2008.10.028

## Acknowledgments

This study was financially supported by the Austrian Research Fund (FWF) through a Lise Meitner grant (project no. M 2503-N29), and the European Union and the State of Hungary, co-financed by the European Regional Development Fund in the project of GINOP-2.3.2.-15-2016-00009 "ICER." PL was partially funded by the Helmholtz Climate Initiative REKLIM (regional climate change; https://www.reklim.de/en, last access: 31 December 2021). TS was funded by the Swedish Research Council (VR grant 2017-03888). The authors thank N. Tepe for ICP-MS analyses, K. Fitzsimmons for donating loess samples from Central Asia and SE Europe, C. Zdanowicz for Yukon loess samples, and V. Ponomareva and L. Bazanova for samples from Kamchatka. The contribution by C.D.B and P.F. is part of the REDDUST project supported by LEFE/INSU and OSU-EFLUVE. The paper benefited from the constructive and insightful comments of three anonymous reviewers.
Biscaye, P. E. (1965). Mineralogy and sedimentation of recent deep-sea clay in the Atlantic Ocean and adjacent seas and oceans. *The Geological Society of America Bulletin*, **76**(7), 803–832. https://doi.org/10.1130/0016-7606(1965)76[803:MASORD]2.0.CO;2

Biscaye, P. E., Grousset, F. E., Revel, M., Van der Gaast, S., Zielinski, G. A., Vaars, A., & Kukla, G. (1997). Asian provenance of last glacial maximum dust in the GISP-2 ice core, Summit, Greenland. *Journal of Geophysical Research*, 102(C12), 26765–26781. https://doi.org/10.1029/97JC01249

Blichter-Toft, J., Chauvel, C., & Albaréde, F. (1997). Separation of Hf and Lu for high-precision isotope analysis of rock samples by magnetic sector-multiple collector ICP-MS. *Contributions to Mineralogy and Petrology*, 127(3), 248–260. https://doi.org/10.1007/s004100050278

Bond, G., Broecker, W., Johnsen, S., McManus, J., Labeyrie, L., Jouzel, J., & Bonani, G. (1993). Correlations between climate records from North Atlantic sediments and Greenland ice. *Nature*, **365**(6442), 143–147. https://doi.org/10.1038/365143a0

Bory, A. J.-M., Biscaye, P. E., & Grousset, F. E. (2003). Two distinct seasonal Asian source regions for mineral dust deposited in Greenland (NorthGRIP). *Geophysical Research Letters*, 30(4), 1167. https://doi.org/10.1029/2002GL016446

Bory, A. J.-M., Biscaye, P. E., Piotrowski, A. M., & Steffensen, J. P. (2003). Regional variability of ice core dust composition and provenance in Greenland. *Geochemistry, Geophysics, Geosystems*, 4(12), 1107. https://doi.org/10.1029/2003GC000627

Bourne, A. J., Abbott, P. M., Albert, P. G., Cook, E., Pearce, N. J. G., Ponomareva, V., et al. (2016). Underestimated risks of recurrent long-range ash dispersal from northern Pacific Arc volcanoes. *Scientific Reports*, 6(1), 29837. https://doi.org/10.1038/srep29837

Bouvier, A., Vervoort, J. D., & Patchett, P. J. (2008). The Lu–Hf and Sm–Nd isotopic composition of CHUR: Constraints from unequilibrated chondrites and implications for the bulk composition of terrestrial planets. *Earth and Planetary Science Letters*, 273(1–2), 48–57. https://doi.org/10.1016/j.epsl.2008.06.010

Braconnot, P., Harrison, S. P., Kageyama, M., Bartlein, P. J., Masson-Delmotte, V., Abe-Ouchi, A., et al. (2012). Evaluation of climate models using palaeoclimatic data. *Nature Climate Change*, 2(6), 417–424. https://doi.org/10.1038/nclimate1456

Brass, G. W. (1975). The effect of weathering on the distribution of strontium isotopes in weathering profiles. *Geochimica et Cosmochimica Acta*, 39(12), 1647–1653. https://doi.org/10.1016/0016-7037(75)90086-1

Burton, G. R., Rosman, K. J. R., Candelone, J.-P., Burn, L. J., Boutron, C. F., & Hong, S. (2007). The impact of climatic conditions on Pb and Sr isotopic ratios found in Greenland ice, 7–150 ky BP. *Earth and Planetary Science Letters*, 259(3–4), 557–566. https://doi.org/10.1016/j.epsl.2007.05.015

Capo, R. C., Stewart, B. W., & Chadwick, O. A. (1998). Strontium isotopes as tracers of ecosystem processes: Theory and methods. *Geoderma*, 82(1–3), 197–225. https://doi.org/10.1016/S0016-7061(97)00102-X

Capron, E., Rasmussen, S. O., Popp, T. J., Erhardt, T., Fischer, H., Landais, A., et al. (2021). The anatomy of past abrupt warmings recorded in Greenland ice. *Nature Communications*, 12(1), 2106. https://doi.org/10.1038/s41467-021-22241-w

Chen, J., Li, G., Yang, J., Rao, W., Lu, H., Balsam, W., et al. (2007). Nd and Sr isotopic characteristics of Chinese deserts: Implications for the provenances of Asian dust. *Geochimica et Cosmochimica Acta*, 71(15), 3904–3914. https://doi.org/10.1016/j.gca.2007.04.033

Cheng, L. Q., Song, Y. G., Wu, Y., Liu, Y. G., Liu, H. F., Chang, H., et al. (2021). Drivers for asynchronous patterns of dust accumulation in central and eastern Asia and in Greenland during the Last Glacial Maximum. *Geophysical Research Letters*, 48(5), e2020GL091194. https://doi.org/10.1029/2020GL091194

Chin, M., Savoie, D. L., Huebert, B. J., Bandy, A. R., Thornton, D. C., Bates, T. S., et al. (2000). Atmospheric sulfur cycle simulated in the global model GOCART: Comparison with field observations and regional budgets. *Journal of Geophysical Research*, 105(D20), 24689–24712. https://doi.org/10.1029/2000JD900385

Chu, N.-C., Taylor, R. N., Chavagnac, V., Nesbitt, R. W., Boella, R. M., Milton, J. A., et al. (2002). Hf isotope ratio analysis using multi-collector inductively coupled plasma mass spectrometry: An evaluation of isobaric interference corrections. *Journal of Analytical Atomic Spectrometry*, 17(12), 1567–1574. https://doi.org/10.1039/B206707B

Chylek, P., Lesins, G., & Lohmann, U. (2001). Enhancement of dust source area during past glacial periods due to changes of the Hadley circulation. *Journal of Geophysical Research*, 106(D16), 18477–18485. https://doi.org/10.1029/2000JD900583

Clark, P. U., Dyke, A. S., Shakun, J. D., Carlson, A. E., Clark, J., Wohlfarth, B., et al. (2009). The Last Glacial Maximum. *Science*, 325(5941), 710–714. https://doi.org/10.1126/science.1172873

Cline, R. M. L., Hays, J. D., Prell, W. L., Ruddiman, W. F., Moore, T. C., Kipp, N. G., et al. (1984). The last interglacial ocean. *Quaternary Research*, 2, 123–224. https://doi.org/10.1016/0033-5894(84)90098-X

Cook, E., Portnyagin, M., Ponomareva, V., Bazanova, L., Svensson, A., & Garbe-Schönberg, D. (2018). First identification of cryptotephra from the Chumchatka Peninsula in a Greenland ice core: Implications of a widespread marker deposit that links Greenland to the Pacific northwest. *Quaternary Science Reviews*, 181, 200–206. https://doi.org/10.1016/j.quascirev.2017.11.036

Craig, H. (1961). Isotopic variations in meteoric waters. *Science*, 133(3465), 1702–1708. https://doi.org/10.1126/science.133.3465.1702

Dansgaard, W. (1964). Stable isotopes in precipitation. *Tellus*, 16(4), 436–468. https://doi.org/10.1111/j.2153-3490.1964.tb00181.x

Dansgaard, W., Johnsen, S. J., Clausen, H. B., Dahl-Jensen, D., Gundestrup, N. S., Hammer, C. U., et al. (1993). Evidence for general instability of past climate from a 250-kyr ice-core record. *Nature*, 364(6434), 218–220. https://doi.org/10.1038/364218a0

Dasch, E. J. (1969). Strontium isotopes in profiles, deep-sea sediments, and sedimentary rocks. *Geochimica et Cosmochimica Acta*, 33(12), 1521–1552. https://doi.org/10.1016/0016-7037(69)90153-7

de Angelis, M., Steffensen, J. P., Legrand, M., Clausen, H., & Hammer, C. (1997). Primary aerosol (sea salt and soil dust) deposited in Greenland ice during the last climatic cycle: Comparison with East Antarctic records. *Journal of Geophysical Research*, 102(C12), 26681–26698. https://doi.org/10.1029/97JC01298

Delmonte, B., Basile-Doelsch, I., Petit, J.-R., Maggi, V., Revel-Rolland, M., Michard, A., et al. (2004). Comparing the EPICA and Vostok dust records during the last 220, 000 years: Stratigraphical correlation and provenance in glacial periods. *Earth-Science Reviews*, 66(1–2), 63–87. https://doi.org/10.1016/j.earscirev.2003.10.004

Erhardt, T., Capron, E., Rasmussen, S. O., Schöpbach, S., Bigler, M., Adolphi, F., & Fischer, H. (2019). Decadal-scale progression of the onset of Dansgaard–Oeschger warming events. *Climate of the Past*, 15(2), 811–825. https://doi.org/10.5194/cp-15-811-2019

Faure, G. (2001). *Origin of igneous rocks: The isotopic evidence*. Springer-Verlag.

Faure, G., & Mensing, T. M. (2005). *Isotopes: Principles and applications* (3rd ed.). John Wiley & Sons Inc.

Feng, J.-L., Hu, Z.-G., Cui, J.-Y., & Zhu, L.-P. (2010). Distributions of lead isotopes with grain size in aeolian deposits. *Terra Nova*, 22, 257–263. https://doi.org/10.1111/j.1365-3121.2010.00941.x

Feng, J.-L., Zhu, L.-P., Zhen, X.-L., & Hu, Z.-G. (2009). Grain size effect on Sr and Nd isotopic compositions in aeolian dust: Implications for tracing dust provenance and Nd model age. *Geochemical Journal*, 43(2), 123–131. https://doi.org/10.2343/geochemj.1.0007

Fischer, H., Siggaard-Andersen, M.-L., Ruth, U., Röthlisberger, R., & Wolff, E. (2007). Glacial/Interglacial changes in mineral dust and sea salt records in polar ice cores: Sources, transport, deposition. *Reviews of Geophysics*, 45(1), RG1002. https://doi.org/10.1029/2005RG000192
Fu, Y., Hao, Q. Z., Peng, S., Markovic, S. B., Gao, X., Han, L., et al. (2021). Clay mineralogy of the Stari Slankamen (Serbia) loess-paleosol sequence during the last glacial cycle d Implications for dust provenance and interglacial climate. *Quaternary Science Reviews*, **263**, 106990.
https://doi.org/10.1016/j.quascirev.2021.106990

Fuhrer, K., Wolff, E. W., & Johnsen, S. J. (1999). Timescales for dust variability in the Greenland Ice Core Project (GRIP) ice core in the last 100, 000 years. *Journal of Geophysical Research*, **104**(D24), 31043–31052. https://doi.org/10.1029/1999JD0900929

Ginoux, P., Chin, M., Tegen, I., Prospero, J. M., Holben, B., Dubovik, O., & Lin, S.-J. (2001). Sources and distributions of dust aerosols simulated with the GOCART model. *Journal of Geophysical Research*, **106**(D17), 20255–20273. https://doi.org/10.1029/2000JD000053

Griffin, J. J., Windom, H., & Goldberg, E. D. (1968). The distribution of clay minerals in the World ocean. *Deep-Sea Research and Oceanographic Abstracts*, **15**(4), 433–459. https://doi.org/10.1016/0011-7471(68)90051-X

Griffin, W. L., Pearson, N. J., Belousova, E. A., & Saeed, A. (2006). Comment: Hf-Isotope heterogeneity in zircon 91500. *Chemical Geology*, **233**(3–4), 358–363. https://doi.org/10.1016/j.chemgeo.2006.03.007

Grousset, F. E., & Biscaye, P. E. (2005). Tracing dust sources and transport patterns using Sr, Nd and Pb isotopes. *Chemical Geology*, **222**(3–4), 149–167. https://doi.org/10.1016/j.chemgeo.2005.05.006

Han, C. H., Hur, S. D., Han, Y. C., Lee, K. H., Hong, S. M., Erhardt, T., et al. (2018). High-resolution isotopic evidence for a potential Saharan provenance of Greenland glacial dust. *Scientific Reports*, **8**(1), 15582. https://doi.org/10.1038/s41598-018-33859-0

Hansson, M. E. (1994). The Renland ice core. A Northern Hemisphere record of aerosol composition over 120, 000 years. *Tellus*, **46B**(5), 390–418. https://doi.org/10.1034/j.1600-0889.1994.t01-4-00005.x

Hines, K. M., Bromwich, D. H., Bai, L., Bitz, C. M., Powers, J. G., & Manning, K. W. (2015). Sea ice enhancements to Polar WRF. *Monthly Weather Review*, **143**(6), 2363–2385. https://doi.org/10.1175/MWR-D-14-00344.1

Hoskin, P. W. O., & Schaltegger, U. (2003). The composition of zircon and igneous and metamorphic petrogenesis. *Reviews in Mineralogy and Geochemistry*, **53**(1), 27–62. https://doi.org/10.2113/0530027

Huber, C., Leuenberger, M., Spahni, R., Flückiger, J., Schwander, J., Stocker, T. F., et al. (2006). Isotope calibrated Greenland temperature record over Marine Isotope Stage 3 and its relation to CH$_{4}$. *Earth and Planetary Science Letters*, **243**(3–4), 504–519. https://doi.org/10.1016/j.epsl.2006.01.002

Jungclaus, J., Giorgetta, M., Reick, C., Legutke, S., Brovkin, V., Crueger, T., et al. (2012). CMIP5 simulations of the Max Planck Institute for Meteorology (MPI-M) based on the MPIESM-P model: The LGM experiment. *ESGF, WDCC at DKRZ*. MXEPlg. https://doi.org/10.1594/WDCC/CMIP5

Kageyama, M., Harrison, S. P., Kapsch, M.-L., Löfverström, M., Lora, J. M., Mikolajewicz, U., et al. (2021). The PMIP4 Last Glacial Maximum experiments: Preliminary results and comparison with the PMIP3 simulations. *Climate of the Past*, **17**(3), 1065–1089. https://doi.org/10.5194/cp-17-1065-2021

Kang, S. G., Roberts, H. M., Wang, X. L., An, Z. S., & Wang, M. (2015). Mass accumulation rate changes in Chinese loess during MIS 2, and asynchrony with records from Greenland ice cores and North Pacific Ocean sediments during the Last Glacial Maximum. *Aeolian Research*, **19**, 251–258. https://doi.org/10.1016/j.aeolia.2015.05.005

Kapp, P., Pullen, A., Pelletier, J. D., Russell, J., Goodman, P., & Cai, F. (2015). From dust to dust: Quaternary wind erosion of the Mu Us desert and Loess Plateau, China. *Geology*, **43**(9), 835–838. https://doi.org/10.1130/G36724.1

Kindler, P., Guillevic, M., Baumgartner, M., Schwander, J., Landais, A., & Leuenberger, M. (2014). Temperature reconstruction from 10 to 120 kyr b2k from the NGRIP ice core. *Climate of the Past*, **10**(2), 887–902. https://doi.org/10.5194/cp-10-887-2014

Kohfeld, K. E., & Harrison, S. P. (2001). DIRTMAP: The geological record of dust. *Earth-Science Reviews*, **54**(1–3), 81–114. https://doi.org/10.1016/S0012-8252(01)00042-3

Lambert, F., Tagliablue, A., Shaffer, G., Lamy, F., Winckler, G., Farias, L., et al. (2015). Dust fluxes and iron fertilization in Holocene and Last Glacial Maximum climates. *Geophysical Research Letters*, **42**(14), 6014–6023. https://doi.org/10.1002/2015GL064250

Lenton, T. M., Rockström, J., Gaffney, O., Rahmstorf, S., Richardson, K., Steffen, W., & Schellnhuber, H. J. (2019). Climate tipping points—Too risky to bet against. *Nature*, **575**(7784), 592–595. https://doi.org/10.1038/d41586-019-03595-0

Li, C., & Born, A. (2019). Coupled atmosphere-ice-ocean dynamics in Dansgaard-Oeschger events. *Quaternary Science Reviews*, **203**, 1–20. https://doi.org/10.1016/j.quascirev.2018.10.031

Licht, A., Pullen, A., Kapp, P., Abell, J., & Giesler, D. (2016). Eolian cannibalism: Reworked loess and fluvial sediment as the main sources of the Chinese Loess Plateau. *The Geological Society of America Bulletin*, **128**(5–6), 944–956. https://doi.org/10.1130/B31375.1

Löfverström, M. (2020). A dynamic link between high-intensity precipitation events in southwestern North America and Europe at the Last Glacial Maximum. *Earth and Planetary Science Letters*, **534**, 116081. https://doi.org/10.1016/j.epsl.2020.116081

Löfverström, M., Caballero, R., Nilsson, J., & Messori, G. (2016). Stationary wave reflection as a mechanism for zonalising Atlantic winter jet at the LGM. *Journal of the Atmospheric Sciences*, **73**(8), 3329–3342. https://doi.org/10.1175/JAS-D-15-0295.1

Ludwig, P., Gavrilov, M. B., Markovic, S. B., Újvári, G., & Lehmkuhl, F. (2021). Simulated regional dust cycle in the Carpathian Basin and the Adriatic Sea region during the Last Glacial Maximum. *Quaternary International*, **581–582**, 114–127. https://doi.org/10.1016/j.quaint.2020.09.048

Ludwig, P., Gomez-Navarro, J. J., Pinto, J. G., Raible, C. C., Wagner, S., & Zorita, E. (2019). Perspectives of regional paleoclimate modeling. *Annals of the New York Academy of Sciences*, **1436**(1), 54–69. https://doi.org/10.1111/nyas.13865

Ludwig, P., Pinto, J. G., Raible, C. C., & Shao, Y. (2017). Impacts of surface boundary conditions on regional climate model simulations of European climate during the Last Glacial Maximum. *Geophysical Research Letters*, **44**(10), 5086–5095. https://doi.org/10.1002/2017GL073622

Ludwig, P., Schaffernicht, E. J., Shao, Y., & Pinto, J. G. (2016). Regional atmospheric circulation over Europe during the Last Glacial Maximum and its links to precipitation. *Journal of Geophysical Research: Atmospheres*, **121**, 2130–2145. https://doi.org/10.1002/2015JD024444

Lugmair, G. W., & Carlson, R. W. (1978). The Sm–Nd history of KREEP. *Proceedings of the Lunar and Planetary Science Conference* 9th, **1**, 689–704.

Lunt, D. J., & Valdes, P. J. (2001). Dust transport to Done C, Antarctica, at the last glacial maximum and present day. *Geophysical Research Letters*, **28**(2), 295–298. https://doi.org/10.1029/2000GL012170

Lupker, M., Aciego, S. M., Bourdon, B., Schwander, J., & Stocker, T. F. (2010). Isotopic tracing (Sr, Nd, U and Hf) of continental and marine aerosols in an 18th century section of the Dye-3 ice core (Greenland). *Earth and Planetary Science Letters*, **295**(1–2), 277–286. https://doi.org/10.1016/j.epsl.2010.04.010

Lynch-Stieglitz, J., Adkins, J. F., Curry, W. B., Dokken, T., Hall, I. R., Herguera, J. C., et al. (2007). Atlantic Meridional overturning circulation during the Last Glacial Maximum. *Science*, **316**(5821), 66–69. https://doi.org/10.1126/science.1137127

Mahowald, N. (2011). Aerosol indirect effect on biogeochemical cycles and climate. *Science*, **334**(6057), 794–796. https://doi.org/10.1126/science.1207374
Mahowald, N., Albani, S., Engelstaedter, S., Winckler, G., & Goman, M. (2011). Model insight into glacial-interglacial paleodust records. *Quaternary Science Reviews*, **30**(7–8), 832–854. https://doi.org/10.1016/j.quascirev.2010.09.007

Mahowald, N., Kohfeld, K., Hansson, M., Balkanski, Y., Harrison, S. P., Prentice, I. C., et al. (1999). Dust sources and deposition during the Last Glacial Maximum and current climate: A comparison of model results with paleodata from ice cores and marine sediments. *Journal of Geophysical Research*, **104**(D13), 15895–15916. https://doi.org/10.1029/1999JD090084

Mahowald, N. M., Muhs, D. R., Levis, S., Rasch, P. J., Yoshioka, M., Zender, C. S., & Luo, C. (2006). Change in atmospheric mineral aerosols in response to climate: Last glacial period, preindustrial, modern, and doubled carbon dioxide climates. *Journal of Geophysical Research*, **111**, D10202. https://doi.org/10.1029/2005JD006653

Martinez-Lamas, R., Toucanne, S., Debret, M., Riboulot, V., Deloffre, J., Boissier, A., et al. (2020). Linking Danube River activity to Alpine Ice-Sheet fluctuations during the last glacial (ca. 33-17 ka BP): Insights into the continental signature of Heinrich Stadials. *Quaternary Science Reviews*, **229**, 106136. https://doi.org/10.1016/j.quascirev.2019.106136

Mayewski, P. A., Meeker, L. D., Whitlow, S., Twickler, M. S., Morrison, M. C., Bloomfield, P., et al. (1994). Changes in atmospheric circulation and ocean ice cover over the North Atlantic during the last 41, 000 years. *Science*, **263**(5154), 1747–1751. https://doi.org/10.1126/science.263.5154.1747

Meyer, I., Davies, G. R., & Stuut, J.-B. W. (2011). Grain size control on Sr-Nd isotope provenance studies and impact on paleoclimate reconstructions: An example from deep sea sediments offshore NW Africa. *Geochemistry, Geophysics, Geosystems*, **12**(3), Q03005. https://doi.org/10.1029/2010GC003355

Meyer, K. W., Aciego, S., & Koorneef, J. M. (2017). Radiogenic isotopic compositions of low concentration dust and aerosol from the GISP2 ice core. *Chemical Geology*, **472**, 31–43. https://doi.org/10.1016/j.chemgeo.2017.09.021

Moine, O., Antoine, P., Hatté, C., Landais, A., Mathieu, J., Prud'homme, C., & Rousseau, D.-D. (2017). The impact of last glacial climate variability in west-European loess revealed by radiocarbon dating of fossil earthworm granules. *Proceedings of the National Academy of Sciences USA*, **114**(24), 6209–6214. https://doi.org/10.1073/pnas.1614751114

Moore, D. M., & Reynolds, R. C. (1997). X-ray diffraction and the identification and analysis of clay minerals (2nd ed.). Oxford University Press.

Patchett, P. J. (1983). Importance of the Lu-Hf isotopic system in studies of planetary chronology and chemical evolution. *Geochimica et Cosmochimica Acta*, **47**(1), 81–91. https://doi.org/10.1016/0016-7037(83)90092-3

Pinto, J. G., & Ludwig, P. (2020). Extratropical cyclones over the North Atlantic and Western Europe during the Last Glacial Maximum and implications for proxy interpretation. *Climate of the Past*, **16**(2), 611–626. https://doi.org/10.5194/cp-16-611-2020

Prömmel, K., Cubasch, U., & Kaspar, F. (2013). A regional climate model study of the impact of tectonic and orbital forcing on African precipitation and vegetation. *Palaeogeography, Palaeoclimatology, Palaeoecology*, **369**, 154–162. https://doi.org/10.1016/j.palaeo.2012.10.015

Rasmussen, S. O., Bigler, M., Blockley, S. P., Blunier, T., Buchardt, S. L., Clausen, H. B., et al. (2014). A stratigraphic framework for abrupt climatic changes during the last glacial period based on three synchronized Greenland ice-core records: Refining and extending the INTI-MATE event stratigraphy. *Quaternary Science Reviews*, **106**, 14–28. https://doi.org/10.1016/j.quascirev.2014.09.007

Reader, M. C., Fung, I., & McFarlane, N. (1999). The mineral dust aerosol cycle during the Last Glacial Maximum. *Journal of Geophysical Research*, **104**(D8), 9381–9398. https://doi.org/10.1029/1999JD900033

Rousseau, D.-D., Antoine, P., & Sun, Y. (2021). How dusty was the last glacial maximum over Europe? *Quaternary Science Reviews*, **254**, 106775. https://doi.org/10.1016/j.quascirev.2020.106775

Ruth, U., Bigler, M., Röthlisberger, R., Siggaard-Andersen, M.-L., Kipfstuhl, S., Goto-Azuma, K., et al. (2007). Ice core evidence for a very tight link between North Atlantic and east Asian glacial climate. *Geophysical Research Letters*, **34**(3), L03706. https://doi.org/10.1029/2006GL027876

Ruth, U., Wagenbach, D., Steffensen, J. P., & Bigler, M. (2003). Continuous record of microparticle concentration and size distribution in the central Greenland NGRIP ice core during the last glacial period. *Journal of Geophysical Research*, **108**, 4098. https://doi.org/10.1029/2002JD002376

Savin, S. M., & Epstein, S. (1970). The oxygen and hydrogen isotope geochemistry of clay minerals. *Geochimica et Cosmochimica Acta*, **34**(1), 25–42. https://doi.org/10.1016/0016-7037(70)90149-3

Savin, S. M., & Hsieh, J. C. C. (1998). The hydrogen and oxygen isotope geochemistry of pedogenic clay minerals: Principles and theoretical background. *Geoderma*, **82**(1–3), 227–253. https://doi.org/10.1016/S0016-7061(97)00103-1

Schaffernicht, E. J., Ludwig, P., & Shao, Y. (2020). Linkage between dust cycle and loess of the Last Glacial Maximum in Europe. *Atmospheric Chemistry and Physics*, **20**(8), 4969–4986. https://doi.org/10.5194/acp-20-4969-2020

Scheuvens, D., Schütz, L., Kandler, K., Ebert, M., & Weinbruch, S. (2013). Bulk composition of northern African dust and its source sediments—A compilation. *Earth-Science Reviews*, **116**, 170–194. https://doi.org/10.1016/j.earscirev.2012.08.005

Schüpbach, S., Fischer, H., Bigler, M., Erhadt, T., Gfeller, G., Leuenberger, D., et al. (2018). Greenland records of aerosol source and atmospheric lifetime changes from the Eemian to the Holocene. *Nature Communications*, **9**(1), 1476. https://doi.org/10.1038/s41467-018-03924-3

Seguinot, J., Ivy-Ochs, S., Jouvet, G., Huss, M., Funk, M., & Preusser, F. (2018). Modelling last glacial cycle ice dynamics in the Alps. *The Cryosphere*, **12**(10), 3265–3285. https://doi.org/10.5194/tc-12-3265-2018

Serno, S., Winckler, G., Anderson, R. F., Maier, E., Ren, H., Gersonde, R., & Haug, G. H. (2015). Comparing dust flux records from the Subarctic North Pacific and Greenland: Implications for atmospheric transport to Greenland and for the application of dust as a chronostratigraphic tool. *Paleoceanography*, **30**(6), 583–600. https://doi.org/10.1002/2014PA002748

Shao, Y. (2004). Simplification of a dust emission scheme and comparison with data. *Journal of Geophysical Research*, **109**, D10202. https://doi.org/10.1029/2003JD004372

Sharp, Z. D., Atudorei, V., & Durakiewicz, T. (2001). A rapid method for determining the hydrogen and oxygen isotope ratios from water and solid hydrous substances. *Chemical Geology*, **178**(1–4), 197–210. https://doi.org/10.1016/S0009-2541(01)00262-5

Sheppard, S. M. F., & Gilg, H. A. (1996). Stable isotope geochemistry of clay minerals. *Clay Minerals*, **31**, 1–24. https://doi.org/10.1180/claymin.1996.031.1.01

Simonsen, M. F., Baccolo, G., Blunier, T., Borunda, A., Delmonte, B., Frei, R., et al. (2019). East Greenland ice core dust record reveals timing of Greenland ice sheet advance and retreat. *Nature Communications*, **10**(1), 4494. https://doi.org/10.1038/s41467-019-12546-2

Singer, A., Dultz, S., & Argamann, E. (2004). Properties of the non-soluble fractions of suspended dust over the Dead Sea. *Atmospheric Environment*, **38**(12), 1745–1753. https://doi.org/10.1016/j.atmosenv.2003.12.026

Skamarock, W. C., Klemp, J. B., Dudhia, J., Gill, D. O., Liu, Z. Q., Berner, J., et al. (2019). A description of the advanced research WRF version 4. In *NCAR technical notes (No. NCAR/TN-556+STR)* (pp. 1–148). https://doi.org/10.5065/1dfh-6p97

Sokolik, I. N., & Toon, O. B. (1996). Direct radiative forcing by anthropogenic airborne mineral aerosols. *Nature*, **381**(6584), 681–683. https://doi.org/10.1038/381681a0

Sprenger, M., & Wernli, H. (2015). The LAGRANTO Lagrangian analysis tool – Version 2.0. *Geoscientific Model Development*, **8**, 2569–2586. https://doi.org/10.5194/gmd-8-2569-2015
Steffensen, J. P. (1997). The size distribution of microparticles from selected segments of the Greenland Ice Core Project ice core representing different climatic periods. *Journal of Geophysical Research*, 102(C12), 26755–26763. https://doi.org/10.1029/97JC01490

Stevens, T., Buylaert, J.-P., Lu, H., Thiel, C., Murray, A., Frechen, M., et al. (2016). Mass accumulation rate and monsoon records from Xifeng, Chinese Loess Plateau, based on a luminescence age model. *Journal of Quaternary Science*, 31(4), 391–405. https://doi.org/10.1002/jqs.2848

Stevens, T., Buylaert, J. P., Thiel, C., Újvári, G., Yi, S., Murray, A. S., et al. (2018). Ice-volume-forced erosion of the Chinese Loess Plateau global Quaternary stratotype site. *Nature Communications*, 9(1), 983. https://doi.org/10.1038/s41467-018-03329-2

Svensson, A., Biscaye, P. E., & Grousset, F. E. (2000). Characterization of late glacial continental dust in the Greenland Ice Core Project ice core. *Journal of Geophysical Research*, 105(D4), 4637–4656. https://doi.org/10.1029/1999JD0910193

Takemura, T., Egashira, M., Matsuzawa, K., Ichijo, H., O'ishi, R., & Abe-Ouchi, A. (2009). A simulation of the global distribution and radiative forcing of soil dust aerosols at the Last Glacial Maximum. *Atmospheric Chemistry and Physics*, 9, 3061–3073. https://doi.org/10.5194/acp-9-3061-2009

Tegen, I., Lacis, A. A., & Fung, I. (1996). The influence on climate forcing of mineral aerosols from disturbed soils. *Nature*, 380(6573), 419–422. https://doi.org/10.1038/380419a0

Tjallingii, R., Claussen, M., Stuut, J. B., Fohlmeister, J., Jahn, A., Bickert, T., et al. (2008). Coherent high- and low-latitude control of the north-west African hydrological balance. *Nature Geoscience*, 1(10), 670–675. https://doi.org/10.1038/ngeo289

Újvári, G. (2021). Clay mineralogy, Sr-Nd-Hf-H isotope data of NGRIP ice core dust and potential source area samples with mixing calculations and Monte Carlo simulations [Dataset]. PANGAEA. https://doi.org/10.1594/PANGAEA.936106

Újvári, G., & Klötzli, U. (2015). U-Pb ages and Hf isotopic composition of zircons in Austrian last glacial loess: Constraints on heavy mineral sources and sediment transport pathways. *International Journal of Earth Sciences*, 104(5), 1365–1385. https://doi.org/10.1007/s00531-014-1139-x

Újvári, G., Klötzli, U., Horschinegg, M., Wegner, W., Hippler, D., Kiss, G. I., & Palcsu, L. (2021). Rapid decomposition of geological samples by ammonium bifluoride (NH₄HF₂) for combined Hf-Nd-Sr isotopic analyses. *Rapid Communications in Mass Spectrometry*, 35(11), e9081. https://doi.org/10.1002/rcm.9081

Újvári, G., Kovács, J., Varga, G., Raucsik, B., & Marković, S. B. (2010). Dust flux estimates for the Last Glacial Period in East Central Europe based on terrestrial records of loess deposits: A review. *Quaternary Science Reviews*, 29(23–24), 3157–3166. https://doi.org/10.1016/j.quascirev.2010.07.005

Újvári, G., Stevens, T., Molnár, M., Demény, A., Lambert, F., Varga, G., et al. (2017). Coupled European and Greenland last glacial dust activity driven by North Atlantic climate. *Proceedings of the National Academy of Sciences USA*, 114(50), 10632–10638. https://doi.org/10.1073/pnas.1712651114

Újvári, G., Stevens, T., Svensson, A., Klötzli, U. S., Manning, C., Németh, T., et al. (2015). Two possible source regions for central Greenland last glacial dust. *Geophysical Research Letters*, 42(23), 10399–10408. https://doi.org/10.1002/2015GL066153

Újvári, G., Wegner, W., Klötzli, U., Horschinegg, M., & Hippler, D. (2018). Sr-Nd-Hf isotopic analysis of <10 mg dust samples: Implications for ice core dust source fingerprinting. *Geochemistry, Geophysics, Geosystems*, 19(1), 60–72. https://doi.org/10.1002/2017GC007136

van de Flierdt, T., Goldstein, S. L., Hemming, S. R., Roy, M., Frank, M., & Halliday, A. N. (2007). Global neodymium-hafnium isotope systematics—Revisited. *Earth and Planetary Science Letters*, 259(3–4), 432–441. https://doi.org/10.1016/j.epsl.2007.05.003

Viers, J., & Wasserburg, G. J. (2004). Behavior of Sm and Nd in a lateritic soil profile. *Geochimica et Cosmochimica Acta*, 68(9), 2043–2054. https://doi.org/10.1016/j.gca.2003.10.034

Werner, M., Tegen, I., Harrison, S. P., Kohfeld, K. E., Prentice, I. C., Balkanski, Y., et al. (2002). Seasonal and interannual variability of the mineral dust cycle under present and glacial climate conditions. *Journal of Geophysical Research*, 107, 4744. https://doi.org/10.1029/2002JD002365

Yokoo, Y., Nakano, T., Nishikawa, M., & Quan, H. (2004). Mineralogical variation of Sr-Nd isotopic and elemental compositions in loess and desert sand from the central Loess Plateau in China as a provenance tracer of wet and dry deposition in the northwestern Pacific. *Chemical Geology*, 204(1–2), 45–62. https://doi.org/10.1016/j.chemgeo.2003.11.004

Zdanowicz, C., Hall, G., Vaive, J., Amelin, Y., Percival, J., Girard, I., et al. (2006). Asian dustfall in the St. Elias Mountains, Yukon, Canada. *Geochimica et Cosmochimica Acta*, 70(14), 3493–3507. https://doi.org/10.1016/j.gca.2006.05.005

Zhao, W., Balsam, W., Williams, E., Long, X., & Ji, J. (2018). Sr-Nd-Hf isotopic fingerprinting of transatlantic dust derived from North Africa. *Earth and Planetary Science Letters*, 486, 23–31. https://doi.org/10.1016/j.epsl.2018.01.004

Zhao, W., Sun, Y., Balsam, W., Zeng, L., Lu, H., Otgonbayar, K., & Ji, J. (2015). Clay-sized Hf-Nd-Sr isotopic composition of Mongolian dust as a fingerprint for regional to hemispherical transport. *Geophysical Research Letters*, 42(13), 5661–5669. https://doi.org/10.1002/2015GL064357

## References From the Supporting Information

Bosq, M., Bertran, P., Degeai, J.-P., Queffelec, A., & Moine, O. (2020). Geochemical signature of sources, recycling and weathering in the Last Glacial loess from the Rhône Valley (southeast France) and comparison with other European regions. *Aeolian Research*, 42, 100561. https://doi.org/10.1016/j.aeolia.2019.100561

Chang, Q., Mishima, T., Yabuki, S., Takahashi, Y., & Shimizo, H. (2000). Sr and Nd isotope ratios and REE abundances of moraines in the mountain areas surrounding the Taklimakan Desert, NW China. *Geochemical Journal*, 34(6), 407–427. https://doi.org/10.2343/geochemj.34.407

Chauvel, C., Garcon, M., Bureau, S., Besnault, A., Jahn, B.-M., & Ding, Z. L. (2014). Constraints from loess on the Hf-Nd isotopic composition of the upper continental crust. *Earth and Planetary Science Letters*, 388, 48–58. https://doi.org/10.1016/j.epsl.2013.11.045

Chen, S.-H., & Sun, W.-Y. (2002). A one-dimensional time dependent cloud model. *Journal of the Meteorological Society of Japan*, 80(1), 99–118. https://doi.org/10.2151/jmsj.80.99

De Maesschalck, R., Jouan-Rimbaud, D., & Massart, D. L. (2000). The Mahalanobis distance. *Chemometrics and Intelligent Laboratory Systems*, 50, 1–18. https://doi.org/10.1016/S0169-7439(99)00047-7

Ferrat, M., Weiss, D., Strekopytov, S., Dong, S., Chen, H., Najorka, J., et al. (2011). Improved provenance tracing of Asian dust sources using rare Earth elements and selected trace elements for palaeomonsoon studies on the eastern Tibetan Plateau. *Geochimica et Cosmochimica Acta*, 75(21), 6374–6399. https://doi.org/10.1016/j.gca.2011.08.025

George, R., Turner, S., Hawkesworth, C., Morris, J., Nie, C., Ryan, J., & Zheng, S.-H. (2003). Melting processes and fluid and sediment transport rates along the Alaska-Aleutian arc from an integrated U-Th-Ra-Be isotope study. *Journal of Geophysical Research*, 108(B5), 2252. https://doi.org/10.1029/2002JB001916

Hao, Q. Z., Guo, Z. T., Qiao, Y. S., Xu, B., & Oldfield, F. (2010). Geochemical evidence for the provenance of middle Pleistocene loess deposits in southern China. *Quaternary Science Reviews*, 29(23–24), 3317–3326. https://doi.org/10.1016/j.quascirev.2010.08.004
Honda, M., Yabuki, S., & Shimizu, H. (2004). Geochemical and isotopic studies of aeolian sediments in China. Sedimentology, 51(2), 211-230.
https://doi.org/10.1111/j.1365-3091.2004.00618.x

Hong, S.-Y., Noh, Y., & Dudhia, J. (2006). A new vertical diffusion package with an explicit treatment of entrainment processes. Monthly Weather Review, 134(9), 2318-2341. https://doi.org/10.1175/MWR3199.1

Iacono, M. J., Delamere, J. S., Mlawer, E. J., Shephard, M. W., Clough, S. A., & Collins, W. D. (2008). Radiative forcing by long-lived greenhouse gases: Calculations with the AER radiative transfer models. Journal of Geophysical Research, 113, D13103. https://doi.org/10.1029/2008JD009944

Jeong, G. Y. (2020). Mineralogy and geochemistry of Asian dust: Dependence on migration path, fractionation, and reactions with polluted air. Atmospheric Chemistry and Physics, 20(12), 7411-7428. https://doi.org/10.5194/acp-20-7411-2020

Jimenez, P. A., Dudhia, J., Gonzalez-Rouco, J. F., Navarro, J., Montávez, J. P., & García-Bustamante, E. (2012). A revised scheme for the WRF surface layer formulation. Monthly Weather Review, 140(3), 898-918. https://doi.org/10.1175/MWR-D-11-00056.1

Kain, J. S. (2004). The Kain-Fritsch convective parameterization: An update. Journal of Applied Meteorology and Climatology, 43(1), 170-181.
https://doi.org/10.1175/1520-0450(2004)043<0170:TKCPAU>2.0.CO;2

Muhs, D. R., Budahn, J. R., Prospero, J. M., & Carey, S. N. (2007). Geochemical evidence for African dust inputs to soils of western Atlantic islands: Barbados, Bahamas, and Florida. Journal of Geophysical Research, 112, F02009. https://doi.org/10.1029/2005JF000445

Pourmand, A., Prospero, J. M., & Sharifi, A. (2014). Geochemical fingerprinting of trans-Atlantic African dust based on radiogenic Sr-Nd-Hf isotopes and rare Earth element anomalies. Geology, 42(8), 675-678. https://doi.org/10.1130/G35624.1

Tewari, M., Chen, F., Wang, W., Dudhia, J., Lemone, M. A., Mitchell, K. A., et al. (2004). Implementation and verification of the unified NOAH land surface model in the WRF model. In 20th conference on weather analysis and forecasting/16th conference on numerical weather prediction (pp. 11-15). American Meteorological Society. Retrieved from http://n2t.net/ark:/85065/d7fb523p

Yang, S., Ding, F., & Ding, Z. (2006). Pleistocene chemical weathering history of Asian arid and semi-arid regions recorded in loess deposits of China and Tajikistan. Geochimica et Cosmochimica Acta, 70(7), 1695-1709. https://doi.org/10.1016/j.gca.2005.12.012

Yang, X. P., Zhu, B. Q., & White, P. D. (2007). Provenance of aeolian sediment in the Taklamakan Desert of western China, inferred from REE and major-elemental data. Quaternary International, 175(1), 71-85. https://doi.org/10.1016/j.quaint.2007.03.005