Montagne Pelée volcano (Martinique,in the French Lesser Antilles) hydrogeological system revealed by high-resolution helicopter-borne electromagnetic imagery

Montagne Pelée, on the French island of Martinique, eastern Caribbean Sea, has been one of the deadliest volcanoes in the world, with 30,000 victims following the 1902 eruption. Thousands of people still live nearby, and this volcano is a strategic “water tank” for Martinique Island, providing 40% of the island’s water supply. This research aimed to better understand its hydrogeological functioning and the relationship with its complex volcanological evolution, taking advantage of a high-resolution helicopter-borne geophysical survey correlated with hydrogeological data from the boreholes and springs databases. Electromagnetic data, correlated with hydrogeological data, allowed for the identification of unsaturated zones, aquifers, and seawater intrusions, as well as the main geological units. In addition, data synthesised from pumping tests revealed that the older the unconsolidated pyroclastic deposits, the lower their hydraulic conductivity. The structural asymmetry between the northeastern and southwestern volcano flanks impacts its hydrogeological functioning. Consequently, the Montagne Pelée hydrogeological conceptual model is marked by several distinguishable aquifers. The upper perched aquifer within recent lava domes is directly involved in, and impacted by, phreatic eruptions, and it supports low flowrate springs. The remaining effective rainfall infiltrates to depth and recharges the hydrothermal system through vertical fractures. The other aquifers are categorized into three groups: northeastern, southeastern and southwestern flank aquifers. This research is a new step toward a better understanding of the Lesser Antilles volcanoes and more broadly of the central and proximal parts of the andesitic active volcanoes.

     

A new method to estimate the oxidation state of basaltic series from microprobe analyses

The oxygen fugacity and therefore the iron redox state of a melt is known to have a strong influence on the liquid line of descent of magmas and thus on the composition of the coexisting melts and crystals. We present a new method to estimate this critical parameter from electron probe microanalyses of two of the most common minerals of basaltic series, plagioclase and clinopyroxene. This method is not based on stoichiometric calculations, but on the different partitioning behaviour of Fe³⁺ and Fe²⁺ between both minerals and a melt phase: plagioclase can incorporate more Fe³⁺ than Fe²⁺, while clinopyroxene can incorporate more Fe²⁺ than Fe³⁺. For example, the effect of oxidizing a partly molten basaltic system (Fe³⁺ is stabilized with respect to Fe²⁺) results in an increase of FeO_total in plagioclase, but a decrease in the associated clinopyroxene. We propose an equation, based on published partition coefficients, that allows estimating the redox state of a melt from these considerations. An application to a set of experimental and natural data attests the validity of the proposed model. The associated error can be calculated and is on average < 1 log unit of the prevailing oxygen fugacity.In order to reduce the number of different variables influencing the Fe²⁺/Fe³⁺ mineral/melt equilibrium, our model is restricted to basaltic series with SiO₂ < 60% that have crystallized at intermediate to low pressure (< 0.5 GPa) and under relatively oxidizing conditions (?FMQ > 0; where FMQ is the fayalite–magnetite–quartz oxygen buffer equilibrium), but it may be parameterized for other conditions. A spreadsheet is provided to assist the use of equations, and to perform the error propagation analysis.     

The significance of stretching lineations in terms of progressive deformation and finite strain

Abstract

For a given increment of deformation, stretching lineations can be defined either by old grains or by new grains.

Progressive strain accumulated in a matrix causes reorientation of old grains. The velocity field acting within a grain depends on the degree of non-coaxiality of the progressive strain in the matrix, and on the competence contrast between grain and matrix (e.g. Lister & Williams, 1983). Deformed populations of old grains do not systematically track the principal directions of finite stretch; this depends on the mechanical behavior of the grains. On the other hand, these populations may be good markers of finite strain magnitude, provided that the initial distribution is known.

In contrast, the statistical shape orientation distribution of population of new grains cannot be related in a simple way to the finite strain magnitude, but their preferred orientation is considered to be the best indicator of the principal directions of finite stretch.

An example is considered from the schistes lustrés near Kenestrelle (Val Chisone, Western Alps), where the observed mineral lineation is attributed to a late E-W stretching event in the alpine orogenic history.     

Apparent flow-dimension approach to the study of heterogeneous fracture network systems

Hydrogeology Journal - The generalized radial flow (GRF) model in well-test analysis employs noninteger flow dimensions to represent the variation in flow area with respect to radial distance from...     

A comparison of approaches to include outcrop information in overburden thickness estimation

