EMRP1.14

Sedimentary crustal porous rocks span a wide range of grain size distributions – from monodisperse to highly polydisperse. The distribution of grain size depends on the location and conditions of rock formation, the chemico-physical processes at play, and is influenced by subsequent geological processes. Well-sorted granular rocks, with a grain size distribution close to monodisperse, and granular rocks with a more polydisperse grain size distribution, have repeatedly been subjected to laboratory experiments. And yet the natural variability from sample to sample and structural heterogeneity within single natural samples all conspire to prevent us from constraining the effect of grain size polydispersivity. While a few studies have focused on the influence of grain size, the control of grain size distribution on the mechanical behavior of rocks has scarcely been studied, especially in the laboratory. In this study, we address this knowledge-gap using synthetic samples prepared by sintering glass beads with controlled polydisperse grain size distributions. When heated above the glass transition temperature, the beads act as viscous droplets and sinter together. Throughout viscous sintering, a bead pack evolves from an initial granular discontinuous state into a solid connected porous state, at which the microstructural geometries and final porosity are known. Variably polydisperse individual samples were prepared by mixing glass beads with diameters of 0.2, 0.5, and 1.15 mm in various proportions, which were sintered together to a final porosity of 0.25 or 0.35. Hydrostatic and triaxial compression experiments were performed for each combination of polydispersivity. The samples were water-saturated, deformed at room temperature, and deformed under drained conditions (with a fixed pore pressure of 10 MPa). Triaxial experiments were conducted at a constant strain rate at effective pressure corresponding to the ductile (compactive) regime. Our mechanical data provide evidence that polydispersivity exerts a significant control on the compactive behavior of porous rocks. Insights into the microstructure were gained using scanning electron microscopy on thin sections prepared from samples before and after deformation. These data allow for the observation of the different deformation features, and by extension the deformation micro-mechanisms, promoted by the different type and degree of polydispersivity. Overall, our data show that, at a fixed porosity, increasing polydispersivity decreases the stress required for compactant failure.

How to cite: Carbillet, L., Heap, M., Wadsworth, F., Baud, P., and Reuschlé, T.: From monodisperse to polydisperse: the influence of grain size distribution on the mechanical behavior of porous synthetic rocks, EGU General Assembly 2021, online, 19–30 Apr 2021, EGU21-2629, https://doi.org/10.5194/egusphere-egu21-2629, 2021.

The localisation of structural damage, in the form of faults and fractures, along a distinct and emergent fault plane is the key driving mechanism for catastrophic failure in the brittle Earth. However, due to the speed at which stable crack growth transitions to dynamic rupture, the precise mechanisms involved in localisation as a pathway to fault formation remain unknown. Understanding these mechanisms is critical to understanding and forecasting earthquakes, including induced seismicity, landslides and volcanic eruptions, as well as failure of man-made materials and structures. We used time-resolved synchrotron x-ray microtomography to image in-situ damage localisation at the micron scale and at bulk axial strain rates down to 10^-7 s^-1. By controlling the rate of micro-fracturing events during a triaxial deformation experiment, we deliberately slowed the strain localisation process from seconds to minutes as failure approached. This approach, originally established to indirectly image fault nucleation and propagation with acoustic emissions, is completely novel in synchrotron x-ray microtomography and has enabled us to image directly processes that are normally too transient even for fast synchrotron imaging methods. Here, we first present the experimental apparatus and control system used to acquire the data, followed by damage localisation and shear zone development in a sample of Clashach sandstone viewed in unprecedented detail. Time-resolved microtomography images demonstrate a strong intrinsic correlation between shear and dilatant strain in the localised zone, with bulk shear strain accomodated by the nucleation and rotation of en-echelon tensile microcracks within a grain-scale shear band. Rotation is accompanied by antithetic to synthetic shear sliding of neighbouring crack surfaces as they rotate. The evolving 4D strain field, measured with incremental digital volume correlation between pairs of recorded x-ray tomographic volumes, independently confirm the correlation between shear and dilatant strain and show how strain localises spontaneously, first through exploration of several competing shear bands at peak stress before transitioning to failure along the optimally-oriented final fault plane. In order to ‘ground-truth’ inferences made from bulk measurements and seismic waves (the primary method of detecting deformation at the field-scale where direct imaging of the subsurface is impossible), we (a) compare rupture energy estimates from local slip measurements with those from bulk slip data, and (b) use AE source location estimates to identify individual cracks and other local changes in the microstucture that may explain the AE source.

