Data report: permeability measurements under confining pressure, Expeditions 315 and 316, Nankai Trough.

International audiencePermeability of six samples from sites C0001 and C0006 were measured in a triaxial cell under effective hydrostatic confining pressure from 1 to 30 MPa. Our results indicate that the initial permeability at 1 MPa of effective confining pressure ranges from 4.6e-18 to 1.8e-19 m2 depending on depth. Actually permeability decreases with increasing depth also corresponding to a decrease of porosity from 62 to 43%. The permeability vs. depth trend is similar for both sites. When the effective confining pressure is increased from 1 to 30 MPa, the permeability decreases for all samples, a decrease interpreted by microfracture closure. However this trend shows some variability indicating a finer microstructural control depending on the lithological origin of the sample

The permeability and elastic moduli of tuff from Campi Flegrei, Italy: implications for ground deformation modelling

The accuracy of ground deformation modelling at active volcanoes is a principal requirement in volcanic hazard mitigation. However, the reliability of such models relies on the accuracy of the rock physical property (permeability and elastic moduli) input parameters. Unfortunately, laboratory-derived values on representative rocks are usually rare. To this end we have performed a systematic laboratory study on the influence of pressure and temperature on the permeability and elastic moduli of samples from the two most widespread lithified pyroclastic deposits at the Campi Flegrei volcanic district, Italy. Our data show that the water permeability of Neapolitan Yellow Tuff and a tuff from the Campanian Ignimbrite differ by about 1.5 orders of magnitude. As pressure (depth) increases beyond the critical point for inelastic pore collapse (at an effective pressure of 10–15 MPa, or a depth of about 750 m), permeability and porosity decrease significantly, and ultrasonic wave velocities and dynamic elastic moduli increase significantly. Increasing the thermal stressing temperature increases the permeability and decreases the ultrasonic wave velocities and dynamic elastic moduli of the Neapolitan Yellow Tuff; whereas the tuff from the Campanian Ignimbrite remains unaffected. This difference is due to the presence of thermally unstable zeolites within the Neapolitan Yellow Tuff. For both rocks we also find, under the same pressure conditions, that the dynamic (calculated from ultrasonic wave velocities) and static (calculated from triaxial stress-strain data) elastic moduli differ significantly. The choice of elastic moduli in ground deformation modelling is therefore an important consideration. While we urge that these new laboratory data should be considered in routine ground deformation modelling, we highlight the challenges for ground deformation modelling based on the heterogeneous nature (vertically and laterally) of the rocks that comprise the caldera at Campi Flegrei

The Influence of Grain Size Distribution on Mechanical Compaction and Compaction Localization in Porous Rocks

The modes of formation of clastic rocks result in a wide variety of microstructures, from poorly-sorted heterogeneous rocks to well-sorted and nominally homogeneous rocks. The mechanical behavior and failure mode of clastic rocks is known to vary with microstructural attributes such as porosity and grain size. However, the influence of the grain size distribution, in particular the degree of polydispersivity or modality of the distribution, is not yet fully understood, because it is difficult to study experimentally using natural rocks. To better understand the influence of grain size distribution on the mechanical behavior of porous rocks, we prepared suites of synthetic samples consisting of sintered glass beads with polydisperse grain size distributions. We performed hydrostatic compression experiments and found that, all else being equal, the onset of grain crushing occurs much more progressively and at lower pressure in polydisperse synthetic samples than in monodisperse samples. We conducted triaxial experiments in the regime of shear-enhanced compaction and found that the stress required to reach inelastic compaction was lower in polydisperse samples compared to monodisperse samples. Further, our microstructural observations show that compaction bands developed in monomodal polydisperse samples while delocalized cataclasis developed in bimodal polydisperse samples, where small grains were systematically crushed while largest grains remained intact. In detail, as the polydispersivity increases, microstructural deformation features appear to transition from localized to delocalized through a hybrid stage where a compaction front with diffuse bands propagates from both ends of the sample toward its center with increasing bulk strain

U–Th and ¹⁰Be constraints on sediment recycling in proglacial settings, Lago Buenos Aires, Patagonia

