Bertrand Villeneuve et ses collègues (UR MALY, Irstea, Laboratoire d'hydro-écologie quantitative ; Laboratoire interdisciplinaire des environnements continentaux, UMR 7360 CNRS—Université de Lorraine) viennent de publier une nouvelle recherche à propos des impacts des différentes pressions sur les invertébrés aquatiques en rivière. Leur problématique est ainsi énoncée : "Le but de notre approche était de prendre en compte les échelles spatiales imbriquées pilotant le fonctionnement des rivières dans la description des liens pressions / état écologique, en analysant les résultats d'un modèle hiérarchique. Le développement de ce modèle nous a permis de répondre aux questions suivantes: La prise en compte des liens indirects entre les pressions anthropiques et l'état écologique des cours d'eau modifie-t-elle la hiérarchie des types de pression impactant les invertébrés benthiques? Les différentes échelles imbriquées jouent-elles des rôles différents dans la relation pressions anthropiques / statut écologique? Ce modèle permet-il de mieux comprendre le rôle spécifique de l'hydromorphologie dans l'évaluation de l'état écologique des cours d'eau?"
Les chercheurs ont sélectionné des hydro-écorégions de plaines et piémont (moins de 450 m), en divisant l'échantillon en deux critères géologiques : calcaire (roches sédimentaires) ou non calcaire. Une autre division, hydrologique, distingue les cours d'eau de petite dimension (ordre 1 à 3 de Strahler) et de moyenne dimension (ordre 4 à 6). Au final, quatre groupes sont analysés à partir de mesures réalisées sur 5 ans (2007-2012) : petites rivières non calcaires (160 sites, 638 échantillonnages), moyennes rivières non calcaires (127, 492), petites rivières calcaires (228, 817) et moyennes rivières calcaires (128, 460). Soit 643 sites et plus de 2400 campagnes d'échantillonnage.
La donnée biologique étudiée est l'indice invertébrés I2M2 (Mondy et al 2012), qui mesure les macro-invertébrés benthiques par la taxonomie et les traits des espèces, produisant un score de qualité en comparant des zones à faible et fort impact humain.
Les données de contexte et pression sont réparties en 3 échelles spatiales emboitées: le bassin versant (8 descripteurs, par exemple urbanisation, agriculture intensive, érosion, irrigation, etc.), le tronçon (13 descripteurs, par exemple digues, barrages, rectification, densité d'arbre en ripisylve), le site (8 descripteurs dont les données physico-chimiques type nitrates et phosphores, les matières en suspension, la mosaïques des substrats du lit).
La méthode statistique utilisée est l'approche PLS (Partial Least Squares, régression partielle moindres carrés), une modélisation qui cherche à mettre en correspondance des données manifestes et des données latentes, en tenant notamment en compte les influences indirectes (quand une variable A modifie une variable B par une action intermédiaire sur une variable C).
Les données du modèle utilisé par Villeneuve et al 2018, art cit, droit de courte citation. Cliquer pour agrandir.
Poids direct et indirect des facteurs de variation de l'I2M2, selon les types de rivière, art cit, droit de courte citation. Cliquer pour agrandir.
Les chercheurs soulignent : "En se focalisant sur l'effet total (direct + indirect) des variables latentes sur les valeurs I2M2, si les variables à impact majeur restent les concentrations en nutriments et en matière organique pour les petits cours d'eau non calcaires, les contributions relatives des effets indirects modifient l'ordre hiérarchique des impacts des autres variables latentes pour les autres types de flux. En effet, pour les petits cours d'eau non calcaires, la contribution décroissante à la variation des valeurs I2M2 expliquée par le modèle sont: nutriments et matière organique (42%), usages des sols du bassin versant (21%), altérations hydromorphologiques à l'échelle du bassin versant (16%), les mosaïques du substrat (15%) et les altérations hydromorphologiques à l'échelle du tronçon (6%). En revanche, dans les cours d'eau de taille moyenne non calcaire, l'ordre décroissant d'importance des impacts a été modifié en: usages des sols du bassin versant (29%), altérations hydromorphologiques à l'échelle du tronçon (29%), nutriments et matière organique (18%), mosaïques (17%) et altérations hydromorphologiques à l'échelle du bassin (7%). Pour les petits cours d'eau calcaires, cet ordre était: usages es sols du bassin hydrographique (33%), mosaïques du substrat (25,0%), nutriments et matières organiques (19%), altérations hydromorphologiques du tronçon (15%) et bassin versant (8%). Pour les cours d'eau de taille moyenne, il s'agissait de: usage des sols (35%), nutriments et matière organique (24%), mosaïque de substrat (18%), altérations hydromorphologiques du tronçon (18%) et du bassin (5%)".
Concernant l'hydromorphologie, les chercheurs observent : "le total des contributions directes des variables hydromorphologiques, aussi bien à l'échelle du bassin versant que du tronçon, représentait de 13% à 23% de la variance totale des valeurs I2M2 expliquées par les modèles. Ces contributions sont passées de 6% (cours d'eau de petite taille) à 13% (cours d'eau de taille moyenne non calcaire) en tenant compte également des effets indirects du bassin versant et de l'hydromorphologie sur les caractéristiques physico-chimiques et les mosaïques des sites, fournissant des contributions totales de 22% à 36% (selon les types de rivière) de la variance expliquée dans les valeurs I2M2."
Poids direct et indirect des échelles spatiales du site, du tronçon ou du bassin, art cit, droit de courte citation. Cliquer pour agrandir.
Discussion
Concernant le modèle proposé dans cette publication, plusieurs points nous sembleraient intéressants à approfondir :
- les scores I2M2 viennent en général du réseau de surveillance de la directive européenne sur l'eau, dont les sites d'implantation ne sont pas représentatifs de tout leur bassin. Il peut y avoir un biais d'échantillonnage (de même qu'il existe une marge d'erreur dans le calcul des scores invertébrés, voir la thèse Wiederkehr 2015 dont l'impact sur la significativité des variations est peu étudiée en général). Des travaux sur un moins grand nombre de rivières, mais avec des mesures plus réparties sur leurs cours pourraient affiner les résultats du modèle;
- une absente dans les entrées du modèle reste la pollution des cours d'eau par les substances autres que les nutriments. Il a été montré par la recherche que les populations européennes d'invertébrés terrestres sont en déclin tendanciel depuis plusieurs décennies, parfois prononcé, avec en ce cas de fortes suspicions sur le rôle des pesticides de synthèse. Mais bien d'autres micropolluants terminent dans les eaux, sans que l'on connaisse leur impact. Ce critère est certes approché par les taux d'urbanisation et d'agriculture, mais les pratiques d'épuration et les types d'agriculture sont assez variables. On soulignera à ce sujet que si le modèle prédit bien la direction de variation des scores I2M2, il est loin d'expliquer toute la variance observée (le R2 est de 33% pour les petits cours d'eau non calcaires, 50% pour les cours d'eau non calcaires, 44% pour les petits cours d'eau calcaires et 40% pour les cours d'eau calcaires);
- le modèle donne une photographie instantanée des bassins. Il serait intéressant de développer une approche plus dynamique où, pour les scores invertébrés dont on dispose de séries pluridécennales et homogènes, l'évolution observée des insectes est mise en lien avec celle des pressions (voir le travail de Van Looy et al 2016). Comme il s'agit ici d'écologie appliquée, avec des avis donnés aux gestionnaires, cette approche dynamique contribuerait à séparer plus efficacement les mesures qui risquent d'avoir peu d'effets (voire des effets négatifs) et les autres;
- enfin, une attente forte réside dans l'évaluation détaillée de l'impact morphologique. Cette catégorie regroupe beaucoup de pressions différentes, comme on l'a vu ci-dessus dans la description du modèle. Le gestionnaire tient depuis une quinzaine d'années un discours sur l'égale importance de la pollution et de la morphologie. Cela peut s'entendre, mais le problème est que la morphologie désigne en réalité l'usage des eaux, des berges et des sols sur toutes les échelles (site, tronçon, bassin) : à ce niveau de généralité, on n'est guère avancé dans la décision! Il faudrait donc hiérarchiser plus finement le poids des pressions morphologiques, tant pour l'atteinte problématique des objectifs DCE 2027 que pour la bonne information du débat démocratique sur la rivière.
Ce modèle de Bertrand Villeneuve et ses collègues montre la complexité et la difficulté d'une étude d'impact sur les bassins versants. Certains lecteurs nous demandent parfois pourquoi nous jugeons les documents de la littérature grise (comme les états des lieux des SDAGE) insuffisants comme outils de décision : cette étude leur apportera un début de réponse. Une politique écologique commence par des mesures in situ et des modélisations pour comprendre le bassin versant sur lequel on investit de l'argent public. Par exemple, dépenser beaucoup d'argent à échelle de sites sans prendre en compte des altérations à échelle supérieure (tronçon, bassin) risque de produire des résultats modestes ou nuls, ce qui est d'ailleurs souvent observé en analyse avant-après d'interventions sur la morphologie (voir par exemple Morandi 2014, Lorenz et al 2018 en Allemagne, cette synthèse 2005-2015). Comme les variations de populations d'insectes aquatiques ne sont pas vraiment la priorité des citoyens, le choix d'interventions peu sensées et le risque de résultats insignifiants vont altérer un consentement à payer déjà assez modeste.
On notera à ce sujet que l'I2M2, comme tous les indicateurs DCE fondés sur l'état de référence, souffre d'une certaine circularité dans sa construction. On s'attend à ce qu'une rivière et un bassin modifiés par l'homme ne produisent pas les mêmes assemblages d'espèces que d'autres très peu modifiés. Le fait de nommer "dégradation" ce changement est un jugement de valeur davantage qu'un jugement scientifique, et il revient in fine à dire (par des moyens un peu complexes) qu'une "bonne" rivière serait une rivière sur laquelle l'homme intervient un minimum (paradigme de la "nature sans l"homme" comme référence idéalisée). Il reste encore à expliquer pourquoi, c'est-à-dire en vertu de quel jugement social partagé certains assemblages d'insectes (dans le cas de l'I2M2) sont préférables à d'autres. On peut douter que l'ovoviviparité ou le polyvotinisme d'une métapopulation d'invertébrés motive en soi un grand nombre de riverains à agir ! Si l'insecte témoigne d'une pollution également dommageable à l'être humain, sa variation a davantage de sens. S'il témoigne d'une évolution physique ou chimique de l'eau sans effet notable sur l'homme, l'enjeu paraît plus difficile à justifier en coût des politiques publiques. Les scores plus aisés à partager seraient certainement la biomasse et la biodiversité (perçues comme signes d'une nature diverse et en bonne santé), mais ils n'apparaissent plus dans le score unique agrégé des indices composites de type I2M2. Dans l'étude ici commentée, les variations des sous-scores composant l'I2M2 ne sont pas détaillées, ce qui est dommage.
Enfin, bien qu'il s'adresse au gestionnaire dans sa conclusion, ce travail apporte surtout des éléments d'intérêt à une théorie de l'évolution écologique des communautés aquatiques, selon une approche systémique. Il est en cela très utile, puisqu'il améliore l'intelligibilité de la dynamique des écosystèmes sous influence anthropique. Sa valeur restera en revanche limitée pour informer les débats très locaux, qui concernent l'écologie de sites particuliers. Par exemple, notre association fait régulièrement observer aux syndicats ou autres gestionnaires que la disparition des ouvrages hydrauliques transversaux (comme suppression d'un impact morphologique en lit mineur) fait dans le même temps disparaître divers micro-habitats secondaires, aquatiques ou humides, qui profitent eux aussi à certaines espèces, dont des invertébrés. Le score I2M2 de la station en lit mineur et la réalité de la biodiversité du site avant-après ne sont pas les mêmes mesures (en l'occurrence, la biodiversité ordinaire perçue par les riverains n'est pas celle de l'I2M2). Les modèles hiérarchisés de bassin versant peuvent donc nourrir des réflexions pour les politiques publiques répondant à des objectifs cadres (comme ceux de la DCE), mais ils ne permettront pas l'économie d'échanges bien plus précis sur la rivière et la biodiversité que veulent les riverains. Et cela ne se décide pas dans les bureaux de Bruxelles...
Référence : Villeneuve B et al (2018), Direct and indirect effects of multiple stressors on stream invertebrates across watershed, reach and site scales: A structural equation modelling better informing on hydromorphological impacts, Science of the Total Environment, 612, 660–671