Cette étude a bénéficié d’un financement de la région Occitanie et de l'Université Fédérale Toulouse Midi-Pyrénées. Les auteurs remercient Johann Blanpied et Tiberiu Groparu pour leur aide lors de l’inventaire de terrain. Sont également remerciés les deux relecteurs anonymes pour leurs suggestions ayant contribué à l’amélioration du manuscrit.
1Les rivières à fond rocheux sont caractérisées par une charge alluviale minimale et/ou transitoire dans un contexte où la capacité de transport des sédiments est excédentaire par rapport à la fourniture sédimentaire sur une période de temps long (Howard, 1994 ; Whipple et al., 2000). Ce type de rivières est présent dans différents contextes géodynamiques (baisse du niveau de base, surrection tectonique par faille ou plissement, etc.) et géomorphologiques (montagne, piémont, plaine). Dans le cas des rivières à fond rocheux de plaine, le rôle des impacts anthropiques dans la réduction des flux solides peut fréquemment être évoqué à l’image des rivières Ropa en Pologne (Rinaldi et al., 2005), Semois en Belgique (Gob et al., 2005), Rambla de la Viuda en Espagne (Calle et al., 2017) ou encore Garonne en France (Jantzi et al., 2018a). Les extractions de granulats représentent la cause majeure de leur déficit sédimentaire. Fréquemment dans ce type de cours d’eau, les affleurements rocheux sont discontinus et peuvent prendre la forme de seuils ou de platiers rocheux. Ces formes, immergées une partie de l’année hydrologique, présentent alors un modelé d’érosion diversifié dont la taille peut varier entre quelques centimètres à plusieurs centaines de mètres (Bradwell, 2005 ; Richardson et Carling, 2005). La formation et l’évolution de ces rivières répondent à une combinaison de facteurs liés à l’hydraulique (type d’écoulement, etc.), à la lithologie (nature de la roche, organisation et homogénéité/hétérogénéité, pendage, présence de failles, etc.) et aux processus d’érosion et d’altération (mécanique, biologique et chimique) (Whol, 1998 ; Hall et al., 2012 ; Inoue et al., 2017 ; Gomez-Heras et al., 2019). Les mécanismes d’érosion façonnant la morphologie des rivières à fond rocheux au travers de formes d’érosion ont été déclinés en six types qui sont (Hancock et al., 1998 ; Whipple et al., 2000 ; Richardson et Carling, 2005) :
-
l’abrasion (i.e., abrasion) produit une usure de la surface et une perte de matière sous l’impact de sédiments en mouvement transportés par charriage et/ou en saltation.
-
l’arrachage (i.e., plucking) entraîne le soulèvement puis le glissement de blocs délimités par des fractures ou diaclases dans la roche ou l’écaillage de roche saine.
-
la cavitation (i.e., cavitation) peut produire une onde de choc liée à la formation de bulles de vapeur dans un fluide lorsque la pression du fluide chute localement en dessous de la pression de vapeur. Elle est susceptible de provoquer un éclatement de la roche.
-
la dissolution (i.e., dissolution) concerne prioritairement des roches solubles, en premier lieu les carbonates, sous l’effet de l’acidité de l’eau. Ce processus peut affecter de manière plus marginale d’autres roches (dissolution de la silice par exemple).
-
le calage hydraulique (i.e., hydraulic wedging) contribue à l’ouverture de fractures préexistantes par des particules sédimentaires faisant office de coin (sables, graviers, galets). Les fissures vont alors s’ouvrir progressivement sous l’effet de la pression du fluide sur les particules contribuant ainsi à disloquer la roche en plusieurs blocs.
-
l’altération mécanique (i.e., physical weathering) entraînant un affaiblissement physique de la roche avec la propagation de fissures et sa désagrégation par différents processus de fragmentation (thermoclastie, hydroclastie) favorisant son érosion par les processus cités aux points précédents. Les roches stratifiées sont généralement particulièrement sensibles à ce type de mécanisme.
2À ces différents mécanismes il faut rajouter le rôle de la végétation qui peut s’installer temporairement (quelques mois à quelques années) mais aussi le rôle des biofilms et des alternances humectation/dessiccation qui ont fait jusqu’à présent l’objet de peu de travail (Seiffert et al., 2016 ; Wieler et al., 2019).
3Parmi ces différents processus d’érosion, l’abrasion et l’écaillage sont considérés comme les plus importants avec pour le premier une prédominance pour les roches fracturées présentant des discontinuités et pour le second une prédominance dans les roches massives faiblement fracturées (Hancock et al., 1998 ; Whipple et al., 2000 ; Chatanantavet et Parker, 2009 ; Lamb et Fonstad, 2010). Les processus de dissolution et de météorisation ont été peu étudiés mais peuvent représenter des facteurs importants d’érosion (Büdel, 1982 ; Lamb et Fonstad, 2010 ; Pelletier et Baker, 2011). La dégradation de la molasse par dissolution en moyenne Garonne a fait l’objet d’une analyse et a montré que le processus contribuait de manière significative à l’érosion des seuils rocheux (Jantzi, 2018). Quant au processus de cavitation, sa contribution dans l’érosion des formes, si elle apparaît comme probable, ne saurait être précisé en raison d’un manque d’étude à ce sujet et des difficultés à observer ce type de processus in situ (Whipple et al., 2000 ; Carling et al., 2009).
4Les formes d’érosion observées dans les lits mineurs des rivières à fond rocheux sont d’une très grande diversité et ont fait l’objet de travaux visant à construire une typologie à partir d’approches morphométriques (Allen, 1971 ; Richardson et Carling, 2005). Les trous d’eau représentent l’une des formes les plus répandues, communément observées et qui ont retenu l’attention (Alexander, 1932 ; Allen, 1982 ; Zen et Prestegaard, 1994 ; Lorenc et al., 1994 ; Richardson et Carling, 2005 ; Springer et al., 2005, 2006 ; Ortega et al., 2014 ; Lima et Binda, 2015 ; Pelletier et al., 2015). Il s’agit de formes en creux cylindriques ou hémisphériques, pouvant mesurer de quelques centimètres à plusieurs mètres de diamètre et de profondeur (Whipple et al., 2000 ; Springer et al., 2005). Elles se développent aussi bien dans des substrats tendres comme l’argile ou résistants comme le granite. Les processus érosifs à l’origine des trous d’eau ne sont pas complètement compris à ce jour, mais les mécanismes d’arrachage, d’abrasion et de cavitation contribueraient au façonnement de ces formes (Hancock et al., 1998 ; Sklar et Dietrich, 1998 ; Dubinski et Wohl, 2013). L’abrasion est cependant mise en avant comme le principal processus de leur formation (Hancock et al., 1998 ; Whipple et al., 2000). Le trou d’eau se développe par érosion de sa paroi et de son fond sous l’effet du déplacement de sédiments piégés à l’intérieur (sables, graviers ou galets). Ces derniers entraînés par un écoulement tourbillonnant, vont abraser la roche et/ou arracher des fragments par leur impact (Springer et al., 2005). Si le charriage fin et grossier est perçu comme étant en grande partie responsable de la croissance des trous d’eau, Whipple et al. (2000) ont montré que l’abrasion par les sédiments en suspension pouvait également jouer un rôle important dans leur développement. Ces formes sont ainsi caractéristiques des rivières à haute énergie (Wohl, 1993 ; Montgomery et al., 1996 ; Whipple et al., 2000).
5L’existence de ces formes pose plusieurs questions d’intérêt pour la compréhension des phénomènes d’évolution des chenaux fluviaux : (i) quel est leur rôle dans le processus général d’incision du chenal à fond rocheux ? (ii) à quelle vitesse et sur quelle durée se développent-elles ? (iii) les différents types de formes sont-elles ubiquistes, uniquement influencées par la lithologie ou sont-elles organisées spatialement en réponse à la variabilité des conditions d’écoulement, ce qui permettrait de les considérer comme des indicateurs du fonctionnement et de l’évolution du chenal.
6Leur étude peut permettre d’identifier d’une part l’organisation des écoulements et d’autre part les processus d’érosion intervenant dans l’évolution des lits fluviaux rocheux et notamment des seuils rocheux. Pour cela, l’analyse morphométrique des formes est une étape importante afin de définir leur origine (Gomez-Heras et al., 2019). Elle représente le principal outil permettant de discriminer les formes, dont la diversité et la complexité sont liées à de multiples facteurs (caractéristiques hydrauliques de l’écoulement, lithologie de la roche, présence de fractures, etc.). Dans ce travail, nous proposons d’étudier les formes d’érosion développées sur un substrat de type molassique en moyenne Garonne. À ce jour, peu de travaux ont été menés sur des roches sédimentaires stratifiées comme les molasses.
7La Garonne est un fleuve situé dans le sud-ouest de la France drainant un bassin versant de 56 000 km2 (fig. 1). Les trois sites d’études (Portet-sur-Garonne, Beauzelle et Grenade) se situent dans la partie médiane de la Garonne, entre la confluence avec l’Ariège à Portet-sur-Garonne en amont et la confluence avec le Tarn à Moissac en aval, soit un linéaire d’environ 80 km (fig. 1). Sa surface drainée est de 15 700 km2 (soit 28 % du bassin versant de la Garonne) pour une pente moyenne de 0,9 m.km-1 et une puissance spécifique moyenne estimée à 70 W.m-2 (tab. 1) (Steiger et al., 2000). Le cours d’eau présente sur ce tronçon un régime hydrologique de type pluvio-nival (Pardé, 1928 ; Lambert, 1989). Le débit moyen annuel est de 190 m3.s-1 à Portet-sur-Garonne (Banque Hydro) et l’année hydrologique se caractérise par des basses eaux de juillet à décembre et des hautes eaux de janvier à juin, période durant laquelle se produit l’essentiel des crues.
Fig. 1 – Carte de localisation des sites d’études
Fig. 1 – Location map of the three study sites
(BD Alti IGN, Photos SMEAG)
(BD Alti IGN, Photos SMEAG).
8Les sites, caractérisés par la présence de seuils rocheux, s’inscrivent dans la plaine holocène (David et al., 2016) large de quelques kilomètres où le chenal s’est incisé dans un substratum sédimentaire molassique de type fluvio-lacustre (Mouline, 1989 ; David et al., 2015). Ce dernier est issu de l’érosion de la chaîne des Pyrénées. Cette formation molassique, datée entre l’Oligocène et le Miocène, est constituée de marnes et d’argiles plus ou moins sableuse où s’intercalent localement des bancs calcaires indurés (Sunyach, 1984). Cette différence de dureté entre strates, induit une érosion différentielle créant un étagement en paliers des seuils rocheux. En périodes d’étiage, ces derniers peuvent couvrir une surface de plusieurs hectares et engendrer une rupture de pente de plus de 2 m.
9Ces formes ont évolué rapidement au cours de la seconde moitié du 20e siècle, essentiellement sous l’effet d’interventions anthropiques dans le lit mineur entre les années 1960 et 1980. La mise en place de barrages hydro-électriques en amont, d’ouvrages de fixation du chenal (digues, enrochements, etc.), ainsi que les extractions de granulats, ont contribué au tarissement des apports solides et à un déficit sédimentaire en moyenne Garonne. Cela s’est traduit par un enfoncement moyen du lit de 1 m avec des valeurs pouvant atteindre localement 2,7 m (Jantzi et al., 2017a ; Jantzi, 2018). Parmi ces activités, les extractions qui ont eu lieu sur chacun des sites, représentent le principal facteur d’apparition des seuils rocheux. Ces dernières, ont entrainé un décapage de la charge alluviale par érosion régressive, laissant apparaître des affleurements molassiques. Les processus d’incision changent alors de nature et passent de l’évacuation des alluvions à l’érosion du substratum. Dans le cas du seuil de Beauzelle, l’incision a atteint 2,8 m en moyenne et jusqu’à 4,5 m localement.
Tab. 1 – Caractéristiques géomorphologiques des sites d’études
Tab. 1 – Geomorphological settings of the study sites
L : longueur du chenal, Wm : largeur moyenne du chenal, SI : sinuosité, Zm : profondeur moyenne du chenal, Sm : pente moyenne du chenal, Ω : puissance spécifique
L: channel length, Wm: mean channel width, SI: sinuosity, Zm: mean channel depth, Sm: mean channel slope, Ω: specific stream power.
10Un inventaire des formes d’érosion pour chaque site a été effectué en 2017 afin d’établir une typologie de ces dernières. Les relevés ont été réalisés au mois de juillet c’est à dire en étiage, afin de permettre la prospection des surfaces les plus étendues possibles. Environ 800 formes ont été relevées à partir d’un échantillonnage systématique. La méthode a consisté en un relevé des formes le long de transects perpendiculaires à l’axe du chenal et matérialisés par un fil tendu tous les 10 m. Chaque forme intersectée a fait l’objet d’une série de relevés qui comporte (i) sa localisation (x, y, z) par dGPS (Trimble Geo 7X ± 1 cm en X et Y), (ii) son orientation (azimut, mesuré à la boussole à bain d’huile Chaix), (iii) la présence de fissures et leur orientation, (iv) la présence de particules (galets, graviers et sables) et la mesure de leur largeur (axe b), (v) sa géométrie (largeur, longueur, profondeur), et (vi) son rattachement à un type de formes sur la base de la typologie de Richardson et Carling (2005) qui, dans certains cas, a été précisée et/ou adaptée ultérieurement à partir de photographies prises lors des relevés.
11La localisation des formes d’érosion relevées par dGPS a été intégrée dans un SIG (sous ArcGis) afin de permettre une analyse de leur répartition spatiale. Dans un premier temps, une approche descriptive a été adoptée à partir d’une analyse cartographique de la répartition des différents types de formes permettant de rendre compte du contexte hydro-géomorphologique dans lequel elles s’inscrivent. Dans un second temps, une approche quantitative de la répartition des formes a été adoptée. Une première analyse a consisté en un calcul de densité par noyau (kernel) permettant de produire une cartographie des densités de formes.
12En complément, une analyse par quadrats a permis une cartographie des densités sous la forme d’un maillage et de déterminer le type de distribution spatiale des formes (aléatoire, agrégative, homogène). Les sites d’étude ont été découpés selon un maillage carré régulier (10 m x 10 m) pour dénombrer les formes et déterminer leur répartition. Le mode de distribution est appréhendé par un test d’hypothèse (Ho) consistant à comparer une distribution observée avec une distribution théorique correspondant à une hypothèse selon laquelle la répartition des individus est le résultat d’un processus particulier (Taylor, 1977). Dans le cas présent, il s’agit de comparer la distribution observée par rapport à une distribution aléatoire (distribution de Poisson, hypothèse Ho) pour laquelle la distribution observée est le résultat d’un processus stochastique. Si Ho est rejetée, il peut être admis que la distribution présente une tendance significative à l’agrégation et/ou à la dispersion des individus.
13Le type de distribution spatiale a été déterminé à partir de l’indice de dispersion (ID) et l’indice de Morisita (IM), qui sont calculés de la manière suivante :
14où :
-
s2 correspond à la variance de la densité,
-
m à la densité moyenne d’individus par maille,
-
k au nombre total de mailles,
-
à la somme des individus par quadrat,
-
à la somme des carrés des individus par quadrat.
15Pour les deux indices, la distribution est considérée comme aléatoire si le résultat est égal à 1, régulière s’il est inférieur à 1 ou agrégative s’il est supérieur à 1. Afin de valider les résultats des indices ID et IM et pouvoir rejeter Ho, le test du Chi-2 (X2) et le test Z ont été réalisés.
16Les trous d’eau étant la forme la plus commune et la mieux définie sur le terrain, une analyse statistique de leurs caractéristiques géométriques a été réalisée. La taille et la forme des trous d’eau échantillonnés étant diversifiées (Jantzi, 2018), seuls les trous d’eau simples et individuels (sans structure secondaire et lien avec des formes voisines), avec une ouverture circulaire ou elliptique et dont la surface interne est relativement lisse, ont été sélectionnés pour cette analyse (fig. 2).
Fig. 2 – Exemple d’un trou d’eau circulaire (A) et d’un trou d’eau elliptique (B)
Fig. 2 – Example of circular pothole (A) and elliptical pothole (B)
Notez la présence de galets (A) et de sables grossiers (B) dans le fond des trous d’eau.
Note the presence of pebbles (A) and coarse sand (B) in the bottom of potholes.
17Afin d’obtenir un jeu de données plus conséquent que celui obtenu par l’échantillonnage systématique, un échantillonnage exhaustif ciblé sur ce type de formes a été réalisé en complément. Toutefois pour cette analyse, seuls les sites de Portet-sur-Garonne et Grenade sont pris en compte en raison d’un nombre de formes insuffisant sur le site de Beauzelle pour obtenir un résultat significatif.
18La caractérisation géométrique de chaque trou d’eau a été réalisée à partir de la mesure du diamètre d’ouverture sur deux axes perpendiculaires D1 et D2, puis de la profondeur Z. A partir de ces mesures, plusieurs indicateurs ont été calculés qui sont (i) le rapport d’aspect d’ouverture (Ro = D1/D2), (ii) le diamètre moyen d’ouverture (D), et (iii) le rapport diamètre-profondeur (Y = D/Z) qui forme une relation indépendante de l’échelle spatiale (Pelletier et al., 2015). Ce rapport permet de caractériser en première approche la dynamique des trous d’eau (Abbott et Pottratz, 1969). La profondeur mesurée correspond à la profondeur maximale de la forme après l’avoir vidée des sédiments. La mesure de l’axe b de tous les sédiments présents dans les trous d’eau étant impossible, seule la plus grosse particule a été mesurée pour chaque forme contenant des sédiments.
19Si le nombre de formes échantillonnées entre Portet-sur-Garonne (nb = 370) et Grenade (nb = 356) est quasi-similaire, l’échantillon pour Beauzelle est bien moins important (nb = 86). Ce nombre relativement faible s’explique : (i) par une surface d’affleurement rocheux bien moins importante par rapport aux autres sites, et (ii) par une quantité importante d’algues dans les zones immergées rendant les formes difficiles à identifier. Cependant, les densités de formes entre sites sont similaires avec 0,03 formes par mètre carré à Portet-sur-Garonne et Grenade et 0,02 formes par mètre carré à Beauzelle. Il existe également une hétérogénéité de la présence des formes entre les seuils. En effet certaines formes, comme les trous d’eau et les sillons sont faiblement représentées à Beauzelle.
20L’inventaire par échantillonnage systématique a permis d’identifier sept groupes de formes qui sont (fig. 3A) :
-
les trous d’eau, qui correspondent à une dépression généralement fermée, de forme circulaire ou elliptique en plan, à paroi plus ou moins verticale (Richardson et Carling, 2005). Le diamètre et la profondeur peuvent aller de quelques centimètres à 50 cm environ.
-
les sillons longitudinaux, qui sont des formes allongées et orientées parallèlement à la direction de l’écoulement local et de manière plus générale à celle du chenal (Richardson et Carling, 2005).
-
les flûtes, qui se caractérisent par un aspect en plan parabolique (Allen, 1971). La partie supérieure (amont) présente généralement un profil concave à forte pente tandis que la partie inférieure (aval) présente une pente plus douce, la profondeur étant plus importante près de l’extrémité amont (Richardson et Carling, 2005). Leur bordure est marquée en zone proximale (extrémité amont), soit de manière arrondie ou angulaire, pour s’effacer vers l’aval et se confondre avec la surface générale (Richardson et Carling, 2005). Leur profondeur peut atteindre une dizaine de centimètres pour une longueur et une largeur n’excédant pas 50 cm.
-
les flaques, qui correspondent à une dépression peu marquée de faible profondeur, inférieure à 10 cm, avec une forme en plan variée (elliptique, allongée, multilobée, etc.). Le fond est relativement plat et les bords sont peu marqués, en pente douce et fuyant de profil convexe. Leur taille est variable avec des dimensions allant de plusieurs dizaines de centimètres à quelques mètres.
-
les dépressions, qui sont des formes marquées de profondeur supérieure à 10 cm et pouvant atteindre plusieurs dizaines de centimètres. La forme en plan présente un aspect varié (elliptique, allongée, multilobée, etc.). Leur fond peut être plat mais également constitué de formes secondaires (trou d’eau, sillon) et se caractérise par des bords arrondis ou vifs bien marqués, de profil concave avec des parois abruptes. Leur dimension peut aller de quelques décimètres à plusieurs mètres.
-
les échancrures, qui s’apparentent à un renfoncement le long d’un rebord d’érosion et dont la forme, en V ou en U, en plan peut être large ou étroite avec des dimensions variables allant de quelques dizaines de centimètres à plusieurs mètres en largeur et en longueur. La partie amont de la forme se caractérise par une rupture de pente marquée, au profil concave, et des bords pouvant être arrondis ou vifs avec des parois abruptes. Le plus souvent, les échancrures s’observent dans le sens de l’écoulement (ouverture vers l’aval) mais peuvent également se développer à contre-sens de ce dernier (ouverture vers l’amont).
-
les buttes, sont des reliefs positifs qui se caractérisent généralement par une surface arrondie et plus ou moins lisse avec des dimensions en largeur et en longueur d’ordre métrique à plurimétrique ainsi qu’une hauteur pouvant être de plusieurs dizaines de centimètres. Elles sont la relique d’un ancien niveau pour lequel l’intensité de l’érosion est moins importante en raison d’un différentiel de dureté entre une couche supérieure indurée (encroûtement calcaire) et une couche inférieure plus tendre (marno-sableuse).
21À cela s’ajoute trois autres groupes identifiés à partir de l’échantillonnage aléatoire mais qui n’ont pas été intégrés dans l’analyse quantitative (fig. 3B) :
-
les formes ondulantes, qui correspondent à des surfaces d’aspect bosselé, relativement arrondies et lisses, avec une orientation longitudinale dans le sens de l’écoulement principal ou non-directionnelle. Elles s’observent en groupe et essentiellement à proximité du chenal principal avec des dimensions variables allant de quelques dizaines de centimètres à plusieurs mètres.
-
les fosses de dissolution, qui sont de petites dépressions d’aspect varié (circulaire, allongé, irrégulier). La taille des formes peut aller de quelques millimètres à plusieurs décimètres. Elles peuvent être isolées, en groupes ou agglomérées par la fusion de plusieurs formes.
-
les formes concaves avec surplomb, qui s’inscrivent sur un plan vertical ou approchant et forment un surplomb. Elles peuvent être profondes avec un surplomb horizontal de la paroi supérieure ou moins profondes avec un surplomb courbé entre la paroi supérieure et la paroi du fond de la forme.
Fig. 3 – A : Effectif et distribution des formes d’érosion pour chaque groupe identifié par échantillonnage systématique. B : Groupes de formes identifiés par échantillonnage aléatoire
Fig. 3 – A: Number and distribution of bedforms for each group identified by systematic sampling. B: Groups of bedforms identified by random sampling
1. Portet-sur-Garonne (nb = 369) ; 2. Beauzelle (nb = 86) ; 3. Grenade (nb = 356) ; 4. Total (nb = 811)
1. Portet-sur-Garonne (#369); 2. Beauzelle (#86); 3. Grenade (#356); 4. Total (#811).
22Ces 10 groupes de formes peuvent se décliner en 52 sous-groupes de formes différentes. Le détail de la présentation des groupes et sous-groupes est décrit par Jantzi (2018). Pour l’ensemble des sites, les trous d’eau et les sillons sont les formes les plus présentes avec respectivement 38 % et 41 % du total. Les flaques et les dépressions représentent respectivement 6 % et 9 %, tandis que les autres catégories (flûtes, échancrures et buttes) représentent 6 % (fig. 3A).
23La comparaison inter-sites met en évidence une variabilité dans la répartition de ces formes (fig. 3A). Les sillons, par exemple, forment le groupe le plus présent sur le site de Grenade avec 61 % tandis qu’à Portet-sur-Garonne et Beauzelle, leur part est deux fois moins importante. A l’inverse, les trous d’eau sont plus fréquents à Portet-sur-Garonne avec 49 % contre 33 % à Grenade et deux fois moins à Beauzelle. Ces variations entre groupes et sites sont interprétées comme le résultat des variations dans les faciès sédimentaires de la molasse et des conditions d’affleurement du substratum d’une part, et des conditions d’érosion d’autre part, mais sans qu’il soit possible de déterminer l’importance respective de chaque facteur.
24L’approche descriptive par analyse cartographique de la répartition des types de formes met en évidence : (i) une localisation préférentielle selon les formes en lien avec les conditions hydro-géomorphologiques dans lesquelles elles s’inscrivent ; (ii) les secteurs où se produisent les processus d’érosion dominant qui interviennent dans l’évolution des seuils. Il est ainsi possible de distinguer les observations générales suivantes.
25Les trous d’eau se présentent généralement dans des zones de convergence d’écoulements et/ou dans des secteurs où le relief présente des irrégularités marquées par une rugosité importante favorable à la formation de courants turbulents. La présence de sédiments dans ces cavités est également courante. Ainsi, il peut être suggéré que les trous d’eau sont entretenus principalement par l’action de flux tourbillonnants, leur morphologie favorisant le phénomène, et par l’abrasion.
26Les sillons se localisent majoritairement près d’un écoulement secondaire. C’est le cas notamment pour le site de Portet-sur-Garonne et Grenade. Ces formes ont pour effet, à faible niveau d’eau, de canaliser l’écoulement avec une hauteur d’eau plus importante et une vitesse plus élevée par rapport à l’écoulement sur les surfaces adjacentes, accentuant l’effet érosif de l’écoulement. Par ailleurs, les sillons s’accompagnent très fréquemment de sédiments grossiers suggérant également une érosion par abrasion.
27Les flaques et dépressions sont essentiellement localisées dans les secteurs au relief faiblement disséquée, où de l’eau stagnante est présente régulièrement. Par ailleurs, on note une absence de piégeage sédimentaire dans le cas des flaques. De forme étalée, à la surface peu rugueuse et de faible profondeur, les flaques semblent se développer essentiellement sous l’effet de la dissolution. Il peut en être de même pour les dépressions, qui semblent correspondre à un stade d’évolution plus avancé de flaques.
28Les échancrures sont principalement observées sur les rebords d’érosion qui correspondent à des ruptures de pentes abruptes, à l’interface entre deux niveaux pouvant être une alternance entre un faciès dur et tendre. Les échancrures se forment généralement dans le sens de l’écoulement et peuvent traduire un processus d’érosion régressive, notamment sur le site de Beauzelle.
29Les résultats de l’approche quantitative illustrés par les cartes de densités, obtenues par la méthode des noyaux (fig. 4A), mettent en évidence une distribution hétérogène avec des secteurs à la densité contrastée pour chacun des sites. Cette observation est également confirmée par les densités issues de la méthode des quadrats (fig. 4B). Globalement, la répartition des formes répond aux observations suivantes.
Fig. 4 – Méthodes de détection et cartographie des formes d’érosion
Fig. 4 – Methods of detection and mapping of erosional forms
A : Carte de densité obtenue par la méthode du noyau avec un rayon d’analyse de 15 m. B : Carte de densité obtenue par la méthode des quadrats avec un maillage de 10 m. La flèche bleue représente la direction de l’écoulement principal
A: Density map obtained by the kernel method with an analysis radius of 15 m. B: Density map obtained by quadrat method with a grid size of 10 m. The blue arrow represents the direction of the main flow
30À Portet-sur-Garonne (fig. 4-5A), trois secteurs sont identifiés. Un secteur amont qui présente une densité de formes élevée à modérée, avec des concentrations plus fortes observées principalement en première moitié de zone et en bordure de berge. Dans la partie centrale, les densités les plus fortes sont observées en bord de chenal, tandis qu’en bord de berge la densité est plus faible avec quelques zones à la densité modérée mais très localisée. Il est à noter la présence d’une bande transversale au centre du secteur sur toute sa largeur, dont la densité est faible. Le troisième secteur (en aval), se caractérise par une première moitié amont à la densité élevée à modérée. A l’inverse, la seconde moitié du secteur, présente une densité faible.
31À Beauzelle (fig. 4-5B), trois secteurs peuvent également être distingués. En amont, la densité est faible avec quelques zones à densité modérée en partie supérieure du secteur. En partie centrale, la densité élevée à modérée peut s’expliquer par la proximité du chenal principal d’écoulement et la situation en dévers induisant des vitesses d’écoulement plus élevées. En aval, la densité est modérée avec néanmoins en bordure de berge une présence assez élevée de formes qui peut s’expliquer par la topographie irrégulière et marquée induisant localement un écoulement turbulent. Comme pour Portet-sur-Garonne, les secteurs à densité élevée correspondent à une proximité au chenal d’écoulement et à des altitudes plus faibles. En plus d’une durée d’immersion plus importante lors de hautes eaux, les contraintes hydrauliques (hauteur d’eau et vitesse) y sont plus importantes. De plus, la partie centrale se trouve à un point de convergence d’écoulements de directions différentes.
32À Grenade (fig. 4-5C), le seuil se présente en deux parties distinctes séparées par un chenal d’écoulement. En rive gauche, le secteur en bord de berge présente une densité faible notamment en amont qui peut s’expliquer par son éloignement au chenal d’écoulement, par une altitude plus élevée, suggérant des contraintes hydrauliques moins importantes, et par la présence de végétation en amont réduisant la vitesse d’écoulement. À l’inverse, la partie centrale montre une densité élevée mais qui tend à diminuer vers l’aval. Cette zone étant plus proche du chenal, les contraintes hydrauliques sont supposées plus importantes. Par ailleurs, de nombreux sillons se trouvent dans l’axe d’un écoulement secondaire empruntant une rupture dans le front du seuil en période de hautes eaux. En bord de chenal, la densité est modérée mais relativement homogène à l’inverse de la partie centrale. Le nombre moins important de formes peut s’expliquer par le fait que ce secteur ne se trouve pas directement dans l’axe d’un écoulement, malgré une altitude plus faible et des contraintes hydrauliques potentiellement plus élevées. En rive droite, une similarité dans la répartition spatiale des formes ainsi qu’une certaine symétrie de situation par rapport au secteur en rive gauche, sont observées. Ainsi, trois secteurs peuvent être distingués à l’image de la partie en rive gauche avec un secteur en bord de chenal où la densité est modérée, un secteur central où la densité est plus importante et un secteur en bord de berge où la densité est plus faible.
Fig. 5 – Croquis des principales unités géomorphologiques et MNT des seuils à Portet-sur-Garonne (A), Beauzelle (B) et Grenade (C)
Fig. 5 – Diagram of the main geomorphological units and DEM of the knickpoints at Portet-sur-Garonne (A), Beauzelle (B) and Grenade (C)
1. Chenal ; 2. Affleurements rocheux ; 3. Bancs de galets ; 4. Bancs végétalisés ; 5. Terrasse holocène ; 6. Talus (h ≥ 4 m) ; 7. Sens d’écoulement
1. Channel; 2. Bedrock outcrops; 3. Gravel bars; 4. Vegetated bars; 5. Holocene terrace; 6. Embankment (h ≥ 4 m); 7. Flow direction
33La caractérisation du type de distribution des formes donne pour chaque site un indice de dispersion (ID) supérieur à 1, ce qui indique une répartition de type agrégative (tab. 2). Les résultats des tests du X2 et ZD indiquent des valeurs significativement supérieures aux valeurs seuils au risque α 0,05 pour Portet-sur-Garonne (X² = 5,76x103 > 21,03 ; 12 d.d.l, p < 0,05 et ZD = 4,54 > 1,96), Grenade (X² = 1,32 x 107 > 25,00 ; 15 d.d.l, p < 0,05 et ZD = 5,29 > 1,96) et Beauzelle (X² = 33,7 < 15,51 ; 8 d.d.l, p < 0,05 et ZD = 0,90 < 1,96). Cependant pour Beauzelle, ZD ne permet pas de rejeter l’hypothèse Ho d’après ce test. Les résultats de l’indice de Morisita (IM) vont dans le même sens avec des valeurs significativement supérieures à 1 pour l’ensemble des sites indiquant une distribution spatiale agrégative (tab. 2). Les valeurs du test du X2 confirment la distribution et permettent de rejeter l’hypothèse Ho. Cette dernière est donc rejetée à partir des résultats des indices pour les trois sites, confirmant la distribution agrégative des formes d’érosion et une liaison entre forme et hydro-dynamisme local.
Tab. 2 – Paramètres statistiques de l’Indice de Dispersion (ID) et de l’Indice de Morisita (IM) pour chaque site
Tab.2 – Statistical parameters of the Dispersion Index (ID) and the Morisita Index (IM) for each site
Dm : densité moyenne, s2 : écart-type, d.d.l. : degrés de liberté, X2 : test du Khi²
Dm: mean density, s2: standard deviation, d.d.l.: degrees of freedom, X2: Chi-square test.
34Les trous d’eau échantillonnés montrent une géométrie hétérogène avec une variation importante de leur dimension, tant sur le plan horizontal que vertical, pouvant aller de quelques centimètres à plusieurs dizaines de centimètres (tab. 3, fig. 6). Les valeurs moyennes sont respectivement de 15,53 cm et 15,71 cm à Portet-sur-Garonne et Grenade pour le diamètre ainsi que de 15,07 cm et 10,41 cm pour la profondeur. Les distributions des valeurs de diamètre et de profondeur présentent une dissymétrie positive avec des coefficients d’asymétrie compris entre 1,77 et 1,98. Les résultats du rapport D1/D2 indiquent que les trous d’eau ont généralement une ouverture elliptique avec des valeurs moyennes de 1,17 à Portet-sur-Garonne et 1,32 à Grenade, une valeur de 1 représentant une ouverture parfaitement circulaire. Le rapport D/Z présente des valeurs moyennes de 1,09 à Portet-sur-Garonne et 1,53 à Grenade. Ces résultats indiquent, que pour le second site, le développement latéral semble plus important que le développement vertical tandis que, pour le premier site, les dynamiques d’incision et d’érosion des parois sont comme similaires.
Fig. 6 – Fréquence d’occurrence des paramètres géométriques pour chaque site
Fig. 6 – Frequency of occurrence of geometrical parameters for each site
Z : profondeur, D : diamètre
Z: depth, D: diameter
Tab. 3 – Caractéristiques géométriques des trous d’eau pour chaque site d’étude
Tab. 3 – Geometrical characteristics of potholes for each study site
n : taille de l’échantillon, M : moyenne, Mo : mode, Me : médiane, SD : écart-type
n: sample size, M: mean, Mo: mode, Me: median, SD: standard deviation.
35L’analyse des données géométriques sur les deux sites met en évidence une forte corrélation entre le diamètre (D) et la profondeur (Z) (fig. 7), se traduisant par une relation de type linéaire pouvant être décrite par :
36ou
37Les résultats montrent dans les deux cas un bon ajustement avec des valeurs R2 de 0,66 à Portet-sur-Garonne et 0,75 à Grenade, confirmant l’observation graphique (fig. 7). Le coefficient ɛ correspond à l’amplitude du diamètre pour une profondeur donnée tandis que β correspond au taux d’augmentation du diamètre en fonction de la profondeur. Les valeurs β obtenues à partir des régressions sont pour chaque site inférieures à 1, ce qui signifie que la profondeur augmente plus vite que le diamètre. Ainsi, plus β est faible, plus l’allongement des formes dans leur composante verticale est important par rapport au développement de leur ouverture. C’est ce qui est traduit graphiquement (fig. 7) avec une pente de régression se démarquant d’une droite de régression où le diamètre et la profondeur augmenteraient au même rythme. Il est à noter que comparé à Portet-sur-Garonne (β = 0,88), la profondeur augmente plus rapidement que le diamètre à Grenade (β = 0,50).
Fig. 7 – Diagramme du diamètre (D) des trous d’eau en fonction de leur profondeur (Z)
Fig. 7 – Plots of potholes diameter (D) versus depths (Z)
38L’observation de trous d’eau en présence d’irrégularités de surface (fractures, micro-dépressions, écaillage, schistosité, etc.) met en évidence une assez bonne corrélation avec, à Portet-sur-Garonne, 75 % des trous d’eau qui s’inscrivent sur une surface irrégulière et 91 % à Grenade (fig. 8). Ces irrégularités peuvent ainsi représenter des faiblesses dans la roche, favorisant la formation de cavités. 42 % des trous d’eau du premier site et 43 % du second site présentent des fractures visibles.
Fig. 8 – Photographies montrant des trous d’eau formés à l’intersection de fissures (A-B) et sur surface irrégulière (A-B-C)
Fig. 8 – Photographs showing potholes developed at joint intersection (A-B) and on an irregular surface (A-B-C)
Les pointillés blancs représentent des fissures
The white dotted lines represent cracks
39Le relevé des sédiments dans les trous d’eau indique que 56 % à Portet-sur-Garonne et 83 % à Grenade piègent des sédiments. La granulométrie des particules les plus grosses se caractérise par un diamètre de grain médian de 58,2 mm pour le premier site et de 40,8 mm pour le second site (tab. 4). La relation entre le diamètre des plus grosses particules et la géométrie des trous d’eau ne permet pas de mettre en évidence une corrélation que cela soit pour le diamètre ou la profondeur (fig. 9A). À l’inverse, la répartition des effectifs montre que les formes avec des diamètres et des profondeurs inférieurs à environ 22 m présentent la plus grande quantité de sédiments piégés (fig. 9B). Il est à noter cependant que la distribution du nombre de particules en fonction de D/Z est plus importante pour un ratio d’environ 1,1 ± 0,2, et que les particules de plus grande taille s’observent pour un ratio d’environ 1,4 ± 0,5. De plus, la taille moyenne ne montre pas de tendance, cette dernière étant homogène sur l’ensemble de la gamme D/Z (fig. 9C).
Tab. 4 – Caractéristiques de la granulométrie des sédiments les plus gros
Tab. 4 – Characteristics of the largest grain size for each site
n : taille de l’échantillon, M : moyenne, Mo : mode, Me : médiane, SD : écart-type, D50 : diamètre médian
n: sample size, M : mean, Mo: mode, Me: median, SD: standard deviation, D50: median diameter
Fig. 9 – Caractéristiques sédimentologiques
Fig. 9 – Sedimentological characteristics
A : Graphique de la taille des particules en fonction du diamètre et de la profondeur. B : Graphique des effectifs des particules en fonction du diamètre et de la profondeur des trous d’eau pour les deux sites. C : Graphique de la répartition des effectifs, de la taille moyenne et maximale des particules en fonction de D/Z pour les deux sites
A: Graph of particle size versus diameter and depth for each site. B: Graph of particle size as a function of diameter and depth of potholes for both sites. C: Graph of the distribution of count, average and maximum particle size according to D/Z for both sites.
40L’inventaire des formes d’érosion a permis d’identifier 10 groupes et 52 sous-groupes de formes différentes. Cela traduit une diversité morphologique importante ainsi que des processus d’érosion et de contrôle impliqués dans l’évolution de ces formes.
41L’analyse de la distribution spatiale des formes d’érosion a mis en évidence pour chaque site une répartition agrégative, soit une distribution non homogène avec des densités contrastées selon les secteurs. Ce résultat suggère que certaines zones présentent une probabilité plus importante d’occurrence de formes avec une distribution pouvant être mise en relation avec la situation hydro-géomorphologique (relief, rugosité de surface, organisation spatiale des seuils, hauteur d’eau, vitesse d’écoulement, etc.). Globalement, il apparaît que les fortes densités de formes en creux s’observent préférentiellement à proximité du chenal d’écoulement principal et en tête de seuil, où les contraintes hydrauliques sont potentiellement les plus élevées (hauteur d’eau et vitesse d’écoulement importantes en crue), tandis que les secteurs les moins denses se localisent préférentiellement en bordure de berges et en queue de seuil, où les conditions hydrauliques sont moins favorables au développement de formes en raison de contraintes plus faibles. Les conditions locales d’écoulement (flux secondaires, rugosité de surface, etc.), en lien avec la micro-topographie, jouent également un rôle important pouvant expliquer la présence très localisée de « noyau » à forte concentration de formes au sein d’un espace faiblement dense. L’analyse par catégorie de formes montre également une localisation préférentielle selon le type, en lien avec le contexte hydro-géomorphologique et pouvant être associés à des processus d’érosion (dissolution, abrasion, etc.) différenciés.
42Les trous d’eau, se présentent généralement dans des zones de convergence d’écoulements et/ou dans des secteurs où la topographie présente des irrégularités marquées par une rugosité importante favorable à la formation de courants turbulents. La présence de sédiments dans ces cavités est également courante. Ainsi, il peut être suggéré que les trous d’eau sont entretenus principalement par l’action de flux tourbillonnants, leur morphologie favorisant le phénomène et l’action des sédiments par abrasion.
43Les sillons, coïncident majoritairement avec la présence d’un écoulement secondaire. C’est le cas notamment sur les sites de Portet-sur-Garonne et Grenade. Ces formes ont pour effet, à faible niveau d’eau, de canaliser l’écoulement avec une hauteur d’eau plus importante et une vitesse plus élevée par rapport à l’écoulement sur les surfaces adjacentes, accentuant l’effet érosif de l’écoulement. Par ailleurs, les sillons s’accompagnent très fréquemment de sédiments grossiers suggérant également une érosion par abrasion.
44Les flaques et dépressions, sont essentiellement localisées dans les secteurs au relief faiblement disséqué, où de l’eau stagnante est présente régulièrement. Par ailleurs, on note une absence de piégeage sédimentaire dans le cas des flaques. De forme étalée, à la surface peu rugueuse et de faible profondeur, les flaques semblent se développer essentiellement sous l’effet de la dissolution. Il peut en être de même pour les dépressions, dont on peut penser qu’elles correspondent à un stade d’évolution plus avancé de flaques.
45Les échancrures, sont principalement observées au niveau de rebords d’érosion qui correspondent à des ruptures de pentes abruptes, à l’interface entre deux niveaux pouvant être une alternance entre un faciès dure et tendre. Les échancrures se forment généralement dans le sens de l’écoulement et peuvent traduire un processus d’érosion régressive. Cela s’observe notamment sur le site de Beauzelle.
46Ainsi, les facteurs locaux semblent déterminants dans l’apparition, l’évolution et la morphologie des formes mineures d’érosion du lit.
47L’analyse morphométrique des trous d’eau cylindriques sur les sites de Portet-sur-Garonne et de Grenade, a mis en évidence une forte corrélation de type linéaire entre le diamètre et la profondeur dans les deux cas. Cette relation traduit une tendance à une augmentation plus rapide de la profondeur par rapport au diamètre d’ouverture. Néanmoins, le taux d’augmentation (β) du diamètre en fonction de la profondeur diffère sensiblement entre les sites (0,88 à Portet-sur-Garonne, 0,50 à Grenade). Cette augmentation plus rapide de la profondeur par rapport à la largeur sur le site de Grenade, peut s’expliquer par une présence moindre de calcaires indurés, les trous d’eau s’inscrivant principalement dans les faciès marno-sableux tendres. De plus, Springer et al. (2006), montrent que la perte de matériaux des trous d’eau lors de leur croissance est plus importante sur la paroi que sur le fond, lorsque β est supérieur à 0,5. Pour ces auteurs, cela indique le rôle principal des sédiments en suspension comme agent d’érosion agissant au-dessus du fond, par rapport aux galets retenus au fond des formes. Les valeurs de β dans la présente étude, suggèrent donc pour le site de Portet-sur-Garonne, (i) une érosion plus importante de la paroi par rapport au fond et (ii) un rôle prépondérant des sédiments en suspension dans l’érosion des trous d’eau par rapport au site de Grenade. Cependant, si les observations sur le terrain semblent aller dans ce sens (cf. infra), cela reste toutefois à confirmer quantitativement.
48La corrélation entre diamètre et profondeur a également été mis en évidence antérieurement par de multiples auteurs (Abbott et Pottratz, 1969 ; Lorenc, 1994 ; Springer et al., 2005, 2006 ; Wang et al., 2009 ; Álvarez-Vázquez et Uña-Álvarez, 2017 ; Dhali et Biswas, 2017 ; Ji et al., 2019). Le Tableau 5 reprend à titre de comparaison des résultats de travaux précédents par rapport à ceux obtenus dans cette étude pour des trous d’eau formés en rivière.
49Selon Springer et al. (2005, 2006), une corrélation importante entre le diamètre et la profondeur suggère que les trous d’eau peuvent être amenés à évoluer sous l’effet d’un processus de coalescence, un phénomène fréquemment observé sur les deux sites (fig. 10A). Pour ces auteurs, si les trous d’eau étaient amenés à s’effacer par l’incision du lit, le diamètre et la profondeur montreraient une dispersion. En effet, en l’absence d’incision la profondeur et le diamètre tendent à augmenter. A l’inverse, en cas d’incision du lit, la profondeur diminuerait tandis que le diamètre continuerait à augmenter pendant le processus d’effacement de la forme, amenant à une réduction de la relation diamètre vs. profondeur (fig. 10B).
Fig. 10 – Caractéristiques morphométriques de formes coalescentes
Fig. 10 – Morphometric characteristics of the coalescent shape forms
A : Forme coalescente résultant de la fusion de plusieurs trous d’eau (ligne en pointillés blancs) à la suite de la rupture de leurs parois (étoiles blanches) durant leur développement latéral. B : Schéma théorique de la relation entre le diamètre et la profondeur si les trous sont amenés à disparaitre par incision du chenal (modifié d’après Springer et al., 2006)
A: Coalescent shape resulting from the fusion of several potholes (dotted lines) following the rupture of their walls (white stars) during their lateral development. B: Theoritical diagram of the relationship between diameter and depth if the potholes were removed by channel incision (Modified from Springer et al., 2006).
50Les constatations visuelles suggèrent également une corrélation entre la présence de trous d’eau et celle d’irrégularités de surface (fractures, micro-dépressions, écaillage, schistosité, etc.). Cette corrélation est vérifiée par une correspondance de 75 à 91 % entre trous d’eau et irrégularités de surface. Dans le sud du Brésil, sur substrat basaltique, Lima et Binda (2015) font la même observation avec une corrélation entre formes d’érosion et irrégularités de surface. Springer et Wohl (2002) et Richardson et Carling (2005) notent le rôle fondamental des petites irrégularités dans la formation des formes et leur développement. La présence de fractures joue un rôle similaire dans la présence des trous d’eau, 42 % des formes montrant cette liaison en moyenne pour chacun des deux sites. L’importance des fractures dans l’initiation des écoulements turbulents et une érosion différentielle (Lima et Binda, 2015), est aussi mise en évidence sur des basaltes par Sengupta et Kale (2011) et Lima et Binda (2015), sur des granites par Lorenc (1994) et Ortega et al. (2014), sur des shales sableux par Wang et al. (2009) ou encore sur des quartzites et des gneiss par Springer et al. (2005, 2006).
Tab. 5 – Exemples de données rapportées pour les trous d’eau suivant une relation linéaire entre le diamètre et la profondeur selon Ji et al. (2018)
Tab. 5 – Examples of reported data for stream potholes following a linear trend between diameter and depth according to Ji et al. (2018)
D : diamètre, Z : profondeur, β : taux d’augmentation du diamètre en fonction de la profondeur, ɛ : l’amplitude du diamètre pour une profondeur donnée, R² : coefficient de détermination
D: diameter, Z : depth, β: diameter increase rate as a function of depth, ɛ: diameter magnitude for a given depth R²: coefficient of determination.
51Le piégeage des sédiments, s’effectue majoritairement dans les formes aux dimensions inférieures à environ 22 cm en largeur et en profondeur. Ceci va à l’encontre des observations faites par Ji et al. (2019) pour lesquels les sédiments sont majoritairement présents dans des trous d’eau plus larges et profonds. Ji et al. (2019), ainsi que Pelletier et al. (2015), font également le constat d’une relation linéaire entre la taille des plus gros sédiments et le diamètre des trous d’eau, traduisant un piégeage épisodique de sédiments dont la taille est comparable à la moitié du diamètre. Dans notre cas, les résultats ne permettent pas de mettre en évidence une telle relation entre la taille maximale des sédiments et le diamètre (R² = 0,24 à Portet-sur-Garonne ; R² = 0,46 à Grenade). Cela suggère, d’une part, que les sédiments les plus gros ne sont pas plus piégés lorsque les trous d’eau sont larges et profonds, et d’autre part, qu’ils ne sont pas un agent d’abrasion plus efficace. Sklar et Dietrich (2001) ont ainsi montré, à partir d’une expérimentation, que le taux d’érosion par rapport à la taille des sédiments présente une dépendance non linéaire, et que l’abrasion est plus efficace pour des particules de taille intermédiaire. En effet, les sédiments doivent être, d’une part, suffisamment petits pour être mis en mouvement, et d’autre part, suffisamment grands pour se déplacer en contact avec le fond.
52Le taux d’occurrence des sédiments dans les trous d’eau est nettement différent entre Portet-sur-Garonne (56 %) et Grenade (83 %). Cela peut suggérer qu’à Grenade l’érosion des trous d’eau par les sédiments grossiers pourrait être prépondérante, à l’inverse de Portet-sur-Garonne. Pour ce dernier site, l’action des sédiments fins en suspension pourrait jouer un rôle important dans le développement des formes en l’absence de charge grossière. La présence de charge grossière dans les trous d’eau ne serait donc pas une condition à leur développement. En effet, s’il a été mis en évidence par certains auteurs que les sédiments grossiers sont plus efficaces que les sédiments fins en suspension dans l’érosion d’un substrat rocheux (Sklar et Dietrich, 2001 ; Wohl, 2015), d’autres ont montré l’importance de l’action des sédiments en suspension dans le développement des trous d’eau (Tinkler, 1993 ; Zen et Prestegaard, 1994 ; Hancock et al., 1998 ; Whipple et al., 2000). Toutefois, ce constat est à nuancer car il s’agit de l’observation d’un état final. En effet, l’analyse de la mobilité de la charge grossière, à l’aide de PIT tags (Passive Integrated Transponder), a montré que même piégés dans des formes d’érosion (sillon, trou d’eau, etc.), le piégeage des galets pouvait être temporaire (Jantzi, 2018). À Portet-sur-Garonne, 38 % des particules injectées dans des formes ont été éjectées hors de la forme et 52 % à Grenade. Il est ainsi possible que le taux d’occurrence des sédiments dans les trous d’eau soit plus important en raison de conditions de décrue plus favorables au piégeage.
53Les résultats de l’étude des formes d’érosion sur substrat molassique en moyenne Garonne peuvent être synthétisés de la manière suivante :
54(i) Les formes d’érosion sur substrat molassique présentent une grande diversité morphologique (10 groupes et 52 sous-groupes différents). Toutefois, ce recensement n’est pas exhaustif en raison de la méthode d’échantillonnage. La répartition des formes montre que tous les types ne sont pas présents dans chaque site, à l’image des échancrures et des flûtes présentes uniquement à Beauzelle et à Portet-sur-Garonne respectivement, ou alors de manière marginale comme les flaques à Grenade et les buttes à Portet-sur-Garonne (fig. 3A). A l’inverse, les trous d’eau et les sillons représentent les formes les plus communément observées sur l’ensemble des sites.
55(ii) La distribution spatiale des formes d’érosion montre une répartition de type agrégative (fig. 4, tab. 2) suggérant que certains secteurs présentent une plus forte probabilité d’occurrence de forme par rapport à d’autres. Cette observation a été mise en relation avec les conditions hydrodynamiques stationnelles qui semble déterminantes dans l’apparition et l’organisation spatiale des formes. Cependant, cette relation devra être précisée par des travaux complémentaires sur les dynamiques d’écoulements dans ces zones lors des hautes eaux. La lithologie représente également un facteur important avec la présence d’irrégularités de surface (fractures, micro-dépressions liée à l’impact des sédiments, schistosité, etc.) qui affaiblissent la roche et favorisent la formation de formes.
56(iii) La formation des seuils dans les années 1960, à la suite de l’évacuation de la charge alluviale sous l’effet de l’incision, suggère des vitesses d’évolution rapides des formes même s’il n’est pas exclu que certaines puissent être héritées. De plus, en raison de différents processus érosifs (abrasion, dissolution, etc.) intervenant dans l’évolution des formes, des caractéristiques lithologiques de la roche (fractures, couche calcaire induré vs. couche marno-sableuse tendre, etc.) et des conditions hydrodynamiques, il est probable que les vitesses soient contrastées. Cela confirme les observations sur les dynamiques d’évolution des seuils qui sont hétérogènes.
57(iv) La morphométrie des trous d’eau tend à montrer que l’érosion des affleurements rocheux n’est pas uniquement due à un processus d’incision du lit. La relation entre le diamètre et la profondeur des trous d’eau (fig. 7) suggère que la fusion des formes entre elles (coalescence) est un facteur important de l’érosion des seuils (fig. 10A). L’observation de nombreuses formes issues d’une coalescence entre des formes de types différents (trous d’eau avec des sillons, etc.) conforte les résultats de l’analyse morphométrique.
58Si l’étude des formes d’érosion a permis d’apporter des éléments de compréhension dans la dynamique d’évolution des seuils rocheux en moyenne Garonne, certains points restent néanmoins à éclaircir concernant les facteurs d’érosion. Des études complémentaires devront ainsi s’attacher à (1) préciser l’influence des conditions d’écoulement, et (2) quantifier le rôle respectif entre charge de fond et charge en suspension.
* Auteur correspondant : Tel : +33 (0)5 61 50 36 26
hugo.jantzi@univ-tlse2.fr (H. Jantzi)