How to cite: Cartwright-Taylor, A., Main, I. G., Butler, I. B., Fusseis, F., Mangriotis, M.-D., Curtis, A., Bell, A., Ling, M., Andò, E., Rizzo, R., Marti, S., Leung, D., and Magdysyuk, O.: Seeing and hearing quasi-static shear band localization in a sandstone, EGU General Assembly 2021, online, 19–30 Apr 2021, EGU21-9048, https://doi.org/10.5194/egusphere-egu21-9048, 2021.

The Jharia region of lower Gondwana in India is one of the largest Underground Coalmine Fires (UCF) affected coalfield in the world. The UCF induced small scale as well as large-scale surface fracturing often creates the life-threatening conditions to coal miners and local surroundings. So, there is a need to understand the thermomechanical behaviour of coal measures rocks to predict the land disturbances in thermo-environmental conditions. It will provide an insight into the UCF induced subsidence mechanism and its preventive measures. The Jharia coal field predominantly consists of sandstone (75-80% by volume) and rest is composed of coal, shale and carbonaceous shale. The present study focuses on thermo-mechanical behaviour of Barakar sandstone (BS) under elevated temperatures. The cores of BS sample were prepared according to the ISRM standards. Further, samples were grouped and thermally treated in temperature range of 25°C, 100°C, 150°C, 300°C, 400°C, 500°C, 600 °C, 700°C and 800°C at a heating rate of 5°C/min for 24 hours in furnace.  Then, these thermally treated BS samples were subjected to laboratory test for stress-strain characteristics. In the process of deformational characteristics evaluation, effect of mineralogical changes and mode of fracture pattern were also studied at the mentioned elevated temperature. Based on the obtained results, the deformational behaviour of thermally treated BS specimens can be grouped into three zones, viz., zone 1 (25-300°C), zone 2 (300-500°C) and zone 3 (500-800°C). In zone 1, the characteristics of the stress-strain curve is similar to those under air dried sandstone specimen. However, small increment in stiffness were observed upto 300°C. The stress-strain curves in this zone shows dominantly brittle fracturing. The increment in stiffness may be related to evaporation of pore water that increases the cohesion between the mineral grains resulting higher stiffness value. In zone 2, the deformation pattern again shows brittle fracturing with continuous decrement in stiffness. The reduction in stiffness may be related to thermally induced porosity and increased microcrack density. In zone 3, the stress strain curve is observed to be concave upward. It indicates the pseudo-ductile behaviour of the thermally treated BS specimens. The observed results suggest a typical behaviour of deformation pattern under UCF induced rock fracturing which may be useful in predicting the land subsidence in UCF affected areas. Present research outcome may be used to design the support measures to reduce the associated hazards.

How to cite: Tripathi, A., Gupta, N., Singh, A. K., Rai, N., and Pain, A.: Deformational characteristics of thermally treated sandstone from an underground coalmine fire region, India, EGU General Assembly 2021, online, 19–30 Apr 2021, EGU21-14332, https://doi.org/10.5194/egusphere-egu21-14332, 2021.

Not all rocks are perfect. Frequently heterogeneities will be present, either in the form of pre-existing fractures, or in the form of sealed fractures. To date, investigation of sample heterogeneity, specifically tensile strength and strength anisotropy has focused on layered rocks, such as shales, sandstones and gneisses. Data is lacking on the effect of single planar heterogeneities, such as pre-existing fractures or stylolites, even though these frequently occur in geo-energy settings.

