1. Introduction

1L’étude de l’impact des usages anciens des sols sur la végétation forestière actuelle constitue depuis plusieurs années un axe de recherche important en écologie du paysage (Peterken et Game 1984, Hermy 1994, Honnay et al. 1998, Honnay et al. 2004). De nombreux travaux montrent que la composition floristique de la végétation forestière est influencée par la dynamique récente de l’occupation du sol. Les espèces végétales caractérisées par une faible capacité de colonisation apparaissent inféodées aux espaces continuellement occupés par la forêt depuis la fin de l’époque moderne. Les espèces plus compétitives, notamment les rudérales, sont en revanche plus fréquentes dans les forêts ayant plus ou moins récemment succédé à des terroirs anciennement cultivés dont les sols ont été amendés (Dupouey et al. 2002b).

2Les études développées dans ce domaine se fondent pour la plupart sur l’analyse régressive de documents cartographiques et restent par conséquent tributaires d’un seuil chronologique relativement tardif imposé par la documentation disponible. L’intégration de données archéologiques relatives aux forêts actuelles constitue une approche complémentaire susceptible de contribuer à l’étude de la dynamique de l’occupation du sol sur la longue durée. Depuis le début des années 1980, la question des relations entre structures archéologiques et végétation forestière a fait l’objet de nombreux travaux, notamment en France. Les recherches conduites dans ce domaine, qui ont particulièrement été développées en région Limousin (Ghestem 2002, 2006) et en Picardie (Vieille 2001, Decocq et al. 2002, Renaux 2003), ont privilégié l’étude d’occupations gallo-romaines et de sites castraux médiévaux. Ces approches se sont pour la plupart focalisées sur l’identification d’espèces végétales potentiellement indicatrices de structures archéologiques. La démarche se fonde sur l’étude de zones présentant une relative homogénéité écologique sur plusieurs dizaines d’hectares, privilégiant ainsi une approche à moyenne échelle. L’acquisition des données floristiques repose sur la réalisation de relevés phytosociologiques situés, d’une part, sur des secteurs comportant des structures archéologiques visibles et, d’autre part, sur des zones dépourvues de vestiges apparents utilisées comme témoins. La mise en œuvre de la méthode sur plusieurs sites et la constitution d’un ensemble de relevés appariés permet l’analyse comparative des relevés « in situ » et des relevés « hors-site ».

3Ces travaux, qui ont principalement été menés dans des contextes plus ou moins acidiphiles, montrent que les espèces indicatrices de sols relativement riches en bases et en éléments nutritifs sont significativement plus fréquentes dans les relevés effectués « in situ » que sur les relevés réalisés « hors site ». Les tendances s’inversent concernant la répartition des espèces acidiphiles, caractéristiques de la végétation potentielle. Ces différences de composition floristique reflètent des disparités sur le plan des propriétés physico-chimiques des horizons superficiels du sol introduites par l’occupation des sites. L’apport de matériaux exogènes, notamment de carbonates de calcium, provoque en particulier une augmentation du pH du sol et favorise le développement d’espèces caractéristiques d’horizons saturés en bases.

4Les approches phyto-écologiques développées sur sites archéologiques se sont pour la plupart inscrites dans une perspective visant à élaborer une méthode de prospection spécifique au milieu forestier. Elles ont dans le même temps négligé l’analyse, à travers une approche fine menée à grande échelle, de la structuration de la végétation sur des sites fossilisés par la forêt. Des travaux ont récemment été conduits dans cette perspective en Lorraine sur des occupations antiques associées à des structures parcellaires (Dupouey et al. 2002a). Ces études, fondées sur l’analyse d’un réseau de plusieurs dizaines de placettes, révèlent que la mise en valeur agricole, au cours de l’Antiquité, de terroirs postérieurement reconquis par la forêt peut entraîner une modification durable des propriétés physiques et chimiques des sols et que les effets rémanents de pratiques agricoles antiques conditionnent la structuration des communautés végétales sur les sites étudiés. Les zones anciennement cultivées ayant bénéficié des transferts de fertilité effectués du saltus vers l’ager se distinguent des secteurs mis en valeur de façon plus extensive par des sols plus profonds, moins acides et plus riches en éléments nutritifs, notamment en azote et en phosphore, ainsi que par une plus grande fréquence d’espèces végétales exigeantes sur le plan trophique comme les nitrophiles.

