The day-to-day weather conditions in the midlatitudes are strongly affected by the passage of extratropical cyclones, which are typically associated with precipitation, strong winds and changes in temperature and cloudiness. Cyclones also play a major role in the water cycle and the redistribution of momentum and energy in the climate system (Hoskins and Valdez, 1990; Chang et al., 2002). The assessment of cyclone activity, notably to analyse their paths, characteristics and impacts, is thus key to determine the day-to-day weather conditions, the regional mean climate and its variability on multiple timescales. In fact, there is a wide range of literature analysing case studies of extreme cyclones (e.g. Wernli et al., 2002; Ludwig et al., 2015), the mean cyclone activity in the midlatitudes in the recent past (e.g. Hoskins and Hodges, 2002; Ulbrich et al., 2009) and possible changes under future climate conditions (e.g. Bengtsson et al., 2009; Ulbrich et al., 2009). On the other hand, studies analysing the structural characteristics of extratropical storms from a climatological perspective are less frequent (e.g. Catto et al., 2010; Rudeva and Gulev, 2011; Dacre et al., 2012; Hewson and Neu, 2015; Sinclair et al., 2020). While some general concepts are available on how warmer climate conditions will affect the intensity and structure of cyclones, there are still several open questions, particularly regarding how dominant the increased latent heating may become compared to other physical processes like low-level and upper-level baroclinicity (see Catto et al., 2019; their Fig. 2).

The availability of studies addressing the characteristics of cyclone activity outside of the period extending from the mid-19th century to the end of the 21st century decreases sharply. Raible et al. (2018) analysed variations of cyclone statistics in a very long simulation with a fully coupled Earth system model from 850 to 2100 CE. While they identified variations on multiple timescales, they found no evidence for an external forcing imprint before 1850. Moreover, Pfahl et al. (2015) analysed cyclone activity in idealised aquaplanet simulations covering a wide range of possible climate conditions (from 270 to 316 K global mean temperatures). While the structure of the majority of the cyclones reveals only small changes on average, larger differences were identified for intense cyclones. For example, cross-front temperature differences are expected to be higher (lower) for considerably colder (warmer) climates (Pfahl et al., 2015, their Fig. 10), whereas the associated precipitation is expected to decrease (increase; their Fig. 12) due to the strong limiting effect of temperature on the atmospheric moisture content.

One important issue preventing non-recent or non-21st-century cyclone analysis is the availability of climate model output with sufficient spatial and temporal resolution to enable identification, tracking and characterisation of such cyclones. For example, model data from the PMIP3 project (Braconnot et al., 2012) are only archived 6-hourly for short (30-year) time slices and typically at a low resolution (approximately 200–300 km). Most pre-20th-century studies consider aggregated measures of cyclone (or synoptic) activity (e.g. Kageyama et al., 1999; Laine et al., 2009; Hofer et al., 2012; Ludwig et al., 2016), thus not enabling a detailed comparison to regional temperature and precipitation variability. One period of particular interest is the Last Glacial Maximum (LGM; Clark et al., 2009), when the European climate was characterised by colder and mostly drier conditions (Bartlein et al., 2011; Annan and Hargreaves, 2013; Újvári et al., 2017; Cleator et al., 2019). Large parts of northern Europe were covered by permanent ice sheets and surrounded by polar desert conditions (Ray and Adams, 2001). Western, central and eastern Europe were largely characterised by open shrublands and grasslands (steppe tundra), while in southern Europe steppe with embedded forest (forest steppe) dominated (Ray and Adams, 2001). Under these conditions, dust storms triggered by strong winds must have been common in Europe, as documented by the major loess deposits found primarily around 50∘ N over western and central Europe and over large parts of eastern Europe (Antoine et al., 2009, 2013; Sima et al., 2013; Újvári et al., 2017).

