Le premier auteur tient à remercier ses co-auteurs ainsi que Sylvie Coutard, Emilie Goval et Jean-Luc Locht pour l’aide qu’ils lui ont apportée sur le terrain, notamment lors du prélèvement des échantillons malacologiques. Les auteurs remercient également Nicole Limondin-Lozouet et Jef Vandenberghe pour leurs commentaires constructifs ayant permis d’améliorer la version initiale de cet article.

1 – Introduction

1Les plaines du nord de la France se distinguent par la présence d’une couverture subcontinue de dépôts lœssiques de la dernière glaciation (Weichselien). Au sein de ces dépôts, le Pléniglaciaire supérieur est représenté par un ensemble de lœss calcaires atteignant localement jusqu’à 6 m d’épaisseur, et incluant plusieurs horizons repères périglaciaires présents au niveau régional. Sur la base de corrélations stratigraphiques de proche en proche d’environ une centaine de profils, étudiés notamment à l’occasion de sondages et de fouilles archéologiques, une séquence synthétique comportant un maximum de trois unités principales séparées par des niveaux repères de type gley de toundra cryoturbé ou géliflué a été définie puis progressivement complétée (Antoine, 1990 ; Antoine et al., 1998, 2003) (fig. 1).

2Au Pléniglaciaire supérieur deux domaines malaco-biogéographiques, nommés "ouest" et "central", ont été différenciés en Europe de l’ouest sur la base des malacofaunes des dépôts lœssiques (Moine, 2008). D’après la valence écologique de ses taxons, le domaine "ouest", qui s’étend sur les régions périphériques à la Manche dont la moitié ouest du nord de la France, est plus humide et à végétation pauvre. Le domaine "central", qui s’étend du Bassin Parisien à l’est de l’Allemagne, est quant à lui plus sec et à végétation plus diversifiée (fig. 2).

3D’après les datations numériques disponibles, les dépôts lœssiques du Pléniglaciaire supérieur se sont mis en place entre environ 35 ka BP (âge moyen du dernier sol brun arctique/boréal du Pléniglaciaire moyen weichselien) et environ 17 ka BP (âge moyen d’arrêt de la sédimentation lœssique en Europe) (fig. 1) (Antoine et al., 2001 ; Vandenberghe et al., 2004). Entre ces deux âges limites, la plupart des enregistrements malacologiques sont discontinus, et restent mal datés (absence de datations numériques). La reconstitution précise de l’évolution paléoenvironnementale spatiale et temporelle n’est donc actuellement pas possible pour cette période.

4Ainsi, dans le domaine "central", seuls les enregistrements malacologiques de Hainerberg et Gräselberg (Wiesbaden, Allemagne) (Remy, 1969), Schwalbenberg (Remagen, Allemagne) (Schiermeyer, 2002) et Nussloch (Allemagne) (Moine et al., 2008) ont fait l’objet d’un échantillonnage selon une colonne continue avec une résolution suffisante pour aborder la dynamique environnementale à l’échelle millénaire. Dans le domaine "ouest", la plupart des enregistrements malacologiques incluent moins de cinq échantillons prélevés isolément, et au mieux dans les différentes unités stratigraphiques (Moine, 2008). A Halling (Royaume-Uni) (Kerney, 1971) et à Maisières-Canal (Belgique) (de Coninck, 1973), 25 et 12 échantillons ont bien été prélevés respectivement en colonne continue, et en colonne discontinue. Toutefois, les alternances sédimentaires de ces deux sites ne sont pas typiques des dépôts lœssiques régionaux en raison de conditions de dépôt particulières (fond de vallée à Maisières-Canal et dépôts de versants crayeux à matrice lœssique à Halling).

5La dynamique malacologique et environnementale millénaire du Pléniglaciaire supérieur weichselien reste donc mal connue, notamment dans le domaine "ouest". Les nouveaux enregistrements à haute résolution réalisés dans le cadre de sondages liés à des opérations d’archéologie préventive ont donc fait l’objet de l’étude qui va suivre.

2 – Contexte sédimentaire et paléoenvironnemental

6Les dépôts lœssiques du Pléniglaciaire supérieur présentent les plus forts taux de sédimentation de la dernière glaciation (Frechen et al., 2003) et sont formés d’alternances de niveaux de lœss et de gleys de toundra visibles dans la plupart des profils du nord-ouest de l’Europe continentale (Semmel, 1968 ; Bibus, 1980 ; Haesaerts, 1980 ; Antoine, 1988 ; Huijzer, 1993 ; Frechen, 1999 ; Frechen et al., 1999 ; Schirmer, 2000 ; Antoine, 2002 ; Vandenberghe et al., 2004 ; Antoine et al., 2009). Dans la majorité des sites étudiés, deux horizons repères sont presque toujours présents. A la base des séquences, le « sol de Saint-Acheul / Villiers-Adam », aussi appelé « sol des Vaux » en Belgique et « Lohner Boden » en Allemagne, représente le dernier sol brun arctique/boréal du Pléniglaciaire moyen weichselien (fig. 1). Le terme « complexe de sols de Saint-Acheul / Villiers-Adam » est utilisé dans les profils où ce sol ne peut pas être distingué des autres horizons du Pléniglaciaire moyen. Vers le sommet des séquences, « l’horizon à langues de Nagelbeek », aussi appelé « Erbenheimer Boden E4 » en Allemagne, est le dernier gley de toundra précédant les dépôts de lœss homogène souvent épais du Dernier Maximum Glaciaire (lœss Brabantien ; Haesaerts et al., 1981). Entre ces deux niveaux repères, plusieurs horizons du type gley de toundra alternent avec des niveaux de lœss lité à microfentes de cryodessication liés à une dynamique nivéo-éolienne (Antoine et al., 2002) et (ou) localement, dans les coupes situées à proximité d’une importante source de sables, à un litage purement éolien comme à Nussloch en rive droite de la vallée du Rhin (Antoine et al., 2009).

7La dénomination "gley de toundra" désigne des horizons hydromorphes résultant de la saturation en eau de la couche active d’un pergélisol riche en glace au cours de périodes de dégel estival. Ils sont caractérisés par une teinte générale grisâtre à verdâtre due à la réduction du fer dans le sédiment, ainsi que par la présence de bandes irrégulières et de tâches de couleur rouille résultant de processus d’oxydation localisés au niveau des traces de racines et de bioturbations diverses (tubules à auréoles oxydées de quelques millimètres de diamètre). On note également la présence d’un léger enrichissement en carbone organique dû à un enracinement plus dense, une activité biologique plus intense et un léger enrichissement en argile (+ 4 à 6 %) par rapport aux lœss encaissants (Antoine et al., 2009). Enfin, une faible décalcification peut redistribuer une partie des carbonates à la base de l’horizon sous la forme de petites concrétions calcaires (poupées du lœss) (Huijzer, 1993 ; Antoine et al., 2001).