5Ces travaux montrent que l’utilisation des outils de bio-indication est susceptible de contribuer à l’analyse fonctionnelle des sites archéologiques actuellement localisés en forêt. Les études que nous développons actuellement sur une série d’occupations antiques du massif forestier de Rambouillet s’inscrivent dans cette problématique. L’originalité de la démarche réside dans l’acquisition de données écologiques (données floristiques, pédochimiques et pédobiologiques) selon un maillage systématique des zones étudiées, permettant ainsi la mise en œuvre d’analyses spatiales conjointement à la réalisation d’analyses statistiques multivariées.

2. Site et méthodes

6Le massif de Rambouillet (20000 ha), localisé dans le sud du département des Yvelines, recouvre une zone de plateaux dont le soubassement formé par un important dépôt de sables stampiens (Sables de Fontainebleau) est partiellement surmonté par des argiles à meulières. Dans cette région, où la couverture limoneuse est très irrégulière, la végétation forestière recouvre des secteurs où les dépôts loessiques sont sporadiques et peu épais (Crahet 1981). Les argiles à meulières et les sables stampiens y déterminent des sols aux potentialités agricoles limitées. Le substrat argileux donne des sols hydromorphes plus ou moins acides et localement caractérisés par une charge pierreuse importante. À partir des Sables de Fontainebleau se sont développés des sols plus ou moins podzolisés caractérisés par une acidité importante et par de faibles capacités de rétention en eau. Étant donné ces conditions édaphiques, le massif forestier de Rambouillet regroupe des formations végétales majoritairement acidiphiles. À l’exception des stations au sol inondé ou engorgé de façon permanente, la végétation se rattache à différents types de chênaies (Bournérias 1972, Bailly et al. 2001). La chênaie à Molinie (Molinia caerulea) s’étend sur les secteurs de plateaux mal drainés et dont les sols sont développés à partir des formations superficielles reposant sur les argiles à meulières. Les secteurs caractérisés par des couvertures limono-sableuses plus épaisses sont occupés par des formations végétales reflétant des niveaux trophiques plus élevés et de meilleures conditions de drainage (chênaie - charmaie à Ronce et Chèvrefeuille des Bois, chênaie - charmaie à Jacinthe des Bois). Sur les sols issus des sables stampiens, la végétation est caractéristique de sols drainants, acides et pauvres en éléments minéraux. Il s’agit principalement de la chênaie sessiliflore à Callune et Leucobryum glauque, d’une part, et de la chênaie sessiliflore à Fougère aigle, d’autre part. Il est à noter, par ailleurs, que les stations acidiphiles associées aux Sables de Fontainebleau présentent des faciès de dégradation ayant régressé au stade de boulaies à Callune (Bailly et al. 2001).

2.1. Le site de la Mare aux Buis

2.1.1. Contexte topographique, pédologique et floristique

7L’établissement antique de la Mare aux Buis est localisé sur la commune de Sonchamp, à environ 2,5 km au Sud-Est du centre-ville de Rambouillet. Situé dans la partie orientale du massif forestier de Rambouillet, il s’étend de part et d’autre de l’allée forestière marquant la limite entre les parcelles XXIII-11 et XXIII-22 de la forêt domaniale (Fig. 1). Localisé dans un vallon très ouvert, sur la rive droite d’un ruisseau, le site est établi sur un léger versant exposé au Sud-Ouest. Ce secteur géographique est caractérisé par des sols hydromorphes plus ou moins acides développés à partir d’une couverture limoneuse superficielle reposant sur des argiles à meulières. Sur l’emprise présumée du site, les formations végétales se répartissent en deux principaux types de stations : chênaie sessiliflore – boulaie acidiphile à Molinie, d’une part, et chênaie pédonculée – charmaie, d’autre part.

