3.1 Dust cycle hypothesis

In line with previous modelling (COHMAP Members, 1988; Ludwig et al., 2016) and fieldwork studies (Dietrich and Seelos, 2010; Krauß et al., 2016; Römer et al., 2016), we hypothesize that east sector winds (i.e. northeasters, easterlies, and southeasters) dominated the mineral dust cycle over central Europe during the LGM (Fig. 2). This hypothesis also implies a linkage of dust sources in central and eastern Europe during the LGM and the loess deposits in Europe. It is suggested here that a greater proportion of all LGM dust deposits in central and eastern Europe comes more from sources in central and eastern Europe than from sources in the English Channel. The east sector winds likely contributed substantially to the formation of the European loess belt in central Europe. Among them, the northeasters and easterlies originated most likely from dry winds that flowed down the slopes of the southern and eastern EIS margins, where they picked up and turned gradually into northeasters and easterlies. By blowing over the bare proglacial EIS areas, they generated dust emissions and carried the dust westwards, implying dust depositions in areas west of the respective dust sources.

3.2 East sector winds and cyclones over central Europe

In agreement with this hypothesis, glacial simulations for 90 kyr ago evidenced katabatic winds over the EIS (Krinner et al., 2004), and global climate model (GCM) simulations for the LGM indicate prevailing east sector winds over central and eastern Europe (COHMAP Members, 1988; Ludwig et al., 2016). In Germany, several aeolian sediment records that are dated to the LGM originated from more eastern sources (Dietrich and Seelos, 2010; Krauß et al., 2016; Römer et al., 2016). The CWT frequencies for the present (not shown) and the PI era are very similar; therefore it is possible to use the term present day to refer to both the PI and the actual present-day frequencies. In contrast to the dominant present-day anticyclones and west sector winds (southwesters, westerlies, and northwesters), east sector winds (36 %) and cyclones (22 %) prevailed over central Europe during the LGM (Table 2). The east sector winds are associated with a strong EIS high (Fig. 2a and COHMAP Members, 1988). The increased frequency of cyclones over central Europe is consistent with the analysis of the LGM storm tracks, which deviated from their present-day course (Hofer et al., 2012), running either along central Europe, the Mediterranean, or the Nordic Seas (Florineth and Schlüchter, 2000; Luetscher et al., 2015; Ludwig et al., 2016). Their Mediterranean course is consistent with the Alpine, western, and southern European climate proxies (Luetscher et al., 2015). In addition, the proxies indicate a storm track branch split-off over the Adriatic that ran past the Eastern Alps to central Europe (Florineth and Schlüchter, 2000; Luetscher et al., 2015; Újvári et al., 2017). These proxy-based findings are in line with the more frequent cyclones in central Europe during the LGM (Table 2). This, in turn, can be related to the stronger and southward-shifted jet stream (Luetscher et al., 2015; Ludwig et al., 2016) and the missing Scandinavian cyclone tracks, which were deflected southwards by the blocking EIS high. As a result, their frequency increased over central Europe (Table 2), consistent with susceptibility- and grain-size-based results that suggest more frequent storms over western Europe. The east sector winds, which more than doubled in frequency in comparison to today (36 % compared to 17 %, Table 2), need to be incorporated to establish a more complete understanding of the main drivers of the dust cycle in Europe during the LGM (Fig. 9a). These winds are also evidenced by northern-central European grain size records for the Late Pleniglacial (Bokhorst et al., 2011). Sediment layers attributed to east wind dated to 36–18 kyr BP are abundant in the Dehner Maar sediments (Eifel, Germany, 50.3^∘ N, 6.5^∘ E; Dietrich and Seelos, 2010). Their provenance showed that up to every fifth dust storm over the Eifel came from the east (Dietrich and Seelos, 2010).