International audienceThe estimation of sediment transfer times remains a challenge to our understanding of sediment budgets and the relationships between erosion and climate. Uranium (U) and thorium (Th) isotope disequilibria offer a means of more robustly constraining sediment transfer times. Here, we present new uranium and tho-rium disequilibrium data for a series of nested moraines around Lago Buenos Aires in Argentine Patagonia. The glacial chronology for the area is constrained using in situ cosmogenic 10 Be analysis of glacial outwash. Sediment transfer times within the periglacial domain were estimated by comparing the deposition ages of moraines to the theoretical age of sediment production, i.e., the comminution age inferred from U disequilibrium data and recoil loss factor estimates. Our data show first that the classical comminution age approach must include weathering processes accounted for by measuring Th disequilibrium. Second, our combined data suggest that the pre-deposition history of the moraine sediments is not negligible, as evidenced by the large disequilibrium of the youngest moraines despite the equilibrium of the corresponding glacial flour. Monte Carlo simulations suggest that weathering was more intense before the deposition of the moraines and that the transfer time of the fine sediments to the moraines was on the order of 100-200 kyr. Long transfer times could result from a combination of long sediment residence times in the proglacial lake (recurrence time of a glacial cycle) and the remobilization of sediments from moraines deposited during previous glacial cycles. 10 Be data suggest that some glacial cycles are absent from the preserved moraine record (seemingly every second cycle), supporting a model of reworking moraines and/or fluctuations in the extent of glacial advances. The chronological pattern is consistent with the U-Th disequilibrium data and the 100-200 kyr transfer time. This long transfer time raises the question of the proportion of freshly eroded sediments that escape (or not) the proglacial environments during glacial periods

Interpretation of porosity and LWD resistivity from the Nankai accretionary wedge in light of clay physicochemical properties: Evidence for erosion and local overpressuring

International audienceIn this study, we used porosity to assess the compaction state of the Nankai accretionary wedge sediments and any implications for stress and pore pressure. However, hydrous minerals affect porosity measurements, and accounting for them is essential toward defining the interstitial porosity truly representative of the compaction state. The water content of sediments was measured in core samples and estimated from logging data using a resistivity model for shale. We used the cation exchange capacity to correct the porosity data for the amount of water bound to clay minerals and to correct the porosity estimates for the surface conductivity of hydrous minerals. The results indicate that several apparent porosity anomalies are significantly reduced by this correction, implying that they are in part artifacts from hydrous minerals. The correction also improves the fit of porosity estimated from logging-while-drilling (LWD) resistivity data to porosity measured on cores. Low overall porosities at the toe of the accretionary wedge and in the splay fault area are best explained by erosion, and we estimated the quantity of sediments eroded within the splay fault area by comparing porosity-effective stress relationships of the sediments to a reference curve. Additionally, a comparison of LWD data with core data (resistivity and P wave velocity) obtained at Site C0001 landward of the mega-splay fault area, suggested a contribution from the fracture porosity to in situ properties on the formation

Multidisciplinary constraints of hydrothermal explosions based on the 2013 Gengissig lake events, Kverkfjöll volcano, Iceland

