Evidence of an early alteration process driven by magmatic fluid in Mururoa volcano

Evidence for a deuteric alteration process induced by a magmatic fluid has been found in the feeder zone of the Mururoa volcano (French Polynesia). Within the dikes, where basaltic glass does not show any evidence of pervasive alteration, vesicles are filled with dioctahedral smectites and calcite, while olivine phenocrysts are replaced by dioctahedral smectites, ankerite and calcite.The ¹³C signature of carbonates, the carbon and H₂O content of the whole rocks and their impoverishment in deuterium are compatible with the presence of magmatic CO₂ during the crystallization of intruding lavas and exclude contamination by seawater. Mass balance calculations on selected thin sections photographs of partly filled up vesicles and replaced olivine crystals, constrain, assuming a closed system interaction, the chemical composition of the initial fluid and the respective amounts of the initial solid phases involved in the alteration process. Thermodynamic modelings using the EQ3/6 software package correctly predict the mineralogic, chemical and isotopic exchanges accompanying alteration, thus validating the closed system assumption. The model which allows prediction of the influence of CO₂ on the alteration products, shows that, above a 0.25 CO₂ mole fraction in the initial fluid, the alteration is entirely controlled by the chemical composition of the initial solid phases. The presence of CO₂ implies the precipitation of dioctahedral smectites and carbonates instead of the magnesian smectites commonly observed in CO₂-free systems.The Mururoa feeder zone shows alteration features typical of a closed system interaction between the basaltic rock and a magmatic fluid in which seawater did not take part.     

Permeability-porosity Relationships in Rocks Subjected to Various Evolution Processes

— It is well known that there is no “universal” permeability-porosity relationship valid in all porous media. However, the evolution of permeability and porosity in rocks can be constrained provided that the processes changing the pore space are known. In this paper, we review observations of the relationship between permeability and porosity during rock evolution and interpret them in terms of creation/destruction of effectively and non-effectively conducting pore space. We focus on laboratory processes, namely, plastic compaction of aggregates, elastic-brittle deformation of granular rocks, dilatant and thermal microcracking of dense rocks, chemically driven processes, as a way to approach naturally occurring geological processes. In particular, the chemically driven processes and their corresponding evolution permeability-porosity relationships are discussed in relation to sedimentary rocks diagenesis.     

Displacements and transformations of nitrate-rich and nitrate-poor water masses in the tropical Pacific during the 1997 El Niño

A Lagrangian analysis was applied to the outputs of a coupled physical-biogeochemical model to describe the redistribution of nitrate-rich and nitrate-poor surface water masses in the tropical Pacific throughout the major 1997 El Niño. The same tool was used to analyze the causes of nitrate changes along trajectories and to investigate the consequences of the slow nitrate uptake in the high nutrient low chlorophyll (HNLC) region during the growth phase of the event. Three patterns were identified during the drift of water masses. The first mechanism is well known along the equator: oligotrophic waters from the western Pacific are advected eastward and retain their oligotrophic properties along their drift. The second concerns the persistent upwelling in the eastern basin. Water parcels have complex trajectories within this retention zone and remain mesotrophic. This study draws attention to the third process which is very specific to the HNLC region and to the El Niño period. During the 1997 El Niño, horizontal and vertical inputs of nitrate decreased so dramatically that nitrate uptake by phytoplankton became the only mechanism driving nitrate changes along pathways. The study shows that because of the slow nitrate uptake characteristic of the tropical Pacific HNLC system, nitrate in the pre-El Niño photic layer can support biological production for a period of several months. As a consequence, the slow nitrate uptake delays the gradual onset of oligotrophic conditions over nearly all the area usually occupied by upwelled waters. Owing to this process, mesotrophic conditions persist in the tropical Pacific during El Niño events.     

An automated salt‐tracing gauge for flow‐velocity measurement

This article introduces the SVG (salt‐velocity gauge), a novel automated technique for measuring flow velocity by means of salt tracing. SVG allows a high measuring rate (up to one every 2 seconds), short control section length (down to 10 cm), high accuracy (+[sol ]?1·5 cm s^?1), and unbiased calculation of the mean velocity in experimental conditions with turbulent, supercritical flow. A few cubic centimetres of saturated salt solution (NaCl) are injected into the flow at regular time intervals using a programmable solenoid valve. The tracer successively passes two conductivity probes placed a short distance downstream. The transformation of the signal between the two probes is modelled as a one‐dimensional diffusion wave equation. Model calibration gives an estimation of the mean velocity and the diffusion for each salt plume. Two implementations of the SVG technique are described. The first was an outdoors simulated rainfall experiment in Senegal (conductivity probes at 40 cm apart, 8 Hz measurement rate, salt injections at 10 second intervals). Mean velocity was estimated to range between 0·1 and 0·3 m s^?1. The second was a laboratory‐based flume experiment (conductivity probes at 10 cm apart, 32 Hz, salt injections at 2 second intervals). Another SVG with probes at 34 cm apart was used for comparison. An acoustic Doppler velocimeter (ADV) was also used to give an independent assessment of velocity. Using the 10 cm salt gauge, estimated mean velocity ranged from 0·6 to 0·9 m s^?1 with a standard deviation of 1·5 cm s^?1. Comparisons between ADV, 10 cm SVG and 34 cm SVG were consistent and demonstrated that the salt‐tracing results were unbiased and independent of distance between probes. Most peaks were modelled with r² > 90 per cent. The SVG technology offers an alternative to the dye‐tracing technique, which has been severely criticized in the literature because of the wide interval of recommended values for the correction factor α to be applied to the timings. This article demonstrates that a fixed value of α is inappropriate, since the correction factor varies with velocity, diffusion and the length of the control section. Copyright © 2005 John Wiley & Sons, Ltd.     

