espèces sélectionnées pour la consommation tels *Unio* et *Anodonta*. Les autres espèces issues du même milieu, mais non collectées, ont été introduites involontairement sur le site archéologique, comme les épiphytes (*Planorbis*, *Planorbarius*, *Lymnaea*) ou les epilithes (*Theodoxus*, *Dreissena*). Ces espèces sont très peu nombreuses. Il convient cependant de souligner leur présence car elles nous offrent des informations ponctuelles sur la diversité des milieux aquatiques. On note uniquement ici les espèces qui préfèrent les eaux à fort courant, bien oxygénées, avec un fond dur comme l'exige le gastéropode *Theodoxus* sp., en association avec les bivalves *Unio crassus*, taxons actuellement très rares dans le delta et présents uniquement dans des bras à fort courant. Les espèces les plus fréquentes sont celles rencontrées dans le Danube et son bras (*Unio tumidus*, *Unio pictorum*, *Anodonta* sp. *Viviparus* sp. et *Dreissena* sp.), là où le courant est modéré, le fond sableux ou même vaseux, mais également dans les lacs alimentés de manière pérenne par des cours d'eau.

Les espèces épiphytes sur végétation submergée sont également présentes à Taraschina, car cette composante du milieu est la plus répandue.

Le rapport *Unio/Anodonta* a été utilisé pour caractériser le lieu de collecte, en comparaison avec les données actuelles sur le site de Hârșova (Radu, 2011) enregistrées pour le Danube, les bras secondaires et lacs alimentés par le fleuve. Les résultats déduits de l'analyse

des divers échantillons montrent que le fleuve et ses bras actifs représentent des lieux potentiels de collecte. Aujourd'hui, dans le delta du Danube, parmi les espèces d'Unio sont uniquement représentés Unio tumidus et Unio pictorum, mais ils sont assez rares (Sárkány-Kiss, Sîrbu, 1998). Ils ne peuvent pas, aujourd'hui, soutenir le besoin d'une population comme cela était le cas durant le 5e millénaire avant notre ère. Les seuls endroits actuels riches en bivalves restent les grands lacs (Sárkány-Kiss, 1996).

Les dimensions de ces individus traduisent une consommation intense, une sur-collecte observée toute au long de l'occupation Gumelnița. Selon nous, ce phénomène a été rendu possible par la conjonction de plusieurs facteurs : l'intérêt des humains pour la consommation des bivalves ; l'existence d'un lieu de collecte assez proche du site ; le niveau peu élevé du fleuve permettant aux bivalves de se concentrer dans certains endroits.

La collecte des mollusques est une activité praticable tout au long de la saison estivale, mais elle est étroitement dépendante des fluctuations du niveau du fleuve, de l'accessibilité aux ressources et des stratégies économiques (qui concernent également les animaux) développées à chaque saison par la communauté de Taraschina.

Durant l'occupation chalcolithique, les interventions anthropiques affectent les populations de mollusques. Elles contribuent, avec d'autres activités (pâturage, déforestation), à l'anthropisation du milieu.

L'analyse des échantillons de mollusques de Taraschina, et principalement les bivalves (dominance des *Unio* et présence d'espèces marines Cardium et Hypanis), place le site chalcolithique à proximité d'un bras important du Danube, non loin d'une lagune. Ce type de paysage est encore visible sur les cartes anciennes (Antipa, 1914 ; Mihăilescu, 1989 ; Carozza et al., 2012a ; Filip, Giosan, 2014). D'autres données viennent soutenir ces conclusions, comme la présence et la fréquence élevée de poissons qui privilégient lacs et lagunes (esturgeons, sandre et gardon de la mer Noire). De surcroît, dans le spectre des mammifères, la dominance des bétails place le site sur un terrain exondé offrant des espaces de pâturages (Bălășescu, Radu, 2011). Les dates obtenues par divers carottages autour du site (Carozza et al., 2012b, 2013 et ce volume) confirment que l'occupation a été fondée sur le large plateau de læss placé aujourd'hui très au nord.

Les changements environnementaux ont des causes multiples : géologiques, climatiques, évolution du niveau de la mer Noire etc. (Panin, 2004 ;

<figure><img src="image_60.png" /><figcaption>**Figure 4**

Comparaisons entre la hauteur des valves d'Unio tumidus actuels du Danube (Fleuve 1, N = 264 ; Fleuve 2, N =152) et les valves d'Unio tumidus archéologiques de Taraschina (Z1E2, N = 47 ; Z2E2, N=170 ; Z1E3, N = 30 ; Z2E3, N=92 ; Z1E4, N = 188 ; Z2E4, N=94 ; Z1E5, N = 24 ; Z2E5, N=101).</figcaption></figure>
Giosan et al., 2006; Lericolais et al., 2009; Carozza et al. 2011c; Weninger, Harper, 2015). Ils laissent une certaine empreinte qui est enregistrée par la communauté des mollusques (réponse biotique). Celle-ci peut être parfois décelée par l'étude des assemblages issus des sites archéologiques.

Si ces données constituent un repère temporel – principalement pour l'évolution et la caractérisation du paysage local –, elles contribuent, en association avec d'autres analyses, à une compréhension globale des phénomènes qui affectent le milieu et le fonctionnement des écosystèmes.

Références bibliographiques

**Antipa 1914**: ANTIPA (G.) – Câteva probleme ştiinţifice şi economice privitoare la Delta Dunării, *Analele Academiei Române*, seria II, tom XXXVI, Memoriile Secţiunii Ştiinţifice no 6, Bucureşti, Librăriile SOCEC &Comp. şi C. SFETEA.

**Bălăşescu, Radu 2011**: BĂLĂŞESCU (A.), RADU (V.) – Paléo-économie animalière et reconstitution de l'environnement, in L. Carozza, C. Bem and C. Micu, (eds.), *Société et environnement dans la zone du Bas Danube durant le 5ᵉ millénaire avant notre ère*, Iaşi, Editura Universității «Alexandru Ioan Cuza», p. 385-408.

**Carozza et al. 2011a**: CAROZZA (L.), BEM (C.), MICU (C.) eds. – Société et environnement dans la zone du bas Danube durant le 5ᵉ millénaire avant notre ère, Iaşi, Editura Universității «Alexandru Ioan Cuza», 430 p.

**Carozza et al. 2011b**: CAROZZA (L.), MICU (C.), CAROZZA (J.-M.), BĂLĂŞESCU (A.), HAITĂ (C.), RADU (V.), BURENS (A.), FURESTIER (R.), MIHAIL (F.), AILINCĂI (S.), FLOREA (M.) – L’habitat chalcolithique de Taraschina et le peuplement ancien du delta du Danube : implications paléogéographiques, in L. Carozza, C. Bem and C. Micu, (eds.), *Société et environnement dans la zone du Bas Danube durant le 5ᵉ millénaire avant notre ère*, Iaşi, Editura Universității «Alexandru Ioan Cuza», p. 409-413.

**Carozza et al. 2011c**: CAROZZA (L.), MICU (C.), CAROZZA (J.-M.) – Préhistoire du delta du

Danube: L'habitat chalcolithique de Taraschina, in L. Carozza, C. Bem and C. Micu, (eds.), *Société et environnement dans la zone du Bas Danube durant le 5ᵉ millénaire avant notre ère*, Iaşi, Editura Universității «Alexandru Ioan Cuza», p. 215-226.

**Carozza et al. 2012a**: CAROZZA (J.-M.), CAROZZA (L.), RADU (V.), LEVÊQUE (F.), MICU (C.), OPREANU (G.), HAITA (C.), DANU (M.) – Aftermath the flooding: geomorphological evolution of the Danube delta after the black sea-mediterranean reconnection and its implication on eneolithic settlements, in P. Gästescu, W. Lewis, P. Breţcan (eds.), *Water Resources and Wetlands, Conference Proceedings, 14-16 September 2012 – Tulcea – Romania*, p. 458-465.

**Carozza et al. 2012b**: CAROZZA (J.-M.), MICU (C.), MIHAIL (F.), CAROZZA (L.) – Landscape change and archaeological settlements in the lower Danube valley and Deta from early neolithic to chalcolitic time: A review, *Quaternary International* 261, p. 21-31.

**Carozza et al. 2013**: CAROZZA (J.-M.), CAROZZA (L.), RADU (V.), LEVÊQUE (F.), MICU (C.), BURENS (A.), OPREANU (G.), HAITĂ (C.), DANU (M.) – Après le déluge : évolution géomorphologique du delta du Danube après la reconnexion mer Noire - Méditerranée et ses implications sur le peuplement Enéolithique / Chalcolithique, *Quaternaire* 24 (4), p. 503-512.

**Chapman 2010**: CHAPMAN (J.) – Houses, Households, villages, and Proto-Cities in southeastern Europe, in D. W. Anthony and J. Y. Chi (eds.),

*The Lost World of Old Europe: The Danube Valley, 5000-3500 BC*, New York and Princeton: Institute for the Study of the Ancient World at New York University and Princeton University Press, p. 74-89.

**Filip, Giosan 2014**: FILIP (F.), GIOSAN (L.) – Evolution of Chilia lobes of the Danube delta: Reorganization of deltaic processes under cultural pressures, *Anthropocene* 5, p. 65-70.

**Giosan et al 2006**: GIOSAN (L.), DONNELLY (J.P.), CONSTANTINESCU (S.), FILIP (F.), OVEJANU (I.), VESPREMEANU-STROE (A.), VESPREMEANU (E.), DULLER (G. A. T.) – Young Danube delta douments stable Black Sea level since the middle Holocene: morphodynamic, paleogeographic, and archaeological implications, *Geology* 34, p. 757-760.

**Grossu 1957**: GROSSU (A. V.) 1957 – Consideraţii paleoecologice şi biogeografice asupra moluşelor găsite în săpături arheologice si în straturi de loess, *Materiale si cercetari arheologice* IV, p. 373-376.

**Grossu 1962**: GROSSU (A. V.) – *Fauna Republicii Populare Romîne. Vol. III, Mollusca. Fascicula. 3, Bivalvia (Scoici). Bucureşti, Editura Academei Republicii Populare Romîne*, 426 p.

**Grossu 1970**: GROSSU (A. V.) – Unele observaţii asupra gasteropode-lor descoperite în straturile romanello-aziliene de la Cuina-Turcului, *Studii şi Cercetări de Istorie Veche* 21, 1, p. 45.

**Grossu 1976**: GROSSU (A. V.) –
Considerații asupra faunei malacolo-
gice din așezarea de la Ripiceni-Izvor,
Studii și Cercetări de Istorie Veche și
Arheologie 27, 1, p. 17.

Grossu 1993 : GROSSU (A. V.) –
Gasteropodele din România. Melci
marini, de uscat și apă dulce.
Compendium. București, 412 p.

Lericolais et al. 2009 : LERICOLAIS
(G.), BULOIS (C.), GILLET (H.),
GUICHARD (F.) – High frequency
sea level fluctuations recorded in the
Black Sea since the LGM, Global
and Planetary Change 66 (1-2 ),
p. 65-75.

Mărgărit, Radu 2014 : MĂRGĂRIT
(M.), RADU (V.) – The use of
autochthonous aquatic resources in
the technologies og Gumelnița
communities/Utilizarea resurselor
acvatice autohtone în tehnologiile
comunităților Gumelnița, in
M. Mărgărit, G. Le Dosseur, A.
Averbouch (eds.), An overview of the
exploitation of hard animal materials
during the Neolithic and
Chalcolithic/O rivire asupra ex-
ploatării materiilor dure animale
de-a lungul neoliticului și calcoli-
ticului, Editura Cetatea de Scaun,
Târgoviște, p. 221-240.

Mărgărit 2016 : MĂRGĂRIT (M.)
–Testing the endurance of prehistoric
adornments: Raw materials from
the aquatic environment, *Journal of
Archaeological Science* 70, p. 66-81.

Micu et al. 2011 : MICU (C.),
AILINCĂI (S.), STĂNICĂ (A.) –
L'occupation ancienne du delta du

Danube : l’apport de l’archéologie,
in L. Carozza, C. Bem and C. Micu,
(eds.), Société et environnement dans
la zone du Bas Danube durant le 5e
millénaire avant notre ère, Iași,
Editura Universității «Alexandru
Ioan Cuza», p. 235-243.

Mihăilescu 1989 : MIHĂILESCU
(N.) – The evolution of the fluvia-
tile network of the Danube Delta in
Pleistocene and Holocene, *Travaux
du Museum d’Histoire Naturelle
Grigore Antipa* XXX, p. 355-366.

Panin, 2004 : PANIN (N.) – The Danube Delta. Geomorphology and Holocene evolution: a synthesis, *Geomorphologie: Relief, Proceesus, Paysage* 9 (4), p. 247-262.

Radu 2011 : RADU (V.) – Le
*Néolithique de la Roumanie*
*Méridionale : Exploitation
des ressources aquatiques dans
les cultures Boian et Gumelnița*,
Sarrebruck, Editions Universitaires
Européennes, 439 p.

Radu et al. 2016 : RADU (V.),
POPOVICI (D. N.), CERNEA (C.),
CERNĂU (I.), BĂLĂȘESCU (A.) –
Harvesting molluscs in the Eneolithic:
a study of freshwater bivalve accu-
mulations from the tell settlements of
Bordușani-Popina and Hârșova
(Romania, 5th millennium BC), *Envi-
ronment Archaeology* 21 (4),
p. 334-350.

Radu (sous presse) : RADU
(V.) – The Lower Danube mollusk
assemblages: zooarchaeological
and paleoenvironmental data from

Eneolithic settlements Hârșova-tell,
Bordușani-Popina and Taraschina
(România, 5th millennium BC),
in proceeding paper *The
archaeology of wetlands, the
landscape, the man and
his environment: Danube Valley in
Prehistory* (2016).

Sárkány-Kiss 1996 :
SÁRKÁNY-KISS (A.) – Contribuții
la cunoașterea faunei de moluște
acvatice din Delta Dunării,
Analele Știintifice ale
Institutului Delta Dunării 1,
p. 5-8.

Sárkány-Kiss, Sîrbu 1998 :
SÁRKÁNY-KISS (A.),
SÎRBU (I.) – Contribuţii la
cunoaşterea asociaţiilor de moluşte
acvatice din lacurile: Roşca,
Rotundu, Belciug, Saraturi-
Murighiol şi Merhei
(Rezervaţia Biosferei Delta
Dunării, România), *Analele
Știintifice ale Institutului de
Cercetare și Proiectare Delta
Dunării* 6 (1), p. 63-70.

Weninger, Harper 2015 :
WENINGER (B.), HARPER (T.) –
The Geographic Corridor for Rapid
Climate Change in Southeast, in S.
Hansen, P. Raczky, A. Anders, and
A. Reingruber (eds.), *Neolithic and
Copper Age between the Carpathians
and the Aegean Sea. Chronologies
and technologies from the 6th to the
4th Millennium BCE*; International
Workshop Budapest 2012
(Archäologie in Eurasien 31), Bonn,
Habelt, p. 475-505.
<figure><img src="image_61.png" /></figure>
# Les sédiments comme révélateurs des activités anthropiques
dans le delta du Danube - phytolithes

Mihaela DANU, Jean-Michel CAROZZA et Erwan MESSAGER

## 1. Introduction

La fouille du tell chalcolithique de Taraschina questionne sur les raisons de l'implantation d'un habitat dans un milieu ouvert, formé de plateaux bordés de chenaux et de zones humides, paysages très différents de ceux de l'actuel delta du Danube (Carozza et al., 2011b ; Carozza et al., 2012). Dans un tel contexte, il est fondamental de comprendre les caractéristiques définissant le paysage végétal durant la période chalcolithique et de savoir comment ces populations ont géré leur environnement, à proximité de leur habitat. Parmi les différents proxys botaniques, les phytolithes se sont avérés être de bons marqueurs pour mettre en évidence l'histoire des formations végétales locales. Ces vestiges, d'origine végétale, ont également permis de mieux comprendre les interactions homme-environnement, en particulier d'évaluer l'intensité des activités agricoles au moment de l'occupation chalcolithique du tell de Taraschina. Les assemblages des phytolithes, issus des dépôts inférieurs et moyens de deux séquences sédimentaires carottées à proximité et sur la bordure du tell (Logs 10 et 12 ; fig. 2), ne reflètent pas des dépôts naturels de phytolithes, mais plutôt une accumulation de vestiges issus du traitement de céréales (Danu et al., 2019).

Ces premiers résultats nous ont encouragé à poursuivre l'analyse sur d'autres carottes. Les premières ont été prélevées à la périphérie du site (Log 21 et Log 22), pour vérifier si le signal enregistré dans la séquence du Log 10 était reproductible. Il s'agissait également de confirmer la nature des activités agraires réalisées par les populations chalcolithiques. Une autre séquence a été collectée dans des secteurs plus éloignés de l'habitat pour rechercher une séquence sédimentaire potentiellement moins affectée par les activités anthropiques (Log 20, fig. 2).

## 2. Le site de Taraschina
et son environnement proche

Le site de Taraschina est situé dans l'actuel delta du Danube, sur la commune de Maliuc, à environ 23 km en amont de l'embouchure de Sulina (fig. 1). La nature exceptionnelle du site de Taraschina est essentiellement due à sa position dans le delta. Le site est attribué à la culture chalcolithique Gumelnița, sur la base de l'étude céramologique. Il est daté du 5e millénaire BC par une série de datations radiocarbone.

<figure><img src="image_62.png" /><figcaption>**Figure 1**
Localisation du site de Taraschina (DAO : L. Carozza).</figcaption></figure>
La plus ancienne occupation découverte en contexte alluvial dans le delta du Danube avait été attribuée à une période de transition entre le Chalcolithique et l'âge du Bronze (Vasiliu, 1995).

La chronologie de l'occupation du site, déduite de l'analyse des niveaux du tell, couvre une période allant de 4800 à 4100 cal. BC (Carozza et al., 2011a, 2011b). Ces occupations sont principalement liées à la culture Gumelnița. Les indices archéologiques indiquent que le site a été progressivement abandonné, vers 4150-4050 cal. BC, probablement en raison d'un changement de fonction (habitat temporaire ou spécialisé) et de modifications de l'environnement. Cependant, le haut de la séquence montre que

les occupations du tell se poursuivent entre 4300 et 4050 cal. BC. Pour l'instant, aucune occupation chalcolithique postérieure à 4100-4050 cal. BC n'est avérée. Cette date pourrait être considérée comme celle de l'abandon définitif du site. Cette hypothèse est renforcée par les assemblages de phytolithes enregistrés dans le Log 10 (Danu et al., 2019). En effet, les phytolithes caractérisant les activités agricoles sont très rares (ou absents) dans le dernier mètre de la séquence.

La cartographie du tell de Taraschina montre que la zone d'habitat s'étend sur plus de 15 000 m² ; ce dernier est composé d'au moins deux principaux sites (les sites 1 et 2, fig. 2). Chacun a livré les restes de bâtiments d'habitation effondrés, associés à des fosses.

<figure><img src="image_63.png" /><figcaption>Figure 2</figcaption></figure>

Localisation des carottes 20, 21 et 22 dans la zone du site de Taraschina (DAO : L. Carozza).
Sur un plan géomorphologique, cette partie du delta correspond à la partie supérieure et ancienne du système deltaïque. Elle est aussi nommée le « delta fluvial interne » (Panin, 1972) jusqu'aux rides de Letea-Caraorman, la structure géomorphologique la plus élevée du delta.

Le delta de fond de baie forme une aire de forme triangulaire, délimitée au nord par le bras de Chilia et au sud par celui de Sulina-Saint Georges. Dans la zone de Taraschina, le nombre de lacs est réduit, en comparaison avec d'autres zones du delta. Les carottages entrepris dans le voisinage du site montrent que la terrasse de læss apparaît sous la surface actuelle de la plaine d'inondation, entre 1 et 4 m. La séquence alluviale holocène, qui recouvre partiellement le site, se compose essentiellement de limons et d'argiles organiques. Les sédiments plus grossiers (sables fins et grossiers) sont rares, ils ont été mis au jour dans la partie nord du site. La reconstitution paléogéographique montre que le site est probablement situé en bordure d'une pente d'environ 2 m (terrasse alluviale ou incision fluviale), surplombant une zone de marécage (fig. 2).

# 3. Matériaux et méthode

## 3.1. Carottage, description des carottes et datation

Dans cette étude, nous présenterons trois séquences : les Logs 20, 21 et 22 (fig. 2). Ces carottes ont été extraites par l'équipe de fouille, à l'extérieur du site, à l'aide d'un carottier à percussion de type Cobra TT. Les séquences ont été corrélées par un certain nombre de couches sédimentaires repères et à l'aide de datations radiocarbone.

### *Log 20 :*

La séquence, longue de 3,00 m a atteint le substratum lœssique à la profondeur de 2,25 m. Cette séquence n'a pas fait l'objet de travaux analytiques de laboratoire et son interprétation repose uniquement sur la caractérisation macroscopique des faciès sédimentaires (texture, couleur, structures sédimentaires et inclusions biologiques et anthropiques). Elle comporte les unités suivantes :

- l'unité 1 correspond au substratum pré-deltaïque constitué de læss de couleur jaune pâle. Il apparaît vers -1,0 m bsl ;

- l'unité 2 est un ensemble stratifié dominé par des limons grossiers et des limons sableux. Le milieu de la séquence, vers 0,50 m bsl, montre une intercalation

sableuse qui traduit des apports de chenaux. Elle a probablement été déposée en bordure de chenal, soit sur une berge, soit en arrière d'une levée de berge ;

- l'unité 3, épaisse de 1,20 m, est constituée d'une alternance de limons fins et d'argiles limono-organiques. Elle traduit une diminution des apports détritiques et un renforcement de la composante organique de la sédimentation, probablement en lien avec l'éloignement du chenal.

### *Log 21 :*

La séquence 21 a été prélevée au sud-ouest du tell, à environ 100 m de celui-ci. La séquence, d'une longueur de 4 m, atteint le substratum lœssique à 2.90 m. Cette séquence n'a pas fait l'objet d'analyses sédimentologiques et son interprétation repose uniquement sur la description des faciès sédimentaires et sur 2 datations ¹⁴C qui y ont été réalisées (fig. 3). Elle se compose des unités suivantes :

- l'unité 1, à la base, correspond au substratum lœssique pré-deltaïque qui a été recoupé sur près d'un mètre d'épaisseur. L'altitude de l'apparition des læss est donc voisine de -1,50 m bsl ;

- l'unité 2 est une séquence limono-organique, présentant localement des enrichissements en sables. Sa base est très riche en éléments coquilliers. Elle a probablement été mise en place dans un environnement de plaine marécageuse distale ;

- l'unité 3 débute par un niveau plus sableux puis évolue rapidement vers un milieu dominé par la sédimentation organo-minérale. Il s'agit probablement d'un marais distal ;

- l'unité 4 marque un retour à des conditions plus connectées au chenal et à une augmentation de la fraction détritique fine, argileuse et limoneuse fine.

### Interprétation générale de la séquence 21

Cette séquence semble se mettre en place dans un contexte d'aggradation rapide, en bordure de chenal, comme l'attestent les récurrences sablo-limoneuses intercalées dans le remplissage organo-minéral.

### Chronologie

La chronologie de la séquence ne repose que sur 2 dates 14C (fig. 3). La datation de la base de la séquence indique un début d'accrétion verticale vers cca. 4450-4310 BC. La séquence est probablement discontinue et marquée par un hiatus et/ou un très fort ralentissement des apports au sommet de l'unité 2. La reprise de la sédimentation, qui débute par un épisode grossier, n'est pas directement renseignée.
Elle est antérieure à la date de 1430-1625 AD donnée par la datation 14C supérieure. Ici aussi, près de la moitié de la séquence est donc historique et vient re-couvrir, après un long hiatus, une séquence ancienne.

*Log 22 :*

Le carottage a traversé une séquence sédimentaire de 3,00 m d'épaisseur et a atteint le *substratum* pré-deltaïque à seulement 1,75 m sous la surface. Là également, la description de cette séquence repose uniquement sur des éléments visuels macroscopiques. Elle se décompose suivant l'organisation suivante :

- l'unité 1 est associée au *substratum* pré-deltaïque lœssique. Il a été recoupé sur une épaisseur de 1,10 m. Le sommet de cette unité se trouve autour de 0,25 - 0,30 m bsl, ce qui atteste de l'existence d'un paléo-relief ;

- l'unité 2 est dominée par des formations limono-organiques très rythmées qui attestent d'une mise en place dans une zone marginale, probablement en bordure de la plaine alluviale. La concentration en malacofaune est importante ;

- l'unité 3 est très homogène et se caractérise par une matrice limoneuse. Là aussi, la présence importante de malacofaune est un trait caractéristique. Cette unité pourrait avoir été mis en place en contexte de plaine alluviale distale soumise à des crues débordantes récurrentes et chargées en limons en suspension.

## 3.2. Échantillonnage, extraction des phytolithes, comptages et classification

### 3.2.1. Préparation et protocole de laboratoire

En vue d'une analyse des phytolithes, 7 échantillons ont été prélevés dans le Log 21. Huit autres échantillons ont été prélevés dans les Logs 20 et 22. L'extraction du contenu phytolithique a été réalisée au CEPAM Nice (UMR 7264), à partir d'un échantillon de sédiment d'environ 1 g, selon le protocole de routine du laboratoire : décarbonatation des sédiments ; défloculation des argiles, filtration, décantations

