Introduction

1Une attention grandissante est portée aux écoulements hydro-sédimentaires de versant dans les hydrosystèmes du Nord-Ouest de la France en raison de leurs conséquences négatives : pertes en terre, pollution, diminution de la fertilité des sols, inondations, crues turbides, etc. Cependant, quantifier l’impact d’une haie et/ou d’un fossé sur le parcours de ces écoulements est bien délicat, sauf si l'on se trouve à proximité du réseau linéaire au moment même où le processus est en cours. En mesurant les débits entrants et sortants, l’observateur pourrait caractériser le fonctionnement hydrologique de ce linéaire (Ouvry et al., 2012) mais cette connaissance, aussi fine soit-elle, ne pourra pas être dupliquée simultanément sur tous les linéaires du bassin versant étudié. Connaitre l’impact de chaque linéaire sur la réponse hydrologique globale d’un bassin devient alors impossible (Baudry et Jouin, 2003). Si les conditions de genèse du ruissellement dans les parcelles agricoles sont aujourd’hui bien établies (Boiffin, 1984 ; Auzet, 1987 ; Boiffin et al., 1988 ; Delahaye, 1992 ; Souchère et al., 1998, Le Bissonnais et al., 2002), le parcours des écoulements de surface (dès que ces derniers quittent la parcelle) est bien plus difficile à appréhender, notamment en contexte bocager où la structure paysagère tend à complexifier ce phénomène (Mérot et al., 1998 ; Gascuel-Odoux et al., 2009 ; Viel et al., 2014). Plusieurs travaux ont déjà mis en évidence le rôle des réseaux linéaires d’origines anthropiques tels que les réseaux de haies, de fossés et de routes (Jones et al., 2000 ; Duke et al., 2003 ; Carluer et De Marsily, 2004 ; Viaud, 2004 ; Branger, 2007 ; Wohlfahrt et al., 2010 ; Viel, 2012 ; Viel et al., 2014). Mérot et al. (1998) montrent par exemple que les haies sur talus et perpendiculaires à la pente tendent à bloquer les flux de surface car ils favorisent l’infiltration en pied de talus. Carluer et De Marsilly (2004) et Moussa (2002) ont appréhendé les réseaux de fossés comme des vecteurs additionnels d’écoulement, tout comme Duke et al. (2003) l’ont fait pour les réseaux routiers. Le Gouée et Delahaye (2008), puis Viel et al (2014), proposent une typologie des espaces selon la connectivité au réseau hydrographique :

• espaces directement connectés grâce à la proximité à un cours d’eau ;

• espaces déconnectés par les haies sur talus ;

• ou espaces indirectement connectés via les réseaux viaires et de fossés.

2Ces travaux mettent en évidence la variabilité de l’efficacité hydrologique du bocage et soulignent le rôle fondamental de la densité de drainage (qu’il soit naturel ou anthropique) qui permet une connexion rapide du ruissellement au cours d’eau et des entrées de champs qui permettent le passage d’un écoulement de surface vers un écoulement de réseau. Ainsi, pour comprendre l'ensemble de ces dynamiques hydro-sédimentaires depuis les parcelles jusqu’au cours d’eau, il faut donc intégrer l’ensemble des entités composant la structure paysagère (Viel et al., 2014). Au-delà des connaissances empiriques acquises, les gestionnaires et/ou les chercheurs ont aujourd’hui besoin d’outils capables de reproduire l’ensemble de ces dynamiques hydro-sédimentaires et de mesurer le rôle de la structure paysagère (actuelle, passée ou même future) sur les transferts de versant (Baudry et Jouin, 2003 ; Gascuel-Odoux et al., 2009 ; Wohlfart et al., 2010). Plusieurs modèles intégrant les réseaux linéaires d’origine anthropique ont été développés ces dernières années : ANTHROPOG (Carluer et De Marsilly, 2004), STREAM (Souchère, 1998, Cerdan et al., 2001) ; RIDEM (Duke, 2003), SACADEAU (Gascuel et al., 2009) ; MHYDAS (Moussa, 2002). Certains sont d’ores et déjà opérationnels et tous ont prouvé leur efficacité dans leur contexte d’application (Moussa, 2002 ; Duke et al., 2003 ; Carluer et De Marsily, 2004 ; Tortrat, 2005 ; Bocher, 2005 ; Doray, 2008 ; Charlier, 2007 ; Levavasseur, 2012). Si ces modèles intègrent le rôle hydrologique des réseaux linéaires et proposent des simulations proches de la réalité (débits de matières à l’exutoire par exemple), ils restent peu adaptés pour mettre en évidence la dynamique des processus d’écoulement à travers les échelles spatiales (depuis les parcelles jusqu'au cours d'eau : groupes de parcelles, tête de bassin, etc.). Or c’est précisément cette dynamique spatiale, souvent fortement influencée par les réseaux linéaires, que l’on cherche à appréhender.

