3 Méthode et outils

Le travail s’est organisé en trois axes : un premier, bibliographique, constitué d’un état de l’art des connaissances actuelles du fonctionnement de l’étang du Curnic et du fonctionnement global des systèmes lagunaires et des problématiques dystrophiques qui leur sont associées ; un deuxième constitué d’analyses statistiques et, en filigrane des deux axes précédents, un travail d’acquisition de données sur le terrain.

3.1 Travail bibliographique et état de l’art des connaissances sur l’étang du Curnic

Un premier travail bibliographique a porté sur la connaissance des systèmes lagunaires dans leur globalité. Ainsi des recherches ont été conduites sur les aspects hydrodynamiques, écologiques, biologiques, physico-chimiques, biogéochimiques et de gestion des lagunes. L’outil de recherche bibliographique Scopus a été utilisé.

Un état de l’art des connaissances concernant l’étang du Curnic a également été mené.

3.2 Les protocoles d’acquisition de données

Le comité de pilotage (COPIL) du site Natura 2000 du 9 mai 2022 a permis de rappeler les questionnements spécifiques de l’étude, nous les reproduisons ci-dessous :

  • Doit-on modifier des franchissements en faveur de la Loutre d’Europe ?

  • De quelles espèces est constituée la population piscicole ? Et dans quel état est-elle ?

  • La tourbière et les prairies humides ont-elles un rôle épuratoire vis-à-vis du bassin versant ?

  • L’Alanan est-il l’apport majeur des nitrates vers l’étang ?

  • L’étang joue-t-il vraiment un rôle épuratoire conséquent ?

  • Quel volume d’eau l’étang peut-il accueillir ?

  • La roselière peut-elle permettre un reméandrage de l’Alanan ? Cela permettra-t-il de renforcer le processus épuratoire ?

  • Comment fluctuent la température et la conductivité dans l’étang ? Existe-il un fonctionnement par secteurs ?

  • Comment circulent les eaux (eau douce et eau de mer) dans l’étang ?

Afin de répondre à ces questionnements, des protocoles ont été mis en place et des études ont été conduites. Pour certaines problématiques, il a été décidé de faire appel à un bureau d’étude ou à des collaborateurs.

Le tableau ci-dessous récapitule les protocoles, outils et études pour chaque thématique abordée dans le cadre de ce projet ainsi que les indicateurs associés.

Thématiques Outils Indicateurs
Propriétés physico-chimique de l’eau
  • Sonde autonome

  • Sonde multiparamètre mobile

  • Analyse en laboratoire

  • Bandelette nitrate

  • Azote global

  • Azote organique et ammoniacale

  • Concentration en nitrates

  • Concentration en nitrites

  • Concentration en orthophosphate

  • Phosphore global

  • Salinité

  • Conductivité (et résistivité)

  • Température de l’eau

  • TDS (Total Disolved Solid)

  • Pression

Hydrodynamique
  • Sonde autonome

  • Vélocimètre

  • Niveau d’eau

  • Débit

Contexte physique de l’étang
  • Scop AQUABIO
  • Bathymétrie
  • Topographie des sédiments
Contexte physique du marais
  • Collaboration avec le laboratoire LETG
  • Topographie
Contexte climatique
  • Pluviomètre à la station d’épuration de Guissény / Kerlouan
  • Pluviométrie
Diagnostic piscicole
  • Scop AQUABIO
  • Espèces en présence

  • État de santé de la population

La figure 3.1 présente l’ensemble des points de mesures/prélèvement par protocoles.

Figure 3.1: Localisation de l’ensemble des points de mesure-prélèvement. Il s’agit d’une carte interactive, vous pouvez choisir le fond de carte de votre choix et décocher les couches que vous souhaitez ; il est possible de cliquer sur les points, une infobulle s’affichera

3.2.1 Contexte climatique

Nous avons présenté les données utilisées dans la section 2.1.4.