8Au cours des multiples alternances gel-dégel qui affectent ces horizons, des processus de cryoturbation (cryoinjection) affectent les gleys de toundra formés en contexte de plateau ou de faible pente et produisent des petits polygones de toundra (maille 20-30 cm) bien visibles en décapage horizontal. Lors des phases de dégradation particulièrement marquées du pergélisol, les involutions subissent des processus de fluage lent (cryoreptation ; Van Vliet-Lanoë, 1985) qui produisent une déformation des structures dans le sens de la pente et la formation des « horizons à langues ». Dans les sites situés sur pente, ce processus débouche sur un fluage généralisé de plus grande ampleur de la couche active (gélifluxion) et génère une érosion des niveaux sous-jacents.

9Dans les séquences du nord de la France, de grandes fentes en coin de type "ice-wedge", de 1 à 2 mètres de profondeur pour 0,5 à 0,6 mètres d’ouverture, sont systématiquement associées aux horizons de gley de toundra (Antoine, 1990 ; 1991). Le fluage de la partie supérieure de la couche active du gley de toundra à l’intérieur de la fente en cours de dégradation (fonte de la glace du coin) démontre une relation intime entre le fonctionnement des réseaux de grande fente et des gleys de toundra, comme c’est actuellement le cas dans les environnements à sols polygonaux (French, 2007), ou comme Haesaerts et Van Vliet-Lanoë (1974) l’ont observé dans les dépôts lœssiques de Belgique.

10Les alternances lœss-gley peuvent être très semblables entre les différents profils d’une même région, comme par exemple dans la vallée du Rhin (Semmel, 1968 ; Bibus, 1980) ou dans le nord de la France (Antoine, 1991). Les niveaux de lœss et les gleys de toundra ont alors les mêmes positions stratigraphiques et parfois des épaisseurs et des signatures sédimentaires très similaires. Une telle homogénéité des signatures stratigraphiques implique que la formation des gleys de toundra dans les dépôts lœssiques résulte de changements climatiques d’envergure continentale et en aucun cas de facteurs locaux.

11Dans le nord de la France comme dans la vallée du Rhin, la présence systématique de faciès lités nivéo-éoliens à microfentes dans la partie inférieure du Pléniglaciaire supérieur confirme la mise en évidence d’une tonalité globalement plus humide (rôle de la couverture nivale) qu’au moment de la mise en place des lœss calcaires homogènes postérieurs à l’horizon à langues de Nagelbeek qui traduisent un environnement particulièrement sec et aride. La partie inférieure des dépôts de lœss lité nivéo-éolien du Pléniglaciaire supérieur inclut généralement plusieurs gleys de toundra épais (~0,5 m) et cryoturbés. Ils sont beaucoup plus ténus et faiblement déformés dans la partie supérieure, excepté l’horizon à langues de Nagelbeek, qui présente un faciès typique à langues grisâtres humifères dans une zone côtière allant de la Normandie à la Belgique, et qui peut-être considéré comme le plus récent des gleys de toundra pléniglaciaires (Haesaerts, 1985 ; Van Vliet-Lanoë, 1992 ; Huijzer, 1993 ; Vandenberghe et al., 1998).

12Parallèlement à cette règle générale, il existe une variabilité importante concernant la puissance des dépôts lœssiques du Pléniglaciaire supérieur, ainsi que le nombre et l’aspect des gleys de toundra. Par exemple, on a d’une part le cas unique de Nussloch, où un maximum de huit gleys de toundra (dont cinq principaux et deux biphasés) ont été observée dans une séquence de plus de treize mètres d’épaisseur (Antoine et al., 2009), et d’autre part de nombreux profils ne montrant pas plus de trois ou quatre gleys de toundra dans une séquence dépassant rarement trois ou quatre mètres d’épaisseur (Antoine, 1990, 1991). Malgré un forçage climatique d’envergure continentale, à l’origine de la formation des gleys de toundra et du dépôt des différentes unités de lœss, des différences significatives ont été notées entre les séquences. Ces dernières peuvent résulter de variations importantes du taux de sédimentation en liaison avec les conditions géomorphologiques locales et la distances aux sources, et (ou) de phases d’érosion conduisant à la présence de hiatus.

13Lorsque deux phases de formation de gley de toundra ne sont séparées que par un mince dépôt de lœss, ce dernier est affecté sur toute son épaisseur par la formation du gley de toundra supérieur. Les deux gleys de toundra apparaissent alors sous la forme d’une seule unité complexe incluant deux phases de cryoturbation et de gleyification en partie superposées et visuellement indissociables. C’est le cas du gley de toundra G2 à Nussloch dans les profils P2 et P3 (Antoine et al., 2002, 2009) et très vraisemblablement dans les séquences du nord de la France au niveau de l’Horizon cryoturbé du Santerre (Antoine, 1991 ; Antoine, obs. pers.).

14Par ailleurs, des phases d’érosion ont pu conduire à la disparition de tout ou partie des sédiments préalablement déjà déposés. Des datations précises sont alors requises pour évaluer les hiatus potentiels. Dans les coupes du bassin de la Somme comme Saint-Sauflieu ou Renancourt par exemple, la sédimentation au Pléniglaciaire supérieur débute directement par des lœss lités sans aucune trace de gley de toundra (Antoine, 1990 ; Fagnart et al., sous presse). Toujours à Saint-Sauflieu, l’horizon à langues de Nagelbeek est remplacé par un niveau d’érosion souligné par un cailloutis cryoturbé dont la mise en place est liée à un contexte plus pentu et à l’affleurement d’un substrat caillouteux sur le versant surplombant le site.

15L’interprétation paléoenvironnementale des alternances lœss-gley du Pléniglaciaire supérieur a récemment progressé grâce à des analyses granulométriques (Antoine et al., 2001 ; Rousseau et al., 2002, 2007 ; Antoine et al., 2009), géochimiques (Hatté et al., 1999, 2001) et malacologiques (Moine et al., 2008) à haute résolution réalisées dans les dépôts lœssiques de Nussloch. Ces travaux ont permis de proposer des corrélations entre les alternances lœss-gley les plus marquées et les variations climatiques cycliques millénaires de la dernière glaciation, i.e. cycles de Dansgaard-Oeschger, qui ont débouché sur le modèle conceptuel suivant (fig. 3) :

16A- Stade : taux de sédimentation lœssique élevé, taille élevée des particules, malacofaunes peu abondantes – Environnement froid et sec.