répétées pour éliminer les plus petites particules (argiles), oxydation de la matière organique, dessiccation du culot, séparation par densimétrie (à l'aide du polytungstate de sodium de densité calibrée à 2,35), rinçage et filtration des culots, séchage et montage entre lame et lamelle.

### 3.2.2. Identification, comptage et construction du diagramme

L'observation des lames au microscope optique à transmission, à des grossissements de 450 fois, permet d'identifier les différents morphotypes des phytolithes présents sur la lame. Les phytolithes sont très bien préservés dans les 3 séquences et, en règle générale, un effectif de 200 particules microscopiques d'opale (SiO₂, nH₂O) d'origine végétale a été compté.

La dénomination des morphotypes s'efforce de suivre les recommandations de l'ICPT (ICPT, 2019).

Les diagrammes de phytolithes ont été réalisés à l'aide du logiciel GpalWin.

### 3.2.3. Classification

Quatorze catégories de phytolithes ont été observées et donc utilisées pour essayer d'apporter de nouvelles informations sur le paléo-milieu et les activités anthropiques sur le secteur de Taraschina.

Les phytolithes de type ELONGATE ENTIRE, ELONGATE DENDRITIC, ACUTE BULBOSUS et BULLIFORM LABELLATE sont essentiellement formés dans l'épiderme des graminées (Twiss et al., 1969 ; Piperno, 1988 ; Mulholland, 1989), mais ils peuvent aussi être produits par d'autres groupes (Piperno, 1988 ; Strömberg et al., 2002). Le type ACUTE BULBOSUS 2 comprend des phytolithes atypiques de trichome, mais n'a pas de valeur taxonomique car il peut être produit tant par des graminées que des Dicotylédones.

Les morphotypes RONDEL, CRENATE, BILOBATE, POLYLOBATE et SADDLE sont produits dans des cellules épidermiques courtes de graminées. Ils peuvent être utilisés pour identifier les principales sous-familles de Poaceae répertoriées dans l'assemblage phyto-

<table><caption>**Figure 3**
Tableau des dates radiocarbone.</caption><thead><tr><td>Site</td><td>ensemble daté</td><td>matériel</td><td>Date BP</td><td>cal 68,2%</td><td>cal 95,4%</td><td>n° labo</td></tr></thead><tbody><tr><td>Taraschina carotte</td><td>Log 21 (-2,09/2,05)</td><td>charbon</td><td>405 ± 30 BP</td><td>1442-1601 AD</td><td>1433-1624 AD</td><td>Poz- 56451</td></tr><tr><td>Taraschina carotte</td><td>Log 21 (-2,80/2,82)</td><td>charbon</td><td>5545 ± 35 BP</td><td>4446-4348</td><td>4452-4341</td><td>Poz- 56452</td></tr></tbody></table>
lithique (Twiss et al., 1969 ; Fredlund et Tieszen, 1994). La forme RONDEL est fréquemment associée à la sous-famille des *Pooideae* (Mulholland, 1989 ; Piperno and Pearsall, 1998), mais on la trouve aussi dans d'autres sous-familles des *Poaceae*, comme les *Arundinoideae* (Barboni, Bremond, 2009). Le morphotype CRENATE peut être considéré comme un bon traceur de la sous-famille des *Pooideae* (Barboni et al., 2007). Le morphotype BILOBATE se retrouve principalement chez les sous-familles des *Panicoidae* et des *Arundinoideae*. Le morphotype de SADDLE est produit en grande proportion par les *Chloridoideae*.

Le type *Papillate* correspond à des papilles épidermiques silicifiées, essentiellement présentes dans les inflorescences de Poaceae (Ball et al., 2001 ; Berlin et al., 2003).

Le morphotype BLOCKY n'a pas été assigné à un groupe taxonomique dans cette étude. Dans les travaux antérieurs, cette classe a été attribuée à *Artemisia* ou à différents genres de *Pinaceae* (Pearsall and Piperno, 1993 ; Blinnikov et al., 2002).

Le type TRACHEARY est généralement attribué aux dicotylédones herbacées ou ligneuses, mais peut également être produit en petites quantités par les conifères et les fougères (Piperno, 1988 ; Runge, 1999).

La forme *silica skeleton* (phytolithes en connexion), correspondant à des fragments d'épiderme siliciés de *Poaceae*, a également été identifiée (Rosen, 1992). Elle correspond à des fragments de tiges, de feuilles et d'inflorescence de plantes de *Poaceae*. Au sein de ce groupe, deux classes ont été individualisées : 1. *silica skeleton* formé de phytolithes *ELONGATE DENDRITIC* provenant d'inflorescences, et 2. *silica skeleton* formé de phytolithes *ELONGATE ENTIRE* provenant de tiges et de feuilles.

De nombreux spicules et diatomées ont été observés, en particulier dans les échantillons supérieurs des séquences.

L'abondance relative de chaque morphotype a été calculée en pourcentage de la somme des phytolithes classifiés. Cette somme, utilisée pour le calcul des pourcentages, était basée sur le nombre total de phytolithes comptés par échantillon (plus de 200 par échantillon). Les *silica skeleton*, les diatomées et les spicules ont été exclus de la somme totale des phytolithes afin d'éviter la surreprésentation de ces groupes. La représentation de chaque groupe a été exprimée en pourcentage de la somme des phytolithes plus la somme du groupe [i.e. % des diatomées = N diatomées / (Σ phytolithes + N diatomées) x 100].

## 4. Résultats

Les phytolithes sont très bien préservés dans les 3 séquences (fig. 4-6). Les phytolithes de graminées (*Poaceae*) dominent tous les assemblages. Au sein de ce groupe, le groupe « *Pooideae* » (représenté par la classe CRENATE et une partie de la classe RONDEL) est la sous-famille la mieux enregistrée. *Pooideae* est le principal groupe de *Poaceae* dans les régions tempérées, mais nous ne pouvons pas exclure qu'une partie de la classe RONDEL puisse provenir aussi de la sous-famille des *Arundinoideae* (par exemple du roseau - Phragmites).

### 4.1. Log 21, périphérie du site

Si l'on excepte l'échantillon de surface de cette séquence, on note une bonne représentation des *ELONGATE DENDRITIC*, des phytolithes issus d'inflorescences (glumes et glumelles) de *Poaceae* (Ball et al., 2001). Une autre caractéristique intéressante de ces spectres est la présence très importante du groupe *silica skeleton*. Les épidermes de graminées (dont la plupart sont coupés) sont très bien enregistrés (20-30 %).

La forme la plus couramment rencontrée dans cette séquence est RONDEL et des pourcentages relativement élevés de CRENATE ont été enregistrés, ce qui souligne l'importance du groupe des *Pooideae*.

À 160 cm de profondeur, on observe une augmentation considérable des spicules d'éponges, une diminution de la valeur de *silica skeleton* et l'absence des phytolithes de type *ELONGATE DENDRITIC*.

L'échantillon de surface, situé à 15 cm de profondeur, montre l'abondance des spicules (70 %). On observe aussi la présence de diatomées (plus de 10 %), associées à de nombreux morphotypes de *BULLIFORM LABELLATE*. Les formes *ELONGATE DENDRITIC* et les épidermes sont absents. Cet échantillon témoigne d'une remise en eau et le retour d'une végétation hydrophile et non anthropisée.

### 4.2. Log 22 (4 échantillons), périphérie du site

Les assemblages issus de cette séquence sont assez similaires aux assemblages provenant de la séquence du Log 21. Les phytolithes RONDEL et CRENATE enregistrent dans cette séquence des valeurs très importantes, la sous-famille *Pooideae* étant donc bien représentée. On observe également une bonne représentation des *ELONGATE DENDRITIC* et on note que les épidermes de graminées (attestés par les *silica skeleton*) sont très bien représentés (plus de 25 %). En ce qui
<figure><img src="image_64.png" /><figcaption>**Figure 4**
Assemblages de phytolithes (%) à partir du Log 20.</figcaption></figure>

<figure><img src="image_65.png" /><figcaption>**Figure 5**
Assemblages de phytolithes (%) à partir du Log 21.</figcaption></figure>
concerne l'échantillon de surface, comme dans le cas du log 21, les spicules et les diatomées sont très bien enregistrés.

## 4.3. Log 20 (4 échantillons), plus éloigné du site

L'échantillon le plus profond (-180 cm) de cette séquence n'a fourni qu'un très faible effectif de phytolithes. Les morphotypes RONDEL et CRENATE n'ont pas été enregistrés dans cet échantillon. L'analyse des échantillons situés à -158 cm et -150 cm a livré des assemblages tout à fait similaires à ceux observés dans les enregistrements des Log 21 et 22. Dans le Log 20, comme dans les autres séquences, l'échantillon du sommet présente un assemblage dominé par les spicules, révélateur du changement de conditions édaphiques.

## 5. Discussion

Les assemblages phytolithiques issus des enregistrements sédimentaires collectés hors-site ont permis de produire de nouvelles données sur l'exploitation du milieu et sur la gestion des ressources végétales à Taraschina. Les recherches menées sur ces trois nouvelles séquences ont permis de retracer l'évolution du couvert végétal et de mettre en évidence une intense activité agricole. L'analyse des phytolithes a révélé des assemblages d'origine anthropique dans lesquels les graminées jouent un rôle majeur.

Les séquences des Logs 21 et 22 sont caractérisées par de pourcentages importants de phytolithes ELONGATE DENDRITIC. Comme nous l'avons vu, l'abondance de cette classe de phytolithes marque le traitement des céréales (Berlin *et al.*, 2003 ; Portillo, Albert, 2011). Ainsi, la présence de nombreux phytolithes ELONGATE DENDRITIC est interprétée comme le signe de la présence de céréales à Taraschina. L'abondance de la forme *silica skeleton* issue des inflorescences confirme l'abondance de ces dernières sur le site. Les pourcentages de *silica skeleton* provenant de tiges et de feuilles de *Poaceae* reflètent par ailleurs une accumulation de résidus de battage. Ce signal phytolithique permet de mettre en évidence la transformation des céréales dans la zone de Taraschina. La consommation de céréales a été attestée à Taraschina par des graines carbonisées (Danu *et al.*, 2019). Ces dernières se composent en particulier d'orge commune (*Hordeum vulgare*) et d'amidonnier (*Triticum dicoccum*).

Les spectres de phytolithes observés dans ces trois séquences « hors site » (caractérisées par de nombreux ELONGATE DENDRITIC et des *silica skeleton*) présentent des similitudes significatives avec les spectres provenant du Log 10 et ceux provenant des dépôts archéologiques, Log 12 (Danu *et al.*, 2019).

Les assemblages phytolithiques des Logs 21, 20 et 22 montrent, à l'instar des deux précédentes séquences étudiées (Logs 10 et 12, Danu *et al.*, 2019), que les

<figure><img src="image_66.png" /><figcaption>**Figure 6**
Assemblages de phytolithes (%) à partir du Log 22.</figcaption></figure>
populations chalcolithiques cultivaient des céréales
dans cette zone. L'hypothèse d'une culture locale est
renforcée par l'abondance des sous-produits de battage
retrouvés dans les trois séquences (essentiellement
la forte quantité de phytolithes de la classe *silica skeleton*). Cela démontre que les plants de céréales
entiers (paille comprise) ont été apportés à proximité
immédiate ou directement sur le site.

Bien que la séquence du Log 21 ait été carottée à
environ 100 m de distance du site, la teneur des sédi-
ments en phytolithes montre que l'influence humaine
est encore très importante. Ces assemblages de phy-
tolithes, retrouvés dans toutes les séquences – y com-
pris les plus éloignées du site –, portent le signal de
traitement des céréales. Ils reflètent l'envergure des
pratiques agraires des communautés chalcolithiques
dans ce secteur. Les assemblages de phytolithes
issus des séquences « hors site » pourraient refléter
des processus de battage, tandis que les assemblages
des dépôts archéologiques (Log 12, Danu et al., 2019)
montrent la pratique du vannage et du dépiquage.

Ces nouveaux résultats confirment que les assemblages
de phytolithes témoignent de la pratique, sur le site, de
différentes étapes de transformation des céréales.

L'origine de ces fortes concentrations en phytolithes
dans les séquences « hors site » reste une question dé-
battue ; plusieurs hypothèses peuvent être envisagées :

1. l'analyse de ces 5 logs révèle des activités anthro-
piques potentiellement liées à une zone de culture,
où les sols auraient été enrichis par fumure (déchets
ménagers et/ou fumier). Le fumier pourrait alors être
composé de sous-produits de transformation, utilisés
comme fourrage pour le bétail (Harvey et Fuller,
2005 ; Lancelloti, Madella, 2012) ;

2. une autre hypothèse pourrait être l'existence
d'une immense aire de battage. La prédominance
de la forme *silica skeleton* – provenant des tiges et
des feuilles – pourrait étayer ce scénario. L'hypo-
thèse d'une aire de battage, associée à l'utilisation du
 traîneau de battage (*Tribulum*), a déjà été proposée
 par Skakun sur les sites chalcolithiques de Bulgarie
 (Skakun, 1993 ; Anderson, 1998 ; Anderson *et al.*,
 2006). Plusieurs exemples de la forme *silica skeleton*,
 présentant des traces de coupe caractéristiques du
 *Tribulum* (Anderson, 2003, 2006 ; Anderson *et al.*,
 2006), ont été observés dans les assemblages des
 phytolithes de Taraschina (Logs 10 et 12, mais aussi
 20, 21 et 22) ;

3. Comme la paille légère peut être utilisée comme
matériau de construction, par exemple comme dégrais-
sant (Hillman, 1984), les phytolithes correspondants

auraient pu être libérés après la destruction et la
dissolution des matériaux de construction, et mélangés
avec les sédiments environnants. Au sein de l'habitat
de Taraschina il semble cependant acquis que la terre
à bâtir qui forme les constructions soit formée d'un
mélange de læss argileux et de fibres végétales issues
des roseaux (tiges et feuilles) ;

4. une autre hypothèse pouvant expliquer la grande
quantité de phytolithes présents dans les dépôts en-
vironnants le site serait le produit de l'accumulation
d'une grande quantité de restes de battage. Ces restes
résulteraient de la pratique intense du traitement des
céréales sur le site même et dans ses environs. Ils ont
pu être utilisés autour du site par les habitants (com-
bustible, bâtiment, artisanat, fourrage ; Anderson,
2003), ou simplement dispersés autour du site.

Les échantillons issus des sommets de séquences sé-
dimentaires des Log 21, 20 et 22 révèlent un déclin
progressif des phytolithes caractérisant le traitement
des céréales. Simultanément, les marqueurs du mi-
lieu aquatique (tels que les spicules d'éponge et/ou
les diatomées) marquent une augmentation significa-
tive, probablement causée par les variations hydro-
logiques. Cet aspect est également observé dans les
Logs 10 et 12. Ces assemblages de silice biogénique
(phytolithes + spicules d'éponge + diatomées) re-
flètent donc un changement environnemental signifi-
catif dans la région de Taraschina, marqué par la mo-
dification des conditions édaphiques autour du site et
le déclin de l'influence humaine.

6. Conclusion

Les assemblages de phytolithes constituent des
marqueurs fiables pour connaître la composition des
communautés végétales présentes sur le tell Taraschina
depuis la période chalcolithique. Grâce à ce proxy ar-
chéobotanique, il a été possible de mettre en évidence
des activités humaines dans le delta du Danube pen-
dant la Préhistoire. Les résultats des Logs 20, 21 et
22 révèlent la prépondérance des marqueurs d'un pa-
léo-milieu fortement anthropisé. La présence massive
des formes caractéristiques du traitement des céréales
dans ces 3 séquences atteste de la transformation des
plantes dans l'environnement proche du site. Cette
étude confirme que les populations chalcolithiques de
Taraschina ont cultivé des céréales au cœur du delta,
il y a plus de 6000 ans. Les données phytolithiques
indiquent ensuite un déclin progressif des marqueurs
de traitement des céréales et l'augmentation simul-
tanée des marqueurs du milieu aquatique, donc de la
montée des eaux et du changement environnemental
qui en a résulté.
Références bibliographiques

**Anderson 1998 : ANDERSON**
(P.C.) – History of Harvesting and
Threshing Techniques for Cereals in
the Prehistoric Near East, in
A.B. Damania, J. Valkoun, G.
Willcox, C.O. Qualset Eds., *The
origins of agriculture and crop
domestication (International Center
for Agricultural Research in the Dry
Areas, Aleppo, Syria)*, p. 145-159.

**Anderson 2003 : ANDERSON**
(P.C.) – Observations on the
threshing sledge and its products in
ancient and present-day Mesopotamia,
in P.C. Anderson, L.S. Cummings,
T.K. Schippers, B. Simonel Eds., *Le
traitement des récoltes : un regard
sur la diversité du Néolithique au
présent*, p. 417-438.

**Anderson 2006 : ANDERSON**
(P.C.) – Premiers Tribulums,
Premières Tractionns Animales au
Proche-Orient vers 8000-7500 BP ?
in P. Pétrequin, R.-M. Arbogast,
A.-M. Pétrequin, S. van Willigen,
M. Bailly Eds., *Premiers chariots,
premières araires. La diffusion de la
traction animale en Europe pendant
les IIe et IIIe millénaires avant notre
ère*. CNRS Editions, Monographies
du CRA 29, p. 299-316.

**Anderson et al. 2006 :** ANDERSON
(P.C.), GEORGES (J.M.), VARGIOLU
(R.), ZAHOUANI (H.) – Insights
from a tribological analysis of the
tribulum, *Journal of Archaeological
Science* 33 (11), p. 1559-1568.

**Ball et al. 2001 :** BALL (T.B.),
GARDNER (J.S.), ANDERSON (N.)
– An approach to identifying
inflorescence phytoliths from
selected species of wheat and barley,
in J.D. Meunier, F. Colin Eds.,
*Phytoliths: Applications in Earth
Sciences and Human History*. A.A.
Balkema Publishers, Exton, PA,
p. 289-302.

**Barboni et al. 2007 :** BARBONI
(D.,) BRÉMOND (L.),
BONNEFILLE (R.) – Comparative
study of modern phytolith
assemblages from inter-tropical
Africa, *Palaeogeography Palaeo-
climatology Palaeoecology* 246,
p. 454-470.

**Barboni, Bremond 2009 :** BARBONI
(D.), BREMOND (L.) – Phytoliths of
East African grasses: An assessment
of their environmental and taxonomic
significance based on floristic data,
*Review of Palaeobotany and
Palynology* 158 (1-2), p. 29-41.

**Berlin et al. 2003 :** BERLIN (A.),
BALL (T.B.), THOMPSON (R.),
KITTLESON (D.), HERBERT (S.C.)
– Ptolemaic Agriculture, “Syrian
Wheat”, and Triticum aestivum,
*Journal of Archaeological Science*
30, p. 115-121.

**Blinnikov et al. 2002 :**
BLINNIKOV (M.S.), BUSACCA
(A.), WHITLOCK (C.) –
Reconstruction of the Late Pleistocene
grassland of the Columbia basin,
Washington, USA, based on phytolith
records in læss, Palaeogeography,
Palaeoclimatology, *Palaeoecology*
177 (1-2), p. 77-101.

**Carozza et al. 2011a :** CAROZZA,
(L.), BEM (C.), MICU (C.) – Société
et environnement dans la zone du bas
Danube durant le 5ᵉ millénaire avant
notre ère, *Iași, Editura Universității
«Alexandru Ioan Cuza»*.

**Carozza et al. 2011b :** CAROZZA,
(L.), MICU (C.), CAROZZA (J.-M.),
HAITĂ (C.), BĂLĂȘESCU (A.),
RADU (V.), BURENS (A.), MIHAIL
(F.), AILINCĂI (S.), FLOREA (M.)
– L’habitat chalcolithique de
Taraschina à Mila 23 et le
peuplement ancien du delta du
Danube durant la première moitié
du 5ᵉ millenaire avant notre ère, in
S. Mills, P. Mirea Eds., *The lower
Danube in prehistory: landscape
changes and human environement
interactions - Proceedings of the
International Conference - Alexandria
3-5 November 2010, Ed.
Renaissance, București*, p. 31-47.

**Carozza et al. 2012 :** CAROZZA
(J.-M.), MICU (C.), MIHAIL (F.),
CAROZZA (L.) – Landscape change
and archaeological settlements in the
lower Danube valley and delta from
early Neolithic to Chalcolithic time:
a review, *Quaternary International*
261, p. 21-31.

**Danu et al. 2019 :** DANU (M.),
MESSAGER (E.), CAROZZA,
(J.-M.), CAROZZA (L.), BOUBY,
(L.), PHILIBERT (S.), ANDERSON
(P.), BURENS (A.), MICU (C.) –
Phytolith evidence of cereal
processing in the Danube Delta
during the Chalcolithic period,
*Quaternary International* 504,
p. 128-138.

**Fredlund et al. 1994 :** FREDLUND
(G.), TIESZEN (L.) – Modern
phytoliths assemblages from the
North American Great Plains, *Journal
of Biogeography* 21, p. 321-335.

**Harvey, Fuller 2005 :** HARVEY
(E.L.), FULLER (D.Q.) –
Investigating crop processing through
phytolith analysis: the case of rice
and millets, *Journal of Archaeological
Science* 32, p. 739-752.

**Hillman 1984 :** HILLMAN (G.) –
Interpretation of archaeological plant
remains: the application of
ethnographic models from Turkey,
in W. van Zeist, W.A. Casparie Eds.,
Plants and Ancient Man, *Studies in
Palaeoethnobotany*,
Rotterdam, Balkema, p. 1-41.

**ICPT :** Neumann (K.), Albert
(R.M.), Ball *et al.* – International
Code for Phytolith Nomenclature
(ICPN) 2.0. *Annals of Botany* 124
(2), p. 189-199.

**Lancelotti, Madella 2012 :**
LANCELOTTI (C.), MADELLA,
(M.M.) – The ‘invisible’ product:
developing markers for identifying
dung in archaeological contexts,
*Journal of Archaeological Science*
39, p. 953-963.

**Mulholland 1989 :** MULHOLLAND
(S.) – Grass opal phytolith
production: A basis for archaeological
interpretation in the northern plains.
Archaeobotany through phytolith
analysis symposium, Annual meeting
of the Soc. Amer. Archaeology,
abstracts, The Phytolitarian Newsletter
6 (1), p. 4.

**Panin 1972 :** PANIN (N.) – Histoire
quaternaire du Delta du Danube,
Essai d’interpretation des facies des
depots deltaiques, *Cercetări Marine*
4, p. 5-15.

**Pearsall, Piperno 1993 :**
PEARSALL (D.M.), PIPERNO (D.R.) – *Current research in phytolith analysis: applications in archaeology and paleoecology*, University of Pennsylvania Press, Philadelphia.

**Piperno 1988 :** PIPERNO (D.) – *Phytolith Analysis: An Archaeological and Geological Perspective*, Academic Press, San Diego.

**Piperno, Pearsall 1998 :** PIPERNO (D.R.), PEARSALL (D.M.) – The Silica Bodies of Tropical American Grasses: Morphology, Taxonomy, and Implications for Grass Systematics and Fossil Phytolith Identification, *Smithsonian contribution to botany*, 85, p. 1-22.

**Portillo, Albert. 2011 :** PORTILLO (M.), ALBERT (R.M.) – Husbandry Practices and Livestock Dung at the Numidian Site of Althiburos (El Médéina, Kef Governorate, Northern Tunisia): the Phytolith and Spherulite Evidence, *Journal of Archaeological Science* 38 (12), p. 3224-3233.

**Rosen 1992 :** ROSEN (A.M.) – Preliminary identification of silica skeletons from Near Eastern archaeological sites: An anatomical approach, in G.R. Rapp, Mulholland, S.C. Eds., *Phytolith Systematic. Plenum Press, New York*, p. 129-147.

**Rovner 1971 :** ROVNER (I.) – Potential of Opal Phytoliths for use in Paleocological reconstructions, *Quaternary Research* 1, p. 343-359.

**Runge 1999 :** RUNGE (F.) – The opal phytolith inventory of soils in

central Africa quantities, shapes,
classification and spectra, *Review of*
*Palaeobotany and Palynology* 107,
p. 23-53.

**Skakun 1993 :** SKAKUN (N.N.) – Agricultural implements in the Neolithic and Eneolithic of Bulgaria, in P. Anderson, S. Beyries, M. Otte, H. Plisson Eds., *Traces et fonction : les gestes retrouvés*. ERAUL 50, Liege, Belgium, p. 361-368.

**Twiss et al. 1969 :** TWISS (C.), SUESS (E.), SMITH (R.) – *Morphology classification of grass phytoliths. Proceedings of the Soil Science Society of America* 33, p. 109-115.

**Vasiliu 1995 :** VASILIU (I.) – *Cercetări arheologice în Delta Dunării*. Mormintele cu ocru de la Chilia-Veche, *Peuce XI*, p. 49-87.
# Une première vision des changements environnementaux dans le delta du Danube : palynologie de la carotte 14 prélevée le long du canal Șontea

Mihaela DANU et Erwan MESSAGER

Dans le secteur nord-ouest du bassin de la mer Noire, les enregistrements palynologiques et paléoécologiques proviennent essentiellement du sud de la Dobroudja, au nord-est de la Bulgarie : le complexe du lac Durankulak (Bozilova, Ivanov, 1985 ; Bozilova, Tonkov, 1985, 1998, 2002 ; Marinova, 2003 ; Marinova, Atanassova, 2006 ; Tonkov *et al.*, 2014), le lac Shabla-Ezeretz (Filipova, 1985 ; Filipova-Marinova, Bozilova, 1990), le lac Varna (Bozilova, Beug, 1994 ; Marinova *et al.*, 2016), le lac Bolata (Tonkov *et al.*, 2011), le lac Arkutino (Bozilova, Beug, 1992) et le lac Oltina (Feurdean *et al.*, 2020). Ces travaux ont permis de reconstituer l'histoire de la végétation pour ces zones de steppes et de forêts-steppes (Bohn, Golub, Hettwer 2000). L'étude de séquences marines a également livré une dynamique régionale de la végétation de toute la région côtière de Dobroudja (Popes-

cu *et al.*, 2020 ; Filipova, Marinova *et al.*, 2013). Au sein de cette région, la zone du delta du Danube constitue un îlot écologique avec des formations végétales propres au contexte deltaïque. Bien que des lacs y soient abondants, les données palynologiques holo-cènes y sont extrêmement rares, car les dépôts sédi-mentaires (lacustres, lagunaires comme fluviatiles) ne semblent pas favorables à la préservation ou à l'enregistrement du matériel sporo-pollinique (Danu *et al.*, 2019). Après de nombreuses tentatives d'ex-traction polliniques sur plusieurs carottes du delta, nous avons enfin réussi à identifier une séquence sédimen-taire (Log. 14) ayant conservé un signal pollinique. La carotte 14, prélevée dans un paléo-chenal du ca-nal de Șontea s'est avérée peu propice pour tenter une reconstruction du paléo-environnement contemporain de l'habitat chalcolithique de Taraschina (fig. 1 et 2).

<figure><img src="image_67.png" /><figcaption>**Figure 1**
Localisation du Log 14 et de l'ensemble des carottages réalisés en marge du tell de Taraschina.</figcaption></figure>
En effet, l'étude géomorphologique (voir Carozza et al., ce volume) révèle l'existence de hiatus chronologiques et de phases hydrologiques actives (dépôts de sables), peu propices à la conservation du matériel sporo-pollinique (fig. 2). Toutefois, la possibilité de disposer d'une carotte prélevée dans la zone centrale du delta du Danube et livrant du pollen nous a conduit à tenter une reconstitution des paléoenvironnements de cette zone à l'aide des pollens et des micro-fossiles non polliniques (MNP).

Ce travail permet une première description des assemblages polliniques contemporains de l'âge du Bronze. Il apporte de nouvelles données sur la connaissance des paléoenvironnements dans le delta du Danube, où peu d'informations sont disponibles.

# 1. Matériel et méthode

La carotte 14 est située 450 m à l'ouest du tell (fig. 1). Les échantillons ont été collectés avec un pas d'échantillonnage de 5 cm lorsque cela était possible. Pour chaque échantillon, un micro-prélèvement de 2 g a été soumis au protocole d'extraction des spores et des grains de pollen correspondant aux méthodes définies par Faegri, Iversen (1989), adaptées par A. Mazuy

(Assistant Ingénieur-chimiste au laboratoire CEPAM Nice). Chaque échantillon de sédiment a fait l'objet d'attaques chimiques acido-basiques (acide chlorhydrique 50 %, acide fluorhydrique 70 %, hydroxyde de potassium 10 %) afin d'éliminer l'ensemble de la matrice sédimentaire minérale et organique. Par la suite, les échantillons ont fait l'objet d'une lévigation en liqueur dense, adaptée à l'analyse pollinique des sédiments d'origine minérale, comme les sédiments archéologiques ou marins. Les préparations ont été montées entre lame et lamelle puis observées au moyen d'un microscope "Zeiss standard™" (X650, X1000) (fig. 3). La détermination des spores et des grains de pollen a été effectuée à l'aide des ouvrages et atlas de référence (Reille, 1992; Beug, 2004), ainsi que de la collection de référence palynologique.

Le pollen était bien préservé dans la partie moyenne de la carotte. En revanche, les échantillons correspondant aux niveaux les plus profonds sont stériles (ni pollen, ni microfossiles non polliniques), à l'exception d'un seul échantillon situé à 505 cm.

Les dates 14C réalisées sur cette séquence montrent que seuls les dépôts inférieurs sont contemporains de la fin de l'occupation chalcolithique de Taraschina. Or le matériel pollinique n'est malheureusement pas bien

Log14

<figure><img src="image_68.png" /><figcaption>**Figure 2**
Stratigraphie et modèle d'âge de la séquence sédimentaire LOG 14.</figcaption></figure>
préservé dans ces dépôts. Les données palynologiques nous permettent donc de documenter essentiellement l'environnement contemporain de la période qui succède à l'occupation chalcolithique du tell : l'âge du Bronze, puis de la période postérieure aux XIII-XIVe siècles en raison d'un important hiatus entre 4000 et 600 cal BP (fig. 1).

Le diagramme palynologique a été réalisé à l'aide du logiciel Tilia. Une analyse CONISS a été employée pour identifier quatre zones polliniques : PAZ (Pollen Assemblage Zones) (fig. 4).

## 2. Résultats et interprétation

La séquence est dominée par les herbacées (entre 50 et 80 %) : Chenopodiaceae, Poaceae, *Artemisia*, Asteraceae, Asteroideae, Apiaceae. Les assemblages indiquent un milieu ouvert dans lequel on observe des taux importants de Chenopodiaceae (entre 20 et 60 %). L'abondance de ces taxons est très probablement liée à la progression des espaces d'eau saumâtre. De nombreuses espèces de cette famille sont des plantes halophiles (par exemple : *Salicornia*

<figure><img src="image_69.png" /><figcaption>Figure 3</figcaption></figure>

Spores et grains de pollen extraits du Log 14 – Taraschina. a : *Cercophora* ; b : *Chenopodiaceae* ; c : *Nymphoides* ; d : *Fagus* ; e : *Pediastrum* ; f : *Pinus*.
<figure><img src="image_70.png" /><figcaption>Figure 4</figcaption></figure>

Diagramme palynologique de la séquence sédimentaire LOG 14.
sp., *Suaeda* sp., *Bassia hirsuta*, *Salsola soda* etc.). Les familles Asteraceae, Caryophyllaceae, Poaceae, Brassicaceae peuvent aussi inclure des espèces halophiles.

## 2.1. Évolution des assemblages selon les 4 zones polliniques locales (PAZ) identifiées

### PAZ1 [505 cm]

Cette zone ne correspond qu'à un seul échantillon, situé à la profondeur de -505 cm. L'assemblage est dominé par les Poacées et les Chénopodiacées (environ 40 % ensemble). Le taux d'arbre est faible (AP < 25 %). On mentionne la présence de *Quercus*, *Carpinus*, *Tilia*. Le diagramme montre aussi la présence de quelques grains de *Fagus*.

L'occurrence de spores de Glomus (HdV-207), un bon indicateur de l'érosion des sols (van Geel, 2001), accompagnées des spores de champignons coprophiles (*Cercophora* - HdV-112, *Coniochaeta cf. ligniaria* - HdV-172, *Sordaria* - HdV-55), pourrait être interprétée comme un indice d'activité anthropique. Le reste de l'assemblage ne donne aucun signe en ce sens : pas de céréales, pas de marqueurs polliniques de pastoralisme. D'après le modèle d'âge, les sédiments analysés (505 cm de profondeur) sont sans doute déjà postérieurs à l'occupation chalcolithique (postérieurs à 4100 cal BC) et ils semblent donc témoigner de l'abandon de cet espace pour l'agriculture.

On observe aussi la présence importante (20 %) de taxons aquatiques et hygrophiles (*Sparganium*, *Alisma plantago-aquatica*, *Nymphoides*). Les zygospores d'algues comme *Pediastrum*, *Mougeotia*, *Spyrogira* sont les marqueurs des variations hydrologiques et trophiques de la zone de dépôt. *Spyrogira* et *Pediastrum* peuvent indiquer dans ce cas l'élévation du niveau de l'eau, de façon durable et importante.

### PAZ2 [398 - 370 cm]

Dans cette zone, les arbres sont encore moins bien représentés (10%<AP<20%) que dans la PAZ1. *Quercus* reste le taxon forestier majoritaire. Les assemblages sont nettement dominés par les Chénopodiacées (jusqu'à 60%). Par rapport à la PAZ1, les Poacées et les Asteraceées enregistrent une hausse. Quelques grains de pollen de plantes messicoles (*Polygonum*, *Rubiaceae*) et des plantes de lieux piétinés et pâturés (*Plantago*, *Artemisia*, *Rumex*, etc.) peuvent être considérés comme de légers indices polliniques d'anthropisation. Un pic de *Cercophora* est enregistré au cours de cette phase. Il faut mentionner que ce genre regroupe plus

de 50 espèces, dont seule une moitié est coprophile et l'autre non (Doveri, 2004 ; Cugny, 2011).

On peut dès lors supposer qu'il s'agit d'une période plus favorable aux activités agro-pastorales, d'autant que les taxons hygrophiles et aquatiques sont en baisse par rapport à la PAZ1, mais les indices d'anthropisation sont ici encore très ténus. On observe aussi que *Pediastrum* et *Botryococcus* sont presque absents.

### PAZ3 [369 - 243 cm]

Le sommet de la séquence est marqué par un accroissement de la représentation des arbres (AP) et une nette diminution des Chénopodiacées. Des formations forestières sont attestées et sont mieux représentées que dans les deux premières zones (PAZ1 et PAZ2). On observe une forêt caducifoliée (chêne, charme, tilleul) et on enregistre aussi des taxons caractéristiques de la ripisylve, tels que *Alnus* et *Salix*.

Les Poacées présentent une relative stabilité. Un seul grain de pollen de *Cerealia* a été identifié. On observe également des pollens de taxons messicoles et de plantes des lieux piétinés et pâturés.

Les spores des champignons coprophiles sont présentes, mais leurs valeurs restent faibles.

### PAZ4 [244 cm]

Les pourcentages de taxons hygrophiles et aquatiques sont en hausse. Dans cette zone, on observe le signal d'une modification des conditions édaphiques, les valeurs de *Pediastrum* et *Botryococcus* étant en hausse. Bien qu'il s'agisse d'algues d'eau douce, ces deux taxons peuvent également tolérer de faibles valeurs de salinité. Cet assemblage est marqué par la hausse des conifères (*Pinus*, *Abies* et *Picea*).

## 2.2. Tentative de reconstruction de l'histoire de la végétation

Les études palynologiques antérieures montrent que, au début de l'Holocène, le sud de la Dobroudja (en Bulgarie) était recouvert d'une végétation steppique constituée de Poaceae, d'Asteraceae, de Chenopodiaceae et d'Artemisia, parsemée de groupements arborés de *Quercus*, *Carpinus*, *Corylus*, *Ulmus*, *Alnus* (Bozilova, Filipova, 1986). À partir de 8000 cal. BC, les arbres se répandent lentement et vers 6200 cal. BC, la région était couverte par une « sylvo-steppe » (Bozilova, Filipova, 1986 ; Bozilova, Tonkov, 1998 ; Filipova-Marinova, Bozilova, 1990). En ce qui concerne le delta du Danube, il est
encore difficile à ce jour de retracer l'évolution de
la végétation pour la période Chalcolithique, faute de
données polliniques contemporaines. Les assem-
blages polliniques du Log 14 mettent en évidence,
pour la période correspondant à la fin de l'occu-
pation chalcolithique, la prépondérance des Poacées
et des Chénopodiacées, ainsi que la faible présence
des arbres. Les taxons aquatiques et hygrophiles,
comme les spores d'algues, suggèrent l'élévation du
niveau de l'eau, de façon durable et importante.

Tandis que les données palynologiques attestent dans
le nord-est de la Bulgarie la pratique de l'agricultu-
re depuis le Néolithique et l'Énéolithique, dans la
séquence sédimentaire du Log 14, pour la période
Chalcolithique, il n'a pas été enregistré de grains de
pollen de Cerealia.

Pour la période correspondant à l’âge du Bronze, les données palynologiques issues de l’analyse du log 14 permettent de reconstruire la dynamique d’une végétation de sylvo-steppe dans le delta du Danube. Ainsi, plusieurs formations forestières sont attes- tées : des forêts de chênes (Quercus) accompagnées de charmes (Carpinus betulus), de charmes d’Orient (Carpinus orientalis) et de tilleuls (Tilia). On remarque également la présence du hêtre (Fagus). On note l’existence des forêts riveraines, composées d’aulnes (Alnus), de bouleau (Betula), de saule (Salix). Le nois- setier (Corylus), l’orme (Ulmus) et le frêne (Fraxinus) sont également enregistrés. Des taxons aquatiques, comme Sparganium, Potamogeton, Butomus, Alisma reflètent le caractère très faiblement saumâtre des eaux. Les signes d’une activité agro-pastorale sont très ténus dans cet enregistrement, ce qui montre que cet espace ne faisait pas l’objet d’une exploitation à cette période de l’âge du Bronze. Un seul grain de pollen de Ceralia a d’ailleurs été identifié.

Le dernier assemblage révélé par l'étude pollinique (daté de 600-700 cal BP) reflète un environnement très différent. Les arbres décidus y sont beaucoup plus rares. Les taux plus importants d'épicéa (Picea), de pin (Pinus) et de sapin (Abies) s'expliquent par des apports lointains, probablement issus des Subcarpates de la Courbure. Comme au cours de l'âge du Bronze, nous enregistrons un léger signal d'activité agro-pastorale avec quelques pollens des plantes messicoles et plusieurs spores de champignons co-prophiles.

Pour cette période, les données palynologiques et les microfossiles non polliniques indiquent un nouveau changement des paramètres édaphiques et hydrologiques avec un accroissement des plantes aquatiques et une raréfaction des Chenopodiaceae.

3. Conclusions

Les données palynologiques obtenues dans le cadre
de cette étude ne sont pas très abondantes et couvrent
une période restreinte de l'Holocène, mais dans un
secteur où aucune donnée n'existait jusqu'à présent,
elles ont le mérite de livrer une première image
des formations végétales présentes durant l'âge du
Bronze, puis au cours des derniers siècles.

Dans cet enregistrement pollinique, nous n'obser-
vons aucun témoignage direct de céréaliculture (à
l'exception d'un unique grain de pollen de céréale),
bien que, pour la période Chalcolithique, l'étude des
phytolithes ait révélé une intense pratique agricole à
Taraschina, au cœur du delta (Danu et al., 2019). La
quasi-absence de signaux d'anthropisation dans l'en-
registrement pollinique de la séquence 14 pourrait
s'expliquer par le fait que l'enregistrement pollinique
débute juste à la fin de l'occupation chalcolithique et
renseigne un environnement contemporain de l'âge
du Bronze. C'est à partir de cette période que le sec-
teur du delta du Danube semble se remettre en eau
comme en témoignent les données sédimentologiques.
Il est donc très probable que, à la fin de la période
chalcolithique, la périphérie du site de Taraschina soit
devenue un lieu peu propice à la céréaliculture, en
raison de la progression des espaces d'eau, comme
en témoigne l'augmentation des taxons hygrophiles.
Ceci expliquerait la quasi-absence des signes d'agri-
culture dans l'enregistrement pollinique.

Les données polliniques obtenues permettent de
décrire les formations forestières qui s'étendaient
dans le secteur à l'âge du Bronze : ripisylve (aulne,
saule) et forêt caducifoliée (chêne, charme, tilleul).

Une hausse des valeurs de pollen de pin, de sapin et d'épicéa est enregistrée pour la fin du Moyen-Âge. Les arbres semblent ainsi mieux représentés dans la partie finale de la séquence, mais il est actuellement difficile de reconstituer leur distribution passée dans le delta du Danube, car l'abondance de ces conifères est le reflet d'apports polliniques lointains.
## Rermerciements

Nous adressons nos plus sincères remerciements au CEPAM NICE (UMR 7264) et à GEODE (UMR 5602) qui nous ont permis de réaliser ce travail au sein de leur laboratoire. Mihaela Danu a bénéficié d'une bourse de postdoctorat financée par le Ministère des Affaires Étrangères (France).

## Références bibliographiques

**Beug 2004 :** BEUG (H.J.) – *Leitfaden der Pollenbestimmung für Mitteleuropa und angrenzende Gebiete*, München: Pfeil.

**Bohn, Gollub, Hettwer 2000 :** BOHN (U.), GOLLUB (G.), HETTWER (C.) – *Karte der natürlichen Vegetation Europas, Bundesamt für Naturschutz Federal Agency for Nature Conservation*, (Bohn-Bad Godesberg).

**Bozilova, Beug 1992 :** BOZILOVA (E.), BEUG (H.J.) – On the Holocene history of vegetation in SE Bulgaria (Lake Arkutino, Ropotamo region), *Veget. Hist. Archaeobot*. 1, p. 19-32.

**Bozilova, Beug 1994 :** BOZILOVA (E.) BEUG (H.J.) – Studies on the vegetation history of Lake Varna region, northern Black Sea coastal area of Bulgaria, *Veget. Hist. Archaeobot*. 3, p. 143-154.

**Bozilova, Filipova 1986 :** BOZILOVA (E.), FILIPOVA (M.) – Palaeoecological environment in NE Black Sea area during Neolithic, Eneolithic and Bronze Periods, *Studia Praehistorica* 8, p. 160-165.

**Bozilova, Ivanov 1985 :** BOZILOVA (E.), IVANOV (I.) – Palaeoenvironment in the area of the Varna Lake during Eneolithic and Bronze Age on the basis of palynological, paleoethnobotanical and archaeological evidence, *Bull Museu Nat. Varna* 21, p. 43-78.

**Bozilova, Tonkov 1985 :** BOZILOVA (E.), TONKOV (S.) – Palaeoecological studies in Lake Durankulak. *Ann Univ Sofia Fac Biol (Bot)* 76, p. 25-30.

**Bozilova, Tonkov 1998 :** BOZILOVA (E.), TONKOV (S.) – Towards the vegetation and settlement history of

the southern Dobrudza coastal region, north-eastern Bulgaria: a pollen diagram from Lake Durankulak, *Veget. Hist. Archaeobot*. 7, p. 141-148.

**Bozilova, Tonkov 2002 :** BOZILOVA (E.), TONKOV (S.) – Paleoecological evidence on the vegetation history and human occupation in the coastal area of Lake Durankulak, Northeastern Bulgaria, in H. Todorova (ed.), *Durankulak. Band II. Die Prahistorischen Graberfelder Von Durankulak, Teil 1*, Sofia, p. 309-312.

**Cugny 2011 :** CUGNY (C.) – *Apports des microfossiles non-polliniques à l'histoire du pastoralisme sur le versant nord Pyrénéen : entre référentiels actuels et reconstitution du passé*, Géographie. Université Toulouse le Mirail - Toulouse II. Français. https://www.theses.fr/2011TOU20089

**Danu et al. 2019 :** DANU (M.), MESSAGER (E.), CAROZZA (J.-M.), CAROZZA (L.), BOUBY (L.), PHILIBERT (S.), ANDERSON (P.), BURENS (A.), MICU (C.) – Phytolith evidence of cereal processing in the Danube Delta during the Chalcolithic period, *Quaternary International* 504, p. 128-138.

**Doveri 2004 :** DOVERI (F.) – *Fungi Fimicoli Italici. A Guide to the recognition of Basidiomycetes and Ascomycetes living on faecal material*. (Ed), Associazione Micologia Bresadola, Trento, Italy, 1104 p.

**Fægri, Iversen 1989 :** FÆGRI (K.), IVERSEN (J.) – *Textbook of pollen analysis, Chichester*: John Wiley and Sons.

**Feurdean et al. 2020 :** FEURDEAN (A.), GRINDEAN (R.), FLORESCU (G.), TANȚĂU (I.), NIEDERMEYER (E.), DIACONU (A.-C.), HUTCHINSON (S.), NIELSEN (A.), SAVA (T.), PANAIT (A.), BRAUN (M.), HICKLER (T.) – The transformation of the forest steppe in the lower Danube Plain of south-eastern Europe: 6000 years of vegetation and land use dynamic, *Biogeosciences* 18 (3), p. 1081-1103.

**Filipova 1985 :** FILIPOVA (M.) – Palaeoecological investigations of Lake Shabla- Ezeretz in NE Bulgaria, *Ecol Mediterr*. 11, p. 147-158.

**Filipova-Marinova, Bozilova 1990 :** FILIPOVA-MARINOVA (M.), BOZILOVA (E.) – Palaeoecological investigation of Lake Shabla-Ezeretz on the Bulgarian Black Sea coast, in T. Krastev (ed) *Geological evolution of the western part of the Black Sea in Neogen-Quaternary*, Publishing House Bulg Acad Sci, Sofia, p. 41-87.

**Filipova-Marinova et al. 2013 :** FILIPOVA-MARINOVA (M.), PAVLOV (D.), COOLEN (M.), GIOSAN (L.) – First high-resolution marinopalynological stratigraphy of Late Quaternary sediments from the central part of the Bulgarian Black Sea area, *Quaternary International*, 293, p. 170-183. doi:10.1016/j.quaint.2012.05.002

**Filipova-Marinova, Pavlov, Giosan 2016 :** FILIPOVA-MARINOVA (M.), PAVLOV (D.), GIOSAN (L.) – Multi-proxy records of Holocene palaeoenvironmental changes in the Varna Lake area, western Black Sea coast, *Quaternary International*, 401, p. 99-108. doi:10.1016/j.quaint.2015.05.009
**Marinova 2003 :** MARINOVA (E.) – The new pollen core Lake Durankulak-3: a contribution to the vegetation history and human impact in Northeastern Bulgaria, in S. Tonkov (ed) *Aspects of palynology and palaeoecology*. Pensoft, Sofia-Moscow, p. 257-268.

**Marinova, Atanassova, 2006 :**
MARINOVA (E.), ATANASSOVA (J.) – Anthropogenic impact on vegetation and environment during the Bronze Age in the area of Lake Durankulak, NE Bulgaria: Pollen, microscopic charcoal, non-pollen palynomorphs and plant macrofossils, *Review of Palaeobotany and Palynology* 141 (1-2), p. 165-178.

**Mudie et al. 2015 :** MUDIE (P.), MUDRYK (I.), GERASIMENKO (N.), ROHOZIN (E.) – *Palynology of the Danube Delta and Shelf - Comparison of Methods*. 2015 GSA Annual Meeting, At

Baltimore, Maryland, USA, 47 (7), 204 p.

**Popescu et al. 2020 :** POPESCU (S.-M.), JIMÉNEZ-MORENO (G.), KLOTZ (S.), LERICOLAIS (G.), GUICHARD (F.), NAMIK ÇAĞATAY (M.), GIOSAN (L.), CALLEJA (M.), FAUQUETTE (S.), SUC (J.-P.) – *Late Quaternary vegetation and climate of SE Europe – NW Asia according to pollen records in three offshore cores from the Black and Marmara seas, Palaeobiodiversity and Palaeoenvironments*, 1-25, Doi: 10.1007/s12549-020-00464-x

**Reille 1992 :** REILLE (M.) — *Pollen et spores d'Europe et d'Afrique du Nord. Marseille, Laboratoire de Botanique Historique et Palynologie*, 543 p.

**Tonkov et al. 2011 :** TONKOV (S.), BEUG (H.-J.), BOZILOVA (E.), FILIPOVA-MARINOVA (M.),

JUNGNER (H.) – Palaeoecological studies at the Kaliakra area, northeastern Bulgarian Black Sea coast: 6000 years of natural and anthropogenic change, *Vegetation History and Archaeobotany*, 20(1), p. 29-40. doi:10.1007/s00334-010-0244-4

**Tonkov et al. 2014 :** TONKOV (S.), MARINOVA (E.), FILIPOVA-MARINOVA (M.), BOZILOVA (E.) – Holocene palaeoecology and human-environmental interactions at the coastal Black Sea Lake Durankulak, northeastern Bulgaria, *Quaternary International*, 328-329, p. 277-286. doi:10.1016/j.quaint.2013.12.004

**van Geel 2001 :** VAN GEEL (B.) – Non-pollen palynomorphs, in JP Smol, H.J.B. Birks, W.M. Last (eds) *Tracking environmental change using lake sediments (terrestrial, algal and silicaceous indicators)*, vol 3. Kluwer, Dordrecht, p. 99-119.
Chapitre 2

*Le tell submergé de Taraschina :
approches géophysiques et géo-archéologiques*
<figure><img src="image_71.png" /></figure>
# Le tell de Taraschina topographie et paysages

Laurent CAROZZA et Cristian MICU

Le tell de Taraschina se trouve au cœur de l'actuel delta du Danube, dans une zone de lacs naturels et de roselières, zébrée de canaux artificiels, témoins des aménagements réalisés durant la seconde moitié du 20e siècle, pendant la période communiste. La découverte du site archéologique remonte au début des années 1980, à l'appui de ramassages de surface. Parmi les cartes anciennes, seule celle établie par Grigore Antipa au début du 20e siècle figure l'emplacement exact du site chalcolithique (Antipa, 1914). Cet érudit n'a pas identifié le site, mais la cartographie qu'il a dressée

des lieux montre une forme ovalaire se détachant nettement dans un paysage de marais, de lacs et de bras morts.

L'analyse de différentes cartes montre que le site archéologique forme une butte artificielle qui émerge des eaux du delta à une altitude d'environ 2 m ALS. Le tell de Taraschina apparaît pour la première fois sur la carte du delta du Danube dressée par Grigore Antipa (fig. 1), en annexe de son ouvrage « Câteva probleme ştiinţifice şi economice privitoare la Delta

<figure><img src="image_72.png" /><figcaption>**Figure 1**
Extrait de la carte établie par Grigore Antipa au début du 20e siècle.</figcaption></figure>
Dunărei » édité en 1914 (Antipa, 1914). Cette carte, d'une précision remarquable, montre l'existence, à la sortie du lac Fortuna, sur la rive gauche du bras Șontea, d'une butte témoin isolée de forme ovalaire (fig. 1). Les cartographes ont matérialisé ce relief isolé, pourtant de taille modeste, qu'ils ont interprété comme une ride ou une levée de berge et non pas comme un édifice artificiel ou un site archéologique (sur la figure 1, les couleurs orangées matérialisent toutes les formes de relief, tel que des levées de berge ou des rides sableuses). Les cartes récentes et les images aériennes que nous avons compilées montrent que cette zone a été largement artificialisée, avec d'une part la fixation du bras de Șontea devenu canal navigable, et d'autre part avec le creusement de

multiples canaux destinés à exploiter les ressources de cette zone (pêche, roseaux, poldérisation...). Sur la carte dressée par Antipa, la taille du tell semble exagérée. Par ailleurs le cartographe a situé Taraschina à proximité du bras Șontea. Si différents documents montrent que le tracé du canal Șontea a été modifié, le positionnement du site de Taraschina sur la carte de 1914 n'est pas suffisamment précis pour en déduire des informations paléo-géographiques.

Aujourd'hui, le site archéologique de Taraschina forme une petite surface émergée d'environ 1 m de la surface des roseaux du delta du Danube. Durant les périodes de crues, cette zone est totalement submergée (fig. 2) et le petit îlot de Taraschina disparaît sous les eaux du Danube. Le site archéologique se trouve à l'est du lac Fortuna (fig. 1 et 3), dans un environnement dont la structure actuelle est fortement artificialisée. Le creusement d'un petit chenal, aménagé dans les années 1970, qui lie le bras Gârla Șontea et le canal Olguța, a recoupé la frange orientale de l'habitat chalcolithique. Les coupes laissent apparaître du mobilier archéologique et une partie des déblais retjetés à l'est du chenal masque les contours du site archéologique.

Mais le caractère exceptionnel et unique du site de Taraschina tient avant tout à sa position très avancée dans le delta actuel du Danube (fig. 4). Jusqu'alors, les témoignages les plus anciens découverts en contexte alluvial dans le delta étaient attribués à l'âge du Bronze. On savait toutefois que Puiu Hașotti

<figure><img src="image_73.png" /><figcaption>**Figure 2**

Vue du site de Taraschina en cours de fouille, recouvert par les eaux du Danube.</figcaption></figure>

<figure><img src="image_74.png" /><figcaption>**Figure 3**

Localisation du tell chalcolithique de Taraschina dans la zone centrale du delta du Danube. Le site se trouve entre les canaux Olguța et Șontea, au nord du bras de Sulina.</figcaption></figure>
<figure><img src="image_75.png" /><figcaption>Figure 4</figcaption></figure>

Localisation du tell chalcolithique de Taraschina et des principaux sites étudiés dans le cadre de la mission «Archéologie du delta du Danube» (Roumanie).

1 : Taraschina ; 2 : Luncaviţa ; 3 : Baia tell «Boruz» ; 4 : Ceamurlia de Jos ; 5 : Baia Goloviţa ; 6 : Lunca ; 7 : Borduşani ; 8 : Hârşova ; 9 : Vităneşti ; 10 : Pietrele Măgura Gorgana ; 11 : Sultana Malu Roşu.
avait mentionné la découverte de quelques éclats en silex sur la ride Caraorman ; l'un d'eux lui paraissant présenter des traces d'utilisation (Haşotti, 1997, p. 13, note 8). Toutefois, ces données n'ont pas pu être vérifiées en étudiant ce mobilier.

La mise en évidence, sur le site de Taraschina, de vestiges datant de l'époque chalcolithique sur une superficie d'environ 1 ha (fig. 5), témoigne de l'ampleur de l'établissement. La densité du mobilier céramique ainsi que la présence d'industries lithique, osseuse

<figure><img src="image_76.png" /><figcaption>**Figure 5**
Localisation des zones de fouille sur le tell de Taraschina à Maliuc.</figcaption></figure>
et de pièces de broyage accentuent l'impression du caractère pérenne de l'établissement (Carozza et al., 2011). Cette découverte bouleverse la perception que l'on pouvait avoir de la paléo-géographie de l'em-bouchure du Danube durant les périodes Néolithique et Chalcolithique. L'intégration de données de peuplement ancien dans le delta du Danube permet de préciser la chronologie de l'évolution du delta, jusqu'ici principalement appuyée sur des informations radio-chronologiques et sur des reconstitutions paléogéographiques.

# 1. Topographie de la surface du tell de Taraschina

La micro-topographie du tell de Taraschina a été réalisée en 2015, à la faveur d'un étiage particulièrement marqué en fin d'été, période durant laquelle les basses eaux du delta nous ont permis d'accéder aux zones périphériques du tell. Cette tâche a pu être réalisée en couplant des données Lidar qui ont été mises à disposition par l'INDD (Institut National de Recherche sur le delta du Danube) dans le cadre d'un partenariat scientifique, à des mesures de terrain réalisées à l'aide d'une station totale. Ce travail de terrain a permis de « débruiter » les données LIDAR, car la présence d'une végétation luxuriante au moment de l'acquisition aéroportée ne permettait pas de disposer d'un modèle suffisamment homogène pour restituer la microtopographie du site. L'ensemble des données acquises dans la zone de Taraschina a été rattaché au système altimétrique absolu (ASL).

La micro-topographie complète du tell de Taraschina permet de distinguer deux sites archéologiques accolés (fig. 5 et 6). Les carottages et les prospections effectués sur le site ont montré que nous étions en présence de deux habitats chalcolithiques, pour partie contemporains. Le site 2, bien que moins élevé, est chronologiquement postérieur à l'occupation du tell principal (site 1 ; fig. 5).

Le tell principal, dit site 1, présente une forme quadrangulaire marquée sur les flancs sud et nord du site. L'analyse des ortho-photographies et des cartes anciennes montre également que le tracé du canal qui lie le lac Taraschina au canal Şontea n'a pas été établi de manière aléatoire, mais de façon à accéder en bateau sur le site. La surface du tell chalcolithique a, en effet, été utilisée durant la période communiste (dans les années 1970) comme plateforme de stockage des roseaux récoltés dans la zone de Taraschina. Le coude formé par le canal tire parti de l'ancienne topographie du tell, bordé sur sa marge orientale par le canal. Le site 1 se développe actuellement sur une

longueur de 90 m pour une largeur de 85 m. Si l'on prend en compte la partie du tell détruite par le creusement du canal, on peut estimer la longueur originelle du tell à 100 m. On peut par conséquent évaluer la surface du tell principal à 0,85 ha.

Le site 2 présente des dimensions nettement plus réduites et un plan grossièrement quadrangulaire de 35 m de long pour une largeur de 30 m. Il forme une excroissance adossée au nord du tell principal.

Au total, nous pouvons estimer à près de 1 ha la surface de la zones habitée et aménagée durant le 5ᵉ millénaire avant notre ère à Taraschina.

Par sa superficie, le tell chalcolithique de Taraschina constitue sans nul doute l'un des sites les plus importants de Dobroudja. À titre de comparaison, le tell de Lunca Pensiune, situé sur les rives du lac Razim, présente une superficie d'environ 3 000 m² (long. : 53 m ; larg. 65 m). Sa puissance stratigraphique est également plus réduite qu'à Taraschina, avec une épaisseur de 2 m au centre de l'édifice. Les données

<figure><img src="image_77.png" /><figcaption>**Figure 6**
Relevé topographique d'un sol chalcolithique du tell de Taraschina.</figcaption></figure>
disponibles pour le tell de Baia indiquent que ce site présentait un plan ovalaire (65 m de long sur 75 m de large) et une superficie estimée à 4 000 m², pour une puissance stratigraphique maximale estimée à 2,80 m. Quant au tell de Luncaviţa, il présente un plan légèrement ovalaire de 100 m de long sur 75 m de large, soit une superficie d'environ 7600 m², pour une puissance stratigraphique de 4 m dans sa partie la plus développée.

Taraschina, malgré une stratigraphique moindre qu'à Luncaviţa, est doté d'une surface analogue à celle de certains tells de la plaine roumaine, comme celui de Borduşani Popină (d'un diamètre d'environ 100 m) ou celui de Sultana Malu Roşu. Les tells de Dobroudja du nord présentent des stratigraphies bien

moins développées, comme à Sultana-Magura Malu Roşu (diamètre de 100 m pour une hauteur de 4 m), à Pietrele Măgura Gorgana (95 m de diamètre pour une hauteur de 6,30 m) ou à Vităneşti Măgurice (95 m de diamètre pour 6 m de haut).

Situé dans l'actuel delta du Danube, le tell de Taraschina n'apparaît pas intégralement dans le paysage, ce qui en fait une spécificité au regard des tells contemporains. Le site est noyé dans les sédiments du Danube et seul le tiers supérieur du site émerge des eaux du delta. La base du tell et la morphologie du paysage sont totalement soustraits à la vue. Cette particularité, qui nous empêche également de fouiller les horizons inférieurs du site de Taraschina, nous a contraint à mettre en œuvre des approches innovantes.

Références bibliographiques

**Antipa 1914** : ANTIPA (G.) – *Delta Dunarii*. Bucarest, 1914.

**Carozza, Bem, Micu 2011** :
CAROZZA (L.), BEM (C.), MICU (C.) – Société et environnement dans la zone du bas Danube durant le 5ᵉ millénaire avant notre ère. Iasi, Editura Universităţii Alexandru Ioan

Cuza, 414 p.

**Carozza, Micu, Haită, Ailincăi, Burens, Mihail, Carozza 2014** :
CAROZZA (L.), MICU (C.), HAITĂ (C.), AILINCĂI (S.), BURENS (A.), MIHAIL (F.), CAROZZA (J.-M), 2014 – Reconnaissance par carottages de la stratigraphie des

habitats pluri-stratifiés de Taraschina dans le delta du Danube. *Dacia*, N.S. LVIII, p. 13-27.

**Haşotti 1997** : HAŞOTTI (P.) – *Epoca neolitica în Dobrogea, Bibliotheca Tomitana I, Constanţa*.
# Le tell de Taraschina
## premières recherches, premiers résultats et développements

Laurent CAROZZA, Cristian MICU et Albane BURENS

L'un des principaux objectifs de cette monographie est de présenter le résultat des recherches menées dans le delta du Danube depuis 2011, date de la publication des premiers travaux conduits sur le tell submergé de Taraschina (Carozza, Bem, Micu, 2011). Nous avions alors précisé l'attribution chrono-culturelle des occupations des niveaux supérieurs du tell, défini les contextes stratigraphiques et présenté les premières hypothèses paléo-économiques.

Depuis, la poursuite des recherches de terrain, conduites à l'appui d'une grande diversité d'approches telles que les prospections géophysiques de surface, les carottages géo-archéologiques *intra* et hors site, les reconstructions paléo-environnementales (...) nous permet aujourd'hui d'apporter de nouvelles données susceptibles de faire évoluer les hypothèses initialement formulées.

À Taraschina, les travaux de terrain ont débuté en 2009 par la conduite de prospections de surface (fig. 1). Parallèlement, des sondages avaient été réalisés dans des zones à forte concentration de mobilier. En 2010, date de la création de la mission archéologique et de son programme pluriannuel de recherche, une première zone de fouille (zone 1) de 72 m² a été ouverte dans les carrés W-R/50-52 (fig. 1). Cette petite zone de fouille a permis de reconnaître la séquence stratigraphique sur une puissance de 1,5 m, la nappe phréatique n'ayant pas permis d'aller plus bas. Les données acquises ont contribué à caractériser les modalités d'abandon de l'habitat chalcolithique et d'en appréhender l'organisation spatiale.

En septembre 2011, une seconde zone de fouille (zone 2) a été implantée dans un secteur où les prospections avaient identifié une forte anomalie magnétique, traduisant la présence de terre à bâtir rubéfiée et pouvant s'apparenter à l'emplacement d'une maison incendiée. Au printemps 2013, en raison des fortes

crues du Danube, la poursuite de la fouille de la zone 2 n'a pas été possible en raison de la submersion totale de la surface du site (fig. 2). Les campagnes de terrain ont repris en octobre 2013.

<figure><img src="image_78.png" /><figcaption>Figure 1</figcaption></figure>

Le tell chalcolithique de Taraschina : microtopographie du site et localisation des différentes zones de fouille, des sondages et des zones prospectées.
En 2015, une fois la fouille de la zone 2 achevée, la réalisation d'une campagne de prospection géophysique extensive a permis de dresser une carte des anomalies magnétiques sur l'ensemble de la surface du site et ainsi de déterminer l'emplacement d'une nouvelle zone de fouille (zone 3). Les conditions environnementales de la prospection réalisée en 2015 ont été optimales et l'équipe a pu réaliser un débroussaillage intégral de la surface du tell. La méthode mise en œuvre en 2015 se caractérise par un couplage entre capteurs et station totale en mode « tracker » (fig. 3). Ce dispositif, calé sur Modèle numérique de Terrain (MNT) du tell de Taraschina, a permis de disposer d'une cartographie à très haute résolution spatiale.

<figure><img src="image_79.png" /><figcaption>**Figure 2**
La zone 2 du site de Taraschina entièrement submergée.</figcaption></figure>

<figure><img src="image_80.png" /><figcaption>**Figure 3**
Prospections géophysiques sur le site de Taraschina (LIENSs).</figcaption></figure>

Parallèlement à ces prospections, nous avons réalisé des mesures géophysiques *in situ*, par diagraphie, à l'intérieur des trous de carottage. On doit à François Lévêque et l'équipe Lienss de l'Université de La Rochelle d'avoir développé les matériels et la mise en œuvre de ce dispositif de mesure innovant. Le couplage entre cartographie de surface des anomalies magnétiques et évaluation de leur nature – nature des unités stratigraphiques, composition, texture... – a ainsi permis de caractériser les anomalies et de les interpréter. La systématisation des carottages tubés, prélevés sous forme de deux transects, a livré des informations essentielles tant sur la morphologie du milieu dans lequel la communauté chalcolithique de Taraschina s'était installée, que sur la dynamique de formation du site. La description systématisée des unités sédimentaires et leur interprétation fonctionnelle nous donne une image très précise de la dynamique de l'habitat et des étapes de l'édification du tell. Ces carottes fournissent également des informations paléo-écologiques permettant d'identifier des activités contribuant à mieux appréhender l'économie des communautés agro-pastorales de Taraschina. Elles constituent également une archive conservée pour partie en chambre froide, qui permettra dans les années à venir de poursuivre les investigations au fur et à mesure de l'évolution de nos questionnements.

Sur la base des données géophysiques, nous avons pu identifier, sur la surface du tell de Taraschina, plusieurs bâtiments incendiés, d'apparence non perturbés par les ré-occupations du site. La problématique liée aux habitations incendiées sur les tells requiert de disposer de structures d'habitat très bien préservées, permettant l'étude des contextes au sein desquels sont intervenus ces incendies (feux volontaires *versus* feux accidentels). Ainsi avons-nous ouvert une vaste zone de fouille (zone 3) de près de 180 m², dans le but de circonscrire l'extension de ces constructions, mais également de les mettre en relation avec des espaces connexes, telles que des zones de circulation, des zones fonctionnelles (battage, stabulation ...) et de définir la présence d'éventuels autres bâtiments non incendiés.

Les données acquises en 2015 (cartographie géophysique et carottages) ont montré la faisabilité de la fouille de cet horizon et son importance pour la connaissance des modalités de l'abandon du tell. D'un point de vue stratigraphique, ces bâtiments incendiés se rencontrent au sommet de la séquence stratigraphique et appartiennent probablement à la dernière étape d'occupation du tell (postérieure à 4350 BC). La fouille en cours de cet ensemble dépasse l'enjeu d'appréhender l'architecture et le fonctionnement d'unités domestiques de la période Gumelnița.
Il s'agit avant tout de comprendre dans quel contexte socio-environnemental ces bâtiments ont pris feux et de mettre en lien l'analyse fonctionnelle de ces structures avec le modèle de transformation du tell de Taraschina, tel que nous l'avons établi à l'aide des données déjà acquises. Il s'agit notamment d'observer s'il existe un lien chronologique entre les phases de déplacement de l'habitat, les possibles hiatus de l'occupation, les changements socio-économiques et la récurrence d'événements comme l'incendie de

bâtiments. Cette approche requiert de disposer, à l'échelle des bâtiments incendiés mis au jour dans la zone 3, des données archéologiques suffisamment précises pour appréhender le phénomène de la destruction de bâtiments par incendie. C'est pourquoi la fouille revêt un caractère essentiel pour permettre de valider (ou non) le modèle d'abandon/réoccupation du tell à la fin du 6e millénaire avant notre ère, sous l'influence de changements environnementaux, tel que nous l'avons proposé.

Références bibliographiques

**Carozza, Bem, Micu 2011 : CA-ROZZA (L.), BEM (C.), MICU (C.)**
– Société et environnement dans la zone du bas Danube durant le 5e millénaire avant notre ère. , Iași, Editura Universității Alexandru Ioan Cuza, 414 p. (ISBN : 978-973-703-725-5).

**Carozza et al. 2012 : CAROZZA**
(J.-M.), MICU (C.), MIHAIL (F.),
CAROZZA (L.) – Landscape change

and archaeological settlements in
the lower Danube valley and delta
from Early Neolithic to Chalcolithic
time. *Quaternary International*,
IGCP 521: Caspian-Black Sea –
Mediterranean Corridors during
the last 30 ka: sea level change and
human adaptive strategies Selected
Papers, Volume 261, 30 May 2012,
p. 1–31, http://dx.doi.org/10.1016/j.
quaint.2010.07.017

**Carozza et al. 2014 : CAROZZA**
(L.), MICU (C.), HAITĂ (C.),
AILINCĂI (S.), BURENS (A.),
MIHAIL (F.), CAROZZA (J.-M) –
Reconnaissance par carottages de la
stratigraphie des habitats pluri-
stratifiés de Taraschina dans le delta
du Danube. *Dacia*, 2014, p. 13-27.
<figure><img src="image_81.png" /></figure>
# Voir le tell submergé de Taraschina sous un autre angle :
prospections géophysiques de surface et diographie

François LÉVÊQUE

## 1. Introduction

L'imagerie géophysique de la proche surface permet d'observer des objets présents sous la surface du sol, invisibles à l'œil, en un temps de prospection modeste. La qualité de la détermination de la géométrie des structures identifiées dépend d'une part du contraste de la propriété physique mesurée entre l'objet et l'encaissant et d'autre part de la densité spatiale des mesures réalisées.

L'objet d'étude, le tell de Taraschina, étant en zone humide, la proximité du toit de la nappe d'eau limite fortement la profondeur d'investigation des méthodes fondées sur la propagation d'un signal électrique ou électromagnétique. Les courants électriques empruntant alors le chemin le plus court, ils restent en surface de la nappe d'eau. En revanche, les méthodes magnétiques sont particulièrement bien adaptées à la présence de terre chauffée, voir cuite et sont par conséquent particulièrement bien appropriées à la problématique de ce site, spécifiquement dans ce contexte de zone humide. Les murs des maisons incendiées constituent des sources magnétiques de forte intensité, d'autant que la terre crue, utilisée en tant que matériaux de construction, est formée de sédiments prélevés à proximité.

Le fer, un des éléments majeurs, gouverne les propriétés magnétiques de la matière naturelle. Du fait des conditions physico-chimiques du milieu, le fer est présent dans les sédiments des zones humides de type deltaïque sous des formes minéralogiques peu magnétiques. L'omniprésence de l'eau et l'absence de l'oxygène qui en résulte sont les moteurs des transformations minéralogiques impliquant le fer. En effet, la pososité du sédiment étant saturée par l'eau, les échanges avec l'air, en particulier avec l'oxygène, sont limités.

L'oxygène libre présent dans le sédiment est rapidement consommé par des réactions d'oxydation, par exemple dans les processus de minéralisation de la matière organique (MO). Le milieu devient réducteur. Le fer est alors réduit, passant de fer ferrique (Fe³⁺) à fer ferreux (Fe²⁺). Ce dernier a une solubilité élevée, favorisant son lessivage. Le changement de couleur du sédiment, passant d'un ton ocre-brun en surface à gris-bleu en profondeur, révèle cette évolution.

Si la teneur en MO est élevée, alors en l'absence d'oxygène disponible, le fer en solution s'associe au souffre libéré par la minéralisation de la MO. Ainsi se forment des sulfures de plus en plus riches en souffre jusqu'à la phase ultime, la pyrite (FeS₂), qui n'est pas magnétique (paramagnétique). Une des phases transitoires précurseurs de la pyrite est la greigite (Fe₃S₄). Contrairement à la pyrite, la greigite est une phase magnétique (ferrimagnétique), apte à porter une aimantation en dehors de celle induite par le champ magnétique local, mais elle est métastable à l'air. Elle ne peut donc être présente que dans les sédiments organiques et disparaît si ces sédiments sont mis à l'air. Sa contribution doit donc être négligeable, du moins en surface. Dans la partie de battement de nappe, l'apport modéré d'oxygène favorise la formation de goethite (FeOOH) un oxyhydroxide de fer ferrique, qui est un pigment de couleur jaune. Ce pigment (antiferromagnétique canté, c'est-à-dire avec les atomes de fer répartis dans deux sous réseaux cristallins avec des moments magnétiques imparfaitement anti-parallèles) possède une faible capacité à porter une aimantation permanente. En conséquence, les sédiments présents dans les zones humides du delta, non impactées par des activités humaines, devraient être peu magnétiques et être caractérisés par la présence de pyrites et de goethites à l'état de trace.
Les parois des murs des maisons sont réalisées avec les sédiments prélevés à proximité du site ; ils contiennent donc potentiellement des traces de goethite et de pyrite, associées à des débris végétaux. L'évolution de ce matériau composite, à l'air libre, favorise la formation de magnétite, tout comme cela se produit dans l'horizon A d'un sol qui devient plus magnétique que le substrat sur lequel il se développe. La magnétite (Fe₃O₄, ferrimagnétique) est l'oxyde de fer le plus magnétique. Néanmoins, sa teneur dans la terre à bâtir des murs reste modeste à ce stade. En comparaison avec les murs des maisons actuelles du même secteur, la susceptibilité magnétique est au mieux doublée. En revanche, si le mur est chauffé à des températures dépassant 300 à 350°C alors, du fait de la teneur en eau du milieu et de la présence de matière organique, la vapeur d'eau ou la combustion de la MO engendre un milieu réducteur qui favorise la transformation d'une partie de la goethite en magnétite, très magnétique, plutôt qu'en hématite (oxyde de fer ferrique ocre, Fe₂O₃, anti-ferromagnétique cancérate), peu magnétique tout comme la goethite. Ainsi, la rubéfaction générée par le feu (dûe aux pigments d'hématite), est elle associée à une augmentation des capacités d'aimantation (dûe à la fraction de magnétite produite) qui, dans le cas présent, peut aller jusqu'au centuple [susceptibilité magnétique passant de 0.25 - 0.30 10⁻³ SI à 25 - 30 10⁻³ SI], du moins cumulé avec l'effet de la transformation de la pyrite.

En effet, de manière analogue à la déshydratation de la goethite, la pyrite subissant une oxydation modérée à haute température peut former de la pyrrhotite (appellation qui regroupe plusieurs formes de type Fe₁-xS (x < 0.2) avec la forme Fe₇S₈ la plus remarquable du point de vue magnétique). Ce minéral, stable à l'air contrairement à la greigite, possède aussi une capacité d'aimantation élevée. Lorsque la chauffe se fait en atmosphère confinée – le cœur d'un bloc massif par exemple – l'oxydation est limitée. Au lieu d'une rubéfaction, il y a alors noircissement. Ces conditions sont, dans le principe, favorables à la formation de tels sulfures magnétiques.

L'incendie d'une maison constitue donc la source d'un enrichissement en minéraux magnétiques des matériaux ayant atteint une température supérieure à 300-350°C. Après l'incendie, les matériaux des murs cuits, ou simplement chauffés, acquièrent en refroidissant une aimantation, dite thermorémanente, d'intensité proportionnelle à celle du champ magnétique local, corrélativement à la capacité d'aimantation des matériaux. La direction de cette aimantation correspond à celle du champ magnétique ambiant au moment de l'incendie. Cette aimantation permanente

qui a enregistré le champ magnétique passé, associée à celle induite dans les matériaux par le champ magnétique local actuel, constitue les sources d'une anomalie locale du champ magnétique. L'enrichissement en minéraux magnétiques contribue à accroître l'anomalie locale du champ magnétique engendrée par l'impact du feu sur les matériaux. L'aimantation thermorémanente acquise lors du refroidissement est en effet proportionnelle à la teneur en minéraux magnétiques ; elle est d'autant accrue que l'enrichissement est élevé.

En résumé, la teneur élevée en minéraux magnétiques des vestiges des murs ayant chauffés tranche avec la faible teneur des sédiments de l'environnement du tell. De même, l'anomalie locale du champ magnétique engendrée par le feu est elle amplifiée par cet enrichissement en minéraux magnétiques. De ce fait, les méthodes magnétiques, fondées sur les propriétés magnétiques intrinsèques aux matériaux (susceptibilité ou viscosité magnétique) ou à leur conséquence sur le champ géomagnétique local (champ produit par les aimantations rémanente et induite), paraissent constituer les outils de la géophysique les plus adaptés pour déterminer la géométrie des maisons incendiées sur le tell.

## 2. L'information magnétique

L'idéal, pour retrouver la géométrie des vestiges des maisons brûlées, serait de déterminer la variation de la concentration en minéraux magnétiques, ou des aimantations naturelles permanentes qu'ils portent, dans un volume sous la surface du sol. En effet, comme nous venons de le décrire, les valeurs élevées d'aimantation ou de concentration en minéraux magnétiques marquent l'emplacement des matériaux ayant chauffés.

La concentration peut être estimée par la capacité d'aimantation d'un volume élémentaire de sol, par exemple après qu'il a été soumis à un champ magnétique, apte à le réaimanter, jusqu'à atteindre la saturation de son aimantation (aimantation rémanente isotherme). De même, la mesure de l'aimantation naturelle se réalise sur un échantillon de volume élémentaire et d'orientation connue. Pour procéder à ces mesures, il faut donc prélever un très grand nombre d'échantillons, puis retourner au laboratoire pour les réaliser. La démarche apparaît caduque, tant par son caractère fortement invasif que par son approche fastidieuse. L'information idéale n'est donc pas accessible dans une démarche non invasive et dans un temps d'acquisition réduit. Pour atteindre ces objectifs, il faut accepter une information approchante.
Sur le terrain, la mesure de la concentration en minéraux magnétiques peut être approchée par la mesure de la susceptibilité magnétique. Cette grandeur, qui correspond au rapport entre une aimantation et le champ magnétique qui la produit, est inféodée à la nature des constituants non magnétiques (inaptes à porter une aimantation permanente, mais portant une aimantation induite par le champ magnétique environnant). Dans le cas d'une concentration en minéraux magnétiques négligeable, la susceptibilité magnétique traduira la proportion de quartz et de calcite (diamagnétique) par rapport aux argiles (paramagnétique). Dans le cas contraire, cette grandeur reflètera la teneur en minéraux les plus magnétiques, telle que la magnétite, avec, comme nous le verrons, un rôle non négligeable de la taille des grains.

Les susceptibilimètres de terrain, ou Kappameter, sont des instruments de contact explorant des volumes pouvant être exprimés en cm³. La mesure est très fortement dépendante de la rugosité de la surface mesurée. L'alternative est l'utilisation de conductivimètres, explorant de beaucoup plus grands volumes, pouvant alors être exprimés en m³. Un conductivimètre mesure une susceptibilité magnétique relative, dite apparente. Les valeurs sont difficilement transposables en vraies valeurs de susceptibilité magnétique telles qu'elles sont mesurées sur échantillon ou avec des capteurs de contact (Kapameter). Ces conductivimètres ne nécessitent pas de contact ; de par leur important volume d'investigation, ces instruments permettent de dégrossir l'étude d'un site. Une grande surface peut être couverte en un temps réduit, mais à une résolution spatiale limitée. Les informations rapidement acquises permettent alors de cibler les zones d'intérêt à étudier à l'aide de méthodes de plus haute résolution, mais nécessitant un temps d'investigation plus long. Dans le cas d'un tell, cette étape de changement d'échelle n'est pas nécessaire compte-tenu de la surface limitée du site.

Les profondeurs d'investigation permises par les Kappameters étant très limitées (quelques centimètres au mieux (Lecoanet et al., 1999), ces outils sont utilisables pendant une fouille archéologique, sur une surface décapée. Leur mise en œuvre reste fastidieuse si l'objectif est de relever une information transposable sous forme cartographique, chaque mesure nécessitant plusieurs secondes. De plus, la mesure réalisée est fortement influencée par la distribution de la taille des particules magnétiques nanométriques. En effet, la réponse du matériau en contact avec le capteur n'est explorée que sur un domaine temporel très bref, inférieur au millième de seconde. Selon leur taille, les particules ont une réaction d'amplitude variable.

Leur temps de réponse à l'excitation par le champ auquel elles sont soumises est aussi dépendant de leur taille. Ainsi, certaines particules, dites superparamagnétiques, très petites, ont une réponse bien plus élevée que des particules légèrement plus grosses (on reste dans le domaine nanométrique, les valeurs de tailles critiques dépendent du minéral considéré). Leur temps de réponse (ou plus exactement leur temps de relaxation) est inférieur aux durées du domaine temporel exploré ; elles se comportent comme des particules paramagnétiques, mais avec une réponse bien plus élevée. De ce fait, une valeur de susceptibilité magnétique ne peut être transposée en une concentration en minéraux magnétiques.

Sur le terrain, la concentration en minéraux magnétiques n'est donc pas directement accessible par de simples mesures à la surface du sol encore végétalisé. En revanche, la variation spatiale de cette concentration en minéraux magnétiques contribue à celle du champ magnétique local, tout comme les aimantations permanentes que portent les minéraux mais aussi l'ensemble des propriétés magnétiques des matériaux environnant, qu'ils soient magnétiques ou non. Ainsi, en accédant à cette variation spatiale du champ géomagnétique local, nous obtenons des informations sur la géométrie d'objets présents sous la surface du sol à partir des contrastes de propriétés magnétiques avec leur environnement.

Essayons de préciser la nature de l'information. Comme nous l'avons vu, une aimantation induite par le champ magnétique local est créée au sein des matériaux, proportionnelle à leur susceptibilité magnétique et à l'intensité du champ qui la produit. Cette aimantation induite s'ajoute à l'aimantation rémanente portée par les matériaux qui est, quant à elle, indépendante du champ magnétique ambiant. Ces aimantations rémanente et induite portées par un matériau créent un champ magnétique autour de lui, qui s'estompe en s'en éloignant. Ce champ d'origine locale s'ajoute au champ magnétique régional, créant une anomalie géomagnétique locale.

La détermination de la variation spatiale de l'intensité du champ magnétique terrestre permet donc d'observer ces anomalies. C'est l'objet de la prospection géomagnétique. Cette méthode de prospection géophysique ne nécessite pas de contact. Grâce aux cadences de mesure élevées des magnétomètres actuels, il est aisé d'obtenir une densité de mesures élevée en un temps de prospection modeste.

Dans le cadre de cette étude, l'information mesurée est l'intensité du vecteur champ magnétique. Cette intensité varie aux environs d'une source en fonction
de sa profondeur, sa distance, sa forme, l'intensité des
aimantations permanentes qu'elle porte, le contraste
de susceptibilité magnétique avec l'encaissant, la
direction du nord magnétique par rapport à celle de
la source, l'inclinaison et l'intensité du champ ma-
gnétique régional, la topographie... Parmi ces para-
mètres, les premiers sont les inconnues. Les derniers,
bien que mesurables, sont communément déterminés
à partir de modèles géomagnétiques globaux. Des
modèles mathématiques permettent une reconstitution
de la géométrie des sources, au prix de nombreuses
hypothèses simplificatrices. Malheureusement, dans
le cadre de cette étude, la distance réduite entre les
sources et le plan dans lequel les mesures sont réali-
sées, la non coaxialité des aimantations rémanentes et
induites, la variation topographique non négligeable
au regard de la dimension du volume exploré... sont
des paramètres difficiles à intégrer dans les modèles.
De ce fait, une reconstitution de la géométrie des
sources, en particulier concernant leur profondeur,
est très incertaine.

Toutefois, afin de pouvoir préciser la position des ves-
tiges des maisons incendiées sur le tell, une prospection
géomagnétique originale 3D a été mise en œuvre. Le
protocole développé vise plus particulièrement à pré-
ciser la profondeur des sources, surtout pour les plus
profondes. Cette prospection se décompose en une
prospection 2D horizontale (2DH) à résolution spatiale
de l'ordre de 20 mesures par m², et une prospection
1D verticale (1DV) correspondant à une diagraphie
en puits de carottage, à une résolution infracentimé-
trique. La fusion spatiale des mesures est réalisée
grâce au géopositionnement 3D (x, y, z) de toutes
les mesures (2DH et 1DV). Pour cela, un réflecteur
360° est solidarisé aux différents capteurs de mesure
du champ magnétique utilisés, par l'intermédiaire
d'un dispositif spécifique à chaque prospection. La
poursuite cinématique du réflecteur, concomitante à
la mesure en continue de l'intensité du champ ma-
gnétique, est réalisée à l'aide d'une station totale
motorisée. La fusion entre les mesures d'intensité
de champ magnétique et leur géopositionnement est
opérée sur la base de l'heure d'acquisition. La syn-
chronisation des instruments est établie a posteriori à
partir d'un évènement enregistré simultanément par
les deux instruments.

Dans un premier temps, le tell a été débroussaillé par
les membres de l'équipe afin de disposer d'un es-
pace dégagé. Par la suite, la prospection 2DH a été
conduite selon des profils parallèles distants d'environ
1 m, en maintenant le capteur bas proches du sol, à
une hauteur la plus régulière possible, de l'ordre de
0,2 m. L'implantation des sites de carottages a ensuite

été établie sur la base des anomalies géomagnétiques
observées. Une fois les carottages réalisés, une
diographie a été effectuée dans chaque puits, à l'aide
d'un dispositif original, permettant ainsi d'accéder à
une information verticale, sous la surface.

3. Matériels et méthodes

La prospection 2DH est réalisée avec un magnétomètre
GSMP-35G GEM Systems Inc, pourvu de deux capteurs,
programmé à une cadence de mesure de 20 Hz. Ce ma-
gnétomètre, fondé sur la technologie dite à pompage
optique, fait partie des magnétomètres absolus. Son
erreur d'étalonnage ne dépasse pas 0.1 nT. Les deux
capteurs sont disposés verticalement, l'un au-dessus
de l'autre, à 1 mètre de distance. Ces capteurs sont
placés à l'extrémité d'un châssis de portage, fait mai-
son, réalisé en tube de fibre de verre. La batterie Li-
ion ainsi que la console de mesure sont placées en
contrepoids à l'extrémité opposée, à environ 2,5 m,
afin d'éloigner ces sources de perturbations électro-
magnétiques des capteurs (fig. 1). L'exactitude de la
mesure obtenue dans cette configuration est meilleure
que 0.1 nT, alors qu'avec le montage compact proposé
par le constructeur, l'anomalie générée par la proximité
des capteurs avec la console et la batterie génère une
anomalie d'une intensité maximale de 0.5 nT, selon
la direction par rapport au nord magnétique. Pour le
géopositionnement des mesures par poursuite laser,
un prisme 360° est fixé au-dessus du capteur haut.
La surface prospectée par l'opérateur est couverte
en parcourant l'espace en continu, à une vitesse de
l'ordre de 1 m/s, selon des profils allers- retours
parallèles. La trajectoire est établie visuellement à
l'aide de repères, matérialisant l'extrémité de chaque
profil, espacés tous les 1 m. Dans la mesure du possible,
le capteur inférieur est maintenu par l'opérateur à une
distance de l'ordre de 0,2 m au-dessus de la surface du sol.

Contrairement à la démarche habituelle de la mesure
en configuration gradiomètre verticale, avec les deux
capteurs superposés, la correction de la variation
temporelle du champ magnétique n'est pas établie
par la détermination du pseudogridient selon cette
verticale (différence des signaux capteur bas - capteur
haut). La variation temporelle est estimée par l'ajus-
tement d'une fonction polynomiale d'ordre 3 à 6 –
selon la nécessité – du signal temporel acquis par le
capteur haut. Le premier ordre de variation de ce signal
temporel correspond à la période des allers-retours.
L'intérêt de cette approche est de conserver les va-
riations de basse fréquence spatiale enregistrées par
le capteur haut, perdu dans le calcul du pseudo-
gradient, pour privilégier les sources profondes qui
sont à l'origine de ces basses fréquences spatiales. Ce
gain est au détriment du non retrait des fluctuations temporelles locales du champ magnétique (pulsations continues et intermittentes) de période courte. Dans le cas de cette étude, compte-tenu de la durée pour réaliser deux profils successifs en sens aller puis retour, les fluctuations de période inférieure à 2 ou 3 minutes ne sont pas retirées. Au pire, en période de forte activité géomagnétique, ces fluctuations ne dépassent pas une amplitude supérieure à quelques nT. En période de faible activité, ces fluctuations ont une amplitude inférieure à 1 nT. La position du site, éloignée de toute zone urbaine, permet de se soustraire des sources de perturbations anthropiques généralement observées près des réseaux électriques.

Concernant la diographie géomagnétique (1DV), le diamètre du carottage ainsi que la longueur de câble nécessaire entre le capteur et la console d'enregistrement limitent fortement la technologie du magnétomètre utilisable. Avec la technologie actuelle, la miniaturisation nécessaire du capteur implique d'utiliser un magnétomètre de technologie fluxgate. Ces instruments sont vectoriels et ne peuvent être considérés comme absolus, présentant entre autres une dérive thermique. Celui mis en œuvre est un fluxgate triaxe Mag 03MSSL70 Bartington connecté à un convertisseur 24 bits DLM24Mag Cx STL Systemtechnik Ludwig GmbH, piloté par un ordinateur portable. Un étalonnage scalaire du capteur triaxe est réalisé à chaque série de diagraphies selon la procédure décrite par Olsen et al., 2001. Bien que la précision de chaque axe soit meilleure que 6 pT à 1 Hz, après étalonnage scalaire, la précision de l'intensité du champ magnétique est de l'ordre de 1 nT. Cette imprécision élevée au regard de la précision

de chaque axe de mesure du capteur s'explique par l'impossibilité d'établir un étalonnage dans un champ magnétique parfaitement constant. Dans la pratique, l'information de diographie exploitée doit avoir une bien meilleure précision relative puisque le capteur reste pratiquement dans la même orientation vis-à-vis du vecteur magnétique local au cours de la diographie. Les erreurs d'inter-étalonnage entre les 3 axes du capteur ont alors peu d'incidence en mesure relative. La cadence de mesure programmée est de 357 Hz, pour permettre un éventuel filtrage des signaux 50 et 150 Hz provenant des réseaux électriques et pour s'affranchir des problèmes d'échantillonnage d'un multiple de ces périodes. Une telle correction n'est pas apparue nécessaire.

Le capteur immersible est inséré dans un étui protecteur en PVC et le câble de liaison est rigidifié par le recouvrement par un jeu de cornières en U emboîtées, en aluminium non magnétique, solidarisé à l'étui de protection (fig. 2). Une première cornière en U recouvre le câble, une seconde ferme le premier U en couvercle, la troisième vient perpendiculairement verrouiller le premier emboîtement. La structure est découpée en segments. La jonction de chaque couche est décalée afin de maintenir la rigidité de la structure. Le prisme 360° est fixé à l'extrémité haute de la perche ainsi constituée. La procédure de la diographie consiste à insérer la perche dans le puits de carottage à une vitesse de l'ordre de 1 m/s, puis à remonter le dispositif à la même vitesse en maintenant la verticalité de la perche. Cette opération, qui ne dure que quelques secondes, est reproduite 2 ou 3 fois par puits afin de pouvoir vérifier la répétabilité du signal de la diographie.

<figure><img src="image_82.png" /><figcaption>Figure 1</figcaption></figure>

Structure de portage du magnétomètre construit avec des tubes emboîtés en fibre de verre, en situation lors d'un test préliminaire à gauche et lors de la prospection du 9 mars 2015 à droite (avec station totale implantée en mode tracking sur le tell de Taraschina).
Le géopositionnement de toutes les mesures d'intensité de champ magnétique est réalisé, que ce soit pour la prospection 2DH ou 1DV, par poursuite laser du prisme réflecteur 360° par une station totale S8 Trimble. Cette dernière ne possède pas d'horloge interne. Elle est connectée à une tablette PC et pilotée par le logiciel Trimble Survey Manager TM. Les mesures d'angles sont réalisées à une cadence de 20Hz

tandis que la distance n'est déterminée qu'à une cadence de 10 Hz, et interpolée pour doubler la fréquence jusqu'à 20 Hz. La capture des mesures se fait par l'intermédiaire du logiciel HyperThermal qui va lire directement l'information que la station totale transmet à la Tablette PC sur le port USB. En raison de l'absence d'horloge interne, l'heure enregistrée n'est pas générée par la station totale, mais par la tablette.

<figure><img src="image_83.png" /><figcaption>**Figure 2**
Dispositif de diagraphie (Taraschina, 2015).</figcaption></figure>
Les processus de transfert des données numériques n'étant pas parfaitement en temps réel, une distorsion du temps enregistrée apparaît. Ce phénomène est révélé par les variations de pas de temps calculé entre deux mesures successives. Les mesures sont, a priori, réalisées à une cadence constante de 0.05 s. Des mesures sont absentes, soit parce que la poursuite n'a pas permis d'établir une position fiable, soit parce que le dispositif d'étalonnage interne de la station totale a bloqué la mesure. Ainsi, la distribution des pas de temps entre deux mesures successives devrait être sur des multiples de 0.05 s. Pour autant, la distribution constatée se rapproche d'une loi de Poisson, avec le mode sur 0.05 s avec des valeurs dès 0.01 s. Afin de disposer d'une information temporelle plus proche de la réalité de la mesure que celle obtenue, une procédure de correction est appliquée. La démarche consiste à attribuer une valeur temps qui suit un pas de 0,05s (pour maintenir la périodicité du phénomène, le pas est très faiblement inférieur à 0,05s) ou un multiple en cas de non mesure.

Le passage du géopositionnement local au géopositionnement global est réalisé en post-traitement à partir du relevé des bornes géoréférencées. Le système de projection utilisé est Stereo 70.

La représentation des données constitue une étape importante de l'exploitation des mesures. Généralement, les résultats d'une prospection géophysique sont représentés sur forme « image », c'est-à-dire de cartes où l'espace est divisé en unités élémentaires (la maille ou le pixel), auxquelles est attribuée une valeur (selon une échelle de couleur, par exemple). Une telle démarche nécessite d'interpoler les mesures existantes, de distribution spatiale non régulière, pour attribuer une valeur à chacun de ces éléments. La valeur attribuée est déterminée à partir des mesures voisines, y compris pour les unités élémentaires non couvertes par une mesure. La loi d'interpolation doit être robuste pour ne pas engendrer une distorsion de l'information. Cette approche nécessaire est d'autant plus approximative que la surface étudiée est sous-échantillonnée. Compte tenu des cadences de mesures élevées des instruments utilisés et des distances réduites entre deux profils, la densité spatiale des mesures tend vers une couverture exhaustive de l'espace. L'étape d'interpolation pour combler les vides n'étant plus nécessaire, il devient alors plus pertinent de représenter le nuage de points en 3D avec l'intensité du champ magnétique en 4e dimension selon une échelle de couleur. Cette démarche présente l'intérêt de conserver une information sur la variation de distance entre la surface du sol et le capteur avec l'altitude (z). Une alternative, beaucoup plus consom-

matrice en temps de traitement numérique, consiste à combler les vides entre chaque mesure en attribuant la valeur de la mesure la plus proche à l'espace considéré. Pour cela, la surface est découpée en polygones élémentaires dont chaque arrête constitue la médiatrice entre deux points de mesure voisins. Ainsi, chaque mesure est-elle associée à un polygone. Cette représentation, nommée diagramme de Voronoï, présente l'avantage de ne nécessiter aucune hypothèse d'interpolation. En interdisant les polygones de trop grande surface, donc non représentatifs car provenant de secteurs sous échantillonnés, seules les données ayant une réelle représentativité sont exploitées.

Les données de diagraphie viennent naturellement s'intégrer dans la représentation 4D. Les graphiques sont réalisés avec le logiciel Voxler 4.0 de Golden Software et Matlab pour la représentation du diagramme de Voronoï.

## 4. Prospection géomagnétique 2DH

L'étude du Tell a débuté en 2009 avec les premières prospections de surface et les premiers sondages. Dès 2011, Radu Dimitriu (GeoEcoMar) a conduit une prospection géomagnétique sur la moitié nord-est du tell (Dimitriu *et al.*, 2011). Du fait de l'utilisation de magnétomètres à proton, nécessitant une mesure statique, la prospection devient rapidement fastidieuse si l'on cherche à obtenir une haute résolution spatiale. Ainsi, lors de cette prospection, les mesures ont-elles été réalisées au nœud d'un carroyage métrique avec deux capteurs disposés respectivement à 0,75 et 2,75 m au-dessus de la surface du sol, soit 1 mesure par m² pour chaque capteur. Cela représente un peu plus de 4000 mesures sur une surface de l'ordre de 4000 m² pour chaque hauteur.

Pour comparaison, lors de la prospection que nous avons réalisée en 2015, la surface couverte en moins d'une journée avec le dispositif utilisé dépasse 7000 m² et a permis d'acquérir plus de 163 000 mesures par capteur. Cette rapidité d'exécution a permis une couverture quasi exhaustive du tell, mais a nécessité un défrichage préalable, plus laborieux, réalisé par toute l'équipe présente. Les zones basses périphériques n'ont pas été couvertes de manière régulière du fait de l'impossibilité de se déplacer sur ces surfaces humides sans portance. Dans la zone régulièrement couverte, plus de 20 mesures par m² ont été réalisées.

Lors de la prospection réalisée en 2011 par Radu Dimitriu, la correction de la variation temporelle du champ magnétique a été réalisée à l'aide de l'enregistrement de cette variation par un second magnétomètre mis
en base (statique) à une cadence de 1 mesure toutes les 15 secondes. Cette faible cadence, imposée par les capacités de mémorisation du magnétomètre, ne permet pas de corriger les variations pseudopériodiques de haute fréquence, c'est-à-dire de période inférieure à 30 secondes. Cette limite est similaire à celle du protocole utilisé pour le traitement des données de notre étude, utilisant les tendances de variation sur les cycles d'aller-retour. De ce fait, bien que réalisée avec des instruments et des protocoles différents, la correction de la variation temporelle entre la prospection de 2011 et celle de 2015 est de même nature. De ce fait, le niveau d'information est comparable entre les deux prospections.

Néanmoins, la hauteur de prospection au-dessus du sol diffère. En 2011 elle est constante, à 0,75 m du sol, tandis qu'en 2015 elle est plus faible, de l'ordre de 0,2 m, pour se rapprocher des sources, mais flucutue avec le balancement produit par le mouvement de la marche et l'évitement des amas de végétaux couchés au sol, imposant un relèvement du dispositif de portage à leur passage. En effet, la dynamique des signaux décroît en s'éloignant de la surface, de manière inversement proportionnelle à la distance à la source, élevée à une puissance entre 2 et 3 en général. Ainsi, pour la zone couverte en 2011, la carte des mesures réalisées à 2,75 m de la surface nécessite une échelle de moindre dynamique, environ 15 nT, par rapport à celle utilisée pour les mesures réalisées à une distance réduite à 0,75 m, qui dépasse 22 nT. L'échelle nécessaire pour représenter les données de la campagne de 2015, réalisée avec une distance à la surface réduite à 0,2 m est pratiquement décuplée, couvrant 150 à 200 nT de variation (bien qu'écrittant les valeurs extrêmes produites par les masses ferreuses tels que les piquets des bornes de géomètre). Pour autant, la proximité accrue des sources lors de la campagne de 2015 ne peut expliquer, à elle seule, l'augmentation de la dynamique du signal. Il s'avère que la couverture quasi exhaustive du tell en 2015 comprend des zones de très fortes anomalies, non explorées en 2011, origine essentielle de l'augmentation de la dynamique du signal.

Une autre différence entre les approches de 2011 et 2015 concerne la logique de représentation des données. L'échelle de couleur, non linéaire, utilisée pour représenter la variation de l'intensité du champ magnétique des mesures de 2011 (Dimitriu et al., 2011) vise à donner le même poids à chaque mesure. À l'opposé, l'échelle utilisée pour représenter les mesures de la campagne de 2015, présentées dans cet article, est fondée sur une répartition linéaire entre plusieurs seuils de couleurs. Elle a pour but de mettre

en valeur les anomalies extrêmes (rouge et bleu), potentiellement associées à des zones de feu, tout en discriminant les structures de plus faible contraste au sein des zones de variations plus modestes (en échelle de gris).

Quelles que soient ces différences, la géométrie des objets remarquables imagés à la fois en 2011 (à 0,75m du sol) et en 2015 (environ 0,2 m du sol) est comparable, démontrant la robustesse du signal géomagnétique. Il est à noter qu'il est nécessaire d'introduire une translation de l'ordre de 5 m vers l'ouest selon la direction de l'axe de carroyage du relevé de 2011 pour permettre une superposition des anomalies remarquables (telle que l'anomalie de la borne de géomètre B-02 dans l'angle sud-ouest par exemple).

La fouille du secteur 2 recoupe en partie nord une des anomalies majeures observées en 2011, non observable en 2015 du fait de la fouille archéologique alors réalisée (excavation et déblais). Cette fouille a mis au jour des blocs de terre cuite déstructurés pouvant être interprétés comme l'amoncellement volontaire de matériaux cuits en position secondaire. L'anomalie générée serait donc produite par les aimantations induites. Les aimantations rémanentes, de direction aléatoire, se compensant quant à elles dès que l'on s'éloigne de l'amoncellement. Il semble donc que la couverture de 2011 ne révèle pas de maison incendiée, puisque l'anomalie majeure observée ne correspondrait pas à un tel objet.

À l'opposé, la campagne de 2015 révèle 3 anomalies majeures dans le quart nord-est du tell, dans la zone non couverte en 2011 (fig. 3). Ces anomalies présentent une structure dipolaire avec une partie positive (figuée en blanc et rouge) plus vaste que la partie négative (figurée en noir et bleu) située au nord. L'intensité de ces anomalies apparaît au moins deux ou trois fois plus importante que l'anomalie majeure imagée en 2011. Deux de ces trois anomalies se trouvent en périphérie de la zone prospectée et ne sont pas imagées dans leur totalité. La première (A1), à l'ouest, est située sur le talus extérieur du tell principal correspondant au site 2. Cette zone plus basse, difficile à défricher à cause de la présence d'eau, n'a été que grossièrement couverte. La partie négative de l'anomalie, au nord, apparaît incomplète.

La seconde anomalie (A2) est située sur la marge nord du tell, dans la pente. La bordure nord de cette anomalie apparaît aussi incomplète (son pôle négatif étant partiellement tronqué), pour les mêmes raisons : cette zone est en limite de l'eau. Cette partie n'a pu être couverte par des profils aller-retour. L'opérateur
a avancé en crabe le long de la bordure émergé du tell, tout en maintenant les capteurs au-dessus de la zone ennoyée et en leur faisant faire un mouvement oscillant de gauche à droite pour couvrir au mieux la zone à portée du bras de la structure de portage. Les bordures est et ouest de l'anomalie positive de cette anomalie A2 apparaissent linéaires et parallèles. Une telle géométrie pourrait révéler des murs d'une construction incendiée.

La troisième anomalie (A3) est couverte dans son intégralité. Elle est implantée plus en hauteur, en bordure ouest du plateau du tell, en regard de la première anomalie située le plus à l'ouest. Cette troisième anomalie présente une structure principale rectangulaire (anomalie positive) à laquelle est adossée une seconde anomalie de moindre intensité, au sud/sud-ouest. Ces anomalies positives dessinent le plan d'un habitat, révélé par les vestiges des murs brûlés lors

<figure><img src="image_84.png" /><figcaption>Figure 3</figcaption></figure>

Prospections géomagnétiques de résolution inframétrique réalisée en 2015 avec l'emprise des sondages et zones de fouilles archéologiques antérieures et postérieures ainsi que l'emprise estimée de la prospection de 2011 (R. Dimitriu et al., 2011). Carte pré-traitée du capteur à 0,2 m au-dessus de la surface du sol. L'emprise de la prospection de 2011 est établie par corrélation des cartes représentant l'information proche de la surface. La différence d'altitude des prospections peut engendrer un léger décalage selon l'axe nord-sud des structures.
de la fouille (secteur 3). Les limites est et ouest de la partie positive de cette troisième anomalie présentent une structure crénelée. Cette apparence ne correspond pas à une réalité mais à un biais du protocole facilement explicable. Tout d'abord, cette forte anomalie est localisée en bordure du plateau et la prospection est réalisée selon des profils qui s'avèrent être sub-perpendiculaires à la structure. De plus, le dispositif de portage est conçu pour que les capteurs soient sur une même verticale avec le réflecteur 360°, qui sert de référence pour le géopositionnement des mesures.

Ainsi, un profil réalisé dans le sens aller, de l'est vers l'ouest, sera conduit en descente au niveau de cette anomalie. Si l'opérateur conserve l'horizontalité du dispositif, alors les capteurs s'éloignent de la source dans la descente et, de ce fait, mesurent un signal atenué. Si l'opérateur corrige la pente en rapprochant les capteurs du sol, alors c'est le géopositionnement qui est affecté, décalant la mesure artificiellement vers l'ouest.

À l'inverse, dans le sens retour de l'ouest vers l'est, la pente en montée impose de relever le dispositif pour ne pas talonner le sol avec le capteur bas. L'inclinaison du dispositif engendre un décalage des mesures vers l'ouest, tout comme dans le sens aller. L'attention de l'opérateur étant captée par la visée de la trajectoire de marche, il est difficile pour lui d'anticiper les changements de topographie du parcours. Arrivé en haut de la pente, puis sur le replat après le changement de pente, l'opérateur redresse le dispositif, en retard par rapport aux capteurs qui sont devant lui. Ainsi, les capteurs s'éloignent du sol pendant un laps de temps. Ces changements de distance capteur-source et de décalage du réflecteur sont la source des artefacts constatés. La figure 4 est une représentation en perspective 3D des données, corrigées de ces effets, mises en vis-à-vis de l'altitude de chaque mesure. En plus du retrait de chevronnage induit par le changement de pente, deux autres corrections sont appliquées :

1- la prospection a été réalisée en deux temps, la partie sud du site a été prospectée durant la matinée ;

<figure><img src="image_85.png" /><figcaption>**Figure 4**</figcaption></figure>

Représentation en vue oblique du nuage de points des mesures d'intensité du champ magnétique corrigée de la variation temporelle et des effets de décalage spatial des mesures et de la variation d'altitude des mesures. Les courbes de niveau sont celles produites par l'étude topographique du site, indépendamment de notre étude. Elles ont été extrapolées à partir d'un maillage de points topographiques.
la partie nord du site dans le courant de l'après-midi. La limite est visible au sud du secteur 2 (fig. 3). La variation temporelle du champ magnétique change de tendance au passage du soleil au zénith. Sur la figure 3, la correction de cette variation temporelle est réalisée sans tenir compte de ces deux périodes. Elle s'avère être imparfaite. Pour corriger ce décalage, une constante est introduite dans la valeur du champ mesurée en seconde période ;

2- de manière analogue, des profils obliques ont été réalisés pour passer d'un secteur à l'autre. Ils recoupent les profils subparallèles, par exemple visibles entre les deux bornes du quart sud-ouest sur la figure 3. Des positions de mêmes coordonnées sont donc couvertes deux fois, mais à des temps de passage différents. La correction de la variation temporelle appliquée ne corrigeant pas les variations de haute fréquence, son imperfection induit un décalage des valeurs. Les profils obliques, qui représentent les mesures réalisées en doublon, brouillent ainsi l'image. Ils sont supprimés.

Depuis cette prospection, un nouveau dispositif de portage a été conçu pour maintenir la verticalité des capteurs avec le dispositif de géopositionnement

(réflecteur 360° ou antenne GPS) et les repères visuels mis en place sont plus nombreux afin de faciliter le suivi du cheminement et permettre ainsi à l'opérateur de suivre plus fidèlement les variations de micro-topographie.

La texture de l'image magnétique obtenue (fig. 4 et 5) permet de distinguer plusieurs zones. Globalement, l'image peut être scindée en deux secteurs, grossièrement séparés par une limite d'axe NNE-SSW. À l'ouest de cette limite, les anomalies ont des dimensions caractéristiques supérieures à celles situées à l'est. Cette délimitation peut aussi être mise en relation avec la topographie, la partie la plus élevée correspondant assez bien avec la partie est. L'élévation de ce secteur et la présence des anomalies mineures structurées selon un axe NNE-SSW pourraient être la manifestation d'un remblaiement par des blocs de terre cuite provenant de maisons incendiées démantelées, tel qu'observé lors de la fouille du secteur 1.

À l'ouest, les trois structures d'habitat incendiées dominent la texturation de l'image. Néanmoins, les anomalies positives mineures situées sur la diagonale entre les deux structures incendiées du site 1 sont de

<figure><img src="image_86.png" /><figcaption>**Figure 5**</figcaption></figure>

Diagramme de Voronoï de la prospection 2DH. Chaque mesure est figurée par un point. Si la surface des polygones délimitée par les médiatrices avec les points voisins est inférieure à 0,25 m², alors l'espace autour du point est considéré de même valeur. Au contraire, si la surface est supérieure à 0,25 m², cas des relevés périphériques par exemple, alors l'espace est considéré être échantillonné à une trop faible densité spatiale de mesure pour permettre une extrapolation à l'espace contigu à chaque point de mesure.
dimensions comparables à celles des structures brûlées. Pour autant, il est difficile de pouvoir affirmer que ces anomalies correspondent à des murs d'habitat non incendiés, d'autant que la dimension d'une anomalie en surface dépend de la profondeur, de la dimension et de l'intensité de sa source. Les reliques d'une construction non incendiée située à une profondeur analogue à celle d'un bâtiment incendié devraient apparaître avec des dimensions plus modestes, ce qui semble le cas.

À l'opposé, la texture des limites nord et sud du secteur est du tell montrent des anomalies de petites dimensios, de contrastes élevés, mais d'intensité plus modeste que les zones incendiées. Ces anomalies, marquées par des variations spatiales rapides, traduisent la présence de sources magnétiques dans la proche surface. Compte-tenu du contexte, ces matériaux magnétiques peuvent correspondre à des tessons, des fragments ou des amas de terre cuite, ou tout autres matériaux magnétiques récents. La prospection géophysique a été réalisée après une campagne de prospection archéologique ciblée, l'ouverture de sondages et la réalisation de fouilles plus étendues.

Ces actions apportent des informations. Elles ont aussi pu laisser leur emprunte magnétique, telle les anomalies magnétiques produites par les tranchés de diagnostiques observables sur les images des prospections géomagnétiques réalisées à la suite de la remise en état du terrain.

La découverte de vestiges d'une habitation sub-contemporaine au sud, matérialisée par la présence de quelques briques de terre cuite et de fragments de restes métalliques (poêles à bois), apporte une explanation probable à l'origine de ces anomalies à fortes variations spatiales, présentes dans le secteur sud.

Quant aux anomalies analogues situées au nord-est du tell, en l'absence d'indices, seuls des blocs de terre cuite présents en sub-surface ou l'apport localisé de terre contenant des résidus de terre cuite pourraient expliquer une telle texture de l'image. L'âge de ces apports, contemporains, sub-contemporains ou associés à une éventuelle phase de ré-haussement du tell par des matériaux provenant de constructions incendiées, ne peut être déterminé par l'image produite.

L'imagerie géomagnétique effectuée depuis la surface ne donne que peu d'indications sur la profondeur des objets. Les carottages réalisés pour l'étude du tell ont permis de recourir à un dispositif de diographie, conçu pour acquérir une information dans les puits de carottage, juste après leur réalisation, avant qu'ils ne se referment naturellement. Au vu de la carte

prétraitée produite (fig. 3), deux transects sub-perpendiculaires ont été définis pour implanter les carottages (fig. 6, figuré bleu).

L'emplacement des carottages sur chaque transect a été déterminé de manière à être situé au milieu des minimums et maximums locaux des anomalies magnétiques. Quatre autres carottages ont été implantés sur la structure incendiée (A3) la plus centrale afin de déterminer les limites nord et sud des murs (fig. 6, figuré vert).

## 5. Prospection géomagnétique 1DV

Le dispositif de diographie, fondé sur une électronique de laboratoire et non de terrain, n'a été testé pour cette application qu'une seule et unique fois avant cette étude. Le paramétrage des cadences de mesure a été défini pour pallier à des pollutions provenant du réseau électrique sans savoir si cette pollution pouvait être significative. De même, la répétabilité du signal reste une inconnue. Pour cette raison, les diagraphies ont été répétées 2 ou 3 fois pour chaque puits, chaque diographie correspondant à un enregistrement décente puis remontée. Il s'avère que les interruptions d'enregistrement de la station totale constituent la source majeure d'absence de reproductibilité du signal, par absence de géopositionnement des mesures. Si des interruptions de l'enregistrement du déplacement par la station totale se reproduisent au cours des diagraphies successives d'un même puits, alors la diographie produite correspond à une courbe composite construite à partir des répliquas. Cette construction d'une diographie composite n'a pas été possible pour un seul puits, par absence d'enregistrement des coordonnées des mesures.

Pour vérifier la répétabilité, les diagraphies du transect WNW-ESE ont été refaites à un jour d'intervalle. Les différences significatives constatées concernent des puits dans lesquels des niveaux à blocs indurés (terre cuite) ont été traversés par le carottier. Dans ces cas, on observe que, lors de la descente ou de la monté du magnétomètre, son étui de protection (partie de plus gros diamètre du dispositif) a été accroché au niveau des blocs indurés. Cela se manifeste sur l'enregistrement de la station totale par un arrêt de déplacement. Des blocs ont donc très certainement été déplacés lors du déblocage du dispositif. Ces blocs étant très magnétiques et portant une aimantation rémanente élevée, leur déplacement est la source de modifications du champ magnétique dans leur voisinage (telles celles observées). Pour cette raison, les diagraphies du premier jour sont celles retenues, considérant que ce type de perturbations peut affecter le signal
après le premier passage. Quant aux autres puits, les erreurs de répétabilité traduisent un jeu de quelques centimètres selon l'axe vertical du capteur dans son étui, trop long de 1 ou 2 cm. Cet effet est détectable en comparant les digraphies réalisées en descente par rapport à celles réalisées en monté. En effet, le changement de contrainte sur le dispositif de diographie entre la descente – selon la gravité mais freiné par les parois (compression) –, et la traction pour la remontée (extension), peuvent entrainer une modification du câlage du magnétomètre dans son étui, n'entrainant au pire qu'un décalage de 1 ou 2 cm, non significatif.

La figure 7 présente les enregistrements successifs de l'une des digraphies réalisées dans un puits montrant des variations d'intensité de champ magnétique de faible dynamique (log 39, extrémité sud-est du tell). Les différences entre les 4 relevés sont mineures. De ce fait, sachant que ce relevé fait partie de ceux de faible variation et que la durée d'acquisition n'est que d'une dizaine de seconde, la correction de la variation temporelle du champ magnétique est totalement superflue et n'est pas réalisée. De même, la recherche

d'un signal associé au réseau électrique n'a pas été entreprise compte-tenu de l'absence de marqueur d'un tel phénomène.

Les courbes obtenues sont présentées dans la figure 8. Elles sont regroupées selon 3 transects en fonction de la position relative de chacun des puits telle que définie dans la figure 5. Les courbes provenant des zones de forte anomalie, qu'elles soient positives ou négatives, présentent les plus fortes amplitudes. Il faut noter que par l'effet de vide du puits de carottage, le sens de l'anomalie apparaît inversé par rapport à l'intensité d'aimantation relative de la source qui l'a produit. Ainsi, le passage dans un niveau très magnétique se traduira par une forte décroissance de l'intensité du champ magnétique, donc par une anomalie négative (Kristjánsson, 1975 ; Daly, Tabbagh, 1988). Afin de faciliter la lecture du signal, les abscisses sont en sens décroissant vers la droite, pour pouvoir traduire intuitivement le signal directement en terme d'aimantation des sources. Toutefois, dans le cas d'une maison brûlée en terre crue, les sources ne peuvent être considérées comme une couche magnétique

<figure><img src="image_87.png" /><figcaption>Figure 6</figcaption></figure>

Implantation des carottages sur le tell en fonction des anomalies magnétiques observées sur la carte prétraitée, interpolée, réalisée les jours suivants la prospection. En bleu : localisation des deux premiers transects ESE-WNW et WSW-ENE ; en vert : complément réalisé sur la structure incendiée (A3) la plus centrale.
d'extension infinie. Il faut aussi considérer des sources ponctuelles. Tout comme pour la prospection géomagnétique en surface (2DH), une source magnétique ponctuelle produira, dans un puits voisin de la source, une anomalie dipolaire. La forme et l'amplitude de cette anomalie du profil de diographie dépend de la distance, mais aussi de la position relative entre le puits et la source par rapport au nord magnétique. Les Log 33 et 43 présentent une courbe d'intensité du champ magnétique de forme générale très similaire. Celles-ci peuvent être interprétées comme le passage dans une couche très magnétique avec toutefois un épaulement en sens opposé sous l'anomalie principale. Cet épaulement traduit très certainement l'hétérogénéité latérale des couches traversées. Cette forme se retrouve aussi avec la courbe du log 36, mais sa partie supérieure apparaît beaucoup plus influencée par l'anomalie positive observée en 2DH. La localisation de ce puits, en position topographique basse dans la partie nord du tell, est peut-être l'explication de cette singularité (effet négatif du creux topographique). Ces 3 courbes proviennent des puits de carottage centrés sur les 3 anomalies majeures interprétées comme des bâtiments incendiés. Le pic correspond alors au niveau cuit.

Les courbes provenant du transect de la structure incendiée (fig. 8) ne reproduisent pas toutes ce schéma. Bien que localisée à l'extrémité sud/sud-ouest de l'anomalie positive, la courbe du log 44 ne présente qu'une faible dynamique. Toutefois, le pic le plus important se situe à la même profondeur que des pics majeurs des puits voisins de ce transect, situé entre 0,9 et 1,1 m de profondeur. Ce pic mineur pourrait donc traduire soit la présence de matériaux chauffés, remobilisés après leur chauffe (retombées de débris de paroi chauffée) soit un faible impact thermique. À l'extrémité opposée, les deux puits nord/nord-est montrent une plus forte complexité. En particulier, le log 41, le plus au nord-est, présente une structure dipolaire à une profondeur de l'ordre de 0,7 m au-dessus du niveau produisant la forte anomalie reconnue dans chaque puits. Ce puits a été implanté pour être situé en bordure externe de la zone incendiée. Il s'avère qu'il recoupe un niveau chauffé entre 1m et 1,1 m de profondeur approximativement, mais de modeste intensité. L'anomalie dipolaire moins profonde correspond probablement à un bloc basculé, tel un parement de mur cuit qui se serait effondré, en position déportée par rapport au puits.

<figure><img src="image_88.png" /></figure>

<figure><img src="image_89.png" /><figcaption>**Figure 7**</figcaption></figure>

Enregistrement de la profondeur relative et de l'intensité du champ magnétique au cours des deux descentes et remontées successives en fonction du temps dans le puits du log 39 (extrémité sud-est du tell). L'échelle de couleur utilisée pour représenter un point de mesure est proportionnelle à la durée entre deux mesures de profondeurs successives déterminées par la poursuite laser. Ce puits est situé dans une zone de faible dynamique de variation spatiale du champ magnétique, au milieu d'une faible anomalie positive. L'intensité légèrement plus élevée autour de 50 cm de profondeur pourrait révéler la source modeste de cette anomalie de surface.
<figure><img src="image_90.png" /><figcaption>Figure 8</figcaption></figure>

Diagraphies d'intensité du champ magnétique des puits de carottage. La couleur des courbes est celle de l'anomalie observée en surface, aux coordonnées du carottage.
La profondeur des niveaux incendiés par rapport à la surface est analogue pour les 3 structures. Sous ces niveaux incendiés, les pics observables sont de beaucoup plus faible amplitude, de l'ordre de 200 nT, à comparer aux 2000 nT du pic principal. Il serait tentant de rechercher une continuité entre ces pics dans l'idée d'identifier un niveau d'occupation du tell remarquable, tel un paléosol, ou un horizon correspondant à un épandage de matériaux associé à une utilisation du feu notamment. Cet exercice semble plus qu'hasardeux, sur la seule base des données de diagraphie géomagnétique, tant il n'est pas possible d'identifier une séquence qui se reproduise d'une courbe à l'autre entre les puits. Seule la partie inférieure, en dessous de 3,2 m de profondeur par rapport à la surface, montre une structure monotone, reflet d'un signal sédimentaire non anthropisé à attribuer aux læss. Ce niveau monotone (Logs 28 et 37) est surmonté d'une augmentation de l'intensité des propriétés magnétiques. Par comparaison avec les séquences des læss quaternaires de Chine, ce niveau correspondrait à un paléosol (un paléosol étant marqué par un enrichissement magnétique pédogénétilique). Les phases d'optimum climatique de la période des dépôts des læss ont favorisé le développement de paléosols. Parfois, le paléosol apparaît dédoublé (Log 30), démontrant que le plateau de læss ne présente pas une simple sédimentation tabulaire.

Tout comme il semble délicat d'interpréter les signaux des diagraphies de faible amplitude dans la partie surmontant les niveaux de læss, la relation entre les anomalies 2DH et 1DV des zones de faible dynamique n'est pas évidente. En effet, les carottages réalisés dans les secteurs de faible amplitude de la variation d'intensité du champ magnétique en surface (2DH, couleur du noir au blanc) présentent des signaux de diagraphie qui ne permettent pas d'identifier l'origine des anomalies de surface par la présence de matériaux plus ou moins magnétiques sous-jacents.

## 6. Prospection géomagnétique 2DH+ 1DV

La représentation 3D de l'ensemble des données permet une vue globale de la variation spatiale de l'information géomagnétique (fig. 9). Les sources magnétiques des 3 structures incendiées, identifiées en prospection de surface, sont bien recoupées en diagraphie. Toutefois, du fait de l'inclinaison du champ magnétique, la position des anomalies positives observées en surface est décalée vers le sud par rapport à la source (cas de l'hémisphère nord). Une démarche mathématique, dénommée la réduction au pôle, permet de recentrer l'anomalie positive à la