3.2.2 Propriété physico-chimique de l’eau

L’objectif est d’évaluer la qualité de l’eau et l’évolution des concentrations/flux en nutriments entre le bassin versant, le réseau du marais et l’étang au cours d’une année ; caractériser l’évolution des températures et de la salinité selon les conditions de marées et les saisons et mieux comprendre les éléments biotiques et abiotiques qui influencent ces fluctuations physico-chimiques.

Concernant la qualité de l’eau, deux méthodes d’acquisition de données sont utilisées : l’analyse en laboratoire et les bandelettes nitrate.

3.2.2.1 Analyse d’eau par laboratoire

Afin d’évaluer les teneurs de différentes formes de l’azote et du phosphore, des prélèvements d’eau sont réalisés à l’embouchure du drain principal du marais dans l’étang (point 2) et à l’exutoire de ce dernier (point 1) (fig. 3.2).

Dans la mesure du possible, les prélèvements sont réalisés le même jour afin de parer au biais engendré par les variations temporelles des concentrations. L’heure est quant à elle sélectionnée en fonction de la marée de telle sorte à prélever dans des conditions similaires. Les lendemains de fortes précipitations sont également évités afin de ne pas obtenir de valeurs extrêmes. Les échantillons sont ensuite transmis au laboratoire Labocea où les éléments suivants sont mesurés : azote total, nitrates, nitrites, l’azote kjeldahl, orthophosphates et phosphore total.

Dans le cadre du suivi de la qualité de l’eau réalisé par le Syndicat Mixte du Bas-Léon, le service environnement de la mairie de Guissény a produit une demande, en passant par la Communauté Lesneven Côte des Légendes (CLCL), pour accorder le nombre de prélèvements, la rythmicité et les éléments analysés qu’il réalise au lieu-dit Cléguer (point 5) avec ceux réalisés dans le cadre de cette étude.

Localisation des points de prélèvement et de mesure

Figure 3.2: Localisation des points de prélèvement et de mesure

3.2.2.2 Bandelette nitrates

Pour compléter la connaissance des flux en nitrate, desde bandelettes nitrate sont utilisées pour les points 1, 2, 3, 4 et 5 (fig. 3.4).

Carte *Deltares* et bandelette nitrate

Figure 3.3: Carte Deltares et bandelette nitrate

Le système utilisé est celui de la marque Deltares médiatisée par l’application Nitrate app. L’incertitude calculée est de 8 % (laboratoire géosciences de Rennes).

Localisation des points de mesure des bandelettes nitrate

Figure 3.4: Localisation des points de mesure des bandelettes nitrate

Deux campagnes (04/07/2022 et 04/08/2022) sur l’ensemble du bassin versant ont également été effectuées afin de posséder une première image des concentrations en nitrates sur celui-ci (fig. 3.1).

3.2.2.3 Sonde multiparamètre

La sonde multiparamètre permet d’acquérir des données de conductivité (et conductivité absolue), de salinité, de pression dans l’eau et dans l’air, de TDS (Total Disolved Solid) ainsi que de température. La sonde utilisée est le modèle HI98199 de la marque Hanna Instruments. Les points mesurés sont situés à la fois dans l’étang et dans le marais (fig. 3.6).

Mesure avec la sonde multiparamètre

Figure 3.5: Mesure avec la sonde multiparamètre

L’objectif est de caractériser le gradient de salinité horizontale dans l’étang afin d’identifier l’ampleur de l’influence de la marée et du coin salé dans le système.

Les mesures ont été effectuées une fois par mois depuis le 14 octobre 2021.

Localisation des points de mesure de la conductivité

Figure 3.6: Localisation des points de mesure de la conductivité

Un protocole a été ajouté au cours de l’étude afin de caractériser plus finement l’influence de la marée. Dix points ont été ajoutés aux initiaux, ils sont cartographiés dans la figure (3.6). Un échantillonnage stratifié22 a été utilisé. Au total, quatre sessions ont été effectuées pour des contextes de marées basses et hautes par coefficient de mortes-eaux et de vives-eaux.

