1Les dynamiques de peuplement des basses plaines littorales ont été fortement influencées par les évolutions paléogéographiques qui ont affecté ces espaces au cours de la seconde moitié de l’Holocène. Dans un premier temps l’absence de sites connus a été imputée à des causes sociales (« répulsivité » des milieux, type d’organisation des territoires et/ou des réseaux de peuplement, habitat de hauteur ; Py, 1993). Ce n’est que récemment que les biais liés à la paléogéographie (colmatage récent de zones lagunaires ou palustres par exemple ou encore exploitation privilégiée de certains types de milieux) ou aux dynamiques sédimentaires (biais taphonomique de recouvrement) ont été pris en compte afin de corriger l’image du peuplement, principalement acquis par prospection de surface (Berger, 2011). L’intégration de ces données a modifié en profondeur la vision des relations entre les sociétés du passé et l’environnement.
2Dans la plaine du Roussillon (fig. 1), la compréhension de l’organisation territoriale des sites et la présence d’espaces vides de toute occupation avant la période Moderne a été l’un des motifs de la mise en place d’un Programme Collectif de Recherche (PCR) (Carozza, 2005). La faible occupation humaine pré-Moderne de la plaine, notamment à l’aval du méridien de Perpignan, malgré de nombreux diagnostics ou prospections, a fait l’objet de différentes explications. Les hypothèses déterministes liées aux milieux ont été privilégiées par l’archéologie, sans toutefois s’appuyer sur des éléments analytiques précis. Marichal et al. (1997) ont suggéré que l’érosion post-Néolithique pouvait également avoir influencé la préservation des artefacts et la répartition des sites connus. Calvet et al. (2002) ont mis en évidence l’importance de l’alluvionnement récent lié au Petit Age Glaciaire et l’effet de masque qu’il pouvait induire. Ces auteurs soulignent toutefois la forte variabilité spatiale de ce recouvrement qui régulariserait une paléotopographie plus différenciée. Celle-ci pourrait être pour partie héritée d’une phase d’incision médiévale entre le xe et le milieu du xiiie siècle (Carozza et al., 2010), qui aurait pu conduire au déstockage d’une partie des sédiments accumulés dans la plaine. Enfin, la préservation différentielle liée au mode d’accumulation et d’empilement des dépôts pourrait également avoir joué un rôle important (Carozza, 2011). Ces différentes hypothèses ne sont pas nécessairement contradictoires et rendent avant tout compte d’échelles d’observation différentes.
3Dans ce travail, nous montrons que la répartition des sites archéologiques à l’échelle de la plaine est avant tout contrôlée par l’architecture sédimentaire de la plaine et par le mode d’édification des deltas à haute énergie, plus particulièrement dans le cas de la basse plaine de l’Agly.
4Afin de reconstituer l’évolution de la basse plaine de la Salanque, une démarche multi-sources a été mise en œuvre par la combinaison de 3 méthodes :
5- une approche de la géométrie des dépôts à petite échelle, réalisée par dépouillement systématique de la Banque de Données du Sous-Sol (BDGSS) du BRGM. Les données recueillies sont principalement les faciès sédimentaires décrits sommairement sur les logs de terrain, plus rarement des données de diagraphie. La chronostratigraphie de ces formations reste sommaire et distingue généralement le substratum pliocène, les terrains pléistocènes et holocènes. Pour l’Holocène, une distinction entre Holocène marin (sableux ou argileux) et continental est parfois proposée, basée sur la révision de Duvail (2004). Toutefois, certaines attributions restent marquées par une forte incertitude et des confusions entre formations marines, lagunaires et paléosols continentaux épais. Les données chrono-stratigraphiques sont absentes.
6- une approche géoarchéologique à partir du suivi d’opérations d’archéologie préventive sur des grands linéaires (RD83 notamment) ou d’interventions ponctuelles (Torreilles) et de carottages mécaniques réalisés à proximité de sites connus. Les carottages ont été effectués à l’aide d’une plateforme APAGEO. Dans les deux cas, l’approche se base sur un relevé systématique de logs lithostratigraphiques et l’identification des principaux faciès pédosédimentaires in situ, incluant les contenus biologiques et archéologiques. Les principales phases d’arrêt ou de ralentissement de la sédimentation enregistrées par des paléosols associés ou non à des occupations ont fait l’objet d’une attention particulière. Un ensemble de datations, issues soit des données archéologiques soit de datations radiocarbones réalisées sur charbons pour le domaine continental ou sur matière végétale pour les formations laguno-marines a également été réalisé afin d’améliorer la résolution chronologique (tab. 1).
7- enfin une approche historique, par le dépouillement systématique des données sur la localisation des cours d’eau a été entreprise. Cette approche permet de couvrir les périodes les plus récentes (huit derniers siècles) et notamment d’intégrer l’évolution enregistrée au cours de la transition entre l’Anomalie Climatique Médiévale et le Petit Age Glaciaire (fin xiie - début xive siècle). Il s’agit principalement de repérer dans les sources les mentions de « confonts » (limites en vis à vis) entre des propriétés et un cours d’eau qui donne un terminus post quem pour l’abandon des chenaux, ou bien les mentions de chenaux anciens (Agly Vieil ou Vell) qui fournissent alors un terminus ante quem de l’arrêt du fonctionnement. Ici, la résolution chronologique peut atteindre une précision inférieure à la dizaine d’année, parfois même annuelle (Puig & Carozza, 2012).
8L’ensemble de ces informations a ensuite été mis en séquences à différentes échelles spatiales.
Fig. 1 : Localisation de la zone d’étude
9La plaine du Roussillon (fig. 1) correspond à une zone de fort élargissement des basses plaines littorales en relation avec le développement d’un bassin sédimentaire néogène de forme trapézoïdale (Calvet, 1994 ; Carozza, 1998). La plaine est organisée en deux ensembles distincts, séparés par les hautes collines de Perpignan à Canet, qui délimitent la Salanque au nord et la plaine d’Elne au sud. La Salanque est principalement drainée par l’Agly et la Têt, qui constituent les deux hydrosystèmes les plus importants du Roussillon. Du fait d’une alimentation par des bassins versants morphologiquement contrastés, le régime des deux cours d’eau diffère significativement. Ils présentent toutefois comme caractéristique commune de conserver une pente aval forte et un débit à plein bord élevé qui en font des systèmes à haute énergie.
10D’un point de vue fonctionnel, la plaine a été divisée en deux entités (Benech & Mussot, 1995). Dans la partie amont, la plaine fonctionne « à l’érosion » (i.e. les formes alluviales y sont étagées). A l’aval, elle fonctionne « à l’accumulation » (i.e. les formes alluviales aggradent ou progradent). La transition entre les deux domaines est matérialisée par l’intersection entre les profils en long de la terrasse pléistocène la plus récente et des basses plaines holocènes (Carozza, 2011). La terrasse pléistocène passe au-dessous des formations holocènes. Ce point de segmentation correspond généralement à l’apex du système de cônes deltaïques holocènes. De ce fait, les formations holocènes de la Salanque sont limitées au nord par les formations de la Crau de Rivesaltes attribuées au Pléistocène supérieur (Fy1, Clauzon et al., 1989). Au centre de la Salanque, l’interfluve Agly-Têt est bien matérialisé par les terrasses du Pléistocène moyen à supérieur (Fw et Fy) de Saint-Estève et Pia. Vers l’aval, cet interfluve s’abaisse et la limite entre les deux systèmes est alors construite et difficile à identifier d’un point de vue topographique. Par contre, la signature pétrographique des apports de l’Agly par les « sables noirs » remaniés des marnes albiennes, constitue un critère fiable de délimitation de l’influence entre les deux systèmes.
11Le maximum transgressif holocène reste localisé en aval du point de segmentation et sa position dépend principalement de l’importance de l’incision pléniglaciaire. En amont, se développe la plaine alluviale, dominée par les processus fluviaux. Vers l’aval se développent les systèmes progradants et aggradants des basses plaines sensu stricto.
12La mobilité des cours d’eau dans la plaine de la Salanque joue un rôle majeur dans l’édification de la plaine et est liée au maintien d’une puissance spécifique forte y compris à l’aval des systèmes (> 40 w.m-2). Elle a fait l’objet de nombreuses hypothèses, notamment en ce qui concerne l’Agly. Sur la base de données de photo-interprétation, Guy (in Buscail-Martin, 1978) a suggéré l’existence d’un tracé orienté sud-nord débouchant au niveau de l’anse de la Roquette. Ce tracé aurait été en relation avec l’occupation protohistorique du Port à Salses-le-Château. L’attribution de ce tracé à l’Holocène pose problème (Carozza et al., 2010) et il pourrait correspondre à une position tardiglaciaire en relation avec la migration de l’Agly jusqu’à son axe holocène. Plus récemment Serrat (2000) a proposé un modèle d’évolution de la Salanque en deux phases. L’Agly aurait d’abord édifié un lobe septentrional et se serait jeté dans l’étang de Leucate jusqu’au ive siècle (lobe 1) avant de prendre une direction proche de l’actuel (lobe 2). Ce schéma partiel s’appuie sur peu d’éléments de datation directs.
Tab. 1 : Datations 14C.
13Deux opérations séparées d’environ 500 m et localisées à l’est de Torreilles ont permis d’appréhender les dynamiques sédimentaires au sud de la zone d’étude (fig. 2).
14Sur le site de Gorg d’en Guillem (Corban, 2007), un sondage profond permis de recouper la séquence holocène sur une épaisseur de 6 m (fig. 3). Elle débute par 2,5 m d’argile sombre qui n’ont pu être observés que de manière indirecte. Les niveaux inférieurs ont livré des fragments de roseaux qui ont fait l’objet d’une datation radiocarbone. Au-dessus, le niveau compris entre 3,50 et 2,90 m montre un enrichissement en sables grossiers et renferme un niveau à malacofaune lagunaire. Cette séquence est couronnée par un niveau argileux sombre à structure prismatique qui a livré des fragments de céramique non tournée. Cette unité est localement ravinée par des apports sableux à sablo-graveleux grossiers, parfois légèrement chenalisés. Cette séquence grano-décroissante se termine par un niveau argilo-limoneux sombre qui renferme des fragments de céramique épars. Ce niveau est recouvert par un limon massif de près de 0,60 m qui a livré à son sommet des artefacts antiques sensu lato. Le sommet de la séquence entre 1,50 m de profondeur et la surface montre un retour à des conditions de dépôt plus dynamiques, dominées par les sables et les sables grossiers. Cette dernière unité renferme un paléosol peu marqué au niveau duquel s’ouvrent des fosses de plantation (paléosol 4 ; fig. 3).
15Les travaux menés sur le site de La Madraguère, ont permis une observation plus précise de la séquence supérieure (fig. 4). Ils montrent une bonne corrélation avec les données de Gorg d’en Guilhem. La séquence peut être divisée en 3 unités. A la base, une formation détritique grossière incisée dans des niveaux sableux occupe la partie nord du site. Elle se raccorde latéralement à des niveaux plus fins qui ont livré des éléments de foyers attribués à la Protohistoire (Vial et al., 2009). Le sommet de la séquence n’est préservé que localement. Il a en effet été partiellement érodé par la mise en place de la deuxième unité sédimentaire qui est associée à un paléochenal et une levée de berge à intercalation de niveaux sableux grossiers, fins et limoneux très rythmés, qui témoignent de l’existence d’un chenal confiné. Le milieu de cette séquence renferme un paléosol peu marqué, associé à un niveau d’épandage de matériel antique et qui peut être corrélé au paléosol 3 de la séquence du Gorg d’en Guillem (fig. 3). La séquence supérieure correspond à un ensemble sablo-limoneux, parfois faiblement chenalisé, qui vient sceller un niveau d’ouverture fosses de plantation médiévales.
16Ces deux sites permettent de reconstituer l’évolution du secteur sud de la Salanque. Cette zone était occupée entre env. 3 000 av. J.-C. et jusque vers env. 1 000 av. J.-C par une lagune. L’émersion se produit peu avant le Bronze final III, soit vers 900-850 av. J.-C., en liaison avec une augmentation des apports détritiques. L’apparition de faciès de chenaux puis de chenaux et de levées de berge attestent de la mise en place d’un cône distal sableux. Les apports se poursuivent jusqu’au ive siècle ap. J.-C., période à laquelle la vitesse de sédimentation moyenne décroît rapidement. Ce changement de dynamique alluviale (passage d’une sédimentation de sables grossiers associée à des environnements de chenal-plaine alluviale proximale) à une sédimentation sablo-limoneuse de plaine distale est le résultat d’une défluviation de l’Agly, documentée par les travaux de la RD83.
Fig. 2 : Localisation des opérations dans la plaine de la Salanque
Fig. 3 : Coupe transversale de la tranchée profonde de Gorg d’en Guillem.
Fig. 4 : Séquence alluviale sur le site de Torreilles - La Madraguère.
17Deux opérations successives ont permis l’observation de l’organisation sédimentaire sur un linéaire total de plus de 4 km (Vignaud, 2008 ; Passarius et al., 2010). Elles ont permis l’identification de quatre séquences alluviales majeures (fig. 5). La base de la séquence est représentée par des formations lagunaires argilo-sableuses, fréquemment riches en coquilles et en fragments végétaux. Cette séquence est identifiée dans la partie orientale de la zone diagnostiquée, mais s’étend peut-être vers l’ouest où l’épaississement des formations supérieures n’a pas permis son observation. Elle est recouverte par un ensemble de formations sableuses grossières dont l’épaisseur diminue vers l’est. Elles permettent de mettre en évidence l’existence d’un chenal principal peu incisé bordé par deux bourrelets de berge, dont l’un porte un site attribué à la période wisigothique (entre le ve et le viiie siècle ap. J.-C.). L’ensemble de cette topographie est colmaté par une formation fluviale à l’ouest et lagunaire à l’est. Celle-ci permet de localiser un cordon interne à « shell hash layer » qui sépare espace fluvial et lagunaire. Cette formation est couronnée par un paléosol faiblement esquissé.
18Les données chronostratigraphiques obtenues confirment ce découpage. Elles permettent de dater les formations lagunaires entre le ier et le iiie siècle ap. J.-C. La mise en place de la deuxième séquence alluviale est antérieure à la fin du viiie siècle ap. J.-C. comme l’indiquent quatre dates radiocarbones obtenues sur des niveaux supérieurs de cette formation. Cette attribution est cohérente avec les données archéologiques et notamment l’identification au niveau de la levée de berge d’une occupation wisigothique (ve-viiie siècle ap. J.-C.). Celle-ci se localise donc au débouché du tracé alto-médiéval de l’Agly dans la lagune de Leucate (fig. 6). La troisième séquence alluviale à dominante sableuse dans la partie ouest peut être attribuée au Moyen Âge central. Elle se raccorde à l’est à un cordon lagunaire attribué à la fin du xve siècle ap. J.-C.
19L’étude complémentaire des sources historiques permet de localiser et de dater les tracés anciens de l’Agly entre la fin du xiie et la fin du xiiie siècle ap. J.-C. Elles montrent que sur cet intervalle de temps, l’Agly a adopté un tracé méridien et parcouru un cône alluvial enraciné à l’ouest de Saint-Laurent-de-la-Salanque. En particulier, le secteur qui porte le toponyme d’Agly Vell (fig. 6) et préserve dans le paysage un tracé méandriforme a pu être rattaché à cette étape de fonctionnement. Ce système a cessé de fonctionner à la fin du xiiie siècle ap. J.-C., probablement entre 1279 et 1280 comme l’atteste la première mention de chenal abandonné à cette dernière date (Puig, 2009).
20Sur la base de ces observations directes et d’informations complémentaires issues de la Banque de Données du Sous-Sol et de travaux antérieurs, il est possible de reconstituer la dynamique d’alluvionnement dans la Salanque entre la fin de l’Âge du Bronze et le début de la période Moderne.
21La zone nord de la Salanque constitue jusqu’à la fin de l’Âge du Bronze un vaste espace lagunaire qui prolonge vers le sud la lagune de Leucate. L’émersion du secteur de Torreilles, en relation avec les apports de l’Agly est clairement identifiée sur les sites de La Madraguère et de Gorg d’en Guillem. Un premier lobe alluvial commence donc à s’édifier au sud de la Salanque vers ~1 000 av. J.-C. L’occupation attribuée au Bronze final sur le site de La Madraguère s’inscrit bien dans ce schéma et est également cohérente avec la découverte d’un site de l’Âge du Bronze IIIa (1050-900 av. J.-C.) enfoui au lieu-dit Eixugador (Kotarba et al., 2008). Non loin, la présence en sub-surface d’une occupation antique datée entre le ier siècle av. J.-C et le ive siècle ap. J.-C. au lieu-dit Les Parrudes s’inscrit également dans ce schéma. L’abandon du site sur lequel des plombs de pêche ont été retrouvés pourrait être une conséquence du déplacement de l’Agly. Vers le ive siècle ap. J.-C., la vitesse de sédimentation se réduit drastiquement sur le site, en relation avec la migration de l’Agly vers le nord. Le cours d’eau débute alors l’édification d’un second lobe au niveau de Saint-Laurent-de-la-Salanque. Cette date de mise en place est confirmée par une datation radiocarbone obtenue sur des formations sableuses de prodelta dans la lagune de Leucate par Buscail-Martin (1978) entre la fin du iiie et le viie siècle ap. J.-C. (recalibrée par nos soins). Ce lobe est fonctionnel durant le Haut Moyen Âge et probablement jusqu’à la fin du ixe siècle ap. J.-C. L’occupation tardo-antique à wisigothique de Figuera Molla (ve -viiie siècle ap. J.-C.) est installée sur une levée de berge en bordure de l’exutoire de l’Agly dans la lagune de Leucate. L’habitat est recouvert par 1,30 m de sédiments attribués à la seconde partie du Moyen Age. A cette date, l’Agly débute l’édification du lobe de Saint-Hippolyte plus au nord. Les sources historiques (Puig, 2009) permettent d’identifier deux phases de construction : une phase précoce à l’est autour du lieu-dit Conangle (xe -fin xiie siècle ap. J.-C. ?) puis une phase tardive entre le xiie siècle et la fin du xiiie siècle ap. J.-C. L’abandon de ce lobe se fait aux environs de cette période et plus probablement autour de 1280 ap. J.-C. L’Agly abandonne définitivement son tracé septentrional et construit alors son lobe actuel (lobe du Barcarès). La délimitation de l’espace lagunaire fini-médiéval (RD83 tranchée 15) suggère que le colmatage de la partie sud de la lagune est tardif et synchrone des premiers effets du Petit Age Glaciaire (fin xiiie - début xive siècle ap. J.-C.) et se prolonge durant une partie du xve - xvie siècle ap. J.-C. L’endiguement de l’Agly amorcé dans la première moitié du xive siècle ap. J.-C. (Puig, 2009), a probablement favorisé la stabilisation de la zone amont du lobe du Barcarès et le colmatage rapide de la partie sud de la lagune de Leucate.
Fig. 5 : Séquence alluviale sur le diagnostic de la RD 83.
Fig. 6 : Carte des lobes alluviaux de l’Agly entre l’an 1000 av. J.-C. et l’actuel.
22L’évolution de la plaine de la Salanque apparaît, à la lumière des travaux géoarchéologiques récents, comme beaucoup plus complexe que ne le suggéraient les travaux antérieurs. L’utilisation couplée de données géomorphologiques, archéologiques et historiques permet la reconstitution sur les trois derniers millénaires des dynamiques sédimentaires et des transformations des paysages avec une précision parfois séculaire. Loin d’une simple progradation régulière des formations continentales sur le littoral et les espaces lagunaires, l’architecture de la basse plaine montre une organisation complexe qui doit être abordée à deux échelles pour rendre compte de la relation entre processus de construction de la plaine et répartition/préservation de l’information archéologique.
23La première, qui a principalement été tentée dans ce travail, est celle de l’organisation des lobes alluviaux à l’échelle de la plaine. Au cours des trois derniers millénaires, le déplacement de l’Agly est marqué par deux périodes. La première correspond à un glissement progressif du sud vers le nord, associée à l’édification successive des lobes de Torreilles, Saint-Laurent-de-la-Salanque et de Saint-Hippolyte. Cette tendance avait déjà été observée par Calvet et al. (2002) qui l’ont interprété en terme de contrôle tectonique de la faille de la Têt. Cette interprétation semble contredite par la mise en place du dernier lobe du Barcarès au cours de la deuxième phase d’évolution récente. Il semble que le facteur déterminant de l’architecture sédimentaire à cette échelle soit les processus autocycliques de compensation latérale (lateral shift, Aslan et al., 2005). La chronologie de la mise en place des lobes contrôle à cette échelle le potentiel archéologique.
24La deuxième échelle de travail est celle de l’architecture interne du lobe et de son mode d’édification. Les défluviations secondaires, mais aussi les vitesses d’aggradation sont un facteur déterminant de la préservation de l’information archéologique. Ainsi, toujours sur le lobe de Torreilles, les sites proto-historiques sont-ils enfouis sous une épaisseur comprise entre 2,40 et 0,80 m. qui correspond à deux épisodes de dépôts différents : dépôts de chenaux-plaine d’inondation proximale protohistorique et antique d’une part, et de plaine distale du bas Moyen Âge d’autre part. Si le second épisode est favorable dans son ensemble à la préservation de l’information archéologique du fait de la faible énergie du milieu, il n’en va pas de même pour les occupations proto-historiques qui semblent occuper, comme cela a déjà été observé dans le delta du Rhône par Arnaud-Fassetta (2000) ou dans le delta du Nil par Welsby et al. (2002) préférentiellement dans les zones de levées de berges. Il est alors difficile de trancher entre préservation différentielle liée à un biais taphonomique de type stratinomique (i.e. lié à l’énergie des milieux de dépôts ; Lombard, 1978) et des formes d’adaptation aux contraintes du milieu. Un schéma similaire se dessine pour les occupations du Haut Moyen Âge sur le lobe de Saint-Laurent-de-la-Salanque. Seuls des travaux complémentaires, notamment sur les formations à haute énergie pourront permettre de trancher entre les deux mécanismes.
25Dans les basses plaines littorales, sur les deltas et cônes associés à des systèmes géomorphologiques à haute énergie, les modalités de construction sont fortement influencées par les défluviations qui commandent la durée de fonctionnement des lobes élémentaires, leur architecture et leur potentiel de préservation de l’information archéologique.
26Dans la Salanque, l’Agly a construit quatre lobes principaux successifs par entre 1 000 av. J.-C. et aujourd’hui (lobes de Torreilles, de Saint-Laurent-de-la-Salanque, de Saint-Hippolyte et du Barcarès). Ils colmatent une vaste zone lagunaire qui prolonge l’étang de Leucate vers le sud et explique l’absence d’occupations anciennes reconnues en surface dans ce secteur. Le vide de la carte archéologique s’explique donc pour partie par un effet paléogéographique. Toutefois, la mise en place de ces différents lobes est systématiquement associée à des occupations qui tendent à montrer l’appropriation rapide par les sociétés du passé de ces terres nouvelles conquises par les dynamiques alluviales sur les espaces lagunaires. Ces occupations sont aujourd’hui recouvertes par les formations mises en place à partir de la seconde moitié du xiiie siècle AD qui dépassent généralement 0,60 m en moyenne et peuvent atteindre 3 m d’épaisseur. L’absence de sites de surface est donc également le résultat d’un biais taphonomique d’enfouissement.