Sediments of Lake Vens (SW European Alps,France) record large-magnitude earthquake events

We studied sediment cores from Lake Vens (2,327 m asl), in the Tinée Valley of the SW Alps, to test the paleoseismic archive potential of the lake sediments in this particularly earthquake-sensitive area. The historical earthquake catalogue shows that moderate to strong earthquakes, with intensities of IX–X, have impacted the Southern Alps during the last millennium. Sedimentological (X-ray images, grain size distribution) and geochemical (major elements and organic matter) analyses show that Lake Vens sediments consist of a terrigenous, silty material (minerals and organic matter) sourced from the watershed and diatom frustules. A combination of X-ray images, grain-size distribution, major elements and magnetic properties shows the presence of six homogenite-type deposits interbedded in the sedimentary background. These sedimentological features are ascribed to sediment reworking and grain sorting caused by earthquake-generated seiches. The presence of microfaults that cross-cut the sediment supports the hypothesis of seismic deposits in this system. A preliminary sediment chronology is provided by ²¹⁰Pb measurement and AMS ¹⁴C ages. According to the chronology, the most recent homogenite events are attributable to damaging historic earthquakes in AD 1887 (Ligure) and 1564 (Roquebillière). Hence, the Lake Vens sediment recorded large-magnitude earthquakes in the region and permits a preliminary estimate of recurrence time for such events of ~400 years.     

New palaeontological data from the organic-rich layers of the Bathonian of the Grands Causses (France)

The palaeontological content (charophytes and vertebrates) of organic-rich layers from the Middle Jurassic (Bathonian) of the Grands Causses at Mostuéjouls (Aveyron, France) has been examined. It is rich in gyrogonites of Porochara douzensis. The vertebrates of the lower layer include a single hybodontiform tooth, “semionotiform” teeth and scales, and one pycnodontiform tooth whereas the upper layer has yielded one hybodontiform tooth, a variety of actinopterygian remains (mostly Caturus sp. and indeterminate pycnodontiforms) and a few possible reptile remains. The variation in fauna between the two layers is ascribed to a different degree of marine influence.     

Groundwater-driven temperature changes at thermal springs in response to recent glaciation: Bormio hydrothermal system,Central Italian Alps

Thermal springs are widespread in the European Alps, with hundreds of geothermal sites known and exploited. The thermal circulation and fluid outflows were examined in the area around Bormio (Central Italian Alps), where ten geothermal springs discharge from dolomite bodies located close to the regional Zebrù thrust. Water is heated in deep circulation systems and upwells vigorously at a temperature of about 40 °C. Heat and fluid transport is explored by steady and transient three-dimensional finite-element simulations taking into account the effect of the last glaciation, which in the study area was recognized to end around 11,000–12,000 years ago. The full regional model (ca. 700 km²) is discretized with a highly refined triangular finite-element planar grid. Numerical simulations suggest a reactivation of the system following the end of the Last Glacial Maximum. Results correctly simulate the observed discharge rate of ca. 2,400 L/min and the spring temperatures after ca. 13,000 years from deglaciation, and show a complete cooling of the aquifer within a period of approximately 50,000 years. Groundwater flow and temperature patterns suggest that thermal water flows through a deep system crossing both sedimentary and metamorphic lithotypes along a fracture network associated with the thrust system. This example gives insights into the influences of deep alpine structures and glaciations on groundwater circulation that control the development of many hydrothermal systems not necessarily associated with convective heat flow.     

Anhydrite–gypsum transition in the argillites of flooded salt workings in eastern France

This study aims at understanding the physico-chemical interactions between the saturated brine and the rocks enclosing the underground salt workings in Lorraine (eastern France). These anhydrite-rich and argillaceous rocks were characterized in terms of mineralogy, micro-texture and connected porosity. Then, the two main lithofacies, massive anhydrite and anhydrite-rich argillite, were immersed in brine during more than 1 year. During this batch experiment, the argillites were affected by macroscopic splitting, contrarily to the massive anhydrite. Micro-texture and brine chemical analyses clearly show the swelling due to the hydration of anhydrite into gypsum inside the argillites, whereas hydration occurs superficially on the massive anhydrite, due to its very low permeability. Anhydrite–gypsum transformation is promoted by the presence of dissolved strontium and potassium in saturated brine. The low activity of water in saturated brine does not allow the clay fraction to swell significantly during the experiment. Thus, the expansion resulting from the hydration of anhydrite into gypsum might be responsible of the splitting of argillite in a saturated brine environment. The superficial anhydrite hydration on massive anhydrite can be explained by the low amount of connected porosity (less than 1%).     

Effects of different crystal faces on the surface charge of colloidal goethite (α-FeOOH) particles: an experimental and modeling study

The surface charge of colloidal particles is usually determined by potentiometric titration. These acid-base titrations make it possible to measure the pH of point-of-zero charge (pzc) for oxide minerals. This macroscopic property is the most important parameter used in surface complexation modeling to reproduce experimental data. The pzc values of goethite reported in the literature vary between 7.0 and 9.5. Carbonate adsorption and/or surface morphology are thought to account for this wide range.We demonstrate a procedure for the removal of the carbonate ions that initially adsorb on goethite and strongly affect the titration curves and pzc determination. We also investigated the crystal-face-specific reactivity of two morphologically different goethites. The z-profiles obtained from atomic force microscopy (AFM) images showed that the goethite with the smallest specific surface area (S = 49 m²/g, denoted G49) exhibits 70% of the (001) face, whereas this value is only 30% for the goethite with largest specific surface area (S = 95 m²/g, denoted G95). This morphologic difference results in slightly different pzc values: 9.0 for G49 goethite and 9.1 for G95 geothite. These experimental pzc values have been correlated with multisite complexation calculations using both the full-site and the 1-pK approaches. We used the full-site approach to consider all of the configurations of hydrogen bond interactions with surface site. The resulting mean charges gave estimated pzc values of 8.9 and 9.2 for the (001) and (101) faces, respectively. Considering these theoretical pzc values for individual faces and the face distributions obtained from AFM analysis, the calculated pzc values are in full agreement with the experimental pzc values. However, this morphologic difference is more expressed in surface charge values than in the pzc values. Indeed, the surface charge of G49 goethite is much higher than that of G95 goethite, and the 1-pK calculations make it possible to fit the titration data satisfactorily.     

Numerical resolution of well-balanced shallow water equations with complex source terms

This paper presents a well-balanced numerical scheme for simulating frictional shallow flows over complex domains involving wetting and drying. The proposed scheme solves, in a finite volume Godunov-type framework, a set of pre-balanced shallow water equations derived by considering pressure balancing. Non-negative reconstruction of Riemann states and compatible discretization of slope source term produce stable and well-balanced solutions to shallow flow hydrodynamics over complex topography. The friction source term is discretized using a splitting implicit scheme. Limiting value of the friction force is derived to ensure stability. This new numerical scheme is validated against four theoretical benchmark tests and then applied to reproduce a laboratory dam break over a domain with irregular bed profile.