Table 2Circulation weather type occurrence frequencies (%) for central Europe (centred at 17.5^∘ E and 47.5^∘ N) during the LGM and the pre-industrial period (PI). The frequencies are based on the LGM and the PI simulation of the Max Planck Institute Earth System Model. The circulation weather type classes are cyclonic (C), anticyclonic (A), northeasterly (NE), and easterly (E) followed by the remaining standard wind directions.

Download Print Version | Download XLSX

Table 3Seasonal CWT occurrence frequencies (%) for central Europe (centred at 17.5^∘ E and 47.5^∘ N) during the LGM. The frequencies are based on MPI-LGM simulation. The CWT classes are cyclonic (C), anticyclonic (A), northeasterly (NE), and easterly (E) followed by the remaining standard wind directions. Sum E is the sum of the east sector winds (NE, E, SE). The seasons are labelled DJF (winter), MAM (spring), JJA (summer), and SON (autumn).

Download Print Version | Download XLSX

Our findings are in agreement with fieldwork-based results of Römer et al. (2016), who found evidence for strong east sector winds over northern, central, and western Germany for 23 to 20 kyr ago. Also, loess in the Harz Foreland indicates a shift to prevailing east sector winds for the LGM (Krauß et al., 2016). The location of aeolian ridges along rivers in northeastern Belgium and a core transect near Leuven also support our finding by evidencing northeasters for the Late Pleniglacial (Renssen et al., 2007). In addition, northerlies, northeasters, and easterlies were inferred from loess deposits west of the Maas (Renssen et al., 2007). Also, for Denmark, wind-polished boulders evidence dominant easterlies and southeasters in the period of 22 to 17 kyr ago (Renssen et al., 2007). The CWT frequency distribution for the LGM (Table 2) contradicts the finding (Renssen et al., 2007) of prevailing west sector winds during the LGM in central Europe (40–55^∘ N, 0–30^∘ E). The distribution also contrasts with the finding (Sima et al., 2013) of prevailing winds from the west-northwest in eastern central Europe, in particular for the area around Stayky (50^∘ N, 31^∘ E). More precisely, the CWT-W and CWT-NW regimes occurred in eastern central Europe in total less than 10 % of the times during the LGM (Table 2), which is even less than the expectation value for a single weather type in the case of a uniform CWT frequency distribution. Conversely, the significant role of the east sector winds (Table 2) is consistent with the deposits on the west bank of the Dnieper (Sima et al., 2013), which are also the loess deposits closest to Stayky. In addition, sandy soil texture and sand dunes indicate prevailing northerlies and northeasters over Dobrudja (44.32^∘ N, 28.18^∘ E), the eastern Wallachian Plain (both located in Romania), and Stary Kaydaky (Ukraine, 48.37^∘ N, 35.12^∘ E; Buggle et al., 2008). The northerlies over Ukraine originated from katabatic winds descending from the EIS (Buggle et al., 2008). The high aridity and grain size variations in the Surduk (Serbia, Table S1) and Stari Bezradychy (Ukraine, Table S1) records evidence prevailing dry and periodically strong east sector winds (Antoine et al., 2009a; Bokhorst et al., 2011).

3.3 Dust emissions from the Eurasian ice sheet margin

The model-simulated dust emission (Fig. 3) indicates that most dust in Europe was emitted from the less elevated corridor between the Alps, the Black Sea, and the EIS (45–55^∘ N). This finding is consistent with loess-based dust flux estimates (Újvári et al., 2010). The highest emission rates (>10⁵ $\mathrm{g}{\mathrm{m}}^{-\mathrm{2}}{\mathrm{yr}}^{-\mathrm{1}}$) occurred along the southern EIS margin (51–53^∘ N, 15–18^∘ E, Fig. 3). This location is in line with the location of the highest emissions found in the Greenland stadial GCM simulation of Sima et al. (2013); yet our simulation indicates a larger upper limit for the emission rates (1000 $\mathrm{g}{\mathrm{m}}^{-\mathrm{2}}{\mathrm{yr}}^{-\mathrm{1}}$). Our results also show high emissions in the dried-up English Channel and the German Bight (Fig. 3). For the latter, they compare well with the average emission of 140 and the maximum emission greater than 200 $\mathrm{g}{\mathrm{m}}^{-\mathrm{2}}{\mathrm{yr}}^{-\mathrm{1}}$ based on a glacial climate simulation (Sima et al., 2009).