verticale de la source. Pour autant, cette démarche n'apporte pas une solution univoque. En effet, le résultat dépend des valeurs des variables utilisées, à savoir la déclinaison et l'inclinaison du champ magnétique qui est à la source des anomalies cartographiées. Les structures dipolaires sont numériquement recherchées selon l'axe définit par la déclinaison magnétique et chaque structure dipolaire est redressée selon la valeur d'inclinaison pour centrer l'anomalie. Ainsi, le pôle positif sera déporté en direction du pôle magnétique défini par la déclinaison introduite. Le déport sur cet axe sera d'autant plus important que l'angle de l'inclinaison avec la verticale (c'est-à-dire : 90°-Inclinaison, en référence à l'inclinaison au pôle) sera élevé.

Dans le cas d'une anomalie générée par une aimantation thermorémanente, le champ qui a créé cette aimantation thermorémanente diffère du champ actuel. Cette différence est due à la variation séculaire qui correspond à la dérive locale du champ magnétique. Ainsi, l'aimantation induite par le champ actuel, qui est proportionnelle à la susceptibilité magnétique des matériaux, et l'aimantation rémanente acquise au moment du refroidissement du matériau, ne sont pas strictement coaxiales. D'autre part, l'inclinaison et la déclinaison du champ magnétique local actuel ne sont pas mesurées lors de la prospection magnétique. Les modèles globaux permettent d'approcher cette information. Par exemple, pour le 9 mars 2015, date de la prospection, l'intensité du champ magnétique est de l'ordre de 48816 nT selon le modèle IGRF2020. Cette valeur est comparable à l'intensité moyenne mesurée, ~48780 nT. Selon ce modèle, l'inclinaison et la déclinaison de cette période sont de l'ordre de 62.8° et 5.9° Est, respectivement. L'intensité, l'inclinaison et la déclinaison magnétique du champ contemporain du refroidissement restent des inconnues. Au-delà de ces incertitudes, le modèle mathématique utilisé ne permet pas de distinguer le champ d'induction actuelle de celui de l'époque de refroidissement. Un champ moyen est alors utilisé.