17B- Transition stade-interstade : baisse du taux de sédimentation lœssique, taille décroissante des particules, mise en place d’un pergélisol, développement d’un gley de toundra, ouverture de coins de glace sous la couche active, et malacofaunes peu abondantes, accroissement des proportions en taxons hygrophiles et palustres – Environnement froid et plus humide.

18C- Interstade : taille minimale des particules, approfondissement de la couche active et du gley de toundra, et proportion maximale en taxons hygrophiles et palustres, puis fonte des coins de glace, gelifluxion et (ou) cryoturbation du gley de toundra, malacofaunes très abondantes, et proportion de juvéniles morts très élevée – Environnement humide avec une saison chaude plus douce et plus longue, et des hivers toujours rigoureux.

19D- Transition interstade-stade : dernières phases de gelifluxion et (ou) cryoturbation du gley de toundra, retour rapide à un fort taux de sédimentation lœssique, légère augmentation de la taille des particules, réduction drastique et brutale des populations de mollusques, diminution des proportions de taxons hygrophiles et palustres. – Environnement froid et sec l’hiver, et rapidement de moins en moins humide à la saison chaude qui baisse en intensité et raccourcit en durée.

20Le modèle proposé par Vandenberghe et Nugteren (2001) est proche de celui-ci et implique également l’installation d’un pergélisol en association avec l’ouverture des coins de glace et le développement des gleys de toundra. La principale différence proposée sur la base des nouvelles recherches menées à Nussloch et dans le nord de la France est un développement des gleys de toundra synchrone de la formation des grandes fentes à coins de glace et non postérieur. Les séquences observées montrent en effet un fluage systématique par gélifluxion des gleys de toundra (couche active) sur la bordure et dans les dépressions liées à la fonte des coins de glace associés. Une fois ce processus terminé, ces dépressions sont colmatées lors de la mise en place de lœss typique (fossilisation et formation des pseudomorphoses de coins de glace) dans un environnement à nouveau très froid et aride. Les analyses malacologiques à haute résolution associées aux séquences étudiées ont en outre permis de préciser la dynamique des variations d’humidité et de température associées à ces alternances lœss-gley.

21Dans le nord de la France, la séquence stratigraphique synthétique des dépôts du Pléniglaciaire supérieur est formée depuis la base, d’un gley de toundra géliflué associé à un réseau de grands coins de glace, d’un lœss calcaire homogène d’épaisseur réduite (~ 0,5 à 0,6 m) mais extrêmement constant au niveau régional, d’un gley de toundra épais et cryoturbé (Horizon cryoturbé du Santerre), d’un lœss calcaire lité à microfentes, d’un gley de toundra géliflué grisâtre faiblement humifère (horizon à langues de Nagelbeek) souvent associé à un deuxième réseau de coins de glace dégradés (gélifluxion), et d’un lœss calcaire homogène (loess de couverture) (fig. 1).

22Dans la vallée du Rhin, la base de la stratigraphie synthétique est identique avec un gley de toundra géliflué à grands coins de glace, et un lœss calcaire homogène d’épaisseur réduite. Ensuite, quatre gleys de toundra (Erbenheimer Boden E1, E2, E3 et E4), se succèdent alternant avec d’épais dépôts de lœss lité. Le plus récent de ces gleys de toundra (E4) est recouvert par un dépôt de lœss homogène comme dans le nord de la France. Le gley de toundra E2 est parfois dédoublé et des coins de glace lui sont souvent associés. Enfin, un niveau de tuf volcanique (tuf d’Eltville) est souvent présent entre les gleys de toundra E3 et E4 (fig. 1).

23L’enregistrement sédimentaire du Pléniglaciaire supérieur de la séquence de Nussloch est particulièrement complet. Grâce aux nombreuses datations par luminescence (Lang et al., 2003 ; Rousseau et al., 2007 ; Antoine et al., 2009) et ¹⁴C (Hatté et al., 2001) composant sa chronologie absolue, des corrélations précises ont été proposées entre les alternances lœss-gley de cette séquence et les variations des indicateurs climatiques globaux de la calotte de glace arctique, notamment ceux de la concentration en poussière. Toutefois, mis à part les séquences de l’est des Carpathes où de nombreuses datations ¹⁴C ont été obtenues à partir de charbons de bois provenant d’horizons essentiellement archéologiques (Haesaerts et al., 2003), bien peu de profils stratigraphiques disposent d’un tel cadre chronologique. La datation des alternances lœss-gley et leur corrélation entre profils stratigraphiques d’une même région ou de régions différentes est donc délicate sur la seule base de leur faciès et de leur position stratigraphique, notamment en raison des variations du taux de sédimentation et de la présence de hiatus d’origine érosive.

24L’objectif de cette étude est donc la caractérisation malacologique des changements paléoenvironnementaux associés aux alternances lœss-gley du Pléniglaciaire supérieur dans le nord de la France, et leur comparaison avec ceux de la vallée du Rhin. En parallèle, nous avons cherché à déterminer si les alternances lœss-gley étaient caractérisées par une dynamique des populations de mollusques terrestres propre à chacune, pouvant constituer un élément de corrélation entre profils stratigraphiques du nord-ouest de l’Europe.

3 – Nouvelles données

3.1 - Localisation des sites

25Le matériel malacologique analysé dans ce travail provient de deux séquences lœssiques situées dans le nord de la France à Curgies (Nord), environ sept kilomètres à l’est de Valenciennes) et à Bourlon (Pas-de-Calais), environ huit kilomètres à l’ouest de Cambrai (fig. 2). Les séquences sédimentaires de ces deux sites ont été observées à la faveur de sondages de diagnostics archéologiques opérés par l’Institut National de Recherches Archéologiques Préventives (INRAP).

26A Curgies, la zone sondée se situe en position d’interfluve sur les plateaux les plus septentrionaux du Bassin Parisien, en rive droite de l’Escaut. Le sondage SP2 a fait l’objet d’un levé stratigraphique, de prélèvements malacologiques et de trois échantillons destinés à réaliser des datations par luminescence optiquement stimulée (OSL) (fig. 4). Les résultats des datations n’étaient toutefois pas encore disponibles au moment de la rédaction de cet article.

27A Bourlon, la zone sondée se situe en contexte de plateau, le cours d’eau le plus proche étant le ruisseau de l’Agache. Deux sondages (ZP7 et ZP 8) ont fait l’objet d’un levé stratigraphique et de prélèvements malacologiques (fig. 4). Les prélèvements malacologiques T1 à T4 du profil ZP7 ont été replacés sur le profil ZP8. Les prélèvements malacologiques réalisés à Bourlon n’ayant au départ qu’une valeur de test, aucune datation n’a été envisagée pour ce site.

3.2 - Chronostratigraphie