2.1.2. Historique des recherches

8Étudié à partir de la fin des années 1960 par les érudits locaux (Zuber 1969, Christmann 1970), le site de la Mare aux Buis a fait l’objet de plusieurs campagnes de prospection depuis près d’une quarantaine d’années. Sur l’emprise du site, les prospecteurs ont signalé de nombreux microreliefs se rapportant à la présence sous-jacente de fondations bien conservées, notamment de maçonneries associant petit appareil, principalement en meulière, et mortier à base de chaux. Au début des années 1990, une campagne de relevés topographiques a été effectuée sur le site à l’initiative du Service Archéologique Départemental des Yvelines (SADY), à la suite d’une coupe forestière réalisée sur la parcelle XXIII-11. Ces relevés ont récemment été complétés par des recherches que nous avons entreprises en 2005 et 2006 en parcelle XXIII-22. Une campagne de relevés microtopographiques a d’une part été effectuée sur une surface de près de 3 ha, ce qui a permis l’élaboration d’un modèle numérique de terrain à haute résolution spatiale (Fig. 1). Par ailleurs, simultanément aux études effectuées sur les sols et la végétation forestière, il a été procédé au ramassage systématique du mobilier archéologique sur les placettes utilisées pour l’acquisition des données écologiques (voir paragraphe 2.2.).

2.1.3. Organisation spatiale et chronologie du site

9Les travaux successivement conduits sur le site permettent, d’une part, de préciser son organisation générale. L’établissement antique présente les vestiges d’un enclos orienté Nord-Nord-Ouest – Sud-Sud-Est délimitant un espace mesurant au minimum 190 m sur 150 m (soit une superficie minimale de 2,28 ha). Conservé sur trois côtés et vraisemblablement pourvu d’une d’entrée sur son côté Ouest, il est matérialisé par un talus au sommet duquel des blocs de meulières régulièrement équarris suggèrent la présence sous-jacente d’un mur d’enceinte. Les relevés topographiques réalisés par le SADY ont permis la restitution, à l’intérieur de l’enclos, du plan d’un bâtiment organisé autour d’une cour quadrangulaire d’une superficie d’environ 2700 m². Contrairement à la parcelle XXIII-11 qui regroupe l’essentiel de la partie résidentielle de l’établissement, la parcelle XXIII-22 est caractérisée par une faible densité de vestiges maçonnés, à l’exception d’un secteur délimité par un petit enclos fossoyé d’environ 20 m de côté, à l’intérieur duquel des matériaux de construction sont observés à la surface du sol (tuiles, appareil en meulière et calcaire, mortier à base de chaux). Une mare de plan rectangulaire, d’environ 40 m sur 20 m, est adossée au talus Ouest de l’enclos, ces deux structures se situant dans le prolongement de l’aile Sud du bâtiment sur cour.

Fig. 1. Plan des structures archéologiques et du dispositif d’étude

10Le site a livré un lot de 271 tessons se rapportant à une période chronologique comprise entre le milieu du I^er s. av. J.-C. et la fin du III^e s. ap. J.-C. La majorité de la céramique collectée, qui se répartit sur 21 des 148 placettes du réseau, provient de l’emprise du petit enclos fossoyé et des secteurs situés à l’intérieur de l’enceinte de l’établissement. Ces zones, caractérisées par des humus actifs sur le plan biologique (mull eutrophe et mull mésotrophe), présentent dès le début du printemps une litière discontinue en raison du turnover rapide de la matière organique. Dans ce contexte, la probabilité de retrouver du matériel épandu à la surface du sol est particulièrement élevée, ce qui peut fausser le résultat des ramassages de surface en accentuant les disparités entre zones de fortes concentrations de matériel et zones de vide apparent. Néanmoins, le mode de répartition de la céramique suggère des contrastes importants concernant l’occupation du site et peut rendre compte d’une pression anthropique plus forte à proximité des structures bâties. Le degré de conservation des structures archéologiques suggère que le site n’a pas fait l’objet de labours agricoles depuis son abandon. À l’appui de cette hypothèse, des datations radiocarbone provenant d’échantillons de charbons de bois prélevés sur 7 charbonnières situées en parcelle XXIII-22 révèlent que les activités de charbonnage ont été pratiquées de manière récurrente depuis le Haut Moyen Âge et suggèrent la continuité de la forêt depuis cette période sur la zone étudiée (Tab. 1).