La démarche consiste à produire différentes réductions au pôle, en recherchant une solution pour minimiser la dissymétrie des anomalies une fois réduite, en modifiant les paramètres. Une première étude a été produite avant toute obtention d'information complémentaire, en utilisant les paramètres du champ magnétique actuel fourni par les modèles de champ magnétique IGRF2020. Cette solution apparaît satisfaisante, réduisant notablement les anomalies négatives au NNW de la structure incendiée A3 (fig. 10b). Depuis, la fouille de cette structure centrale a été entreprise et a atteint les structures chauffées. Les relevées archéologiques permettent de localiser
la position des vestiges des matériaux chauffés. La confrontation de ce relevé avec la première réduction au pôle réalisée montre un décalage vers l'ouest de l'anomalie réduite par rapport aux relevés archéologiques. Fort de cette constatation, une nouvelle optimisation des valeurs d'inclinaison et de déclinaison a été entreprise en recherchant un meilleur accord entre le relevé archéologique et l'anomalie réduite. Avec une augmentation de la déclinaison de +4.5°, en conservant l'inclinaison actuelle, il y a un meilleur accord entre le pourtour du relevé archéologique et celui de l'anomalie positive (fig. 10c). Pourtant, l'épaulement négatif de l'anomalie au nord, qui devrait être atténué si la solution était meilleure, est accentué. Avec une augmentation la déclinaison de +9.5°, toujours en conservant l'inclinaison actuelle, l'accord avec la donnée archéologique est encore amélioré (fig. 10d) mais en accentuant le défaut d'atténuation de l'épaulement négatif de l'anomalie au nord. Cet antagonisme entre 1. l'optimisation du paramétrage