28L’interprétation chronostratigraphique des profils de Curgies et de Bourlon s’appuie sur la synthèse des travaux réalisés dans les dépôts lœssiques du nord de la France au cours des vingt dernières années (Antoine et al., 1998, 1999, 2003). En l’absence de datations absolues, elle est principalement fondée sur l’analyse de la succession des unités litho-pédo-stratigraphiques des horizons périglaciaires, leur hiérarchisation et leur intégration dans le canevas stratigraphique régional (fig. 4).

29Dans le profil de Curgies, le limon brun orangé de base (1) a été interprété comme un horizon Bt de sol brun lessivé interglaciaire rapporté à l’Eemien, alors que le limon brun orangé sus-jacent (2), interprété comme un sol forestier, représente quant à lui un bilan condensé de plusieurs pédogénèses du Début-Glaciaire weichselien. Le Début-Glaciaire weichselien est surtout représenté par le limon gris (3), interprété comme un sol steppique, et attribué à l’un des derniers épisodes de pédogenèse de cette longue période de transition. La plupart des artéfacts archéologiques découverts à Curgies proviennent de la partie supérieure de cet horizon humifère. Ce pédocomplexe a été observé à des profondeurs très variables dans les différents sondages réalisés, témoignant d’une paléotopographie nettement différente de l’actuelle. Le niveau repère suivant est le limon argileux brun-orangé feuilleté (9) interprété comme le sol de Saint-Acheul / Villiers-Adam ou la partie supérieure du complexe de sols de Saint-Acheul / Villiers-Adam. L’ensemble de limons jaune et de gleys de toundra calcaires compris entre le complexe de sols du Début-Glaciaire et le complexe de sols de Saint-Acheul / Villiers-Adam (unités 4 à 8) serait donc d’âge Weichselien inférieur ou moyen. Pour confirmer son âge une datation OSL est en cours. L’ensemble de limons jaune et de gleys de toundra calcaires sus-jacents au complexe de sols de Saint-Acheul / Villiers-Adam (unités 10 à 12) est donc logiquement attribuable au Pléniglaciaire weichselien supérieur. Des datations OSL du sol de Saint-Acheul / Villiers-Adam (9) et du limon jaune calcaire (11) remplissant la pseudomorphose de coin de glace sont en cours pour préciser leur âge. Bien qu’étant directement sus-jacent au complexe de sols de Saint-Acheul / Villiers-Adam, aucun élément ne permet de préciser si le gley de toundra de base est corrélé avec le gley de toundra de base (A) ou avec l’horizon cryoturbé du Santerre (B) dans la stratigraphie synthétique des dépôts lœssiques du nord de la France (fig. 4). L’horizon géliflué qui termine cet ensemble (12) est interprété comme étant l’équivalent local de « l’horizon à langues de Nagelbeek ».

30A Bourlon, seule la partie supérieure des sondages ZP7 et ZP8 a été explorée en détail. Le limon brun argileux feuilleté (1), épais de plusieurs mètres, et au sein duquel aucune subdivision n’a été mise en évidence, a été mis en parallèle avec le complexe de sols de Saint-Acheul / Villiers-Adam. L’ensemble lœssique calcaire sus-jacent (unités 2 à 5) est donc logiquement d’âge Weichselien supérieur. Compte tenu du faible développement vertical de la séquence de Bourlon, il est possible de corréler le gley de toundra (2) observé à la base des lœss calcaires soit avec le gley de toundra de base (hypothèse A), soit avec l’horizon cryoturbé du Santerre (hypothèse B) dans la stratigraphie synthétique régionale (fig. 4). En raison de sa position stratigraphique, le gley de toundra supérieur (4) est très certainement l’expression locale de l’horizon à langues de Nagelbeek.

3.3 - Prélèvement et analyse des échantillons malacologiques

31A Curgies, des échantillons de 5 litres de sédiment ont été prélevés en continu dans deux colonnes distinctes (fig. 4). Au-dessous le complexe de sols de Saint-Acheul / Villiers-Adam, la première colonne s’étend de la base du gley de toundra inférieur au sommet du limon jaune supérieur. De la base vers le sommet, les échantillons ont des épaisseurs de respectivement 8, 25 et 25 cm. Au-dessus du sol de Saint-Acheul / Villiers-Adam, la seconde colonne s’étend de la base du gley de toundra épais à coins de glace presque jusqu’à la limite de décarbonatation située dans le limon jaune sus-jacent. Tous les échantillons correspondent à des tranches de sédiment d’une épaisseur de 10 cm.

32A Bourlon, des échantillons de 10 litres de sédiment et de 10 cm d’épaisseur, excepté T4 (20 cm), ont été prélevés ponctuellement (T1 à T4) (fig. 4) à titre de tests dans les différentes structures pédo-sédimentaires visibles dans le sondage ZP7, à savoir le limon brun feuilleté, le gley de toundra de base, le remplissage lœssique d’une grande fente à coin de glace, et le gley de toundra supérieur à petites fentes de gel (réseau polygonal décimétrique). Les échantillons 1 à 7 ont été prélevés en colonne continue dans le sondage ZP8 depuis le sommet du complexe de sols de Saint-Acheul / Villiers-Adam jusqu’à la limite de décarbonatation qui s’enfonçait plus profondément dans ce deuxième sondage.

33Au laboratoire, les échantillons prélevés ont été tamisés à l’eau sur un tamis de 425 μm, triés, et dénombrés selon la méthode de comptage de Ložek (1964) (tab. 1). Le volume des échantillons prélevés à Curgies est inférieur de moitié à celui préconisé par Ložek (1964) et Puisségur (1976) pour les lœss. Les valeurs d’abondance obtenues pour ce site n’ont toutefois pas été doublées pour conserver l’intégrité des données de comptage. La composition des assemblages ne montre en outre aucun changement aberrant suggérant une mauvaise représentativité de la malacofaune dans ces échantillons. Les exigences écologiques des taxons ont été compilées d’après divers ouvrages de référence (tab. 2). Pour l’interprétation des données malacologiques, nous avons aussi utilisé les spectres écologiques de Ložek (1964) et de Puisségur (1976), et les indices de diversité et d’équitabilité de Simpson (1949) et de Shannon (Shannon & Weaver, 1949). Le rapport Juvéniles/Adultes (J/A) proposé par Moine et al. (2005) a été calculé pour les populations des taxons du genre Pupilla. Les courbes obtenues pour P. muscorum et P. alpicola montrant les mêmes variations, leurs effectifs ont été fusionnés pour le calcul de cet indice. En théorie, lors d’une période climatique plus favorable, l’abondance augmente indépendamment de la proportion en individus juvéniles via le renforcement du cycle de reproduction. Or, en période glaciaire, des valeurs très élevées du rapport J/A peuvent accompagner les augmentations d’abondance. La mortalité constatée de nombreux individus juvéniles nés lors des périodes d’amélioration climatique rapides intra-glaciaires peut s’expliquer par une plus grande proportion d’éclosions trop tardives ou trop précoces durant la saison estivale. Les conséquences en sont respectivement une maturité physiologique insuffisante pour permettre l’entrée en hibernation (Moine et al., 2005) ou un effet létal de périodes de gel printanier tardives postérieures au début de la période d’activité printanières des individus juvéniles.