Tab. 1. Datations radiocarbone obtenues sur 7 charbonnières

2.2. L’acquisition des données écologiques : matériel et méthodes

11Les approches écologiques développées sur le site ont exclusivement porté sur la parcelle XXIII- 22, étant donné que la parcelle XXIII-11, qui a fait l’objet de plusieurs coupes forestières successives depuis une quinzaine d’années, présente un faciès au stade de fourrés arbustifs et demeure de ce fait difficilement accessible. L’aire d’étude, mesurant 290 m sur 190 m (soit une superficie de 5,51 ha), comporte deux types de biotopes. Les secteurs comportant les structures archéologiques sont recouverts par une formation végétale de type chênaie - charmaie traitée en taillis sous futaie. La strate arborescente, très recouvrante, favorise les espèces de sous-bois ombragés au niveau de la strate herbacée. Les périphéries Ouest et Nord de la zone étudiée sont recouvertes par une formation de type chenaie - boulaie à Molinie. Le milieu est plus ouvert et favorise les espèces végétales héliophiles (Molinie, Ronce et Chèvrefeuille).

12Un dispositif de 148 placettes de 100 m² a été mis en place selon un maillage systématique de l’aire étudiée basé sur une résolution spatiale de 20 m (Fig. 2). Un relevé phytosociologique a été réalisé sur 133 placettes, les placettes situées à proximité immédiate de la route forestière traversant le site ayant été exclues de l’étude afin de limiter les effets de pollutions éventuelles1. Les relevés, effectués par strate de végétation (arborescente, arbustive et herbacée), se sont limités à la flore vasculaire. Pour chacune des strates, l’abondance relative de chaque espèce a été exprimée au moyen d’un coefficient d’abondance - dominance, selon la méthode sigmatiste (Braun-Blanquet 1951)2. Sur chacune des 148 placettes, trois carottages manuels ont été effectués au niveau de l’horizon hémi-organique du sol (entre 1,5 et 10 cm de profondeur à partir de la surface du sol, selon le type d’humus en présence) de manière à constituer un échantillon composite représentatif de chaque unité étudiée. Après séchage à l’air libre, les échantillons ont été transmis au laboratoire d’analyses des sols de l’INRA où ils ont été soumis à des analyses chimiques (pH, rapport carbone total/azote total – C/N, phosphore assimilable – P2O5 et cations échangeables : aluminium – Al³⁺, calcium – Ca²⁺, potassium – K⁺ et magnésium – Mg²⁺). Des prélèvements de sol et de litière ont par ailleurs été effectués sur un ensemble de 57 placettes pour procéder à des études de pédofaune3. L’humidité relative du sol a été déterminée pour chaque échantillon après séchage complet des prélèvements. Les données concernant la faune du sol, en cours d’analyse, ne sont pas présentées dans la présente contribution.

Fig. 2. Relevés microtopographiques