The Frequency Dependence of Harmonic Fluid Flow Through Networks of Cracks and Pores

—The "dynamic" permeability k(ω) of heterogeneous networks of cracks, tubes and spheres, was determined by numerically simulating the harmonic flow of an interstitial fluid for a wide range of frequencies. For comparison with previous works, this procedure was applied to the 100 network realizations used in Bernabé (1995). In most cases, the calculated frequency dependence of the real and imaginary parts of k(ω) was consistent with the JKD model (Johnson et al., 1987), showing a transition from "viscous", macroscopic flow at low frequencies to "inertial" flow at high frequencies. The viscous skin depth δ _c at the transition was found to be proportional to the critical capillary radius r _c from a capillary invasion (Katz and Thompson, 1986). A simple explanation is that these two length scales arise from the same percolation problem. On the other hand, δ _c was not well correlated with the JKD parameter Λ. The conclusion is that Λ and δ _c (or r _c ?) are two independent parameters, derived from two unrelated approaches (i.e., weighted averaging and percolation theory). Finally, an attempt was made to relax the initial assumptions of a rigid solid matrix and an incompressible fluid. It was observed that the effect of the fluid compressibility could occasionally be very large, especially when networks with large amounts of storage pore space were considered.     

Long-term groundwater resource observatory for Southwestern Madagascar

Madagascar has one of the highest poverty rates in the world and consequently the long-term monitoring of groundwater resources is not a priority for the authorities. However, groundwater is often the only sustainable resource that has a satisfactory quality to supply the population. This is especially true in the south-west of the country, which is a semi-arid region and a global change hot spot (intense land use and climate changes). In response to the lack of data, the Groundwater Resource Observatory for Southwestern Madagascar (GROSoM) was established to monitor piezometry and meteorology over the longer term as part of a humanitarian response. The first site was setup in 2014 in a catchment located over a carbonate plateau; in 2018, a second site was installed in an alluvial setting within a crystalline basement catchment and a third site will be installed in 2020 to monitor groundwater dynamics in a coastal setting. The three sites, located between Toliara and Taolagnaro cities, are complementary and representative of various hydrogeological systems in Southwestern Madagascar. Each site includes a weather station and between 3 and 6 piezometric probes. The monitoring data indicate a strong inter-annual variability in precipitation, which induces a strong variability in aquifers recharge. One of the driest years in 2016 seems to be consistent with strong El Niño – Southern Oscillation (ENSO) effects observed at the global scale, while years with higher recharge appear to be related to cyclones such as Fundi in 2015 and Eketsang in 2019. Preliminary results of cross-disciplinary studies demonstrated a link between groundwater and health issues (i.e., admissions to basic health centres). This observatory aims to produce long-term data and has two objectives: (i) strengthening the early warning system for humanitarian crises in Madagascar; (ii) contributing to a better understanding of the effects of climate change on groundwater resources in this semi-arid region.     

Comparison of LiDAR waveform processing methods for very shallow water bathymetry using Raman,near‐infrared and green signals

