Abridged English version
The end of the Norwegian Margin formation during the opening of the North Atlantic Ocean coincides with the Palaeocene–Eocene transition, about 55 Ma ago [6]. From the Miocene to the Pliocene times, the fine-grained hemipelagic sediments of the Kai Formation are deposited on the margin. The Kai Formation is affected by a polygonal fault system in relation with fluid expulsion [1]. From the Pliocene to the Pleistocene times, the thick prism of sediments of the Naust Formation is deposited on the shelf during several glacial episodes [14]. From 2.5 to 1.3 Ma, the prograding platform fills the Storegga depression. During the last 1.3 Ma, the sedimentation is controlled by glacial and interglacial cycles [11], with a maximum activity at 0.6–0.5 Ma when the Fenno-Scandian ice sheet repeatedly covers the area. The Storegga slides have been mapped over more than 800 km, from the head wall, on the top of the continental slope at about 400 m water depth down to the Norwegian basin abyssal plain at 3600-m water depth [4]. Eight retrogressive slides, triggered over a short time period, have been identified [8]. Whereas the main phase of the Storegga slides was dated at about 7250–7300 years 14C (8150–8200 bp), recent studies on the northern flank of the slides show evidence for more recent sliding episodes: 5000 ± 300 years 14C et 2200–2500 years 14C. Gas hydrates are composed of gas molecules trapped within water molecules cages [12]. Gas hydrates formation is mainly controlled by thermodynamical conditions (high pressure and low temperature) that occur on most of the continental margins at water depths between 500 and 1000 m. The occurrence of marine gas hydrates in the sediments is mainly inferred from the observation on seismic sections of a Bottom Simulating Reflector (BSR). Geophysical studies [9] have shown that the BSR is generated at the transition between sediments containing a variable amount of solid gas hydrate above and sediments containing a small volumetric fraction of free gas below. The stability field for gas hydrate being mainly -dependent, and isotherms approximately parallel to the sea floor, the BSR parallels the seabed, and shows a polarity opposite to the seafloor one and often cuts across sedimentary layer reflections. In the Storegga area, BSRs have been identified north of the slides and extend below the slides in the eastern part of the area [2,3,5]. However, a geotechnical boring drilled over a well-defined BSR area [10] has shown no evidence for gas or gas hydrates.
The HYDRATECH cruise on N/O Le Suroît (June–July 2002) aimed at a high-resolution survey of a limited area of the Storegga slides. Data collected during the cruise include bathymetric and reflectivity EM300 maps (Fig. 1), chirp echo-sounder records (Fig. 3c), and high-resolution (HR) 2D (Fig. 2) and 3D (Fig. 4) seismic data [13]. The study area (Fig. 1), about 3200 km2, extends from the upper slope at the northern boundary of the slides to within the slides, in over 2-km water depth. It crosses slide scars and zones of potential present-day mechanical instability on the northern flank of the slides. It covers fields of pockmarks and various intra-sedimentary structures interpreted as fluid escape structures. Seismic data interpretation confirms that, as proposed in Bouriak et al. [2] and Bünz et al. [5], the existence, location and amplitude of the BSR is controlled by the sediments properties. The main controlling factor would be the sediments ability to favour circulation and trapping of free gas below the base of the gas hydrate stability zone (BGHS). The BSR could then mainly be due to free gas trapped below the BGHS, without any indication of high gas-hydrate concentrations above, in agreement with the low gas-hydrate concentration inferred from refraction seismic data analysis [5].
High-resolution seismic lines point to intense small-scale fracturing and diffuse degassing through the near seafloor sediments in the slides (Fig. 2). In contrast, the slope area immediately north of the edge of the Storegga slides is nearly void of any seismic evidence pointing to seafloor degassing structures. Further upslope, fields of numerous gas chimney-type structures are observed starting at a distance of 10–20 km from the slides, in water depth of less than 1 km (Fig. 1). High-resolution images obtained with a 2–5-kHz chirp profiler illustrate various types of these structures (Fig. 3c). Not all of the chimneys are observed to reach the seafloor. When they do, the chimneys are often associated with seafloor depressions or pockmarks. Bathymetry and acoustic seafloor imagery data (Fig. 3a and b) show that these pockmarks can be up to 500 m large in diameter, and are usually associated with a higher seabed reflectivity. They are shallow, rarely exceeding a few meters in depth. Pull-downs are commonly observed in and below the chimneys, which suggests the presence of gas. In several cases, surrounding reflectors are bent upwards near the chimneys, which may result from sediment deformation or an increase in the velocity in the surrounding sediments. Most of the chimney-like structures are rooted below the BSR. On seismic and chirp profiles, amplitude anomalies are observed in the undisturbed sediments above the chimneys that do not reach the seafloor. These anomalies affect the sediments up to the seafloor and suggest recent fluid expulsion activity of these chimneys.
Bathymetric data (Fig. 1) point to preferential directions in the slide scars and suggest a tectonic control of the sliding process. On the 3D seismic block (Fig. 4), two shallow disturbed areas were identified located at about 150 m below seabed. Those slides affect thin (about 15 m for the upper slide) layers in the Naust sedimentary formation, extend over a large part of the study area and are sealed by posterior sedimentation, which implies early sliding of the sediments after deposition. The upper slide affects layers dated approximately 120 000 years BP. There is no clear link between these slides and gas hydrates' occurrence in the sediments: gas hydrate dissociation that could have had an impact on the strength of the sediments is expected to take place preferentially at the BGHS. However, the slides are nowadays still located well above the BGHS. The sliding layers are identified as glide planes for the −8200-year main Storegga slide event [7], which means that the instability of these layers could have lasted since the last 120 000 years. In places where these layers are undisturbed, the present-day stability is thus still in question.
1 Introduction
La zone de Storegga est caractérisée par l'existence d'hydrates de gaz dans les sédiments de la pente continentale et la présence reconnue de grands glissements. C'est donc une région privilégiée pour l'étude des interactions entre présence de fluides et d'hydrates de gaz et instabilités de pente. La déstabilisation des hydrates de gaz, lors du réchauffement qui a suivi la dernière glaciation, est une hypothèse avancée pour expliquer le déclenchement des glissements de Storegga. Réciproquement, se pose la question de la quantité de méthane libérée vers l'atmosphère lors du rétrécissement de la zone de stabilité des hydrates à la suite des grands glissements de Storegga. Pour apporter des éléments de réponse à ces questions, la campagne Hydratech se proposait d'établir une image sédimentaire précise, à l'aide d'outils géophysiques à haute définition, d'un secteur du flanc nord des glissements de Storegga. Il s'agissait de mettre en lumière la distribution du gaz et des hydrates de gaz dans les sédiments, les différentes structures d'échappement de fluides, les déformations intrasédimentaires, ainsi que les liens éventuels entre ces différents phénomènes. À travers la présentation d'extraits de données acquises sur la zone, cette contribution illustre l'intérêt de combiner différents outils géophysiques dont les données se complètent pour rendre compte des phénomènes physiques différents qui affectent la pente continentale.
2 La zone de Storegga
2.1 Contexte géologique
La fin de la formation de la marge de Norvège lors de l'ouverture de l'Atlantique nord coïncide avec la transition Paléocène–Éocène il y a environ 55 Ma [6]. Du Miocène au Pliocène, des sédiments fins hémipélagiques se déposent sur la marge : les séries de la formation Kai. Cette formation est affectée par un système de failles polygonales dont l'origine est liée à des expulsions de fluides [1]. Du Pliocène au Pléistocène, la marge subit plusieurs épisodes glaciaires, qui déposent sur le plateau un prisme épais de sédiments : la formation Naust [14]. De 2,5 à 1,3 Ma, la progradation de la plate-forme remplit la dépression de Storegga. Lors des 1,3 derniers millions d'années, la sédimentation est conditionnée par l'alternance de cycles glaciaires et interglaciaires [11], avec une activité maximale vers 0,6–0,5 Ma, lorsque la calotte glaciaire fenno-scandienne couvre la zone à répétition.
2.2 Les glissements de Storegga
Les glissements de Storegga, qui affectent les séries actuelles, ont été cartographiés sur plus de 800 km, depuis la cicatrice en tête des glissements sur la partie haute de la pente continentale, vers 400 m de profondeur, jusqu'à la limite des sédiments déplacés dans la plaine abyssale du bassin norvégien, à des profondeurs de 3600 m [4]. À ce jour, huit glissements récents, qui auraient fonctionné de façon rétrogressive et sur une courte durée (quelques heures), ont été identifiés [8]. Alors que le glissement principal de Storegga a été daté de 7250–7300 ans 14C (8150–8200 BP), des datations récentes sur le flanc nord font état de glissements à 5000 ± 300 ans 14C et 2200–2500 ans 14C, ce qui pose la question du risque actuel d'instabilité de cette zone [8].
2.3 Hydrates de gaz et BSR
Les hydrates de gaz sont composés de molécules d'eau organisées en cages, qui piègent des molécules de gaz [12]. Leur formation est en particulier contrôlée par des conditions thermodynamiques particulières, qui sont atteintes sur la plupart des marges continentales, à des profondeurs d'eau comprises entre 500 et 1000 m : fortes pressions, faibles températures. L'observation d'un Bottom Simulating Reflector (BSR) est aujourd'hui encore l'indice géophysique principal pour conclure à la présence d'hydrates de gaz. Ces réflecteurs sismiques, à la limite inférieure de la zone de stabilité des hydrates de gaz [9], suivent les ondulations topographiques du fond de mer, en recoupant les réflecteurs sédimentaires stratigraphiques. Ils sont caractérisés par une polarité négative, qui exprime, suivant les résultats de la plupart des études entreprises, la présence de gaz libre sous les sédiments contenant les hydrates.
Sur la zone de Storegga, des BSRs ont été observés sur des profils sismiques au nord des glissements. Ils s'étendent à l'est sous les glissements [2,3,5]. Un forage géotechnique a été réalisé dans une zone à BSR [10], mais ni hydrates, ni gaz libre n'ont été échantillonnés.
3 La campagne Hydratech : outils utilisés et données
Trois outils géophysiques principaux ont été utilisés. Ils ont permis d'obtenir des données complémentaires de surface et de sub-surface sur la zone.
3.1 Bathymétrie – Imagerie
Le sondeur multifaisceaux EM300 du N/O Le Suroît a permis de relever une carte bathymétrique (Fig. 1) et une carte de réflectivité acoustique du fond sur l'ensemble de la zone d'étude ().
3.2 Imagerie sismique haute résolution
Un outil de sismique haute résolution a été mis en œuvre pour, d'une part, acquérir des profils régionaux (sismique 2D conventionnelle) et, d'autre part, et pour la première fois dans le milieu universitaire, acquérir des données de sismique 3D couvrant une surface limitée.
3.2.1 Sismique haute résolution 2D
Une série de 24 profils haute résolution, enregistrés sur une flûte numérique de 750 m (120 traces de 6,25 m) immergée à 3 m, a été acquise. La source était composée de 2 canons à air mini-GI guns immergés à 1,5 m, qui permettent d'émettre dans une gamme de fréquences de 20 à 300 Hz. Après application d'un filtre passe bande et correction de l'atténuation géométrique, les données ont été sommées et migrées en profondeur, en prenant en compte la navigation réelle des sources et des récepteurs (Fig. 2).
3.2.2 Sismique haute résolution 3D (HR 3D)
Un système d'acquisition original développé à l'Ifremer [13] a été mis en œuvre, pour la première fois, dans le cadre d'une campagne scientifique. Le dispositif comprend deux flûtes de 48 traces chacune (intertrace de 6,25 m) et deux sources fonctionnant alternativement. Quatre lignes « point milieu » sont éclairées par ligne naviguée, d'où un gain de temps d'acquisition considérable par rapport à un dispositif conventionnel composé d'une source et d'une flûte. L'intertir est d'environ 6,25 m. Les lignes naviguées sont espacées de 25 m, ce qui a permis de couvrir une zone de en 12 jours environ. Le bloc de données 3D obtenu présente une résolution horizontale théorique de (dimensions du bin) et une résolution verticale de l'ordre du mètre. La qualité de l'image 3D finale dépend en grande partie du traitement soigné des données de navigation. La phase de prétraitement des collections points de tir comprend, en particulier, filtrage fréquentiel et corrections statiques. L'ensemble des données est ensuite migré en profondeur (migration avant sommation de type Kirchhoff). Le bloc migré présenté a été obtenu en utilisant un gradient de vitesse dans les sédiments, estimé à l'aide d'analyses de vitesse sur les lignes HR2D régionales. Les données sismiques ont été interprétées à l'aide du logiciel Sismage (Total).
3.3 Sondeur chirp
Le Suroît est équipé d'un sondeur de sédiments de coque émettant un signal modulé en fréquence sur une bande de fréquences assez large (2,5–4,5 kHz). Il a été utilisé en mode chirp, avec un signal d'une durée de 10 ms, à 25 % de la puissance maximale du sondeur, pendant 12 h par jour, sur toute la durée de la campagne (24 jours sur zone). Les données sont de bonne qualité, avec une pénétration moyenne de 150 à 200 ms temps double (Fig. 3).
4 Données et interprétations préliminaires
4.1 Hydrates de gaz et Bottom Simulating Reflector (BSR)
Le profil htk2D01–1, présenté sur la Fig. 2, illustre les caractéristiques principales du BSR sur la zone d'étude. On peut noter, en particulier, les propriétés variables de ce réflecteur. Il est continu sur le haut de la pente entre les points milieu commun (CDP) 1000 à 4000. À cet endroit, il est matérialisé par un réflecteur de polarité inverse par rapport au fond de l'eau et presque coïncident avec la stratigraphie. Le BSR disparaît alors (CDP 4500 et 6500) en tant que réflecteur. La base du domaine de stabilité des hydrates de gaz (notée dans la suite BHSZ pour Bottom of Hydrate Stability Zone) est définie par la terminaison d'anomalies d'amplitude. Ensuite, et jusqu'à environ 2 km au sud-ouest de la limite nord des glissements de Storegga, la BHSZ est alternativement marquée par des segments de réflecteur ou des terminaisons d'anomalies d'amplitude (CDP 6500 à 9000). Enfin, sous la zone affectée par les glissements (CDP 9000 à 13 000), anomalies d'amplitude et BSR ne sont pas visibles. Il faut noter que la profondeur du BSR en limite des glissements coïncide avec la base des niveaux sédimentaires affectés par les glissements. Une carte de la répartition du BSR sur la zone d'étude 3D a été réalisée à partir des données de la campagne. Les zones où le BSR apparaît en tant que réflecteur continu à la résolution de la sismique utilisée et celles où il se matérialise par la terminaison d'anomalies d'amplitude ont été cartographiées. Cette carte confirme le caractère discontinu du BSR. En particulier, l'amplitude relative du BSR dépend de la strate intersectée par celui-ci. Ces variations d'amplitude pourraient être dues à des quantités variables de gaz piégées sous le BSR, en fonction de la porosité et de la perméabilité des séries sédimentaires recoupées par le BSR. On note aussi, dans le bas de la zone 3D, un changement de l'amplitude relative du BSR, lié à la présence d'une faille. En première analyse, la présence et l'amplitude du BSR sont donc contrôlées par la nature des séries sous-jacentes, ainsi que par l'existence de failles intra-sédimentaires. Il semble que ce soit principalement la capacité des séries sédimentaires à piéger le gaz au niveau de la BHSZ qui contrôle ces amplitudes. Aucune anomalie d'amplitude qui pourrait être liée à la présence d'hydrate dans les séries au-dessus du BSR n'est observée.
4.2 Sorties de fluides
Les cartes bathymétriques (Figs. 1 et 3b) et d'imagerie (Fig. 3a) mettent en évidence de nombreuses structures sub-circulaires, d'un diamètre moyen de 300 m environ, et qui affectent le fond de l'eau en haut de pente de la zone d'étude. La plupart de ces structures sont des dépressions (pockmarks) sur la bathymétrie et correspondent à des zones de forte réflectivité sur l'imagerie. On peut observer des structures en dômes, mais celles-ci sont plus rares. Ces structures sont corrélées sur les profils de sondeur de sédiment et sur les profils sismiques (Figs. 2–4) à des zones verticales (cheminées), le long desquelles les réflecteurs semblent interrompus. Certaines de ces cheminées atteignent le fond de l'eau, d'autres semblent être scellées. Une carte de la répartition des cheminées, réalisée à partir des données sismiques (Fig. 1) montre qu'elles se développent au-dessus d'une profondeur d'environ 950–900 m à l'est et de 1050–1000 m à l'ouest (Fig. 1). À la périphérie de ces cheminées, les réflecteurs sont rebroussés vers le haut, ou vers le bas, ou encore très peu affectés (Fig. 3c). Sur les profils sismiques (Figs. 2 et 4), les cheminées peuvent être suivies en profondeur jusqu'au BSR. Il est difficile de préciser si elles affectent effectivement les réflecteurs plus profonds ou si elles se terminent au niveau des anomalies d'amplitude, dans les horizons sous le BSR. Une analyse détaillée des profils sismiques et de sondeur de sédiments montre des anomalies d'amplitude qui affectent les séries non déformées à la verticale des cheminées scellées, et ce, jusqu'au fond de l'eau. Ces anomalies suggèrent une activité récente ou actuelle d'expulsion de fluides à travers ces cheminées. Particulièrement bien visibles sur les profils de sondeur de sédiment, mais aussi sur la sismique, des failles normales de faible rejeu (quelques mètres) limitent de petites zones effondrées. Le bloc sismique 3D permet de montrer que ces structures, qui pourraient être assimilées à de petits grabens d'extension sur une ligne 2D, ont en réalité une forme circulaire et pourraient représenter un stade peu évolué de ces cheminées.
4.3 Glissements
Le niveau de détail des cartes bathymétrique (Fig. 1) et d'imagerie apporte des informations essentielles : il permet de définir la limite des sédiments perturbés par les glissements récents ou encore met en évidence, à l'intérieur de la zone glissée, plusieurs plans de glissement situés à des niveaux différents, et dont les limites quasi linéaires semblent orientées dans des directions préférentielles (Fig. 1). Il est possible aussi d'identifier, en bordure de glissement, des zones perturbées, mais dont le déplacement horizontal est très faible. Sur les profils sismiques 2D, ces zones correspondent à des séries litées en surface, au-dessus d'un faciès acoustiquement sourd, avec quelques réflexions discontinues à la base.
Pour l'identification des instabilités sédimentaires anciennes, dans les séries de la zone investiguée, c'est l'analyse du bloc de données sismiques 3D qui se révèle particulièrement intéressante. La Fig. 4 représente la carte en 3D d'un horizon situé à environ 150 m sous le fond. Sur un profil 2D (Fig. 2), cet horizon est continu au nord-est, puis perturbé à partir du CDP 5500, et enfin très discontinu après le CDP 6500. Sur la carte, on montre immédiatement que cet horizon est affecté par un glissement dont l'ampleur du déplacement horizontal s'accroît à partir du CDP 6500 (Figs. 2 et 4). La partie supérieure du glissement (CDP 5500–6500) présente un déplacement faible ; la partie inférieure a plus largement mobilisé les sédiments sur une distance plus importante. Ces deux phases de glissement semblent toutefois scellées par les séries sédimentaires immédiatement postérieures, ce qui implique un déclenchement précoce de ces instabilités. Un second niveau glissé, avec un comportement similaire au précédent, peut être identifié à une profondeur de quelques dizaines de mètres sous le premier glissement. Les séries affectées par ces deux glissements appartiennent à la formation Naust. Plus précisément, le glissement supérieur affecte les séries basales de la formation Naust A2, datées de −250 000 ans BP [7].
5 Discussion
L'analyse des données de la campagne Hydratech, en particulier celle des cartes d'amplitude issues du bloc 3D, confirme le contrôle de la position et de l'amplitude du BSR par la nature des séries sédimentaires qui recoupent la limite inférieure de la BGHS. Comme le proposent Bouriak et al. [2] et Bünz et al. [5], le facteur principal de contrôle du BSR est probablement la capacité de ces séries à permettre la circulation et le piégeage du gaz libre sous la BGHS. Si le BSR est principalement marqué par le toit de séries chargées en gaz libre, sa présence n'est donc pas synonyme de fortes teneurs en hydrates au-dessus du BSR. L'absence de réflecteur susceptible de marquer le toit des hydrates ainsi que la faible teneur en hydrates (inférieure à quelques pour-cent de la porosité) estimée à partir de données de sismique réfraction sur la zone par Bünz et al. [5] corroborent ce raisonnement.
L'analyse conjointe des données de bathymétrie, de sondeur de sédiment et de sismique HR permet d'identifier de nombreuses cheminées de sortie de fluides et de montrer que celles-ci sont actives ou ont été actives très récemment car, même si elles ne sont pas marquées par un relief sur le fond de l'eau, l'amplitude du fond de la mer sur les profils chirp et sismique est affectée à la verticale de ces cheminées. À partir des données acquises, on montre que les cheminées s'enracinent sous le BSR, mais il n'est pas possible de préciser si c'est au niveau des séries chargées en gaz immédiatement sous le BSR, ou plus profondément. Leur répartition limitée au haut de pente pourrait être contrôlée par l'existence de surpressions liées à la présence, mise en évidence par les données sismiques, de gaz libre dans les séries sous le BSR.
Les deux glissements décrits dans le bloc sismique 3D se sont déclenchés peu après le dépôt des séries glissées, comme le prouve le scellement de ces séries par les sédiments postérieurs. Ils sont caractérisés par une faible épaisseur (de l'ordre de 15 m en moyenne pour le glissement supérieur), et par une extension horizontale régionale. À partir des données disponibles, il n'existe pas de lien évident entre la présence de gaz ou d'hydrates de gaz et ces deux glissements superficiels observés sur le bloc 3D : les glissements mis en évidence sont éloignés de la BGHS et l'étaient d'autant plus lors de leurs mises en place. Or, c'est principalement au niveau de la BGHS que la déstabilisation des hydrates de gaz contenus dans les sédiments est susceptible de se produire et d'influencer la stabilité de ces derniers.
Les plans de glissement identifiés ont joué ultérieurement un rôle lors du glissement principal de Storegga [7]. L'instabilité de ces séries aurait donc perduré depuis leur dépôt daté de −250 000 ans [7] jusqu'à −8200 ans, date de déclenchement du glissement de Storegga [8]. Sur son flanc nord, les zones où ces séries sont encore intactes pourraient alors continuer de présenter un risque d'instabilité important.
6 Conclusion
La combinaison des différentes données géophysiques acquises au cours de la campagne Hydratech apporte des informations sur les structures détaillées et la répartition des séries sédimentaires, du BSR, des sorties de fluides et des glissements sur la zone étudiée. Une première analyse des données permet (1) de confirmer le contrôle prépondérant de la position et de l'amplitude du BSR par les séries sédimentaires, (2) de montrer l'activité récente ou actuelle des sorties de fluides dans la zone d'étude, (3) de minimiser l'influence des hydrates de gaz dans le déclenchement des deux glissements superficiels identifiés sur le secteur étudié du flanc nord des glissements de Storegga, (4) de mettre en évidence le caractère durablement instable de certaines séries de la formation Naust, et donc l'instabilité potentielle des zones où ces séries sont encore intactes sur le flanc nord des glissements de Storegga.
En revanche, la question de l'impact de la possible déstabilisation des hydrates de gaz dans le déclenchement des glissements principaux de Storegga n'est pas résolue. La poursuite de l'exploitation des données de la campagne Hydratech est toutefois prometteuse.
Acknowledgments
La campagne Hydratech a été financée par l'Ifremer et le projet européen Hydratech. Nous remercions le commandant Thierry Alix, les officiers et l'équipage du N/O Le Suroît pour leur collaboration.