Role of the bedrock topography in the Quaternary filling of a giant estuarine basin: the Lower St. Lawrence Estuary,Eastern Canada

The geometry of estuarine and/or incised‐valley basins and their protected character compared with open sea basins are favourable for the preservation of sedimentary successions. The Lower St. Lawrence Estuary Basin (LSLEB, eastern Canada) is characterized by a thick (>400 m in certain areas) Quaternary succession. High‐ and very high‐resolution seismic reflection data, multibeam bathymetry coverage completed by core and chronostratigraphic data as well as a 3‐D seismic stratigraphic model are used to document the geometrical relationships between the bedrock and the Quaternary units of the LSLEB. The bedrock geometry of LSLEB is characterized by two large troughs that are interpreted as resulting mainly from repeated (?) periods of glacial overdeepening of a pre‐Quaternary drainage system. However, other mechanisms with complex feedback effects such as differential glacio‐isostatic uplift, erosion, sedimentary supply, and subsidence may have contributed to the formation of bedrock troughs. The two large bedrock troughs are mostly filled by ～200 m thick Wisconsinan (Marine Isotopic Stages 2–4) and possibly older sediments. Overlying units recorded the retreat of the Laurentian Ice Sheet during the Late Wisconsinan (Marine Isotopic Stage 2) and estuarine conditions during the Holocene. The strong correlation existing between the bedrock topography and the thickness of the Quaternary succession is indicative of the effectiveness of the LSLEB as a sediment trap.     

Conditioning of Multiple-Point Statistics Facies Simulations to Tomographic Images

Geophysical tomography captures the spatial distribution of the underlying geophysical property at a relatively high resolution, but the tomographic images tend to be blurred representations of reality and generally fail to reproduce sharp interfaces. Such models may cause significant bias when taken as a basis for predictive flow and transport modeling and are unsuitable for uncertainty assessment. We present a methodology in which tomograms are used to condition multiple-point statistics (MPS) simulations. A large set of geologically reasonable facies realizations and their corresponding synthetically calculated cross-hole radar tomograms are used as a training image. The training image is scanned with a direct sampling algorithm for patterns in the conditioning tomogram, while accounting for the spatially varying resolution of the tomograms. In a post-processing step, only those conditional simulations that predicted the radar traveltimes within the expected data error levels are accepted. The methodology is demonstrated on a two-facies example featuring channels and an aquifer analog of alluvial sedimentary structures with five facies. For both cases, MPS simulations exhibit the sharp interfaces and the geological patterns found in the training image. Compared to unconditioned MPS simulations, the uncertainty in transport predictions is markedly decreased for simulations conditioned to tomograms. As an improvement to other approaches relying on classical smoothness-constrained geophysical tomography, the proposed method allows for: (1) reproduction of sharp interfaces, (2) incorporation of realistic geological constraints and (3) generation of multiple realizations that enables uncertainty assessment.     

Design of sediment detention basins: Scaled model experiments and application

Sediment detention basins are implemented on mountain rivers to trap solid material that may aggravate the flooding of downstream settlements.However,retention structures built in the past may unnecessarily retain sediment during non-hazardous flood events,resulting in high maintenance costs and sediment deficit downstream.In addition,the so-called spontaneous self-flushing of previously retained sediment during floods has occasionally been observed.Recent research suggests to design sediment detention basins for controlling sediment passage with a guiding channel across the deposition area upstream of a hybrid barrier.Such barriers consist of a check dam with a slot orifice and an upstream bar screen with a bottom clearance in order to benefit from a combined mechanical-hydraulic retention control.The present paper enhances this pioneering research with the help of new experimental data,including a wide range of sediment mixtures and large wood,as well as variable barrier heights.Improved design criteria are provided regarding the bar screen and the basin storage capacity.The functionality of the enhanced concept for sediment detention is illustrated by a case study on a physical model:the protection of the Villard-Bonnot village(France)against torrential hazards.     

First-principles study of the structural and isotopic properties of Al- and OH-bearing hematite

A series of hematite structures containing various amounts of aluminum and hydroxyl groups was modeled using first-principles methods based on the density functional theory. Evolution of the lattice parameters was quantified as a function of Al and H concentrations. The a and c lattice parameters decrease with the aluminum content and increase with the water content. This allows explaining experimental data reported for synthetic hematites, in particular the observed deviation from the Vegard’s rule. Investigation of several hydroxyl configurations associated with cationic vacancies suggests that the speciation of water also significantly affects the hematite structure. The ⁵⁷Fe and ¹⁸O reduced partition functions (β-factor) were determined. Results show a linear dependence of the iron and oxygen β-factors on the aluminum content. An incorporation of 18 mole % Al₂O₃ in hematite would increase the iron β-factor of about 0.6‰ and the oxygen β-factor of about 5.5‰ at 0 °C. This effect is sufficiently large to be measurable and to affect the interpretation of natural isotopic compositions. On the other hand, the effect of water is found to be negligible for the hydroxyl configuration investigated.     

