1À la demande du Conseil général des Côtes d’Armor, un suivi topo-morphologique à haute fréquence des deux cordons dunaires du Grand Rocher et de Saint-Michel-en-Grève, situés en baie de Saint-Michel-en-Grève (fig. 1), a été entrepris entre les mois d’octobre 2002 et septembre 2003. Ce travail fut décidé alors que des phénomènes d’érosion affectaient ces derniers (Suanez, 2004 ; Stéphan et Suanez, 2004), afin d’analyser et de comprendre les mécanismes à l’origine de leur recul.
Fig. 1 – Carte de localisation.
Fig. 1 – Location map.
2Cet article présente les résultats obtenus durant cette année de suivi topo-morphologique et examine le rôle des forçages météorologiques dans l’évolution des cordons dunaires étudiés. Les travaux réalisés sur le fonctionnement saisonnier des milieux dunaires insistent sur la complexité des processus alternant entre phase d’érosion et phase d’alimentation (Davidson-Arnott et Law, 1990 ; Pye, 1991 ; Sherman et Bauer, 1993). Comme de nombreux auteurs l’ont montré, les tempêtes, épisodes morphogènes générant une élévation brutale et temporaire du niveau marin jouent un rôle important dans l’érosion des dunes (Sherman et Nordstrom, 1985 ; Carter et al., 1990 ; Pye et Neal, 1994 ; Arens, 1996a). Ces phénomènes érosifs sont renforcés lorsqu’une onde de tempête se combine à une pleine mer de vive-eau (Ruz et Meur-Ferrec, 2004 ; Battiau-Queney, 2004). Toutefois, comme le rappelle P.A. Pirazzoli (2000), aucune situation combinant une très forte tempête et une marée de vive-eau ne s’est encore produite sur le littoral nord atlantique, à l’exception peut-être de la tempête de l’hiver 1989-1990 qui a eu sur le littoral breton des effets spectaculaires et dévastateurs (Fichaut et Hallégouët, 1989). Les travaux réalisés sur les relations entre les niveaux marins extrêmes et l’érosion des dunes apportent un canevas méthodologique intéressant (Ruggiero et al., 1997 ; Dail et al., 2000 ; Ruggiero et al., 2001). Ces auteurs soulignent le rôle important des fortes houles et des phénomènes de run-up qui leur sont associés dans la dynamique de recul des fronts de dune. Pour autant, cette érosion n’est pas totalement irréversible car des phases d’alimentation et de régénération des cordons dunaires peuvent prendre place en période estivale où l’activité morphogène est considérablement réduite (Psuty, 1988 ; Davidson-Arnott et Law, 1990 ; Arens, 1996a, 1997 ; Ruz et Meur-Ferrec, 2004). La réversibilité des effets de l’érosion peut également s’inscrire dans un temps beaucoup plus long, pluriannuel, comme l’ont montré des suivis réalisés sur diverses plages (Thom et Hall, 1991 ; Battiau-Queney et al., 2000). Aussi, une bonne connaissance de ces fonctionnements est une aide indispensable à la gestion rationnelle de ces milieux (Hesp, 2002 ; Ruz, 2004 ; Battiau-Queney, 2004).
3D’un point de vue méthodologique, deux types de levés ont été réalisés afin de mesurer la variabilité saisonnière du bilan sédimentaire des cordons dunaires étudiés. Le premier repose sur la mesure de la superficie dunaire donnant ainsi une vision globale de l’évolution en plan des cordons. Le second s’est attaché à mesurer les changements morphosédimentaires du système cordon dunaire / haut de plage le long de six transects perpendiculaires au trait de côte permettant ainsi d’analyser le comportement du front de dune dans sa dimension verticale. L’étude des conditions de forçages météo-marins a été réalisée à partir de données Météo France pour le vent et la pression barométrique, d’enregistrements marégraphiques fournies par le SHOM (Service Hydrographique et Océanographique de la Marine-Brest) et de données de houle issues d’une modélisation acquises auprès du LNHE-EDF (Laboratoire National d’Hydraulique et Environnement de l’EDF – Chatou). Dans un dernier point, l’analyse des données morphosédimentaires et des conditions de forçages météo-marins s’est faite à partir d’une série d’Analyses factorielles en Composantes Principales (ACP).
4La baie de Saint-Michel-en-Grève correspond à une plage de fond de baie cernée par de hautes falaises dont l’altitude décroît d’ouest en est de quatre-vingts mètres à une cinquantaine de mètres. Elle est ouverte aux flux océaniques d’ouest à nord-ouest et il n’est pas rare d’avoir des régimes de nord-est, principalement au printemps. Toutefois, la présence des pointes de Berg ar Forn à l’est et de Beg Douar à l’ouest confère à la baie un caractère abrité limitant les impacts morphogènes. L’importance du marnage (9,3 m à Locquirec) explique la grande superficie de l’estran (7 km2) qui peut découvrir sur plus de deux kilomètres de distance. Les cordons dunaires occupent un espace extrêmement réduit (25 000 m2 environ, soit 0,04 % du système littoral). Cela s’explique en grande partie par le « blocage » du système littoral par les falaises de fond de baie qui n’a pas permis l’édification d’une zone supralittorale étendue où s’exercerait la dynamique éolienne. À cela s’ajoute également le rôle des aménagements avec la construction de la route littorale D786. Cette dernière a repoussé « artificiellement » le trait de côte vers l’aval du système entraînant ainsi une submersion plus fréquente du haut d’estran et, avec elle, la destruction des premières dunes (Gad et al., 2003). Les seules dunes holocènes que l’on puisse trouver en baie de Saint-Michel-en-Grève se situent dans les rentrants que forment les débouchés des vallées fluviales. Elles sont aujourd’hui totalement masquées par des constructions d’origine anthropique que l’on a édifiées sur leurs revers (habitats, activités touristiques balnéaires, etc.).
5La baie de Saint-Michel-en-Grève connaît un engraissement très important que l’on mesure particulièrement bien depuis le début du xixe siècle (Suanez et al., 2002 ; Stéphan et Suanez, 2005). Depuis 1990, ces apports sédimentaires sableux ont entraîné un exhaussement suffisamment conséquent du haut de plage pour que se mettent en place des dunes bordières (Stéphan et Suanez, 2004). Elles se répartissent en trois cordons dunaires bien distincts. Les deux premiers, les cordons dunaires de Saint-Efflam et du Grand Rocher, se situent sur la commune de Plestin-les-Grèves, le troisième, de Saint-Michel-en-Grève, est situé à l’est-nord-est sur la commune du même nom (fig. 1). D’un point de vue morphologique, ces « micro-reliefs dunaires » correspondent à des placages éoliens très peu élevés dont la partie externe forme un bourrelet inférieur à un mètre d’altitude par rapport au niveau des plus hautes mers. La granulométrie de ces derniers est très homogène d’un secteur à l’autre de la baie. Les mesures effectuées sur l’ensemble des sables dunaires montrent des variations très faibles comprises entre 163 µm et 142 µm, pour un grain moyen égal à 155 µm. C’est en général un sédiment bien trié comme l’indiquent les indices de tri égal à 0,41 et d’asymétrie proche de 0 (Stéphan, 2003). Notre attention portera sur les deux cordons dunaires du Grand Rocher et de Saint-Michel-en-Grève qui font l’objet de cet article.
6Le cordon dunaire du Grand Rocher occupe la partie centrale de la baie. Sa superficie ne dépasse pas 1 900 m2 et sa hauteur reste inférieure à 1,5 m. Il présente une morphologie bosselée constituée de bourrelets peu prononcés parallèles au rivage et séparés les uns des autres par des secteurs plans de type « plate-forme » (fig. 2). Son orientation et sa localisation font qu’il est protégé des flux d’ouest par le promontoire du Grand Rocher. À l’inverse, il est particulièrement bien exposé aux flux d’est à nord-est. Il apparaît extrêmement végétalisé (taux de couverture végétale > 90 %). Les formations végétales qui le couvrent sont des plus communes : dans la partie externe les graminées prédominent (Agropyrum junceinforme et Carex arenaria), la partie interne, proche de l’enrochement, est quant à elle couverte d’associations de pelouse rase (plantin, euphorbe, fucus et des lichens).
Fig. 2 – Morphologie du cordon dunaire et du haut de plage du Grand Rocher.
Fig. 2 – Morphology of the coastal dune and upper beach of Grand Rocher
(cliché : Serge Suanez)
(photo : Serge Suanez).
7Le cordon dunaire de Saint-Michel-en-Grève est de taille plus importante ; il s’étend sur 800 m et couvre environ 14 500 m2. Sa hauteur reste inférieure à 1,5 m. La morphologie de cet ensemble est similaire à celle décrite pour le cordon dunaire du Grand Rocher (fig. 3). À l’inverse, son orientation ouest-nord-ouest l’expose principalement aux flux océaniques d’ouest rentrant dans la baie. Son taux de couverture végétale est également très important (proche de 100 %) et présente les mêmes associations vétégales que celles décrites précédemment.
Fig. 3 – Morphologie du cordon dunaire et du haut de plage de Saint-Michel-en-Grève.
Fig. 3 – Morphology of the coastal dune and upper beach of Saint-Michel-en-Grève
(source : http://www.francevuesurmer.com)
(http://www.francevuesurmer.com).
Fig. 4 – Conditions de vents sur la zone d’étude à partir des données enregistrées pour la station météorologique de Lannion.
Fig. 4 – Wind conditions for the study area obtained from data recorded at the Lannion wind station.
8La dynamique éolienne est en grande partie conditionnée par les flux météo-marins d’ouest. Néanmoins, suivant les périodes de l’année, les vents de secteur est à nord-est jouent un rôle non négligeable. Les données de vent obtenues auprès de la station météorologique de Lannion permettent de décrire une « climatologie » générale et saisonnière des vents dans le secteur de la baie de Saint-Michel-en-Grève (fig. 4).
9À l’échelle annuelle, la majorité des vents les plus forts est de direction ouest à sud-ouest (220° à 260°). Ces vents accompagnent le passage des perturbations atlantiques génératrices de tempêtes sur le littoral. Les circulations d’est à nord-est (40° à 100°) représentent la deuxième principale direction des vents annuels. Ils montrent toutefois des vitesses moins importantes que pour les vents d’ouest. À l’échelle saisonnière, une variation non négligeable des directions et des forces de vents est observée. Durant l’hiver, du mois de janvier au mois de mars, les vents de sud-ouest (200° à 260°) sont les plus fréquents et les plus virulents (13 % du régime total compris entre 4 et 8 m/s et 8 % > à 8 m/s). Ces vents forts sont souvent associés à des hauts niveaux marins où peuvent se combiner marée de vive-eau et onde de tempête (surcote). Ces conditions sont propices à l’érosion des dunes. À partir du mois d’avril, les vents d’ouest sont moins intenses. Ils prennent une direction ouest-nord-ouest (260° à 300°). Dans le même temps, on observe une nette recrudescence des vents d’est à nord-est (40° à 80°) toutefois peu vigoureux (majorité du régime < à 4 m/s). Cette période printanière associant des vents morphogènes et une diminution des pluies est plutôt favorable à l’engraissement des dunes, aussi bien de Saint-Michel-en-Grève que du Grand Rocher. En été, du mois de juillet au mois de septembre, les vents faiblissent quelle que soit la direction. Durant cette période « sèche », les vents les plus morphogènes soufflent de l’ouest. Cette circulation est alors propice à l’alimentation des dunes de Saint-Michel-en-Grève. À partir du mois d’octobre, les perturbations atlantiques reviennent et s’accompagnent d’une augmentation de la force des vents de sud-ouest (200° à 260°). Ces flux océaniques sont les plus représentés durant la période automnale. Ils s’accompagnent de hauts niveaux marins (marée d’équinoxe et ondes de tempêtes) capables de créer de graves dommages sur les cordons dunaires.
Tableau 1 – Inventaire des levés topo-morphologiques effectués entre les mois d’octobre 2002 et septembre 2003
Table 1 – List of topographic surveys realised between October 2002 and September 2003
(* à l’exception du profil 03 ; ** à l’exception du profil 06).
(*except profile 03 ; ** except profile 06).
10Entre les mois d’octobre 2002 et septembre 2003, nous avons réalisé onze campagnes de mesures topo-morphologiques des cordons dunaires du Grand Rocher et de Saint-Michel-en-Grève (tab. 1). La fréquence des relevés a été particulièrement élevée durant la période hivernale afin de cerner tous les épisodes de forte activité morphogène. À partir du printemps, lorsque les conditions anticycloniques ont prédominé, la fréquence est passée à un relevé par mois. Ces levés de terrain ont été effectués au tachéomètre laser (Nicon DTM 300 et Leica TCR303). Ils sont composés d’un levé du front de dune afin de mesurer le recul ou l’avancée de ce dernier et de mesurer les surfaces dunaires perdues ou gagnées. La limite retenue pour délimiter le front de dune est clairement matérialisée par une rupture de pente marquant le contact entre la dune et le haut de plage. Cet escarpement dessine également la limite de la végétation dunaire. L’analyse quantitative de ces relevés (calcul des surfaces) a été réalisé sous SIG (ArcView 8.3) suivant une procédure classique d’intégration des données. Les levés de terrain sont aussi composés de mesures effectuées le long de six profils transversaux disposés perpendiculairement au trait de côte englobant la dune et le haut de plage : quatre profils sur le cordon dunaire de Saint-Michel-en-Grève (fig. 3) et deux sur celui du Grand Rocher (fig. 2). À partir de ces données, il a été possible de quantifier les volumes de sable perdus ou gagnés par la dune en utilisant la méthode des « surfaces verticales ». La superposition de deux profils délimite des aires en accrétion et en érosion entre chaque nœud d’intersection. Le calcul de la « surface verticale » revient à mesurer la superficie de ces aires par un calcul d’intégrale. En ramenant ces superficies au mètre linéaire, on obtient des valeurs en m3/m.l. Des marges d’erreur inhérentes à la mesure tachéométrique, calculées par méthode statistique, sont comprises entre 5 et 7 cm en X et Y et restent inférieures à 3 cm en Z (Suanez et al., 2005).
11Nous séparerons les résultats acquis sur l'évolution du front de dune de ceux portant sur les profils de dune et du haut de plage.
12Pour le cordon dunaire du Grand Rocher, le bilan effectué au terme de la période d'observation montre un recul important du front de dune, la surface perdue atteignant 332 ± 5 m2 (fig. 5). Cette évolution s’inscrit dans une tendance au recul qui est observée depuis 1999 (Suanez, 2004). Toutefois, c’est principalement entre les mois d’octobre et de mars 2003 que le recul s’effectue. Durant cette période, le cordon dunaire perd près de 618 ± 5 m2, soit un quart de sa superficie initiale du mois d’octobre 2002 (2 725 ± 5 m2). Si l’érosion affecte l’ensemble de l'édifice, c’est toutefois la partie ouest qui reste la plus touchée. À certains endroits, la dune a totalement disparu découvrant ainsi l’enrochement qui est désormais battu par la houle. Dans le détail, les mesures montrent que le recul s’est effectué par à-coups, les mois de novembre 2002 et janvier 2003 enregistrant les pertes les plus importantes (respectivement -330 m2 et -185 m2). À partir du mois de mars 2003, la position du front de dune se stabilise dans un premier temps avant de connaître une phase de progradation importante durant l’été 2003. La surface gagnée durant cette période atteint plus de 365 ± 5 m2.
Fig. 5 – évolution du cordon dunaire du Grand Rocher entre le mois d'octobre 2002 et le mois de septembre 2003.
Fig. 5 – Dunes of Grand Rocher evolution between October 2002 and September 2003.
Evolution du front de dune du Grand Rocher entre (A) les mois d'octobre 2002 et de février 2003 et entre (B) les mois de février 2003 et de septembre 2003.
Dune of Grand Rocher evolution between (A) October 2002 and February 2003 ; between (B) February and September 2003.
Fig. 6 – évolution du cordon dunaire de Saint-Michel-en-Grève entre le mois d'octobre 2002 et le mois de septembre 2003.
Fig. 6 – Evolution of the coastal dunes of Saint-Michel-en-Grève between October 2002 and September 2003.
13À l’inverse du cordon dunaire du Grand Rocher, une grande stabilité du front de dune de Saint-Michel-en-Grève est observée (fig. 6). Entre les mois d’octobre 2002 et de mai 2003, les modifications restent mineures. L’évolution des superficies dunaires illustre ces changements peu importants. Au mois de mai 2003, la surface du cordon dunaire équivaut à celle du mois d’octobre 2002. Le levé effectué au mois de juin 2003 montre une augmentation très nette d’environ 1 450 ± 28 m2. Cette tendance se poursuit durant la période estivale car les mesures effectuées en septembre 2003 indiquent une nouvelle augmentation de la surface d’environ 1 776 ± 28 m2. Sur l’ensemble de la période, le cordon dunaire de Saint-Michel-en-Grève a gagné 3 363 ± 28 m2.
14L’évolution du système « cordon dunaire / haut de plage » du Grand Rocher montre une évolution contrastée entre les profils 1 et 2 (fig. 7). Les bilans obtenus sur l’ensemble de la période montrent un déficit très important au niveau du profil 1 situé à l’ouest (-4m3/m.l), alors qu’au niveau du profil 2 situé à l’est, le système s’est engraissé (1,5 m3/m.l).
Fig. 7 – évolution du bilan sédimentaire (m3/m.l) du cordon dunaire du Grand Rocher mesurée à partir des profils 1 et 2.
Fig. 7 – Evolution of the sediment budget (m3/l.m) of the dunes of Grand Rocher based on the surveys of profiles 1 and 2.
Fig. 8 – évolution du bilan sédimentaire (m3/m.l) du cordon dunaire de Saint-Michel-en-Grève mesurée à partir des profils 3, 4, 5 et 6.
Fig. 8 – Evolution of the sediment budget (m3/l.m) of the dunes of Saint-Michel-en-Grève based on the surveys of profiles 3, 4, 5 and 6.
15Dans le détail, on s’aperçoit que l’érosion a prédominé durant la période hivernale entre les mois d’octobre 2002 et de mars 2003. Une phase particulièrement érosive enregistrée au mois de janvier 2003 s’illustre par une déperdition de matériel sédimentaire atteignant respectivement -2,5 et -1,5 m3/m.l pour les profils 1 et 2 (fig. 7). À partir de cette période, l’ensemble du système connaît un engraissement continu, particulièrement net dans la partie orientale du cordon dunaire (profil 2). Au niveau du profil 1, on note encore une phase érosive bien marquée durant le printemps 2003.
16Le cordon dunaire de Saint-Michel-en-Grève connaît un engraissement sédimentaire continu tout au long de la période d’observation (fig. 8). Si quelques phases de démaigrissement apparaissent, à l’exemple du profil 3 à la date du 18 février 2003, la tendance observée au niveau des quatre profils montre que le cordon s’engraisse de façon importante. Le bilan sédimentaire global calculé à partir des quatre radiales reste largement positif. Il est compris entre +1 et +2,5 m3/m.l suivant les secteurs.
17Dans le but de cerner l’impact des forçages météorologiques sur l’évolution morphosédimentaire des cordons dunaires, une analyse des données météo-marines sur la période d'observation a été effectuée afin de distinguer, d’une part, les événements tempétueux intervenant dans la dégradation des dunes et, d’autre part, les périodes d’accalmie généralement propices à la régénération des systèmes dunaires. Quatre types de données ont été retenus : les données de vent tri-horaires (vitesse et direction) pour la station de Lannion, la pression barométrique (valeurs tri-horaires pour la station de Lannion), la marée enregistrée à Roscoff et des données de houle obtenues à partir de simulation numérique effectuée avec le logiciel TOMAWAC (fig. 9).
18Les données de houle pour l’ensemble de la période d’étude du 1er septembre 2002 au 30 septembre 2004 ont été acquises par simulation numérique à partir du modèle d’états de mer océanique du logiciel TOMAWAC dans le cadre du projet de construction d’un atlas numérique de houle le long des côtes Atlantique-Manche-Mer du Nord (collaboration entre Météo France, le Centre d’Etudes Techniques Maritimes et d’Études Fluviales (CETMEF), et LNHE-EDF). Les données de vent utilisées pour la modélisation proviennent de mesures satellitaires effectuées dans le cadre du projet ISLSCP (International Satellite Land Surface Climatology Project) NOAA-NCEP/Reanalysis2 (champs de vent à 10 m sur un pas de temps de 6 h et un pas d’espace de 1,875°). Les caractéristiques de houle ont été calculées en un point situé à l’entrée NW de la baie de Saint-Michel-en-Grève (coordonnées 3,7368° de longitude W et 48,9294° de latitude N). Trois paramètres réduits du spectre directionnel de houle ont été calculés toutes les 3 h : la hauteur significative spectrale Hmo en mètre, la direction moyenne de propagation ¦m en degré, la période de pic (Tpic) par la méthode de Read à l’ordre 5 en seconde.
19Les données obtenues sur l’ensemble de la période d’étude montrent des hauteurs maximales d’environ 9 m associées à des périodes de plus de 21 secondes. La hauteur moyenne est égale à 1,75 m pour une période de 10 secondes (tab. 2). La prise en compte des valeurs modales de H1/3 = 2 m et de T1/3 = 9 s, calculées à partir des mesures effectuées au large de Ouessant entre 1985 et 2001, et disponibles en libre accès sous la base de données Candhis du CETMEF (http://www.cetmef.equipement.gouv.fr/donnees/candhis/), permet de cerner la période de forte agitation qui s’étend du mois d’octobre 2002 au mois d’avril 2003 (fig. 9).
Tableau 2 – Caractéristiques morphodynamiques des secteurs du Grand Rocher et de Saint-Michel-en-Grève et calcul du run-up à partir de l’équation de P. Ruggiero et al. (2001).
Table 2 – Morphodynamic characteristics of the Grand Rocher and Saint-Michel-en-Grève zones and runup calculation using the P. Ruggiero et al.’s formula (2001).
Fig. 9 – Caractéristiques des conditions météo-marines enregistrées durant la période d’observation (octobre 2002 – septembre 2003).
Fig. 9 – Meteorological and marine conditions recorded during the period of survey (October 2002 – September 2003).
20À l’exception des oscillations du plan d’eau liées à la marée astronomique, trois éléments principaux interviennent dans l’élévation « temporaire » du niveau marin à la côte définissant ainsi la notion de niveaux extrêmes et jouant dans l’érosion des littoraux, que ce soient les côtes à falaise ou les littoraux meubles (Shih et Komar, 1994 ; Ruggeiro et al., 1997, 2001).
21Les surcotes, le plus souvent associées aux passages des perturbations atlantiques, s’expliquent par deux phénomènes (Bouligand et Pirazzoli, 1999 ; Pirazzoli, 2000 ; Betts et al., 2004 ; Pirazzoli et al., 2004) : une diminution de la pression barométrique (par définition < 1 015 hpa) et l’intervention de vents de mer violents poussant la masse océanique à la côte.
22L’élévation du niveau moyen de la mer liée à l’agitation, encore appelée le setup ou ηmax, (Komar, 1998). Cette déformation du plan d’eau forme un angle dont la pente orientée vers le rivage trouve son origine au niveau de la zone de déferlement.
23Enfin, les processus qui commandent la zone de swash interviennent également dans l’élévation du niveau marin au rivage. Cette dernière appelée run-up est définie comme étant la somme de deux composantes : les variations du plan d’eau liées à l’agitation (setup = ηmax) et l’action du jet de rive (wave uprush = R), soit RT = ηmax + R (Komar, 1998 ; Ruggiero et al., 2001 ; Masselink et Hughes, 2003).
24Dans le cas présent, l’analyse des niveaux extrêmes repose sur la quantification de deux paramètres, les surcotes et le run-up. Les surcotes ont été calculées par simple différence entre la marée observée et la marée prédite enregistrée par le marégraphe de Roscoff pour l’ensemble de la période considérée (fig. 9). Ces données ont été acquises auprès du SHOM (Brest) dans le cadre d’une convention de partenariat. Le run-up a été calculé à partir des équations issues des travaux de R.A. Holman (1986) et P. Nielsen et D.J. Hanslow (1991). Ces derniers insistent sur la relation entre l’élévation verticale maximale du run-up et les deux paramètres que sont la hauteur de la houle au large Hsig ou Hmo et le nombre d’Iribarren (surf similarity parameter) :
TTomomoξoT2%run-up25Plus récemment, les travaux de P. Ruggiero et al. (2001) portant sur l’analyse du run-up dans un contexte de plage de type dissipant (Agate beach, Oregon) ont montré que la relation entre la pente S, la longueur d’onde Lo et le run-up était bien plus significative. L’expression mathématique que ces auteurs proposent est la suivante :
ßomoo2π22pic26Dans le cadre de cette étude, l’équation 4 a été utilisée en tenant compte du caractère « hyper-dissipant » de notre zone d’étude. Les résultats obtenus sont présentés dans le tableau 2. Pour des pentes tan ß respectives de 0,012 et de 0,025 pour les secteurs de Grand Rocher et de Saint-Michel-en-Grève, les valeurs tri-horaires calculées sur l’ensemble de la période donnent un RT2% moyen égal à 0,48 m pour le secteur du Grand Rocher et de 0,69 m pour celui de Saint-Michel-en-Grève, les valeurs maxima atteignant respectivement 1,92 m et 2,77 m.
27Les données de vents (vitesse et direction) et de pressions barométriques proviennent des enregistrements tri-horaires Météo France de la station de Lannion. L’une des difficultés est de déterminer une vitesse de vent critique à partir de laquelle se produit un événement morphogène tempétueux. Les nombreuses études réalisées sur les tempêtes, en domaine méditerranéen ou atlantique, suggèrent des vitesses de vent très variées (Lemasson et Régnauld, 1997 ; Bruzzi, 1996, 1998 ; Audran, 1999 ; Clabaut et al., 2000 ; Meur-Ferrec et Ruz, 2002 ; Hénaff, 2004). Les seuils fixés par les auteurs dépendent de l’échelle spatio-temporelle à laquelle ils travaillent. Les travaux réalisés sur de grands domaines géographiques à des échelles de temps pluri-décennales ne retiennent que les événements très exceptionnels où les vitesses de vent sont supérieures à 15 m/s. À l’inverse, les études menées sur l’impact des événements épisodiques sur des secteurs plus restreints utilisent un seuil beaucoup plus bas, compris généralement entre 8 et 11 m/s. Dans le cadre de cette étude, nous considérons que, au-delà de 10 m/s, un événement potentiellement morphogène peut prendre place.
Tableau 3 – Caractéristiques hydrodynamiques et météo-marines des principaux épisodes morphogènes enregistrés durant la période d’observation.
Table 3 – Atmospheric and hydrodynamic characteristics of the main geomorphic events recorded during the period of survey.
28La figure 9 synthétise l’ensemble des données météo-marines enregistrées sur l’ensemble de la période d’étude. Entre le 7 octobre 2002 et le 15 avril 2003, on note une forte activité morphogène. Treize épisodes définis comme étant potentiellement morphogènes ont été recensés (tab. 3). Ils présentent des caractéristiques hydrodynamiques et météo-marines très variables. Certaines tempêtes, comme celles du 13 octobre 2002, 29 décembre 2002, 28 février 2003, 19 avril 2003 et du 2 mai 2003 coïncident avec de fortes marées (niveaux marégraphiques observés > 3,5 m) et/ou de fortes surcotes pouvant dépasser 30 cm. À l’inverse, les hauteurs de houle Hmo enregistrées durant ces événements ont été peu élevées (entre 3 et 4 m). Dans d’autres cas, des caractéristiques inverses sont observées. Les épisodes morphogènes du 27 octobre 2002, du 2 décembre 2002 et du 30 janvier 2003 sont marqués par des hauteurs de houles très importantes (Hmo compris entre 6 et 9 m) ; à l’inverse les niveaux d’eau enregistrés durant ces événements sont peu élevés (faible marée observée et/ou faible surcote). Notons enfin que, à l’exception des épisodes morphogènes (toutefois de faible intensité) du 19 avril et du 2 mai 2003, une longue période d’accalmie prend place dès le mois de mars 2003. Ces conditions climatiques « exceptionnelles » marquées par un régime anticyclonique de beau temps se sont accompagnées le plus souvent de vents de secteurs est à nord-est.
29L’analyse des forçages météo-marins et des dynamiques morphosédimentaires repose sur un traitement statistique permettant de corréler l’évolution diachronique des bilans sédimentaires aux conditions de forçage météo-marines enregistrées entre chaque campagne de mesure. À cet effet, l’Analyse factorielle en Composante Principale (ACP) a été utilisée. L’ACP mesure et représente séparément les relations entre les individus et celles entre les variables, sous la forme de deux nuages de points, dans un espace à n-1 dimensions défini par n-1 axes factoriels (Chamussy et al., 1994). Les fortes relations s’expriment selon la distance entre les points qui, dans le cas de l’ACP, est une distance euclidienne. Les axes factoriels sont hiérarchisés de telle sorte qu’ils comportent une part décroissante d’information, encore appelée « inertie ». On considère qu’au-delà du deuxième axe factoriel, la part d’information retenue par le modèle n’est plus significative.
30Dans le cas qui nous occupe, deux dynamiques morphosédimentaires sont analysées : les phases d’érosion et de régénération des cordons dunaires. Pour les phases d’érosion, nous recherchons les relations les plus fortes qui permettent d’identifier le rôle des différents agents météo-marins dans la diminution du bilan sédimentaire des cordons dunaires : les superficies dunaires (SUPDU1 pour le cordon dunaire du Grand Rocher ; SUPDU2 pour celui de Saint-Michel-en-Grève) et les surfaces verticales calculées sur les six profils (SUPROF1, SUPROF2, SUPROF3, etc.). Notre attention a été portée sur l’importance des niveaux extrêmes comme facteur explicatif de l’érosion des dunes. Ainsi, cinq variables ont été retenues, la hauteur de la houle (Hmo), la marée observée (MAROBS), qui fait la somme du niveau marin donné par la marée prédite et de la surcote, le run-up (RUNUP1 pour le cordon dunaire du Grand Rocher ; RUNUP2 pour celui de Saint-Michel-en-Grève), la vitesse du vent (VITV) et la pression barométrique (PRESS).
31Les valeurs ayant été retenues pour chaque variable correspondent à celles qui caractérisent l’épisode le plus morphogène enregistré entre chaque campagne de mesure. Cette discrimination comporte une part d’incertitude car, comme nous l’avons vu précédemment, deux à trois épisodes tempétueux ont pu avoir lieu entre deux relevés topo-morphologiques, comme cela s’est produit durant les périodes d’octobre à novembre, de janvier à février et d’avril à mai (tab. 3). Dans tous les cas, notre choix s’est fait en fonction de deux paramètres essentiels, les plus importantes hauteurs d’eau de marée observée et/ou les plus fortes hauteurs de houle. Le tableau 4 consigne les épisodes morphogènes retenus. Pour la période allant du mois de mai à septembre 2003, marquée par une absence de tempête, nous avons sélectionné « l’épisode » caractérisé par une hauteur de houle et/ou un niveau marégraphique le(s) plus important(s).
Tableau 4 – Dates, heures et caractéristiques des épisodes morphogènes retenus pour l’analyse factorielle en composante principale.
Table 4 – Dates, hours, and characteristics of the geomorphic storm events selected for the PCA statistical analysis (Principal Component Analysis).
Tableau 5 – Calcul théorique du transit éolien orienté vers les cordons dunaires durant la période d’observation à partir des formules de R.A. Bagnolds (Eq. 6) et de R. Kawamura (Eq. 7). Total des précipitations (mm) tombées dans l’intervalle de temps et précipitations (mm) ramenées à la journée pour chaque période.
Table 5 – Theoretical calculation of eolian sand transport to the dunes during the period of survey using R.A. Bagnolds’s (Eq.6) and R. Kawamura’s formulas (Eq.7). Total rainfall recorded for each interval of time and rainfall per day for each period.
32Pour l’analyse des phases de régénération des cordons dunaires, nous avons privilégié une approche théorique en travaillant sur les périodes durant lesquelles les vents de mer ont pu générer un engraissement du système haut de plage / cordon dunaire. Une quantification des volumes de sables transportés par le vent a été réalisée pour chaque intervalle de temps séparant deux levés (tab. 5). Pour ce faire, une discrimination des vents morphogènes a été effectuée à partir de deux filtres : la vitesse d’arrachement des sables et la direction qui est fonction de l’orientation et de la localisation des cordons dunaires : 280° à 80° nord pour le cordon dunaire du Grand Rocher et 220° à 20° nord pour celui de Saint-Michel-en-Grève. Les vitesses d’arrachement des sables ont été retenues dès lors que la vitesse du vent mesurée V* > V*c (vitesse critique d’arrachement des sables). Dans le cas présent, la vitesse critique d’arrachement V*c a été calculée à partir de la formule de R.A. Bagnold (1941) :
33V*c = k1 [(ρs / ρa) (gd)]0,5 (5)
où k1 est une constante = 0,1 ; ρs : densité des sables (2,7 g/cm3) ; ρa : densité de l’air (0,0023 g/cm3) ; d : diamètre du grain moyen en cm ; g : accélération de la pesanteur (9,81 m/s2). Pour un grain moyen de 155 μm, constituant les sables vifs du système haut de plage / dune, la vitesse critique d’arrachement V*c est égale à 5,7 m/s. Un seuil de 6 m/s a été choisi pour filtrer les données brutes.
34Les formules de quantification ayant été retenues sont celles de R.A. Bagnolds (1941) et de R. Kawamura (1964). Ces équations sont les plus fréquemment utilisées dans les études portant sur le transit sédimentaire éolien (Svasek et Terwindt, 1974 ; Sarre, 1988 ; Arens, 1994 ; Meur-Ferrec et Ruz, 2002 ; Alcantara-Carrio et Alonso, 2002) :
35équation de R.A. Bagnolds (1941)
Q = C (d/D)0,5 V*c3 (ρa/g) (6)
36équation de R. Kawamura, (1964)
Q = K(ρa/g) (V*-V*c) (V*+V*c)2 (7)
où Q : quantité de sable transporté (kg/m/s) ; V* : vitesse du vent mesuré en cm/s ; V*c : vitesse critique d’arrachement en cm/s ; C : coefficient sans dimension qui dépend du tri du sédiment = 1,5 très bon tri, 1,8 bon tri, 2,8 mauvais tri, 3,5 éléments non mobilisables ; K : constante de 2,7 pour des sables compris entre 100 et 800 μm ; ρa : densité de l’air (1,23 kg/m3) ; g : accélération de la pesanteur (9,81 m/s2) ; d : le diamètre moyen des grains de sables (en cm) et D : valeur du grain standard (0,025 cm).
37La figure 10 montre une bonne relation entre les fortes houles de tempête, les hauts niveaux de submersion qui leur sont associés (run-up) et l’érosion des cordons dunaires. Ces phénomènes sont visibles sur le terrain (fig. 11). Cette relation structure pleinement le premier facteur se rapportant à l’analyse des superficies dunaires (fig. 10A). Elle illustre également la forte saisonnalité du fonctionnement morphosédimentaire en opposant la période allant de la fin du mois d’octobre au début du mois de décembre 2002, marquée, comme on l’a vu, par des tempêtes fortement morphogènes à celle allant de la mi-mai au début du mois de septembre 2003, caractérisée en revanche, par des conditions météo-marines très calmes.
Fig. 10 – Analyse statistique en composantes principales visant à corréler le bilan sédimentaire des cordons dunaires et les conditions météo-marines enregistrées pendant la période d’observation.
Fig. 10 – Statistical analysis (PCA) of dune sediment budget and meteo-marine conditions recorded during the monitoring period.
Fig. 11 – Effets du run-up et submersion du cordon dunaire de Saint-Michel-en-Grève (cliché : 15/11/ 2002, S. Suanez).
Fig. 11 – Effects of run-up on the submergence of the coastal dune of Saint-Michel-en-Grève dune (photo : 15/11/2002, S. Suanez).
Fig. 12 – Correspondance entre la marée observée, les niveaux d’eau extrêmes (marée observée + run-up) et l’altitude du pied de dune pour les six profils, sur l’ensemble de la période d’observation.
Fig. 12 – Correspondence between observed tide, high water level (observed tide water level + run-up) and the dune toe height for six profiles over the whole survey period.
A : 1 : pleine mer de marée observée (marégraphe de Roscoff) ; 2 : pleine mer de marée observée + run-up pour le secteur de Saint-Michel-en-Grève ; 3 : altitude du pied de dune au niveau du profil 3 ; 4 : altitude du pied de dune au niveau du profil 4 ; 5 : altitude du pied de dune au niveau du profil 5 ; 6 : altitude du pied de dune au niveau du profil 6 ; B : 1 : pleine mer de marée observée (marégraphe de Roscoff) ; 2 : pleine mer de marée observée + run-up pour le secteur du Grand Rocher ; 3 : altitude du pied de dune au niveau du profil 1 ; 4 : altitude du pied de dune au niveau du profil 2.
A : 1 : measured high tide level (Roscoff maregraph) ; 2 : measured high tide level + run-up at Saint-Michel-en-Grève site ; 3 : height of the dune toe for profile 3 ; 4 : height of the dune toe for profile 4 ; 5 : height of the dune toe for profile 5 ; 6 : height of the dune toe for profile 6 ; B : 1 : measured high tide level (Roscoff maregraph) ; 2 : measured high tide level + run-up at Grand Rocher site ; 3 : height of the dune toe for profile 1 ; 4 : height of the dune toe for profile 2.
38L’analyse portant sur les profils montre des résultats moins tranchés. Si ce schéma semble à peu près fonctionner pour le cordon dunaire du Grand Rocher (fig. 10B), aucune relation allant clairement dans ce sens ne s’exprime pour le cordon dunaire de Saint-Michel-en-Grève (fig. 10C). Cela s’explique par le fait que l’érosion induite par les phénomènes de run-up n’affecte pas l’ensemble du cordon dunaire mais uniquement certaines parties de ce dernier. Comme le montre la figure 12, l’altitude du pied de dune conditionne l’importance des phénomènes érosifs dès lors qu’elle est inférieure à la hauteur du niveau d’eau extrême (marée observée + run-up). Sur l’ensemble de la période d’observation, ces hauts niveaux marins ont atteint le pied de dune du Grand Rocher bien plus fréquemment que celui de Saint-Michel-en-Grève (fig. 12A et 12B). Comme on l’a vu précédemment, cela s’est traduit par un engraissement continu de ce dernier tout au long de la période d’observation. On note que les hauts niveaux marins correspondant à la marée observée stricto sensu ne jouent pas dans l’érosion des dunes, ce qu’illustre la figure 12. Dans la majorité des cas, cette variable structure le deuxième axe factoriel mais en tant qu’élément isolé. Cela montre que les niveaux de pleine mer de marée observée sont inférieurs à l’altitude des pieds de dunes et que, comme l’a indiqué P.A. Pirazzoli (2000), les plus fortes marées de vive-eau ont rarement coïncidé avec un épisode tempéteux de très forte intensité.
39Enfin, comme on l’a vu précédemment, la saisonnalité s’exprime particulièrement bien au niveau des périodes 2 et 3 (fin du mois d’octobre au début du mois de décembre 2002) qui, dans tous les cas, s’opposent fortement aux périodes 10 et 11 (mi-mai au début du mois de septembre 2003) (fig. 13). Pour le reste des individus on note une forte dispersion qui brouille ce fonctionnement saisonnier. Cela peut s’expliquer par l’ensemble des processus agissant durant les périodes de régénération des cordons dunaires qui, comme on l’a dit précédemment, peuvent gommer l’effet des tempêtes. L’estimation du transit éolien, couplé aux précipitations, illustre bien ces mécanismes (tab. 5). Durant les périodes 1 et 2 marquées par de fortes tempêtes, les valeurs théoriques des volumes sédimentaires transportés vers les cordons dunaires sont considérables. Toutefois, si l’on tient compte des précipitations importantes enregistrées durant ces périodes (> 1 mm/jour), on peut penser que les déplacements sédimentaires ont été quasiment nuls, en raison du rôle limitant du taux d’humidité dans le transport éolien (Svasek et Terwindt, 1974 ; Sarre, 1988 ; Arens, 1996b ; Meur-Ferrec et Ruz, 2002). De ce fait, la régénération des cordons dunaires durant ces périodes a été inexistante, d’où la forte signature des épisodes morphogènes sur la morphologie. À l’inverse, pour les autres périodes hivernales et printanières (périodes 4, 6, 7 et 8), l’impact morpho-sédimentaire des tempêtes a été en partie gommé par un engraissement éolien plus ou moins important favorisé par un temps sec.
Fig. 13 – Exemples de morphologie du haut de plage et du front dunaire, dune du Grand Rocher.
Fig. 13 – Examples of the morphology of the upper beach and dune front, Grand Rocher dune.
A : présence d’une falaise d’érosion (cliché : 23/10/2002, S. Suanez). B : accumulation en pied de dune et végétation pionnière (cliché : 12/09/2003, S. Suanez).
A : erosional dune scarp (photo : 23/10/2002, S. Suanez). B : sand accumulation on the dune toe and pioneer plants (photo : 12/09/2003, S. Suanez).
40D’une façon générale, les résultats obtenus sont en accord avec le schéma classique du fonctionnement saisonnier des cordons dunaires décrit dans la littérature (Davidson-Arnott et Law, 1990 ; Sherman et Bauer, 1993 ; Carter et al., 1990 ; Hesp, 1981 ; Hallégouët, 1981). À la période automnale-hivernale de forte activité morphogène marquée par un stade érosif succède une phase de régénération, particulièrement efficace durant l’été. Toutefois, lorsque l’on compare le fonctionnement des deux secteurs étudiés, des différences apparaissent. Elles rendent compte de spécificités morphosédimentaires et hydrodynamiques propres à chacun des sites. Ces dernières permettent de replacer le fonctionnement de ces deux cordons dunaires au sein des modèles conceptuels (Short et Hesp, 1982 ; Psuty, 1992 ; Hesp, 2002 ; Ruz, 2004).
41Si l’on se fonde sur le modèle de A.D. Short et P.A. Hesp (1982) reposant sur une classification morphodynamique des environnements côtiers, les cordons dunaires du Grand Rocher et de Saint-Michel-en-Grève devraient être très étendus et bénéficier d’une alimentation sédimentaire très importante. En suivant ce schéma, d’une part le caractère « hyper-dissipant » de la baie de Saint-Michel-en-Grève minimiserait les phénomènes érosifs par un fort amortissement de la houle, d’autre part, la grande largeur de l’estran très plan expliquerait un transport éolien de la plage vers la dune très important. Si pour le cordon dunaire de Saint-Michel-en-Grève ce modèle peut s’appliquer, pour celui du Grand Rocher marqué par une forte érosion, le fonctionnement apparaît bien plus complexe.
Fig. 14 – Modèle de formation et d’évolution des cordons dunaires fondé sur le bilan sédimentaire de la plage et de la dune (d’après N.P. Psuty, 1992).
Fig. 14 – Model of construction and evolution of coastal dunes based on beach and dune sediment budget (from Psuty, 1992).
42Le modèle proposé par N.P. Psuty (1992) reposant sur la notion de bilan sédimentaire du couple plage / dune semble plus approprié (fig. 14). Ainsi, on peut classer le cordon dunaire de Saint-Michel-en-Grève dans le scénario (B) décrivant une situation où le bilan sédimentaire de la plage et de la dune est positif. Dans ce contexte, l’auteur souligne l’existence de conditions propices à l’édification de dunes bordières successives au fur et à mesure de la progression du rivage. Cette disposition rappelle à bien des égards la morphologie décrite précédemment. À l’inverse, le cordon dunaire du Grand Rocher se situerait plutôt entre les scénarios (D) et (C) décrivant respectivement des situations où le bilan sédimentaire du couple plage / dune est déficitaire, déficitaire pour la plage et équilibré à positif pour la dune. Dans ce contexte on assiste à une dégradation du cordon dunaire liée à l’action des phénomènes de submersion (onde de tempête et/ou PM de vive-eau) avec, durant les périodes de faible activité morphogène, des courtes phases de régénération. Ces notions portant sur la relation entre le bilan sédimentaire et la dynamique des cordons dunaires ont en partie été analysées et décrites par S. Suanez (2004).
43Enfin, le modèle proposé par P.A. Hesp (1988, 2002) reposant sur une classification morpho-écologique permet de faire la synthèse entre les dynamiques agissant à court terme et l’évolution à plus long terme des cordons dunaires. Ainsi, le cordon dunaire de Saint-Michel-en-Grève appartiendrait au stade 2 correspondant à une dune bordière bien constituée et bien végétalisée (70 à 90 % de couverture végétale), pouvant être précédée de dunes embryonnaires. Dans ce contexte, des formes mineures d’érosion (siffle-vent) peuvent affecter la morphologie du versant externe ; toutefois, sur le long terme, ce cordon évoluerait vers un stade 2A. Ce dernier caractérise les côtes stables ou en accumulation et se distingue par un développement de la dune par exhaussement et/ou formation de crêtes dunaires successives. Sur du court terme, on peut avoir un stade 2D décrivant une érosion en micro-falaise par les vagues de tempête. Néanmoins, pour le cordon dunaire de Saint-Michel-en-Grève, la réversibilité de ces phénomènes s’illustre par un ensevelissement rapide de ces formes d’érosion et une revégétalisation du versant externe.
44Le cordon dunaire du Grand Rocher reste difficilement classable tant les processus de recul qui l’affectent sont rapides. Ainsi, les cycles d’érosion décrits par l’auteur ont du mal à s’appliquer dans ce cas. Du point de vue morpho-écologique et de la dynamique, il appartiendrait à un stade intermédiaire entre les stades 2B et 3B, stade marqué par un recul continu sur le long terme. Sur le court terme, la dynamique serait liée à l’action incessante de l’érosion entaillant en micro-falaise le versant externe qui jamais ne retrouverait sa position initiale malgré des phases temporaires de réalimentation.
45Le suivi à haute fréquence réalisé dans le cadre de cette étude a permis d’appréhender de façon satisfaisante les changements morphosédimentaires rapides des cordons dunaires. L’expérience montre que les deux méthodes de mesure mises en œuvre apparaissent complémentaires et utiles. Si la mesure le long des profils permet de cerner les modifications sur l’ensemble du système plage / dune, elle reste toutefois très dépendante de la localisation du profil. Cela introduit une part importante de subjectivité dans les résultats obtenus qui sont alors difficilement généralisables à l’ensemble du cordon dunaire. La mesure du front de dune permet d’éliminer cette incertitude et de donner des résultats plus pertinents à partir de l’évolution des surfaces dunaires. Toutefois, elle ne peut s’appliquer que pour des cordons dunaires présentant une morphologie assez simple et dont la limite frontale est facilement identifiable, ce qui est le cas des cordons dunaires de la baie de Saint-Michel-en-Grève.
46Le recours à l’Analyse factorielle en Composante Principale (ACP) apparaît efficace et relativement simple à mettre en œuvre. Elle a permis d’une part de montrer le rôle important que pourrait jouer les phénomènes de run-up associés aux fortes houles de tempêtes dans le recul du cordon dunaire du Grand Rocher, confirmant ainsi la grande fragilité de ce secteur. D’autre part, elle souligne le rôle limité des hauts niveaux de marée de vive-eau qui, dans la plupart des cas, n’atteignent jamais le pied de dune.
47Enfin, l’analyse des forçages météo-marins et des dynamiques morphosédimentaires montre un fonctionnement saisonnier bien marqué entre l’hiver et l’été. Toutefois, les mécanismes associés aux courtes périodes de faible activité morphogène venant s’intercaler aux épisodes tempétueux brouillent quelque peu ce schéma. Comme l’ont montré plusieurs auteurs, cela souligne les interactions entre les processus agissant à différentes échelles temporelles (Davidson-Arnott, 1988 ; Thom et Hall, 1991 ; Battiau-Queney et al., 2000 ; Battiau-Queney, 2004), c'est-à-dire saisonnière et épisodique dans le cas présent.