3Ainsi, face à ces différents constats, un outil de modélisation a été conçu et développé sous système multi-agents (SMA) pour recréer les dynamiques de ruissellement. Ce type de modélisation a été choisi pour sa capacité à appréhender un système complexe par l’implémentation de règles d’interactions simples et faciles à comprendre (Langlois, 2002 ; Caillault, 2011). Le développement de plateformes facilement accessibles telle que NetLogo (Wilensky, 1999) ou MAGéo (Blanpain et al., 2013) a rendu populaire l’utilisation des SMA en sciences humaines et sociales, et notamment en géographie environnementale (Caillault, 2011 ; Rebaudo et al., 2011 ; Vannier, 2011 ; Schwarz, 2012 ; Caillault et al., 2013 ; Anselme et al., 2014). Dans l’étude des processus physiques des écoulements de surface, l’approche par système multi-agents est encore assez récente mais trouve son intérêt dans sa capacité à offrir de l’information sur l’origine des dynamiques spatiales des écoulements et sur leur variabilité temporelle (Servat, 2000 ; Reaney, 2008, Rakotoarisoa, 2014). L'approche "agent" permet, en effet, d'obtenir des informations détaillées sur les dynamiques hydrologiques à travers le paysage au cours d'un épisode pluvieux : trajectoires d'écoulements, interaction infiltration/écoulement, destin de chaque particule d'eau ou encore, origine de l'eau alimentant l'hydrogramme à l'exutoire (Reaney, 2008). Elle permet également de mettre en place, assez facilement, des scénarios d'aménagements et d'évaluer leurs conséquences sur les écoulements. Dans ce domaine, seuls trois travaux peuvent être cités :

1. les travaux de David Servat et plus globalement du projet Rivage sur le développement d’une approche multi-agents pour simuler les interactions entre les acteurs hétérogènes de l’infiltration et du ruissellement d’eau sur une surface de sol ;

2. ceux de Sim Reaney sur le couplage d’un modèle physique avec des agents représentant des gouttes d’eau ;

3. ceux de Mahefa Rakotoarisoa sur le développement d’un système multi-agent couplé à une fonction de production de type SCS_Curve Number.

Développements méthodologiques

5Le modèle a été créé sous Netlogo (Wilensky, 1999, 2004), version 5.0.3. La présentation du modèle suit ici le protocole ODD (Overview, Design Concept and Details) mis en place afin de standardiser la description des modèles individus-centrés et multi-agents (Grimm et al., 2006, 2010).

Vue d’ensemble (Overview)

Objectif

6L'objectif du modèle est de mettre en évidence le rôle des entités linéaires sur les écoulements de surface dans les bassins versants agricoles fortement anthropisés. En plus de l’intérêt scientifique que présente le développement de ce système (pertinence des SMA pour mieux comprendre les dynamiques spatiales du ruissellement de surface), l’objectif est également d’ordre opérationnel avec la possibilité d’identifier les points stratégiques dans un bassin (haies limitant le ruissellement, fossés qui connectent au cours d’eau, espaces connectés au cours d’eau…) tout en suivant en tout point de l’espace la dynamique, tant spatiale que temporelle, des processus de ruissellement de surface.

7Pour cela, le modèle s'appuie sur une reproduction des écoulements de surface par une approche agent, c'est-à-dire par une discrétisation de l’élément dynamique, l’eau, en plusieurs agents (Treuil et al., 2008). Les agents se déplacent en fonction de la pente jusqu’à atteindre l’exutoire de l’environnement de simulation. Au cours de leur déplacement, les agents sont amenés à interagir avec différents éléments de la structure paysagère que sont les haies et les fossés. Il n’y a pas de prise en compte des vitesses d’écoulement : les agents se déplacent en fonction des itérations dites "machine". D’autres simplifications ont été opérées pour avoir un modèle le plus simple possible (pas de prise en compte de l’occupation du sol, de la pédologie, de la géologie ou des flux sédimentaires par exemple).

Entités

9Agentgoutte : Ce sont des entités individuelles qui se déplacent dans l’environnement de la simulation. Chaque "agentgoutte" contient un volume d’eau discret : autrement dit, les écoulements sont discrétisés en agents. Ces agents sont comparables aux "boules d’eau" (ou paquet d’eau) de Servat (2000) dans le sens où, à la différence des gouttes de pluie, leur existence se poursuit au-delà même du contact avec le sol. Les "agentgouttes" se déplacent dans un environnement qui est constitué de patches et avec lesquels ils interagissent. Le déplacement est conditionné par la variable "altitude" des patches (issue d’un modèle numérique de terrain). Les agents se déplacent itérativement de proche en proche en suivant la plus grande pente sauf s’ils se situent sur un patch identifié comme "fossé" ; dans ce cas, ils devront suivre la pente imposée par l’élément linéaire. Une fois arrivé à l’exutoire du linéaire l’"agentgoutte" se déplace sur le patch voisin de plus petite altitude. Si ce patch correspond à un fossé, alors l’agent reste dans un écoulement de réseau ; si ce n’est pas le cas, l’agent reprend un écoulement topographique. Chaque "agentgoutte" dispose de variables spécifiques mais qui peuvent ensuite évoluer au cours de la simulation :