Under the influence of the continental ice sheets and extended sea ice, the PMIP3 general circulation models (GCMs) show stronger meridional temperature gradients, leading to a southward-displaced, more intense and less variable North Atlantic jet than under current climate conditions (Löfverström et al., 2014, 2016; Merz et al., 2015; Wang et al., 2018). These differences have been related, e.g. to more dominant cyclonic Rossby wave breaking near Greenland (Riviére et al., 2010), stationary wave packets trapped in the midlatitude wave guide (Löfverström, 2020) and enhanced meridional eddy momentum flux convergence over the North Atlantic (Wang et al., 2018). In line with a southward-displaced and stronger jet stream, several studies show a more intense and southward-shifted North Atlantic storm track compared to today's climate (e.g. Hofer et al., 2012; Luetscher et al., 2015; Ludwig et al., 2016). However, other studies display reduced storm track activity over the North Atlantic in spite of the enhanced baroclinicity (e.g. Donohoe and Battisti, 2009; Riviére et al., 2010; Löfverström et al., 2016). Riviére et al. (2018) discusses a reduced baroclinic conversion as a possible reason for this apparent discrepancy, arguing that the eddy heat fluxes are less well aligned with the mean temperature gradient for LGM than for the pre-industrial (PI) conditions. Other arguments for the reduced storminess include model resolution, parameterisations and boundary conditions (e.g. Donohoe and Battisti, 2009; Riviére et al., 2018). Thus, the intensity differences between LGM and PI North Atlantic storm track activity may be model dependent.

The PMIP3 models show indications that while Europe was largely drier than today, this was not the case for some regions, notably for the Iberian Peninsula (e.g. Hofer et al., 2012; Beghin et al., 2016; Ludwig et al., 2016). However, wetter conditions over the Iberian Peninsula are not in line with (most of) the proxy data, which themselves are often associated with considerable uncertainties (e.g. Bartlein et al., 2011; Moreno et al., 2014; Cleator et al., 2019). Nevertheless, the substantial misrepresentation of the regional climate for the LGM in PMIP3 models compared to proxies is regarded as a general issue (Harrison et al., 2015). In some cases, such caveats can be partially traced back to shortcomings of the GCMs and/or their boundary conditions. For example, Ludwig et al. (2017) implemented more realistic boundary conditions in terms of the North Atlantic sea surface temperatures (SSTs), land use types and vegetation cover in a regional climate model (RCM) to simulate the regional climate under LGM conditions. Their results in terms of LGM temperature, precipitation and the permafrost margin are in better agreement with the proxies than without the implemented boundary conditions. Still, further studies are needed, notably at the regional scale (e.g. Ludwig et al., 2018), in order to further our confidence in the modelling capabilities and our understanding of the palaeoclimate conditions for Europe in key periods like the LGM (Harrison et al., 2015, 2016; Ludwig et al., 2019).

The present work aims to advance our understanding of the LGM climate over the North Atlantic and Europe through a more detailed analysis of the cyclonic activity and its associated impacts, notably in terms of precipitation, temperature and wind speed. The LGM cyclones are first identified and tracked on a simulation with the coupled MPI-ESM-P model, for which data with high temporal resolution were archived. Secondly, a subsample of extreme cyclones is downscaled with a RCM to analyse possible changes in LGM cyclone characteristics compared to their modern counterparts at high spatial and temporal resolution. The identified characteristics of LGM extreme cyclones are discussed in terms of the available proxies for LGM climate across western Europe. In addition to the precipitation and temperature, the importance of the dominant land cover conditions and the frequent occurrence of dust storms is evaluated. The final section presents the summary and main conclusions.