Stochastic simulations of fault networks in 3D structural modeling

3D structural modeling is a major instrument in geosciences, e.g. for the assessment of groundwater and energy resources or nuclear waste underground storage. Fault network modeling is a particularly crucial step during this task, for faults compartmentalize rock units and plays a key role in subsurface flow, whether faults are sealing barriers or drains. Whereas most structural uncertainty modeling techniques only allow for geometrical changes and keep the topology fixed, we propose a new method for creating realistic stochastic fault networks with different topologies. The idea is to combine an implicit representation of geological surfaces which provides new perspectives for handling topological changes with a stochastic binary tree to represent the spatial regions. Each node of the tree is a fault, separating the space in two fault blocks. Changes in this binary tree modify the fault relations and therefore the topology of the model.     

Preferred orientation distribution of shock‐induced planar microstructures in quartz and feldspar

Shocked quartz and feldspar grains commonly exhibit planar microstructures, such as planar fractures, planar deformation features, and possibly microtwins, which are considered to have formed by shock metamorphism. Their orientation and frequency are typically reported to be randomly distributed across a sample. The goal of this study is to investigate whether such microstructures are completely random within a given sample, or whether their orientation might also retain information on the direction of the local shock wave propagation. For this work, we selected samples of shatter cones, which were cut normal to the striated surface and the striation direction, from three impact structures (Keurusselkä, Finland, and Charlevoix and Manicouagan, Canada). These samples show different stages of pre‐impact tectonic deformation. Additionally, we investigated several shocked granite samples, selected at different depths along the drill core recovered during the joint IODP‐ICDP Chicxulub Expedition 364 (Mexico). In this case, thin sections were cut along two orthogonal directions, one parallel and one normal to the drill core axis. All the results refer to optical microscopy and universal‐stage analyses performed on petrographic thin sections. Our results show that such shock‐related microstructures do have a preferred orientation, but also that relating their orientation with the possible shock wave propagation is quite challenging and potentially impossible. This is largely due to the lack of dedicated experiments to provide a key to interpret the observed preferred orientation and to the lack of information on postimpact orientation modifications, especially in the case of the drill core samples.     

Numerical modelling of stream–aquifer interaction: Quantifying the impact of transient streambed permeability and aquifer heterogeneity

Stream–aquifer interaction plays a vital role in the water cycle, and a proper study of this interaction is needed for understanding groundwater recharge, contaminants migration, and for managing surface water and groundwater resources. A model‐based investigation of a field experiment in a riparian zone of the Schwarzbach river, a tributary of the Rhine River in Germany, was conducted to understand stream–aquifer interaction under alternative gaining and losing streamflow conditions. An equivalent streambed permeability, estimated by inverting aquifer responses to flood waves, shows that streambed permeability increased during infiltration of stream water to aquifer and decreased during exfiltration. Aquifer permeability realizations generated by multiple‐point geostatistics exhibit a high degree of heterogeneity and anisotropy. A coupled surface water groundwater flow model was developed incorporating the time‐varying streambed permeability and heterogeneous aquifer permeability realizations. The model was able to reproduce varying pressure heads at two observation wells near the stream over a period of 55 days. A Monte Carlo analysis was also carried out to simulate groundwater flow, its age distribution, and the release of a hypothetical wastewater plume into the aquifer from the stream. Results of this uncertainty analysis suggest (a) stream–aquifer exchange flux during the infiltration periods was constrained by aquifer permeability; (b) during exfiltration, this flux was constrained by the reduced streambed permeability; (c) the effect of temporally variable streambed permeability and aquifer heterogeneity were found important to improve the accurate capture of the uncertainty; and (d) probabilistic infiltration paths in the aquifer reveal that such pathways and the associated prediction of the extent of the contaminant plume are highly dependent on aquifer heterogeneity.     

Spatial assessment of fishing effort around European marine reserves: implications for successful fisheries management

We examined the spatial dynamic of artisanal fishing fleets around five European marine protected areas (MPAs) to derive general implications for the evaluation of MPAs as fisheries management tools. The coastal MPAs studied were located off France, Malta and Spain and presented a variety of spatial designs and processes of establishment. We developed a standardized methodology to define factors influencing effort allocation and to produce fishing effort maps by merging GIS with geostatistical modelling techniques. Results revealed that in most cases the factors “distance to the no-take”, “water depth”, and “distance to the port” had a significant influence on effort allocation by the fishing fleets. Overall, we found local concentration of fishing effort around the MPA borders. Thus, neglecting the pattern of fishing effort distribution in evaluating MPA benefits, such as spillover of biomass, could hamper sound interpretation of MPAs as fisheries management tools.