The Grand Brûlé exploration drilling: New data on the deep framework of the Piton de la Fournaise volcano. Part 2: Secondary minerals

The Grand Brûlé borehole intersects a thick pile of basaltic lavas, down to 1010 m, and a basic-ultrabasic intrusive complex, from 1010 to 3003.50 m.The lavas are, in general, unaltered except in two fractured zones, where hydrothermal fluids circulated at temperatures not exceeding 350 ° C. The main secondary minerals there are pyroxene, feldspar, epidote, actinolite and chlorite.The entire thickness of the intrusive body intersected contains secondary minerals representing three stages of cooling:

1. 1. A late magmatic episode (600–900°C) characterized by biotite, kaersutite, edenite and pargasite.

2. 2. A hydrothermal episode (T<- 350°C) characterized by epidote, albite, biotite, actinolite and chlorite.

3. 3. A phase of serpentinization (T≤ 350°C).

It is very likely that the two later events occurred simultaneously, with physical and chemical interference.     

Towards Improved Remedial Outcomes in Categorical Aquifers with an Iterative Ensemble Smoother

Categorical parameter distributions consisting of geologic facies with distinct properties, for example, high-permeability channels embedded in a low-permeability matrix, are common at contaminated sites. At these sites, low-permeability facies store solute mass, acting as secondary sources to higher-permeability facies, sustaining concentrations for decades while increasing risk and cleanup costs. Parameter estimation is difficult in such systems because the discontinuities in the parameter space hinder the inverse problem. This paper presents a novel approach based on Traveling Pilot Points (TRIPS) and an iterative ensemble smoother (IES) to solve the categorical inverse problem. Groundwater flow and solute transport in a hypothetical aquifer with a categorical parameter distribution are simulated using MODFLOW 6. Heads and concentrations are recorded at multiple monitoring locations. IES is used to generate posterior ensembles assuming a TRIPS prior and an approximate multi-Gaussian prior. The ensembles are used to predict solute concentrations and mass into the future. The evaluation also includes an assessment of how the number of measurements and the choice of the geological prior determine the characteristics of the posterior ensemble and the resulting predictions. The results indicate that IES was able to efficiently sample the posterior distribution and showed that even with an approximate geological prior, a high degree of parameterization and history matching could lead to parameter ensembles that can be useful for making certain types of predictions (heads, concentrations). However, the approximate geological prior was insufficient for predicting mass. The analysis demonstrates how decision-makers can quantify uncertainty and make informed decisions with an ensemble-based approach.     

Cretaceous length of day perturbation by mantle avalanche

These last 10 years, numerical models of mantle convection have emphasized the role of the 670 km endothermic phase change in generating avalanches that trigger catastrophic mass transfers between upper and lower mantle. On the other hand, scientists have emphasized the concomitance of large-scale worldwide geophysical and tectonic events, which could find their deep thermal roots in the huge mass transfers induced by the avalanches. In particular, the paleontological records show two periods of length of day (l.o.d.) shortening between 420 and 360, and 200 and 80 Myr BP. This last event is synchronous with a strong true polar wander and a global warming of the upper mantle. In order to study the potential effects of the avalanche on the main component of the Earth’s rotation, the Liouville equation has been solved and the l.o.d. evolution has been calculated from the perturbations of the inertia tensor. The results show that the inertia tensor of the Earth’s is mainly sensitive to the global transfers through the 670 km discontinuity. The l.o.d. perturbations will be synchronous with the global thermal effects of the avalanche. These theoretical results allow proposing a self-consistent physical mechanism to explain periods of the Earth’s rotation acceleration. Within this context, the l.o.d. shortening during the Cenozoic and Cretaceous brings one more clue to the possible participation of a mantle avalanche in generating the concomitant large scale events which have occurred during this very particular period of the Earth’s history.     

Eruptive history of the Colima volcanic complex (Mexico)

The evolution of the Colima volcanic complex can be divided into successive periods characterized by different dynamic and magmatic processes: emission of andesitic to dacitic lava flows, acid-ash and pumice-flow deposits, fallback nuées ardentes leading to pyroclastic flows with heterogeneous magma, plinian air-fall deposits, scoriae cones of alkaline and calc-alkaline nature. Four caldera-forming events, resulting either from major ignimbrite outbursts or Mount St. Helens-type eruptions, separate the main stages of development of the complex from the building of an ancient shield volcano (25 × 30 km wide) up to two summit cones, Nevado and Fuego.The oldest caldera, C1 (7–8 km wide), related to the pouring out of dacitic ash flows, marks the transition between two periods of activity in the primitive edifice called Nevado I: the first one, which is at least 0.6 m.y. old, was mainly andesitic and effusive, whereas the second one was characterized by extrusion of domes and related pyroclastic products. A small summit caldera, C2 (3–3.5 km wide), ended the evolution of Nevado I.Two modern volcanoes then began to grow. The building of the Nevado II started about 200,000 y. ago. It settled into the C2 caldera and partially overflowed it. The other volcano, here called Paleofuego, was progressively built on the southern side of the former Nevado I. Some of its flows are 50,000 y. old, but the age of its first outbursts is not known. However, it is younger than Nevado II. These two modern volcanoes had similar evolutions. Each of them was affected by a huge Mount St. Helens-type (or Bezymianny-type) event, 10,000 y. ago for the Paleofuego, and hardly older for the Nevado II. The landslides were responsible for two horseshoe-shaped avalanche calderas, C3 (Nevado) and C4 (Paleofuego), each 4–5 km wide, opening towards the east and the south. In both cases, the activity following these events was highly explosive and produced thick air-fall deposits around the summit craters.The Nevado III, formed by thick andesitic flows, is located close to the southwestern rim of the C3 caldera. It was a small and short-lived cone. Volcan de Fuego, located at the center of the C4 caldera, is nearly 1500 m high. Its activity is characterized by an alternation of long stages of growth by flows and short destructive episodes related to violent outbursts producing pyroclastic flows with heterogeneous magma and plinian air falls.The evolution of the primitive volcano followed a similar pattern leading to formation of C1 and then C2. The analogy between the evolutions of the two modern volcanoes (Nevado II–III; Paleofuego-Fuego) is described. Their vicinity and their contemporaneous growth pose the problem of the existence of a single reservoir, or two independent magmatic chambers, after the evolution of a common structure represented by the primitive volcano.     