3.4 - Malaco-biogéographie

34La composition relativement pauvre de leurs malacofaunes des profils de Curgies et de Bourlon (tab. 1) confirme leur appartenance au domaine malaco-biogéographique "ouest" au Pléniglaciaire supérieur (fig. 2). L’enregistrement malacologique de Curgies est le premier à inclure la totalité des taxons caractérisant le domaine "ouest" dans plusieurs échantillons d’une même séquence. Cette caractéristique peut être due à la proximité de la limite avec le domaine "central" où la malacofaune est plus diversifiée, ou à des conditions micro-environnementales favorables à la diversification de la malacofaune. En outre, on notera la présence dans les deux sites de Pupilla alpicola, identifiée pour la première fois dans des dépôts lœssiques du Pléniglaciaire supérieur du nord de la France (fig. 5). Cette espèce caractéristique des environnements palustres (tab. 2) est parfois difficile à différencier de Pupilla muscorum qui est plutôt une espèce polyvalente de milieux ouverts secs ou humides. En raison de la grande variabilité morphologique de Pupilla muscorum résultant de sa plasticité écophénotypique (Rousseau & Laurin, 1984 ; Rousseau, 1997), ses morphotypes de milieux très froids et très humides à la spire haute et à l’ouverture arrondie et peu épaissie peuvent être très semblables à ceux de Pupilla alpicola. Une révision du matériel malacologique du domaine "ouest" en grande partie antérieur à ces analyses morphométriques permettrait de rechercher la présence de Pupilla alpicola et donc l’occurrence de phases palustres dans les séquences lœssiques concernées.

35

4 – Résultats

4.1 - Définition des malacozones de Curgies et de Bourlon

36Les deux enregistrements malacologiques (fig. 6 & 7, tab. 1) ont été chacun divisés en plusieurs malacozones sur la base des combinaisons de variations les plus significatives et les plus diagnostiques de différents paramètres des malacofaunes (tab. 3).

4.2 - Interprétations paléoenvironnementales de la séquence de Curgies

37MZ1 « Zone à T. hispidus » : une aridité importante et un couvert végétal très peu développé (très faible nombre de granules de vers de terre) peuvent expliquer la quasi-absence de mollusques terrestres et la faible épaisseur des deux gleys de toundra. Aucun indice ne suggère que les coquilles aient pu être détruites postérieurement à leur enfouissement.

38MZ2 « Zone à T. hispidus et limaces » : dans la moitié inférieure du gley de toundra, les proportions des taxons et des spectres écologiques sont relativement équilibrées, et suggèrent la présence d’une micro-mosaïque environnementale à l’aspect ouvert, très humide et peu végétalisé. Le milieu devient plus humide (proportion croissante du taxon hygrophile Succinella oblonga et apparition en faibles effectifs du taxon palustre Oxyloma elegans et de Columella columella, espèce typique des stations de montagne humides à couvert végétal), mais les températures notamment estivales restent froides (apparition du taxon boréo-alpin Columella columella, et valeurs faibles du rapport J/A et de l’abondance). Ces conditions ont conduit à la formation de grandes fentes en coin en liaison avec un pergélisol riche en glace.

39MZ3 « Zone à S. oblonga et O. elegans » : dans la moitié supérieure du gley de toundra, l’humidification du milieu s’accélère au détriment de la végétation (doublement des effectifs de l’espèce hygrophile Succinella oblonga, disparition des limaces, de Columella columella et de Trochulus hispidus, taxons nécessitant un couvert végétal plus développé), mais la micro-mosaïque environnementale perdure (maintien simultané de Pupilla muscorum et des taxons palustres Pupilla alpicola et Oxyloma elegans), soit dans l’espace avec des flaques d’eau éparses persistantes alternant avec des points hauts plus secs, soit dans le temps avec une succession saisonnière d’environnements humides puis secs permettant le développement successif des populations des espèces présentes. L’occurrence d’une succession environnementale temporelle a notamment été suggérée par Davies (1996) pour expliquer la présence de Pupilla muscorum au sein d’assemblages malacologiques dominés par des cortèges d’espèces aquatiques, palustres et de milieux humides dans des dépôts alluviaux holocènes. Au sommet du gley de toundra, l’humidité est maximale et la végétation très appauvrie : on observe une stagnation des effectifs de P. muscorum et P. alpicola, et une explosion démographique des populations hygrophiles et palustres de S. oblonga et O. elegans. Par ailleurs l’augmentation importante du rapport J/A et des abondances de S. oblonga et O. elegans indique un radoucissement sensible des températures. Dans le cas de Nussloch, la position de la séquence en sommet de greda (ride de lœss allongée) a favorisé le drainage des eaux de fonte, permettant à tous les taxons présents d’accroître leur population sous l’effet de l’amélioration climatique. A Curgies, un environnement très humide résulte du mauvais drainage des eaux de fonte, et l’intensification du cycle de reproduction ne profite qu’aux taxons hygrophiles et palustres se contentant d’un couvert végétal pauvre. Une fois passé l’optimum thermique, les phases de regel hivernal provoquent la cryoturbation du gley de toundra par mise en pression de la couche active entre la table du pergélisol et la surface à nouveau gelée du sol (Van Vliet-Lanoë, 1985, 1987).

40MZ4 « Zone à O. elegans » : cette dernière malacozone est caractérisée par le retour brutal d’une sédimentation lœssique typique, par la fin des processus de cryoturbation du gley de toundra, et par des températures estivales froides (chute drastique de l’abondance de tous les taxons). Les transitions gley-lœss s’accompagnent aussi généralement d’un milieu plus sec (augmentation des proportions en taxons de milieu ouvert et en individus mésophiles). Dans le cas présent, le milieu reste très humide et peu favorable à la diversification végétale (maintien de l’espèce palustre Oxyloma elegans, disparition de l’espèce palustre Pupilla alpicola, disparition progressive de Succinella oblonga, présence sporadique de Pupilla muscorum et Trochulus hispidus, et richesse et diversité minimales), probablement en raison d’une dépression générée par la fonte du coin de glace qui a favorisée le maintien d’eau à cet endroit. Cette situation n’est pas unique, il arrive parfois de trouver des taxons aquatiques d’eau douce en grande quantité dans des sédiments lœssiques (Mazenot, 1951 ; Moine et al., 2005) qui témoignent d’environnements inondés plusieurs mois par an. Des malacofaunes similaires, bien qu’un peu plus riches, ont été observées dans les dépôts du Dryas Récent du nord de la France et ont conduit à une interprétation similaire (Limondin & Antoine, 2001).

