Habitat heterogeneity influences cold‐seep macrofaunal communities within and among seeps along the Norwegian margin. Part 1: macrofaunal community structure

Håkon Mosby mud volcano (HMMV) is one of the most active and most studied seep sites in European waters. Many authors have described its thermal activity, dynamic of mud flows, and geochemical and microbial processes. It is characterised by a concentric zonation of successive biogenic habitats related to an activity and geochemical gradient from its centre to its periphery. Around the central area covered by mud flows, white and grey microbial mats occur among areas of bare sediment, whereas siboglinid tubeworm fields of Sclerolinum contortum and/or Oligobrachia haakonmosbiensis colonise the peripheral areas. The meiofaunal community is known to be structured among habitats, but the macrofauna has rarely been investigated and has never been sampled in situ. As part of the European project HERMES, using the ROVs Victor 6000 and Quest 4000, we sampled quantitatively the different habitats of the volcano for macrofauna sensus lato, retained on a 250‐ or 500‐μm sieve. We also sampled a newly discovered pockmark on Storegga slide (cne 5.6) and two pockmarks (G11, G12) in the Nyegga area. Macrofauna was identified and counted from phylum to family level. Our results on HMMV showed a gradient of increasing density and diversity from the volcano centre (1–3 taxa; 260 ind·m^?2) to the peripheral siboglinid fields (8–14 taxa, 93,000 ind·m^?2), with an intermediate situation for microbial mats. For macrofauna ≥500 μm, non‐siboglinid polychaetes dominated the communities of the central mud volcano area, white mats and S. contortum fields (83, 89 and 37% of the total, respectively), whereas gastropods dominated grey mats and O. haakonmosbiensis fields (89 and 44% of the total, respectively). Polychaete families followed the same pattern of diversity according to habitats within HMMV. Of 23 polychaete families identified, only one occurred in the centre, and three in the microbial mats. Capitellidae and Dorvilleidae (typical of organically and sulphide‐enriched areas) occurred at remarkably high densities in white microbial mats and in O. haakonmosbiensis fields. The S. contortum fields were the most diverse habitat with 12 polychaete families. The 250‐μm fraction showed similar taxa dominating the habitats, but taking meiofauna into account, nematodes became the major taxon in white mats and in S. contortum fields, where they were particularly large in size, whereas copepods dominated in other habitats. Meiofauna and macrofauna did not show the same patterns of density according to habitats. Using principal components analysis the habitats at HMMV were clearly distinct, and clustered according to dominant species of siboglinids and type of microbial mats. Pockmarks at Nyegga showed a similar concentric pattern of habitats around fluid sources as on the volcano, which seemed similarly to influence macrofauna composition, but at a much smaller scale. Total taxa and polychaete diversity are high in the S. contortum fields in these pockmarks as well. Regional‐scale comparisons including HMMV and Storegga suggested a higher influence of habitat‐type than seep‐site on the community structure.     

Alternate bars in a sandy gravel bed river: generation,migration and interactions with superimposed dunes

A field study was carried out to investigate the development of alternate bars in a secondary channel of the Loire River (France) as a function of discharge variations. We combined frequent bathymetric surveys, scour chains and stratigraphical analysis of deposits with measurements and modelling of flow dynamics. The channel exhibited migrating bars, non‐migrating bars and superimposed dunes. Possible mechanisms of bar initiation were found to be chutes associated with changes of bank direction and instability resulting from interactions between existing bars during the fall in water level after floods. We propose that the reworking of bar sediments during low flows (high width‐to‐depth ratio β), reinforced by high values of the Shields mobility parameter, can explain the formation or re‐generation of new alternate migrating bars during a subsequent flood. The migration pattern of the bars was found to be cyclic and to depend mainly on (i) channel layout and (ii) the dynamics of superimposed dunes with heights and lengths depending on location and discharge value. For instance, the hysteresis affecting the steepness of dunes influences the flow resistance of the dunes as well as the celerity of migrating bars during flood events. We compare the findings from the field with results from theoretical studies on alternate bars. This gives insight in the phenomena occurring in the complex setting of real rivers, but it also sheds light on the extent to which bar theories based on idealized cases can predict those phenomena. Copyright © 2014 John Wiley & Sons, Ltd.     

