Multi-scale modeling for dynamics of structure-soil-structure interactions

This paper presents two cases of multi-scale modeling in the context of dynamic structure soil-structure interaction. The first case concerns the behaviour of reinforced soils. It is shown that such system may involve both shear and bending effect at the leading order, which corresponds to a second gradient material. The second case addresses the seismic response of soils in presence of a densely urbanized city. It appears that the effect of resonance of the whole buildings in interaction actually modify the seismic response. In both cases the theoretical approach is completed by a validation through analogous samples tested on a shaking table.     

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.     

Contrasting responses of oligochaetes (Annelida) and chironomids (Diptera) to the abatement of eutrophication in Lake Neuchâtel

: Lake Neuchâtel (Switzerland) has been recovering from human-induced eutrophication since 1982. Oligochaetes and chironomid larvae were used to monitor the recovery of sediment at a depth of 40 m. In oligochaete communities, the mean proportion of individuals belonging to species indicative of oligotrophic conditions (mostly Spirosperma velutinus and Stylodrilus heringianus) increased from 16 % in 1984 to 33 % in 1992, but decreased to 11 % in 1997. Based on oligochaetes, benthic conditions deteriorated in 1997 whereas, based on the increase of oligotrophic chironomids (mostly Micropsectra), they improved. This difference was attributed to the impact of large planktonic algae whose sedimentation had altered oxygen concentrations at the water-sediment interface. Oligochaetes were more affected than chironomids because they are less mobile and they depend more on the inner sediment for their food and reproduction than the latter. However, based on the combined response of oligochaetes and chironomids, benthic recovery was stalled in 1997 but was not reversed.     

Burial rates of organic matter along the eastern Canadian margin and stable isotope constraints on its origin and diagenetic evolution

Organic carbon (OC) and nitrogen (N) contents and δ¹³C and δ¹⁵N values in total organic matter (OM) were measured in sub-surface sediments (0–30 cm sub-bottom) from 21 cores raised from the Laurentian Channel of the Gulf of St. Lawrence and the Labrador Sea, to document OM fluxes and storage along the eastern Canadian margin. Storage rates as high as 2.5 g m⁻² yr⁻¹ for OC and 0.2 g m⁻² yr⁻¹ for N are observed in the Laurentian Channel, suggesting that the shelf plays a significant role in terms of OM storage (from 1 to 2% of the primary production). Based on the isotopic composition of the essentially marine OM of the Labrador Sea (δ¹³C/V-PDB=−21.9±0.4‰; δ¹⁵N/AIR=7.6±0.6‰; n=12), there is no isotopic evidence for a significant relative input of terrestrial OM along the Laurentian Channel (δ¹³C/V-PDB=−21.9±0.4‰; δ¹⁵N/AIR=8.0±0.9‰; n=10), either due to high relative fluxes of marine OM and/or to the trapping of continental OM in the estuary and upstream. High storage rates of OM are also observed on the continental rise of the Labrador Sea (as high as 1.1 g C m⁻² yr⁻¹ and 0.09 g N m⁻² yr⁻¹). They contrast with one order of magnitude lower rates on the slope, due to low sedimentation rates (SR) and sediment winnowing by the Western Boundary Undercurrent (WBUC). Reduced early diagenetic alteration of OM is observed, particularly in the Laurentian Channel. It results in discrete (i) losses of OC and N, (ii) shifts in C/N ratios, suggesting preferential removal of N-bearing OM also highlighted by losses in total hydrolysable amino acids (HAA). In the Labrador Sea slope records, due to low SR, OM concentration changes linked to long term temporal variations may superimpose on these diagenetic trends, and some influence of the WBUC is noticeable.