13Après transformation des coefficients d’abondance – dominance en valeurs numériques selon la méthode de pondération proposée par van der Maarel (1979)4, les données collectées lors des relevés phytosociologiques ont été étudiées au moyen d’une analyse factorielle des correspondances (AFC), l’analyse ayant été restreinte à la strate herbacée. Pour contribuer à l’interprétation des gradients écologiques exprimés par l’AFC, les variables pédochimiques ont été intégrées à l’analyse en tant que critères additionnels et regroupées au sein d’une matrice passive de présence/absence. Pour ce faire, les valeurs de chacune des variables ont été reparties en trois classes au moyen d’un partitionnement univarié selon la méthode de Fisher (1958)5, les classes obtenues, notées par le suffixe « - », « = » et « + », constituant ainsi de nouvelles variables (Tab. 2a) dont les valeurs ont été codées selon une modalité binaire (0/1). Une série de variables topographiques a été intégrée à l’analyse selon le même principe afin d’apprécier les relations entre le positionnement des placettes par rapport aux structures archéologiques et leur composition floristique (Tab. 2b). La distance entre chaque placette et différentes entités spatiales (enclos principal, enclos fossoyé et mare adjacente, fondations) a été calculée par la mise en œuvre de requêtes effectuées au moyen d’un système d’information géographique. La transformation de ces variables quantitatives en variables binaires s’est fondée sur la même méthode que celle adoptée pour le traitement des variables pédochimiques. Par ailleurs, d’autres critères, initialement codés de façon binaire (présence de matériaux de construction, présence de vestiges de fondations, présence de céramique, placette située sur l’emprise de l’enclos fossoyé et/ou de la mare adjacente, placette située sur l’emprise de l’enclos principal) ont été joints à la matrice de données, sans transformation préalable. Les variables édaphiques (variables pédochimiques et humidité relative du sol) des 57 placettes qui font actuellement l’objet d’études pédobiologiques ont été étudiées au moyen d’une analyse en composantes principales. Chacune des variables pédochimiques a d’autre part fait l’objet d’une analyse géostatistique à travers la réalisation d’un semi-variogramme afin de déterminer si la distribution de chaque paramètre sur l’aire étudiée est aléatoire ou si elle obéit à une structure spatiale. Dans ce dernier cas, les variables concernées – en l’occurrence, pH, P2O5, Al³⁺, Ca²⁺, K⁺ et Mg²⁺– ont fait l’objet d’une interpolation spatiale par kriegeage.

Tab. 2. a-b. Critères additionnels intégrés dans l’analyse factorielle des correspondances

3. Résultats et discussion

14Les espèces contribuant le plus au premier axe du plan factoriel de l’AFC (Fig. 3 et Fig. 4) se répartissent en deux ensembles écologiques bien distincts. On peut distinguer, d’une part, un pôle neutrophile, avec Acer campestre, Mercurialis perennis (espèces neutrocalciphiles), Ajuga reptans, Ranunculus ficaria (neutronitrophiles), Arum maculatum, Carex sylvatica (neutronitroclines), Dryopteris filix-mas, Lamiastrum galeobdolon, Millium effusum (neutroclines à moyenne amplitude) et, d’autre part, un pôle acidiphile avec Deschampsia flexuosa, Hypericum pulchrum et Molinia caerulea. L’analyse de la contri bution des espèces à l’axe F1 montre que la végétation se structure principalement selon un gradient trophique, ce que confirme la position des variables pédochimiques le long de l’axe (Fig. 3). Les placettes dont les scores factoriels sont les plus bas sont caractérisées par les sols les plus acides, les plus désaturés en azote et en minéraux. Il s’agit également des sols dont le complexe absorbant est dominé par l’ion aluminium (Al³⁺), lequel exerce une action toxique à l’égard du développement racinaire des espèces tolérant mal les sols acides. À l’autre extrémité du gradient, les placettes correspondent aux sols les plus proches de la neutralité et présentant les capacités nutritives les plus élevées, notamment pour le calcium et le phosphore. Le premier axe de l’AFC exprime également un gradient de luminosité, le pôle acide regroupant des espèces de pelouses et d’ourlets plus ou moins héliophiles, le pole neutre regroupant des espèces de forêts se développant préférentiellement au sein de sous-bois relativement ombragés.