de la réduction au pôle pour réduire la partie négative de l'anomalie au nord et 2. la superposition du pourtour du relevé archéologique à celui de l'anomalie positive montre que la paléo-déclinaison du champ magnétique divergerait d'une dizaine de degrés avec le champ actuel. Cette amplitude est compatible avec les modèles de variation séculaire. Il faut aussi considérer que des blocs aient pu basculer. L'aimantation thermorémanente de ces blocs modifie localement le champ magnétique, rendant l'anomalie trop complexe pour répondre parfaitement à un modèle de réduction au pôle.

## 7. Conclusion

L'imagerie magnétique constitue un outil à privilégier pour étudier les structures incendiées en zone humide. Le protocole de diographie géomagnétique en puits de carottage, développé spécifiquement pour cette

<figure><img src="image_91.png" /><figcaption>**Figure 9**</figcaption></figure>

Représentation 3D des données géomagnétiques de surface 2DH et de diographie 1DV du tell. Les courbes topographiques, équidistantes de 0,1 m, sont surélevées de 0,5 m afin de délimiter le tell. Les puits sont matérialisés par une ligne blanche de 5 m de profondeur afin de pouvoir appréhender la profondeur relative de chacun d'entre eux.
application, a démontré son intérêt. Il permet
d'apporter une indication précise sur les profon-
deurs des vestiges chauffés, voir cuits, parfaitement
complémentaire à la prospection géomagnétique de
surface. En dehors de ces objets très magnétiques,
le constat est plus mitigé. Pour les zones de faible
anomalie, la confrontation des données d'imageries

