Fault zone geometry of a mature active normal fault: A potential high permeability channel (Pirgaki fault, Corinth rift, Greece)

We present results from petrophysical analysis of a normal fault zone with the aim of defining the flow pathways and their behavior during seismic and interseismic periods. Data are obtained on porosity geometry, strain structure and mineralogy of different domains of a normal fault zone in the Corinth rift. Data point out a close relationship between mineralogy of the clayey minerals, porosity network and strain structures and allow definition of a macroscopic anisotropy of the flow parameters with a strong control by microscopic ultracataclasite structures. The Pirgaki fault zone, developed within pelagic limestone, has a sharp asymmetric porosity profile, with a high porosity volume in the fault core and in the damage zone of the hanging wall. From porosity volumes and threshold measurements, a matrix permeability variation of 6 orders of magnitude could be expected between the protolith and the fault core. Modifications of this pathway during seismic and interseismic phases are depicted. Healing of cracks formed during seismic slip events occurred in the fault core zone and the porous network in the damage zone is sealed in a second step. The lens geometry of the fault core zone is associated with dissolution surfaces and open conduits where dissolved matter could move out of the fault core zone. These elementary processes are developed in particular along Riedel's structures and depend on the orientation of the strain surfaces relative to the local stress and depend also on the roughness of each surface type. P-surfaces are smooth and control shearing process. R-surfaces are rough and present two wavelengths of roughness. The long one controls localization of dissolution surfaces and conduits; the short one is characteristic of dissolution surfaces. The dissolved matter can precipitate in the open structures of the hanging wall damage zone, decreasing the connectivity of the macroscopic conduit developed within this part of the fault zone.     

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.     

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.     

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.     

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.     

Field characterization and data integration to define the hydraulic heterogeneity of a shallow granular aquifer at a sub-watershed scale

Providing a sound basis for aquifer management or remediation requires that hydrogeological investigations carried out to understand groundwater flow and contaminant transport be based on representative data that capture the heterogeneous spatial distribution of aquifer hydraulic properties. This paper describes a general workflow allowing the characterization of the heterogeneity of the hydraulic properties of granular aquifers at an intermediate scale of a few km². The workflow involves characterization and data integration steps that were applied on a 12-km² study area encompassing a decommissioned landfill emitting a leachate plume and its main surface water receptors. The sediments composing the aquifer were deposited in a littoral–sublittoral environment and show evidence of small-scale transitional heterogeneities. Cone penetrometer tests (CPT) combined with soil moisture and electrical resistivity (SMR) measurements were thus used to identify and characterize spatial heterogeneities in hydraulic properties over the study area. Site-specific statistical relationships were needed to infer hydrofacies units and to estimate hydraulic properties from high-resolution CPT/SMR soundings distributed all over the study area. A learning machine approach was used due to the complex statistical relationships between colocated hydraulic and CPT/SMR data covering the full range of aquifer materials. Application of this workflow allowed the identification of hydrofacies units and the estimation of horizontal hydraulic conductivity, vertical hydraulic conductivity and porosity over the study area. The paper describes and discusses data acquisition and integration methodologies that can be adapted to different field situations, while making the aquifer characterization process more time-efficient and less labor-intensive.     

Mineralisation of amethyst-bearing geodes in Ametista do Sul (Brazil) from low-temperature sedimentary brines: evidence from monophase liquid inclusions and stable isotopes

Fluid inclusion studies in combination with hydrogen, oxygen and sulphur isotope data provide novel insights into the genesis of giant amethyst-bearing geodes in Early Cretaceous Paraná continental flood basalts at Amestita do Sul, Brazil. Monophase liquid inclusions in colourless quartz, amethyst, calcite, barite and gypsum were analysed by microthermometry after stimulating bubble nucleation using single femtosecond laser pulses. The salinity of the fluid inclusions was determined from ice-melting temperatures and a combination of prograde and retrograde homogenisation temperatures via the density maximum of the aqueous solutions. Four mineralisation stages are distinguished. In stage I, celadonite, chalcedony and pyrite formed under reducing conditions in a thermally stable environment. Low δ³⁴S_V-CDT values of pyrite (?25 to ?32?‰) suggest biogenic sulphate reduction by organotrophic bacteria. During the subsequent stages II (amethyst, goethite and anhydrite), III (early subhedral calcite) and IV (barite, late subhedral calcite and gypsum), the oxidation state of the fluid changed towards more oxidising conditions and microbial sulphate reduction ceased. Three distinct modes of fluid salinities around 5.3, 3.4 and 0.3 wt% NaCl-equivalent characterise the mineralisation stages II, III and IV, respectively. The salinity of the stage I fluid is unknown due to lack of fluid inclusions. Variation in homogenisation temperatures and in δ¹⁸O values of amethyst show evidence of repeated pulses of ascending hydrothermal fluids of up to 80–90 °C infiltrating a basaltic host rock of less than 45 °C. Colourless quartz and amethyst formed at temperatures between 40 and 80 °C, while the different calcite generations and late gypsum precipitated at temperatures below 45 °C. Calculated oxygen isotope composition of the amethyst-precipitating fluid in combination with δD values of amethyst-hosted fluid inclusions (?59 to ?51?‰) show a significant ¹⁸O-shift from the meteoric water line. This ¹⁸O-shift, high salinities of the fluid inclusions with chloride-sulphate composition, and high δ³⁴S values of anhydrite and barite (7.5 to 9.9?‰) suggest that sedimentary brines from deeper parts of the Guaraní aquifer system must have been responsible for the amethyst mineralisation.     

Fluvial incisions in the North‐Western Pyrenees (Aspe Valley): Dissection of a former planation surface and some tectonic implications

Dissected remnants of a high elevation low relief surface (HELRS) are encountered all along the Pyrenean range. All authors agree that the HELRS was shaped during the post‐shortening evolution of the belt, but doubt remains on its original elevation. Notably, whether a post‐shortening uplift event occurred after the generation of the Pyrenean planation surface is still debated. Based on a geomorphological study of the entrenchment of the Aspe River in the North‐Western Pyrenees, we describe a post‐Oligocene multi‐stage dissection of the HELRS linked to a regional post‐shortening uplift. Each incision stage is recorded by erosional triangular facets and associated stepped remnants of former smooth topographies. Compared analysis of patterns of incision in the Axial Zone and its northern border allows to evidence differential vertical motions. We discuss the forcing mechanism(s) that controlled this morphological evolution.