The loess deposits (Újvári et al., 2010) and the model results are consistent in that the Carpathian Basin was both a dust source and a dust sink (Figs. 3 and 4). Major dust sources surrounding the Carpathians and the Eastern Alps (Fig. 3) are in line with deposits in Serbia and the Carpathian Basin (Újvári et al., 2010; Bokhorst et al., 2011). The dust emissions from the Lower Danube Basin (Fig. 3) are in agreement with plentiful sediment supply, strong winds, and dry conditions inferred from the plateau loess in Urluia, located near the Black Sea in southeastern Romania (Fitzsimmons and Hambach, 2014). Also, the emissions from the western Black Sea littoral (Fig. 3) are consistent with provenance analyses of eastern Dobrogea loess in the Lower Danube Basin (Jipa, 2014). Our results indicate a close relationship between strong dust emissions and low terrains (or basins). This relationship is found for the North Sea Basin and the European plains bordering the EIS, the Caucasus, the Carpathians, or the Massif Central (Figs. 1 and 3). The dust emissions from the EIS margin and from the foothills of the European mountains (Fig. 3) are consistent with the loess-based finding of significant aeolian dust contributions from glaciogenic and orogenic dust sources (Újvári et al., 2010).

3.4 Conforming dust deposition and loess accumulation rates

Compared with the GCMs (Werner, 2002; Mahowald et al., 2006; Hopcroft et al., 2015; Sudarchikova et al., 2015; Albani et al., 2016), the WRF-Chem-LGM dust deposition rates (F_D, Fig. 4) reproduce the MARs (Table S1, Fig. 4a and b) and MAR10 (Table S1, Fig. 4c and d) better, at least by 1 order of magnitude. One factor for this improvement is most likely the higher spatio-temporal resolution (Ludwig et al., 2019) of the WRF-Chem-LGM experiments combined with the provided more highly resolved geographical input data, for example the regional LGM topography, land use, and dynamic (yet monthly prescribed) vegetation cover. The boundary conditions provided by the MPI-LGM could also be a factor for this improvement. Taking into account that the MPI-ESM experiment for the present reproduces the observed atmospheric circulation over Europe better than other GCMs (Ludwig et al., 2016), it is likely that MPI-LGM also reproduces the LGM conditions more realistically. Another factor could be the orography-based estimated fraction of alluvium (Ginoux et al., 2001) combined with the proxy-based reconstructed bare soil fraction (Cline et al., 1984) to calculate the spatial erodibility distribution. Based on this distribution, the WRF-Chem-LGM was able to suppress unrealistic numerical dust emission from areas with low or zero erodibility. Most likely, the improvement also results from selecting the well-tested and observation-confirmed Shao dust emission scheme (Shao, 2004; Kang et al., 2011). For example, this scheme takes into account the dynamic moisture changes at the soil surface. Due to our recent improvement of the Shao dust emission scheme, the effect of snow cover on dust emission has also been taken into account in the WRF-Chem-LGM experiments.

Figure 4Dust deposition rates for the Last Glacial Maximum, comprising particles of up to 20 µm in diameter (F_D20) using (a) dynamic downscaling (F_D20DD) and (b) statistic dynamic downscaling (F_D20SD). Panels (c) and (d) are as (a) and (b), but for particles up to 12 µm (F_D12). Each blue circle size represents one mass accumulation rate (MAR, Table S1 column 5) magnitude. Each red circle size represents one reduced mass accumulation rate (MAR10, Table S1 column 6) magnitude. MAR and MAR10 values compiled in Table S1. The simulation-based (F_D20, F_D12) and the fieldwork-based (MAR, MAR10) rates result from independent data. Delineated are the Danube (light blue), the coastlines (grey; Braconnot et al., 2012), and the ice sheet extents (turquoise; Cline et al., 1984).