géomagnétiques de diographie et de prospection de
surface n'apporte pas d'explication aux anomalies
observées en surface. La densité spatiale des diagra-
phies, déjà élevée, est certainement trop faible au
regard de l'hétérogénéité spatiale des dépôts très for-
tement anthropisés. Ce point constitue très certaine-
ment la limite de la méthode. En effet, augmenter la

<figure><img src="image_92.png" /><figcaption>Figure 10</figcaption></figure>

Verticalisation de l'anomalie géomagnétique par réduction au pôle. Les carrés verts matérialisent les limites de la fouille. Le trait rose correspond à la limite des matériaux chauffés reconnus à la fouille. A : données interpolées avec une maille carrée de 0,2475 m de côté. Les secteurs non couverts de manière régulière ont été exclus ; B : réduction au pôle avec une déclinaison magnétique de 5,9° Est et une inclinaison de 62,8° correspondant aux caractéristiques du champ magnétique local à la période de la prospection ; C : réduction au pôle avec une déclinaison magnétique de 10,4°E et une inclinaison de 62,8° ; D : réduction au pôle avec une déclinaison magnétique de 15,4°E et une inclinaison de 62,8°.
densité spatiale de l'information (ce qui reviendrait à
réaliser des carottages très proches les uns des autres)
irait à l'encontre de la démarche mise en place : la
moins invasive et la moins destructrice possible.
D'autre part, du fait de la profondeur des vestiges,
l'imagerie de surface ne permettrait pas d'augmenter
notablement la précision de l'information, même en
décuplant la densité spatiale des mesures, comme cela a
été développé pour l'étude des sols des grottes préhis-
toriques. Les protocoles d'imagerie 3D à résolution
infra-décimétrique sont transposables en plein-air,
mais sur des surfaces limitées (une centaine de m²).
Ils ont déjà fait leur preuve sur des foyers à galets
(Lévêque et al., 2019), des sites gaulois de production
de sel par méthode ignigène (non publié) ou encore
des épaves antiques (Lévêque et al., 2016). Une in-
formation de meilleure qualité pourrait alors être ob-
tenue, ce qui permettrait d'optimiser les traitements
de réduction au pôle. En revanche, ces protocoles
nécessitent, avec la technologie actuelle, un matériel
conséquent. La mise en œuvre d'une telle démarche,
susceptible de localiser l'implantation d'une fouille,
apparaît démesurée. Par contre, elle peut constituer