Morphological analysis of active Mount Nemrut stratovolcano, eastern Turkey: evidences and possible impact areas of future eruption

Mount Nemrut, an active stratovolcano in eastern Turkey, is a great danger for its vicinity. The volcano possesses a summit caldera which cuts the volcano into two stages, i.e. pre- and post-caldera. Wisps of smoke and hot springs are to be found within the caldera. Although the last recorded volcanic activity is known to have been in 1441, we consider here that the last eruption of Nemrut occurred more recently, probably just before 1597. The present active tectonic regime, historical eruptions, occurrence of mantle-derived magmatic gases and the fumarole and hot spring activities on the caldera floor make Nemrut Volcano a real danger for its vicinity. According to the volcanological past of Nemrut, the styles of expected eruptions are well-focused on two types: (1) occurrence of water within the caldera leads to phreatomagmatic (highly energetic) eruptions, subsequently followed by lava extrusions, and (2) effusions–extrusions (non-explosive or weakly energetic eruptions) on the flanks from fissures. To predict the impact area of future eruptions, a series of morphological analyses based on field observations, Digital Elevation Model and satellite images were realized. Twenty-two valleys (main transport pathways) were classified according to their importance, and the physical parameters related to the valleys were determined. The slope values in each point of the flanks and the Heim parameters H/L were calculated. In the light of morphological analysis the possible impact areas around the volcano and danger zones were proposed. The possible transport pathways of the products of expected volcanic events are unified in three main directions: Bitlis, Guroymak, Tatvan and Ahlat cities, the about 135 000 inhabitants of which could be threatened by future eruptions of this poorly known and unsurveyed volcano.     

An operational method for the real-time monitoring of E. coli numbers in bathing waters

The aim of this study was to investigate the potential application of the beta-d-glucuronidase (GLUase) activity measurement for the routine detection and quantification of E. coli in marine bathing waters. GLUase activity was measured as the rate of hydrolysis of 4-methylumbelliferyl-beta-d-glucuronide. Culturable E. coli were quantified by the most probable number (MPN) microplate method. Both methods were applied to a large set of seawater samples. Significant correlation was found between the log of GLUase activity and the log of culturable E. coli. The mean coefficient of variation (CV) of the GLUase activity was less than 15% at concentrations around the current standards of International regulations whereas the CV of the microplate method was around 30%. When samples were stored at 4 degrees C and 20 degrees C, the mean CV of the GLUase activity remained below 15% up to 6 hours after sample collection whereas the range of variation of the microplate method varied between 10 and 50%. We concluded that the GLUase activity is an operational, reproducible, simple, very rapid and low cost method for the real-time enumeration of E. coli in bathing waters and should be preferred to the microplate method. The GLUase activity method should be routinely applied to the rapid enumeration of E. coli in recreational waters and recommendations for its application were suggested to water quality managers.     

Dose-dependent response of a benthic system to biodeposition from suspended blue mussel (Mytilus edulis) culture

This study reports the results of a field experiment using benthic mesocosms that examined dose-dependent effects of mussel biodeposition on the benthic environment. Mesocosms were placed in the natural sea bottom and subjected to one of eight levels of biodeposition (from 0 to 1400 mussels m^?2). Most analyses indicated non-linear (i.e., threshold) effects. Sediment characteristics changed significantly between 200 and 400 mussels m^?2 as did multivariate community structure. Community structure effects were characterised by changes in abundances of species that are very sensitive or tolerant to organic loading. The multivariate AZTI Marine Biotic Index (M-AMBI) indicated that the benthic status changed from High to Good in all mesocosms receiving biodeposits. Sediments acted as a sink for oxygen (O₂), but results suggest O₂ sediment demand was not sensitive enough to evaluate organic loading impacts. Results from this and improved experiments can be used to determine the environmental carrying capacity of sites for bivalve culture.     

Identification of Groundwater Parameters Using an Adaptative Multiscale Method

The identification of groundwater parameters in heterogeneous systems is a major challenge in groundwater modeling. Flexible parameterization methods are needed to assess the complexity of the spatial distributions of these parameters in real aquifers. In this article, we introduce an adaptative parameterization to identify the distribution of hydraulic conductivity within the large‐scale (4400 km²) Upper Rhine aquifer. The method is based on adaptative multiscale triangulation (AMT) coupled with an inverse problem procedure that identifies the parameters' distributions by reducing the error between measured and simulated heads. The AMT method has the advantage of combining both zonation and interpolation approaches. The AMT method uses area‐based interpolation rather than an interpolation based on stochastic features. The method is applied to a standard 2D groundwater model that takes into account the interactions between the aquifer and surface water bodies, groundwater recharge, and pumping wells. The simulation period covers 204 months, from January 1986 to December 2002. Recordings at 109 piezometers are used for model calibration. The simulated heads are globally quite accurate and reproduce the main dynamics of the system. The local hydraulic conductivities resulting from the AMT method agree qualitatively with existing local experimental observations across the Rhine aquifer.