The MARs and MAR10 (Table S1 and Fig. 4) were reconstructed from samples that were extracted during fieldwork campaigns from loess paleosol sites. The MAR for a specific site was inferred by taking into account all particles found in the respective sample, independent of their diameter. In contrast, the MAR10 for the same site was inferred by taking into account only particles up to 10 µm in diameter. Most of the MARs and MAR10 (Table S1 and Fig. 4) result from sites of the European loess belt. This belt plays a key role in assessing palaeoclimatic dust cycle simulations for Europe (Kukla, 1977; Little et al., 2002; Haase et al., 2007; Sima et al., 2009). During the LGM, it corresponded approximately to the fraction of the European land area that was bounded northwards by the EIS and southwards by the Alps, Dinaric Alps, and Black Sea. The F_D20 in Fig. 4a and b (F_D12 in Fig. 4c and d) denotes the WRF-Chem-LGM deposition rates caused only by particles smaller than 20 µm (12 µm) in diameter. To distinguish the deposition rates obtained from the two downscaling methods, the F_D20DD and F_D12DD relate to the dynamic, while the F_D20SD and F_D12SD relate to the statistic dynamic downscaling simulations.

For central Europe, the dynamic (Fig. 4a and c) and statistic dynamic downscaling (Fig. 4b and d) resulted in similar F_D values, confirming the suitability of the statistic dynamic downscaling.

Figure 5Dust emission rates for (a) winter (DJF), (b) spring (MAM), (c) summer (JJA), and (d) autumn (SON) during the Last Glacial Maximum. This reconstruction is based on dynamic downscaling. The Danube (light-blue line) and the extent of the continental ice sheets (white) are shown.

During the LGM, the largest F_D20 (>10⁵ $\mathrm{g}{\mathrm{m}}^{-\mathrm{2}}{\mathrm{yr}}^{-\mathrm{1}}$) occurred in western Poland (Fig. 4a). Slightly lower F_D20 values (10⁴–10⁵ $\mathrm{g}{\mathrm{m}}^{-\mathrm{2}}{\mathrm{yr}}^{-\mathrm{1}}$) were found in adjacent areas, in eastern Germany, for example. F_D20 was 10³–10⁴ $\mathrm{g}{\mathrm{m}}^{-\mathrm{2}}{\mathrm{yr}}^{-\mathrm{1}}$ on the North German Plain, in the dried-up German Bight, eastern England, northern and western France, the Benelux region, and southeast of the Carpathians. Regional deposition maxima of 10³–10⁴ $\mathrm{g}{\mathrm{m}}^{-\mathrm{2}}{\mathrm{yr}}^{-\mathrm{1}}$ occurred along the French LGM coastline (46–48^∘ N), on the eastern side of the Carpathians (44–47^∘ N, including the eastern Romanian Danube Plain), and near the Caucasus (44–45^∘ N, Fig. 4a). They coincide with today's extensive loess derivates along the Atlantic coastline of France and at the European foothills north of 42^∘ N, with the loess thickness maximum in the Romanian Danube Plain (Haase et al., 2007; Jipa, 2014). The quality of the simulations is also recognizable in the Carpathian Basin, which is now half covered with loess and clay of aeolian origin (Varga et al., 2012). There, the simulated F_D20 values of 100–1000 $\mathrm{g}{\mathrm{m}}^{-\mathrm{2}}{\mathrm{yr}}^{-\mathrm{1}}$ (Fig. 4a) are in good agreement with the MARs (200–500 $\mathrm{g}{\mathrm{m}}^{-\mathrm{2}}{\mathrm{yr}}^{-\mathrm{1}}$). In Ukraine and at the eastern margins of the EIS, F_D20 values of 100–1000 $\mathrm{g}{\mathrm{m}}^{-\mathrm{2}}{\mathrm{yr}}^{-\mathrm{1}}$ are in line with the MARs (Fig. 4a). Over Ukraine and consistent with our results, dust transport and deposition by east sector winds are evidenced by loess deposits on the west bank of the Dnieper (Sima et al., 2013).