The estimation of overburden sediment thickness is important in hydrogeology, geotechnics and geophysics. Usually, thickness is known precisely at a few sparse borehole data. To improve precision of estimation, one useful complementary information is the known position of outcrops. One intuitive approach is to code the outcrops as zero thickness data. A problem with this approach is that the outcrops are preferentially observed compared to other thickness information. This introduces a strong bias in the thickness estimation that kriging is not able to remove. We consider a new approach to incorporate point or surface outcrop information based on the use of a non-stationary covariance model in kriging. The non-stationary model is defined so as to restrict the distance of influence of the outcrops. Within this distance of influence, covariance parameters are assumed simple regular functions of the distance to the nearest outcrop. Outside the distance of influence of the outcrops, the thickness covariance is assumed stationary. The distance of influence is obtained thru a cross-validation. Compared to kriging based on a stationary model with or without zero thickness at outcrop locations, the non-stationary model provides more precise estimation, especially at points close to an outcrop. Moreover, the thickness map obtained with the non-stationary covariance model is more realistic since it forces the estimates to zero close to outcrops without the bias incurred when outcrops are simply treated as zero thickness in a stationary model.     

Assessment of the Design of an Experimental Cover with Capillary Barrier Effect Using 4 Years of Field Data

One important step in the design of inclined covers with capillary barrier effect (CCBE) is the determination of the water diversion length (DL). Numerical simulations can predict the DL more precisely than steady-state analytical solutions. Nevertheless, as simplified methods have always been part of engineering design, the application of analytical solutions with conservative boundary conditions, may allow engineers to make reasonable predictions, particularly during the pre-feasibility stage of a project. In this study, a CCBE was designed, constructed and instrumented at the Saint-Tite-des-Caps landfill, Quebec, Canada. This CCBE included a seepage control layer superimposing a sand-gravel capillary barrier. The seepage control layer was made up of deinking by-products (DBP), an industrial by-product that was previously disposed of as waste. The capillary barrier was designed using an adaptation of the Ross analytical solution and the scenario considered was that of steady-state flow during constant seepage flow applied uniformly at the top of the sand-gravel capillary barrier. Although these conditions appear simplistic, they were deemed reasonable because placement of the seepage control layer on the top of the capillary barrier led to very low suctions at the interface, thereby allowing uniform downward seepage rates, limited by the saturated hydraulic conductivity of the DBP. In this paper, a discussion about the behaviour of the cover system based on 4 years of field data from several instruments is presented. The challenge of using DBP, more precisely the settlement of the DBP layer and its impact on k _sat , is also assessed. The DL was reassessed considering the new k _sat . A discussion on the validity of employing analytical solutions to determine DL is also presented. This paper illustrates how certain variables affect the design of inclined CCBEs that include a highly compressible material as seepage control layer.     

Volcanic ash hazard in the Central Mediterranean assessed from geological data

Volcanic ash produced during explosive eruptions can have very severe impacts on modern technological societies. Here, we use reconstructed patterns of fine ash dispersal recorded in terrestrial and marine geological archives to assess volcanic ash hazards. The ash-dispersal maps from nine Holocene explosive eruptions of Italian volcanoes have been used to construct frequency maps of distal ash deposition over a wide area, which encompasses central and southern Italy, the Adriatic and Tyrrhenian seas and the Balkans. The maps are presented as two cumulative-thickness isopach maps, one for nine eruptions from different volcanoes and one for six eruptions from Somma-Vesuvius. These maps represent the first use of distal ash layers to construct volcanic hazard maps, and the proposed methodology is easily applicable to other volcanic areas worldwide.     

Investigating possible links between the North Atlantic Oscillation and rainfall variability in Marrakech (Morocco)

Eighty-two-year rainfall time series have been studied together with climatic patterns of NAO using classical statistical methods. Then, the wavelet approach has been applied to show annual (1 year (1y)) and inter-annual (2–4 years (2-4y), 5–8 years (5-8y), and 8–16 years (8-16y)) modes distributed following four major discontinuities: 1945, 1960, 1975, and 1995. The 1y, 2-4y, and 5-8y powers show high energy during the wet period 1922–1930 and a low one in 1928–1938. After 1945, the annual mode highlights a high energy while the inter-annual modes present low energy. Between 1975 and 1995, powerful modes of 1 and 2–4 years are identified with low power of 5-8y and 8-16y modes. Since 1995, the low power of 5–8y decreases, while the 8-16y mode emphasizes a high variability. The coherence between NAO and Marrakech precipitation is strongly defined for low frequencies with a total contribution of 75 %. This coherence is in phase in the beginning and presents out phase signs since 1945. The change of phase can be associated to a decreasing of coherence especially around 1990. This finding is useful to understand the relationship between the hydrological variability and NAO climate patterns in the southern side of the Mediterranean Sea.