We have performed Brazilian Disc tests on limestone samples containing planar heterogeneities, investigating Brazilian test Strength (BtS) and the effects of orientation on strength. We used prefractured Indiana limestone to represent a planar heterogeneity without cohesion and Treuchtlinger Marmor samples with central stylolites to represent a planar heterogeneity of unknown strength (as an example of a sealed fracture). The planar discontinuity was set at different rotation angles of approximately 0–20–30–45–60–90⁰, where 90⁰ (steep angle) is parallel to the principal loading direction, and 0⁰ (low angle) to the horizontal axis of the sample. All experiments were filmed, and where possible Particle Image Velocimetry was used to determine internal particle motion. Moreover, we used a 2D Comsol model in which we simplified the stylolite surface as a sinusoid. The model was used to qualitatively determine how i) a different period of the sinusoid and ii) relative strength of sinusoid/matrix affect the results.

Our results show that all imperfect samples are weaker than intact samples. The 2D Comsol model indicates that the qualitative results remain unaffected by changing the period (assumed to be representative of roughness) of the cohesive heterogeneity, nor by the relative strength contrast: the location of the first fracture remains unaffected. For both heterogeneity types, the fracture patterns can be divided into four categories, with two clear endmembers, and a more diffusive subdivision in between.

For a cohesion-less heterogeneity:

• steep angles lead to frictional sliding along the interface, and only a small hypothesized permeability increase.
• Intermediate angles lead to a combination of tensile failure of the matrix and sliding along the interface, where for steeper angles more new fractures form which follow the path of the existing fracture.
• Low angles lead to closure of the old fracture and new tensile failure.

For a cohesive heterogeneity of unknown cohesion:

• Steep angles lead to intensive failure of the heterogeneous zone, attributed to the presence of a stress concentrator.
• Intermediate angles lead to partial failure along the heterogeneous zone, and the formation of new fractures in the matrix, potentially instigated by mode II failure to accommodate motion.
• Low angles lead to the formation of a new fracture plus opening within the heterogeneous zone.

These results imply that hydrofracture (i.e. creating tensile stresses) of a stylolite-rich zone will lead to more fractures than fractures in a homogeneous zone, where the orientation of the stylolites and bedding will control the orientation of the permeable pathways.

How to cite: Pluymakers, A., Barnhoorn, A., and Bakker, R.: Effect of a heterogeneity on tensile failure: interaction between fractures in a limestone, EGU General Assembly 2021, online, 19–30 Apr 2021, EGU21-5959, https://doi.org/10.5194/egusphere-egu21-5959, 2021.

A postmortem technique is introduced to investigate the fracture connectivity evolution under elevated confining pressures via a sensitivity analysis. Three Onagawa shale samples are deformed under brittle, ductile, and transition conditions, by increasing the confining pressures. Brittle deformation is characterized by longitudinal splitting of the sample at 3% axial strain, and the onset of transition from brittle to ductile deformation is between 4% ~ 5% axial strain. The ductile deformation is characterized by a distributed conjugate fracture network and strain hardening. In completion of the deformation, the samples are scanned in a commercially available X-ray CT machine. The grayscale values of the primary 2D images were reversed, stacked, and surface rendered to obtain the 3D volume distribution of the fractures. Reversing and surface rendering allowed the acquisition of volume and surface data of the fractures along with their direct visualization. Further, utilizing a residual analysis, the voxel value density distribution that fabricated the fracture network is extracted (Residual histogram). Thresholding of the residual histogram generated volume segments of the final fracture network demonstrating the sensitivity of the fracture network to the choice of threshold. Voxel volumes of fractures alone are obtained by thresholding post-peak voxel values of the residual histogram and consecutive post-peak thresholding shows that the generated volume segments of the fracture network can be utilized to interpret, possible nucleation sites after strain localization, propagation of fractures, and coalescence. Fracture connectivity is quantified by means of relative entropy from information theory, and the relative entropy of size distribution of fracture volumes showed that it is closer to zero with the fractures being well connected. Moreover, the cumulative fracture volume shows a power-law growth towards the failure after a unique threshold to each sample. These results have been validated by previous acoustic emission studies and a 4D tomographic investigation on strain localization of shale. Therefore, despite the postmortem nature of the investigation, the new technique has opened possibilities to investigate the fracture properties and their evolution under elevated confining pressures.

