Modeling dust emission response to MIS 3 millennial climate variations from the perspective of East European loess deposits
1Laboratoire de Météorologie Dynamique, UMR8539, INSU-CNRS, & CERES-ERTI, Ecole Normale Supérieure, 24 rue Lhomond, 75231 Paris cedex 5, France
2Laboratoire des Sciences du Climat et de l'Environnement, UMR8212, CNRS-CEA-UVSQ, CE Saclay, l'Orme des Merisiers, Bât. 701, 91191 Gif-sur-Yvette cedex, France
3Lamont-Doherty Earth Observatory of Columbia University, Palisades, NY 10964, USA
4Laboratoire de Géographie Physique, UMR8591, CNRS – Université Paris I, place A. Briand, 92158 Meudon cedex, France
Abstract. European loess sequences of the last glacial period (~ 100–15 kyr BP) show periods of strong dust accumulation alternating with episodes of reduced sedimentation, favoring soil development. In the western part of the loess belt centered around 50° N, these variations appear to have been caused by the North Atlantic rapid climate changes: the Dansgaard-Oeschger (DO) and Heinrich (H) events. It has been recently suggested that the North-Atlantic climate signal can be detected further east, in loess deposits from Stayky (50° 05.65' N, 30° 53.92' E), Ukraine. Here we use climate and dust emission modeling to investigate this data interpretation. We focus on the areas north and northeast of the Carpathians, where loess deposits can be found, and the corresponding main dust sources must have been located as well. The simulations, performed with the LMDZ atmospheric general circulation model and the ORCHIDEE land-surface model, represent a Greenland stadial, a DO interstadial and an H event respectively. Placed in Marine Isotope Stage 3 (~ 60–25 kyr BP) conditions, they only differ by the surface conditions imposed in the North Atlantic between 30° and 63° N. The main source for the loess deposits in the studied area is identified as a dust deflation band, with two very active spots located west–northwest from our reference site. Emissions only occur between February and June. Differences from one deflation spot to another, and from one climate state to another, are explained by analyzing the relevant meteorological and surface variables. Over most of the source region, the annual emission fluxes in the "interstadial" experiment are 30 to 50% lower than the "stadial" values; they would only be about 20% lower if the inhibition of dust uplift by the vegetation were not taken into account. Assuming that lower emissions result in reduced dust deposition leads us to the conclusion that the loess-paleosol stratigraphic succession in the Stayky area reflects indeed North-Atlantic millennial variations. In the main deflation areas of Western Europe, the vegetation effect alone determined most of the ~ 50% stadial-interstadial flux differences. Even if its impact in Eastern Europe is less pronounced, this effect remains a key factor in modulating aeolian emissions at millennial timescale. Conditions favorable to initiating particularly strong dust storms within a few hundred kilometers upwind from our reference site, simulated in the month of April of the "H event" experiment, support the identification of H events as layers of particularly coarse sedimentation in some very detailed profiles.