3.2.2.4 La sonde autonome

Une sonde autonome a été placée le 1 juin 2021 à l’exutoire de l’étang (point 6, fig. 3.6). Le modèle utilisé est le CTD de la marque OTT. Elle mesure le niveau d’eau, la conductivité, la salinité, les TDS et la température tous les quarts d’heure et moyenne les résultats toutes les demi-heures. Les mesures sont enregistrées à une hauteur comprise entre 15 et 20 cm au-dessus de la surface du sédiment.

L’objectif est de quantifier l’ensemble des influences à la fois marines et continentales sur les indicateurs mesurés à l’exutoire dans l’objectif de modéliser son fonctionnement hydrodynamique et physico-chimique.

La sonde a été déplacée le 28 juillet 2022 au point 12, au sud de l’étang (fig. 3.6).

3.2.3 Outils d’acquisition des données hydrauliques

3.2.3.1 Vélocimètre

Les mesures de débit ont été réalisées au vélocimètre. Le modèle utilisé est le FW01 de la marque COMETEC, la sensibilité théorique minimum de l’instrument est de 0,1 \(m.s^{-1}\).

Les données de débits permettront de calculer le bilan hydrique de l’étang et de déterminer les flux en éléments chimiques.

Trois sections sont mesurées, une à l’exutoire (point 1), une à l’exutoire du drain principal du marais, au sud de l’étang (point 2) et une au lieu-dit “Cléguer” avant son passage sous la route (fig. 3.2).

Le protocole utilisé est celui de l’Office National de l’Eau et des Milieux Aquatiques (Le Coz et al., 2011) (fig. 3.7). Les débits ont été calculés à l’aide de la méthode de la section médiane.

 Exemple de répartition des points de mesure dans la section de jaugeage pour une section large. ONEMA, 2011

Figure 3.7: Exemple de répartition des points de mesure dans la section de jaugeage pour une section large. ONEMA, 2011

Les débits sont mesurés une fois par mois en hiver et deux fois par mois au printemps et en été. Les débits (\(Q\)) sont exprimés en \(m^{3}.s^{-1}\).

3.2.3.2 Sonde autonome

La sonde autonome calcule la variation du niveau d’eau à l’exutoire de l’étang. L’outil est présenté dans la section 3.2.2.4.

3.2.4 Contexte physique de l’étang

3.2.4.1 Bathymétrie et topographie des sédiments

L’objectif est de connaître le volume d’eau accueilli ainsi que la répartition des sédiments dans l’étang. Le volume d’eau permettra de calculer le bilan hydrique de l’étang et le temps de résidence des eaux dans le système.

Un levé bathymétrique a été réalisé par la Scop AQUABIO le jeudi 31 mars 2022. L’ensemble de l’étang ainsi que le canal de l’Alanan ont été levés.

Le relevé consiste à circuler sur le plan d’eau à l’aide d’une embarcation, munie d’un sondeur couplé à un GPS, qui permet d’enregistrer des points pour lesquels sont alors connues les coordonnées GPS et la profondeur (Le Ruyet, 2022).

La topographie du fond a également été mesurée « à l’aide d’une perche graduée équipée d’un système coulissant en fonction de la hauteur de sédiments » (Le Ruyet, 2022).

À l’issue de l’étude, un rapport et deux MTN ont été fournis, un concernant la bathymétrie et un pour la topographie des sédiments.

3.2.5 Contexte physique du marais

3.2.5.1 La topographie

Le levé topographique concerne la zone, en bordure du canal de l’Alanan, où ont été creusés les bassins de pisciculture dans les années 70 ; la zone représente environ 6 ha. Il s’agit d’un protocole prospectif dont la finalité est d’approximer les volumes nécessaires à exporter — et les coûts associés — dans l’hypothèse d’un détournement de l’Alanan dans la roselière afin de favoriser l’abattement de l’azote.