How to cite: Jayawickrama, E., Muto, J., Sasaki, O., and Nagahama, H.: A postmortem approach via X-ray computed tomography and thresholding to investigate fracture network evolution in Onagawa shale, EGU General Assembly 2021, online, 19–30 Apr 2021, EGU21-11553, https://doi.org/10.5194/egusphere-egu21-11553, 2021.

Stromboli volcano has experienced four sector collapses over the past 13 thousand years, resulting in the formation of the Sciara del Fuoco (SDF) horseshoe-shaped depression and an inferred NE / SW striking rift zone across the SDF and the western sector of the island. These events have resulted in the formation of steep depressions on the slopes on the volcano where episodes of instability are continuously being observed and recorded. This study aims to quantify the fracture density inside and outside the rift zone to identify potential damaged zones that could reduce the edifice strength and promote fracturing. In order to do so we have carried out a multiscale analysis, by integrating satellite observations, field work and seismic and electrical resistivity analyses on cm scales blocks belonging to 11 lava units from the main volcanic cycles that have built the volcano edifice, ie. Paleostromboli, Nestromboli and Vancori. 0.5 m resolution Pleiades satellite data has been first used to highlight 23635 distinct linear features across the island. Fracture density has been calculated using Fracpaq based on the Mauldon et al (2001) method to determine the average fracture density of a given area on the basis of the average length of drawn segments within a predetermined circular area. 41.8 % of total fracture density is found around intrusions and fissures, with the summit area and the slopes of SDF having the highest average fracture density of 5.279  . Density, porosity, P- wave velocity in dry and wet conditions and electrical resistivity (in wet conditions) were measured  via an ultrasonic pulse generator and acquisition system (Pundit) and an on purpose built measuring quadrupole on cm scale blocks of lavas collected from both within and outside the proposed rift zone to assess the physical state and the crack damage of the different lava units.  Preliminary results show that P-wave velocity between ~ 2.25 km/s < Vp < 5km/s decreases with porosity while there is high variability electrical resistivity with 21.7 < ρ < 590 Ohm * m. This is presumably due to the lavas texture and the variable content of bubble/vesicles porosity and crack damage, that is reflected by an effective overall porosity between 0 and 9 %. Higher porosity is generally mirrored by lower p-wave velocity values. Neostromboli blocks show the most variability in both P-wave velocity and electrical resistivity. Further work will assess crack density throughout optical analyses and systematically investigate the UCS and elastic moduli. This integrated approach is expected to provide a multiscale fracture density and allow to develop further laboratory testing on how slip surfaces can evolve to a flank collapse at Stromboli.

How to cite: Alcock, T., Vinciguerra, S., Benson, P., and Vagnon, F.: Multiscale fracture density analysis at Stromboli Volcano, Italy: implications to flank stability, EGU General Assembly 2021, online, 19–30 Apr 2021, EGU21-8872, https://doi.org/10.5194/egusphere-egu21-8872, 2021.