The MARs of a few loess sites are higher than the F_D20 in their surroundings. Such an underestimation could be explained by particles larger than 20 µm, which are not taken into account by the F_D20. For some regions, the MARs of closely related sites vary over orders of magnitude, e.g. between 10² and 10⁴ $\mathrm{g}{\mathrm{m}}^{-\mathrm{2}}{\mathrm{yr}}^{-\mathrm{1}}$ near the Rhine and in Belgium (Fig. 4a). This may be due to strong small-scale variability, loess dating uncertainties (Singhvi et al., 2001; Renssen et al., 2007), or age model inaccuracies (Bettis et al., 2003). For western Germany, a transition from higher F_D20 (10³–10⁴ $\mathrm{g}{\mathrm{m}}^{-\mathrm{2}}{\mathrm{yr}}^{-\mathrm{1}}$) in the northeast to lower F_D20 (10²–10³ $\mathrm{g}{\mathrm{m}}^{-\mathrm{2}}{\mathrm{yr}}^{-\mathrm{1}}$) in the southwest was found (Fig. 4a). For a few sites in southwestern Germany, Austria, Ukraine, and along the Danube, F_D20 is an order of magnitude lower than the respective MARs (Fig. 4a). Given the 50 km grid spacing of the WRF-Chem-LGM simulation, this may be attributed to missing local dust sources, such as dried-up riverbeds and floodplains. Possibly, the MARs of these sites are also inferred from particles that were predominantly larger than 20 µm, yet data on particle sizes are not available. The peak deposition locations and the overall shape of the F_D20 and F_D12 patterns are very similar (Fig. 4). The F_D12 values are also consistent with the MAR10 almost everywhere (Fig. 4c and d). Those F_D12 values that overestimate the MAR10 do not contradict the consistency since the F_D12 also takes into account particles that are (by definition) excluded by the MAR10. In summary, high consistency was found between the simulated dust deposition rates and the MARs and MAR10 that were reconstructed from on-site samples.

Figure 6Dust deposition rates for (a) winter (DJF), (b) spring (MAM), (c) summer (JJA), and (d) autumn (SON) during the Last Glacial Maximum. This reconstruction is based on dynamic downscaling. Ice sheet extents (turquoise; Cline et al., 1984), Danube (light-blue line), and coastlines (grey; Braconnot et al., 2012) are delineated.

Figure 7Snow cover (%, left column), vegetation cover (%, centre), and soil moisture (m³ m⁻³, right), resolved for winter (DJF), spring (MAM), summer (JJA), and autumn (SON) for the Last Glacial Maximum. These reconstructions are based on dynamic downscaling.

3.5 Seasonal dust cycle patterns

During the LGM, the strongest emission and deposition in Europe occurred in summer, followed by autumn and spring (Figs. 5 and 6). The areas with the overall highest emission were also those with the highest seasonal emission (Figs. 3 and 5). The spring and winter emissions have the same order of magnitude. The low winter and spring emission rates along the EIS margin were caused by the then extensive snow cover there. During winter, emissions peaked only in northern France, consistent with its little snow cover and the vegetation cover (Fig. 7) that was prescribed to the WRF-Chem-LGM. Major dust emissions occurred from the Carpathian Basin and along the northwest coast of the Black Sea. During spring, slightly attenuated emissions are simulated for France, despite the decreasing snow cover but in accordance with its increasing vegetation cover. Considerably higher emission rates are simulated from along the German and Polish EIS margin where the snow cover had retreated. For eastern Europe, the growing vegetation cover and the slight soil moisture increase account for partly lower spring than winter emission rates. The soil moisture increase possibly resulted from meltwater of the retreating snow cover. The highest emission rates occurred during summer and were located along the German and Polish EIS margin. Slightly lower emissions are found to the east of the EIS. These findings are in coherence with the surface properties of these areas during summer, i.e. they were mostly snow free and the least moist. During autumn, the snow cover increased, causing a decrease in dust emissions, except for the area north of the Black Sea which encountered its annual maximum. This maximum can be attributed to the retreat of the vegetation cover and the dry soil conditions there.