Fig. 3. Plan F1 x F2 de l’analyse factorielle des correspondances (AFC) réalisée sur la composition floristique de 133 placettes (variables pédochimiques traitées en critères additionnels)

Fig. 4. Plan F1 x F2 de l’analyse factorielle des correspondances (AFC) réalisée sur la composition floristique de 133 placettes (variables géomatiques traitées en critères additionnels)

15Le deuxième axe du plan factoriel tend à opposer les espèces à amplitude écologique relativement restreinte aux espèces à large amplitude trophique et hydrique, comme Rubus fruticosus et Lonicera peryclimenum. Les placettes correspondant aux scores factoriels négatifs de l’axe F2 caractérisent des faciès forestiers de transition entre, d’une part, la chênaie – boulaie à Molinie et la chênaie – charmaie, d’autre part. Les cortèges floristiques témoignent de niveaux trophiques intermédiaires : la strate herbacée, largement investie par la Ronce et le Chèvrefeuille, est caractérisée par la présence, en proportions relativement équilibrées, d’espèces à large amplitude, acidiclines et neutroclines (Convallaria maialis, Anemone nemorosa, Euphorbia amygdaloides, Hedera helix).

16Les placettes contribuant le plus au premier axe du plan factoriel de l’AFC (notamment F2, G2, H2, G3, H3 et G4) présentent majoritairement des scores factoriels positifs et correspondent par conséquent aux relevés caractérisés par les niveaux trophiques les plus élevés. Ces placettes où les espèces acidiphiles sont très peu fréquentes hébergent, d’une part, des espèces forestières indicatrices de sols relativement riches en bases et en éléments nutritifs et, d’autre part, des espèces caractéristiques de milieux rudéralisés dont les sols offrent de bonnes réserves en eau (Ajuga reptans, Gallium aparine, Stachys sylvatica). Ces placettes coïncident avec les pics observés pour les taux de calcium échangeable et les teneurs en phosphore assimilable. Celles-ci se concentrent autour du petit fossoyé et à proximité de la mare adjacente. Les placettes ayant livré de la céramique et celles qui présentent des matériaux de construction associés à des vestiges de fondations sont très significativement plus fréquentes parmi les placettes correspondant aux scores factoriels les plus élevés. La distance des placettes vis-à-vis des structures archéologiques détermine par ailleurs leur distribution le long du gradient trophique exprimé par l’axe F1 du plan factoriel (Fig. 4). Les placettes correspondant aux classes de distances les plus réduites vis-à-vis des différentes structures sont très significativement plus fréquentes parmi les placettes ayant les scores factoriels les plus élevés. Inversement, les placettes les plus éloignées des structures sont plus fréquentes parmi les placettes affectées des scores les plus négatifs.

17L’analyse en composantes principales (ACP) réalisée sur les variables édaphiques conforte les tendances dégagées lors de l’analyse factorielle des correspondances (Fig. 5). La première composante principale est fortement corrélée avec les valeurs du pH, les taux de calcium et ceux de magnésium. Elle est d’autre part négativement corrélée avec les taux d’aluminium et les valeurs observées pour le rapport carbone/azote. L’opposition entre ces groupes de caractères montre que le gradient qui s’exprime sur le premier axe de l’ACP correspond à un gradient trophique. Il oppose, à gauche de l’axe, des placettes caractérisées par un sol acide, désaturé et peu actif sur le plan biologique et, à droite, des placettes dont le sol est moins acide, plus riche en éléments nutritifs et marqué par un turnover plus rapide de la matière organique. La seconde composante principale est fortement corrélée à l’humidité relative du sol et exprime un gradient hydrique, les placettes situées dans la partie inférieure du graphique se rapportant à des sols relativement secs, celles placées dans la partie supérieure à des sols moins bien drainés.