un des volets d'une étude cherchant à reconstituer
l'histoire thermique d'un objet. Pour cela, l'imagerie
géomagnétique de surface et en puits doit être cou-
plée à une étude paléomagnétique des matériaux
chauffés, pour déterminer les blocs encore en place de
ceux basculés post refroidissement. L'imagerie à ré-
solution infra- décimétrique permet de faire le lien
entre la prospection géomagnétique conventionnelle
et l'analyse paléomagnétique sur échantillon. Ce lien
est réalisé par un changement d'échelle tout d'abord
– du site à la structure incendiée – mais aussi du bloc
prélevé à l'échantillon qui en est extrait (Lévêque et
al., 2019). En zones humides – comme dans le delta
du Danube – en dehors du contexte des structures in-
cendiées ou des structures de chauffe, les méthodes
magnétiques sont peu propices à l'imagerie de sites
archéologiques, du moins tant qu'il n'y a pas d'ap-
port ou de création de matériaux magnétiques. En
effet, le contraste entre les sols, mêmes anthropisés,
et les sédiments est faible. Il est donc préférable de
privilégier les secteurs d'étude où il y a suspicion de
présence de structures incendiées ou structures de
chauffe, dans le premier mètre sous la surface.

Références bibliographiques

Daly, Tabbagh 1988 : DALY (L.), TABBAGH (A.) – Towards the in situ measurement of the remanent magnetization of oceanic basalts, *Geophysical Journal International*, Volume 95, Issue 3, December 1988, p. 481–489, https://doi.org/10.1111/j.1365-246X.1988.tb06698.x

Dimitriu et al. 2011 :
DIMITRIU (R. G.), CAROZZA (L.),
CONSTANTINESCU (A.M.), MICU
(C.), POP (I.C.), FLOREA (M.) –
*Shallow subsurface geophysical
investigation of a Chalcolithic site
within the Danube Delta.*
Unpublished. https://doi.
org/10.13140/RG.2.1.4610.4408

Kristjánsson 1975 :
KRISTJÁNSSON (L.) –
Geomagnetic measurements in drill
measurements in drill holes through
layered strata, *Geoexploration*, Vo-

lume 13, Issues 1–4 , p. 45-55, doi.
org/10.1016/0016-7142(75)90005-8.

Lecoanet et al., 1999 : LECOANET (H.), LÉVÊQUE (F.), SEGURA (S.) – Magnetic Susceptibility in Environmental Applications: omparison of Field Probes. *Physics of the Earth and Planetary Interiors*, 115, p. 191-204. https://doi.org/10.1016/S0031-9201(99)00066-7

**Lévêque et al. 2016 :** LÉVÊQUE (F.), REMAZEILLES (C.), MINJACQ (M.), PLASSON (N.), MATHE (V.), SANCHEZ (C.), JEZEGOU (M-P.) – Magnetic wood impregnated with corrosion products of ghosts nails: the case of ancient shipwrecks. 15<sup>th</sup> Castle Meeting New trends on Paleo, Rock and Environmental Magnetism, 21-27 august 2016, Dinant, Belgique, p. 91-92.

**Lévêque et al. 2019 :** LÉVÊQUE (F.), SOHN (M.), BURENS (A.), GAY (M.) – Fonction des structures de combustion à galets chauffés néolithiques du sud de la France : quelles informations peuvent apporter les méthodes magnétiques ? *Unpublished in L'archéométrie dans tous ses états : Perspectives inter-disciplinaires sur une discipline en constante évolution*. XXIIe colloque du GMPCA,9 au 12 mai, 2019, Université de Montréal, Montréal, Québec, Canada, Poster 17.

Olsen et al. 2001 : OLSEN (N.), TOFFNER-CLAUSEN (L.), RISBO (T.), BRAUER (P.), MERAYO (J.), PRIMDAHL (F.), SABAKA (T.) – In-flight calibrationmethods used for the Ørsted mission. ESA SP on Space Magnetometer Calibration. doi=10.1.1.41.3567
<figure><img src="image_93.png" /></figure>
# Chronologie des occupations et des grandes étapes de l'édification du tell de Taraschina

Laurent CAROZZA

Depuis 2010 et le début des recherches de terrain, nous avons multiplié les mesures d'âges de manière à ten- ter d'établir une chronologie des grandes étapes de la construction du tell et de ses occupations postérieures à l'époque chalcolithique. Les datations radiocarbone obtenues pour différents ensembles archéologiques mis au jour sur le site nous donnent aujourd'hui la possibilité de mieux appréhender la chronologie de l'occupation chalcolithique du tell, comme celle des différentes occupations postérieures à l'abandon de l'habitat. Ces datations absolues permettent d'identi- fier des rythmes d'occupations en lien avec les chan- gements environnementaux propres à la zone centrale du delta du Danube. Nous proposons ainsi, sur la base de l'analyse des âges mesurés sur les carottes, une synchronisation entre occupations humaines et phases de transformation des hydrostèmes.

## 1. L'échantillon

Le tell de Taraschina a fait l'objet de 26 datations radiocarbone (fig. 1). Compte tenu des particularités du site, les datations ont été réalisées tant sur des échantillons prélevés dans les différentes zones de fouille (21 dates) que dans les carottes prélevées sur le site (5 dates). De plus, nous disposons de 18 da- tations réalisées sur des carottes prélevées hors-site, utilisées pour synchroniser les évolutions de l'habitat et de l'environnement.

La zone 1 a fait l'objet de 6 dates, dont une sur carottage (Log 12) ; la zone 2 a été plus largement échantillonnée avec 14 dates uniquement obtenues sur des échan- tillons prélevés dans des ensembles archéologiques. Deux dates ont à ce jour été obtenues pour la zone 3, sur des horizons de surface. Quatre datations ont été effectuées sur carottage, notamment dans le but d'établir une chronologie de la formation de l'habitat.

L'essentiel des datations a été obtenu sur restes osseux (20 échantillons), d'autres ont été réalisées sur charbons de bois prélevés dans des carottes (5 échantillons) ou sur carporestes (graines de céréales carbonisées). Les restes osseux datés ont été sélec- tionnés en fonction de leur taille (de manière à mini- miser le risque de contamination par migration) et de la nature du taxon. Nous avons privilégié les grands ruminants (bovins, équidés) au détriment des suidés. Chez les suidés, une alimentation partiellement com- posée de mollusques (Unio) ou de poissons aurait pu influer sur le rapport isotopique du carbone mesuré. Cette réserve vaut également pour les datations réa- lisées sur les restes humains mis au jour dans la zone 2. La diète des individus a pu influer sur les résultats obtenus.

À titre exploratoire, nous avons réalisé une datation sur une mandibule de micro-vertébré prélevée dans une structure soumise à des bioturbations (Fs 1134). Le résultat obtenu (1650 ± 30 BP, soit 264-533 AD) nous a conduit à écarter cette datation de l'analyse.