Upscaling of CO2 vertical migration through a periodic layered porous medium: The capillary-free and capillary-dominant cases

We present an upscaled model for the vertical migration of a CO₂ plume through a vertical column filled with a periodic layered porous medium. This model may describe the vertical migration of a CO₂ plume in a perfectly layered horizontal aquifer. Capillarity and buoyancy are taken into account and semi-explicit upscaled flux functions are proposed in the two following cases: (i) capillarity is the main driving force and (ii) buoyancy is the only driving force. In both cases, we show that the upscaled buoyant flux is a bell-shaped function of the saturation, as in the case of a homogeneous porous medium. In the capillary-dominant case, we show that the upscaled buoyant flux is the harmonic mean of the buoyant fluxes in each layer. The upscaled saturation is governed by the continuity of the capillary pressure at the interface between layers. In the capillary-free case, the upscaled buoyant flux and upscaled saturation are determined by the flux continuity condition at the interface. As the flux is not continuous over the entire range of saturation, the upscaled saturation is only defined where continuity is verified, i.e. in two saturation domains. As a consequence, the upscaled buoyant flux is described by a piecewise continuous function. Two analytical approximations of this flux are proposed and this capillary-free upscaled model is validated for two cases of heterogeneity. Upscaled and cell averaged saturations are in good agreement. Furthermore, the proposed analytical upscaled fluxes provide satisfactory approximations as long as the saturation set at the inlet of the column is in a range where analytical and numerical upscaled fluxes are close.     

TIMS measurements of Ra and Ra in the Gulf of Lion, an attempt to quantify submarine groundwater discharge

Submarine groundwater discharge (SGD) is now recognized as an important pathway for water and chemical species fluxes to the coastal ocean. In order to determinate SGD to the Gulf of Lion (France), we measured the activities of ²²⁶Ra and ²²⁸Ra by thermal ionization mass spectrometry (TIMS) in coastal waters and in the deep aquifer waters of the Rhone deltaic plain after pre-concentration of radium by MnO₂. Compared to conventional counting techniques, TIMS requires lower quantities of water for the analyses, and leads to higher analytical precision. Radium isotopes were thus measured on 0.25–2 L water samples containing as little as 20 fg of ²²⁶Ra and 0.2–0.4 fg of ²²⁸Ra with precision equal to 2%. We demonstrate that coastal surface waters samples are enriched in ²²⁶Ra and ²²⁸Ra compared to the samples further offshore. The high precision radium measurements display a small but significant ²²⁶Ra and ²²⁸Ra enrichment within a strip of circa 30 km from the coast. Radium activities decrease beyond this region, entrained in the northern current along the shelf break or controlled by eddy diffusion. The radium excess in the first 30 km cannot be accounted for by the river nor by the early diagenesis. The primary source of the radium enrichment must therefore be ascribed to the discharge of submarine groundwater. Using a mass-balance model, we estimated the advective fluxes of ²²⁶Ra and ²²⁸Ra through SGD to be 5.2 × 10¹⁰ and 21 × 10¹⁰ dpm/d respectively. The ²²⁶Ra activities measured in the groundwater from the Rhone deltaic plain aquifer are comparable to those from other coastal groundwater studies throughout the world. By contrast, ²²⁸Ra activities are higher by up to one order of magnitude. Taking those groundwater radium activities as typical of the submarine groundwater end-member, a minimum volume of 0.24–4.5 × 10¹⁰ l/d is required to support the excess radium isotopes on the inner shelf. This has to be compared with the average rivers water runoff of 15.4 × 10¹⁰ l/d during the study period (1.6 to 29% of the river flow).     

Atlantic multidecadal oceanic variability and its influence on the atmosphere in a climate model