18Le plan factoriel montre que les sols caractérisés par les pH les plus élevés sont relativement secs. Les placettes ayant livré de la céramique antique et/ou présentant des vestiges de fondations correspondent pour la plupart à ces caractéristiques. La présence de structures bâties, ainsi que l’existence de niveaux de démolition, impliquant, au niveau des horizons pédologiques superficiels, une charge importante en éléments grossiers, notamment en matériaux carbonatés, fournit à cet égard une explication possible du caractère à la fois drainant et peu acide du sol. Les placettes situées sur l’emprise de l’enclos principal du site qui sont dépourvues de vestiges de fondation et de matériaux de construction tendent à se répartir en deux groupes très différents : le premier (placettes M2, M3, N1 et N2), situé à proximité de la limite Nord de l’enceinte, est caractérisé par des sols humides et acides ; le second (placettes J1, I2, H1), situé le long du côté Ouest de l’enclos, par des sols secs et modérément acides. Ce constat conduit à émettre l’hypothèse d’affectations différenciées au sein de l’enclos, sans qu’il soit possible, dans l’état actuel des données, de statuer sur les facteurs à l’origine de ces disparités.

Fig. 5. Plan F1 x F2 de l’analyse en composantes principales (ACP) de 8 variables édaphiques réalisée sur 57 placettes

Fig. 6. Cartographie par interpolation spatiale des taux de phosphore assimilable au niveau de l’horizon hémi-organique du sol (résolution de l’interpolation : 5 m)

19Les placettes présentant les valeurs les plus élevées pour les taux de phosphore assimilable correspondent à des sols à tendance hydromorphe qui caractérisent les abords Sud de la mare et de l’enclos fossoyé (placettes G3, E3, D2 et F2). Ces mêmes placettes contribuent par ailleurs fortement à l’inertie du premier axe du plan factoriel de l’AFC et se distinguent par une proportion relativement élevée d’espèces rudérales hygroclines à mésohygrophiles (Ajuga reptans, Galium aparine, Ranunculus ficaria). La cartographie des taux de phosphore assimilable au niveau de l’horizon hémi-organique du sol (Fig. 6) montre que ce secteur s’inscrit au sein d’une zone de forme subrectangulaire d’environ 800 m² caractérisée par des valeurs supérieures à 400 ppm et au sein de laquelle les taux atteignent un pic à 921 ppm, c’est-à-dire un niveau comparable à ceux que l’on peut observer sur des terres actuellement cultivées et soumises de longue date à des pratiques de fertilisation. L’effet de seuil traduit par la variabilité des taux de phosphore permet de poser l’hypothèse de fonctions spécifiques à proximité de l’enclos fossoyé et de la mare attenante. La disponibilité élevée en phosphore du sol peut indiquer des apports en fumure animale. Par ailleurs, le tassement du sol pourrait expliquer l’accentuation de son caractère hydromorphe. Ceci conduit à envisager l’hypothèse d’une aire de stabulation attenante à l’enclos principal de l’établissement antique et aménagée à proximité immédiate d’un point d’eau pour le bétail.

4. Conclusion

20Les recherches interdisciplinaires conduites sur l’établissement antique de la Mare aux Buis montrent que la structuration des communautés végétales sur la zone étudiée et la variabilité des propriétés chimiques du sol répondent en grande partie à des modifications anciennes du milieu introduites par l’occupation du site au cours du Haut-Empire. L’étude des gradients écologiques qui s’expriment sur le site, à travers l’analyse d’indicateurs phyto-écologiques et pédochimiques, permet de poser un certain nombre d’hypothèses concernant les modalités de l’occupation du site et de l’exploitation du sol. L’intégration de données écologiques selon une approche spatiale conduite à grande échelle fournit par ailleurs des éléments d’appréciation de l’emprise du site et de son organisation. Ce type de démarche ouvre de ce point de vue un certain nombre de perspectives méthodologiques concernant les modalités d’utilisation des outils de bio-indication dans le cadre de recherches archéologiques spécifiques aux espaces forestiers.