La distribution stratigraphique des échantillons datés montre que nous avons privilégié la base et le som- met de la séquence stratigraphique. La sélection des échantillons a, de fait, été dictée par la problématique de la recherche, c'est-à-dire la datation des fosses et des niveaux d'occupation formant le sommet de la séquence stratigraphique (établissement de la data- tion de l'abandon), et par la volonté d'approcher la chronologie de la formation du tell en sélectionnant des échantillons à la base des carottes. Nous avons également tenté de répondre à des questions relatives au statut du site, notamment en lien avec la présence de restes humains disséminés dans la partie sommitale de la séquence stratigraphique. Par ailleurs, la datation de restes d'équidés a été dictée par une
problématique en lien avec le développement de cette espèce dans les horizons chalcolithiques de la culture Gumelnița.

## 2. Rythmes des occupations et des abandons du tell de Taraschina

Nous allons dans un premier temps prendre en compte l'ensemble des 25 datations dans le but d'appréhender les rythmes de l'occupation du tell : de sa formation à ses multiples réoccupations.

### 2.1. Chronologie des occupations du tell chalcolithique

L'occupation chalcolithique du tell de Taraschina a été datée via 18 datations radiocarbone (fig. 1). Les bornes chronologiques maximales sont comprises entre 4800 et 4050 BC, soit une durée maximale d'un peu plus de 700 ans. La datation qui fonde la formation du site provient de la carotte 4 ; elle a été obtenue sur un reste carpologique carbonisé recueilli au contact du sol naturel : 5840 ± 40 BP, soit 4798-4584 BC (95,4 %) 4780-4621 (68,2 %). Une seconde date a été réalisée à la base de la carotte 33, à la cote -3,83m : 5690 ± 30 BP, soit 4604-4456 (95,4 %) 4548-4466 (68,2 %). Cette datation se rapproche de celle obtenue pour la base de la carotte 12, réalisée dans la zone 1 : 5680 ± 30 BP, soit 4590-4453 (95,4 %) 4541-4466 (68,2 %). Ces données sont encore trop indigentes pour établir un modèle solide, mais tout porte à croire que la formation du tell est intervenue durant la première moitié du 5e millénaire BC, probablement entre 4800 et 4600 BC.

La partie inférieure de la stratigraphie reste pour l'instant non datée. L'essentiel des datations a été obtenu pour les ensembles archéologiques fouillés dans la moitié supérieure de la séquence stratigraphique du tell (fig. 1) ; elles se rapportent principalement à l'intervalle 4550 et 4250 BC (13 dates). Trois dates obtenues sur les Us 2047 (4520-4356 BC), 2045 (4520-4356 BC) et la structure 2054 (4505-4351 BC) se rapportent à la seconde moitié du 46e siècle avant notre ère et la première moitié du 44e siècle BC. Un lot de 6 autres datations caractérise l'intervalle de temps compris entre 4450 et 4350 BC. Ces dates ont été obtenues dans la zone de fouille 1, sur les Us 1102 (4461-4344 BC), 1107 (4487-4348 BC) et 1206 (4456-4335 BC) qui forment des épandages ou des niveaux de remblais. La datation obtenue pour le silo 1134, structure creusée dans les niveaux précédemment datés, est légèrement postérieure (4415-4282 BC). Dans la zone 2, la datation de l'Us 2038

(4455-4330 BC), formant un niveau d'épandage, est contemporaine de la dernière occupation de la zone 1. Ces occupations sont statistiquement contemporaines du démantèlement d'une unité d'habitation incendiée en cours de fouille dans la zone 3 (Us 3004 - 4450-4335 BC). Ces 6 datations, bien qu'issues de contextes stratigraphiques différents, matérialisent toutes la dernière étape de l'édification du tell, dans le sens où elles ont été obtenues sur des unités stratigraphiques formant des paléosols ou des épandages mis au jour dans la partie sommitale de l'habitat. Ces niveaux sont ponctuellement recoupés par des fosses et des silos creusés dans ces horizons et dont la caractéristique est de présenter une extension large. Dans la zone 2, le silo 2015 (4444-4275 BC) et l'Us 2009 du silo 2010 (4325-4052) sont postérieurs à la formation de ces niveaux du tell principal. La datation obtenue pour l'Us 2044 (4368-4236 BC), unité qui se rapporte à l'ensemble stratigraphique 5, a été rejetée car stratigraphiquement largement antérieure. L'unité se situe à environ 0,85 m sous la surface du tell. L'échantillon daté résulte très vraisemblablement d'une contamination en lien avec le creusement des fosses silos.

Les deux dates obtenues pour le site 2 (qui forme une excroissance adossée au nord du tell principal) se rapportent à la même tranche chronologique. Elles ont été obtenues dans des carottes : la première à 2,18 m sous la surface du sol dans la carotte 33 (4368-4246 BC) ; la seconde à 1,45 m dans la carotte 19 (4336-4076 BC). Si ces deux datations indiquent une concomitance de l'occupation entre le tell principal (site 1) et le site 2, tout porte à croire que l'occupation se serait développée plus longuement sur le site 2. La plate-forme du tell principal n'aurait été quant à elle occupée que de manière épisodique, ou à d'autres fins que celle d'un habitat. Cette donnée soulève la nécessité de s'interroger sur la chronologie de l'abandon général du tell en tant que site d'habitat. Deux datations s'écartent légèrement de ce modèle. Elles ont été obtenues pour une fosse-silo : Fs 2010 : 5350 ± 35 BP, soit 4323-4052 (95,4 %) 4313-4073 (68,2 %) et pour la carotte 19 prélevée sur le site 2 - Log 19 -1,45 m 5350 ± 35 BP, soit 4323-4052 (95,4 %), 4325-4236 (83 %), 4313-4073 (68,2 %). Ces deux datations se singularisent par une large plage temporelle après calibration, due à l'effet plateau qui caractérise les derniers siècles du 5e millénaire avant notre ère. De fait, l'intervalle de temps couvert par ses deux dates est long et compris entre 4325 et 4050 BC. La chronologie de l'abandon du tell peut, d'un point de vue théorique, être estimée entre 4250 et 4050. L'analyse de l'assemblage mobilier de la structure 2010 devrait permettre de préciser cette datation par l'examen des critères typologiques de la céramique.
<figure><img src="image_94.png" /><figcaption>**Figure 1**
Ensemble des dates réalisées sur le site de Taraschina.</figcaption></figure>
## 2.2. 4100-2900 BC : un long hiatus avant une occupation à caractère funéraire ?

On observe un hiatus, long de plus d'un millénaire (4050-2900 BC), entre la fin de l'occupation chalcolithique et la ré-occupation de la surface du tell à des fins funéraires durant l'âge du Bronze. La mise en évidence d'ossements humains, découverts dans les horizons supérieurs de la zone 2 (Constantinescu ce volume), nous a incité à lancer une campagne de datations radiocarbone dans le but de préciser la datation de ces restes épars.

Deux ossements mis au jour dans l'Us 2005 ont été datés : 4135 ± 35 BP, 2873-2588 cal BC (95,4 %) – 2873-2619 cal BC (93%) – 2862-2631 cal BC (68,2%) ; 4190 ± 35 BP 2892-2666 cal BC (95,4 %) /2892-2835 cal BC (24,3%) / 2817-2666 cal BC (71,1%), 2884-2700 cal BC (68,2%). Ces deux dates sont très cohérentes et montrent que des dépôts funéraires, dont d'autres ossements ont été découverts à l'interface entre le premier paléosol chalcolithique et les niveaux de mise en culture Us 2001 et 2002, se rapportent à la première moitié du 3e millénaire avant notre ère. Us 2005, métatarsien 4, gauche a été daté de 3835 ± 40 BP, soit 2460-2151 cal BC (95,4 %) – 2460-2198 cal BC (93,2%) / 2396-2205 cal BC (68,2%) – 2436-2205 cal BC (65,9%). Cette date se rapporte sans conteste à l'âge du Bronze ancien mais correspond à la phase finale de cette période. Elle est cohérente avec la date obtenue sur un prélèvement effectué sur une connexion anatomique d'un ulna et d'un radius mis au jour dans l'Us 2042. Le résultat obtenu se rapporte également au Bronze ancien : Us 2042, ulna et radius gauches d'un même individu. Ce dernier a été daté de 4045 ± 35 BP, soit 2836-2473 cal BC (95,4 %) – 2670-2473 cal BC (91%) – 2621-2491 cal BC (68,2%). Ces ossements correspondent très probablement à une sépulture en fosse, datée du Bronze ancien, mise au jour dans la zone 2. Ce dépôt est situé le long de la berme de la fouille et n'a pu être intégralement fouillé à ce jour. Les quatre datations attestent d'une fréquentation de la zone de Taraschina entre 2900 et 2200 cal BC et de l'implantation de probables monuments funéraires durant le Bronze ancien.

## 2.3. Une réoccupation à la fin de l'âge du Bronze

Nous avons également identifié une probable occupation du site à la fin de l'âge du Bronze. Trois dates ont été réalisées sur des faunes domestiques mises au jour dans l'Us 2030, à la surface du remplissage de la

fosse 2010 et au sein de l'Us 3002. Il s'agit d'éléments intrusifs au sein de structures attribuées à l'occupation chalcolithique, culture de Gumelnița.

Toutefois, ces éléments hors contexte révèlent de manière indéniable que la surface du tell a accueilli des occupations temporaires datant de l'âge du Bronze, bien qu'aucune structure ni aucun mobilier céramique caractéristique n'ait pu être identifié. La première datation est issue d'un ossement mis au jour au sein de l'Us 2030 : 3320 ± 30 BP, soit 1684-1523 cal BC (95,4 %) 1634-1532 cal BC (68,2%). La seconde date a été obtenue pour la structure Fs 2010, 3080 ± 35 BP soit 1427-1260 cal BC (95,4 %) 1406-1296 cal BC (68,2%). L'échantillon daté pour l'Us 3002 a été prélevé au contact de la surface cultivée et des premiers niveaux archéologiques in situ : 3070 ±30 BP, soit 1415-1236 cal BC (95,4 %) 1393-1278 cal BC (68,2%). Ce résultat est compatible avec celui obtenu pour la fosse Fs 2010. Ces trois mesures témoignent de la fréquentation du site durant le Bronze moyen, dans un intervalle de temps compris entre 1750-1200 BC.

## 2.4. Des occupations aux périodes historiques

L'analyse des mobiliers recueillis lors de la fouille des niveaux superficiels et des prospections permet de jeter un regard neuf sur les réoccupations du tell durant les périodes historiques.

Deux monnaies datant du début du second siècle de notre ère (vers 230 AD) ont été mises au jour sur le site (voir les études correspondantes dans Ailincăi et al., et Nuțu ce volume). L'analyse des mobiliers isus des zones 1 et 2 a permis d'exhumer quelques fragments de céramique byzantine. Ces éléments permettent d'attester qu'une fréquentation de la butte artificielle que formait le tell est intervenue entre le 10e et le 11e siècle de notre ère. Par ailleurs, une pipe en terre cuite, de facture ottomane, montre que le site a été visité au 19e siècle, alors qu'aucune structure archéologique contemporaine n'a été identifiée. Enfin, la mise au jour de restes humains dans les niveaux de décapage et en prospection nous rappellent que le toponyme « movila » que l'on peut traduire par cimetière, pourrait indiquer la présence de sépultures de russes Lipovènes dont nous savons qu'ils ont peuplé la zone à la fin du 18e siècle.

Plus anecdotiques, les enquêtes orales indiquent que le tell de Taraschina a été occupé durant la période communiste, dans les années 1960, notamment en tant que plateforme de stockage pour les roseaux.
### 3. Datation des carottes et synchronisation avec les occupations archéologiques

Une dizaine de carottes ont été prélevées dans le proche environnement du tell. Quatre d'entre elles ont fait l'objet de datations ¹⁴C (14 dates). Ces carottes sont, pour certaines, très proches de l'habitat chalcolithique, distant de 25 mètres (carotte 24) à 45 m (carotte 21). Ces séquences sédimentaires, prélevées à l'interface entre l'habitat *stricto sensu* et les milieux dits « hors site », sont soumises à l'influence des occupations humaines. Lors des phases d'abandon du tell, ces mêmes séquences décrivent des dynamiques davantage naturelles. Pour échapper à cet effet, nous avons réalisé une carotte le long de l'actuel canal Şontea, distant de près de 350 m du tell. Ces datations ne peuvent toutefois pas être utilisées de manière directe pour caractériser les occupations humaines du tell de Taraschina (Reingruber, 2015).

Les datations radiocarbone ont été réalisées, pour l'essentiel, sur du matériel carbonisé (graines et micro charbons de bois) ou bien sur restes végétaux imbibés (feuilles et bois). De toute évidence, ces matériaux proviennent pour partie des activités humaines réalisées sur le site et dans son environnement proche. Cet argument nous a parfois été opposé (Carozza et al., 2011) pour rejeter la validité chronologique des séquences analysées. Lorsque l'on met en relation les dates avec chacune des séquences (fig. 3), on observe que toutes les dates présentent une hiérarchisation chronologique cohérente, sans inversion manifeste. Ce seul argument suffit à montrer que s'opèrent multimentionnément, à compter de la seconde moitié du 5e millénaire BC, la construction du tell (occupations chalcolithiques et occupations postérieures) et l'édification de la plaine d'inondation dans l'environnement proche du site archéologique.

La carotte 24, prélevée au sud du tell, a livré une séquence de plus de 3,50 m. La base de la carotte est caractérisée par la présence de niveaux contemporains de l'occupation du tell. La première datation (-3,46 m) est comprise entre 4462 et 4338 cal BC ; la seconde, réalisée sur un échantillon prélevé à -3,20 m, est légèrement postérieure : 4342-4242 cal BC. La troisième date est issue d'un échantillon prélevé à -2,33 m. Son intervalle est compris entre 2009 et 1772 cal BC, moment où aucune occupation du tell n'a été identifiée.

La carotte 10b, située au nord ouest du tell, a fait l'objet de 5 datations (Carozza et al., 2011). La particularité de cette séquence est d'avoir livré 5 dates qui couvrent une plage de temps de moins de 600 ans, comprise entre 4592 et 3997 BC (fig. 2). Parmi ces 5

dates, qui décrivent une dynamique de sédimentation très rapide (-1,35 m en 6 siècles), l'une d'elles s'écarte légèrement du modèle (-1,92/1,90 m - 4504-4351 cal BC). Si l'on rejette cette datation, on observe qu'une grande portion de la séquence est contemporaine de l'occupation du tell. Les carottes 10b et 24 présentent la même caractéristique de voir se développer sur le substrat de loess, durant la seconde moitié du 5e millénaire avant notre ère, de puissants niveaux hydromorphes. Cette chronologie de mise en place d'une plaine d'inondation est confortée par l'obtention d'une autre datation pour la carotte 21 : 5545 ± 35 BP soit 4452-4341 cal BC (95,4 %). L'ensemble de ces données montre que des changements environnementaux très rapides s'opèrent entre le 45e et le 44e siècle avant notre ère.

La carotte 14 a fait l'objet de 6 datations qui permettent de synchroniser changements environnementaux et rythmes de l'occupation humaine. La base de la séquence, à l'instar des autres séquences, montre le développement de niveaux hydromorphes durant la seconde moitié du 5e millénaire BC (Log 14 : - 4,87 m, 4258-4001 cal BC). Cette mesure d'âge est pour partie contemporaine de la fin de l'occupation et du moment de l'abandon du tell chalcolithique. Deux autres dates ont été effectuées dans la tranche basse de la carotte (- 3,73 m 3003-2876 cal BC ; - 3,41 m 2281-2038 cal BC ; - 2,61 m 2131-1912 cal BC). Ces trois mesures encadrent les phases de réoccupation du tell durant l'âge du Bronze. Mais ce qui caractérise la base de la séquence, c'est la dynamique de sédimentation très rapide. On observe en effet que 2,25 m de sédiments se sont accumulés dans cette zone basse en un peu plus de deux millénaires. Un autre fait marquant est la mise en évidence d'un hiatus sédimentaire, long de plus de 3 millénaires, durant lesquels aucun événement géomorphologique n'a été enregistré. La reprise d'une dynamique rapide de formation de la plaine d'inondation s'opère à compter de la seconde moitié du 11e siècle de notre ère (fig. 2). Durant ce hiatus sédimentaire, on observe la présence d'occupations humaines dans la zone de Taraschina au cours de la période protohistorique (second âge du Fer), de l'Antiquité (3e siècle de notre ère) et de la période Byzantine (10e - 11e siècle).

### 4. Chronologie des occupations et des grandes étapes de la construction du tell

Après avoir pris en compte l'ensemble des datations pour observer le rythme des occupations relatives à la période Chalcolithique et à l'âge du Bronze, nous avons réalisé une seconde analyse prenant davantage
<figure><img src="image_95.png" /></figure>

<figure><img src="image_96.png" /></figure>

<figure><img src="image_97.png" /><figcaption>Figure 2</figcaption></figure>

Synchronisation des occupations humaines de la zone proche du tell de Taraschina et des principales séquences sédimentaires prélevées hors site.
en compte les contextes stratigraphiques. Lors de la
fouille, nous avons opéré des regroupements d'uni-
tés stratigraphiques au sein d'ensembles cohérents.
Ces ensembles peuvent réunir des groupes de fosses,
différents horizons formant une unité d'habitation ou
bien des niveaux d'abandon. Nous avons ainsi opéré
des regroupements au sein de chaque zone de fouille,
puis proposé des équivalences qui permettent d'ap-
préhender le site dans sa globalité.

4.1. Hiérarchisation des datations de la séquence chalcolithique

Le traitement des données a été réalisé de manière clas-
sique, en calibrant les dates (intervalle de confiance de
95 %) et en établissant une représentation graphique
qui montre leur hiérarchisation, sans tenir compte des
données stratigraphiques. Nous n'entrerons pas dans
le détail des données, telles que discutées précédem-
ment, mais nous nous intéresserons aux principales
tendances. Les incertitudes relatives à la chrono-
logie de la genèse du tell ne sont pas levées. Nous
disposons d'une seule et unique datation antérieure à
4600 BC pour la base de la carotte 4. Toutes les autres
dates sont postérieures et cette question ne sera pas
tranchée avant obtention de nouveaux résultats. Les
deux dates obtenues pour la base des carottes 12 et
34 sont en revanche très cohérentes. Toutefois, la
moitié inférieure de la séquence stratigraphique, dont
il faut rappeler qu'elle dépasse 3,5 m par endroits, de-
meure méconnue, notamment en raison de l'absence
de matériel disponible dans les carottes pour réaliser
des datations. En l'absence de restes carbonisés et
d'ossements de mammifères, nous nous trouvons
face à une impasse. À l'heure actuelle, nous tentons
de valider la faisabilité de l'utilisation de méthodes
géo-chronologiques alternatives, sachant que les
échantillons proviennent de carottages et ne peuvent
pas être orientés.

La moitié supérieure du tell est, comme nous l'avons
vu, très bien échantillonnée dans l'espace et en fonction
des contextes (fosses, silos, sols, remblais...). Nous
disposons de 15 datations, dont deux se rapportent au
site 2. La hiérarchisation de ces datations nous a per-
mis de distinguer quatre ensembles. Le premier est
constitué de datations obtenues dans la zone 2, à la
base de la séquence fouillée. Trois dates obtenues pour
les Us 2047, Us 2045 et la structure 2054 s'inscrivent
dans une plage de temps large, comprise entre 4500
et 4350 BC. Le groupe suivant est essentiellement
formé de fosses et de sols archéologiques fouillés
dans les niveaux supérieurs des zones 1 et 2. Ces da-
tations sont pour l'essentiel légèrement postérieures

à celles du groupe précédent, avec un intervalle de
temps compris entre 4450 et 4350 BC. L'unique data-
tion obtenue pour la zone 3, au sommet de la couche
d'incendie d'une maison, s'inscrit dans ce modèle
(4450-4335 BC). Seule la datation obtenue pour le
silo Fs 1134 (4415-4282 BC) s'écarte de cette phase
chronologique, bien qu'elle en soit archéologique-
ment contemporaine. Enfin, les dates obtenues pour
le site 2 sont toutes postérieures à 4370 BC. L'oc-
cupation de cet habitat, adossé au tell principal, se
serait développée entre 4350 et 4200 BC si l'on s'en
tient aux deux seules dates réalisées. Les habitats qui
forment cette partie du tell sont pour partie contem-
porains des fosses mises au jour dans la zone 2 du
site 1 (Fs 2015 et 2010).

Sur la base de ces données, nous pouvons proposer l'hypothèse que les niveaux d'habitation mis au jour sur le tell principal (zone 1, Us 1102, Us 1206, Us 3004) datent tous schématiquement d'une période antérieure à 4350 BC. L'occupation du site 2 s'amorcerait quant à elle postérieurement. Nous avons utilisé cet argument pour évoquer la possible existence d'une courte phase d'abandon du site autour de 4350 BC, rupture qui pourrait transparaître par l'absence de datations couvrant nettement l'intervalle de temps compris entre 4350 et 4250 BC (fig. 3). Ce phénomène pourrait également résulter d'un biais dans l'échantillonnage, lié à l'exiguïté des zones fouillées qui représentent actuellement moins de 5 % de la superficie totale du tell.

4.2. Histogramme cumulatif pondéré de l'occupation chalcolithique

Une autre manière de représenter la distribution des
datations radiocarbone, telles qu'elles représentent
la durée de l'occupation du tell, est de construire un
histogramme cumulatif pondéré (fig. 4). L'intérêt de
ce type de représentation est d'obtenir une vision
synthétique d'un ensemble de dates. Toute la diffi-
culté réside dans l'interprétation des tendances et de
la distribution. Trop souvent, nous avons tendance
à interpréter ces données en terme de dynamique,
alors qu'elles ne transcrivent que des tendances, tri-
butaires des contextes archéologiques datés. Dans le
cas spécifique du tell de Taraschina, il nous a semblé
interessant de voir où se situaient les ruptures dans la
distribution des données. Si l'on travaille avec un in-
tervalle de confiance à 95%, on observe une très forte
densité autour de la période comprise entre 4500 et
4300 cal BC. Ces données consolident ce que nous
observions précédemment, notamment les change-
ments qui s'opèrent autour de 4350-4300 BC dans
<figure><img src="image_98.png" /><figcaption>**Figure 3**
Hiérarchisation des dates de la séquence chalcolithique de Taraschina.</figcaption></figure>
l'organisation générale de l'habitat. La distribution de ces mêmes données avec une crédibilité de 68 % tend à accentuer cette même tendance, bien que l'on observe que la rupture s'opère plus tardivement, autour de 4250 BC. Ces données sont compatibles avec les modèles acquis récemment dans la plaine Roumaine, notamment sur le site de Pietrele, qui montrent que le déclin de l'occupation des tells s'opèrerait vers 4250-4200 cal BC (Hansen, 2007 ; Reingruber, 2015 ; Weninger et al., 2015). À Taraschina, l'hétérogénéité des données ne permet toutefois pas encore d'assurer le scénario historique selon lequel le déclin de l'occupation s'amorcerait vers 4350 BC, avec pour conséquence une recomposition de l'espace habité et de son environnement proche. Pour poursuivre cette analyse, nous avons tenté de développer une modélisation bayésienne des données, sachant que cette approche produit une vision « idéalisée ».

## 4.3. Traitement bayésien de l'ensemble des dates

Nous avons réalisé, dans une seconde étape, un traitement statistique des données en utilisant un modèle bayésien (logiciel chronomodel V1.3.4). L'intérêt de ce type de traitement réside dans l'introduction de contraintes dans le modèle, et notamment l'opportunité offerte par le logiciel d'insérer des liens stratigraphiques lorsque les unités datées (Us) posédent un lien d'antériorité direct ou déduit. Nous

avons en conséquence pris en compte, pour chaque zone de fouille, les liens stratigraphiques observés entre les unités stratigraphiques datées. Dans le cas d'un lien direct d'antériorité, nous avons introduit une contrainte. Dans le cas d'une succession déduite, nous avons organisé les événements en fonction de ce lien supposé d'antériorité. L'analyse stratigraphique de chaque zone de fouille a permis de mettre en lumière des ensembles cohérents (niveaux de sols, de remblais, groupe de fosses...) correspondant à des étapes de la construction du tell. Des contraintes stratigraphiques ont ainsi pu être établies entre certains de ces ensembles, notamment dans la zone 2. On observe toutefois que le lien stratigraphique introduit entre la fosse 1134 et l'Us 1102 tend à vieillir la datation de cette structure (credibility interval 95,0 % -4619-4333 BC). De la même manière, l'estimation de la durée de l'Us 3004 – correspondant au démantèlement de la maison incendiée de la zone 3 – a tendance à être légèrement plus longue (credibility interval 95,0 % -4510-4253 BC).

## 4.4. Modèle chronologique bayésien par phases

Bien que le nombre de datations disponibles pour caractériser la séquence du tell de Taraschina ne soit pas très important, il nous a semblé intéressant de présenter un modèle chronologique fondé sur le phasage chronologique tel qu'observé à la fouille en regroupant, lorsque cela était possible, les dates se rapportant à une phase (fig. 6). Nous avons ainsi intégré dans le modèle des contraintes stratigraphiques entre les ensembles 5-6, 4 et 2, déduits de l'analyse chrono-stratigraphique de la zone 2.

### - La phase initiale

La phase de fondation du tell a été déterminée à l'aide de 3 dates, toutes situées au contact entre des niveaux anthropiques et le substratum. Seules deux dates sont associées à des niveaux archéologiques. La date obtenue pour le Log 34 a été réalisée sur un micro-charbon de bois présent dans une couche naturelle, de type paléosol ou remplissage d'une dépression naturelle.

D'un point de vue statistique, nous pouvons avancer l'hypothèse que le début de cette phase initiale pourrait débuter entre 4893 et 4560 cal BC et s'achever entre 4595 et 4458 cal BC. Seule la date la plus ancienne, obtenue à la base de la carotte 4, situerait dans cette zone le début de l'occupation entre le début du 48e siècle BCE et la fin du 46e siècle BCE. Les deux autres datations sont plus cohérentes et situent entre 4600 et 4450 les premières occupations au nord du tell.

<figure><img src="image_99.png" /><figcaption>**Figure 4** Histogramme cumulatif pondéré de l'occupaton chalcolithique.</figcaption></figure>
**Duration (posterior distrib.) :**

MAP : 190 Mean : 209 Std
deviation : 89
Q1 : 146 Q2 (Median) : 204 Q3 : 262
HPD Region (95,0 %) : [34; 381] (95,0 %)
Credibility Interval (95,0 %) : [37, 380]

**Begin (posterior distrib.) :**

MAP : -4709 Mean : -4718 Std
deviation : 83
Q1 : -4765 Q2 (Median) : -4713 Q3 : -4662
HPD Region (95,0 %) : [-4884; -4553] (95,0 %) BC/AD
Credibility Interval (95,0 %) : [-4881, -4554]

**End (posterior distrib.) :**

MAP : -4508 Mean : -4509 Std
deviation : 38
Q1 : -4531 Q2 (Median) : -4508 Q3 : -4484
HPD Region (95,0 %) : [-4587; -4437] (95,0 %) BC/AD
Credibility Interval (95,0 %) : [-4585; -4438]

## - La phase 2 : ensembles 5 et 6 de la zone 2

La seconde phase correspond à des datations obtenues dans la zone 2, issues des ensembles 5 et 6, correspondant à une unité d'habitation et à ses structures connexes. D'un point de vue archéologique, cette phase est très cohérente. Parmi les quatre dates réalisées pour cet ensemble, trois témoignent de cette forte homogénéité (Us 2045 - 4488-4347 cal BC ; Us 2050 (Fs 2054) - 4505-4351 cal BC ; Us 2045 -4520-4356 cal BC). La quatrième date obtenue pour l'Us 2044 est discordante et nettement trop jeune pour pouvoir être retenue (5460 ±40 BP, soit 4368-4236 cal BC). D'un point de vue archéologique il est probable que cet échantillon recueilli dans l'ensemble 5 provienne d'une contamination liée aux creusements réalisés dans la partie supérieure de la séquence (fosses de l'ensemble 2). Si l'on conserve cette datation dans le traitement bayésien, elle tend à rajeunir la datation de cette phase (début 4547-4430 BC, fin 4483-4385 BC). Le traitement des trois dates se révèle beaucoup plus cohérent et situerait le début de la phase vers 4564-4417 cal BC et sa fin entre 4477 et 4367 BC avec un intervalle de crédibilité de 95 %.

Les faits archéologiques datés peuvent être attribués à une période comprise entre la seconde moitié du 46e siècle avant notre ère et le tout début du 44e siècle BCE. La phase 2 est bien postérieure à la phase à laquelle elle succède, bien que les plages de probabilités se chevauchent entre la fin de la phase 1 (4595-4458 BC) et le début de la phase 2 (4547-4430).

**Duration (posterior distrib.) :**

MAP : 44 Mean : 65 Std deviation : 37
Q1 : 37 Q2 (Median) : 60 Q3 : 86
HPD Region (95,0 %) : [2; 136] (95,0 %)
Credibility Interval (95,0 %) : [3, 133]

**Begin (posterior distrib.) :**

MAP : -4479 Mean : -4486 Std deviation : 38
Q1 : -4508 Q2 (Median) : -4483 Q3 : -4460
HPD Region (95,0 %) : [-4562; -4413] (95,0 %) BC/AD
Credibility Interval (95,0 %) : [-4564, -4417]

**End (posterior distrib.) :**

MAP : -4418 Mean : -4421 Std deviation : 29
Q1 : -4440 Q2 (Median) : -4420 Q3 : -4401
HPD Region (95,0 %) : [-4477; -4365] (95,0 %) BC/AD
Credibility Interval (95,0 %) : [-4477, -4367]

## - La phase 3 : les ensembles 4 et 7 de la zone 2

Cette phase n'est représentée que par deux datations obtenues dans la fouille de la zone 2. La première date retenue pour l'Us 2038 correspond à un probable remblai ou, plus vraisemblablement, à une zone de rejets domestiques (5530 ±40 BP, 4455-4330 cal BC). La seconde date a été réalisée pour l'Us 2045 (5610 ±40 BP, 4520-4356 cal BC) qui forme un épannage de l'ensemble 7, antérieur à la formation de l'Us 2038. L'intérêt de ces deux dates est de caractériser, d'un point de vue chronologique, un changement dans le statut de l'occupation de cette zone. Le début de la phase peut être estimé à l'intervalle de temps compris entre 4451 et 4343 cal BC ; la fin de la phase est datée entre 4428 et 4316 cal BC.

**Duration (posterior distrib.) :**

MAP : 8 Mean : 31 Std deviation : 25
Q1 : 11 Q2 (Median) : 24 Q3 : 43
HPD Region (95,0 %) : [0; 80] (94,9 %)
Credibility Interval (95,0 %) : [0, 79]

**Begin (posterior distrib.) :**

MAP : -4390 Mean : -4396 Std deviation : 29
Q1 : -4414 Q2 (Median) : -4395 Q3 : -4376
HPD Region (95,0 %) : [-4453; -4342] (95,0 %) BC/AD
Credibility Interval (95,0 %) : [-4451, -4343]

**End (posterior distrib.) :**

MAP : -4362 Mean : -4366 Std deviation : 28
Q1 : -4383 Q2 (Median) : -4365 Q3 : -4348
HPD Region (95,0 %) : [-4426; -4312] (95,0 %) BC/AD
Credibility Interval (95,0 %) : [-4428, -4316]

## - La phase 4 : incendie maison de la zone 3

Cette phase correspond à une maison incendiée en cours de fouille de la zone 3. Pour l'heure, nous ne disposons que d'une seule datation caractérisant le sommet de la couche d'incendie (5530 ±35 BP, 4450-4335). Il ne s'agit pas d'une phase *stricto sensu*, mais il nous a semblé intéressant d'observer la manière dont cet événement contribuait au modèle chronologique (contribution au modèle avec intervalle de confiance de 95% 4512- 4255 BC).
<figure><img src="image_100.png" /><figcaption>**Figure 5**

Traitement bayésien des dates de Taraschina.</figcaption></figure>
**Begin (posterior distrib.) :**

MAP: -4362 Mean: -4379 Std deviation: 67
Q1: -4421 Q2 (Median): -4379 Q3: -4343
HPD Region (95,0%): [-4511; -4249] (95,0%) BC/AD
Credibility Interval (95,0%): [-4512, -4255]

**End (posterior distrib.) :**

MAP: -4362 Mean: -4379 Std deviation: 67
Q1: -4421 Q2 (Median): -4379 Q3: -4343
HPD Region (95,0%): [-4511; -4249] (95,0%) BC/AD
Credibility Interval (95,0%): [-4512, -4255]

- La phase 5 – ensemble 2, fosses et niveaux de la zone 1

Cette phase réunit quatre dates obtenues pour des structures (Fs 1107 et 1134) et des niveaux de sols (Us 1102 et 1206) fouillés dans la zone 1. Ces faits archéologiques se rapportent à l'ensemble stratigraphique 2 (ENS 2) communs aux zones de fouille 1 et 2. Le début de la phase est estimé à l'intervalle 4522-4366 BC ; la fin de la phase est comprise entre 4428 et 4181 cal BC.

**Duration (posterior distrib.) :**

MAP: 84 Mean: 117 Std deviation: 70
Q1: 70 Q2 (Median): 102 Q3: 145
HPD Region (95,0%): [8; 259] (95,0%)
Credibility Interval (95,0%): [14, 256]

**Begin (posterior distrib.) :**

MAP: -4441 Mean: -4442 Std deviation: 40
Q1: -4464 Q2 (Median): -4440 Q3: -4417
HPD Region (95,0%): [-4522; -4363] (95,0%) BC/AD
Credibility Interval (95,0%): [-4522, -4366]

**End (posterior distrib.) :**

MAP: -4356 Mean: -4326 Std deviation: 67
Q1: -4368 Q2 (Median): -4342 Q3: -4302
HPD Region (95,0%): [-4433; -4179] (95,0%) BC/AD
Credibility Interval (95,0%): [-4428, -4181]

**Observation :**

on remarque que la dernière occupation mise au jour dans la zone 1 est pour partie contemporaine des niveaux d'abandon du bâtiment incendié de la zone 3. De la même manière, on observe une probable antériorité de cette occupation avec la phase d'abandon du site (phase 6).

- La phase 6 : ensemble 2 et fosses de la zone 2

La phase 6 réunit des datations obtenues pour des fosses-silos fouillées dans la zone 2. Ces structures ont été creusées dans les derniers niveaux du tell ; elles matérialisent la dernière occupation chalcolithique du tell dans cette zone. Le début de la phase peut être estimé à l'intervalle 4384-4267 BC ; la fin de la phase entre 4342 et 4108 BC.

Observation : La plage de temps probable couverte par cette phase montre une contemporanéité forte avec l'occupation du site 2. La date obtenue pour la structure Fs 2010 contribue à placer la date de l'abandon de cette partie du tell entre 4365-4256 (HPD à 95%), soit durant la première moitié du 43e siècle avant notre ère.

**Duration (posterior distrib.) :**

MAP: 57 Mean: 87 Std deviation: 59
Q1: 46 Q2 (Median): 74 Q3: 111
HPD Region (95,0%): [0; 211] (95,0%)
Credibility Interval (95,0%): [2, 211]

**Begin (posterior distrib.) :**

MAP: -4339 Mean: -4332 Std deviation: 30
Q1: -4351 Q2 (Median): -4334 Q3: -4314
HPD Region (95,0%): [-4387; -4268] (95,0%) BC/AD
Credibility Interval (95,0%): [-4384, -4267]

**End (posterior distrib.) :**

MAP: -4268 Mean: -4245 Std deviation: 61
Q1: -4283 Q2 (Median): -4257 Q3: -4221
HPD Region (95,0%): [-4346; -4107] (95,0%) BC/AD
Credibility Interval (95,0%): [-4342, -4108]

- La phase 7 : chronologie des occupations du site 2

La phase 7 correspond à l'occupation du site 2. Seules deux dates, issues de carottages, montrent que la phase a pu débuter dans l'intervalle 4482-4226 et s'achever entre 4354 et 4059 BC. La plage de temps couverte est pour partie contemporaine de la fin de l'occupation du site 1, et plus particulièrement de certaines des fosses creusées au sommet de la séquence stratigraphique.

**Duration (posterior distrib.) :**

MAP: 45 Mean: 117 Std deviation: 92
Q1: 45 Q2 (Median): 91 Q3: 157
HPD Region (95,0%): [0v; 292] (94,9%)
Credibility Interval (95,0%): [0, 288]

**Begin (posterior distrib.) :**

MAP: -4337 Mean: -4342 Std deviation: 74
Q1: -4365 Q2 (Median): -4334 Q3: -4300
HPD Region (95,0%): [-4482; -4219] (95,0%) BC/AD
Credibility Interval (95,0%): [-4482, -4226]

**End (posterior distrib.) :**

MAP: -4259 Mean: -4229 Std deviation: 77
Q1: -4281 Q2 (Median): -4242 Q3: -4186
HPD Region (95,0%): [-4356; -4057] (95,0%) BC/AD
Credibility Interval (95,0%): [-4354, -4059]

## 5. Conclusion

Le traitement des datations radiocarbone nous a permis de phaser et de dater les étapes successives de l'édification du tell. Ce travail a pris en compte les ensembles
archéologiques mis au jour dans les différentes zones du site et permet de proposer leur datation. Il a égale-ment permis d'établir l'existence probable d'un hia-tus ou d'un abandon du site vers 4350 BC, marquant un changement de statut et le déplacement total ou

partiel de l'habitat vers le site 2 (le site 1 étant occupé de manière secondaire). Le traitement bayésien montre bien la manière dont chacune des phases de succède. En revanche il s'avère peu efficient pour déterminer la chronologie des moments de ruptures.

<figure><img src="image_101.png" /><figcaption>**Figure 6**

Taraschina : modèle chronologique bayésien par phases.</figcaption></figure>
<table><thead><tr><th>SITE</th><th>ensemble daté</th><th>matériel</th><th>date BP</th><th>cal 68,2%</th><th>cal 95,4%</th><th>n° laboratoire</th></tr></thead><tbody><tr><td>TARASCHINA site 1</td><td>Fs 2010 - Us 1116</td><td>charbon</td><td>3080 ± 35 BP</td><td>1406-1296</td><td>1427-1260</td><td>Poz-42871</td></tr><tr><td>TARASCHINA site 1</td><td>Us 2030 N54</td><td>ossement</td><td>3320 ± 30 BP</td><td>1634-1532</td><td>1684-1523</td><td>Poz-56751</td></tr><tr><td>TARASCHINA site 1</td><td>Us 2005 N55</td><td>ossement</td><td>3835 ± 30 BP</td><td>2341-2207</td><td>2457-2200</td><td>Poz-76945</td></tr><tr><td>TARASCHINA site 1</td><td>Us 2042 M53ab</td><td>ossement</td><td>4045 ± 35 BP</td><td>2621-2491</td><td>2836-2473</td><td>Poz-76944</td></tr><tr><td>TARASCHINA site 1</td><td>Us 2005 N57</td><td>ossement</td><td>4135 ± 35 BP</td><td>2863-2630</td><td>2873-2578</td><td>Poz-56752</td></tr><tr><td>TARASCHINA site 1</td><td>Us 2005 N57 BIS</td><td>ossement</td><td>4190 ± 35 BP</td><td>2884-2700</td><td>2892-2640</td><td>Poz-57710</td></tr><tr><td>TARASCHINA site 1</td><td>Silo 1134</td><td>micro-faune</td><td>1650 ± 30 BP</td><td>351-424 AD</td><td>264-533 AD</td><td>Beta-304710</td></tr><tr><td>TARASCHINA site 1</td><td>Fs 2010 - Us 2009</td><td>ossement</td><td>5350 ± 35 BP</td><td>4313-4072</td><td>4323-4052</td><td>Poz-49883</td></tr><tr><td>TARASCHINA site 1</td><td>Log 19 (-1,45)</td><td>charbon</td><td>5390 ± 30 BP</td><td>4325-4236</td><td>4336-4076</td><td>Beta-350742</td></tr><tr><td>TARASCHINA site 1</td><td>Us 2044 N54a</td><td>ossement</td><td>5460 ± 40 BP</td><td>4350-4264</td><td>4368-4236</td><td>Poz-76943</td></tr><tr><td>TARASCHINA site 1</td><td>Log 33 (-218)</td><td>ossement</td><td>5470 ± 35 BP</td><td>4354-4267</td><td>4368-4246</td><td>Beta-76942</td></tr><tr><td>TARASCHINA site 1</td><td>Fs 2015 - Us 2016</td><td>ossement</td><td>5480 ± 35 BP</td><td>4358-4269</td><td>4444-4275</td><td>Poz-49884</td></tr><tr><td>TARASCHINA site 1</td><td>Silo 1134</td><td>ossement</td><td>5510 ± 30 BP</td><td>4369-4427</td><td>4282-4415</td><td>Beta-309452</td></tr><tr><td>TARASCHINA site 1</td><td>Us 2038</td><td>ossement</td><td>5530 ± 40 BP</td><td>4445-4340</td><td>4445-4330</td><td>Poz-65695</td></tr><tr><td>TARASCHINA site 1</td><td>Us 1206</td><td>ossement</td><td>5540 ± 40 BP</td><td>4446-4345</td><td>4456-4335</td><td>Poz-50175</td></tr><tr><td>TARASCHINA site 1</td><td>Us 1102</td><td>ossement</td><td>5570 ± 35 BP</td><td>4448-4363</td><td>4461-4344</td><td>Poz-46585</td></tr><tr><td>TARASCHINA site 1</td><td>Fosse 1107</td><td>ossement</td><td>5580 ± 35 BP</td><td>4449-4369</td><td>4487-4348</td><td>Poz-36217</td></tr><tr><td>TARASCHINA site 1</td><td>Us 2048</td><td>ossement</td><td>5580 ± 40 BP</td><td>4449-4369</td><td>4488-4347</td><td>Poz-65697</td></tr><tr><td>TARASCHINA site 1</td><td>Fs 2054 Us 2060</td><td>ossement</td><td>5600 ± 40 BP</td><td>4459-4369</td><td>4505-4351</td><td>Poz-65694</td></tr><tr><td>TARASCHINA site 1</td><td>Us 2045</td><td>ossement</td><td>5610 ± 40 BP</td><td>4486-4370</td><td>4520-4356</td><td>Poz-65696</td></tr><tr><td>TARASCHINA site 1</td><td>Us 2047</td><td>ossement</td><td>5610 ± 40 BP</td><td>4486-4370</td><td>4520-4356</td><td>Poz-65968</td></tr><tr><td>TARASCHINA site 1</td><td>Log 12 (-2,56)</td><td>charbon</td><td>5680 ± 30 BP</td><td>4541-4466</td><td>4590-4453</td><td>Beta-350741</td></tr><tr><td>TARASCHINA site 1</td><td>Log 34 -383</td><td>charbon</td><td>5690 ± 30 BP</td><td>4548-4466</td><td>4604-4456</td><td>Beta-442135</td></tr><tr><td>TARASCHINA site 1</td><td>Log 4 (base)</td><td>carporeste</td><td>5840 ± 340 BP</td><td>4780-4621</td><td>4798-4584</td><td>Poz-30313</td></tr><tr><td>TARASCHINA site 1</td><td>Us 3004-1 N34</td><td>ossement</td><td>5530 ± 35 BP</td><td>4444-4340</td><td>4450-4335</td><td>Poz-85830</td></tr></tbody></table>

Figure 7

Ensemble des dates réalisées dans les zones 1 et 2 du site de Taraschina.