Nos travaux ont été développés dans le cadre du projet SEDILOIR inscrit dans la Plan Loire et financé par la région Pays de la Loire, le FEDER et l'Agence de l'eau Loire-Bretagne (n ° PRESAGE 38212). Nous remercions les relecteurs de l’article pour leurs remarques utiles, qui ont permis d’approfondir la réflexion et d’améliorer la présentation de nos résultats.
1La thématique du transport sédimentaire et des changements morphologiques des systèmes fluviaux est très largement étudiée, en particulier pour les systèmes fluviaux caractérisés par une forte énergie (puissance spécifique supérieure à 30 W/m²). C’est le cas des rivières de montagne (Piton et Recking, 2016) et des grands cours d’eau (Claude et al., 2012 ; Arnaud et al., 2015). La littérature est peu abondante quand il s’agit des rivières à faible énergie (Petit et al., 2015 ; Lespez et al., 2015 ; Melun, 2012). Cependant, depuis quelques années ces rivières intéressent les chercheurs et les gestionnaires des cours d'eau. Elles sont, comme tout autre cours d'eau, l'objet de réaménagements dans le cadre de la DCE (Directive Cadre sur l’Eau de 2000 ; ONEMA, 2013) qui prescrit le libre écoulement de l’eau afin d’assurer les continuités écologiques et sédimentaires.
2Les travaux récents mettent le plus souvent en évidence, outre la faiblesse générale de la dynamique fluviale, la grande hétérogénéité des configurations qui se révèle au fur et à mesure des recherches (Harnischmacher, 2007 ; Csiki et Rhoads, 2014 ; Rollet et al., 2014 ; Dufour et al., 2015 ; Depret et al., 2015 ; Kuo et al., 2015). Elles sont liées aux conditions géomorphologiques des cours d'eau ou tronçons de cours d'eau étudiés (tracé, pente du lit fluvial, caractéristique des sédiments disponibles pour le transport solide, etc.), à la diversité de l'impact des conditions hydrologiques sur les morphologies fluviales (rôle des crues, du débit à pleins bords), et également des situations hydrologiques intermédiaires au cours desquelles les sédiments qui se déposent seront évacués lors de la crue suivante, plus ou moins forte). Le développement de technologies adaptées performantes a permis de la révéler. Il s'agit des mesures par bathymétrie et par LiDAR terrestre 3D (Light Detection And Ranging) en particulier (Thoma et al., 2001 ; Nasermoaddeli et Pasche, 2008 ; Merwade, 2009 ; Melun, 2012 ; Foerst et Rüther, 2013 ; Hamel et al., 2013) ainsi que par profileur acoustique Doppler ADCP (acoustic Doppler current profiler) permettant d'évaluer la structure de l’écoulement (Le Coz et al., 2009 ; Dongsu et Muste, 2012 ; Claude, 2012).
3Par ailleurs, une grande partie des caractéristiques des cours d'eau de faible énergie est reliée aux aménagements des lits, souvent importants et anciens de plusieurs siècles (seuils transversaux, digues, chenalisation, etc.). Ils ont contribué à modifier l'écoulement en agissant sur la vitesse de l’eau, les volumes écoulés et leurs variations au cours de l'année. Ils ont également contribué à modifier les caractéristiques du transport sédimentaire en réduisant la mobilité de la charge la plus grossière et en donnant ainsi la faveur aux éléments fins tels que les limons et les sables (Talbot et Lapointe, 2002 ; James et Marcus, 2006 ; Gregory, 2006 ; Ikeda, 1984 ; Jaeggi, 1984 ; Hohensinner et al., 2014 ; Ollero, 2010). Le rôle des barrages sur les cours d’eau, largement étudié, est assez bien connu, du moins dans les grandes lignes (Williams et Wolman, 1984 ; Hadley et Emmett, 1998 ; Schmidt et Wilcock, 2008 ; Zahar et al., 2008 ; Csiki et Rhoads, 2010 ; Baker et al., 2011 ; Lobera et al., 2015 ; Rinaldi et al., 2013). On admet généralement que la stabilisation des lits fluviaux a permis que les paysages se modifient (fixation des lits fluviaux, modification de la végétation des berges, transformation de la fonction des rivières et des aménagements associés, etc.) (Germaine et Barraud, 2013). Ainsi, la situation actuelle des cours d'eau est-elle issue de la combinaison des conditions naturelles de l'érosion fluviale et de l'ajustement des cours d'eau à des modifications anthropiques dont les effets se superposent ou se juxtaposent aux conditions naturelles.
4Les travaux conduits sur le Loir permettent de préciser l’ampleur du rôle de ces transformations anthropiques dans le fonctionnement actuel de la rivière et de confirmer ou non les caractéristiques du transport sédimentaire qui concerne surtout des limons et des sables fins. Malgré la faible énergie du cours d’eau (de l’ordre de 5 W/m²), des secteurs fortement érodés se sont développés ; nos mesures visent à montrer quels sont les processus qui y participent. Le terrain d’étude prend en compte un secteur de Loir long de 60 km (de Poncé sur le Loir à La Flèche), qui sert à déterminer les conditions générales de la dynamique de la rivière (énergie fluviale et état des berges). Ce secteur comprend trois sites implantés à proximité de barrages (Ruillé, Bruant, et Cherré) où sont réalisées des mesures de l’érosion du lit fluvial. Les résultats obtenus devraient contribuer à la compréhension du fonctionnement des rivières de faible énergie souvent modifiées par les interventions anthropiques.
5Dans sa partie aval, le Loir est incisé de 60 à 80 m dans les plateaux sédimentaires de l’ouest du bassin de Paris (fig. 1). La rivière draine un fond de vallée large en moyenne de 2 000 m à l'intérieur duquel le lit majeur est dominé par une très basse terrasse culminant entre +1 à +3 m (fig. 2). Les matériaux du fond alluvial, épais de 4 à 5 m, ont été décrits par Piana et al. (2011) et Despriée et al. (2003). Ils sont composés de bas en haut par des cailloutis très grossiers surmontés par des graviers (diamètre 5 cm) disposés en niveaux sub-horizontaux, épais d’une cinquantaine de centimètres et séparés par des niveaux de sables et de graviers. Au-dessus, on observe 1,50 à 2 m de sédiments limono-argileux récents.
Fig. 1 – Le bassin versant du Loir.
Fig. 1 – Loir’s watershed.
A. Situation. B. Grands ensembles lithologiques. 1. Limite du bassin versant du Loir ; 2. Loir et principaux affluents ; 3. Ville principale ; 4. Limite de département ; 5. Faille ; 6. Faille avec compartiment affaissé ; 7. Sable, argile et graviers de fond de vallée ; 8. Calcaire de Beauce ; 9. Argiles à silex des collines du Perche ; 10. Craie du Turonien et du Sénonien ; 11. Sables et grès du Cénomanien ; 12. Calcaire en affleurement du Jurassique ; 13. Tronçon fluvial étudié.
A. Situation. B. Major lithological units. 1. Catchment area; 2. Loir river and main tributaries; 3. Main cities; 4. Administrative limit (« département »); 5. Fault; 6. Fault and downlifted zone; 7. Sand, clay, and gravel of the bottom valley; 8. Lacustrine limestone of beauce; 9. Flint clay of the perche hills; 10. Turonian and senonian chalks; 11. Cenomanian sands and sandstones; 12. Jurassic outcrops limestone; 13. Fluvial section studied.
Fig. 2 – Partie de la vallée du Loir comprenant les sites de mesure.
Fig. 2 – Part of the Loir’s valley where measurement sites are localized.
1. Vestige de seuil hydraulique ; 2. Seuil hydraulique ruiné ; 3. Seuil hydraulique fonctionnel ; 4. Point coté ; 5. Principales communes ; 6. Le Loir ; 7. Affluent ; 8. Bras de dérivation ; 9. Ancien chenal ; 10. Canal exploitant un ancien chenal ; 11. Site de mesures ; 12. Principaux sapements ; 13. Lit majeur ; 14. Basse terrasse ; 15. Limite du fond de vallée.
1. Remnants of ancient sill; 2. Removed sill; 3. Sill; 4. Altitude; 5. Main villages; 6. Loir River; 7. Tributary; 8. Spillway; 9. Paleochannel; 10. Canal located in a paleochannel; 11. Survey sites; 12. Main lateral erosion sites; 13. Major bed/floodplain; 14. Terrace; 15. Floodplain (including terraces).
6Le Loir, profond de 2,5 m et large de 40 m environ, décrit un tracé sinueux (coefficient de sinuosité de 1,3) dont les boucles les plus serrées (rayon de courbure de l’ordre de 150 m) sont héritées de tracés anciens (Piana et al., 2011). La faible pente du lit fluvial est comprise entre 0,2 ‰ et 0,4 ‰. Le lit mineur est délimité par des berges hautes de 0,5 à 1,5 m, fixées par une végétation arbustive et buissonnante, à l’exception de quelques endroits où elles sont soumises à l'érosion (fig. 3-4).
Fig. 3 – Berge limono-sableuse à Cherré (juin 2015).
Fig. 3 – Loamy sand riverbank at Cherré (June 2015).
Fig. 4 – Matériaux grossiers en pied de berge à Cherré (octobre 2014).
Fig. 4 – Coarse materials at the riverbank base visible at Cherré (October 2014).
7Le Tableau 1 donne les caractéristiques hydrologiques principales du Loir. Les valeurs atteintes sont le fait (i) des apports pluviométriques atteignant 650 mm en moyenne par an, et (ii) de la lithologie du substratum (calcaire et grès), qui soutiennent les écoulements en période de basses eaux. Les crues sont caractéristiques des régions océaniques tempérées de plaine et plateau, avec une montée et une descente de crue en 4 à 7 jours (Larue, 1989). Des aménagements au fil de l’eau ont été réalisés pour permettre différents usages. De discrètes marques semblent liées aux travaux effectués pour le flottage du bois (en aval de la Flèche, entre la Forêt de Bercé et la Loire) et la navigation arrêtée dès le début du XIXe siècle (secteurs de lit rectifiés, anciens seuils noyés, chenal dédoublé, réalisations anciennes encore suggérées par la toponymie, port, gué, etc.). Les aménagements les plus visibles sont les travaux de canalisation du lit dans la traversée des bourgades (Le Lude, Château du Loir), et les barrages. Ces derniers, d’une hauteur comprise entre de 1,5 et 1,8 m, comprennent une chaussée souvent oblique et des clapets mobiles (Corbonnois et Techkpo, 2013). Le long du tronçon d'étude, ils servent à la production d'hydroélectricité, au maintien de zones humides et surtout à la gestion des niveaux d'eau. Plusieurs sont en ruine (fig. 2). Le Loir est concerné par la directive cadre européenne sur l'eau (2000) qui prescrit le rétablissement des continuités écologiques et sédimentaires qui devra intervenir avant 2027 ; ici, l’objectif principal est de permettre que s’améliore la qualité de l’eau et des milieux aquatiques (SAGE Loir, 2016).
Tab. 1 – Débits moyens annuels et mensuels sur le secteur d’étude (données issues de la Banque Hydro consultée le 22 mai 2016).
Tab. 1 – Annual and monthly mean rates measured in the studied area (Data from Hydro database accessed May 22nd 2016).
8Les trois sites (Ruillé, Bruant, et Cherré) retenus pour l’étude détaillée des conditions de l’érosion fluviale concernent un tronçon de lit sinueux de 100 à 300 m de long, caractérisé par des berges en partie érodées (fig. 2). Les caractéristiques géomorphologiques des trois sites ont été cartographiées à partir du terrain et de photographies aériennes. La Figure 5 représente celles du site de Bruant (SEDILOIR, 2016).
Fig. 5 – Caractéristiques géomorphologiques du Loir au site de Bruant.
Fig. 5 – Geomorphological characteristics of the Loir River at Bruant.
A. Fond de vallée. 1. Berge sapée ; 2. Berge dégradée ; 3. Berge stable végétalisée ; 4. Lit mineur ; 5. Ancien passage du loir ; 6. Très basse terrasse ; 7. Basse terrasse ; 8. Limite du fond de vallée ; 9. Barrage. B. Lit du loir. 1. Ancien déversoir ; 2. Berge protégée (maçonnée, palplanche…) ; 3. Berge stable, végétalisée ; 4. Berge boisée ; 5. Arbre isolé ; 6. Berge dégradée ; 7. Berge sapée ; 8. Zone de levé bathymétrique ; 9. Secteur de berge « scannée » ; 10. Profil sédimentaire ; 11. Banc vif d’alluvions ; 12. Banc d’alluvions végétalisé ; 13. Chenal de crue ; 14. Ligne principale de courant ; 15. Seuil et niveau graveleux ; 16. Passées graveleuses dans les matériaux des berges ; 17. Rampe de mise à l’eau ; 18. Station de pompage.
A. Valley floor. 1. Eroded bank; 2. Degraded bank; 3. Vegetated stable bank; 4. Active channel/river bed; 5. Paleochannels; 6. Lower terrace; 7. Terrace; 8. Floodplain (including terraces); 9. Dam; b. Loir’s channel. 1. Ancient spillway; 2. Ripraped bank; 3. Vegetated stable bank; 4. Vegetated stable bank (dense riparian tree cover); 5. Isolated tree; 6. Degraded bank; 7. Eroded bank; 8. Bathymetrical survey; 9. Terrestrial lidar survey; 10. Sediment architecture and grain size survey; 11. Bar; 12. Vegetated bar; 13. Flood channel; 14. Main stream; 15. Threshold and gravel level; 16. Coarse sediment beds (gravels) in fine bank sediments; 17. Boat access area ; 18. Pumping station.
9À l’exception de 106 profils transversaux du lit du Loir levés en 2004 tous les 500 m entre la confluence de la Braye près de Poncé-sur-le Loir et l’amont de la Flèche (fig. 2) (EPL, 2004 ; HYDRATEC, 2008), toutes les données utilisées proviennent de mesures effectuées dans les trois sites entre l’automne 2013 et mai 2015 (fig. 6, tab. 2), pour différentes conditions hydrologiques afin d’en saisir les impacts sur le façonnement des formes fluviales. Des traitements spécifiques ont permis de restituer les actions ténues de l'érosion enregistrées par les berges et le fond du lit. Les données des trois sites ont été traitées selon des méthodes similaires. À titre d’exemple, la Figure 6 représente pour le site de Bruant l'hydrogramme journalier couvrant la période du suivi, ainsi que les dates et les types de mesures réalisées.
Fig. 6 – Hydrogramme de la station de Flée (juin 2013 à mai 2015) et mesures effectuées sur le site de Bruant.
Fig. 6 – Hydrograph acquired at Flée gauging station (from June 2013 to May 2015) and measurements carried out at Bruant.
Tab. 2 – Mesures réalisées sur les trois sites d’étude.
Tab. 2 – Measurements carried out on the three sites examined.
10Les 106 profils transversaux ont permis de calculer (i) le débit à pleins bords du Loir à partir des données de la station de Flée ; il est de l'ordre de 140 m3/s (selon la loi de Gumbel) et varie longitudinalement en fonction des configurations morphologiques du lit, (ii) les puissances spécifiques (Petit et al., 2005 ; Schmitt et al., 2001). Pour toute la section fluviale concernée, ces dernières sont généralement faibles, atteignant le plus fréquemment quelques watts par m². Les plus élevées (une vingtaine de valeurs) dépassent 15 W/m² mais restent inférieures à 30 W/m². Parfois, la pente du lit trop faible pour être évaluée avec une précision suffisante, conduit à obtenir des valeurs de puissances spécifiques proches de zéro.
11L’auscultation des berges réalisée grâce à un Scanner Laser 3D (modèle Leica Scan Station C10) permet un relevé précis de la topographie des berges sur toute leur hauteur (0,7 à 2 m). L’érosion est évaluée à partir de la comparaison des relevés effectués à chaque campagne de mesures. Les données collectées sur le terrain se présentent sous la forme de nuages de points denses (fonction de la fréquence d’échantillonnage spatial : ici au minimum 1 point par 5 cm²). Une classification hiérarchique descendante permet d’isoler les points représentant les matériaux qui composent les berges de ceux correspondant à la végétation rivulaire (sur 0,2 à 0,3 m sous la surface du sol). La méthodologie appliquée est identique à celle proposée par Brodu et Lague (2013), soit la sélection d’échantillons caractéristiques de l’environnement ausculté et la création d’une classification multi-échelles. Seuls les points classifiés ayant un indice de confiance supérieur ou égal à 90 % ont été sauvegardés (des mesures supplémentaires ont été faites à Bruant et à Cherré en mars et avril 2015 pour le traitement des parties supérieures de berges avant que la végétation annuelle se développe). Pour évaluer le recul des berges d'une campagne à l'autre, les nuages de points ont été consolidés et géoréférencés dans le système géodésique français RGF93.
12Le calcul des distances de recul des berges et des volumes érodés est fait à partir de deux algorithmes dénommés respectivement « Distance Computation – According to a Given Direction » (DC-AGD) et « Volume Computation – According to a Given Direction » (VC-AGD), et développés au sein du Laboratoire Géomatique et Foncier du Conservatoire national des arts et métiers (Fadili, 2015). S’agissant du calcul du recul de berge, l’algorithme DC-AGD s’inspire de l’algorithme M3C2 développé par Lague et al. (2013). La différence principale réside dans la prise en compte de la direction de l’évolution du recul de la berge au cours du temps (perpendiculaire ou oblique) et non selon la normale comme le proposent les deux auteurs. Ceci permet d'obtenir une meilleure adéquation avec les dynamiques observées. Une distance de recul point par point est alors calculée. Pour le Calcul du volume érodé, l’algorithme VC-AGD permet des cubatures comprises entre deux surfaces et selon une direction donnée. Le principe est de projeter une surface maillée, qui correspond ici à la campagne la plus ancienne, sur le nuage de point de la campagne la plus récente. Afin d’évaluer les marges d’incertitude liées au calcul, une simulation de nuages fictifs a été réalisée là où les distances de recul et les volumes érodés sont connus. L’incertitude est inférieure à 0,005 m pour les distances ; elle est de 0,01 m3 pour les volumes calculés. À titre d’exemple, la Figure 7 représente pour le site de Bruant, la distribution statistique d des distances de recul entre les différentes campagnes calculées à l’aide de l’algorithme DC-AGD. Ainsi le 1er et le 9ème décile, la médiane, ainsi que le 1er et 3ème quartile, sont représentés afin de réaliser une comparaison diachronique des résultats obtenus.
Fig. 7 – Distribution statistique des distances de recul calculées à l'aide de l'algorithme DC-AGD selon la direction Est-Ouest sur le site de Bruant.
Fig. 7 – Statistical distribution of recession distances from east to west on the site of Bruant, calculated using DC-AGD algorithm.
13Le logiciel Hypack a été utilisé pour l’acquisition des données sur le terrain (Hypack Survey). Le géoréférencement des levés bathymétriques est réalisé par GPS en mode temps réel RTK (Real-Time Kinematic) utilisant le réseau de stations permanentes Orphéon. L’incertitude sur la hauteur est estimée à 0,1 mètre. Par ailleurs, du fait des conditions du terrain (pied de berge difficile à mesurer et végétation en fond de lit en période estivale) les mesures ont été contrôlées par des mesures manuelles et mono-faisceau (sondeur bathymétrique Tritech PA500 kHz, gamme de profondeurs : 0,1 à 10 m). Différents Modèles Numériques de Terrain, avec un pas de 0,5 mètre ont été produits à l’aide du logiciel Surfer 3D. Les résultats du traitement ont permis de repérer de manière satisfaisante la topographie du fond du lit (fig. 8), puis de mettre en évidence son évolution et de calculer des bilans sédimentaires pour la période de temps considérée.
Fig. 8 – Levé bathymétrique à Bruant le 28 janvier 2014.
Fig. 8 – Bathymetric survey conducted at Bruant on January 28th 2014.
14Pour ces mesures, nous avons mis en œuvre un ADCP M9 (SonTek) depuis une embarcation légère de type Neveu Secu 12. Le suivi planimétrique du profil théorique a été réalisé par le pilote du bateau grâce au logiciel Hypack 2013A et par l’asservissement d’un DGPS Magelan Proflex 500 en mode RTK équipé du système différentiel (l’incertitude altimétrique est estimée à 10 cm). Le traitement des données (5 profils par site, soit au total 15 profils de mesures) a permis de déterminer dans les sections mouillées, les champs de vitesses du courant. Nous utiliserons les résultats obtenus durant la campagne de janvier 2014 alors que le débit du Loir était proche du débit à pleins bords (fig. 9). Les codes MatLab développées par Claude (2012) pour exploiter ces informations ont été mobilisés dans ce travail.
Fig. 9 – Exemple de représentation des vitesses longitudinales en m/s issues du traitement réalisé à l’aide du développement MatLab de Claude (2012 ; profil 9 mesuré à Bruant)
Fig. 9 – Cross section of the Loir showing the longitudinal stream velocities in m/s provided by Claude’s (2012) MatLab code (profile n° 9 acquired at Bruant).
15Les échantillons de sédiments ont été prélevés manuellement. Dans les berges dénudées (fig. 3-4), les prélèvements ont été faits tous les 20 à 30 cm de la surface du sol jusque sous le niveau de l’eau. Dans le fond du lit, les matériaux ont été échantillonnés au droit des 15 sections de mesure ADCP, réparties dans les 3 sites, à raison de 3 à 5 échantillons par section (fig. 10). La granulométrie de la fraction grossière (supérieure à 2 mm) a été obtenue par tamisage humide sur une colonne de tamis à mailles comprises entre 2 à 50 mm (NF ISO 3310-1 et 3310-2). La fraction fine (inférieure à 2 mm) a été analysée par granulométrie laser (Mastersizer 3000 - ISO 13320).
Fig. 10 – Localisation des prises d’échantillons de sédiments (berges et fond du lit). Site de Bruant.
Fig. 10 – Location of sediment sampling sites (riverbanks and riverbed bottom) at Bruant.
16Dix-huit profils en travers ont été mesurés sur les secteurs de Bruant et Cherré dans le but de réaliser une modélisation hydraulique 1D en régime permanent non uniforme à l’aide du logiciel HEC-Ras 4.1. Le levé de la partie terrestre a été effectué par GPS en mode NRTK (Near Real-Time Kinematic) et celui de la partie aquatique par levés réalisés avec un mono-faisceau Tritech PA500 couplé à un GPS en mode NRTK. L’incertitude des mesures de ce levé est inférieure à 0,05 m dans les trois composantes. Cette modélisation sert à mettre en évidence, à l’échelle du site, l’évolution des paramètres qui influencent la dynamique de la rivière et à établir des relations avec les résultats de nos mesures.
17Dans la partie inférieure du Loir, la densité des mesures apparaît insuffisante pour que des corrélations significatives puissent être établies entre la puissance spécifique calculée pour chaque profil transversal du lit et les conditions du milieu (surface des sections mouillées, pente, aménagements, débit), en particulier la présence de barrages (SEDILOIR, 2016). Pour ceux-ci, aucune influence forte n’est démontrée comme cela a également été mis en évidence par les travaux de Csiki et Rhoads (2014) pour des rivières américaines. Tout juste peut-on établir d’après la Figure 11, qui prend en compte la dispersion des puissances spécifiques, des tendances ; les barrages successifs contribuent à la faible énergie dans les secteurs 2 et 4. Dans les autres, les liens ne peuvent pas être établis. Les puissances spécifiques, toujours inférieures à 30 W/m², expliquent que l’érosion fluviale généralement soit peu active. D’après les prospections de terrain, 85 % de la longueur du lit du Loir sont stables et présentent des berges végétalisées. Les 15 % restant comprennent des marques d’érosion, en particulier des sapements de berges (fig. 2, 12-13) situés dans les boucles les plus serrées du tracé ou encore associés localement à des désordres provoqués par les interventions anthropiques. Nos sites de mesures sont localisés dans les secteurs 3, 4, et 5 en des endroits où la puissance spécifique est de l’ordre de 5 W/m², mais où l’érosion fluviale est bien visible.
Fig. 11 – Puissances spécifiques d du Loir, de l'aval (La Flèche) à l'amont (Poncé-sur-le-Loir).
Fig. 11 – Variation of the unit stream power of the Loir River from the Braye confluence to upstream of La Flèche.
Fig. 12 – État des berges (site de Bruant).
Fig. 12 – Riverbank situation (Bruant).
Fig. 13 – Barrage de Cherré.
Fig. 13 – The dam of Cherré.
18À l’échelle des trois sites d’étude et pour la période d'observation, les résultats des mesures montrent que les berges reculent très irrégulièrement. Sur le site de Bruant (47°39’ N ; 0°20’E), quatre campagnes de LIDAR terrestre ont été réalisées en juin 2013, juillet 2014, octobre 2014 et mars 2015. Au cours de l'ensemble de la période d'observation, la station de Flée a enregistré neuf pics de crue entre juin 2013 et juillet 2014 (débits supérieurs à 80 m3/s pour un maximum de 143 m3/s) et trois entre octobre 2014 et mars 2015 (débits entre 45 m3/s et 68,8 m3/s). Le maximum enregistré sur la période correspond au débit de crue de fréquence 0,5 (fig. 14). L'érosion a fait reculer certaines portions de berge sur une distance de plus de 1,50 m (fig. 15). Mais ce recul reste inférieur à 0,15 m pour 75 % des points entre juillet 2014 et octobre 2014 (saison estivale) et inférieur à 0,30 m entre octobre 2014 et mars 2015 (saison hivernale). Les reculs plus importants, mesurés pour l’ensemble de la période allant de juin 2013 à juillet 2014, sont liés à la plus longue durée de cette période (13 mois) mais aussi aux conditions hydrologiques plus favorables à l'érosion (3ème quartile égal à 0,99 m). Ces reculs, très irréguliers (fig. 15), indiquent que l’essentiel se fait par l’approfondissement d’encoches d’érosion puis par l’ablation des « promontoires » qui les séparent (fig. 16). Ces actions surviennent principalement en période hivernale comme le montre la Figure 14 qui représente pour le site de Bruant, les volumes de sédiments érodés en m3 par m² de surface de berge scannée, en regard de la variation de débits.
Fig. 14 – Hydrogramme des débits moyens journaliers à la station de Flée et bilan sédimentaire des berges à Bruant (calculé à l'aide de l'algorithme VC-AGD).
Fig. 14 – Hydrograms of mean daily discharges at Flée station and sedimentary budget of riverbanks calculated using DC-AGD algorithm at Bruant.
Fig. 15 – Distance de recul point par point de la berge entre Juin 2013 et Juillet 2014 sur le site de Bruant calculée à l’aide de l’algorithme DC-AGD selon la direction Est-Ouest.
Fig. 15 – Point to point recession distance of the riverbank between June 2013 and July 2014 at Bruant calculated by means of DC-AGD algorithm in an east-west direction.
Fig. 16 – Distance de recul moyennée par tranche de 0,5 mètre de long sur le site de Bruant à l'aide de l'algorithme DC-AGD.
Fig. 16 – Recession distances acquired from east to west on Bruant and averaged in slices of 0.5 m in width, calculated using DC-AGD algorithm.
19Les levés bathymétriques montrent que le lit du Loir a une forme irrégulière, en adéquation d'abord avec le tracé de la rivière (Bravard et Petit, 1997 ; Malavoi et Bravard, 2010). Il comprend (i) de rares bancs d'alluvions fines, souvent végétalisés ou surmontés par une mince tranche d'eau ; ils sont allongés contre la berge ou présents dans les zones où la vitesse du courant est très faible, (ii) des hauts fonds en aval des points d'inflexion des sinuosités, formés de graviers et petits galets, le plus souvent immergés, (iii) des mouilles atteignant 1 à 3 m de profondeur, creusées dans le fond du chenal, et localisées en rive concave d’une sinuosité ou bien dans l'axe du chenal, (iv) des traces d'aménagements divers ayant modifié la forme du lit surtout aux abords des barrages (fond plat à l'amont de barrages, surcreusement en aval).
20Nos mesures montrent que ces formes se sont peu modifiées pendant la période d’observation, de l'automne 2013 au printemps 2015. Des variations de la morphologie du fond, de l'ordre du décimètre, sont liées aux dépôts et reprises de sédiments qui affectent indifféremment toutes les unités. La Figure 17 indique pour le site de Bruant les modifications intervenues entre janvier et octobre 2014. Celles-ci concernent surtout les bancs de vase mobile, situés sous l'eau. Ailleurs il n’y a pas de lien avec la topographie du fond du lit. Il apparaît ainsi que ce sont surtout les matériaux fins qui sont déplacés, en faibles volumes. À Bruant (budget sédimentaire : -0,32.10-5 m³m²j-1 entre le 10 octobre 2013 et le 28 janvier 2014, et +0,38.10-5 m3m²j-1 entre le 28 janvier 2014 et le 20 octobre 2014), on estime qu’ils sont mis en mouvement à partir d’un débit de 70 m3/s survenu une dizaine de fois pendant la période d'observation. La situation est la même dans les autres sites.
Fig. 17 – Écart d'altitudes du fond du lit entre le 28 janvier et le 20 octobre 2014 (site de Bruant).
Fig. 17 – Differences of river bottom elevation measured between January 28th and October 20th 2014 at Bruant.
21Les résultats montrent que les matériaux issus du recul des berges sont majoritairement évacués. Les dépôts dans le fond du lit sont quasi-négligeables (tab. 3). À Ruillé, où le site de mesure est directement en aval d'un ancien barrage qui dérive vers le canal d’un moulin une partie du débit en crue (ce barrage est percé par plusieurs brèches), les berges livrent le moins de matériaux. Sur le fond, les sédiments proviennent principalement de vase déposée là où la vitesse du courant est très faible et facile à remobiliser lors des crues. L’érosion est plus importante dans les deux autres sites. À Bruant, les berges livrent autant de matériaux que le chenal en évacue. À Cherré, les apports de berges sont moins efficacement entraînés vers l’aval. En ces deux sites, la puissance spécifique du Loir, calculée pour un débit à pleins bords à partir de profils transversaux serrés mesurés en mars 2015, est faible ; ce résultat est en cohérence avec ceux obtenus ailleurs le long du Loir (fig. 11). Elle atteint au maximum 12,5 W/m² à Bruant (fig. 18) et 20 W/m² à Cherré. Cette énergie suffit à évacuer les fines fournies par les berges. Elle est insuffisante pour prendre en charge la fraction grossière qui par ailleurs a peu bougé pendant la période de mesure, comme le montrent les résultats de la bathymétrie.
Tab. 3 – Résultats des mesures de l'érosion.
Tab. 3 – Results of erosion measurements.
Fig. 18 – Puissances spécifiques calculées sur les 18 profils réalisés sur le site de Bruant.
Fig. 18 – Unit stream power computed on the 18 cross sections carried out at Bruant.
22Nos résultats issus de l'interprétation des analyses granulométriques montrent que la fraction grossière présente dans le fond du lit du Loir a une taille supérieure à celle des graviers et galets qui affleurent dans les berges (échantillons BR 02, 03, 04 prélevés dans le profil B1). Les galets (> 50 mm) sont les plus abondants à Bruant où leur part dépasse 40 % du poids des échantillons dans la moitié des prélèvements (contre 1/5 à Cherré et 1/6 à Ruillé). Le reste des matériaux est majoritairement constitué par des graviers et granules de 20 à 50 mm et de 5 à 20 mm de diamètre. Les fractions inférieures à 5 mm sont peu abondantes sauf en des points particuliers, comme les pieds de berge alimentés par des éboulements de paquets de limons sableux et des bancs d'alluvions vaseuses (échantillons Bru 11 et 13 sur le site de Bruant) (fig. 19). Ainsi, si la faible capacité de transport du Loir limite le transfert amont-aval des sédiments, les profils ADCP mesurés à Bruant en janvier 2014 pour un débit proche du débit à pleins bords (fig. 9-10, 19), montrent que les champs de vitesses, assez variés d'une section à l'autre et à l'intérieur d'une section mouillée, semblent suffisamment contrastés pour générer des turbulences capables de mouvoir des éléments grossiers dans certaines parties du lit mineur. L’analyse des résultats de nos campagnes de mesures permet d’établir des liens entre diamètres de particules et vitesses du courant (tab. 4, fig. 20). Ainsi est-il possible de proposer plusieurs hypothèses relatives à la dynamique des sédiments grossiers dans le fond du lit du Loir. Il comprend (i) des galets quasi-immobiles en pied de berges concaves (diamètre supérieur à 50 mm) ; ils constituent ou dépassent 70 % du poids des échantillons (Rui 14, Bru 6, 7, 15, Che 8, 9, 10). Ils sont hérités (recoupement par le Loir de niveaux grossiers remplissant le fond alluvial) et forment un pavage qui bloque l'approfondissement actuel du lit, (ii) des graviers de 20 à 50 mm de diamètre (fraction prédominante dans la plupart des échantillons) et dans une moindre mesure, de granules de 5 à 20 mm de diamètre. Mobilisables par les courants les plus forts, ils forment une armure qui protège de l'incision sauf en conditions exceptionnelles, (iii) des matériaux fins en transit dans le chenal et dont la mobilité au cours des saisons hivernales a été montrée à partir des mesures de bathymétrie.
Fig. 19 – Principales fractions granulométriques des échantillons prélevés sur le site de Bruant exprimé en % (Bru. Matériaux du fond du lit ; Br. Matériaux des berges).
Fig. 19 – Main grain size classes (in %) of the samples taken at Bruant (Bru. Riverbed materials ; Br. Riverbank materials).
Tab. 4 – Dimensions des particules et classes de vitesses du courant à l’intérieur d’une section mouillée du Loir (janvier 2014 pour un débit proche du débit à pleins bords).
Tab. 4 – Particle sizes and stream velocity classes on a wetted cross section of the Loir River (on January 2014 given a flow rate close to the bankfull discharge).
Fig. 20 – Champ de vitesse du Loir. Exemple du profil N° 7 à Bruant (débit proche du débit à pleins bords).
Fig. 20 – Velocity field of the Loir River. Example from the profile #7 acquired at Bruant.
23L’hypothèse généralement admise est que les aménagements des cours d'eau comme le Loir ont contribué à la modification d’un paysage fluvial dont on ne connaît pas l’état l’antérieur du fait de l’ancienneté des installations. Ce paysage fait référence davantage à la végétation des berges et aux usages du cours d’eau qu’à l’érosion fluviale. Les ouvrages hydrauliques ont fixé un lit naturellement peu mobile compte tenu de la faiblesse de la pente et de sa faible capacité à mobiliser une charge sédimentaire abondante. La puissance spécifique moyenne, calculée à partir des 106 profils transversaux mesurés sur le Loir sur une longueur de 60 km, est de 11 W/m² ; elle est donc proche des 12 W/m² mesurée par Depret et al. (2015) sur le Cher. Dans 70 % des cas, cette puissance est inférieure à 5 W/m² et dans 10 % des cas seulement, supérieure à 20 W/m². Ces caractéristiques sont en cohérence avec l’état du lit dont les berges très peu mobiles sont végétalisées sur 85 % du linéaire de la rivière. Ailleurs, la stabilité du lit est ou a été rompue du fait d’aménagements parfois très anciens ou en lien avec le tracé du Loir (rives concaves des boucles les plus serrées). Le processus principal de l’érosion des berges là où elle s’est développée, est le sapement. Il agit sur les matériaux cohésifs souvent disloqués par la météorisation (variations d’humidité et de température). Les éboulements de nature gravitaire observés sur d’autres cours d’eau sont peu fréquents (Vandenberghe et al., 2012), et n’interviennent que lorsque des lentilles de graviers affleurent en pied de berge (Nasermoaddelli et Pasche, 2008).
24Sur les sites de Bruant et de Cherré, une érosion notable a entraîné le recul des berges d’une dizaine de centimètres en moyenne aux endroits des mesures. Le recul est plus important dans les encoches (jusqu’à 1,50 m). Les sédiments évacués par la rivière sont des fines (fraction inférieure à 5 mm). Les processus interviennent pour une très large gamme de débits et ne nécessitent pas de fortes crues. Comme l’ont montré d’autres recherches dont celles de Nardi et al. (2013), de Hooke (2015) ou de Palmer et al. (2014), il n’y a pas de corrélation directe entre les débits, en particulier la magnitude des crues, et l’érosion des berges pour des rivières de faible énergie. Le déblaiement des sédiments qui conduit au recul des berges a besoin d’une phase antérieure de préparation des matériaux lorsque le niveau limnimétrique est bas (période estivale ou entre deux crues), avant leur évacuation lors d’évènements hydrologiques favorables (remontée et/ou variation du niveau d’eau).
25Sur le Loir, le déblaiement ne concerne que la fraction fine inférieure à 5 mm de diamètre. Les mesures réalisées sur le site de Bruant (fig. 16), montrent que les modifications du trait de berge sont surtout visibles après la période des hautes eaux. Au cours de la période de basses eaux, les matériaux affaiblis par des fentes de dessiccation s'éboulent en pied de berge d'où ils seront repris au moment des crues. Mais ces éboulements ne concernent que de faibles volumes de sédiments.
26Ces modifications du lit interviennent malgré la faible énergie du cours d’eau. À Bruant et à Cherré, la puissance spécifique moyenne est respectivement de 6 et 4 W/m², Il faut donc prendre en considération d’autres moteurs. Dans le détail, la forme du lit du Loir est irrégulière ; le fond comprend des mouilles profondes (jusqu'à 2,5 à 3 m) et les berges sont échancrées par des encoches qui accroissent leur rugosité et contribuent à accentuer les forces de cisaillement. Ces caractéristiques sont propices à la formation d'un écoulement plus turbulent qui doit amplifier l'érosion. La turbulence à l’intérieur du chenal a été étudiée en laboratoire sur canal expérimental (Yue et al., 2015) et dans le cadre d’approches hydro-écologiques pour la définition des habitats (Kronvang et al., 2012). Dans nos sites, elle permet que les sédiments fins prélevés dans le lit soient évacués. Sont ainsi mis à nu des matériaux délaissés par une érosion insuffisamment puissante ; les conditions hydro-climatologiques actuelles ne permettent pas la reprise de ces matériaux et l'incision du lit. On peut émettre l’hypothèse que leur mobilisation est très peu fréquente. La précision actuelle de nos mesures n’est pas suffisante pour évaluer les paramètres hydrauliques qui pourraient le confirmer (contraintes de cisaillement et forces tractrices). La granulométrie des éléments grossiers hérités des conditions hydro-climatiques du Quaternaire et leur distribution sont variables tant dans le fond du lit que dans les berges, comme cela a été montré ; ces caractères n’ont pas permis de proposer un diamètre médian de galets représentatif pour un tronçon de cours d’eau. Par ailleurs, la très faible inclinaison des lignes d’eau nécessite des mesures complémentaires pour saisir leur modification lorsque les conditions d’écoulement varient (conditions hydrologiques ou influence anthropique liées aux manipulations associées aux barrages). Ainsi, l'essentiel de l'énergie de la rivière, en particulier celle qui est déterminée par les vitesses les plus fortes à l'intérieur des sections mouillées, se dissipe latéralement, ce qui permet aux berges de reculer principalement par sapement et localement seulement, par éboulement.
27Le fonctionnement du Loir répond ainsi à des règles générales révélées par les longs tronçons de cours d’eau, où la stabilité ne peut être remise en cause que par des évènements catastrophiques. Le tracé de la rivière bouge peu. À plus grande échelle, l’érosion se manifeste de manière très irrégulière et en des secteurs de berges courts et discontinus correspondant à des sites spécifiques. Dans ce fonctionnement hydro-géomorphologique, le rôle des héritages géomorphologiques sur la dynamique fluviale est aussi grand que celui des conditions actuelles, tant du climat que des actions anthropiques.
28Il est probable que l'arasement des barrages envisagé par l'application de la loi sur l'eau ne change pas fondamentalement les conditions de la morphogenèse compte tenu de la faiblesse des pentes. Les changements dans les faciès d’écoulement liés à l’abaissement de la ligne d’eau ainsi que les caractéristiques de la morphologie du fond n'auront probablement qu'un impact modéré. On peut émettre cependant l’hypothèse qu’une partie de l’activité morphogène mise en évidence sur le site de Bruant soit influencée par l’érosion régressive associée au percement du barrage situé 500 m en aval (percement naturel intervenu au cours des années 1980). Mais dans tous les cas, la dynamique fluviale est sous le contrôle de l'énergie, insuffisante pour mobiliser les matériaux grossiers.
29Ainsi sera-t-il nécessaire de compléter les mesures que nous venons d’analyser, d’une part par une approche spatio-temporelle se rapportant à de plus longs tronçons de rivière et suivis au cours d'évènements hydrologiques plus nombreux. Et d’autre part par des mesures de ligne d’eau de précision sub-centimétrique. Il s'agit ainsi de mieux prendre en compte l'hétérogénéité des conditions de l’écoulement et de ses implications hydro-géomorphologiques.
30Les résultats obtenus sur le Loir contribuent à la compréhension du fonctionnement des cours d'eau de faible énergie, du moins dans ses grandes lignes. S’ils ne sont valables que pour les conditions dans lesquelles nos mesures ont été faites (pas de forte crue survenue au cours de la période de mesure, débit à pleins bords atteint une seule fois) ils sont cependant révélateurs des processus qui participent à la dynamique fluviale : très grande discontinuité des manifestations de l’érosion, tant spatiales que temporelles, prédominance de l’ablation sur les berges et formes construites très rares, forte influence des héritages et de l’aménagement des lits fluviaux. Par ailleurs, les travaux réalisés sur le Loir ont un caractère exploratoire. Ils ont permis que soient mises au point des méthodes de traitement des données de terrain adaptées aux caractéristiques de l’érosion des rivières de faible énergie. La chaîne d’acquisition pour quantifier les transports sédimentaires et l’évolution géomorphologique du Loir a combiné pour partie la chaîne de mesures développée dans le cadre des travaux de thèse de doctorat de Nicolas Claude sur les processus et flux hydro-sédimentaires en rivière sablo-graveleuse (Claude, 2012) et une méthode d’auscultation des berges par scanner 3D développée au sein du Laboratoire Géomatique et Foncier du CNAM (traitements importants tant au niveau du nettoyage des données que du développement de méthodes de quantification). Une approche pluridisciplinaire a été appliquée afin d’appréhender au mieux le fonctionnement hydro-sédimentaire et morphologique du Loir. Mais il est encore nécessaire d’une part de continuer à affiner les mesures à grande échelle afin d’avoir une connaissance plus fine des phénomènes d’érosion, et d’autre part de développer une méthode pour passer de l’échelle du site à celle du tronçon et du cours d’eau afin de comprendre et de quantifier plus largement les processus du transport sédimentaire sur le Loir.
▪ some local shoreline recessions ranging from 0.80 to 1.40 m in length, which appeared in the course of a few months during the high-water period (in the absence of major floods),
▪ some minor changes of the riverbed floor, including a deepening of 1.5 to 1.8 m depth between permanently submerged shallow areas.