The Fracture Induced Electromagnetic Radiation (FEMR) technique has gradually progressed in the past decade as a useful geophysical tool to determine the direction and magnitude of recent crustal stresses, visualize the modification and realignment of stresses inside tunnels thus proving to be an important precursor for geohazards, earthquake forecasting, as well as delineate landslide-prone slip planes in unstable regions. Its working principle is based on the generation of geogenic electromagnetic radiation emanating from the brittle rock bodies that are fractured being subjected to an incremental increase of the differential stress in the near-surface of the Earth’s crust. The “Process zone” at the fractured crack tip contains numerous microcracks which subsequently creates dipoles due to the polarization of charges on such microcrack tips which rapidly oscillates emitting FEMR waves of frequencies between KHz to MHz range. The coalescence of the microcracks eventually leads to a macro failure dampening the amplitude of the FEMR pulses. The attenuation of FEMR pulses is comparatively lesser than seismic waves making it a more efficient precursor to potential tectonic activities indicating an upcoming earthquake a few hours/days before the actual event. In the current study, we have attempted to exploit this technique to identify the locations of the potential active faults across the tectonically active Narmada-Son Lineament (NSL), Central India. Although the first tectonic stage involved rifting and formation of the NSL during the Precambrian time, the rifting continued at least till the time of Gondwana deposition. Later, tectonic inversion took place as a result of the collision between the Indian and the Eurasian plate resulting in reverse reactivation of the faults. Episodic reverse movement along NSL caused recurrent earthquakes and linear disposition of the sediments that were deposited at the foothills of the Satpura Horst. Although the origin of East-West trending NSL dates back to the Precambrian time, it is very much tectonically active as manifested by recent earthquakes. The study has been conducted by taking linear FEMR readings across 3 traverses along the NSL which on analysis provides an idea about the potential active faults, their locations, and frequency of occurrence. The accumulation of strain in the brittle rocks that can eventually lead to a macro failure is demarcated as an anomalous increase in the amplitude of the FEMR pulses indicative of an upcoming tectonic episode in the region. To further corroborate the analysis, we have attempted to determine the neo-tectonic activity in the region by calculating the morphometric parameters across the Khandwa-Itarsi-Jabalpur region, Central India. Finally, we attempt to comment on the tectonic evolution of Central India in the recent past. We also encourage researchers to adapt the novel technique of FEMR which is swift, affordable, and feasible compared to conventional techniques deployed to survey the active tectonics of a region.

How to cite: Das, S. and Mallik, J.: Delineating active faults across the Narmada-Son Lineament (NSL), India using the technique of Fracture Induced Electromagnetic Radiation (FEMR), EGU General Assembly 2021, online, 19–30 Apr 2021, EGU21-2370, https://doi.org/10.5194/egusphere-egu21-2370, 2021.

Rock mass damaging has become a topic of great interest in the engineering-geology research community during the last decades as it can significantly influence slopes stability. In this sense, the study of mechanics and dynamics of jointed rock masses represents a challenge because it will allow to better understand how external continuous and transient stressors can influence the short- to long-term stability controlling their pre-failure behavior. Consequently, the detection of permanent changes in physical and mechanical parameters, due to periodic or transient stressors, is an important target to mitigate the related geological risk as it can potentially lead rock masses to failure, especially when infrastructures and natural or cultural heritages are exposed elements. In this framework, the Acuto field laboratory (Central Italy) has been designed and implemented in 2016 within an abandoned quarry by employing an integrated geotechnical and geophysical monitoring system, with the aim of investigating how natural and anthropic conditioning factors could lead fractured rock masses to failure. The integrated monitoring system, which is installed on a potentially unstable 20-m³ jointed rock block, is composed of several strain devices (i.e., strain gauges -SG- and jointmeters -JM-), one fully equipped weather station, one rock thermometer, eight high-sensitivity microseismic uniaxial accelerometers and optical and InfraRed Thermal cameras. The acquisition of long-term monitoring time-series, coupling multimethodological approaches, allowed to establish cause-to-effect relationships among different environmental stressors and induced strain effects, highlighting the continuous action of thermal stresses on rock mass deformations both at the daily and seasonal timescales. In fact, while the analysis of thermal and strain monitoring data allowed to characterize the cyclic contraction and relaxation response of major rock fractures and microcracks to temperature fluctuations, the microseismic monitoring array was able to detect during thermal transient (i.e., freezing conditions) the occurrence of microseismic emissions potentially related to the genesis or progressive growth of pre-existing cracks.