La réalisation de ce levé est à envisager hors période de nidification des oiseaux d’eau et de roselière. La fin de l’hiver ou bien l’automne est donc à privilégier.

Les levés sont réalisés avec la collaboration du laboratoire LETG, IUEM (interlocuteur, Serge Suanez).

3.2.6 Diagnostic piscicole

L’objectif du diagnostic piscicole est d’établir un état des lieux de la population piscicole de l’étang. L’ichtyofaune23 est l’un des indicateurs biologiques requis par la Directive Cadre sur l’Eau (DCE) pour évaluer l’état écologique des masses d’eau.

La Scop AQUABIO a réalisé ce diagnostic entre le 13 et le 14 juin 2022.

L’étang du Curnic étant saumâtre, l’emploi d’engins passifs est indispensable. Le prestataire a utilisé des verveux doubles sur quatre stations réparties en différents points de l’étang.

3.3 Traitement des données

Dans un premier temps, un travail de gestion de données a été effectué. Les données préexistantes ont été normalisées et stockées dans un drive24 et enrichies par celles acquises au cours du stage. Ces données ont été importées directement dans l’environnement RStudio à l’aide package googlesheets4 (Bryan & RStudio, 2021). Les scripts ont été rédigés de telle sorte à ce que leurs sorties soient mises à jour en fonction de la base de données ; ainsi, les membres du service environnement pourront les réutiliser sans avoir à les modifier si la base continue à être alimentée.

Le travail statistique proprement dit a commencé par une analyse descriptive (univariée, bivariée et multivariée) des données issues des protocoles présentés plus hauts.
Dans un second temps, des modélisations stochastiques25 de séries temporelles26 ont été réalisées sur les signaux enregistrés par la sonde autonome.

En parallèle de cette analyse statistique, les débits, le bilan hydrologique, le temps de résidence, les flux, les bilans et les valeurs de l’abattement ont été calculés. Dans un souci de lisibilité, les méthodologies propres à chaque résultat seront présentées en amont de chacun d’entre eux.

L’ensemble des traitements ont été effectués sous les logiciels Microsoft Excel, RStudio et QGIS.

References

Bryan, J., RStudio, 2021. Googlesheets4: Access Google Sheets Using the Sheets API V4.
Le Coz, J., Camemen, B., Dramais, G., Ribot-Bruno, J., Ferry, M., Rosique, J.-L., 2011. Contrôle Des Débits Réglementaires. Application de l’article L. 214-18 Du Code de l’environnement. ONEMA, Cemagref.
Le Ruyet, O., 2022. Étude Bathymétrique et Évaluation de Volumes de Sédiments. Étang Du Curnic et Alanan Aval à Guissény (29). AQUABIO.
Raymond, S., 2011. Incertitudes des flux transportés par les rivières (Matière en suspension, nutriments, sels dissous) Vers un système expert d’optimisation des méthodes de calcul. thèse de doctorat, Université François Rabelais - Tours.

  1. « L’échantillonnage stratifié consiste à prélever à des moments privilégiés. » (Raymond, 2011).↩︎

  2. L’ichtyofaune est la partie de la faune rassemblant les poissons. (wikipedia, cons. 05/08/2022 12:09)↩︎

  3. Un drive est un dossier sur le site googledrivepartageable à des tiers.↩︎

  4. La notion de stochasticité fait référence à l’aléatoire dans un système donné ; les effets “imprévisibles”. La modélisation stochastique s’oppose à la modélisation déterministe qui n’admet pas d’effet aléatoire.↩︎

  5. Une série temporelle, ou série chronologique, est une suite de valeurs numériques représentant l’évolution d’une quantité spécifique au cours du temps. (wikipedia, cons. 06/08/2022 12:03)↩︎