41Les périodes de transition gley-lœss sont souvent mal enregistrées dans les sédiments en raison de leur courte durée et des perturbations stratigraphiques générées par les phénomènes de cryoturbation (présence d’une zone de mélange de 15 à 20 centimètres d’épaisseur). En outre, leur caractérisation précise nécessiterait de séparer les lœss piégés entre les involutions des gleys de toundra dans les zones cryoturbées, ce qui est irréalisable compte tenu du volume de prélèvement nécessaire aux analyses malacologiques. De même qu’à Nussloch, cette transition semble plus rapide que la transition inverse lœss-gley.

4.3 - Interprétations paléoenvironnementales de la séquence de Bourlon

42MZ1 « Zone à P. muscorum » : le complexe de sols de Saint-Acheul / Villiers-Adam semble caractérisé par des conditions environnementales sèches (prédominance nette de l’espèce de milieu ouvert Pupilla muscorum), une végétation peu diversifiée (pauvreté de la malacofaune) et des températures très froides (abondances faibles). Ces conditions environnementales peuvent toutefois caractériser la période de mise en place du dépôt lœssique plutôt que sa transformation ultérieure en limon brun. Le sommet du complexe de sols de Saint-Acheul / Villiers-Adam n’est donc pas toujours totalement décalcifié dans le nord de la France. De nouveaux échantillonnages en colonne continue permettront de mieux connaître la composition des malacofaunes et l’évolution paléoenvironnementale au Pléniglaciaire moyen weichselien dans cette région.

43MZ2 « Zone à P. muscorum et S. oblonga » : un réchauffement important intervient à la base (apparente) du gley de toundra (abondances maximales de la majorité des taxons, valeurs constantes et peu élevées du rapport J/A). La mosaïque environnementale est aussi humide qu’en surface d’un gley de toundra (augmentation des proportions de l’espèce hygrophile Succinella oblonga).

44MZ3 « Zone à S. oblonga et P. muscorum » : les températures diminuent mais restent douces (forte baisse des effectifs qui restent néanmoins élevés dans ce contexte, valeurs toujours constantes du rapport J/A). Le milieu est un peu moins humide (légère baisse des proportions en taxons hygrophiles et palustres). Cette phase post-réchauffement est plutôt longue et bien enregistrée contrairement à celle associée à la transition gley-lœss à Curgies. Les proportions légèrement plus élevées de Pupilla muscorum et des limaces au détriment de celles de Succinella oblonga dans l’échantillon T2 par rapport aux échantillons 5 et 4 suggèrent des conditions un peu moins humides peut-être dues à une faible surélévation ou au drainage des eaux de surface aux abords du coin de glace en activité.

45MZ4 « Zone à P. alpicola » : un second réchauffement, moins intense, intervient au sommet du gley de toundra (accroissement soudain de la population de l’espèce palustre Pupilla alpicola). Contrairement à celui de MZ2, il se caractérise par un allongement de la saison de reproduction avec persistance d’hivers toujours rigoureux et/ou de gelées printanières intenses (forte augmentation du rapport J/A indiquant une surmortalité printanière ou automnale importante des juvéniles). L’humidité augmente à nouveau (baisse accrue des effectifs de Pupilla muscorum et des limaces, et arrêt de la baisse des effectifs de Succinella oblonga plus hygrophile) probablement en raison de la fonte saisonnière plus prononcée du pergélisol, ayant également entraîné la cryoturbation du gley de toundra.

46Cet épais gley de toundra pourrait-il être biphasé ? Cette situation a déjà été observée auparavant. D’une part, dans la séquence de lœss de Savy (Aisne, France), un gley de toundra, situé dans la même position stratigraphique que celui de Bourlon, apparaît constitué de deux unités distinctes séparées par un niveau de pseudomorphoses de grandes fentes à coin de glace colmatées par des lœss homogènes (Locht et al., 2006). D’autre part, dans la séquence de lœss de Nussloch, un épais gley de toundra (G2) observé dans les premiers profils analysés (P2 et P3) (Antoine et al., 2001) s’est finalement révélé être constitué de deux unités distinctes (G2a et G2b) dans le dernier profil (P4) mis à jour (Antoine et al., 2002). La surimposition de deux gleys de toundra distincts à Bourlon est donc tout à fait plausible dans un contexte local à faible taux de sédimentation.

47MZ5 « Zone à P. muscorum » : la reprise de la sédimentation lœssique s’accompagne de conditions beaucoup plus froides (réduction drastique des populations, rapport J/A très faible), plus sèches (disparition de l’espèce palustre Pupilla alpicola et faibles proportions en taxons hygrophiles) et d’une légère diversification végétale (réapparition sporadique de Trochulus hispidus). Le coin de glace ayant fondu, la grande fente a pu localement piéger ou retenir plus d’humidité ou d’eau expliquant ainsi les 26 % de Succinella oblonga de l’assemblage T3 dans ce contexte lœssique. Le matériel malacologique du genre Pupilla de cette unité lœssique se caractérise par un déficit en apex, qui est la partie la plus fine, légère et fragile de la coquille (valeurs de plus en plus négatives du rapport Juvéniles/Adultes). La même observation a déjà été faite dans les niveaux les moins abondants des lœss de Nussloch (Moine et al., 2008). En l’absence de traces de destruction physique ou chimique sur les coquilles récoltées, ce manque d’apex peut résulter de la très faible concentration en matériel dans les échantillons, ou de pertes durant les phases d’échantillonnage et de tamisage en raison de cette très faible concentration. La composition similaire des échantillons malacologiques 2 et T3 suggère un remplissage rapide de la grande fente en coin ainsi que leur contemporanéité. Enfin, au niveau de l’échantillon T1 intervient une augmentation d’humidité (formation du gley de toundra à polygones décimétriques) accompagnée d’un appauvrissement de la végétation (absence de Trochulus hispidus, comme dans l’échantillon de base de cette malacozone). Aucun réchauffement ne semble associé à ce gley de toundra peu épais (absence d’augmentations importantes de l’abondance et du rapport J/A) dont la composition de la malacofaune est plus proche de celle du lœss sous-jacent que du gley de toundra épais.

5 - Discussion

5.1 – Implications paléoclimatiques

48La dynamique régionale de formation des gleys de toundra en Europe de l’ouest et le modèle conceptuel établi à Nussloch (voir Contexte et fig. 3) nous permettent, en l’absence de datations absolues pour les profils de Curgies et de Bourlon, de proposer des corrélations entre les changements paléoenvironnementaux déduits de l’analyse des malacofaunes et des sédiments, et les différentes phases des cycles stades-interstades de Dansgaard-Oeschger (tab. 4) à l’instar de ce qui a été proposé pour la séquence malacologique de Nussloch.