The mechanisms controlling the decadal to multidecadal variability of the Atlantic Meridional Overturning Circulation (MOC) and its influence on the atmosphere are investigated using a control simulation with the IPSL-CM4 climate model. The multidecadal fluctuations of the MOC are mostly driven by deep convection in the subpolar gyre, which occurs south of Iceland in the model. The latter is primarily influenced by the anomalous advection of salinity due to changes in the East Atlantic Pattern (EAP), which is the second mode of atmospheric variability in the North Atlantic region. The North Atlantic Oscillation is the dominant mode, but it plays a secondary role in the MOC fluctuations. During summer, the MOC variability is shown to have a significant impact on the atmosphere in the North Atlantic–European sector. The MOC influence is due to an interhemispheric sea surface temperature (SST) anomaly with opposite signs in the two hemispheres but largest amplitude in the northern one. The SST pattern driven by the MOC mostly resembles the model Atlantic Multidecadal Oscillation (AMO) and bears some similarity with the observed one. It is shown that the AMO reflects both the MOC influence and the local atmospheric forcing. Hence, the MOC influence on climate is best detected using lagged relations between climatic fields. The atmospheric response resembles the EAP, in a phase that might induce a weak positive feedback on the MOC.     

Forced regressive deposits of a deglaciation sequence: Example from the Late Quaternary succession in the Lake Saint‐Jean basin (Québec,Canada)

Differentiating between forced regressive deposits from deglacial periods in high latitude domains and forced regressive deposits from the onset of glacial periods in low latitude domains is fundamental for the accurate interpretation of glacial cycles within the geological record and then for the reconstruction of palaeogeography and palaeo‐climate. A forced regressive deglacial sequence is documented from the Lake Saint‐Jean basin (Québec, Canada). In this area, the Late Pleistocene to Holocene sediments have recorded the Laurentide ice sheet retreat accompanied by the invasion of marine waters (Laflamme Gulf) from ca 12·9 cal kyr bp . Subsequently, fluvio‐deltaic and coastal prograding wedges were deposited; they followed the base‐level fall due to glacio‐isostatic rebound. This succession, representing a transition from glacial to post‐glacial periods within a previously glaciated area, was investigated through recent mapping, preserved landforms, facies analysis, and new optical stimulated luminescence and radiocarbon dates. Three basin‐scale geological sections share a common lower part made of isolated ice‐contact fan deposits overlying bedrock. Throughout the entire basin, ice‐contact fans are capped by glacimarine muds. Above, fluvial and coastal prograding systems were deposited and evolved through four steps: (i) deltaic systems progressively increased in width; (ii) coastal influence on sedimentation increased; (iii) hydrographic drainage systems became more organised; and (iv) deltas graded from steep (Gilbert delta) to low‐angle foresets (mouth‐bar delta). Deposited during the base‐level fall from glacio‐isostatic rebound, the complete succession has been designated as a single falling stage system tract referred to as a deglacial falling stage system tract. It is representative of a deglaciation sequence in areas previously covered by ice during glacial periods (i.e. medium to high latitude domains). Diagnostic criteria are provided to identify such a deglacial falling stage system tract in the geological record, which may aid identification of previously unknown glacial cycles.     

Mechanisms of the atmospheric response to North Atlantic multidecadal variability: a model study

The atmospheric circulation response to decadal fluctuations of the Atlantic meridional overturning circulation (MOC) in the IPSL climate model is investigated using the associated sea surface temperature signature. A SST anomaly is prescribed in sensitivity experiments with the atmospheric component of the IPSL model coupled to a slab ocean. The prescribed SST anomaly in the North Atlantic is the surface signature of the MOC influence on the atmosphere detected in the coupled simulation. It follows a maximum of the MOC by a few years and resembles the model Atlantic multidecadal oscillation. It is mainly characterized by a warming of the North Atlantic south of Iceland, and a cooling of the Nordic Seas. There are substantial seasonal variations in the geopotential height response to the prescribed SST anomaly, with an East Atlantic Pattern-like response in summer and a North Atlantic oscillation-like signal in winter. In summer, the response of the atmosphere is global in scale, resembling the climatic impact detected in the coupled simulation, albeit with a weaker amplitude. The zonally asymmetric or eddy part of the response is characterized by a trough over warm SST associated with changes in the stationary waves. A diagnostic analysis with daily data emphasizes the role of transient-eddy forcing in shaping and maintaining the equilibrium response. We show that in response to an intensified MOC, the North Atlantic storm tracks are enhanced and shifted northward during summer, consistent with a strengthening of the westerlies. However the anomalous response is weak, which suggests a statistically significant but rather modest influence of the extratropical SST on the atmosphere. The winter response to the MOC-induced North Atlantic warming is an intensification of the subtropical jet and a southward shift of the Atlantic storm track activity, resulting in an equatorward shift of the polar jet. Although the SST anomaly is only prescribed in the Atlantic ocean, significant impacts are found globally, indicating that the Atlantic ocean can drive a large scale atmospheric variability at decadal timescales. The atmospheric response is highly non-linear in both seasons and is consistent with the strong interaction between transient eddies and the mean flow. This study emphasizes that decadal fluctuations of the MOC can affect the storm tracks in both seasons and lead to weak but significant dynamical changes in the atmosphere.     