Variables propres à chaque agent

start-patch : nom du patch de départ ;
stock_en_eau : quantité d’eau discrète (en mètre précipité, équivalent d’une lame d’eau précipitée) ;
distance_parcourue : distance totale parcourue au cours d’une simulation (en mètre) ;
distance_jusqu'à_un_lineaire : distance parcourue avant de rencontrer un patch de type linéaire (en mètre) ;
distance_parcourue_sur_lineaire : distance parcourue sur un linéaire (en mètre).

1. Les patches "haies" (couleur verte) : ils ont une valeur d’altitude systématiquement supérieure à celles de leurs voisins de types différents (fossé ou neutre) qui leur confère un rôle d’obstacle à l’écoulement (un "agentgoutte" ne se déplaçant que sur un patch voisin d’altitude inférieure) ;

2. • 5 Profondeur moyenne des fossés observées sur le terrain d'étude.

   Les patches "fossés" (couleur rouge) : ils ont toujours une altitude inférieure (environ 50 cm5) à celles de leurs voisins de types différents (haie ou neutre). Cette diminution de l’altitude des patches fossés favorise ainsi leur capacité à concentrer les "agentgouttes". Cette catégorie intègre également le réseau routier qui est considéré ici comme étant encaissé dans le versant ;

3. Les autres patches sont dits "neutres" (aucune occupation du sol spécifique ne leur est attribuée). Majoritaires dans l’environnement de simulation, ces entités peuvent correspondre à toutes les autres entités composant la structure paysagère d’un bassin versant agricole : parcelles agricoles, mares, entrées de champs, zones urbanisées, cours d’eau, etc. Leur intégration future dans le modèle est envisagée.

11Les patches possèdent, eux aussi, des variables qui leur sont propres. Certaines sont issues du travail préparatoire sous SIG (l’altitude par exemple) mais d’autres sont attribuées en début de simulation. Le terme de "patch linéaire" sera utilisé par la suite pour désigner les patches représentant une portion de réseau linéaire.

Variables propres à chaque patch

- altitude : valeur d’altitude issue d’un MNT (obtenue sous SIG)
- cible : nom du patch voisin de plus petite altitude (obtenue sous SIG, valeur égale à direction)
- linear : neutre (valeur numérique : 0), haie (10), fossé (20)
- flux : nombre d’"agentgouttes" passés sur le patch au cours de la simulation
- exutoire : vrai si le patch est identifié comme l’exutoire du bassin versant étudié
- proximite_d'un_lineaire : vrai si l’un des voisins est de type linéaire

+ autres variables envoyées par l’"agentgoutte" en fin de simulation (distance parcourue…)

Echelles spatiales et temporelles

L’espace