Highlights • A multidisciplinary approach to unravel the energetics of hydrothermal explosions. • Pressure failure caused by a lake drainage triggered the hydrothermal explosions. • Bedrock nature controlled the explosion dynamics and the way energy was released. • Approx. 30% of the available thermal energy is converted into mechanical energy. • Released seismic energy as proxy to detect past (and future?) hydrothermal explosions. Hydrothermal explosions frequently occur in geothermal areas showing various mechanisms and energies of explosivity. Their deposits, though generally hardly recognised or badly preserved, provide important insights to quantify the dynamics and energy of these poorly understood explosive events. Furthermore the host rock lithology of the geothermal system adds a control on the efficiency in the energy release during an explosion. We present results from a detailed study of recent hydrothermal explosion deposits within an active geothermal area at Kverkfjöll, a central volcano at the northern edge of Vatnajökull. On August 15th 2013, a small jökulhlaup occurred when the Gengissig ice-dammed lake drained at Kverkfjöll. The lake level dropped by approximately 30 m, decreasing pressure on the lake bed and triggering several hydrothermal explosions on the 16th. Here, a multidisciplinary approach combining detailed field work, laboratory studies, and models of the energetics of explosions with information on duration and amplitudes of seismic signals, has been used to analyse the mechanisms and characteristics of these hydrothermal explosions. Field and laboratory studies were also carried out to help constrain the sedimentary sequence involved in the event. The explosions lasted for 40–50 s and involved the surficial part of an unconsolidated and hydrothermally altered glacio-lacustrine deposit composed of pyroclasts, lavas, scoriaceous fragments, and fine-grained welded or loosely consolidated aggregates, interbedded with clay-rich levels. Several small fans of ejecta were formed, reaching a distance of 1 km north of the lake and covering an area of approximately 0.3 km2, with a maximum thickness of 40 cm at the crater walls. The material (volume of approximately 104 m3) has been ejected by the expanding boiling fluid, generated by a pressure failure affecting the surficial geothermal reservoir. The maximum thermal, craterisation and ejection energies, calculated for the explosion areas, are on the order of 1011, 1010 and 109 J, respectively. Comparison of these with those estimated by the volume of the ejecta and the crater sizes, yields good agreement. We estimate that approximately 30% of the available thermal energy was converted into mechanical energy during this event. The residual energy was largely dissipated as heat, while only a small portion was converted into seismic energy. Estimation of the amount of freshly-fragmented clasts in the ejected material obtained from SEM morphological analyses, reveals that a low but significant energy consumption by fragmentation occurred. Decompression experiments were performed in the laboratory mimicking the conditions due to the drainage of the lake. Experimental results confirm that only a minor amount of energy is consumed by the creation of new surfaces in fragmentation, whereas most of the fresh fragments derive from the disaggregation of aggregates. Furthermore, ejection velocities of the particles (40–50 m/s), measured via high-speed videos, are consistent with those estimated from the field. The multidisciplinary approach used here to investigate hydrothermal explosions has proven to be a valuable tool which can provide robust constraints on energy release and partitioning for such small-size yet hazardous, steam-explosion events

KG2B, a collaborative benchmarking exercise for estimating the permeability of the Grimsel granodiorite - Part 1: Measurements, pressure dependence and pore-fluid effects

Measuring the permeability of tight rocks remains a challenging task. In addition to the traditional sources of errors that affect more permeable formations (e.g. sample selection, non-representative specimens, disturbance introduced during sample acquisition and preparation), tight rocks can be particularly prone to solid–fluid interactions and thus more sensitive to the methods, procedures and techniques used to measure permeability. To address this problem, it is desirable to collect, for a single material, measurements obtained by different methods and pore-fluids. For that purpose a collaborative benchmarking exercise involving 24 laboratories was organized for measuring the permeability of a single low permeability material, the Grimsel granodiorite, at a common effective confining pressure (5 MPa). The objectives of the benchmark were: (i) to compare the results for a given method, (ii) to compare the results between different methods, (iii) to analyze the accuracy of each method, (iv) to study the influence of experimental conditions (especially the nature of pore fluid), (v) to discuss the relevance of indirect methods and models and finally (vi) to suggest good practice for low permeability measurements. In total 39 measurements were collected that allowed us to discuss the influence of (i) pore-fluid, (ii) measurement method, (iii) sample size and (iv) pressure sensitivity. Discarding some outliers from the bulk data set (4 out of 39) an average permeability of 1.11 × 10−18 m² with a standard deviation of 0.57 × 10−18 m² was obtained. The most striking result was the large difference in permeability for gas measurements compared to liquid measurements. Regardless of the method used, gas permeability was higher than liquid permeability by a factor approximately 2 (kgas = 1.28 × 10−18 m² compared to kliquid = 0.65 × 10−18 m²). Possible explanations are that (i) liquid permeability was underestimated due to fluid-rock interactions (ii) gas permeability was overestimated due to insufficient correction for gas slippage and/or (iii) gases and liquids do not probe exactly the same porous networks. The analysis of Knudsen numbers shows that the gas permeability measurements were performed in conditions for which the Klinkenberg correction is sufficient. Smaller samples had a larger scatter of permeability values, suggesting that their volume were below the Representative Elementary Volume. The pressure dependence of permeability was studied by some of the participating teams in the range 1–30 MPa and could be fitted to an exponential law k = ko.exp(–γPeff) with γ = 0.093 MPa−1. Good practice rules for measuring permeability in tight materials are also provided