Stochastically-driven multidecadal variability of the Atlantic meridional overturning circulation in CCSM3

The Atlantic meridional overturning circulation (AMOC) in the last 250?years of the 700-year-long present-day control integration of the Community Climate System Model version 3 with T85 atmospheric resolution exhibits a red noise-like irregular multi-decadal variability with a persistence longer than 10?years, which markedly contrasts with the preceding ~300 years of very regular and stronger AMOC variability with ~20?year periodicity. The red noise-like multi-decadal AMOC variability is primarily forced by the surface fluxes associated with stochastic changes in the North Atlantic Oscillation (NAO) that intensify and shift northward the deep convection in the Labrador Sea. However, the persistence of the AMOC and the associated oceanic anomalies that are directly forced by the NAO forcing does not exceed about 5?years. The additional persistence originates from anomalous horizontal advection and vertical mixing, which generate density anomalies on the continental shelf along the eastern boundary of the subpolar gyre. These anomalies are subsequently advected by the mean boundary current into the northern part of the Labrador Sea convection region, reinforcing the density changes directly forced by the NAO. As no evidence was found of a clear two-way coupling with the atmosphere, the multi-decadal AMOC variability in the last 250?years of the integration is an ocean-only response to stochastic NAO forcing with a delayed positive feedback caused by the changes in the horizontal ocean circulation.     

Stationarity Scores on Training Images for Multipoint Geostatistics

This research introduces a novel method to assess the validity of training images used as an input for Multipoint Geostatistics, alternatively called Multiple Point Simulation (MPS). MPS are a family of spatial statistical interpolation algorithms that are used to generate conditional simulations of property fields such as geological facies. They are able to honor absolute “hard” constraints (e.g., borehole data) as well as “soft” constraints (e.g., probability fields derived from seismic data, and rotation and scale). These algorithms require 2D or 3D training images or analogs whose textures represent a spatial arrangement of geological properties that is presumed to be similar to that of a target volume to be modeled. To use the current generation of MPS algorithms, statistically valid training image are required as input. In this context, “statistical validity” includes a requirement of stationarity, so that one can derive from the training image an average template pattern. This research focuses on a practical method to assess stationarity requirements for MPS algorithms, i.e., that statistical density or probability distribution of the quantity shown on the image does not change spatially, and that the image shows repetitive shapes whose orientation and scale are spatially constant. This method employs image-processing techniques based on measures of stationarity of the category distribution, the directional (or orientation) property field and the scale property field of those images. It was successfully tested on a set of two-dimensional images representing geological features and its predictions were compared to actual realizations of MPS algorithms. An extension of the algorithms to 3D images is also proposed. As MPS algorithms are being used increasingly in hydrocarbon reservoir modeling, the methods described should facilitate screening and selection of the input training images.     

Deltaic plain development and environmental changes in the Petite Camargue, Rhone Delta, France, in the past 2000 years

The deltaic plain of the Petite Camargue which constitutes the western part of the Rhone Delta, began its main progradation around 2000 yr ago. Several delta lobes follow each other and have participated in the deltaic evolution. The deltaic lobes have distinct morphologies which reflect the dynamic fluvial and marine processes under the influence of climatic and human controls. Two delta lobe systems were built by the Daladel and Peccaïs channels, after which a deflected wave-influenced delta lobe was formed by the La Ville and Saint-Roman channels. The latest channel, the Rhone Vif channel, is skewed because this channel was completely canalized and engineered up to its mouth in the beginning of the 16th century. Since the avulsion of this channel about 1550 A.D., the coastline of the Petite Camargue has been especially affected by the influence of waves and currents. The spits replaced the beach ridges which juxtaposed themselves and have migrated westward since the 16th century. The formation of the western part of the delta in the last 2000 yr is affected by not only the fluvial sedimentary fluxes and the coastal dynamics to the mouth but also climatic change and human influence.