49Dans les deux sites étudiés, seule la séquence de Curgies présente l’enregistrement d’une transition stade-interstade (Zone à T. hispidus et limaces). A Nussloch, ce type de phase de transition est associé au développement d’un gley de toundra (Moine et al., 2008) par engorgement de la couche active suite à la formation d’une nappe perchée au sommet du pergélisol (Haesaerts & Van Vliet-Lanoë, 1981). Le gel saisonnier accompagné de la formation de lentilles de glace cisaillant les racines dans le sol aurait conduit à l’appauvrissement de la végétation. Ces phénomènes seraient favorisés par des épisodes d’affaiblissement des hautes pressions de la cellule polaire générant une pénétration plus fréquente des dépressions en provenance de l’Atlantique Nord et une augmentation des précipitations (Hatté & Guiot, 2005), très probablement sous la forme de neige.

50Durant les interstades (MZ3 à Curgies ; MZ2 à Bourlon), l’arrivée d’un peu de chaleur d’origine océanique sur le continent provoque une augmentation de l’intensité de la fonte saisonnière du sommet du pergélisol et débouche sur la fusion des coins de glace (approfondissement de la couche active, épaississement de l’horizon de gley dans les dépôts sous-jacents, fluage du gley sur les bords des grandes fentes après la fonte des coins de glace). Cet apport de chaleur est à relier à des périodes d’affaiblissement maximal des hautes pressions de la cellule polaire qui permettent le passage plus fréquent de dépressions atlantiques sur l’Europe de l’ouest. Toutefois, cette tendance s’inverse rapidement, et la fin des phases interstadiaires se caractérise déjà par un début de baisse des températures, notamment estivales, et le ralentissement des processus de cryoturbation et (ou) de gélifluxion des gleys de toundra (Antoine et al., 2009).

51A Curgies comme à Bourlon, le brusque retour de la sédimentation lœssique et de conditions très froides et sèches au sommet de l’épais gley de toundra central suggèrent un rééquilibrage tout aussi brutal des masses d’air. La brièveté de cette phase de transition interstade-stade ne permet pas de lui associer un enregistrement sédimentaire distinct, excepté peut-être les dépôts de lœss piégés entre les involutions qui affectent le sommet des gleys de toundra cryoturbés. Les lœss et les malacofaunes des malacozones supérieures de chaque site sont quant à eux typiques d’une phase stadiaire.

52Notons enfin deux événements particuliers. Premièrement, MZ3 est une malacozone particulière présente uniquement à Bourlon. Malgré des effectifs élevés dans ce contexte glaciaire, la forte décroissance des abondances et les faibles valeurs du rapport J/A ne permettent pas d’attribuer cette malacozone à une phase strictement interstadiaire comme MZ2 sous-jacente. Une malacozone aux caractéristiques équivalentes a également été distinguée à Nussloch au sein du gley de toundra G2. Deuxièmement, l’unique échantillon T1 prélevé à Bourlon dans le gley de toundra supérieur à polygones décimétriques ne permet pas d’aborder la dynamique de sa formation. Toutefois, l’absence d’augmentation significative de l’abondance en mollusques qui caractérise ce niveau a déjà été observée dans les gleys de toundra peu épais faiblement ou non cryoturbés de la moitié supérieure des dépôts lœssiques du Pléniglaciaire supérieur de Nussloch, et notamment dans le gley de toundra G8 (équivalent local de l’Horizon à langues de Nagelbeek ; Moine et al., 2008). Pourtant, la seule présence des polygones de taille décimétrique et des petites fentes associées, implique au moins une faible augmentation de l’humidité même si la malacofaune ne la reflète pas. Ces phénomènes résulteraient donc d’une phase de transition stade-interstade n’ayant pas débouché sur des conditions climatiques pleinement interstadiaires. Les caractéristiques des gleys de toundra de la seconde moitié des dépôts lœssiques du Pléniglaciaire supérieur de Nussloch (faible épaisseur, quasi-absence de cryoturbations, absence de maxima d’abondance en mollusques), incluant notamment l’équivalent local de l’horizon à langues de Nagelbeek (G8), ont conduit à la même interprétation (Moine et al., 2008).

5.2 - Implications stratigraphiques

53Les enregistrements malacologiques des gleys de toundra de Curgies et de Bourlon montrent une dynamique différente (aucun indice de bi-partition du gley de toundra à Curgies) (fig. 6 & 7). Malgré leur position stratigraphique similaire, il est probable que leur formation ne soit pas contemporaine et qu’elle résulte donc de cycles climatiques stade-interstade différents. Pour tester cette hypothèse, nous les avons comparés avec les enregistrements malacologiques des gleys de toundra de Nussloch en adaptant l’épaisseur des pics d’abondance de Curgies et de Bourlon à ceux de Nussloch (fig. 8).

54La structure bi-phasée de l’enregistrement malacologique du gley de toundra de Bourlon permettait d’envisager une corrélation soit avec le gley de toundra G2, soit avec les gleys de toundra G3 et G4 de Nussloch. Après comparaison des données, la première possibilité s’est révélée la plus cohérente (fig. 8) L’analyse malacologique du gley de toundra G2 a été réalisée dans le profil P3 de Nussloch où les subdivisions stratigraphiques de cet horizon n’étaient pas lisibles (Moine et al., 2008). La comparaison des enregistrements malacologiques des gleys de toundra de Nussloch (G2) et de Bourlon montre un profil d’abondance similaire : maxima à la base, décroissance importante et palier au sommet. Dans les deux sites, le second réchauffement apparaît de moindre intensité, et bénéficie surtout au taxon hygrophile Succinella oblonga. Les deux enregistrements malacologiques ont aussi en commun des valeurs de richesse relativement stables et un rapport J/A plus élevé durant le second réchauffement que durant le premier. Seules les valeurs des indices de diversité et d’équitabilité sont un peu plus variables à Nussloch, sûrement en raison d’une malacofaune plus riche qui inclut des taxons plus sensibles aux changements climato-environnementaux que ceux de Bourlon, qui sont présents dans quasiment toutes les malacofaunes glaciaires européennes (Ložek, 1964).