Airborne light detection and ranging (LiDAR) bathymetry appears to be a useful technology for bed topography mapping of non‐navigable areas, offering high data density and a high acquisition rate. However, few studies have focused on continental waters, in particular, on very shallow waters (<2 m) where it is difficult to extract the surface and bottom positions that are typically mixed in the green LiDAR signal. This paper proposes two new processing methods for depth extraction based on the use of different LiDAR signals [green, near‐infrared (NIR), Raman] of the SHOALS‐1000T sensor. They have been tested on a very shallow coastal area (Golfe du Morbihan, France) as an analogy to very shallow rivers. The first method is based on a combination of mathematical and heuristic methods using the green and the NIR LiDAR signals to cross validate the information delivered by each signal. The second method extracts water depths from the Raman signal using statistical methods such as principal components analysis (PCA) and classification and regression tree (CART) analysis. The obtained results are then compared to the reference depths, and the performances of the different methods, as well as their advantages/disadvantages are evaluated. The green/NIR method supplies 42% more points compared to the operator process, with an equivalent mean error (?4·2 cm verusu ?4·5 cm) and a smaller standard deviation (25·3 cm verusu 33·5 cm). The Raman processing method provides very scattered results (standard deviation of 40·3 cm) with the lowest mean error (?3·1 cm) and 40% more points. The minimum detectable depth is also improved by the two presented methods, being around 1 m for the green/NIR approach and 0·5 m for the statistical approach, compared to 1·5 m for the data processed by the operator. Despite its ability to measure other parameters like water temperature, the Raman method needed a large amount of reference data to provide reliable depth measurements, as opposed to the green/NIR method. Copyright © 2010 John Wiley & Sons, Ltd.     

Oscillating Flow of a Compressible Fluid through Deformable Pipes and Pipe Networks: Wave Propagation Phenomena

Similarly to blood pulse propagation in the artery system, oscillating flow can propagate as a wave in fluid-saturated pipes, networks of pipes or, by extension, in porous media, if the fluid is compressible and/or the pipes are elastically deformable. First, propagation of flow waves generated in a semi-infinite pipe by harmonic pressure oscillations at the pipe entrance is analyzed. The dispersion equation is derived, allowing determination of the phase velocity and quality factor as functions of frequency. Wave reflections at the end of a finite-length pipe and ensuing interferences between forward and backward traveling waves are then examined. Because of fluid storage in the pipe, the amplitude of the AC volumetric fluxes entering and exiting the pipe at its upstream and downstream ends are not equal. Thus, two different, upstream and downstream, frequency-dependent, AC hydraulic conductivities are introduced. Superposed on the classic viscous-inertial flow transition (controlled by the value of the pipe radius), these complex-valued parameters show another transition between an interference-free regime at low frequencies and a strong interference regime above a critical frequency that roughly scales as the pipe length. Because of attenuation, the flow wave interferences tend to gradually weaken with increasing frequencies. Finally, the single pipe model is used to investigate fluid flow waves through pipe networks with results very similar to those described above. The flow waves analyzed here are akin to the Biot’s slow P waves and their propagation properties could affect seismic soundings in some geological settings.     

The snowball Earth aftermath: Exploring the limits of continental weathering processes

Carbonates capping Neoproterozoic glacial deposits contain peculiar sedimentological features and geochemical anomalies ascribed to extraordinary environmental conditions in the snowball Earth aftermath. It is commonly assumed that post-snowball climate dominated by CO₂ partial pressures several hundred times greater than modern levels, would be characterized by extreme temperatures, a vigorous hydrological cycle, and associated high continental weathering rates. However, the climate in the aftermath of a global glaciation has never been rigorously modelled. Here, we use a hierarchy of numerical models, from an atmospheric general circulation model to a mechanistic model describing continental weathering processes, to explore characteristics of the Earth system during the supergreenhouse climate following a snowball glaciation. These models suggest that the hydrological cycle intensifies only moderately in response to the elevated greenhouse. Indeed, constraints imposed by the surface energy budget sharply limit global mean evaporation once the temperature has warmed sufficiently that the evaporation approaches the total absorbed solar radiation. Even at 400 times the present day pressure of atmospheric CO₂, continental runoff is only 1.2 times the modern runoff. Under these conditions and accounting for the grinding of the continental surface by the ice sheet during the snowball event, the simulated maximum discharge of dissolved elements from continental weathering into the ocean is approximately 10 times greater than the modern flux. Consequently, it takes millions of years for the silicate weathering cycle to reduce post-snowball CO₂ levels to background Neoproterozoic levels. Regarding the origin of the cap dolostones, we show that continental weathering alone does not supply enough cations during the snowball melting phase to account for their observed volume.     

Subsurface DC resistivity mapping: approximate 1-D interpretation

Resistivity prospecting is the main tool used to investigate the shallow structure of the ground. A series of new techniques for determining the 2-D and 3-D geometry of the ground is now finding increasing use, but the light and simple Wenner prospecting technique remains a practical and efficient tool for rapidly mapping lateral variations in resistivity. When the resistivity changes are smooth, 1-D modelling can be used to interpret the data, and the criteria governing this approximation can be defined from synthetic data generated by a 3-D slab-model. For a Wenner array, two quadripole configurations can be used, Normal and Dipole-Dipole. For these two configurations the width of the transition zone, the apparent anisotropy effect and the precision of the resistivity values recovered from 1-D inversion differ. However the simultaneous inversion of both sets of data gives better results than for either configuration by itself. Two examples illustrate that in geological contexts where the thickness of the weathered zone causes the changes in the apparent resistivity value, this parameter can be recovered from 1-D inversion.