Figure 8Dust emission rate fractions caused by the (a) northeasters, easterlies, and southeasters; (b) cyclonic regimes; (c) southwesters, westerlies, and northwesters; and (d) southerlies during the Last Glacial Maximum. The simulated emission rates are weighted according to the occurrence frequency of the associated wind regime(s) in the Max Planck Institute Earth System Model (Table 2). Dust particles up to 20 µm in diameter have been considered. The Danube (light-blue line) and the extent of the continental ice sheets (white) are shown.

Figure 9Dust deposition rate fractions caused solely by the (a) northeasters, easterlies, and southeasters; (b) cyclonic regimes; (c) southwesters, westerlies, and northwesters; and (d) the southerlies during the Last Glacial Maximum. The simulated deposition rates are weighted according to the occurrence frequency of the associated wind regime(s) in the Max Planck Institute Earth System Model (Table 2). Dust particles up to 20 µm in diameter have been considered. The ice sheet extents (turquoise; Cline et al., 1984), the Danube (light blue), and the coastlines (grey; Braconnot et al., 2012) are delineated.

The winter CWT distribution indicates prevailing east sector winds (37 %) in contrast to cyclonic regimes, which occurred much less frequently than on an annual average (13 %; Tables 2 and 3). The winter deposition rates northwest of the Alps were considerably above the annual average, while the rates at the central and eastern European EIS margin were below the annual average (Figs. 4 and 6a). In western Europe, the highest deposition rates occurred near the sources, yet a considerable dust fraction was also transported and deposited to the west and northwest of the sources, which requires east sector winds. Low deposition rates were found for southern France; however marked depositions occurred when subjected to cyclonic regimes (Fig. 9b). The deposition pattern for the central Mediterranean area (Italy, the Adriatic) suggests significant dust transport by east sector winds and anticyclonic winds, in sum prevailing 51 % of the time. In eastern Europe, considerable winter deposition rates covered areas south of the dust sources, in particular the western Black Sea and regions south of the Danube. This indicates a significant contribution to the dust transport by northerlies (6 %), northeasters (12 %), and the anticyclonic regimes (14 %).

Also, the spring deposition rates evidence the importance of the east sector winds (42 %, Table 2) for the dust cycle. In western Europe, major deposition areas are to the west and northwest of the sources, while they are to the west and southwest in eastern Europe (Fig. 6b). An increase in the dust transport towards the south in western Europe and towards the north in eastern Europe indicates an increasing role of the cyclonic regimes (27 %) during the spring.

The summer deposition rates are distributed zonally along the EIS margin, suggesting an approximately latitude-parallel dust transport by the west (21 %) and/or east sector (24 %) wind directions. In addition, the northern flanks of cyclonic regimes (24 %) likely contributed to a westward dust transport. Over the northeasternmost part of Europe (62^∘ N, 40^∘ E), the deposition rates suggest east sector winds. The autumn deposition rates over western and central Europe show a westward running plume from the southern EIS margin over Germany and Poland, corroborating the major role of the east sector winds (38 %) in the dust cycle. The high deposition rates in eastern Europe suggest that the cyclonic regimes (19 %) also contributed during autumn.

3.6 Wind-regime-based dust cycle decomposition