The starting point of our analysis is data from the third phase of the Paleoclimate Modeling Intercomparison Project (PMIP3) (Braconnot et al., 2012) (http://pmip3.lsce.ipsl.fr/, last access: 28 March 2020). The simulations were performed according to the PMIP3 21 ka experimental design, which includes the lower sea level and blended ice sheet data (Peltier et al., 2015; Lambeck and Chappell, 2001; Lambeck et al., 2002; Tarasov and Peltier, 2002, 2003), orbital parameters and lower greenhouse gas concentrations (see Table 1). From the PMIP3 GCMs, we have selected the MPI-ESM-P (Stevens et al., 2013; Jungclaus et al., 2013) constant forcing simulations, for which two 30-year time slices with 6-hourly output data are available for PI and LGM conditions. This choice was motivated by the availability of 6-hourly 3-D model level data needed for running the RCM. Ludwig et al. (2016) recently analysed a small ensemble of PMIP3 models in terms of large-scale circulation for Europe, jet stream synoptic activity, precipitation and temperature for the LGM. While the Max Planck Institute for Meteorology Earth System Model in Paleo Mode (MPI-ESM-P) has a slightly different jet structure to some of the other PMIP3 models (cf. Ludwig et al., 2016; their Fig. 2), which also impacts, for example, the storm track and precipitation (their Figs. 4, 6), its main characteristics are generally close to the ensemble average.

Individual extratropical cyclones over the North Atlantic and Europe are identified and tracked based on 6-hourly mean sea level pressure data with a widely used automatic tracking algorithm (Murray and Simmonds, 1991; Pinto et al., 2005). The resulting cyclone statistics provide information on the lifetime of each identified cyclone, thus enabling the computation of mean cyclone statistics like track density, mean maximum intensity, cyclogenesis, cyclolysis, propagation speed and deepening rates. Cyclone statistics obtained with this method compare well with other methodologies (e.g. Neu et al., 2013; Hewson and Neu, 2015). Following Pinto et al. (2009), cyclones are selected based on the following conditions: (a) cyclone lifetime of at least 24 h, (b) a minimum core mean sea level pressure (MSLP) value below 1000 hPa, (c) a maximum vorticity (approximated by the Laplacian of MSLP) value above 0.6 hPa deg. lat.-2 and (d) a maximum deepening rate of 0.3 hPa deg. lat.-2 s-1 is achieved at least once during their lifetime. The method is applied here to the MPI-ESM-P data for the extended (ONDJFM) winter season. In order to analyse the characteristics of the most extreme cyclones affecting Europe in more detail, the most intense 30 cyclones (top 30) for the PI and LGM periods are selected based on their peak intensity in terms of vorticity and their passage within a pre-defined box over the eastern North Atlantic (Fig. 1, dashed box). The selection of the box enables the creation of a cyclone ensemble that impacts western Europe and permits a comparison with terrestrial proxy data, e.g. for precipitation and dust.

The Weather Research and Forecasting (WRF) model (Skamarok et al., 2008) is used in version 3.9.1.1 to simulate the top 30 cyclones (from PI and LGM) with a grid spacing of 12.5 km (including 35 vertical layers up to 30 hPa). To achieve a grid spacing of 12.5 km, a two-step nesting approach is necessary. Cyclones were initially simulated on a 50 km grid forced by MPI-ESM-P data as initial and boundary conditions with an update frequency of 6 h. The final 12.5 km grid spacing is achieved by a second nesting step within the WRF model. An overview of the parameterisation choices is given in Table 2. For the calculation of wind gusts, a gust parameterisation based on 10 m wind speed and friction velocity (Schulz and Heise, 2003; Schulz, 2008) has been implemented into the WRF model. This gust parameterisation shows overall good agreement with observed wind gusts, particularly over flat terrain (Born et al., 2012). For the WRF simulations, global PI and LGM boundary conditions were adapted considering specifications of the PMIP3 protocol (Braconnot et al., 2012; see also Ludwig et al., 2017). These changes encompass orbital parameters, trace gases (see Table 1), the consideration of ice sheets (extent and height), an associated lowering of the sea level and adaption of land use cover (CLIMAP Project Members, 1984).

The top 30 cyclone tracks simulated by WRF were identified manually based on relative vorticity field at 850 hPa. For comparison of the PI and LGM cyclone characteristics, and following the methodology from Catto et al. (2010; their Fig. 3), each track was rotated so that the cyclones were each moving in west–east directions, enabling the generation of composites for different atmospheric variables (cf. Dacre et al., 2012). Composites have been created for peak intensity (0) and 6, 12, 18 and 24 h before peak intensity, and 6 and 12 h afterwards. For brevity, we will primarily discuss the time frames of (i) 12 h before peak intensity and (ii) peak intensity. The variables analysed from the 12.5 km WRF simulations include mean sea level pressure, precipitation, column-integrated water vapour, 850 hPa equivalent-potential temperature, 925 hPa winds, and near-surface wind gusts.

In this section, we analyse the general characteristics of cyclones over the North Atlantic and Europe under LGM conditions and compare them to PI climate conditions. Figure 2 shows the cyclone track density for the extended winter for PI and LGM climate conditions. In spite of the lower spatial resolution of MPI-ESM-P, the cyclone track density for PI is close to cyclone statistics obtained with reanalysis datasets, with a slight southerly shift of cyclonic activity (see Fig. S1 in the Supplement for comparison with NCEP reanalysis data of Kalnay et al., 1996) and CMIP GCMs for recent climate conditions (cf. Pinto et al., 2007; their Fig. 1). Still, some regional shortcomings are identified, notably the limited cyclone activity over the Mediterranean basin. The North Atlantic storm track shows a clear tilt towards northern Europe and the Arctic Ocean for PI, and its location and orientation are closely related with the eddy-driven jet stream (black contours in Fig. 2a) and the associated upper-air baroclinicity (Hoskins and Valdes, 1990; Pinto et al., 2009). A comparison of the jet stream between MPI-ESM-P PI and NCEP reanalysis data shows a slight tilt towards Europe by the MPI model (Fig. S2), in line with the enhanced (reduced) southward (northward) cyclone activity (Fig. S1).

The North Atlantic storm track looks quite different under LGM conditions relative to PI: the cyclone track density is higher over the North Atlantic and more constrained to the ice edge (Figs. 2b, S1). Close to Europe, a bifurcation is found, and cyclones are either deflected northward along the border of the Scandinavian ice sheet or southeastward towards central Europe and the Mediterranean (Fig. 2c). In accordance, the eddy-driven jet is stronger under LGM conditions in the MPI-ESM-P (Fig. S2), thus establishing more favourable conditions for the occurrence of intense storms affecting western and central Europe. For the North Atlantic (35–70∘ N, 70∘ W–0), the total number of cyclones for the analysed 30-year period is about 26 % larger for LGM than for PI conditions (12071 vs. 9541 individual cyclone counts in 30 years, corresponding to roughly 2.2 cyclones d-1 vs. 1.8 cyclones d-1).

Other properties of cyclone activity are depicted in Fig. 3 for LGM in terms of absolute values (left column) and their differences to PI (right column). Cyclogenesis is dominant along the North American east coast for the LGM (Fig. 3a) and much stronger than in PI (Fig. 3b), which is in line with a much stronger upper-level jet stream during glacial conditions. Moreover, the rate of cyclogenesis increased south of Greenland and over western and central Europe. On the other hand, there is more cyclolysis along the borders of the Greenland and Scandinavian ice sheets (Fig. 3c) relative to the PI conditions (Fig. 3d). Mean maximum cyclone intensity is typically attained in a region extending from Newfoundland to Iceland and the British Isles, with a secondary maximum over eastern Europe (Fig. 3e). Compared to the PI cyclones, the LGM cyclones have stronger intensities, particularly in an area extending from the south of Greenland to the British Isles and over most of continental Europe. On the other hand, cyclone intensity close to the North American east coast is considerably lower (Fig. 3f). Deepening rates are stronger for LGM cyclones over the central north Atlantic, as well as their filling rates close to the ice edge/ice sheets (Fig. 3g, h). These results point towards different typical development of LGM cyclones compared to their PI counterparts, which occurs either more zonally at lower latitudes towards central Europe or further downstream closer to the ice edge towards the Arctic.

Figure 4 displays the relative frequency distribution of the cyclone intensity over the North Atlantic area (35–70∘ N, 70∘ W–0), revealing that LGM cyclones are on average more intense (mean (median): 1.58 (1.41) hPa deg. lat.-2 s-1) than their PI counterparts (mean (median): 1.43 (1.28) hPa deg. lat.-2 s-1). In particular, the number of cyclones exceeding 3 hPa deg. lat.-2 s-1 is twice as large for LGM as for PI. In fact, a small number (18) of LGM cyclones attain intensities exceeding the range identified for PI cyclones. The statistics for the region close to Europe (box) are similar (not shown). All these results document a shift towards stronger intensities for LGM cyclones, both in terms of average numbers and extreme values.

To analyse the characteristics of cyclones for the LGM, a subset of 30 extreme cyclones passing over the eastern North Atlantic is selected for both the PI and LGM periods based on the highest values for vorticity (Laplacian of MSLP) within the selected box. The trajectories of the selected MPI-ESM cyclones are depicted in Fig. 5 and key information is given in Table 3. Two important facts are clearly identifiable: LGM extreme cyclone trajectories are more zonally orientated and constrained to a narrower corridor (particularly until 15∘ W) than their PI counterparts, and they achieve higher vorticity (mean Laplacian of MSLP for LGM: 2.80; PI: 2.52) values during lifetime (Table 3).

For each of these cyclones, WRF simulations down to 12.5 km grid spacing are performed along the most intense segment of their lifetimes, thus gaining 3-D data to analyse the cyclones, e.g. in terms of the evolution of their structure, air masses, winds and precipitation. Care was taken that storms agree between the original MPI-ESM-P data and WRF cyclone tracks. Generally, the obtained tracks based on the WRF simulations reveal lower core pressures and higher vorticity than their low-resolution counterparts. Figure 6 shows a comparison between high- and low-resolution data for a selected cyclone. The considered cyclone tracks (MPI-ESM LGM cyclone no. 24) show a good superimposition for MPI-ESM, WRF 50 km and WRF 12.5 km and good agreement of the position of maximum intensity with a slightly further southeastward location for WRF 50 km (Fig. 6a). While the time series, centred at peak intensity, for the development of core pressure show a similar strong pressure drop for all resolutions, the relative vorticity exhibits much stronger values for the high-resolution simulation (WRF 12.5 km) in comparison to the coarser realisations. Note that the WRF simulations do not cover the whole trajectory of the cyclones, and thus the time series are shorter. Figure 6c and d depict the corresponding precipitation patterns at peak intensity, in this case revealing more small-scale structures and higher precipitation values for the 12.5 km WRF simulation.

Based on the top 30 cyclones for PI and LGM and using the composite methodology described in Sect. 2, average cyclone characteristics are investigated for the intensification phase, focussing on the time frames (i) 12 h before peak intensity and (ii) peak intensity. Figure 7 displays the MSLP fields for PI and LGM for peak intensity; the corresponding panels for 12 h before peak intensity can be found in Fig. S3. The anomalies compared to the mean MSLP are displayed in colours, in order to analyse MSLP gradients (Figs. 7a, b; S3a, b). The core MSLP of the LGM has higher values compared to PI (LGM: 982.7 hPa, PI: 975.2 hPa). This can be explained by the lower sea levels (∼120 m) for LGM, causing a difference of global mean sea level pressure of about ∼13 hPa (PI: 1010 hPa, LGM: 1023 hPa). Taking this into account, the LGM cyclones reach deeper MSLP values compared to the global average MSLP, being consistent with the stronger deepening rates for LGM cyclones identified by the tracking algorithm (Fig. 3g, h). Additionally, the closer isobars south of the cyclone core indicate stronger pressure gradients for the LGM cyclones, which is supported by the LGM – PI differences (Fig. 7c). This is particularly the case on the expected location of the frontal areas. Figures 7d, e and S3d, e display the anomalies of equivalent-potential temperature (θE) at 850 hPa for LGM relative to PI. Results show that the cross-frontal gradients are particularly intense across the warm front and that the whole development is displaced forward in the cyclone for LGM conditions, indicating a faster occlusion (Fig. 7f). On the other hand, the total water content is much higher under PI conditions, primarily due to the effect of the higher environmental temperatures (Figs. 7g, h; S3g, h) with differences in the warm sector reaching up to 10 and 6–8 mm close to the cyclone core at peak intensity (Fig. 7i).

The strong difference in available water content again leads to a large difference in terms of accumulated precipitation, which is clearly larger for the PI composites (Figs. 8a, b; S4a, b). While some of the smaller deviations can be potentially attributed to slightly different developments, the total precipitation is considerably lower for the LGM extreme cyclones (up to 1.7 mm h-1) particularly in the area of peak precipitation close to the cyclone centre (Fig. 8c). On the other hand, wind speed at 925 hPa is much higher for LGM cyclones (Figs. 8d, e; S4d, e), where wind speeds over large areas south of the cyclone exceed their PI counterparts by 10–12 ms-1, particularly for 12 h before peak intensity (Figs. 8f; S4f). Strong differences are also revealed for near-surface wind gusts (Figs. 8g–i; S4g–i). In this case, the wind gusts are particularly enhanced along the expected location of the cold front at peak intensity, with deviations exceeding 5 ms-1. The above-described patterns remain true for the cyclone characteristics 12 h after peak intensity (Figs. S5, S6). For example, stronger wind gusts remain dominant in the area south of the cyclone core (Fig. S4i).

In summary, LGM cyclones display steeper MSLP gradients (in agreement with higher intensity in terms of circulation), larger temperature gradients between the air masses, weaker precipitation and stronger wind gusts than their PI counterparts do (see also the summary for wind and precipitation in Table S1 in the Supplement).

The above analysis shows that extreme cyclones under LGM conditions were typically more intense than PI extreme cyclones. Likewise, they were associated with larger frontal temperature gradients, stronger winds and reduced precipitation. In this section, we analyse in how far these cyclone characteristics can help to explain the climate conditions in western and central Europe at the LGM and the frequent occurrence of dust storms in that area.

According to the available proxy records, western and central Europe were colder and largely drier under LGM conditions (Bartlein et al., 2011; Annan and Hargreaves, 2013; Cleator et al., 2019) and were to a large extent covered by open shrublands and grasslands (Ray and Adams, 2001). The colder and generally drier conditions are reproduced by the PMIP3 GCMs, but considerable differences in detail have been identified compared to proxy-based climate reconstructions (Beghin et al., 2016; Ludwig et al., 2016; Cleator et al., 2019). These can be both attributed to shortcomings of the GCMs and/or uncertainties in the reconstructions, which are often not well constrained (Bartlein et al., 2011; Ludwig et al., 2019; Cleator et al., 2019). For example, most PMIP3 GCMs show enhanced precipitation over the Iberian Peninsula compared to PI climate (Beghin et al., 2016; Ludwig et al., 2016), which disagrees with the proxy data (e.g. Bartlein et al., 2011). Later, Ludwig et al. (2017, 2018) show that the wet bias over the Iberian Peninsula and western Europe can strongly be reduced by considering more realistic boundary conditions (particularly in terms of sea surface temperatures, land use and vegetation cover) in high-resolution RCM simulations.

The present results enable a more detailed evaluation of the above hypotheses and interpretations. Under current climate conditions, precipitation in the midlatitudes is largely associated with the passage of extratropical cyclones, where extreme cyclones have a comparatively larger contribution to total/extreme precipitation (Pfahl and Wernli, 2012; Hawcroft et al., 2012). Hawcroft et al. (2016) have provided evidence that GCMs typically have shortcomings in representing the precipitation associated with midlatitude cyclones. It is reasonable to assume that those shortcomings may be exacerbated at lower resolutions, e.g. those typical of the PMIP3 models. An idealised study by Pfahl et al. (2015) has shown that precipitation increases (decreases) disproportionately in considerably warmer (colder) climate conditions. With the help of the high-resolution WRF simulations, we could provide evidence that extreme cyclones under LGM conditions indeed induce considerably less precipitation than their PI counterparts (∼22 % within a 10∘ radius around the cyclone centre). Even though lower cyclone-related precipitation may be partially compensated by enhanced moisture advection embedded in the (stronger) westerly large-scale flow, particularly for areas where orographic precipitation dominates (e.g. upward slopes of mountain ranges/glaciers), it is consistent with the hypothesis of a drier western and central Europe. The view of a drier Europe is also consistent with the dominant land cover types estimated from proxy data (Ray and Adams, 2001) or a statistical reconstruction based on temperature and precipitation (Shao et al., 2018), namely polar desert close to the glaciers, forest steppe over southern Europe and steppe in between. The fact that we have selected an area close to the Iberian Peninsula for cyclone selection leads to the strong assumption that largely drier conditions must also have been found for southwestern Europe, in agreement with the proxies (e.g. Bartlein et al., 2011; Moreno et al., 2014).

Mineral dust plays an important role in our climate system (Shao et al., 2011). Dust emissions are typically initiated by the wind stress on land surfaces with little to no vegetation cover and easily erodible soils (e.g. Prospero et al., 2002). Such areas were very common in Europe under LGM conditions (e.g. Ray and Adams, 2001; Ugan and Byers, 2007), when the global dust cycle is estimated to have been stronger than at present (Maher et al., 2010). A large number of loess deposits over western and central Europe (particularly around 50∘ N; Antoine et al., 2009) indicate that dust storms were a common feature of the European climate. In particular, Antoine et al. (2009) identified cyclic variations in loess deposition between 34 and 17 ka at several sites in France, Germany and Belgium, with particularly high sedimentation rates, and attributed these to numerous and intense dust storms in periods with stronger large-scale flow and reduced precipitation. Furthermore, the high accumulation rates of loess in the middle and lower Danube basin indicate cold, dry and windy conditions during the LGM in southeastern Europe (Fitzsimmons et al., 2012), consistent with increased storm activity over central Europe (Fig. 2c). These and other findings document a more intense (global) dust cycle for LGM conditions (e.g. Albani et al., 2016; Újvári et al., 2017; Albani and Mahowald, 2019). The occurrence of dust storms over western and central Europe has been conceptually associated with the passage of intense extratropical cyclones penetrating deep into the continent (Antoine et al., 2009; their Fig. 12). Following on previous studies (e.g. Laine et al., 2009; Hofer et al., 2012; Ludwig et al., 2016), the present results provide evidence for the first time that individual LGM cyclones would be indeed capable of triggering such dust events: their frequent tracks over western/central Europe and strong wind could easily trigger dust emission and transport over short (for coarse grain) and large (for fine grain material) distances (cf. Shao et al., 2011). As moisture acts to make surface dust particles more cohesive (e.g. Ishizuka et al., 2008), the reduced cyclone precipitation and higher wind speeds in LGM cyclones would have actually been more conducive to generating dust storms. In addition to the role of the westerlies and embedded cyclones into generating dust storms in Europe, there is evidence that situations with persistent easterlies associated with anticyclonic flow triggered by a strong anticyclone over the Scandinavian ice sheet may have also played a significant role for loess deposition not only over eastern Europe but also over central Europe (Újvári et al., 2017; Schaffernicht et al., 2019).

The statistics and characteristics of extratropical cyclones over the North Atlantic and western Europe were analysed for time-slice experiments for PI and LGM conditions. First, the statistics of the climatologies of PI and LGM cyclones are analysed and compared based on global MPI-ESM-P simulations. Second, the characteristics of extreme LGM cyclones over the eastern North Atlantic were analysed in detail based on high-resolution simulations (12.5 km grid spacing) with the RCM WRF and compared to their PI counterparts. The results were discussed with available proxy reconstructions of climate parameters and vegetation types. The main conclusions are as follows:

The North Atlantic storm track was more intense under LGM conditions, featuring more frequent and intense synoptic systems than under PI conditions. One of the downstream branches brought more often extreme cyclones towards western and central Europe and the Mediterranean area.

Even though the added value of RCMs in palaeoclimate applications is still controversially discussed (Armstrong et al., 2019), there is a general call for improvements towards a new generation of reliable regional projections (e.g. Harrison et al., 2015; Kageyama et al., 2018). The present and other studies (see Ludwig et al., 2019 for a review) provide clear arguments for the extended use of RCMs in the scope of palaeoclimate studies, as they can play an important role towards a meaningful joint interpretation of proxies and climate model data. The upcoming new palaeoclimate simulations within CMIP6/PMIP4 (Kageyama et al., 2017, 2018), as well as new proxy-based reconstructions of climates (Cleator et al., 2019), will provide novel possibilities to expand our understanding of past climates and to reduce uncertainties on both the numerical and reconstruction branches.