Starting from 2018, experimental activities at the Acuto field lab are supported by the “Dipartimento di Eccellenza” project of the Italian Ministry of Education Universities and Research funds attributed to the Department of Earth Sciences of the University of Rome “Sapienza”.  In this framework, the Acuto filed laboratory will undergo a structural upgrade that will be aimed at the investigation of two new sectors of the abandoned quarry. These new sectors will be instrumented with innovative thermal profiles probe, fiber Brag grating sensors and traditional SG and JM for detailed stress-strain monitoring, acoustic emission sensors and high-frequency and low-frequency geophones for ambient seismic noise monitoring and microseismic events detection as well as accelerometers for evaluating the rock mass response in the case of seismic shaking. The main goal of such an improvement will be both technical and methodological, and will shed light on the application of integrated geophysical and geotechnical monitoring approaches in investigating the multiscale rock mass damaging process as well as the detection of rock mass failure precursors by using non-conventional combinations and configurations of geotechnical and broad-band geophysical devices.

How to cite: Grechi, G., D'Angiò, D., Fiorucci, M., Iannucci, R., Lenti, L., Marmoni, G. M., and Martino, S.: Integrated geophysical and geotechnical monitoring for multiscale rock mass damaging investigation at the Acuto Field-Lab (Italy), EGU General Assembly 2021, online, 19–30 Apr 2021, EGU21-9403, https://doi.org/10.5194/egusphere-egu21-9403, 2021.

In 2009, a large-scale landslide was triggered by typhoon rainfall and buried an entire village, which named Hsiaolin and located in Taiwan. 
After that, Soil and Water Conservation Bureau (SWCB) has promoted a national project for the prevention work of large-scale landslide. The national project includes with the investigation of potential area, the design of monitoring system, and the design of warning system, etc. 
The investigation of potential large-sclae landslide was based on the  digital elevation model with 1 meter resolution. However, the investigation of the underground was lack and not clear enough. Therefore, the specific landslide's body is hardly to estimate and it causes difficulty in follow-up works. 
This study applied two methods to investigate the scenario of slope failure. The first method is based on the limited equilibrium method, which proposed by Yoshino and Uchida (2019). The method was used to search the specific region of unstable slope based on a series of high-resolution digital elevation models. After the specific region of unstable slope was confirmed, the landslide can be simulated by a numerical model, which this study proposed to represent the entire landslide process from occurrence to post-failure . 
These proposed methods were applied at Baolai area, south Taiwan to track the evolution of the potential area. The failure scenario could be evaluated by the proposed numerical model. By this study, the investigation of underground can be evaluated and these results are very important information for the design of monitoring system.

How to cite: Tsai, Y. and Lee, W.: Application of the high-resolution digital elevation model in an integrated numerical method for the occurrence mechanism and post-failure behavior of the landslide, EGU General Assembly 2021, online, 19–30 Apr 2021, EGU21-10764, https://doi.org/10.5194/egusphere-egu21-10764, 2021.

In the pre-feasibility study stage, only a small amount of borehole data can be obtained. Since the available geological information is insufficient, the engineering geological conditions of the project can only be preliminarily and approximately estimated during this stage. In this study, we attempt to seek a method to make a preliminary analysis and evaluation of the stability of the surrounding rock masses of an underground rock carven project, which makes full and optimum use of the limited borehole data to accomplish the assessment of the investigated site. The basic information on rock fractures is extracted from the borehole Television logging data and the fracture extension directions are also determined. Providing that the cracks detected in the borehole would extend to the cavern area, the cracks with appropriate direction, larger width and larger hydraulic conductivity can be selected. These selected cracks are considered in the numerical model established using ANSYS, and the stability of surrounding rock of cavern is analyzed under this situation. In the absence of large amount of borehole data, this method, which set up an extreme case, can be used to analyze possible failure of rock mass under extreme adverse conduction in advance. In general, the proposed method for stability analysis could contribute to the design and construction practice of a tunnel project constructed in fractured rock masses.

How to cite: Tan, L. and Xiang, W.: Stability analysis of the surrounding rock of cavern under extreme situation based on limited borehole data., EGU General Assembly 2021, online, 19–30 Apr 2021, EGU21-230, https://doi.org/10.5194/egusphere-egu21-230, 2021.