The wind regime occurrence frequency distribution (Table 2) demonstrates the temporal dominance of the east sector winds during the LGM. This temporal dominance likely shaped the dust cycle but the contribution of each wind regime type has so far not been analysed. This analysis is provided here by discussing the dust emission and deposition characteristics associated with different CWTs which reveal that the east sector winds caused by far the largest dust emission and deposition during the LGM (Figs. 8a and 9a). In sum, they generated an average dust emission of 1111 $\mathrm{g}{\mathrm{m}}^{-\mathrm{2}}{\mathrm{yr}}^{-\mathrm{1}}$ (Fig. 8a), which is more than twice the rate generated by cyclonic regimes (494 $\mathrm{g}{\mathrm{m}}^{-\mathrm{2}}{\mathrm{yr}}^{-\mathrm{1}}$, Fig. 8b). The west sector winds contributed on average even less to the dust cycle (375 $\mathrm{g}{\mathrm{m}}^{-\mathrm{2}}{\mathrm{yr}}^{-\mathrm{1}}$, Fig. 8c). Compared to the southerlies (232 $\mathrm{g}{\mathrm{m}}^{-\mathrm{2}}{\mathrm{yr}}^{-\mathrm{1}}$, Fig. 8d), this rate is low for a wind sector that sums the contribution of three wind directions (SW, W, NW).

The cyclonic wind regimes caused the most heterogeneously distributed emissions (Fig. 8b) with four main centres: the largest located in the German–Polish–Czech border region, another in eastern England, and the remaining two near the EIS margin in western Russia. This distribution resembles a subset of the emission distribution of the east sector winds (Fig. 8a). Together with the location of the CWT reference regions, this resemblance could be explained by the fact that all records classified as cyclonic must centre their cyclonic pressure distribution approximately around the central point for the CWT classification (47.5^∘ N, 17.5^∘ E). This implies that the corresponding emissions could have been triggered by easterlies on the northern flanks of the cyclones. Dust was hardly emitted from areas on the southern flanks of the cyclones which are commonly affected by fronts and precipitation (Booth et al., 2018). In addition to the dust emission areas that occurred equally during both regimes (cyclonic and east sector winds), the east sector winds also generated emissions in Austria, Slovakia, Hungary, Ukraine, central Germany, the Danube Basin, and the North Sea Basin. In contrast, the west sector winds produced a more homogeneous distribution of markedly smaller emission rates extending from western Ukraine to the French Atlantic coast. While northwesters with a strong northerly component most likely forced emissions from the German–Polish EIS margin, the west sector winds and the southerlies controlled the emissions from France, southwestern Germany, the English Channel, and the Alps foreland (Fig. 8c and d). The combination of the emission and deposition rate patterns of the east sector winds (Figs. 8a and 9a) indicates major westward dust transport along the southern and eastern EIS margin. The conic shape of the deposition rate distribution in western and central Europe (between 10² and 10³ $\mathrm{g}{\mathrm{m}}^{-\mathrm{2}}{\mathrm{yr}}^{-\mathrm{1}}$) suggests that these depositions can be attributed to emissions from more eastern sources. The east sector winds also deposited considerable amounts of dust in and south of the Danube Basin as well as along the Danube.

The deposition rates of the cyclonic regimes (Fig. 9b) indicate two main dust transport directions: westwards over central and eastern Europe and southwards over western Europe. The shape and location of the emission and deposition areas caused by the west sector winds are almost congruent (Fig. 8c and 9c). This implies that a unique dust transport direction cannot be inferred for this wind regime. Instead, dust may have been transported in various directions. Dust deposition in Ireland, western Great Britain, the Bay of Biscay, and near the eastern margin of the EIS even suggests a westward dust transport (Fig. 9c), implying that east sector winds may have occurred locally while possibly weak west sector winds prevailed over central Europe. More precisely, dust was transported westwards from Poland to eastern and central Germany, while it was carried southwards from eastern England to the English Channel and northwestern France up to the Pyrenees foreland. The depositions caused by southerlies show a northwestward transport over central Europe (Fig. 9d). Considerable amounts of dust (between 10³ and 10⁵ $\mathrm{g}{\mathrm{m}}^{-\mathrm{2}}{\mathrm{yr}}^{-\mathrm{1}}$) were transported from sources in western Poland, eastern Germany, and Czechia to northern Germany, Denmark, southern Sweden, and the North Sea Basin. The deposition pattern also suggests a northwestward transport in France.