55La structure uni-événementielle de l’enregistrement malacologique du gley de toundra de Curgies permet d’envisager trois corrélations possibles avec d’une part les gleys de toundra G3 ou G4, et d’autre part le pic d’abondance G2’ de Nussloch (fig. 8). Ce dernier n’est toutefois pas associé à un gley de toundra clairement exprimé (Moine et al., 2008 ; Antoine et al., 2009). L’abondance maximale à Curgies n’est que deux fois moindre que pour G2’, et les deux pics sont surtout dus à des effectifs croissants des taxons hygrophiles et palustres. Or, une fois adaptée l’épaisseur du pic d’abondance de Curgies à celui de G2’, seules les courbes du rapport J/A montrent une dynamique similaire. La mise en parallèle des deux maximas d’abondance avec ou sans adaptation de l’épaisseur du pic d’abondance de Curgies (non figuré ; léger décalage vers le bas des courbes de Curgies) montre beaucoup plus de cohérence entre les variations des différents paramètres dans les deux sites. On constate que la faible épaisseur du pic d’abondance G2’ à Nussloch, composé de seulement deux échantillons, rend difficile toute comparaison avec d’autres enregistrements malacologiques. Une éventuelle corrélation entre les enregistrements malacologiques du gley de toundra de Curgies et du pic d’abondance G2’ n’est donc pas facile à mettre en évidence, mais elle n’est pas pour autant exclue.

56La comparaison des paramètres malacologiques de Curgies avec les enregistrements malacologiques des gleys de toundra G3 et G4 est beaucoup plus probante (fig. 8). Même si les pics d’abondances à Nussloch sont environ cinq fois plus importants, ils ont une allure plus semblable à celui de Curgies, notamment celui de G4 plus asymétrique. En outre, la richesse et les indices de diversité montrent tous une forte décroissance, les rapports J/A présentent un maxima synchrone, et les indices d’équitabilité montrent une baisse notable suivie d’une augmentation égale dans G3 et un peu plus importante dans G4. Les principales différences sont l’absence d’une augmentation de la proportion du taxon hygrophile Succinella oblonga durant G3, et la désynchronisation de ce signal entre Curgies et Nussloch dans G4. Cette désynchronisation est due à une croissance plus prédominante des effectifs des taxons hygrophiles et palustres à Curgies qu’à Nussloch au sommet de G4. Ceci peut s’expliquer par le contexte topographique plus plat à Curgies qui permet le développement de zones humides, alors qu’à Nussloch la position en sommet de dune de lœss implique des environnements mieux drainés moins favorables aux mollusques hygrophiles. La corrélation du gley de toundra de Curgies avec les gleys de toundra G3 ou G4 de Nussloch apparaît donc plus probable qu’avec le niveau associé au pic d’abondance G2’. Toutefois, même si la dynamique de l’ensemble des indices semble un peu plus favorable à G4, un contrôle chronologique indépendant reste néanmoins absolument nécessaire pour trancher entre ces trois possibilités.

57Cette tentative inédite de corrélation d’enregistrements malacologiques entre les deux domaines "ouest" et "central" suggère que des malacofaunes contemporaines mais distantes peuvent être caractérisées par une dynamique similaire à l’échelle millénaire. Ceci a été montré également par les comparaisons des malacofaunes tardiglaciaires de séries d’Europe de l’ouest et centrale (Royaume-Uni, France, Pays-Bas, Allemagne) (Limondin-Lozouet, 2002). Seule la multiplication de tels enregistrements malacologiques associée à des datations absolues permettra de renforcer cette hypothèse. En cas de confirmation, l’approche spatiale à haute résolution temporelle de l’évolution des paléoenvironnements d’Europe de l’ouest au Pléniglaciaire supérieur, qui constitue à l’heure actuelle l’objectif de ces analyses malacologiques, serait donc envisageable. Cette tentative de corrélation inédite a également des implications stratigraphiques. En effet, malgré leur position stratigraphique similaire, les gleys de toundra de Curgies et de Bourlon ne seraient pas tous deux corrélables avec le gley de toundra de base de la stratigraphie synthétique des dépôts lœssiques du nord de la France si l’on tient compte de leurs malacofaunes. Des hiatus d’âge et de durée différents affecteraient donc les enregistrements sédimentaires du Pléniglaciaire supérieur du nord de la France, notamment dans sa moitié inférieure, ce qui expliquerait la relative simplicité de ces enregistrements sédimentaires par comparaison avec ceux de la vallée du Rhin (voir Contexte). La corrélation entre les gleys de toundra épais échantillonnés en colonnes continues à Curgies et à Bourlon avec la stratigraphie synthétique régionale n’est donc pas évidente. En outre, l’absence d’âge pour l’horizon cryoturbé du Santerre ne permet pas d’établir une corrélation avec les gleys de toundra de Curgies et de Bourlon malgré les corrélations proposées avec les gleys de toundra datés de Nussloch. Si les datations réalisées à Curgies étayent les corrélations établies sur la base de la malacofaune, alors le gley de toundra de Curgies serait plus récent que celui de Bourlon. L’âge absolu de l’horizon cryoturbé du Santerre n’étant pas connu, deux hypothèses de corrélation sont possibles (fig. 4) : (1) le gley de toundra de Curgies serait corrélable avec l’Horizon cryoturbé du Santerre (B), et celui de Bourlon au gley de toundra de base de la stratigraphie synthétique régionale (A) ou à un autre gley de toundra intermédiaire, ou (2) le gley de toundra de Bourlon serait corrélé à l’horizon cryoturbé du Santerre (B) et celui de Curgies à un autre gley de toundra sus-jacent non-figuré. Les datations OSL en cours pour Curgies et celles associées aux analyses malacologiques à venir devraient permettre à l’avenir de lever ces questionnements.

6 – Conclusion

58L’analyse de ces deux premiers enregistrements malacologiques à haute résolution dans le domaine malaco-biogéographique "ouest" a permis malgré le faible nombre d’échantillons composant chaque séquence :

• de renforcer le découpage et la caractérisation biogéographique du nord-ouest de l’Europe au Pléniglaciaire supérieur.

• de décrire pour la première fois la dynamique malacologique et paléoenvironnementale associée aux alternances lœss-gley dans ce domaine biogéographique.

• de mettre en évidence la variabilité des réponses environnementales associées aux alternances lœss-gley.

• de proposer, par analogie avec les résultats précédemment obtenus à Nussloch, une première tentative de corrélation entre les alternances lœss-gley décrites dans le nord de la France et les différentes phases des oscillations climatiques rapides stade-interstade (cycles de Dansgaard-Oeschger) qui jalonnent la dernière période glaciaire.

• de proposer des corrélations entre les gleys de toundra de Curgies et de Bourlon, et certains gley de toundra de Nussloch.

59Les résultats obtenus posent toutefois de nouvelles questions concernant notamment la stratigraphie des dépôts de lœss du nord de la France (corrélation des profils, âge et durée des formations sédimentaires,…) auxquelles seule la multiplication de ce genre d’analyses malacologiques, et surtout la réalisation d’un cadre chronologique absolu, permettront d’apporter de nouveaux éléments de réponse.