12L’environnement de simulation est discrétisé en patches de 10x10 mètres. La discrétisation de l’espace, effectuée sous SIG se base sur les mailles carrées régulières d’un modèle numérique de terrain obtenu auprès de l’Institut National de l’Information Géographique et Forestière (IGN). Les données composant l’espace d’étude sont présentées dans la partie données en entrées. La taille de l’environnement de simulation doit être plus grande que celle du bassin étudié (afin d’éviter tout problème lié aux effets de bordure) et de taille rectangulaire de façon à intégrer l’ensemble du bassin versant (Figure 1). Tous ces prétraitements (définition d'un environnement de travail et exportation des données) sont réalisés sous un système d’information géographique (SIG). La taille des cellules est fixée à 10 mètres, ce qui est un bon compromis entre précision et réalité du terrain (cette résolution présente aussi l’avantage de ne pas trop altérer la qualité des données d’entrée sans engendrer un trop grand nombre de cellules).

Le temps

13Les simulations se fondent sur les itérations de la plateforme NetLogo. Dans cette première version, le modèle est de type événementiel. Aucune fonction de transfert intégrant les vitesses de déplacement n’a été intégrée : à chaque itération, les "agentgouttes" se déplacent d’un patch à un autre situé dans un voisinage de Moore.

Diagramme d’activité

Eléments de conception (Design concept)

14Les grands concepts du modèle sont présentés ci-dessous selon l’ordre établi par Grimm et al. (2006).

Principes

15En comparant les simulations avec et sans l’intégration des réseaux, on cherche à mesurer l’impact des réseaux linéaires sur le déplacement des agents. Les premières simulations sont dites "topographiques" (i.e sans les réseaux) : les "agentgouttes" se déplacent uniquement en fonction des valeurs "altitude" du MNT à partir du sous-modèle Ecoul_Ciblé. Les contraintes à l’écoulement sont ensuite intégrées aux simulations (i.e avec les réseaux) et prises en compte par les "agentgouttes" dans leur déplacement (sous-modèles Ecoul_Topo et Ecoul_Réseau). On part donc du principe que les interactions entre les agents se limitent à deux fonctions : combler les cuvettes topographiques présentes dans l’environnement de simulation et adapter leur déplacement selon les contraintes linéaires. Chaque règle décisionnelle est appliquée aux "agentgouttes" c'est-à-dire à un niveau micro (ou local). De l’ensemble de ces interactions émerge la dynamique globale d’écoulement, et il est alors possible d’évaluer le poids des contraintes paysagères sur cette dynamique à l’aide de différents indices de sensibilité (détaillés ci-après).

Emergence

16Les simulations mettent en évidence les axes préférentiels d’écoulement, les changements induits par les réseaux, les espaces où les "agentgouttes" interagissent pour combler les cuvettes et les endroits où les agents quittent leur circulation de réseau pour un écoulement "topographique". Dans cette première approche, l'émergence produite au cours des simulations reste limitée. Mais l'intégration prochaine aux patches d'attributs "occupation des sols" et d’autres liés à des données "physiques" (pédologie, épaisseur de sol, antécédents hydro-climatiques, etc.) renforcera les interactions "agentgouttes" / patches, desquelles émergeront davantage de phénomènes non prévisibles (zones de ruissellement réel, prairies traversées par un écoulement, fossés saturés, etc.).

Adaptation

17Lorsqu’un "agentgoutte" se trouve sur un patch qui n’a pas d’exutoire (comme par exemple sur une cuvette topographique), il a la possibilité d’augmenter la valeur d’altitude d’un patch sur lequel il se situe grâce à la quantité d’eau qu’il transporte (variable stock_en_eau ). Cette fonctionnalité, déjà utilisé par Darboux et al. (2002), permet de s’affranchir des problèmes des cuvettes fréquemment soulignés dans la littérature. Pour limiter le nombre de cuvettes présentes dans le MNT après l’intégration des réseaux de haies et de fossés, un traitement raster (fillsink) (Rakotoarisoa et al., 2014) est initié en début de simulation et permet de combler artificiellement les cuvettes.

Objectifs

18Les "agentgouttes" ont pour seul objectif de se déplacer en fonction de la plus grande pente. Tant qu’un "agentgoutte" a la possibilité de se déplacer, il le fait. Seuls les agents situés sur le patch exutoire ne peuvent pas se déplacer puisqu’il s’agit du patch ayant la plus basse altitude.

Apprentissage

19Au cours de la simulation, les "agentgouttes" enregistrent un certain nombre de données qui seront utilisées à la fin de chaque simulation pour obtenir différents indices spatiaux (par exemple, la distance parcourue). L’enregistrement de certaines données, comme l’identification du premier patch linéaire rencontré, ne peut se faire que sous certaines conditions dépendantes de l’expérience de l’agent au cours de la simulation. Par exemple, pour qu’un agent obtienne le nom du premier patch linéaire rencontré, deux conditions, dépendantes de son expérience, doivent être réunies : 1/ que celui-ci n'ait jamais rencontré de patch linéaire préalablement, 2/ que le patch sur lequel se trouve l'agent est bien identifié comme étant un patch linéaire.

Perception

20Les agents n’ont aucune connaissance de l’environnement global de simulation mais une vision limitée aux huit patches voisins (voisinage de Moore) et aux "agentgouttes" qui pourraient se trouver dessus.

Interaction

21Les "agentgouttes" interagissent avec les patches au cours de la simulation. Lorsque l’un d’eux se situe sur un patch "fossé", il se déplace préférentiellement sur un patch fossé voisin (d’altitude inférieure), si aucun patch fossé voisin n'a une altitude inférieure alors l’agent quitte le fossé et reprend un écoulement "topographique". Lorsque plusieurs "agentgouttes" se trouvent sur une cuvette ils interagissent entre eux pour la combler (cf adaptation).

Stochasticité

22La stochasticité est rare dans ce modèle car les règles sont principalement déterministes. Seule exception : dans les sous modèles Ecoul_Topo et Ecoul_Réseau, si deux patches voisins ont la même altitude, c’est par un choix aléatoire que l’"agentgoutte" se déplace sur l’un d’eux.

Détails (Details)

Initialisation

23A l’initialisation, il est possible d’intégrer ou non les contraintes linéaires. Chaque patch est ensuite affecté d’un "agentgoutte" : le nombre d’"agentgouttes" est donc identique à celui des patches. La densité d’"agentgouttes" n’est pas modifiable dans cette version. L’utilisateur a cependant la possibilité d’agir sur le stock_en_eau transporté par les agents. Bien que peu réaliste (puisque l’on peut générer de "très grosses gouttes"), cette option présente l’avantage d’éviter de générer un grand nombre d’agents et donc de diminuer la durée des simulations.

Données d’entrée

25Les réseaux linéaires (Figure 4) ont été obtenus par un travail de terrain et par photo-interprétation des orthoimages de l’IGN. S’il est possible d’obtenir facilement les réseaux de haies (données "végétation" de la BDtopo de l’IGN ou par photo-interprétation), il n’existe pas de référentiel sur les réseaux de fossés car ces entités linéaires sont construites à la seule initiative des agriculteurs et donc référencées dans aucun registre (Carluer et De Marsily, 2004). Dans un premier temps, l’ensemble des réseaux est relevé dans toute sa complexité : routes, fossés, route et fossés associés, haies et fossés associés, talus et fossé associé, haie, haies sur talus, routes associées à haies sur talus etc. De cette diversité, une classification simple est opérée. L'échelle de sortie étant le pixel à 10 mètres, un seul élément est considéré. La présence ou l’absence d’un fossé est le seul facteur discriminant. Ainsi seules deux classes sont utilisées en sortie : haie et fossé. La classe haie rassemble tous les réseaux linéaires non associés à un réseau de fossé. La classe fossé correspond, au contraire, au réseau associé à au moins un fossé (Figure 4). Cette typologie permet ensuite de modifier sous Netlogo les valeurs « altitudes » des entités concernées et d’affecter des règles d’interactions spécifiques.

Les sous-modèles

Tableau 1 : Code du sous-modèle Ecoul_Ciblé

Code	Traduction
To Ecoul_Ciblé face cible move-to cible end	Procédure Ecoul_Ciblé L’"agentgoutte" fait face à sa cible L’"agentgoutte" va sur sa cible Fin de la procédure

27Si le modèle Ecoul_Ciblé offre l’avantage d’être rapide dans son exécution, les valeurs cibles utilisées ne peuvent pas être directement modifiées sous NetLogo lors de l’intégration des réseaux linéaires. En effet, l’intégration des réseaux linéaires modifie les valeurs cibles des patches situés à proximité (puisque les altitudes des patches linéaires sont modifiées). Aussi, pour conserver rapidité (d’exécution) et cohérence (hydrologique), le modèle Ecoul_Ciblé a été couplé à Ecoul_Topo, pour intégrer les valeurs d’altitude dans son algorithme de calcul. Ecoul_Topo s’applique à tous les "agentgouttes" situés sur des patches possédant au moins un patch voisin linéaire (en jaune sur la figure 5), tandis que le modèle Ecoul_Ciblé est utilisé par les agents sur les patches non linéaires.

Tableau 2 : Code du sous-modèle Ecoul_Topo

Code

Traduction

To Ecoul_Topo  [let p one-of neighbors with-min [altitude] with-min [distance myself] ifelse [altitude] of p < altitude [ face p move-to p ] [ let profondeur_cuvette [altitude] of p - altitude let eau_disponible sum [stock_en_eau] of turtles-here  ifelse profondeur_cuvette > = eau_disponible [set altitude altitude + stock_en_eau ] [let R eau_disponible - profondeur_cuvette  let distribution_R R / count turtles-here set altitude altitude + ( profondeur_cuvette + 0.00001) set stock_en_eau stock_en_eau - (stock_en_eau -distribution_R) ]] end

Procédure Ecoul_Topo Soit p un patch voisin dont les valeurs altitude et distance sont les plus petites des 8 patches voisins Si l’altitude de p est inférieure à l’altitude de l’"agentgoutte" alors  Faire face à p  Avancer sur p… Sinon Soit profondeur_cuvette : la profondeur de la cuvette en mètre Soit eau_disponible : quantité totale d’eau transportée par les agents présents sur la cuvette Si la profondeur de la cuvette est supérieure à la quantité d’eau disponible alors L’altitude augmente de la somme des stocks d’eau des agents présents Sinon Soit R : quantité d’eau restant après comblement cuvette Soit distribution_R : R divisé par nombre d’agent présents L’altitude augmente de la valeur de la profondeur de la cuvette et de 0.00001 Le stock d’eau de chaque agent diminue de distribution_R Fin de la procédure

28Un troisième sous-modèle (Ecoul_Réseau) permet de gérer le déplacement des "agentgouttes" dans les réseaux fossés en les forçant à se déplacer en fonction de la pente du fossé.

Tableau 3 : Code du sous-modèle Ecoul_Réseau

Code

Traduction

To Ecoul_Réseau [let a one-of neighbors with [linear = 20] with-min [altitude] with-min [distance myself]  ifelse altitude > [altitude] of a [face a move-to a] [ let a1 one-of neighbors with-min [altitude] with-min [distance myself]  ifelse [altitude] of a1 < = altitude  [face a1 move-to a1] [let profondeur_cuvette [altitude] of a1 - altitude let eau_disponible sum [stock_en_eau] of turtles-here      ifelse profondeur_cuvette > = eau_disponible   [ set altitude altitude + stock_en_eau ] [ let R eau_disponible - profondeur_cuvette    let distribution_R R / count turtles-here   set altitude altitude + (profondeur_cuvette + 0.00001)  set stock_en_eau stock_en_eau - distribution_R  Ecoul_Topo ]]] End

Procédure Ecoul_Réseau Soit a un patch voisin fossé dont les valeurs altitude et distance sont les plus petites des 8 patches voisins Si la valeur d’altitude du patch sur lequel est l’"agentgoutte" est supérieure à l’altitude de a L’"agentgoutte" fait face à a L’"agentgoutte" va sur a Soit a1 un patch voisin dont les valeurs altitude et distance sont les plus petites des 8 patches voisins Si l’altitude de a1 est inférieure à l’altitude de l’"agentgoutte" alors L’"agentgoutte" fait face à a1 L’"agentgoutte" va sur a1 Soit profondeur_cuvette : la profondeur de la cuvette en mètre Soit eau_disponible : quantité totale d’eau transportée par les agents présents sur la cuvette Si la profondeur de la cuvette est supérieure à la quantité d’eau disponible alors L’altitude augmente de la somme des stocks d’eau des agents présents Sinon Soit R : quantité d’eau restant après comblement cuvette Soit distribution_R : R divisé par nombre d’agent présents L’altitude augmente de la valeur de la profondeur de la cuvette et de 0.00001 Le stock d’eau de chaque agent diminue de distribution_R (retour comblement cuvette ?) Sous-modèle Ecoul_Topo Fin de la procédure

29Il existe donc trois sous-modèles, les "agentgouttes" utilisent l’un des trois en fonction du patch sur lequel ils se situent (Figure 5).

Application du modèle et premiers résultats

Présentation du bassin de la Bouderie, support des simulations

31Les linéaires (tous confondus) sont nombreux dans ce bassin (127 mètres linéaires par hectares). Néanmoins, cet indice global, peu pertinent, masque une forte variabilité entre les extrémités du bassin et la partie intermédiaire : la densité varie entre 95 et 113 mètres linéaires par hectares dans les parties amont et aval, tandis qu'elle atteint plus de 160 mètres linéaires par hectare dans la partie centrale (Figure 6). L’occupation du sol, principalement agricole, explique ce contraste : les parties aval et amont présentent de grandes parcelles cultivées alors que dans la partie centrale, le maintien d’herbages permanents sur de petites parcelles favorise la forte densité de linéaires. Par ailleurs, la grande majorité de ces réseaux est constituée de haies sur talus associées à un fossé. La présence quasi systématique de fossés s’explique notamment par la présence de formations argileuses.

Résultats

Résultats des simulations sans linéaire (sous modèle Ecoul_Ciblé)

33En interagissant avec la valeur "cible" des patches, les "agentgouttes" se déplacent de patch en patch jusqu’à atteindre l’exutoire. Une telle simulation est proche de ce que proposent des modèles de simulation d'écoulements tels que STREAM (Cerdan et al., 2001) ou RUICELLS (Delahaye, 2002 ; Douvinet, 2008), à savoir la mise en évidence des axes d'écoulement. Les patches situés à l’amont de l’exutoire sont ainsi identifiés tout comme il est possible de connaitre le nombre d’"agentgouttes" ayant circulé sur chacun des patches (Figure 7). Cette première approche est satisfaisante puisque : 1) les agents interagissent avec les valeurs "cible" des patches pour se déplacer et atteindre l’exutoire du bassin et 2) une superficie de bassin de 162 hectares est obtenue, ce qui correspond à la superficie calculée avec d’autres outils SIG.

Résultats des simulations intégrant les linéaires

35L’impact des réseaux anthropiques sur les écoulements de surface est ensuite exploré en analysant l’évolution des distances à l’exutoire de chacun des patches composant le bassin versant. Autrement dit, on cherche à savoir si les réseaux linéaires tendent à faciliter le transfert des écoulements de surface vers le cours d’eau (ce qui pourrait, dans ce cas précis, diminuer le temps de concentration et raccourcir le temps nécessaire à la gestion des crues) ou si, au contraire, ils augmentent la distance hydrologique des écoulements depuis les versants (accentuation des piégeages des écoulements hydriques et solides, diminution du pic de crue et des pertes en terre).

36La distance à parcourir par un agent qui partirait d’un patch pour rejoindre l’exutoire est fonction de la pente mais aussi des contraintes qu’il rencontrera au cours de son trajet. Ainsi, en calculant l’évolution de la distance parcourue par chaque agent au cours d’une simulation avec et sans linéaire (ratio distance topographique / distance avec linéaires) il est possible de mesurer l’impact des réseaux linéaires sur les distances à l’exutoire de chaque patch et de dresser une typologie des espaces en fonction de l’évolution de ces distances. Trois types d’espaces peuvent ainsi être identifiés (Figure 9) :

• Tout d’abord, certains espaces sont connectés à l’exutoire dès lors que les réseaux linéaires sont intégrés (A). En l’absence de linéaires, les écoulements dont les zones d’émission se situent dans ces espaces ne rejoignent pas l’exutoire identifié mais s’écoulent vers un autre point. Ici, cela concerne 15 hectares qui sont connectés à l’exutoire et dont l’intégration n’est pas prise en compte par les outils classiques de délimitation de bassin versant (ex : outil "délimitation de bassin versant" proposé par les logiciels SIG) ;

• A l’inverse, d’autres espaces (9 hectares) sont déconnectés de l’exutoire (B) lorsque l’on intègre les réseaux linéaires. Cela s’explique par la présence de haies qui bloquent les écoulements ou bien par la présence de réseaux de haies et de fossés qui dévient les écoulements vers des bassins voisins (cas des espaces situés sur les marges du bassin) ;

• Enfin, la grande majorité des espaces du bassin, soit 140 hectares, reste connectée à l’exutoire après intégration des linéaires, mais la distance à parcourir par les agents pour atteindre l’exutoire est modifiée (C). Le gradient de couleur sur la carte (Figure 9) représente l’évolution de ces distances. La couleur rouge met en évidence les espaces dont la distance à l’exutoire a augmenté (la distance à parcourir par un agent est plus grande), en vert, les espaces dont la distance a diminué, enfin en jaune, les espaces dont la distance a peu évolué. Les espaces représentés en rouge ont une distance à l’exutoire qui a augmenté de plus de 10 % (en pourcentage de la distance initiale), avec un maximum de + 25 %. Les espaces en vert ont une distance qui a diminué jusqu'à 19 % de la distance initiale. Globalement, sur le bassin étudié, la majeure partie des patches a une distance à l’exutoire qui augmente du fait de la présence de réseaux linéaires.

37La cartographie de l’évolution des distances à l’exutoire permet donc d’appréhender l’impact de la structure paysagère sur les écoulements de surface. Il est ainsi possible de discriminer les espaces d’un bassin versant en fonction de l’impact des réseaux linéaires sur le parcours des écoulements qu’ils émettent. Les résultats offrent l’avantage d’intégrer explicitement l’orientation des réseaux par rapport aux pentes ainsi que leur connectivité les uns par rapport aux autres. La comparaison de cette cartographie avec celle présentant les densités de linéaires est significative (Figure 10). Sur les focus 1 et 2, de densité équivalente avec plus de 18 mètres de linéaires par 100 m², il apparaît très clairement que l’impact sur les distances à l’exutoire n’est pas le même. Dans le focus 1, les réseaux, principalement de fossés et orientés dans le sens de la pente, diminuent les distances hydrologiques à l'exutoire, alors que dans le focus 2, les distances augmentent du fait de la présence de haies perpendiculaires à la pente qui dévient les écoulements. Cette cartographie illustre les propos évoqués par de nombreux auteurs : plus que la densité, c’est l’organisation et la structuration dans l’espace des réseaux qui influence le plus les écoulements de surface (Mérot, 1998 ; Delahaye, 2002 ; Baudry et Jouin, 2003 ; Douvinet, 2008 ; Viel, 2014).

38Cette cartographie permet ainsi de visualiser un premier indice obtenu grâce à la simulation multi-agents. Basé sur peu de données en entrées (MNT + réseaux linéaires) et transposable, cet indice fournit une première estimation du poids des réseaux sur les écoulements de surface au sein d'un petit bassin versant agricole. Ce dernier, en plus de livrer une première appréciation de la rugosité d'un paysage, est construit à partir d'informations spatiales pertinentes à une échelle inférieure. Il ouvre également des perspectives sur la capacité des systèmes multi-agents à proposer, en tant qu’outil d’aide à la décision, des indices simples de caractérisation de l’efficacité des réseaux sur le parcours des écoulements de surfaces : distance des parcelles à un réseau linéaire (permet d’estimer les temps de concentration dans un réseau), poids de la circulation dans un réseau dans le déplacement des agents, points de passage d'un écoulement de surface à un écoulement de réseau, etc. De plus, les perspectives en matière de recherche sont importantes notamment autour du concept de connectivité. Ce concept, de plus en plus utilisé en hydrologie et géomorphologie, peine encore à être clairement défini et surtout à être caractérisé par des indices simples et transposables. Le modèle présenté ici montre toute son acuité pour être un outil pertinent dans ce sens. Cependant, la version actuelle ne prend pas en compte les vitesses d'écoulement et le type de surface, et ne permet donc pas de pouvoir mesurer finement les effets des réseaux sur les temps de parcours des écoulements. Car si ce premier indice permet d'évaluer l'évolution des distances, d'autres phénomènes ne sont pas pris en compte (accélération des vitesses en fonction de la nature du sol, effet d'étranglement dans les fossés, etc.) alors que leurs effets sur les temps de parcours sont très importants.

39Plusieurs difficultés ont été rencontrées durant la phase de développement du modèle. Nous les présentons ici brièvement pour conforter les choix qui ont été faits :

1. Grâce à différents établissements publics tels que l’Institut National Géographique les données nécessaires au modèle sont, en partie, accessibles (MNT, réseau de haies, occupation du sol…). En revanche, il est compliqué de s’affranchir du terrain pour recenser les réseaux de fossés. Par ailleurs, si les orthoimages de l’IGN offrent de bonnes informations sur l’utilisation du sol, un travail de terrain pour connaître les différents types de cultures en place doit être envisagé. Ainsi, tant que la technologie ne permettra pas d’obtenir rapidement et facilement toutes ces données, la méthode présentée restera semi-automatique. Les outils de la télédétection restent une piste ouverte pour automatiser le travail d’acquisition de données ;

2. • 11 Initialement le MNT ne contient aucune cuvette car "nettoyé" sous SIG. En revanche, après intégrati (...)

   L’intégration des réseaux linéaires dans la simulation modifie l’altitude des patches concernés, ce qui modifie directement la structure du MNT sous-jacent. Ainsi à l’initialisation, le nombre de cuvettes (patches ayant une altitude inférieure à leurs 8 voisins) est de 95011. Pour limiter ce nombre, un algorithme de calcul a été codé directement sous NetLogo (Rakotoarisoa et al., 2014) afin de diminuer le nombre de cuvette à une cinquantaine, soit 0.1 % du nombre total de patches. Même en réduisant à une cinquantaine de cuvettes, le MNT reste hydrologiquement incorrect et justifie le développement du sous-modèle Ecoul_Topo, qui permet aux agents de gérer collectivement le problème des cuvettes (cf : adaptation) ;

3. Dans cette première approche, seuls les réseaux linéaires ont été intégrés. Il apparaît désormais indispensable d’intégrer d’autres entités du paysage (buse enterrée, mares, entrés de champs…) et de l’occupation du sol pour tendre vers des simulations capables de mesurer avec finesse le poids du paysage sur les écoulements de surface.

Conclusion

40La structure paysagère et plus particulièrement les réseaux linéaires jouent un rôle important sur le parcours des transferts hydro-sédimentaires sur les versants en modifiant les chemins d’écoulements. Pour appréhender le rôle de ces réseaux sur les dynamiques d’écoulement un outil de modélisation sous une plateforme multi-agents (NetLogo) a été développé. En discrétisant les écoulements par un ensemble d’agents, dont on modélise le comportement, il est possible de reconstruire les écoulements de surface et de mesurer l’impact des éléments linéaires. Après s’être assuré du bon fonctionnement du modèle, un premier indice d'analyse spatiale a été présenté, à savoir l’impact des réseaux sur les distances à parcourir par les écoulements pour rejoindre l'exutoire d'un bassin. Les premières simulations montrent toute la pertinence des SMA comme outil d’évaluation du rôle de la structure paysagère sur les écoulements de surface. En analysant le rôle de chaque linéaire sur le parcours des écoulements, le modèle permet de démontrer que plus que la densité de linéaire dans l’espace, c’est leur organisation et leur structuration qui doivent être analysées. Cette hypothèse souvent expliquée de manière empirique trouve ici une réponse et offre surtout une visualisation de ces fines dynamiques spatiales. En cela les SMA sont des outils parfaitement adaptés. L’intégration prochaine des vitesses d'écoulement et de l’occupation du sol devrait permettre d’améliorer les résultats actuels notamment en définissant les zones susceptibles d’émettre du ruissellement et celles pouvant les limiter (prairie, bandes enherbées, etc.). A ce titre, le modèle développé est aussi un véritable outil d'aide à la décision, pouvant permettre de guider les changements opérés sur l'occupation du sol, de visualiser l'incidence de ces changements sur le fonctionnement global du bassin versant, ou même, dans une logique expérimentale, d'évaluer le rôle de l'organisation spatiale des réseaux et des aires productrices de ruissellement sur la réponse hydrologique et sédimentaire du bassin versant.